Fidelity of DNA polymerases in DNA amplification

Minimal twister sister-like self-cleaving ribozymes in the human genome revealed by deep mutational scanning

Despite their importance in a wide range of living organisms, self-cleaving ribozymes in the human genome are few and poorly studied. Here, we performed deep mutational scanning and covariance analysis of two previously proposed self-cleaving ribozymes (LINE-1 and OR4K15). We found that the regions essential for ribozyme activities are made of two short segments, with a total of 35 and 31 nucleotides only. The discovery makes them the simplest known self-cleaving ribozymes. Moreover, the essential regions are circular permutated with two nearly identical catalytic internal loops, supported by two stems of different lengths. These two self-cleaving ribozymes, which are shaped like lanterns, are similar to the catalytic regions of the twister sister ribozymes in terms of sequence and secondary structure. However, the nucleotides at the cleavage site have shown that mutational effects on two twister sister-like (TS-like) ribozymes are different from the twister sister ribozyme. The discovery of TS-like ribozymes reveals a ribozyme class with the simplest and, perhaps, the most primitive structure needed for self-cleavage.

Advances in Genotyping Detection of Fragmented Nucleic Acids

Single nucleotide variant (SNV) detection is pivotal in various fields, including disease diagnosis, viral screening, genetically modified organism (GMO) identification, and genotyping. However, detecting SNVs presents significant challenges due to the fragmentation of nucleic acids caused by cellular apoptosis, molecular shearing, and physical degradation processes such as heating. Fragmented nucleic acids often exhibit variable lengths and inconsistent breakpoints, complicating the accurate detection of SNVs. This article delves into the underlying causes of nucleic acid fragmentation and synthesizes the strengths and limitations of next-generation sequencing technology, high-resolution melting curves, molecular probes, and CRISPR-based approaches for SNV detection in fragmented nucleic acids. By providing a detailed comparative analysis, it seeks to offer valuable insights for researchers working to overcome the challenges of SNV detection in fragmented samples, ultimately advancing the accurate and efficient detection of single nucleotide variants across diverse applications.

Fighting the battle against evolution: designing genetically modified organisms for evolutionary stability

Synthetic biology has made significant progress in many areas, but a major challenge that has received limited attention is the evolutionary stability of synthetic constructs made of heterologous genes. The expression of these constructs in microorganisms, that is, production of proteins that are not necessary for the organism, is a metabolic burden, leading to a decrease in relative fitness and make the synthetic constructs unstable over time. This is a significant concern for the synthetic biology community, particularly when it comes to bringing this technology out of the laboratory. In this review, we discuss the issue of evolutionary stability in synthetic biology and review the available tools to address this challenge.

Evidence of vertical transmission of fowl adenovirus 8b in ducks

Fowl adenovirus mainly causes hydropericardium hepatitis syndrome (HHS), inclusion body hepatitis (IBH) and gizzard erosion (GE), etc. In 2015, the first outbreak of HHS was reported in broiler chickens in central China, followed by an outbreak in waterfowl. The first outbreak of HHS in broiler flocks in central China in 2015, followed by outbreaks in waterfowl, has severely restricted the healthy development of the poultry industry. During the investigation, fowl adenovirus was detected in ducklings from a total of seven hatcheries in Shandong, Inner Mongolia and Jiangsu provinces. In addition, the DNA of fowl adenovirus was detected in breeding ducks and their progeny. To test the hypothesis that FAdV can be transmitted vertically, sixty 250-day-old Cherry Valley breeder ducks were divided equally into three groups for experimental infection. FAdV-8b SDLY isolate (duck/Shandong/SDLY/2021, SDLY) preserved in our laboratory was injected intramuscularly into group A and inoculated orally into group B. FAdV-8b DNA was detected in the yolk membranes, embryos and allantoic fluid of duck embryos in the FAdV-infected group after inoculation. In addition, the FAdV-8b hexon gene isolated from yolk membranes, embryos, allantoic fluid and duck eggs was close to 100% nucleotide homology to the FAdV-8b hexon gene isolated from laying duck ovaries, indicating that fowl adenovirus can be transmitted vertically in ducks. These findings provide evidence for the possible vertical transmission of fowl adenovirus from breeder ducks to ducklings.

Pre-PCR Mutation-Enrichment Methods for Liquid Biopsy Applications

Liquid biopsy is having a remarkable impact on healthcare- and disease-management in the context of personalized medicine. Circulating free DNA (cfDNA) is one of the most instructive liquid-biopsy-based biomarkers and harbors valuable information for diagnostic, predictive, and prognostic purposes. When it comes to cancer, circulating DNA from the tumor (ctDNA) has a wide range of applications, from early cancer detection to the early detection of relapse or drug resistance, and the tracking of the dynamic genomic make-up of tumor cells. However, the detection of ctDNA remains technically challenging, due, in part, to the low frequency of ctDNA among excessive circulating cfDNA originating from normal tissues. During the past three decades, mutation-enrichment methods have emerged to boost sensitivity and enable facile detection of low-level mutations. Although most developed techniques apply mutation enrichment during or following initial PCR, there are a few techniques that allow mutation selection prior to PCR, which provides advantages. Pre-PCR enrichment techniques can be directly applied to genomic DNA and diminish the influence of PCR errors that can take place during amplification. Moreover, they have the capability for high multiplexity and can be followed by established mutation detection and enrichment technologies without changes to their established procedures. The first approaches for pre-PCR enrichment were developed by employing restriction endonucleases directly on genomic DNA in the early 1990s. However, newly developed pre-PCR enrichment methods provide higher sensitivity and versatility. This review describes the available pre-PCR enrichment methods and focuses on the most recently developed techniques (NaME-PrO, UVME, and DEASH/MAESTRO), emphasizing their applications in liquid biopsies.

Hairpin Structure Facilitates Multiplex High-Fidelity DNA Amplification in Real-Time Polymerase Chain Reaction

Clinically and biologically, it is essential to detect rare DNA-sequence variants for early cancer diagnosis or drug-resistance mutation identification. Some of the common quantitative polymerase chain reaction (qPCR)-based variant detection methods are restricted in the limit of detection (LoD) because the DNA polymerases used for these methods have a high polymerase misincorporation rate; thus, the detection sensitivity is sometimes unsatisfactory. With the proofreading activity, high-fidelity (HiFi) DNA polymerases have a 50- to 250-fold higher fidelity. However, there are currently no proper probe-based designs functioning as the fluorescence indicator allowing multiplexed HiFi qPCR reactions, thus restricting the application of HiFi DNA polymerases like the variant detection. We presented the occlusion system, composed of a 5'-overhanged primer with a fluorophore modification and a probe with a short-stem hairpin and a 3' quencher modification. We demonstrated that the occlusion system allowed multiplexing HiFi qPCR reaction, and it was compatible with the current variant-enrichment method to improve the LoD up to 10-fold. Thus, the occlusion system satisfactorily functioned as an efficient fluorescence indicator in HiFi qPCR reactions and allowed the application of HiFi DNA polymerases in variant detection methods to improve detection sensitivity.

