DNA Extraction Method Affects the Detection of a Fungal Pathogen in Formalin-Fixed Specimens Using qPCR

Museum collections provide indispensable repositories for obtaining information about the historical presence of disease in wildlife populations. The pathogenic amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has played a significant role in global amphibian declines, and examining preserved specimens for Bd can improve our understanding of its emergence and spread. Quantitative PCR (qPCR) enables Bd detection with minimal disturbance to amphibian skin and is significantly more sensitive to detecting Bd than histology; therefore, developing effective qPCR methodologies for detecting Bd DNA in formalin-fixed specimens can provide an efficient and effective approach to examining historical Bd emergence and prevalence. Techniques for detecting Bd in museum specimens have not been evaluated for their effectiveness in control specimens that mimic the conditions of animals most likely to be encountered in museums, including those with low pathogen loads. We used American bullfrogs (Lithobates catesbeianus) of known infection status to evaluate the success of qPCR to detect Bd in formalin-fixed specimens after three years of ethanol storage. Our objectives were to compare the most commonly used DNA extraction method for Bd (PrepMan, PM) to Macherey-Nagel DNA FFPE (MN), test optimizations for Bd detection with PM, and provide recommendations for maximizing Bd detection. We found that successful detection is relatively high (80–90%) when Bd loads before formalin fixation are high, regardless of the extraction method used; however, at lower infection levels, detection probabilities were significantly reduced. The MN DNA extraction method increased Bd detection by as much as 50% at moderate infection levels. Our results indicate that, for animals characterized by lower pathogen loads (i.e., those most commonly encountered in museum collections), current methods may underestimate the proportion of Bd-infected amphibians. Those extracting DNA from archived museum specimens should ensure that the techniques they are using are known to provide high-quality throughput DNA for later analysis.     

A Bayesian Framework to Identify Methylcytosines from High-Throughput Bisulfite Sequencing Data

High-throughput bisulfite sequencing technologies have provided a comprehensive and well-fitted way to investigate DNA methylation at single-base resolution. However, there are substantial bioinformatic challenges to distinguish precisely methylcytosines from unconverted cytosines based on bisulfite sequencing data. The challenges arise, at least in part, from cell heterozygosis caused by multicellular sequencing and the still limited number of statistical methods that are available for methylcytosine calling based on bisulfite sequencing data. Here, we present an algorithm, termed Bycom, a new Bayesian model that can perform methylcytosine calling with high accuracy. Bycom considers cell heterozygosis along with sequencing errors and bisulfite conversion efficiency to improve calling accuracy. Bycom performance was compared with the performance of Lister, the method most widely used to identify methylcytosines from bisulfite sequencing data. The results showed that the performance of Bycom was better than that of Lister for data with high methylation levels. Bycom also showed higher sensitivity and specificity for low methylation level samples (<1%) than Lister. A validation experiment based on reduced representation bisulfite sequencing data suggested that Bycom had a false positive rate of about 4% while maintaining an accuracy of close to 94%. This study demonstrated that Bycom had a low false calling rate at any methylation level and accurate methylcytosine calling at high methylation levels. Bycom will contribute significantly to studies aimed at recalibrating the methylation level of genomic regions based on the presence of methylcytosines.     

Sequencing Degraded DNA from Non-Destructively Sampled Museum Specimens for RAD-Tagging and Low-Coverage Shotgun Phylogenetics

Ancient and archival DNA samples are valuable resources for the study of diverse historical processes. In particular, museum specimens provide access to biotas distant in time and space, and can provide insights into ecological and evolutionary changes over time. However, archival specimens are difficult to handle; they are often fragile and irreplaceable, and typically contain only short segments of denatured DNA. Here we present a set of tools for processing such samples for state-of-the-art genetic analysis. First, we report a protocol for minimally destructive DNA extraction of insect museum specimens, which produced sequenceable DNA from all of the samples assayed. The 11 specimens analyzed had fragmented DNA, rarely exceeding 100 bp in length, and could not be amplified by conventional PCR targeting the mitochondrial cytochrome oxidase I gene. Our approach made these samples amenable to analysis with commonly used next-generation sequencing-based molecular analytic tools, including RAD-tagging and shotgun genome re-sequencing. First, we used museum ant specimens from three species, each with its own reference genome, for RAD-tag mapping. Were able to use the degraded DNA sequences, which were sequenced in full, to identify duplicate reads and filter them prior to base calling. Second, we re-sequenced six Hawaiian Drosophila species, with millions of years of divergence, but with only a single available reference genome. Despite a shallow coverage of 0.37±0.42 per base, we could recover a sufficient number of overlapping SNPs to fully resolve the species tree, which was consistent with earlier karyotypic studies, and previous molecular studies, at least in the regions of the tree that these studies could resolve. Although developed for use with degraded DNA, all of these techniques are readily applicable to more recent tissue, and are suitable for liquid handling automation.     

