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.

Nondestructive enzymatic deamination enables single-molecule long-read amplicon sequencing for the determination of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution.

The predominant methodology for DNA methylation analysis relies on the chemical deamination by sodium bisulfite of unmodified cytosine to uracil to permit the differential readout of methylated cytosines. Bisulfite treatment damages the DNA, leading to fragmentation and loss of long-range methylation information. To overcome this limitation of bisulfite-treated DNA, we applied a new enzymatic deamination approach, termed enzymatic methyl-seq (EM-seq), to long-range sequencing technologies. Our methodology, named long-read enzymatic modification sequencing (LR-EM-seq), preserves the integrity of DNA, allowing long-range methylation profiling of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) over multikilobase length of genomic DNA. When applied to known differentially methylated regions (DMRs), LR-EM-seq achieves phasing of >5 kb, resulting in broader and better defined DMRs compared with that previously reported. This result showed the importance of phasing methylation for biologically relevant questions and the applicability of LR-EM-seq for long-range epigenetic analysis at single-molecule and single-nucleotide resolution.

Absence of founder effect and evidence for adaptive divergence in a recently introduced insular population of white-tailed deer (Odocoileus virginianus).

Islands are generally colonized by few individuals which could lead to a founder effect causing loss of genetic diversity and rapid divergence by strong genetic drift. Insular conditions can also induce new selective pressures on populations. Here, we investigated the extent of genetic differentiation within a white-tailed deer (Odocoileus virginianus) population introduced on an island and its differentiation with its source mainland population. In response to their novel environmental conditions, introduced deer changed phenotypically from mainland individuals, therefore we investigated the genetic bases of the morphological differentiation. The study was conducted on Anticosti Island (Québec, Canada) where 220 individuals were introduced 120 years ago, resulting in a population size over 160,000 individuals. We used genotyping-by-sequencing (GBS) to generate 8,518 filtered high-quality SNPs and compared patterns of genetic diversity and differentiation between the continental and Anticosti Island populations. Clustering analyses indicated a single panmictic island population and no sign of isolation by distance. Our results revealed a weak, albeit highly significant, genetic differentiation between the Anticosti Island population and its source population (mean FST  = 0.005), which allowed a population assignment success of 93%. Also, the high genetic diversity maintained in the introduced population supports the absence of a strong founder effect due to the large number of founders followed by rapid population growth. We further used a polygenic approach to assess the genetic bases of the divergent phenotypical traits between insular and continental populations. We found loci related to muscular function and lipid metabolism, which suggested that these could be involved in local adaptation on Anticosti Island. We discuss these results in a harvest management context.

A Microbiome DNA Enrichment Method for Next-Generation Sequencing Sample Preparation.

"Microbiome" is used to describe the communities of microorganisms and their genes in a particular environment, including communities in association with a eukaryotic host or part of a host. One challenge in microbiome analysis concerns the presence of host DNA in samples. Removal of host DNA before sequencing results in greater sequence depth of the intended microbiome target population. This unit describes a novel method of microbial DNA enrichment in which methylated host DNA such as human genomic DNA is selectively bound and separated from microbial DNA before next-generation sequencing (NGS) library construction. This microbiome enrichment technique yields a higher fraction of microbial sequencing reads and improved read quality resulting in a reduced cost of downstream data generation and analysis. © 2016 by John Wiley & Sons, Inc.

Selective Depletion of Abundant RNAs to Enable Transcriptome Analysis of Low-Input and Highly Degraded Human RNA.

Ribosomal RNAs (rRNAs) are extremely abundant, often constituting 80% to 90% of total RNA. Since rRNA sequences are often not of interest in genomic RNA sequencing experiments, rRNAs can be removed from the sample before the library preparation step, in order to prevent the majority of the library and the majority of sequencing reads from being rRNA. Removal of rRNA can be especially challenging for low quality and formalin-fixed paraffin-embedded (FFPE) RNA samples due to the fragmented nature of these RNA molecules. The NEBNext rRNA Depletion Kit (Human/Mouse/Rat) depletes both cytoplasmic (5 S rRNA, 5.8 S rRNA, 18 S rRNA, and 28 S rRNA) and mitochondrial rRNA (12 S rRNA and 16 S rRNA) from total RNA preparations from human, mouse, and rat samples. Due to the high similarity among mammalian rRNA sequences, it is likely that rRNA depletion can also be achieved for other mammals but has not been empirically tested. This product is compatible with both intact and degraded RNA (e.g., FFPE RNA). The resulting rRNA-depleted RNA is suitable for RNA-seq, random-primed cDNA synthesis, or other downstream RNA analysis applications. Regardless of the quality or amount of input RNA, this method efficiently removes rRNA, while retaining non-coding and other non-poly(A) RNAs. The NEBNext rRNA Depletion Kit thus provides a more complete picture of the transcript repertoire than oligo d(T) poly(A) mRNA enrichment methods. © 2016 by John Wiley & Sons, Inc.