The Integration of Gold Nanoparticles with Polymerase Chain Reaction for Constructing Colorimetric Sensing Platforms for Detection of Health-Related DNA and Proteins

Polymerase chain reaction (PCR) is the standard tool in genetic information analysis, and the desirable detection merits of PCR have been extended to disease-related protein analysis. Recently, the combination of PCR and gold nanoparticles (AuNPs) to construct colorimetric sensing platforms has received considerable attention due to its high sensitivity, visual detection, capability for on-site detection, and low cost. However, it lacks a related review to summarize and discuss the advances in this area. This perspective gives an overview of established methods based on the combination of PCR and AuNPs for the visual detection of health-related DNA and proteins. Moreover, this work also addresses the future trends and perspectives for PCR-AuNP hybrid biosensors.

On the relation between input and output distributions of scRNA-seq experiments

MOTIVATION

Single-cell RNA sequencing determines RNA copy numbers per cell for a given gene. However, technical noise poses the question how observed distributions (output) are connected to their cellular distributions (input).

RESULTS

We model a single-cell RNA sequencing setup consisting of PCR amplification and sequencing, and derive probability distribution functions for the output distribution given an input distribution. We provide copy number distributions arising from single transcripts during PCR amplification with exact expressions for mean and variance. We prove that the coefficient of variation of the output of sequencing is always larger than that of the input distribution. Experimental data reveals the variance and mean of the input distribution to obey characteristic relations, which we specifically determine for a HeLa dataset. We can calculate as many moments of the input distribution as are known of the output distribution (up to all). This, in principle, completely determines the input from the output distribution.

AVAILABILITY AND IMPLEMENTATION

Source code freely available at https://github.com/danielschw188/InputOutputSCRNASeq.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Error Characterization and Statistical Modeling Improves Circulating Tumor DNA Detection by Droplet Digital PCR

BACKGROUND

Droplet digital PCR (ddPCR) is a widely used and sensitive application for circulating tumor DNA (ctDNA) detection. As ctDNA is often found in low abundance, methods to separate low-signal readouts from noise are necessary. We aimed to characterize the ddPCR-generated noise and, informed by this, create a sensitive and specific ctDNA caller.

METHODS

We built 2 novel complimentary ctDNA calling methods: dynamic limit of blank and concentration and assay-specific tumor load estimator (CASTLE). Both methods are informed by empirically established assay-specific noise profiles. Here, we characterized noise for 70 mutation-detecting ddPCR assays by applying each assay to 95 nonmutated samples. Using these profiles, the performance of the 2 new methods was assessed in a total of 9447 negative/positive reference samples and in 1311 real-life plasma samples from colorectal cancer patients. Lastly, performances were compared to 7 literature-established calling methods.

RESULTS

For many assays, noise increased proportionally with the DNA input amount. Assays targeting transition base changes were more error-prone than transversion-targeting assays. Both our calling methods successfully accounted for the additional noise in transition assays and showed consistently high performance regardless of DNA input amount. Calling methods that were not noise-informed performed less well than noise-informed methods. CASTLE was the only calling method providing a statistical estimate of the noise-corrected mutation level and call certainty.

CONCLUSIONS

Accurate error modeling is necessary for sensitive and specific ctDNA detection by ddPCR. Accounting for DNA input amounts ensures specific detection regardless of the sample-specific DNA concentration. Our results demonstrate CASTLE as a powerful tool for ctDNA calling using ddPCR.

Critical assessment of structure-based approaches to improve protein resistance in aqueous ionic liquids by enzyme-wide saturation mutagenesis

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for green industrial processes involving biocatalysts, but often reduce enzyme activity. Experimental and computational methods are applied to predict favorable substitution sites and, most often, subsequent site-directed surface charge modifications are introduced to enhance enzyme resistance towards aIL. However, almost no studies evaluate the prediction precision with random mutagenesis or the application of simple data-driven filtering processes. Here, we systematically and rigorously evaluated the performance of 22 previously described structure-based approaches to increase enzyme resistance to aIL based on an experimental complete site-saturation mutagenesis library of Bacillus subtilis Lipase A (BsLipA) screened against four aIL. We show that, surprisingly, most of the approaches yield low gain-in-precision (GiP) values, particularly for predicting relevant positions: 14 approaches perform worse than random mutagenesis. Encouragingly, exploiting experimental information on the thermostability of BsLipA or structural weak spots of BsLipA predicted by rigidity theory yields GiP = 3.03 and 2.39 for relevant variants and GiP = 1.61 and 1.41 for relevant positions. Combining five simple-to-compute physicochemical and evolutionary properties substantially increases the precision of predicting relevant variants and positions, yielding GiP = 3.35 and 1.29. Finally, combining these properties with predictions of structural weak spots identified by rigidity theory additionally improves GiP for relevant variants up to 4-fold to ∼10 and sustains or increases GiP for relevant positions, resulting in a prediction precision of ∼90% compared to ∼9% in random mutagenesis. This combination should be applicable to other enzyme systems for guiding protein engineering approaches towards improved aIL resistance.

Phylogeographic and SNPs Analyses of Bemisia tabaci B Mitotype Populations Reveal Only Two of Eight Haplotypes Are Invasive

Simple Summary

The whitefly Bemisia tabaci taxon consists of an undefined number of morphologically identical genetic variants of which only a few, including the B, harbor invasive haplotypes. These haplotypes have potential to upsurge and become important pests and plant virus vectors in irrigated agroecosystems worldwide. In the 1980s, unprecedented outbreaks associated with the B variant were reported worldwide, however, the precise origin(s) of the invasive haplotypes has not been determined. In this study, available B. tabaci mitochondrial gene sequences were examined for patterns of conserved single nucleotide changes (SNPs). The whitefly sequence records represented North Africa-Middle Eastern habitats, the proposed B variant center of origin, and distant locales recently invaded by haplotype(s) of the B variant. Unexpectedly, the analysis revealed eight SNPs groups (haplotypes) demonstrating that the genetic architecture of the B mitoype is more complex than previously recognized. Also, the distribution patterns of the eight B haplotypes were tightly linked to well-defined eco-geographic regions, suggesting the different groups have diversified by geographic isolation. Contrary to claims that collectively, the B variant is invasive, only two of the eight haplotypic groups have established in geographical locations outside of their zone of endemism.