碳酸钙沉淀法回收琼脂糖凝胶中DNA的探讨

采用碳酸钙沉淀法回收琼脂糖凝胶中的DNA,达到分离纯化目的,回收后的DNA可用于重组、PCR等研究。首先将含有目的DNA的琼脂糖凝胶用Nal溶液融解,然后加入cacl2,和NaHCO3,生成CaCO3,沉淀,DNA与cac03形成复合物,通过离心分离出沉淀复合物,利用稀酸溶解沉淀,再用无水乙醇沉降,即可回收目标DNA。利用该方法回收了质粒、毛白杨和转基因羊基因组DNA,同收率为20％～50％,0D260／OD280,为1．7～19,最大回收了21kb片段,最小回收250bp片段,回收后的DNA样品进行了PCR扩增和限制性内切酶反应,PCR可以扩增出目的片段,同时限制性内切酶可以将回收后的DNA切开,表明DNA质量良好。利用碳酸钙沉淀法可以回收琼脂糖凝胶中的DNA,此法简单、易行,较为有效。     

A Method Enabling High-Throughput Sequencing of Human Cytomegalovirus Complete Genomes from Clinical Isolates

Human cytomegalovirus (HCMV) is a ubiquitous virus that can cause serious sequelae in immunocompromised patients and in the developing fetus. The coding capacity of the 235 kbp genome is still incompletely understood, and there is a pressing need to characterize genomic contents in clinical isolates. In this study, a procedure for the high-throughput generation of full genome consensus sequences from clinical HCMV isolates is presented. This method relies on low number passaging of clinical isolates on human fibroblasts, followed by digestion of cellular DNA and purification of viral DNA. After multiple displacement amplification, highly pure viral DNA is generated. These extracts are suitable for high-throughput next-generation sequencing and assembly of consensus sequences. Throughout a series of validation experiments, we showed that the workflow reproducibly generated consensus sequences representative for the virus population present in the original clinical material. Additionally, the performance of 454 GS FLX and/or Illumina Genome Analyzer datasets in consensus sequence deduction was evaluated. Based on assembly performance data, the Illumina Genome Analyzer was the platform of choice in the presented workflow. Analysis of the consensus sequences derived in this study confirmed the presence of gene-disrupting mutations in clinical HCMV isolates independent from in vitro passaging. These mutations were identified in genes RL5A, UL1, UL9, UL111A and UL150. In conclusion, the presented workflow provides opportunities for high-throughput characterization of complete HCMV genomes that could deliver new insights into HCMV coding capacity and genetic determinants of viral tropism and pathogenicity.     

HTSstation: A Web Application and Open-Access Libraries for High-Throughput Sequencing Data Analysis

The HTSstation analysis portal is a suite of simple web forms coupled to modular analysis pipelines for various applications of High-Throughput Sequencing including ChIP-seq, RNA-seq, 4C-seq and re-sequencing. HTSstation offers biologists the possibility to rapidly investigate their HTS data using an intuitive web application with heuristically pre-defined parameters. A number of open-source software components have been implemented and can be used to build, configure and run HTS analysis pipelines reactively. Besides, our programming framework empowers developers with the possibility to design their own workflows and integrate additional third-party software. The HTSstation web application is accessible at http://htsstation.epfl.ch.     

SLAF-seq: An Efficient Method of Large-Scale De Novo SNP Discovery and Genotyping Using High-Throughput Sequencing

Large-scale genotyping plays an important role in genetic association studies. It has provided new opportunities for gene discovery, especially when combined with high-throughput sequencing technologies. Here, we report an efficient solution for large-scale genotyping. We call it specific-locus amplified fragment sequencing (SLAF-seq). SLAF-seq technology has several distinguishing characteristics: i) deep sequencing to ensure genotyping accuracy; ii) reduced representation strategy to reduce sequencing costs; iii) pre-designed reduced representation scheme to optimize marker efficiency; and iv) double barcode system for large populations. In this study, we tested the efficiency of SLAF-seq on rice and soybean data. Both sets of results showed strong consistency between predicted and practical SLAFs and considerable genotyping accuracy. We also report the highest density genetic map yet created for any organism without a reference genome sequence, common carp in this case, using SLAF-seq data. We detected 50,530 high-quality SLAFs with 13,291 SNPs genotyped in 211 individual carp. The genetic map contained 5,885 markers with 0.68 cM intervals on average. A comparative genomics study between common carp genetic map and zebrafish genome sequence map showed high-quality SLAF-seq genotyping results. SLAF-seq provides a high-resolution strategy for large-scale genotyping and can be generally applicable to various species and populations.     