A method for selectively enriching microbial DNA from contaminating vertebrate host DNA.

DNA samples derived from vertebrate skin, bodily cavities and body fluids contain both host and microbial DNA; the latter often present as a minor component. Consequently, DNA sequencing of a microbiome sample frequently yields reads originating from the microbe(s) of interest, but with a vast excess of host genome-derived reads. In this study, we used a methyl-CpG binding domain (MBD) to separate methylated host DNA from microbial DNA based on differences in CpG methylation density. MBD fused to the Fc region of a human antibody (MBD-Fc) binds strongly to protein A paramagnetic beads, forming an effective one-step enrichment complex that was used to remove human or fish host DNA from bacterial and protistan DNA for subsequent sequencing and analysis. We report enrichment of DNA samples from human saliva, human blood, a mock malaria-infected blood sample and a black molly fish. When reads were mapped to reference genomes, sequence reads aligning to host genomes decreased 50-fold, while bacterial and Plasmodium DNA sequences reads increased 8-11.5-fold. The Shannon-Wiener diversity index was calculated for 149 bacterial species in saliva before and after enrichment. Unenriched saliva had an index of 4.72, while the enriched sample had an index of 4.80. The similarity of these indices demonstrates that bacterial species diversity and relative phylotype abundance remain conserved in enriched samples. Enrichment using the MBD-Fc method holds promise for targeted microbiome sequence analysis across a broad range of sample types.

Varicella outbreak among Afghan National Civil Order Police recruits-Herat Regional Military Training Center, Herat, Afghanistan, 2010.

In December 2010, an outbreak of varicella was reported among student recruits enrolled at the Afghan National Civil Order Police Herat Regional Military Training Center. The outbreak had an overall attack rate of 9.8% (31 of 316 recruits) with primary, secondary, and tertiary attack rates of 6.3% (20 of 316), 3.4% (10 of 296), and 0.35% (1 of 286). Fortunately, the outbreak did not lead to any deaths or serious complications. However, it significantly interfered with Afghan National Civil Order Police training by causing a loss of 378 person-days of training. Medical personnel from the Afghan National Police, DynCorp International, Government of the Islamic Republic of Afghanistan Ministry of Public Health, and NATO Training Mission-Afghanistan Herat Joint Medical Operation Cell joined together to control and characterize the outbreak and prepare and disseminate recommendations for preventing future outbreaks. Control measures were quickly implemented, but less than ideal. Varicella vaccine was not available in Afghanistan to immunize exposed recruits. The outbreak was reported to medical authorities through a slow and convoluted process. And the majority of varicella cases did not self-report for care. Rather, medical personnel diagnosed most cases only after recruits were directed to report for a physical examination.

The MspJI family of modification-dependent restriction endonucleases for epigenetic studies.

MspJI is a novel modification-dependent restriction endonuclease that cleaves at a fixed distance away from the modification site. Here, we present the biochemical characterization of several MspJI homologs, including FspEI, LpnPI, AspBHI, RlaI, and SgrTI. All of the enzymes specifically recognize cytosine C5 modification (methylation or hydroxymethylation) in DNA and cleave at a constant distance (N(12)/N(16)) away from the modified cytosine. Each displays its own sequence context preference, favoring different nucleotides flanking the modified cytosine. By cleaving on both sides of fully modified CpG sites, they allow the extraction of 32-base long fragments around the modified sites from the genomic DNA. These enzymes provide powerful tools for direct interrogation of the epigenome. For example, we show that RlaI, an enzyme that prefers (m)CWG but not (m)CpG sites, generates digestion patterns that differ between plant and mammalian genomic DNA, highlighting the difference between their epigenomic patterns. In addition, we demonstrate that deep sequencing of the digested DNA fragments generated from these enzymes provides a feasible method to map the modified sites in the genome. Altogether, the MspJI family of enzymes represent appealing tools of choice for method development in DNA epigenetic studies.

High sensitivity 5-hydroxymethylcytosine detection in Balb/C brain tissue.