Abstract

The Bemisia tabaci cryptic species contains 39 known mitotypes of which the B and Q are best recognized for having established outside their extant endemic range. In the 1980s, previously uncharacterized haplotype(s) of the B mitotype rapidly established in tropical and subtropical locales distant from their presumed center of origin, leading to displacement of several native mitotypes and extreme damage to crops and other vegetation particularly in irrigated agroecosystems. To trace the natural and evolutionary history of the invasive B haplotypes, a phylo-biogeographic study was undertaken. Patterns of single nucleotide polymorphisms (SNPs) and signatures potentially indicative of geographic isolation were investigated using a globally representative mitochondrial cytochrome oxidase I gene (mtCOI) sequence database. Eight haplotype groups within the North Africa-Middle East (NAFME) region were differentiated, NAFME 1–8. The NAFME 1–3 haplotypes were members of the same population that is associated with warm desert climate niches of the Arabian Peninsula and east coastal Africa-Ethiopia. The NAFME 4 and 5 haplotypes are endemic to warm and cold semi-arid niches delimited by the Irano-Turanian floristic region, itself harboring extensive biodiversity. Haplotypes 6 and 7 co-occurred in the Middle East along eastern Mediterranean Sea landmasses, while NAFME 8 was found to be endemic to Cyprus, Turkey, and desert micro-niches throughout Egypt and Israel. Contrary to claims that collectively, the B mitotype is invasive, NAFME 6 and 8 are the only haplotypes to have established in geographical locations outside of their zone of endemism.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

The Diversity of a Polyclonal FluCell-SELEX Library Outperforms Individual Aptamers as Emerging Diagnostic Tools for the Identification of Carbapenem Resistant Pseudomonas aeruginosa

Textbook procedures require the use of individual aptamers enriched in SELEX libraries which are subsequently chemically synthesized after their biochemical characterization. Here we show that this reduction of the available sequence space of large libraries and thus the diversity of binding molecules reduces the labelling efficiency and fidelity of selected single aptamers towards different strains of the human pathogen Pseudomonas aeruginosa compared to a polyclonal aptamer library enriched by a whole-cell-SELEX involving fluorescent aptamers. The library outperformed single aptamers in reliable and specific targeting of different clinically relevant strains, allowed to inhibit virulence associated cellular functions and identification of bound cell surface targets by aptamer based affinity purification and mass spectrometry. The stunning ease of this FluCell-SELEX and the convincing performance of the P. aeruginosa specific library may pave the way towards generally new and efficient diagnostic techniques based on polyclonal aptamer libraries not only in clinical microbiology.

Evidence of vertical transmission of novel astrovirus virus in goose

In 2019, a novel goose astrovirus (GoAstV) epidemiological investigation on geese was conducted in Shandong province, China. During the investigation, a high prevalence of novel GoAstV was observed in symptom-free breeding geese flocks. Moreover, the novel GoAstV-specific RNA was detected in either breeder birds or their progenies. To verify the hypothesis that the novel GoAstV could be transmitted vertically, a total of 42 WuLong breeder geese, aged 335 days, were equally divided into three groups for experimental infection. The SDPY isolate of novel GoAstV (A/goose/Shandong/SDPY/2018, SDPY), preserved in our laboratory, was injected intramuscularly to subjects of group A while orally inoculated to those of group B. After the inoculation, novel GoAstV RNA was detected in vitelline membrane, embryos, and allantoic fluid of goose embryos in novel GoAstV infected groups. Moreover, the ORF2 gene of novel GoAstV from vitelline membrane, embryo, allantoic fluid as well as gosling shared almost 100 % nucleotide homology to a novel GoAstV virus isolated from the goose ovary which produced the egg, suggesting that the novel GoAstV can be vertically transmitted in the goose. Taken together, the findings provide evidence of possible vertical transmission of novel GoAstV from breeding goose to goslings.

Accurate inference of the full base-pairing structure of RNA by deep mutational scanning and covariation-induced deviation of activity

Despite the large number of noncoding RNAs in human genome and their roles in many diseases include cancer, we know very little about them due to lack of structural clues. The centerpiece of the structural clues is the full RNA base-pairing structure of secondary and tertiary contacts that can be precisely obtained only from costly and time-consuming 3D structure determination. Here, we performed deep mutational scanning of self-cleaving CPEB3 ribozyme by error-prone PCR and showed that a library of <5 × 104 single-to-triple mutants is sufficient to infer 25 of 26 base pairs including non-nested, nonhelical, and noncanonical base pairs with both sensitivity and precision at 96%. Such accurate inference was further confirmed by a twister ribozyme at 100% precision with only noncanonical base pairs as false negatives. The performance was resulted from analyzing covariation-induced deviation of activity by utilizing both functional and nonfunctional variants for unsupervised classification, followed by Monte Carlo (MC) simulated annealing with mutation-derived scores. Highly accurate inference can also be obtained by combining MC with evolution/direct coupling analysis, R-scape or epistasis analysis. The results highlight the usefulness of deep mutational scanning for high-accuracy structural inference of self-cleaving ribozymes with implications for other structured RNAs that permit high-throughput functional selections.

Determination of Mutational Spectra Induced by Environmental Toxicants in Complex Human Cell Populations

Polycyclic aromatic hydrocarbons (PAHs) are widely distributed in the environment and have potent mutagenic and carcinogenic activities. Studies of mutations induced by these compounds in human cells can help acquire an understanding of their mutagenic pathways. In this chapter, independent cultures of a human cell line expressing cytochrome P450 CYP1A1 (cell line MCL-5) were treated with benzo(a)pyrene (BaP) or dibenzo(a,l)pyrene (DBP), and mutants at the hypoxanthine phosphoribosyltransferase (HPRT) locus were selected en masse by 6-thioguanine resistance (6TG^R). The kinds and positions of the mutations occurring in the third exon of the HPRT gene were analyzed in the mixed HPRT^R mutant cell populations using a combination of polymerase chain reaction (PCR) and denaturing gradient gel electrophoresis (DGGE). Mutant bands were excised from the gel, amplified using PCR, and sequenced to determine the kinds and positions, or spectrum of mutations.

Fidelity of DNA polymerases in the detection of intraindividual variation of mitochondrial DNA

Here we investigated the consequences of PCR amplification errors in the identification of intraindividual mtDNA variation. The bumblebee Bombus morio was chosen as model for the COI gene amplification tests with two DNA polymerases (Taq and Q5) presenting different error rates. The amplifications using Taq resulted in a significant increase of singleton haplotypes per individual in comparison to Q5. The sequence characteristics indicated that Taq resulted haplotypes are mostly due to amplification errors. Studies focusing on intraindividual variability should address special attention to the DNA polymerase fidelity to avoid overestimation of heteroplasmic haplotypes.