A Method for Preparing DNA Sequencing Templates Using a DNA-Binding Microplate

A DNA-binding matrix was immobilized on the surface of a 96-well microplate and used for plasmid DNA preparation for DNA sequencing. The same DNA-binding plate was used for bacterial growth, cell lysis, DNA purification, and storage. In a single step using one buffer, bacterial cells were lysed by enzymes, and released DNA was captured on the plate simultaneously. After two wash steps, DNA was eluted and stored in the same plate. Inclusion of phosphates in the culture medium was found to enhance the yield of plasmid significantly. Purified DNA samples were used successfully in DNA sequencing with high consistency and reproducibility. Eleven vectors and nine libraries were tested using this method. In 10 μl sequencing reactions using 3 μl sample and 0.25 μl BigDye Terminator v3.1, the results from a 3730xl sequencer gave a success rate of 90–95% and read-lengths of 700 bases or more. The method is fully automatable and convenient for manual operation as well. It enables reproducible, high-throughput, rapid production of DNA with purity and yields sufficient for high-quality DNA sequencing at a substantially reduced cost.     

Detecting DNA Modifications from SMRT Sequencing Data by Modeling Sequence Context Dependence of Polymerase Kinetic

DNA modifications such as methylation and DNA damage can play critical regulatory roles in biological systems. Single molecule, real time (SMRT) sequencing technology generates DNA sequences as well as DNA polymerase kinetic information that can be used for the direct detection of DNA modifications. We demonstrate that local sequence context has a strong impact on DNA polymerase kinetics in the neighborhood of the incorporation site during the DNA synthesis reaction, allowing for the possibility of estimating the expected kinetic rate of the enzyme at the incorporation site using kinetic rate information collected from existing SMRT sequencing data (historical data) covering the same local sequence contexts of interest. We develop an Empirical Bayesian hierarchical model for incorporating historical data. Our results show that the model could greatly increase DNA modification detection accuracy, and reduce requirement of control data coverage. For some DNA modifications that have a strong signal, a control sample is not even needed by using historical data as alternative to control. Thus, sequencing costs can be greatly reduced by using the model. We implemented the model in a R package named seqPatch, which is available at https://github.com/zhixingfeng/seqPatch.     

Proboscidean DNA from Museum and Fossil Specimens: An Assessment of Ancient DNA Extraction and Amplification Techniques

Applications of reliable DNA extraction and amplification techniques to postmortem samples are critical to ancient DNA research. Commonly used methods for isolating DNA from ancient material were tested and compared using both soft tissue and bones from fossil and contemporary museum proboscideans. DNAs isolated using three principal methods served as templates in subsequent PCR amplifications, and the PCR products were directly sequenced. Authentication of the ancient origin of obtained nucleotide sequences was established by demonstrating reproducibility under a blind testing system and by phylogenetic analysis. Our results indicate that ancient samples may respond differently to extraction buffers or purification procedures, and no single method was universally successful. A CTAB buffer method, modified from plant DNA extraction protocols, was found to have the highest success rate. Nested PCR was shown to be a reliable approach to amplify ancient DNA templates that failed in primary amplification.     

Tn-Seq Explorer: A Tool for Analysis of High-Throughput Sequencing Data of Transposon Mutant Libraries

Tn-seq is a high throughput technique for analysis of transposon mutant libraries. Tn-seq Explorer was developed as a convenient and easy-to-use package of tools for exploration of the Tn-seq data. In a typical application, the user will have obtained a collection of sequence reads adjacent to transposon insertions in a reference genome. The reads are first aligned to the reference genome using one of the tools available for this task. Tn-seq Explorer reads the alignment and the gene annotation, and provides the user with a set of tools to investigate the data and identify possibly essential or advantageous genes as those that contain significantly low counts of transposon insertions. Emphasis is placed on providing flexibility in selecting parameters and methodology most appropriate for each particular dataset. Tn-seq Explorer is written in Java as a menu-driven, stand-alone application. It was tested on Windows, Mac OS, and Linux operating systems. The source code is distributed under the terms of GNU General Public License. The program and the source code are available for download at http://www.cmbl.uga.edu/downloads/programs/Tn_seq_Explorer/ and https://github.com/sina-cb/Tn-seqExplorer.     