DNA hydroxymethylation is a long known modification of DNA, but has recently become a focus in epigenetic research. Mammalian DNA is enzymatically modified at the 5(th) carbon position of cytosine (C) residues to 5-mC, predominately in the context of CpG dinucleotides. 5-mC is amenable to enzymatic oxidation to 5-hmC by the Tet family of enzymes, which are believed to be involved in development and disease. Currently, the biological role of 5-hmC is not fully understood, but is generating a lot of interest due to its potential as a biomarker. This is due to several groundbreaking studies identifying 5-hydroxymethylcytosine in mouse embryonic stem (ES) and neuronal cells. Research techniques, including bisulfite sequencing methods, are unable to easily distinguish between 5-mC and 5-hmC . A few protocols exist that can measure global amounts of 5-hydroxymethylcytosine in the genome, including liquid chromatography coupled with mass spectrometry analysis or thin layer chromatography of single nucleosides digested from genomic DNA. Antibodies that target 5-hydroxymethylcytosine also exist, which can be used for dot blot analysis, immunofluorescence, or precipitation of hydroxymethylated DNA, but these antibodies do not have single base resolution.In addition, resolution depends on the size of the immunoprecipitated DNA and for microarray experiments, depends on probe design. Since it is unknown exactly where 5-hydroxymethylcytosine exists in the genome or its role in epigenetic regulation, new techniques are required that can identify locus specific hydroxymethylation. The EpiMark 5-hmC and 5-mC Analysis Kit provides a solution for distinguishing between these two modifications at specific loci. The EpiMark 5-hmC and 5-mC Analysis Kit is a simple and robust method for the identification and quantitation of 5-methylcytosine and 5-hydroxymethylcytosine within a specific DNA locus. This enzymatic approach utilizes the differential methylation sensitivity of the isoschizomers MspI and HpaII in a simple 3-step protocol. Genomic DNA of interest is treated with T4-BGT, adding a glucose moeity to 5-hydroxymethylcytosine. This reaction is sequence-independent, therefore all 5-hmC will be glucosylated; unmodified or 5-mC containing DNA will not be affected. This glucosylation is then followed by restriction endonuclease digestion. MspI and HpaII recognize the same sequence (CCGG) but are sensitive to different methylation states. HpaII cleaves only a completely unmodified site: any modification (5-mC, 5-hmC or 5-ghmC) at either cytosine blocks cleavage. MspI recognizes and cleaves 5-mC and 5-hmC, but not 5-ghmC. The third part of the protocol is interrogation of the locus by PCR. As little as 20 ng of input DNA can be used. Amplification of the experimental (glucosylated and digested) and control (mock glucosylated and digested) target DNA with primers flanking a CCGG site of interest (100-200 bp) is performed. If the CpG site contains 5-hydroxymethylcytosine, a band is detected after glucosylation and digestion, but not in the non-glucosylated control reaction. Real time PCR will give an approximation of how much hydroxymethylcytosine is in this particular site. In this experiment, we will analyze the 5-hydroxymethylcytosine amount in a mouse Babl/C brain sample by end point PCR.

MmeI: a minimal Type II restriction-modification system that only modifies one DNA strand for host protection.

MmeI is an unusual Type II restriction enzyme that is useful for generating long sequence tags. We have cloned the MmeI restriction-modification (R-M) system and found it to consist of a single protein having both endonuclease and DNA methyltransferase activities. The protein comprises an amino-terminal endonuclease domain, a central DNA methyltransferase domain and C-terminal DNA recognition domain. The endonuclease cuts the two DNA strands at one site simultaneously, with enzyme bound at two sites interacting to accomplish scission. Cleavage occurs more rapidly than methyl transfer on unmodified DNA. MmeI modifies only the adenine in the top strand, 5'-TCCRAC-3'. MmeI endonuclease activity is blocked by this top strand adenine methylation and is unaffected by methylation of the adenine in the complementary strand, 5'-GTYGGA-3'. There is no additional DNA modification associated with the MmeI R-M system, as is required for previously characterized Type IIG R-M systems. The MmeI R-M system thus uses modification on only one of the two DNA strands for host protection. The MmeI architecture represents a minimal approach to assembling a restriction-modification system wherein a single DNA recognition domain targets both the endonuclease and DNA methyltransferase activities.

Efficient isolation of targeted Caenorhabditis elegans deletion strains using highly thermostable restriction endonucleases and PCR.

Reverse genetic approaches to understanding gene function would be greatly facilitated by increasing the efficiency of methods for isolating mutants without the reliance on a predicted phenotype. Established PCR-based methods of isolating deletion mutants are widely used for this purpose in Caenorhabditis elegans. However, these methods are inefficient at isolating small deletions. We report here a novel modification of PCR-based methods, employing thermostable restriction enzymes to block the synthesis of wild-type PCR product, so that only the deletion PCR product is amplified. This modification greatly increases the efficiency of isolating small targeted deletions in C.elegans. Using this method six new deletion strains were isolated from a small screen of approximately 400 000 haploid genomes, most with deletions <1.0 kb. Greater PCR detection sensitivity by this modification permitted approximately 10-fold greater pooling of DNA samples, reducing the effort and reagents required for screens. In addition, effective suppression of non-specific amplification allowed multiplexing with several independent primer pairs. The increased efficiency of this technique makes it more practical for small laboratories to undertake gene knock-out screens.