Patterns of PCR Amplification Artifacts of the Fungal Barcode Marker in a Hybrid Mushroom

The polymerase chain reaction (PCR) is widely used in modern biology and medicine. However, PCR artifacts can complicate the interpretation of PCR-based results. The internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster is the consensus fungal barcode marker and suspected PCR artifacts have been reported in many studies, especially for the analyses of environmental fungal samples. At present, the patterns of PCR artifacts in the whole fungal ITS region (ITS1+5.8S+ITS2) are not known. In this study, we analyzed the error rates of PCR at three template complexity levels using the divergent copies of ITS from the mushroom Agaricus subrufescens. Our results showed that PCR using the Phusion^® High-Fidelity DNA Polymerase has a per nucleotide error rate of about 4 × 10^–6 per replication. Among the detected mutations, transitions were much more frequent than transversions, insertions, and deletions. When divergent alleles were mixed as templates in the same reaction, a significant proportion (∼30%) of recombinant molecules were detected. The in vitro mixed-template results were comparable to those obtained from using the genomic DNA of the original mushroom specimen as template. Our results indicate that caution should be in place when interpreting ITS sequences from individual fungal specimens, especially those containing divergent ITS copies. Similar results could also happen to PCR-based analyses of other multicopy DNA fragments as well as single-copy DNA sequences with divergent alleles in diploid organisms.

High-resolution mapping of DNA polymerase fidelity using nucleotide imbalances and next-generation sequencing

DNA polymerase fidelity is affected by both intrinsic properties and environmental conditions. Current strategies for measuring DNA polymerase error rate in vitro are constrained by low error subtype sensitivity, poor scalability, and lack of flexibility in types of sequence contexts that can be tested. We have developed the Magnification via Nucleotide Imbalance Fidelity (MagNIFi) assay, a scalable next-generation sequencing assay that uses a biased deoxynucleotide pool to quantitatively shift error rates into a range where errors are frequent and hence measurement is robust, while still allowing for accurate mapping to error rates under typical conditions. This assay is compatible with a wide range of fidelity-modulating conditions, and enables high-throughput analysis of sequence context effects on base substitution and single nucleotide deletion fidelity using a built-in template library. We validate this assay by comparing to previously established fidelity metrics, and use it to investigate neighboring sequence-mediated effects on fidelity for several DNA polymerases. Through these demonstrations, we establish the MagNIFi assay for robust, high-throughput analysis of DNA polymerase fidelity.

Loop assembly: a simple and open system for recursive fabrication of DNA circuits