A High-Throughput Screening Method to Reengineer DNA Polymerases for Random Mutagenesis

A screening system for directed evolution of DNA polymerases employing a fluorescent Scorpion probe as a reporter has been developed. The screening system has been validated in a directed evolution experiment of a distributive polymerase from the Y-polymerase family (Dpo4 from Sulfolobus solfataricus) which was improved in elongation efficiency of consecutive mismatches. The engineering campaign yielded improved Dpo4 polymerase variants one of which was successfully benchmarked in a sequence saturation mutagenesis experiment especially with regard to the desirable consecutive transversion mutations (>2.5-fold increase in frequency relative to a reference library prepared with Dpo4 WT). The Scorpion probe screening system enables to reengineer polymerases with low processivity and fidelity, and no secondary activities (i.e. exonuclease activity or strand displacement activity) to match demands in diversity generation for directed protein evolution.     

DDBJ Read Annotation Pipeline: A Cloud Computing-Based Pipeline for High-Throughput Analysis of Next-Generation Sequencing Data

High-performance next-generation sequencing (NGS) technologies are advancing genomics and molecular biological research. However, the immense amount of sequence data requires computational skills and suitable hardware resources that are a challenge to molecular biologists. The DNA Data Bank of Japan (DDBJ) of the National Institute of Genetics (NIG) has initiated a cloud computing-based analytical pipeline, the DDBJ Read Annotation Pipeline (DDBJ Pipeline), for a high-throughput annotation of NGS reads. The DDBJ Pipeline offers a user-friendly graphical web interface and processes massive NGS datasets using decentralized processing by NIG supercomputers currently free of charge. The proposed pipeline consists of two analysis components: basic analysis for reference genome mapping and de novo assembly and subsequent high-level analysis of structural and functional annotations. Users may smoothly switch between the two components in the pipeline, facilitating web-based operations on a supercomputer for high-throughput data analysis. Moreover, public NGS reads of the DDBJ Sequence Read Archive located on the same supercomputer can be imported into the pipeline through the input of only an accession number. This proposed pipeline will facilitate research by utilizing unified analytical workflows applied to the NGS data. The DDBJ Pipeline is accessible at http://p.ddbj.nig.ac.jp/.     

Modelling the Evolutionary Dynamics of Viruses within Their Hosts: A Case Study Using High-Throughput Sequencing

Uncovering how natural selection and genetic drift shape the evolutionary dynamics of virus populations within their hosts can pave the way to a better understanding of virus emergence. Mathematical models already play a leading role in these studies and are intended to predict future emergences. Here, using high-throughput sequencing, we analyzed the within-host population dynamics of four Potato virus Y (PVY) variants differing at most by two substitutions involved in pathogenicity properties. Model selection procedures were used to compare experimental results to six hypotheses regarding competitiveness and intensity of genetic drift experienced by viruses during host plant colonization. Results indicated that the frequencies of variants were well described using Lotka-Volterra models where the competition coefficients β_ij exerted by variant j on variant i are equal to their fitness ratio, r_j/r_i. Statistical inference allowed the estimation of the effect of each mutation on fitness, revealing slight (s = −0.45%) and high (s = −13.2%) fitness costs and a negative epistasis between them. Results also indicated that only 1 to 4 infectious units initiated the population of one apical leaf. The between-host variances of the variant frequencies were described using Dirichlet-multinomial distributions whose scale parameters, closely related to the fixation index F_ST, were shown to vary with time. The genetic differentiation of virus populations among plants increased from 0 to 10 days post-inoculation and then decreased until 35 days. Overall, this study showed that mathematical models can accurately describe both selection and genetic drift processes shaping the evolutionary dynamics of viruses within their hosts.     

Sequence to Medical Phenotypes: A Framework for Interpretation of Human Whole Genome DNA Sequence Data

High throughput sequencing has facilitated a precipitous drop in the cost of genomic sequencing, prompting predictions of a revolution in medicine via genetic personalization of diagnostic and therapeutic strategies. There are significant barriers to realizing this goal that are related to the difficult task of interpreting personal genetic variation. A comprehensive, widely accessible application for interpretation of whole genome sequence data is needed. Here, we present a series of methods for identification of genetic variants and genotypes with clinical associations, phasing genetic data and using Mendelian inheritance for quality control, and providing predictive genetic information about risk for rare disease phenotypes and response to pharmacological therapy in single individuals and father-mother-child trios. We demonstrate application of these methods for disease and drug response prognostication in whole genome sequence data from twelve unrelated adults, and for disease gene discovery in one father-mother-child trio with apparently simplex congenital ventricular arrhythmia. In doing so we identify clinically actionable inherited disease risk and drug response genotypes in pre-symptomatic individuals. We also nominate a new candidate gene in congenital arrhythmia, ATP2B4, and provide experimental evidence of a regulatory role for variants discovered using this framework.     

