A combinatorial approach to the restriction of a mouse genome

cuRRBS: simple and robust evaluation of enzyme combinations for reduced representation approaches

DNA methylation is an important epigenetic modification in many species that is critical for development, and implicated in ageing and many complex diseases, such as cancer. Many cost-effective genome-wide analyses of DNA modifications rely on restriction enzymes capable of digesting genomic DNA at defined sequence motifs. There are hundreds of restriction enzyme families but few are used to date, because no tool is available for the systematic evaluation of restriction enzyme combinations that can enrich for certain sites of interest in a genome. Herein, we present customised Reduced Representation Bisulfite Sequencing (cuRRBS), a novel and easy-to-use computational method that solves this problem. By computing the optimal enzymatic digestions and size selection steps required, cuRRBS generalises the traditional MspI-based Reduced Representation Bisulfite Sequencing (RRBS) protocol to all restriction enzyme combinations. In addition, cuRRBS estimates the fold-reduction in sequencing costs and provides a robustness value for the personalised RRBS protocol, allowing users to tailor the protocol to their experimental needs. Moreover, we show in silico that cuRRBS-defined restriction enzymes consistently out-perform MspI digestion in many biological systems, considering both CpG and CHG contexts. Finally, we have validated the accuracy of cuRRBS predictions for single and double enzyme digestions using two independent experimental datasets.

Bidirectional Retroviral Integration Site PCR Methodology and Quantitative Data Analysis Workflow

Integration Site (IS) assays are a critical component of the study of retroviral integration sites and their biological significance. In recent retroviral gene therapy studies, IS assays, in combination with next-generation sequencing, have been used as a cell-tracking tool to characterize clonal stem cell populations sharing the same IS. For the accurate comparison of repopulating stem cell clones within and across different samples, the detection sensitivity, data reproducibility, and high-throughput capacity of the assay are among the most important assay qualities. This work provides a detailed protocol and data analysis workflow for bidirectional IS analysis. The bidirectional assay can simultaneously sequence both upstream and downstream vector-host junctions. Compared to conventional unidirectional IS sequencing approaches, the bidirectional approach significantly improves IS detection rates and the characterization of integration events at both ends of the target DNA. The data analysis pipeline described here accurately identifies and enumerates identical IS sequences through multiple steps of comparison that map IS sequences onto the reference genome and determine sequencing errors. Using an optimized assay procedure, we have recently published the detailed repopulation patterns of thousands of Hematopoietic Stem Cell (HSC) clones following transplant in rhesus macaques, demonstrating for the first time the precise time point of HSC repopulation and the functional heterogeneity of HSCs in the primate system. The following protocol describes the step-by-step experimental procedure and data analysis workflow that accurately identifies and quantifies identical IS sequences.

Focussing reduced representation CpG sequencing through judicious restriction enzyme choice

Current restriction enzyme based reduced representation methylation analyses aim for limited, but unbiased, methylome coverage. As the current best estimate suggests that only ~20% of CpGs are dynamically regulated, we characterised the CpG and genomic context surrounding all suitable restriction enzyme sites to identify those that were located in regions rich in dynamically methylated CpGs. The restriction-site distributions for MspI, BstUI, and HhaI were non-random. CpGs in CGI and shelf+shore could be enriched, particularly in gene bodies for all genomic regions, promoters (TSS1500, TSS200), intra- (1st exon, gene body, 3'UTR, 5'UTR) and inter-genic regions. HpyCH4IV enriched CpG elements in the open sea for all genomic elements. Judicious restriction enzyme choice improves the focus of reduced representation approaches by avoiding the monopolization of read coverage by genomic regions that are irrelevant, unwanted or difficult to map, and only sequencing the most informative fraction of CpGs.