High-efficiency methods for DNA assembly have enabled the routine assembly of synthetic DNAs of increased size and complexity. However, these techniques require customization, elaborate vector sets or serial manipulations for the different stages of assembly. We have developed Loop assembly based on a recursive approach to DNA fabrication. The system makes use of two Type IIS restriction endonucleases and corresponding vector sets for efficient and parallel assembly of large DNA circuits. Standardized level 0 parts can be assembled into circuits containing 1, 4, 16 or more genes by looping between the two vector sets. The vectors also contain modular sites for hybrid assembly using sequence overlap methods. Loop assembly enables efficient and versatile DNA fabrication for plant transformation. We show the construction of plasmids up to 16 genes and 38 kb with high efficiency (> 80%). We have characterized Loop assembly on over 200 different DNA constructs and validated the fidelity of the method by high-throughput Illumina plasmid sequencing. Our method provides a simple generalized solution for DNA construction with standardized parts. The cloning system is provided under an OpenMTA license for unrestricted sharing and open access.

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the "Mediterranean" (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the "African silver-leafing" (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the "ASL" does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies.

Detecting Somatic Mutations in Normal Cells

Somatic mutations have been studied extensively in the context of cancer. Recent studies have demonstrated that high-throughput sequencing data can be used to detect somatic mutations in non-tumor cells. Analysis of such mutations allows us to better understand the mutational processes in normal cells, explore cell lineages in development, and examine potential associations with age-related disease. We describe here approaches for characterizing somatic mutations in normal and non-tumor disease tissues. We discuss several experimental designs and common pitfalls in somatic mutation detection, as well as more recent developments such as phasing and linked-read technology. With the dramatically increasing numbers of samples undergoing genome sequencing, bioinformatic analysis will enable the characterization of somatic mutations and their impact on non-cancer tissues.

Sequence-based nomenclature: a reply to Thines et al. and Zamora et al. and provisions for an amended proposal "from the floor" to allow DNA sequences as types of names

We reply to two recently published, multi-authored opinion papers by opponents of sequence-based nomenclature, namely Zamora et al. (IMA Fungus9: 167-175,2018) and Thines et al. (IMA Fungus9: 177-183, 2018). While we agree with some of the principal arguments brought forward by these authors, we address misconceptions and demonstrate that some of the presumed evidence presented in these papers has been wrongly interpreted. We disagree that allowing sequences as types would fundamentally alter the nature of types, since a similar nature of abstracted features as type is already allowed in the Code (Art. 40.5), namely an illustration. We also disagree that there is a high risk of introducing artifactual taxa, as this risk can be quantified at well below 5 %, considering the various types of high-throughput sequencing errors. Contrary to apparently widespread misconceptions, sequence-based nomenclature cannot be based on similarity-derived OTUs and their consensus sequences, but must be derived from rigorous, multiple alignment-based phylogenetic methods and quantitative, single-marker species recognition algorithms, using original sequence reads; it is therefore identical in its approach to single-marker studies based on physical types, an approach allowed by the Code. We recognize the limitations of the ITS as a single fungal barcoding marker, but point out that these result in a conservative approach, with "false negatives" surpassing "false positives"; a desirable feature of sequence-based nomenclature. Sequence-based nomenclature does not aim at accurately resolving species, but at naming sequences that represent unknown fungal lineages so that these can serve as a means of communication, so ending the untenable situation of an exponentially growing number of unlabeled fungal sequences that fill online repositories. The risks are outweighed by the gains obtained by a reference library of named sequences spanning the full array of fungal diversity. Finally, we elaborate provisions in addition to our original proposal to amend the Code that would take care of the issues brought forward by opponents to this approach. In particular, taking up the idea of the Candidatus status of invalid, provisional names in prokaryote nomenclature, we propose a compromise that would allow valid publication of voucherless, sequence-based names in a consistent manner, but with the obligate designation as "nom. seq." (nomen sequentiae). Such names would not have priority over specimen- or culture-based names unless either epitypified with a physical type or adopted for protection on the recommendation of a committee of the International Commission on the Taxonomy of Fungi following evaluation based on strict quality control of the underlying studies based on established rules or recommendations.

Evolution of a Thermophilic Strand-Displacing Polymerase Using High-Temperature Isothermal Compartmentalized Self-Replication

Strand-displacing polymerases are a crucial component of isothermal amplification (IA) reactions, where the lack of thermal cycling reduces equipment needs and improves the time to answer, especially for point-of-care applications. In order to improve the function of strand-displacing polymerases, we have developed an emulsion-based directed evolution scheme, high-temperature isothermal compartmentalized self-replication (HTI-CSR) that does not rely on thermal cycling. Starting from an algorithm-optimized shuffled library of exonuclease-deficient Family A polymerases from Geobacillus stearothermophilus (Bst LF) and Thermus aquaticus (Klentaq), we have applied HTI-CSR to generate a more thermostable strand-displacing polymerase variant that performs well in loop-mediated isothermal amplification and rolling circle amplification, even after thermal challenges of up to 95 °C that lead to better primer annealing. The new enzyme (v5.9) is also capable of a variety of new reactions, including isothermal hyperbranched rolling circle amplification. The HTI-CSR method should now prove useful for evolving additional beneficial phenotypes in strand-displacing polymerases.

Retrieval of a million high-quality, full-length microbial 16S and 18S rRNA gene sequences without primer bias

Small subunit ribosomal RNA (SSU rRNA) genes, 16S in bacteria and 18S in eukaryotes, have been the standard phylogenetic markers used to characterize microbial diversity and evolution for decades. However, the reference databases of full-length SSU rRNA gene sequences are skewed to well-studied ecosystems and subject to primer bias and chimerism, which results in an incomplete view of the diversity present in a sample. We combine poly(A)-tailing and reverse transcription of SSU rRNA molecules with synthetic long-read sequencing to generate high-quality, full-length SSU rRNA sequences, without primer bias, at high throughput. We apply our approach to samples from seven different ecosystems and obtain more than a million SSU rRNA sequences from all domains of life, with an estimated raw error rate of 0.17%. We observe a large proportion of novel diversity, including several deeply branching phylum-level lineages putatively related to the Asgard Archaea. Our approach will enable expansion of the SSU rRNA reference databases by orders of magnitude, and contribute to a comprehensive census of the tree of life.

Locked Nucleic Acid Technology for Highly Sensitive Detection of Somatic Mutations in Cancer

The molecular diagnosis of the cancer mutational status is essential for modern clinical laboratory medicine. Mutations in EGFR, KRAS, BRAF, and PIK3CA genes are widely analyzed in solid tumors such as lung cancer, colorectal cancer, breast cancer, and melanoma. The allele-specific polymerase chain reaction, high-resolution melting, and Sanger sequencing are used for detecting and identifying gene mutations in many clinical laboratories. The locked nucleic acid (LNA) is a class of nucleic acid analogs that contain a methylene bridge connecting the 2' oxygen and 4' carbon in the ribose moiety. This methylene bridge locks the ribose group into a C3'-endo conformation. LNA, including an oligonucleotide, increases the thermal stability of hybrid strands. The use of LNA technology in molecular diagnostic methods improves the specificity and sensitivity of assays. This review describes routinely analyzed mutations and molecular diagnostic methods used in the clinical laboratory along with the performance improvement of mutational analysis with LNA.

A systematic comparison of error correction enzymes by next-generation sequencing

Gene synthesis, the process of assembling gene-length fragments from shorter groups of oligonucleotides (oligos), is becoming an increasingly important tool in molecular and synthetic biology. The length, quality and cost of gene synthesis are limited by errors produced during oligo synthesis and subsequent assembly. Enzymatic error correction methods are cost-effective means to ameliorate errors in gene synthesis. Previous analyses of these methods relied on cloning and Sanger sequencing to evaluate their efficiencies, limiting quantitative assessment. Here, we develop a method to quantify errors in synthetic DNA by next-generation sequencing. We analyzed errors in model gene assemblies and systematically compared six different error correction enzymes across 11 conditions. We find that ErrASE and T7 Endonuclease I are the most effective at decreasing average error rates (up to 5.8-fold relative to the input), whereas MutS is the best for increasing the number of perfect assemblies (up to 25.2-fold). We are able to quantify differential specificities such as ErrASE preferentially corrects C/G transversions whereas T7 Endonuclease I preferentially corrects A/T transversions. More generally, this experimental and computational pipeline is a fast, scalable and extensible way to analyze errors in gene assemblies, to profile error correction methods, and to benchmark DNA synthesis methods.

Spiked Genes: A Method to Introduce Random Point Nucleotide Mutations Evenly throughout an Entire Gene Using a Complete Set of Spiked Oligonucleotides for the Assembly

In vitro mutagenesis methods have revolutionized biological research and the biotechnology industry. In this study, we describe a mutagenesis method based on synthesizing a gene using a complete set of forward and reverse spiked oligonucleotides that have been modified to introduce a low ratio of mutant nucleotides at each position. This novel mutagenesis scheme named "Spiked Genes" yields a library of clones with an enhanced mutation distribution due to its unbiased nucleotide incorporation. Using the far-red fluorescent protein emKate as a model, we demonstrated that Spiked Genes yields richer libraries than those obtained via enzymatic methods. We obtained a library without bias toward any nucleotide or base pair and with even mutations, transitions, and transversion frequencies. Compared with enzymatic methods, the proposed synthetic approach for the creation of gene libraries represents an improved strategy for screening protein variants and does not require a starting template.

Limitations and challenges of genetic barcode quantification

Genetic barcodes are increasingly used to track individual cells and to quantitatively assess their clonal contributions over time. Although barcode quantification relies entirely on counting sequencing reads, detailed studies about the method’s accuracy are still limited. We report on a systematic investigation of the relation between barcode abundance and resulting read counts after amplification and sequencing using cell-mixtures that contain barcodes with known frequencies (“miniBulks”). We evaluated the influence of protocol modifications to identify potential sources of error and elucidate possible limitations of the quantification approach. Based on these findings we designed an advanced barcode construct (BC32) to improved barcode calling and quantification, and to ensure a sensitive detection of even highly diluted barcodes. Our results emphasize the importance of using curated barcode libraries to obtain interpretable quantitative data and underline the need for rigorous analyses of any utilized barcode library in terms of reliability and reproducibility.

Mapping DNA polymerase errors by single-molecule sequencing

Genomic integrity is compromised by DNA polymerase replication errors, which occur in a sequence-dependent manner across the genome. Accurate and complete quantification of a DNA polymerase's error spectrum is challenging because errors are rare and difficult to detect. We report a high-throughput sequencing assay to map in vitro DNA replication errors at the single-molecule level. Unlike previous methods, our assay is able to rapidly detect a large number of polymerase errors at base resolution over any template substrate without quantification bias. To overcome the high error rate of high-throughput sequencing, our assay uses a barcoding strategy in which each replication product is tagged with a unique nucleotide sequence before amplification. This allows multiple sequencing reads of the same product to be compared so that sequencing errors can be found and removed. We demonstrate the ability of our assay to characterize the average error rate, error hotspots and lesion bypass fidelity of several DNA polymerases.

Assembly of long DNA sequences using a new synthetic Escherichia coli-yeast shuttle vector

Synthetic biology is a newly developed field of research focused on designing and rebuilding novel biomolecular components, circuits, and networks. Synthetic biology can also help understand biological principles and engineer complex artificial metabolic systems. DNA manipulation on a large genome-wide scale is an inevitable challenge, but a necessary tool for synthetic biology. To improve the methods used for the synthesis of long DNA fragments, here we constructed a novel shuttle vector named pGF (plasmid Genome Fast) for DNA assembly in vivo. The BAC plasmid pCC1BAC, which can accommodate large DNA molecules, was chosen as the backbone. The sequence of the yeast artificial chromosome (YAC) regulatory element CEN6-ARS4 was synthesized and inserted into the plasmid to enable it to replicate in yeast. The selection sequence HIS3, obtained by polymerase chain reaction (PCR) from the plasmid pBS313, was inserted for screening. This new synthetic shuttle vector can mediate the transformation-associated recombination (TAR) assembly of large DNA fragments in yeast, and the assembled products can be transformed into Escherichia coli for further amplification. We also conducted in vivo DNA assembly using pGF and yeast homologous recombination and constructed a 31-kb long DNA sequence from the cyanophage PP genome. Our findings show that this novel shuttle vector would be a useful tool for efficient genome-scale DNA reconstruction.

Polymerase specific error rates and profiles identified by single molecule sequencing

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

Rapid Determination of HIV-1 Mutant Frequencies and Mutation Spectra Using an mCherry/EGFP Dual-Reporter Viral Vector

The high mutation rate of human immunodeficiency virus type-1 (HIV-1) has been a pivotal factor in its evolutionary success as a human pathogen, driving the emergence of drug resistance, immune system escape, and invasion of distinct anatomical compartments. Extensive research has focused on understanding how various cellular and viral factors alter the rates and types of mutations produced during viral replication. Here, we describe a single-cycle dual-reporter vector assay that relies upon the detection of mutations that eliminate either expression of mCherry or enhanced green fluorescent protein (EGFP). The reporter-based method can be used to efficiently quantify changes in mutant frequencies and mutation spectra that arise due to a variety of factors, including viral mutagens, drug resistance mutations, cellular physiology, and APOBEC3 proteins.

Conservation of nucleotide sequences for molecular diagnosis of Middle East respiratory syndrome coronavirus, 2015

•

Protocols for the laboratory diagnosis of Middle East respiratory syndrome coronavirus (MERS-CoV) were developed in 2012.

•

Nucleotide sequences of the viruses detected between 2012 and 2015 were analyzed.

•

Mismatches between the sequences of primers/probes and viruses were found.

•

Such mismatches could lead to a lower sensitivity of the assay.

Infection due to the Middle East respiratory syndrome coronavirus (MERS-CoV) is widespread. The present study was performed to assess the protocols used for the molecular diagnosis of MERS-CoV by analyzing the nucleotide sequences of viruses detected between 2012 and 2015, including sequences from the large outbreak in eastern Asia in 2015. Although the diagnostic protocols were established only 2 years ago, mismatches between the sequences of primers/probes and viruses were found for several of the assays. Such mismatches could lead to a lower sensitivity of the assay, thereby leading to false-negative diagnosis. A slight modification in the primer design is suggested. Protocols for the molecular diagnosis of viral infections should be reviewed regularly after they are established, particularly for viruses that pose a great threat to public health such as MERS-CoV.

Sources of PCR-induced distortions in high-throughput sequencing data sets

PCR permits the exponential and sequence-specific amplification of DNA, even from minute starting quantities. PCR is a fundamental step in preparing DNA samples for high-throughput sequencing. However, there are errors associated with PCR-mediated amplification. Here we examine the effects of four important sources of error—bias, stochasticity, template switches and polymerase errors—on sequence representation in low-input next-generation sequencing libraries. We designed a pool of diverse PCR amplicons with a defined structure, and then used Illumina sequencing to search for signatures of each process. We further developed quantitative models for each process, and compared predictions of these models to our experimental data. We find that PCR stochasticity is the major force skewing sequence representation after amplification of a pool of unique DNA amplicons. Polymerase errors become very common in later cycles of PCR but have little impact on the overall sequence distribution as they are confined to small copy numbers. PCR template switches are rare and confined to low copy numbers. Our results provide a theoretical basis for removing distortions from high-throughput sequencing data. In addition, our findings on PCR stochasticity will have particular relevance to quantification of results from single cell sequencing, in which sequences are represented by only one or a few molecules.

Extensive Mitochondrial Heteroplasmy in Natural Populations of a Resurging Human Pest, the Bed Bug (Hemiptera: Cimicidae)

Homoplasmy, the occurrence of a single mitochondrial DNA haplotype within an individual, has been the accepted condition across most organisms in the animal kingdom. In recent years, a number of exceptions to this rule have been reported, largely due to the ease with which single nucleotide polymorphisms can be detected. Evidence of heteroplasmy-two or more mitochondrial variants within a single individual-has now been documented in a number of invertebrates; however, when present, heteroplasmy usually occurs at low frequencies both within individuals and within populations. The implications of heteroplasmy may be far reaching, both to the individual in relation to its health and fitness, and when considering the evolutionary dynamics of populations. We present novel evidence for frequent mtDNA heteroplasmy in the bed bug, Cimex lectularius L. (Hemiptera: Cimicidae). Our findings show that heteroplasmy is common, with 5 of 29 (17%) populations screened exhibiting two mitochondrial variants in a ∼1:2 ratio within each individual. We hypothesize that the mechanism underlying heteroplasmy in bed bugs is paternal leakage because some haplotypes were shared among unrelated populations and no evidence for nuclear mitochondrial DNA sequences was detected.

Targeted ultradeep next-generation sequencing as a method for KIT D816V mutation analysis in mastocytosis

Next-generation sequencing (NGS) is becoming increasingly used for diagnostic mutation analysis in myeloid neoplasms and may also represent a feasible technique in mastocytosis. However, detection of the KIT D816V mutation requires a highly sensitive method in most patients due to the typically low mutation levels. In this study, we established an NGS-based KIT mutation analysis and analyzed the sensitivity of D816V detection using the Ion Torrent platform. Eighty-two individual NGS analyses were included in the study. All samples were also analyzed using highly sensitive KIT D816V mutation-specific qPCR. Measurements of the background level in D816V-negative samples supported a cutoff for positivity of 0.2% in three different NGS panels. Clinical samples from patients with SM that tested positive using qPCR with a D816V allele burden >0.2% also tested positive using NGS. Samples that tested positive using qPCR with an allele burden <0.2% tested negative using NGS. We thereby demonstrate that caution should be taken when using the potentially very sensitive NGS technique for KIT D816V mutation analysis in mastocytosis, as many patients with SM have D816V mutation levels below the detection limit of NGS. A dedicated and highly sensitive KIT D816V mutation analysis therefore remains important in mastocytosis diagnostics.

Investigation into the annotation of protocol sequencing steps in the sequence read archive

Background

The workflow for the production of high-throughput sequencing data from nucleic acid samples is complex. There are a series of protocol steps to be followed in the preparation of samples for next-generation sequencing. The quantification of bias in a number of protocol steps, namely DNA fractionation, blunting, phosphorylation, adapter ligation and library enrichment, remains to be determined.

Results

We examined the experimental metadata of the public repository Sequence Read Archive (SRA) in order to ascertain the level of annotation of important sequencing steps in submissions to the database. Using SQL relational database queries (using the SRAdb SQLite database generated by the Bioconductor consortium) to search for keywords commonly occurring in key preparatory protocol steps partitioned over studies, we found that 7.10%, 5.84% and 7.57% of all records (fragmentation, ligation and enrichment, respectively), had at least one keyword corresponding to one of the three protocol steps. Only 4.06% of all records, partitioned over studies, had keywords for all three steps in the protocol (5.58% of all SRA records).

Conclusions

The current level of annotation in the SRA inhibits systematic studies of bias due to these protocol steps. Downstream from this, meta-analyses and comparative studies based on these data will have a source of bias that cannot be quantified at present.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0064-7) contains supplementary material, which is available to authorized users.

Enhancing the speed of morpholino-DNA biosensor by electrokinetic concentration of DNA in a microfluidic chip

Electrokinetic methods that conveniently concentrate charged analytes by orders of magnitude are highly attractive for nucleic acid assays where they can bypass the complexity and costs of enzyme-based amplification. The present study demonstrates an electrokinetic concentration device incorporating charge-neutral morpholino (MO) probes: as DNA analyte is concentrated in a microfluidic channel using ion concentration polarization (ICP) it is simultaneously hybridized to spots of complementary MO probes immobilized on the channel floor. This approach is uniquely favored by the match between the optimum buffer ionic strength of approximately 10mM for both MO-DNA surface hybridization and electrokinetic concentration. The simple and easily scalable poly(dimethylsiloxane) (PDMS) microfluidic device was fabricated using soft lithography and contact printing of a conductive polymer, poly(3,4-ethylenedioxythiophene)-polystyrene sulfonate (

PEDOT

PSS) as a cation-selective membrane material. Using the microfluidic concentrator, we could increase the concentration of DNA by three orders of magnitude in less than 5 min at an electric field of 75 Vcm(-1). The 1000-fold increase in concentration of DNA led to an increase in the speed of MO-DNA hybridization by two orders of magnitude and enabled a detection sensitivity of ~1 nM within 15 min of concentration. Using the proposed microfluidic concentrator, we also demonstrated a rapid hybridization with a binary DNA mixture, containing a fully complementary and a non-complementary sequence to mimic molecular backgrounds present in real DNA samples.

Multi-template polymerase chain reaction

PCR is a formidable and potent technology that serves as an indispensable tool in a wide range of biological disciplines. However, due to the ease of use and often lack of rigorous standards many PCR applications can lead to highly variable, inaccurate, and ultimately meaningless results. Thus, rigorous method validation must precede its broad adoption to any new application. Multi-template samples possess particular features, which make their PCR analysis prone to artifacts and biases: multiple homologous templates present in copy numbers that vary within several orders of magnitude. Such conditions are a breeding ground for chimeras and heteroduplexes. Differences in template amplification efficiencies and template competition for reaction compounds undermine correct preservation of the original template ratio. In addition, the presence of inhibitors aggravates all of the above-mentioned problems. Inhibitors might also have ambivalent effects on the different templates within the same sample. Yet, no standard approaches exist for monitoring inhibitory effects in multitemplate PCR, which is crucial for establishing compatibility between samples.

The spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β

To provide molecular-level insights into the spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β (polβ), we report four crystal structures of polβ complexed with dG•dTTP and dA•dCTP mismatches in the presence of Mg2+ or Mn2+. The Mg(2+)-bound ground-state structures show that the dA•dCTP-Mg2+ complex adopts an 'intermediate' protein conformation while the dG•dTTP-Mg2+ complex adopts an open protein conformation. The Mn(2+)-bound 'pre-chemistry-state' structures show that the dA•dCTP-Mn2+ complex is structurally very similar to the dA•dCTP-Mg2+ complex, whereas the dG•dTTP-Mn2+ complex undergoes a large-scale conformational change to adopt a Watson-Crick-like dG•dTTP base pair and a closed protein conformation. These structural differences, together with our molecular dynamics simulation studies, suggest that polβ increases replication fidelity via a two-stage mismatch discrimination mechanism, where one is in the ground state and the other in the closed conformation state. In the closed conformation state, polβ appears to allow only a Watson-Crick-like conformation for purine•pyrimidine base pairs, thereby discriminating the mismatched base pairs based on their ability to form the Watson-Crick-like conformation. Overall, the present studies provide new insights into the spontaneous replication error and the replication fidelity mechanisms of polβ.

Determining lower limits of detection of digital PCR assays for cancer-related gene mutations

•

Digital PCR offers very high assay sensitivity and limit of detection.

•

An approach for calculating limit of detection is demonstrated for two EGFR assays.

•

Assay LoDs have been evaluated for eighteen cancer targets.

Digital PCR offers very high sensitivity compared to many other technologies for processing molecular detection assays. Herein, a process is outlined for determining the lower limit of detection (LoD) of two droplet-based digital PCR assays for point mutations of the epidermal growth factor receptor (EGFR) gene. Hydrolysis probe mutation-detection assays for EGFR p.L858R and p.T790M mutations were characterized in detail. Furthermore, sixteen additional cancer-related mutation assays were explored by the same approach. For the EGFR L8585R assay, the assay sensitivity is extremely good, and thus, the LoD is limited by the amount of amplifiable DNA that is analyzed. With 95% confidence limits, the LoD is one mutant in 180,000 wild-type molecules for the evaluation of 3.3 μg of genomic DNA, and detection of one mutant molecule in over 4 million wild-type molecules was achieved when 70 million copies of DNA were processed. The measured false-positive rate for the EGFR L8585R assay is one in 14 million, which indicates the theoretical LoD if an unlimited amount of DNA is evaluated. For the EFGR T790M assay, the LoD is one mutant in 13,000 for analysis of a 3.3 μg sample of genomic DNA, and the dPCR assay limit sensitivity approaches one mutant in 22,000 wild-type molecules.

Error Rate Comparison during Polymerase Chain Reaction by DNA Polymerase

As larger-scale cloning projects become more prevalent, there is an increasing need for comparisons among high fidelity DNA polymerases used for PCR amplification. All polymerases marketed for PCR applications are tested for fidelity properties (i.e., error rate determination) by vendors, and numerous literature reports have addressed PCR enzyme fidelity. Nonetheless, it is often difficult to make direct comparisons among different enzymes due to numerous methodological and analytical differences from study to study. We have measured the error rates for 6 DNA polymerases commonly used in PCR applications, including 3 polymerases typically used for cloning applications requiring high fidelity. Error rate measurement values reported here were obtained by direct sequencing of cloned PCR products. The strategy employed here allows interrogation of error rate across a very large DNA sequence space, since 94 unique DNA targets were used as templates for PCR cloning. The six enzymes included in the study, Taq polymerase, AccuPrime-Taq High Fidelity, KOD Hot Start, cloned Pfu polymerase, Phusion Hot Start, and Pwo polymerase, we find the lowest error rates with Pfu, Phusion, and Pwo polymerases. Error rates are comparable for these 3 enzymes and are >10x lower than the error rate observed with Taq polymerase. Mutation spectra are reported, with the 3 high fidelity enzymes displaying broadly similar types of mutations. For these enzymes, transition mutations predominate, with little bias observed for type of transition.

Improved multiple displacement amplification (iMDA) and ultraclean reagents

BACKGROUND

Next-generation sequencing sample preparation requires nanogram to microgram quantities of DNA; however, many relevant samples are comprised of only a few cells. Genomic analysis of these samples requires a whole genome amplification method that is unbiased and free of exogenous DNA contamination. To address these challenges we have developed protocols for the production of DNA-free consumables including reagents and have improved upon multiple displacement amplification (iMDA).

RESULTS

A specialized ethylene oxide treatment was developed that renders free DNA and DNA present within Gram positive bacterial cells undetectable by qPCR. To reduce DNA contamination in amplification reagents, a combination of ion exchange chromatography, filtration, and lot testing protocols were developed. Our multiple displacement amplification protocol employs a second strand-displacing DNA polymerase, improved buffers, improved reaction conditions and DNA free reagents. The iMDA protocol, when used in combination with DNA-free laboratory consumables and reagents, significantly improved efficiency and accuracy of amplification and sequencing of specimens with moderate to low levels of DNA. The sensitivity and specificity of sequencing of amplified DNA prepared using iMDA was compared to that of DNA obtained with two commercial whole genome amplification kits using 10 fg (~1-2 bacterial cells worth) of bacterial genomic DNA as a template. Analysis showed >99% of the iMDA reads mapped to the template organism whereas only 0.02% of the reads from the commercial kits mapped to the template. To assess the ability of iMDA to achieve balanced genomic coverage, a non-stochastic amount of bacterial genomic DNA (1 pg) was amplified and sequenced, and data obtained were compared to sequencing data obtained directly from genomic DNA. The iMDA DNA and genomic DNA sequencing had comparable coverage 99.98% of the reference genome at ≥1X coverage and 99.9% at ≥5X coverage while maintaining both balance and representation of the genome.

CONCLUSIONS

The iMDA protocol in combination with DNA-free laboratory consumables, significantly improved the ability to sequence specimens with low levels of DNA. iMDA has broad utility in metagenomics, diagnostics, ancient DNA analysis, pre-implantation embryo screening, single-cell genomics, whole genome sequencing of unculturable organisms, and forensic applications for both human and microbial targets.

Heat-stable enzymes from extremely thermophilic and hyperthermophilic microorganisms

Only in the last decade have microorganisms been discovered which grow near or above 100°C. The enzymes that are formed by these extremely thermophilic (growth temperature 65 to 85°C) and hyperthermophilic (growth temperature 85 to 110°C) microorganisms are of great interest. This review covers the extracellular and intracellular enzymes of these exotic microorganisms that have recently been described. Polymer-hydrolysing enzymes, such as amylolytic, cellulolytic, hemicellulolytic and proteolytic enzymes, will be discussed. In addition, the properties of the intracellular enzymes involved in carbohydrate and amino-acid metabolism and DNA-binding and chaperones and chaperone-like proteins from hyperthermophiles are described. Due to the unusual properties of these heat-stable enzymes, they are expected to fill the gap between biological and chemical processes.

Uncovering Wolbachia diversity upon artificial host transfer

The common endosymbiotic Wolbachia bacteria influence arthropod hosts in multiple ways. They are mostly recognized for their manipulations of host reproduction, yet, more recent studies demonstrate that Wolbachia also impact host behavior, metabolic pathways and immunity. Besides their biological and evolutionary roles, Wolbachia are new potential biological control agents for pest and vector management. Importantly, Wolbachia-based control strategies require controlled symbiont transfer between host species and predictable outcomes of novel Wolbachia-host associations. Theoretically, this artificial horizontal transfer could inflict genetic changes within transferred Wolbachia populations. This could be facilitated through de novo mutations in the novel recipient host or changes of haplotype frequencies of polymorphic Wolbachia populations when transferred from donor to recipient hosts. Here we show that Wolbachia resident in the European cherry fruit fly, Rhagoletis cerasi, exhibit ancestral and cryptic sequence polymorphism in three symbiont genes, which are exposed upon microinjection into the new hosts Drosophila simulans and Ceratitis capitata. Our analyses of Wolbachia in microinjected D. simulans over 150 generations after microinjection uncovered infections with multiple Wolbachia strains in trans-infected lines that had previously been typed as single infections. This confirms the persistence of low-titer Wolbachia strains in microinjection experiments that had previously escaped standard detection techniques. Our study demonstrates that infections by multiple Wolbachia strains can shift in prevalence after artificial host transfer driven by either stochastic or selective processes. Trans-infection of Wolbachia can claim fitness costs in new hosts and we speculate that these costs may have driven the shifts of Wolbachia strains that we saw in our model system.