A Non-radioactive DNA Sequencing Method Using Biotinylated Dideoxynucleoside Triphosphates and {Delta}Tth DNA Polymerase

We synthesized a set of four biotinylated dideoxynucleosidetriphosphates (biotin-9-ddNTPs) and optimized the reaction conditionsfor non-radioactive cycle sequencing using modified Tth DNApolymerase (Tth) and a chemiluminescent detection system. Theresulting sequencing ladders showed lower background comparedto those with the conventional non-radioactive sequencing methodwhich uses 5'-biotinylated primers, especially when PCR productswere analysed. With our method, DNA sequences can be determinedat any primer positions without preparing 5'-biotinylated primersfor dideoxy chaintermination.     

The Use of High-Throughput DNA Sequencing in the Investigation of Antigenic Variation: Application to Neisseria Species

Antigenic variation occurs in a broad range of species. This process resembles gene conversion in that variant DNA is unidirectionally transferred from partial gene copies (or silent loci) into an expression locus. Previous studies of antigenic variation have involved the amplification and sequencing of individual genes from hundreds of colonies. Using the pilE gene from Neisseria gonorrhoeae we have demonstrated that it is possible to use PCR amplification, followed by high-throughput DNA sequencing and a novel assembly process, to detect individual antigenic variation events. The ability to detect these events was much greater than has previously been possible. In N. gonorrhoeae most silent loci contain multiple partial gene copies. Here we show that there is a bias towards using the copy at the 3′ end of the silent loci (copy 1) as the donor sequence. The pilE gene of N. gonorrhoeae and some strains of Neisseria meningitidis encode class I pilin, but strains of N. meningitidis from clonal complexes 8 and 11 encode a class II pilin. We have confirmed that the class II pili of meningococcal strain FAM18 (clonal complex 11) are non-variable, and this is also true for the class II pili of strain NMB from clonal complex 8. In addition when a gene encoding class I pilin was moved into the meningococcal strain NMB background there was no evidence of antigenic variation. Finally we investigated several members of the opa gene family of N. gonorrhoeae, where it has been suggested that limited variation occurs. Variation was detected in the opaK gene that is located close to pilE, but not at the opaJ gene located elsewhere on the genome. The approach described here promises to dramatically improve studies of the extent and nature of antigenic variation systems in a variety of species.     

一种高效快速的鱼类标本基因组DNA提取方法

以95%酒精保存的黄鳝(M onopterus albus)和斑鳢(Channa maculates)标本为材料,采用先沉降DNA再去除杂质的方法从鱼类标本中提取基因组DNA。基因组DNA的琼脂糖凝胶电泳和紫外分光光度法检测以及PCR扩增结果显示,本方法提取的鱼类基因组DNA的电泳主带清晰明亮;A260/A280值在1.7830-2.0144之间;PCR扩增产物条带清晰明亮,且单一整齐没有拖带,表明本方法可从酒精保存的鱼类标本中提取比较纯净的DNA,能够满足一般分子生物学试验需要。与传统苯酚/氯仿法相比,本方法操作简单快速,避免了苯酚等物质对后续实验的影响,可作为一种常规动物组织DNA提取方法。     

一种快速简便构建DNA病毒测序文库的方法

为了对1株中国棉铃虫核型多角体缺失病毒HZ-9进行基因组测序,采用了一种新的方法,通过超声波振断HaBacHZ9细菌人工染色体质粒（bacterial artificial chromosome plasmid,Bacmid）基因组DNA,用Taq酶在DNA片段两端加腺噤呤A,胶回收后得到预期的1—2kb的DNA片段,然后与pGEM-Teasy载体连接,构建了中国棉铃虫缺失病毒HaBacHZ9的亚克隆文库。结果随机挑选10个克隆子酶切分析,显示9个克隆子有1500bp左右的插入片断,并对HaBaeHZ9进行了全基因组测序。结论成功构建了HaBaeHZ9的DNA测序文库,为HZ-9功能基因组学研究奠定了基础,这是一种简单快速的构建DNA病毒测序文库的方法。