Error Rate Comparison during Polymerase Chain Reaction by DNA Polymerase

Data recovery methods for DNA storage based on fountain codes

Today's digital data storage systems typically offer advanced data recovery solutions to address the problem of catastrophic data loss, such as software-based disk sector analysis or physical-level data retrieval methods for conventional hard disk drives. However, DNA-based data storage currently relies solely on the inherent error correction properties of the methods used to encode digital data into strands of DNA. Any error that cannot be corrected utilizing the redundancy added by DNA encoding methods results in permanent data loss. To provide data recovery for DNA storage systems, we present a method to automatically reconstruct corrupted or missing data stored in DNA using fountain codes. Our method exploits the relationships between packets encoded with fountain codes to identify and rectify corrupted or lost data. Furthermore, we present file type-specific and content-based data recovery methods for three file types, illustrating how a fusion of fountain encoding-specific redundancy and knowledge about the data can effectively recover information in a corrupted DNA storage system, both in an automatic and in a guided manual manner. To demonstrate our approach, we introduce DR4DNA, a software toolkit that contains all methods presented. We evaluate DR4DNA using both in-silico and in-vitro experiments.

Deletion of exons 45 to 55 in the DMD gene: from the therapeutic perspective to the in vitro model

BACKGROUND

Gene editing therapies in development for correcting out-of-frame DMD mutations in Duchenne muscular dystrophy aim to replicate benign spontaneous deletions. Deletion of 45-55 DMD exons (del45-55) was described in asymptomatic subjects, but recently serious skeletal and cardiac complications have been reported. Uncovering why a single mutation like del45-55 is able to induce diverse phenotypes and grades of severity may impact the strategies of emerging therapies. Cellular models are essential for this purpose, but their availability is compromised by scarce muscle biopsies.

METHODS

We introduced, as a proof-of-concept, using CRISPR-Cas9 edition, a del45-55 mimicking the intronic breakpoints harboured by a subset of patients of this form of dystrophinopathy (designing specific gRNAs), into a Duchenne patient's cell line. The edited cell line was characterized evaluating the dystrophin expression and the myogenic status.

RESULTS

Dystrophin expression was restored, and the myogenic defects were ameliorated in the edited myoblasts harbouring a specific del45-55. Besides confirming the potential of CRISPR-Cas9 to create tailored mutations (despite the low cleavage efficiency of our gRNAs) as a useful approach to generate in vitro models, we also generated an immortalized myoblast line derived from a patient with a specific del45-55.

CONCLUSIONS

Overall, we provide helpful resources to deepen into unknown factors responsible for DMD-pathophysiology.

An efficient approach for overproduction of DNA polymerase from Pyrococcus furiosus using an optimized autoinduction system in Escherichia coli

High fidelity DNA polymerase from Pyrococcus furiosus (Pfupol) is an attractive alternative to the highly popular DNA polymerase from Thermus aquaticus. Because this enzyme is in great demand for biotechnological applications, optimizing Pfupol production is essential to supplying the industry's expanding demand. T7-induced promoter expression in Escherichia coli expression systems is used to express recombinant Pfupol; however, this method is not cost-effective. Here, we have effectively developed an optimized process for the autoinduction approach of Pfupol expression in a defined medium. To better examine Pfupol's activities, its purified fraction was used. A 71 mg/L of pure Pfupol was effectively produced, resulting in a 2.6-fold increase in protein yield when glucose, glycerol, and lactose were added in a defined medium at concentrations of 0.05%, 1%, and 0.6%, respectively, and the condition for production in a 5 L bioreactor was as follow: 200 rpm, 3 vvm, and 10% inoculant. Furthermore, the protein exhibited 1445 U/mg of specific activity when synthesized in its active state. This work presents a high level of Pfupol production, which makes it an economically viable and practically useful approach.

Chemical unclonable functions based on operable random DNA pools

Physical unclonable functions (PUFs) based on unique tokens generated by random manufacturing processes have been proposed as an alternative to mathematical one-way algorithms. However, these tokens are not distributable, which is a disadvantage for decentralized applications. Finding unclonable, yet distributable functions would help bridge this gap and expand the applications of object-bound cryptography. Here we show that large random DNA pools with a segmented structure of alternating constant and randomly generated portions are able to calculate distinct outputs from millions of inputs in a specific and reproducible manner, in analogy to physical unclonable functions. Our experimental data with pools comprising up to >1010 unique sequences and encompassing >750 comparisons of resulting outputs demonstrate that the proposed chemical unclonable function (CUF) system is robust, distributable, and scalable. Based on this proof of concept, CUF-based anti-counterfeiting systems, non-fungible objects and decentralized multi-user authentication are conceivable.

Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma

BACKGROUND

Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory.

METHODS

Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples.

RESULTS

dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response.

CONCLUSIONS

Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

Storage-D: A user-friendly platform that enables practical and personalized DNA data storage

Deoxyribonucleic acid (DNA) has been suggested as a very promising medium for data storage in recent years. Although numerous studies have advocated for DNA data storage, its practical application remains obscure and there is a lack of a user-oriented platform. Here, we developed a DNA data storage platform, named Storage-D, which allows users to convert their data into DNA sequences of any length and vice versa by selecting algorithms, error-correction, random-access, and codec pin strategies in terms of their own choice. It incorporates a newly designed "Wukong" algorithm, which provides over 20 trillion codec pins for data privacy use. This algorithm can also control GC content to the selected standard, as well as adjust the homopolymer run length to a defined level, while maintaining a high coding potential of ~1.98 bis/nt, allowing it to outperform previous algorithms. By connecting to a commercial DNA synthesis and sequencing platform with "Storage-D," we successfully stored "Diagnosis and treatment protocol for COVID-19 patients" into 200 nt oligo pools in vitro, and 500 bp genes in vivo which replicated in both normal and extreme bacteria. Together, this platform allows for practical and personalized DNA data storage, potentially with a wide range of applications.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Extracting, filtering and simulating cellular barcodes using CellBarcode tools

Identifying true DNA cellular barcodes among polymerase chain reaction and sequencing errors is challenging. Current tools are restricted in the diversity of barcode types supported or the analysis strategies implemented. As such, there is a need for more versatile and efficient tools for barcode extraction, as well as for tools to investigate which factors impact barcode detection and which filtering strategies to best apply. Here we introduce the package CellBarcode and its barcode simulation kit, CellBarcodeSim, that allows efficient and versatile barcode extraction and filtering for a range of barcode types from bulk or single-cell sequencing data using a variety of filtering strategies. Using the barcode simulation kit and biological data, we explore the technical and biological factors influencing barcode identification and provide a decision tree on how to optimize barcode identification for different barcode settings. We believe that CellBarcode and CellBarcodeSim have the capability to enhance the reproducibility and interpretation of barcode results across studies.

Site-Directed Mutagenesis to Mutate Multiple Residues in a Single Reaction

Site-directed mutagenesis (SDM) is a technique that allows mutation of specific nucleotide(s) in a codon to study its functional implications in a protein. Commercial kits are available, which require high-performance liquid chromatography purified oligos for this purpose. These kits are expensive, and they are not very efficient, so one has to sequence several clones to get a desired one. We present here a simple method that requires only crude oligos, commercially available high-fidelity enzymes, and the success rate is close to 100%. In addition, up to 6 different mutations can be introduced in one reaction without causing any fortuitous change in the vector backbone. Using this strategy, we have introduced 32 S/T➔A substitutions in the N-terminus head and 13 changes in the C-terminus tail domain of vimentin.

Codon optimization of a gene encoding DNA polymerase from Pyrococcus furiosus and its expression in Escherichia coli

BACKGROUND

DNA polymerase is an essential component in PCR assay for DNA synthesis. Improving DNA polymerase with characteristics indispensable for a powerful assay is crucial because it can be used in wide-range applications. Derived from Pyrococcus furiosus, Pfu DNA polymerase (Pfu pol) is one of the excellent polymerases due to its high fidelity. Therefore, we aimed to develop Pfu pol from a synthetic gene with codon optimization to increase its protein yield in Escherichia coli.

RESULTS

Recombinant Pfu pol was successfully expressed and purified with a two-step purification process using nickel affinity chromatography, followed by anion exchange chromatography. Subsequently, the purified Pfu pol was confirmed by Western blot analysis, resulting in a molecular weight of approximately 90 kDa. In the final purification process, we successfully obtained a large amount of purified enzyme (26.8 mg/L). Furthermore, the purified Pfu pol showed its functionality and efficiency when tested for DNA amplification using the standard PCR.

CONCLUSIONS

Overall, a high-level expression of recombinant Pfu pol was achieved by employing our approach in the present study. In the future, our findings will be useful for studies on synthesizing recombinant DNA polymerase in E. coli expression system.

A digital twin for DNA data storage based on comprehensive quantification of errors and biases

Archiving data in synthetic DNA offers unprecedented storage density and longevity. Handling and storage introduce errors and biases into DNA-based storage systems, necessitating the use of Error Correction Coding (ECC) which comes at the cost of added redundancy. However, insufficient data on these errors and biases, as well as a lack of modeling tools, limit data-driven ECC development and experimental design. In this study, we present a comprehensive characterisation of the error sources and biases present in the most common DNA data storage workflows, including commercial DNA synthesis, PCR, decay by accelerated aging, and sequencing-by-synthesis. Using the data from 40 sequencing experiments, we build a digital twin of the DNA data storage process, capable of simulating state-of-the-art workflows and reproducing their experimental results. We showcase the digital twin's ability to replace experiments and rationalize the design of redundancy in two case studies, highlighting opportunities for tangible cost savings and data-driven ECC development.

Manipulative-based Activity Using Pop Beads for Demonstration of Sanger Sequencing

Sanger sequencing, also known as dideoxy sequencing, is a widely used method for DNA sequencing, particularly for cloned plasmids and clinical samples. This technique requires a combination of essential biochemistry skills, such as a chain-termination reaction, gel electrophoresis, and fluorescence detection. Unfortunately, there is a lack of activities that replicate the Sanger sequencing process for students to learn and practice these skills. To address this issue, a manipulative-based Sanger sequencing activity was developed that incorporates colorful pop beads to demonstrate a chain-termination reaction, separation of products, and fluorescence detection. The beads represent deoxynucleotides and dideoxynucleotides, allowing for a visual representation of DNA fragment generation. This kinesthetic learning activity offers a high visual impact for students, aiding in their understanding of the Sanger sequencing process, and can also be used to illustrate polymerase chain reaction (PCR)-based techniques.

Prevalence of Human Norovirus GII.4 Sydney 2012 [P31] between 2019 and 2021 among Young Children from Rural Communities in South Africa

Acute gastroenteritis (AGE) accounts for considerable morbidity and mortality in the paediatric population worldwide, especially in low-income countries. Human norovirus (HNoV), particularly GII.4 strains, are important agents of AGE. This study aimed to detect and characterise HNoV in children with and without AGE. Between 2019 and 2021, 300 stool samples (200 AGE and 100 without AGE) were collected from children below 5 years of age referred to the healthcare facilities of the rural communities of Vhembe District, South Africa. After detection using real-time RT-PCR, HNoV positive samples were subjected to RT-PCR and Sanger sequencing. Partial nucleotide sequences (capsid/RdRp) were aligned using the Muscle tool, and phylogenetic analysis was performed using MEGA 11. The nucleotides' percent identity among HNoV strains was compared using ClustalW software. A significant difference in HNoV prevalence between AGE children (37%; 74/200) and non-AGE (14%; 14/100) was confirmed (p < 0.0001). Genogroup II (GII) HNoV was predominant in AGE children (80%; 59/74), whereas most non-AGE children were infected by the GI norovirus genogroup (64%; 9/14). GII.4 Sydney 2012 [P31] strains were dominant (59%; 19/32) during the study period. A phylogenetic analysis revealed a close relationship between the HNoV strains identified in this study and those circulating worldwide; however, ClustalW showed less than 50% nucleotide similarity between strains from this study and those from previously reported norovirus studies in the same region. Our findings indicate significant changes over time in the circulation of HNoV strains, as well as the association between high HNoV prevalence and AGE symptoms within the study area. The monitoring of HuNoV epidemiology, along with stringent preventive measures to mitigate the viral spread and the burden of AGE, are warranted.

DNA Polymerases for Whole Genome Amplification: Considerations and Future Directions

In the same way that specialized DNA polymerases (DNAPs) replicate cellular and viral genomes, only a handful of dedicated proteins from various natural origins as well as engineered versions are appropriate for competent exponential amplification of whole genomes and metagenomes (WGA). Different applications have led to the development of diverse protocols, based on various DNAPs. Isothermal WGA is currently widely used due to the high performance of Φ29 DNA polymerase, but PCR-based methods are also available and can provide competent amplification of certain samples. Replication fidelity and processivity must be considered when selecting a suitable enzyme for WGA. However, other properties, such as thermostability, capacity to couple replication, and double helix unwinding, or the ability to maintain DNA replication opposite to damaged bases, are also very relevant for some applications. In this review, we provide an overview of the different properties of DNAPs widely used in WGA and discuss their limitations and future research directions.

Enhanced Biomass Production of Recombinant Pfu DNA Polymerase Producer Escherichia coli BL21(DE3) by Optimization of Induction Variables Using Response Surface Methodology

Pfu DNA polymerase is one of the most preferred molecular enzymes that is isolated from the hyperthermophilic Pyrococcus furiosus and used for high-throughput DNA synthesis by the polymerase chain reaction. Therefore, an efficient Pfu DNA polymerase production method is necessary for molecular techniques. In the present study, Pfu DNA polymerase was expressed in recombinant Escherichia coli BL21(DE3) and significant parameters for the biomass production were optimized using the central composite design which is the most popular method of response surface methodology. Induction conditions including cell density prior induction (OD_600nm), post-induction temperature, IPTG concentration, and post-induction time and their interactions on biomass production were investigated. The maximum biomass production (14.1 g/L) in shake flasks was achieved using the following predicted optimal conditions: OD_600nm before induction of 0.4 and the induction at 32 °C for 7.7 h, with 0.6 mM IPTG. Optimized culture conditions were implemented to scale up experiments. 22% and 70% increase in biomass production was achieved in 3 L and 10 L bioreactors, respectively as compared to initial biomass production observed in unoptimized conditions. Similary, a 30% increase of Pfu DNA polymerase production was obtained after the optimization. The polymerase activity of the purifed Pfu DNA polymerase was assessed by PCR amplification and determined as 2.9 U/μl by comparison with commercial Pfu DNA polymerase. The findings of this study indicated that the proposed fermentation conditions will contribute to further scale‑up studies to enhance the biomass for the production of other recombinant proteins.

Mobile and Self-Sustained Data Storage in an Extremophile Genomic DNA

DNA has been pursued as a novel biomaterial for digital data storage. While large-scale data storage and random access have been achieved in DNA oligonucleotide pools, repeated data accessing requires constant data replenishment, and these implementations are confined in professional facilities. Here, a mobile data storage system in the genome of the extremophile Halomonas bluephagenesis, which enables dual-mode storage, dynamic data maintenance, rapid readout, and robust recovery. The system relies on two key components: A versatile genetic toolbox for the integration of 10-100 kb scale synthetic DNA into H. bluephagenesis genome and an efficient error correction coding scheme targeting noisy nanopore sequencing reads. The storage and repeated retrieval of 5 KB data under non-laboratory conditions are demonstrated. The work highlights the potential of DNA data storage in domestic and field scenarios, and expands its application domain from archival data to frequently accessed data.

Characterisation of Primary Human Hippocampal Astrocyte Cell Culture Following Exposure to Hypoxia

Background

The present study aimed to understand the characterisation of human hippocampal astrocyte following hypoxia exposure. Based on the preliminary screening, 15 min was chosen as the time point and the cells were exposed to different oxygen percentages.

Methods

The Trypan blue viability assay used to examine cell death. Immunofluorescence assay, glial fibrillary acidic protein (GFAP) was used to portray the morphology of astrocytes. The hypoxia-inducible factor 1 (HIF-1) staining was performed to confirm hypoxia induced cell death and there was a dramatic expression of HIF-1α displayed in exposed astrocyte cells compared to the control. In molecular level, genes were chosen, such as glyceraldehyde 3-phosphate dehydrogenase (GAPDH), GFAP, HIF-1α and B-cell lymphoma 2 (Bcl-2) and ran the reverse transcription-polymerase chain reaction (RT-PCR).

Results

Microscope revealed a filamentous and clear nucleus appearance in a control whereas the rupture nuclei with no rigid structure of the cell were found in the 3% oxygen. The control and hypoxia cells were also stained with the annexin V-fluorescein isothiocyanate (annexin V-FITC). Fluorescence microscope reveals astrocyte cells after hypoxia showed higher expression of nuclei but not in control. Merging PI and FITC showed the differences of nuclei expression between the control and hypoxia. In the molecular analysis, there were significant changes of GFAP, HIF-1α and Bcl-2 in hypoxia exposed cells when compared to the control group.

Conclusion

Cells that were exposed to hypoxia (3% oxygen for 15 min) clearly showed damage. General view of human hippocampal astrocyte genomic response to hypoxia was obtained.

Easy method for six-fragment Golden Gate Assembly of modular vectors

Efficient cloning techniques are a requirement for synthetic biology. This study provides a simplified cloning method based on Golden Gate Assembly that can be used for rapid vector construction. The building of multiple expression vectors with customizable modules is achieved in a timely manner with minimal hands-on time by removing unnecessary steps in the workflow. The authors constructed a total of 21 mammalian episomal expression vectors and conducted a fluorescence expression comparison for different regulatory region combinations post-transfection in HEK293T and HEPG2 cells. Screening revealed that using the EF-1α promoter in combination with the bovine growth hormone polyadenylation sequence seemed to perform best in the types of cells tested compared with other variants.

Determining On-Target, Off-Target, and Copy Number Status of Transgenic Events After CRISPR/Cas9 Targeted AAVS1 Safe-Harbor Modification of iPSCs Using Double-Control Quantitative Copy Number PCR

Double-control quantitative copy number PCR (dc-qcnPCR) is a recently described tool that can be used to quantify donor DNA insertions in genetically modified monoclonal cell lines. In conjunction with an insert-confirmation PCR, the technique can quickly and easily identify clones containing on-target heterozygous or homozygous donor DNA integrations and exclude off-target insertions. The genetic manipulation of immortal cell lines is a versatile tool to elucidate cellular signaling pathways and protein functions. Despite recent advances in the precision of genetic engineering tools such as CRISPR/Cas9, transcription-activator-like effector nucleases (TALENs), and zinc-finger nucleases (ZFNs), it is still essential to verify the accurate insertion of the sequence of interest (donor DNA) into the targeted genomic DNA (gDNA) locus. This precise integration into a genetic safe harbor, and exclusion of the donor DNA from functionally relevant genes, can ensure normal cellular functionality. Current methods to analyze the specificity of donor DNA insertions either are cost-prohibitive or create dependency on manufacturers for assay design and production. The dc-qcnPCR method is a simple, yet powerful, approach that can be prepared and carried out in any laboratory equipped with standard molecular biology supplies. Here we provide step-by-step instructions to prepare and perform the dc-qcnPCR, and its companion insert-confirmation PCR, to determine donor DNA insertion numbers in monoclonal cell lines genetically modified through CRISPR/Cas9. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Genetic modification at AAVS1 safe harbor in induced pluripotent stem cells (IMR90-4) using CRISPR/Cas9: from plasmid design to monoclonal expansion Support Protocol 1: Measurement of Gaussia luciferase activity to verify reporter protein functionality Support Protocol 2: Verification of monoclonal expansion using immunofluorescence. Basic Protocol 2: Insert-confirmation PCR Basic Protocol 3: Design and preparation of double-control quantitative copy number PCR reagents and quantification of donor DNA integrations in genetically modified monoclonal cells.

Low-level point heteroplasmy detection in human mitogenomes amplified with different polymerases and sequenced on MiSeq FGx platform

Introduction

Massively parallel sequencing of mitogenomes usually requires prior amplification. The PCR step may influence the quality of the data obtained, especially when low-level heteroplasmy detection is applied.

Aim

The aim of this study was to compare the reliability of two different DNA polymerases in detecting homoplasmic and heteroplasmic substitutions in human mitogenomes.

Material and methods

Mitogenomes of five samples were amplified with Long PCR Enzyme Mix from Fermentas or TaKaRa LA Taq DNA Polymerase from TaKaRa. Then, NexteraTM XT DNA libraries were sequenced on MiSeq FGx platform (Illumina). mtDNA substitutions were called for alternative variants above the 1% level.

Results

All homoplasmic substitutions detected in amplicons generated with polymerases studied here and sequenced on MiSeq FGx system were consistently identified as homoplasmies with alternative sequencing methods. TaKaRa LA Taq DNA Polymerase was found to be less accurate in low-level heteroplasmy detection than Long PCR Enzyme Mix enzyme as more false negative and false positive results were observed for minority variants called above the 1% level. Nevertheless, both PCR systems studied can be successfully used to detect authentic mtDNA substitutions, for which minority variants exceed the 3.61% level assuming at least 10,000x coverage and sequencing Nextera XT DNA libraries on MiSeq FGx machine.

Conclusions

The accuracy and sensitivity of point heteroplasmy detection with the MiSeq FGx instrument varies on polymerase used for mtDNA amplification. Therefore, it is recommended to validate the laboratory protocols used for mtDNA substitution detection prior to their implementation for the forensic or medical genetics purposes.

Cost-Effective PCR-Based Identification of Tunga penetrans (Siphonaptera) Larvae Extracted from Soil Samples Containing PCR Inhibitor-Rich Material

Tungiasis is a neglected tropical disease caused by skin-penetrating female Tunga penetrans fleas. Although tungiasis causes severe health problems, its ecology is poorly understood and morphological descriptions of the larvae are unavailable. To identify T. penetrans immature stages and sites where they develop, diagnostic PCRs are required. However, flea larvae feed on soil organic matter rich in PCR inhibitors. Here, three DNA preparation methods, including a soil DNA kit that removes inhibitors, a simple ammonium acetate precipitation approach (AmAcet) and a crude lysate of larvae (CL), were combined with amplification by the highly processive FIREPol^® Taq or the inhibitor-resistant Phusion^® polymerase. Independent of the polymerase used, the frequency of successful amplification, C_q values and PCR efficacies for the low-cost CL and AmAcet methods were superior to the commercial kit for amplification of a 278 bp partial internal transcribed spacer-2 (ITS-2) and a 730 bp pan-Siphonaptera cytochrome oxidase II PCR. For the CL method combined with Phusion^® polymerase, the costs were approximately 20-fold lower than for the methods based on the soil DNA kit, which is a considerable advantage in resource-poor settings. The ITS-2 PCR did not amplify Ctenocephalides felis genomic or Tunga trimammilata ITS-2 plasmid DNA, meaning it can be used to specifically identify T. penetrans.

The mutational signatures of formalin fixation on the human genome

Clinical archives of patient material near-exclusively consist of formalin-fixed and paraffin-embedded (FFPE) blocks. The ability to precisely characterise mutational signatures from FFPE-derived DNA has tremendous translational potential. However, sequencing of DNA derived from FFPE material is known to be riddled with artefacts. Here we derive genome-wide mutational signatures caused by formalin fixation. We show that the FFPE-signature is highly similar to signature 30 (the signature of Base Excision Repair deficiency due to NTHL1 mutations), and chemical repair of DNA lesions leads to a signature highly similar to signature 1 (clock-like signature due to spontaneous deamination of methylcytosine). We demonstrate that using uncorrected mutational catalogues of FFPE samples leads to major mis-assignment of signature activities. To correct for this, we introduce FFPEsig, a computational algorithm to rectify the formalin-induced artefacts in the mutational catalogue. We demonstrate that FFPEsig enables accurate mutational signature analysis both in simulated and whole-genome sequenced FFPE cancer samples. FFPEsig thus provides an opportunity to unlock additional clinical potential of archival patient tissues.

evSeq: Cost-Effective Amplicon Sequencing of Every Variant in a Protein Library

Widespread availability of protein sequence-fitness data would revolutionize both our biochemical understanding of proteins and our ability to engineer them. Unfortunately, even though thousands of protein variants are generated and evaluated for fitness during a typical protein engineering campaign, most are never sequenced, leaving a wealth of potential sequence-fitness information untapped. Primarily, this is because sequencing is unnecessary for many protein engineering strategies; the added cost and effort of sequencing are thus unjustified. It also results from the fact that, even though many lower-cost sequencing strategies have been developed, they often require at least some access to and experience with sequencing or computational resources, both of which can be barriers to access. Here, we present every variant sequencing (evSeq), a method and collection of tools/standardized components for sequencing a variable region within every variant gene produced during a protein engineering campaign at a cost of cents per variant. evSeq was designed to democratize low-cost sequencing for protein engineers and, indeed, anyone interested in engineering biological systems. Execution of its wet-lab component is simple, requires no sequencing experience to perform, relies only on resources and services typically available to biology labs, and slots neatly into existing protein engineering workflows. Analysis of evSeq data is likewise made simple by its accompanying software (found at github.com/fhalab/evSeq, documentation at fhalab.github.io/evSeq), which can be run on a personal laptop and was designed to be accessible to users with no computational experience. Low-cost and easy-to-use, evSeq makes the collection of extensive protein variant sequence-fitness data practical.

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.

m6A demethylation of cytidine deaminase APOBEC3B mRNA orchestrates arsenic-induced mutagenesis

The cytidine deaminase APOBEC3B (A3B) is an endogenous inducer of somatic mutations and causes chromosomal instability by converting cytosine to uracil in single-stranded DNA. Therefore, identification of factors and mechanisms that mediate A3B expression will be helpful for developing therapeutic approaches to decrease DNA mutagenesis. Arsenic (As) is one well-known mutagen and carcinogen, but the mechanisms by which it induces mutations have not been fully elucidated. Herein, we show that A3B is upregulated and required for As-induced DNA damage and mutagenesis. We found that As treatment causes a decrease of N6-methyladenosine (m6A) modification near the stop codon of A3B, consequently increasing the stability of A3B mRNA. We further reveal that the demethylase FTO is responsible for As-reduced m6A modification of A3B, leading to increased A3B expression and DNA mutation rates in a manner dependent on the m6A reader YTHDF2. Our in vivo data also confirm that A3B is a downstream target of FTO in As-exposed lung tissues. In addition, FTO protein is highly expressed and positively correlates with the protein levels of A3B in tumor samples from human non-small cell lung cancer patients. These findings indicate a previously unrecognized role of A3B in As-triggered somatic mutation and might open new avenues to reduce DNA mutagenesis by targeting the FTO/m6A axis.

Prediction and validation of hematopoietic stem and progenitor cell off-target editing in transplanted rhesus macaques

The programmable nuclease technology CRISPR-Cas9 has revolutionized gene editing in the last decade. Due to the risk of off-target editing, accurate and sensitive methods for off-target characterization are crucial prior to applying CRISPR-Cas9 therapeutically. Here, we utilized a rhesus macaque model to compare the predictive values of CIRCLE-seq, an in vitro off-target prediction method, with in silico prediction (ISP) based solely on genomic sequence comparisons. We use AmpliSeq HD error-corrected sequencing to validate off-target sites predicted by CIRCLE-seq and ISP for a CD33 guide RNA (gRNA) with thousands of off-target sites predicted by ISP and CIRCLE-seq. We found poor correlation between the sites predicted by the two methods. When almost 500 sites predicted by each method were analyzed by error-corrected sequencing of hematopoietic cells following transplantation, 19 off-target sites revealed insertion or deletion mutations. Of these sites, 8 were predicted by both methods, 8 by CIRCLE-seq only, and 3 by ISP only. The levels of cells with these off-target edits exhibited no expansion or abnormal behavior in vivo in animals followed for up to 2 years. In addition, we utilized an unbiased method termed CAST-seq to search for translocations between the on-target site and off-target sites present in animals following transplantation, detecting one specific translocation that persisted in blood cells for at least 1 year following transplantation. In conclusion, neither CIRCLE-seq or ISP predicted all sites, and a combination of careful gRNA design, followed by screening for predicted off-target sites in target cells by multiple methods, may be required for optimizing safety of clinical development.

Detecting MUC1 Variants in Patients Clinicopathologically Diagnosed With Having Autosomal Dominant Tubulointerstitial Kidney Disease

Introduction

Autosomal dominant tubulointerstitial kidney disease (ADTKD)-MUC1 is predominantly caused by frameshift mutations owing to a single-base insertion into the variable number tandem repeat (VNTR) region in MUC1. Because of the complexity of the variant hotspot, identification using short-read sequencers (SRSs) is challenging. Although recent studies have revealed the usefulness of long-read sequencers (LRSs), the prevalence of MUC1 variants in patients with clinically suspected ADTKD remains unknown. We aimed to clarify this prevalence and the genetic characteristics and clinical manifestations of ADTKD-MUC1 in a Japanese population using an SRS and an LRS.

Methods

From January 2015 to December 2019, genetic analysis was performed using an SRS in 48 patients with clinically suspected ADTKD. Additional analyses were conducted using an LRS in patients with negative SRS results.

Results

Short-read sequencing results revealed MUC1 variants in 1 patient harboring a cytosine insertion in the second repeat unit of the VNTR region; however, deeper VNTR regions could not be read by the SRS. Therefore, we conducted long-read sequencing analysis of 39 cases and detected MUC1 VNTR variants in 8 patients (in total, 9 patients from unrelated families). With the inclusion of family-affected patients (n = 31), the median age at the development of end-stage kidney disease (ESKD) was 45 years (95% CI: 40–40 years).

Conclusion

In Japan, the detection rate of MUC1 variants in patients with clinically suspected ADTKD was 18.8%. More than 20% of patients with negative SRS results had MUC1 variants detected by an LRS.

Optimized Loopable Translation as a Platform for the Synthesis of Repetitive Proteins

The expression of long proteins with repetitive amino acid sequences often presents a challenge in recombinant systems. To overcome this obstacle, we report a genetic construct that circularizes mRNA in vivo by rearranging the topology of a group I self-splicing intron from T4 bacteriophage, thereby enabling "loopable" translation. Using a fluorescence-based assay to probe the translational efficiency of circularized mRNAs, we identify several conditions that optimize protein expression from this system. Our data suggested that translation of circularized mRNAs could be limited primarily by the rate of ribosomal initiation; therefore, using a modified error-prone PCR method, we generated a library that concentrated mutations into the initiation region of circularized mRNA and discovered mutants that generated markedly higher expression levels. Combining our rational improvements with those discovered through directed evolution, we report a loopable translator that achieves protein expression levels within 1.5-fold of the levels of standard vectorial translation. In summary, our work demonstrates loopable translation as a promising platform for the creation of large peptide chains, with potential utility in the development of novel protein materials.

MuA-based Molecular Indexing for Rare Mutation Detection by Next-Generation Sequencing

Detection of low-frequency mutations in cancer genomes or other heterogeneous cell populations requires high-fidelity sequencing. Molecular barcoding is one of the key technologies that enables the differentiation of true mutations from errors, which can be caused by sequencing or library preparation processes. However, current approaches where barcodes are introduced via primer extension or adaptor ligation do not utilize the full power of barcoding, due to complicated library preparation workflows and biases. Here we demonstrate the remarkable tolerance of MuA transposase to the presence of multiple replacements in transposon sequence, and explore this unique feature to engineer the MuA transposome complex with randomised nucleotides in 12 transposon positions, which can be introduced as a barcode into the target molecule after transposition event. We applied the approach of Unique MuA-based Molecular Indexing (UMAMI) to assess the power of rare mutation detection by shortgun sequencing on the Illumina platform. Our results show that UMAMI allows detection of rare mutations readily and reliably, and in this paper we report error rate values for the number of thermophilic DNA polymerases measured by using UMAMI.

Characterizing the Relationship Between Neutralization Sensitivity and env Gene Diversity During ART Suppression

Although antiretroviral therapy (ART) successfully suppresses HIV-1 replication, ART-treated individuals must maintain therapy to avoid rebound from an integrated viral reservoir. Strategies to limit or clear this reservoir are urgently needed. Individuals infected for longer periods prior to ART appear to harbor more genetically diverse virus, but the roles of duration of infection and viral diversity in the humoral immune response remain to be studied. We aim to clarify a role, if any, for autologous and heterologous antibodies in multi-pronged approaches to clearing infection. To that end, we have characterized the breadths and potencies of antibody responses in individuals with varying durations of infection and HIV-1 envelope (env) gene diversity as well as the sensitivity of their inducible virus reservoir to broadly neutralizing antibodies (bNAbs). Plasma was collected from 8 well-characterized HIV-1+ males on ART with varied durations of active infection. HIV envs from reservoir-derived outgrowth viruses were amplified and single genome sequenced in order to measure genetic diversity in each participant. IgG from plasma was analyzed for binding titers against gp41 and gp120 proteins, and for neutralizing titers against a global HIV-1 reference panel as well as autologous outgrowth viruses. The sensitivity to bNAbs of these same autologous viruses was measured. Overall, we observed that greater env diversity was associated with higher neutralizing titers against the global panel and also increased resistance to certain bNAbs. Despite the presence of robust anti-HIV-1 antibody titers, we did not observe potent neutralization against autologous viruses. In fact, 3 of 8 participants harbored viruses that were completely resistant to the highest tested concentration of autologous IgG. That this lack of neutralization was observed regardless of ART duration or viral diversity suggests that the inducible reservoir harbors 'escaped' viruses (that co-evolved with autologous antibody responses), rather than proviruses archived from earlier in infection. Finally, we observed that viruses resistant to autologous neutralization remained sensitive to bNAbs, especially CD4bs and MPER bNAbs. Overall, our data suggest that the inducible reservoir is relatively resistant to autologous antibodies and that individuals with limited virus variation in the env gene, such as those who start ART early in infection, are more likely to be sensitive to bNAb treatment.

Lentiviral-Driven Discovery of Cancer Drug Resistance Mutations

Identifying resistance mutations in a drug target provides crucial information. Lentiviral transduction creates multiple types of mutations due to the error-prone nature of the HIV-1 reverse transcriptase (RT). Here we optimized and leveraged this property to identify drug resistance mutations, developing a technique we term LentiMutate. This technique was validated by identifying clinically relevant EGFR resistance mutations, then applied to two additional clinical anticancer drugs: imatinib, a BCR-ABL inhibitor, and AMG 510, a KRAS G12C inhibitor. Novel deletions in BCR-ABL1 conferred resistance to imatinib. In KRAS-G12C or wild-type KRAS, point mutations in the AMG 510 binding pocket or oncogenic non-G12C mutations conferred resistance to AMG 510. LentiMutate should prove highly valuable for clinical and preclinical cancer-drug development. SIGNIFICANCE: LentiMutate can evaluate a drug's on-target activity and can nominate resistance mutations before they occur in patients, which could accelerate and refine drug development to increase the survival of patients with cancer.

The Importance of Quality Control of LSDV Live Attenuated Vaccines for Its Safe Application in the Field

Vaccination is an effective approach to prevent, control and eradicate diseases, including lumpy skin disease (LSD). One of the measures to address farmer hesitation to vaccinate is guaranteeing the quality of vaccine batches. The purpose of this study was to demonstrate the importance of a quality procedure via the evaluation of the LSD vaccine, Lumpivax (Kevevapi). The initial PCR screening revealed the presence of wild type LSD virus (LSDV) and goatpox virus (GTPV), in addition to vaccine LSDV. New phylogenetic PCRs were developed to characterize in detail the genomic content and a vaccination/challenge trial was conducted to evaluate the impact on efficacy and diagnostics. The characterization confirmed the presence of LSDV wild-, vaccine- and GTPV-like sequences in the vaccine vial and also in samples taken from the vaccinated animals. The analysis was also suggestive for the presence of GTPV-LSDV (vaccine/wild) recombinants. In addition, the LSDV status of some of the animal samples was greatly influenced by the differentiating real-PCR used and could result in misinterpretation. Although the vaccine was clinically protective, the viral genomic content of the vaccine (being it multiple Capripox viruses and/or recombinants) and the impact on the diagnostics casts serious doubts of its use in the field.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

An error prone PCR method for small amplicons

Error-prone PCR (epPCR) is a commonly employed approach in molecular biology, especially in directed evolution, to generate libraries of DNA molecules with broad mutational spectrums. Though commonly applied to mutagenize protein coding sequences of several hundreds or thousands of basepairs, we found that commonly used protocols were not suitable for small (<100 bp) amplicons. Here we report a modified error-prone PCR protocol utilizing a Touchdown approach and employing only commercially available components, that should be broadly useful for the researcher interested in concentrating mutations into a small region of plasmid DNA. It will also be useful for achieving very high mutational loads on a standard-sized amplicon.

Establishing a gold standard method for the detection of Cherax reovirus using reverse transcriptase, quantitative, polymerase chain reaction

Cherax reovirus infects redclaw crayfish (Cherax quadricarinatus) and it may be involved in mortalities between 5-20 % and stunting of up to 40 % of survivors. The sequence of the RNA-dependent RNA polymerase was used to develop a reverse transcription, quantitative, PCR (RT-qPCR) which was specific against seven other crustacean viruses (Athtab bunyavirus, Chequa iflavirus, Macrobrachium rosenbergii nodavirus, Gill-associated virus, Taura syndrome virus, White spot syndrome virus, and Penaeus stylirostris Penstylhamaparvovirus) although GAV produced a reaction that was easily separated by melt curve analysis. A strong linear correlation (r² = 0.9965) was obtained between viral quantities ranging from 10⁷ to 10 viral copies/reaction with an amplification efficiency of 0.92. This RT-qPCR is 2-times faster and 100 times more sensitive than a standard RT-PCR using agarose gel electrophoresis with the potential to detect the virus down to 7.64 copies/reaction in clinical samples. In clinical crayfish samples, it was able to detect Cherax reovirus in crayfish when the traditional RT-PCR was negative. Its' measurement of uncertainty was less than 2% (0.02-1.9), similar to PCRs for other crustacean viruses. This RT-qPCR is proposed as the gold standard and should be used for the screening of populations of C. quadricarinatus for broodstock before being used in hatcheries or on farms.

An Automated Correction Algorithm (ALPACA) for ddPCR Data Using Adaptive Limit of Blank and Correction of False Positive Events Improves Specificity of Mutation Detection

BACKGROUND

Bio-Rad droplet-digital PCR is a highly sensitive method that can be used to detect tumor mutations in circulating cell-free DNA (cfDNA) of patients with cancer. Correct interpretation of ddPCR results is important for optimal sensitivity and specificity. Despite its widespread use, no standardized method to interpret ddPCR data is available, nor have technical artifacts affecting ddPCR results been widely studied.

METHODS

False positive rates were determined for 6 ddPCR assays at variable amounts of input DNA, revealing polymerase induced false positive events (PIFs) and other false positives. An in silico correction algorithm, known as the adaptive LoB and PIFs: an automated correction algorithm (ALPACA), was developed to remove PIFs and apply an adaptive limit of blank (LoB) to individual samples. Performance of ALPACA was compared to a standard strategy (no PIF correction and static LoB = 3) using data from commercial reference DNA, healthy volunteer cfDNA, and cfDNA from a real-life cohort of 209 patients with stage IV nonsmall cell lung cancer (NSCLC) whose tumor and cfDNA had been molecularly profiled.

RESULTS

Applying ALPACA reduced false positive results in healthy cfDNA compared to the standard strategy (specificity 98 vs 88%, P = 10-5) and stage IV NSCLC patient cfDNA (99 vs 93%, P = 10-11), while not affecting sensitivity in commercial reference DNA (70 vs 68% P = 0.77) or patient cfDNA (82 vs 88%, P = 0.13). Overall accuracy in patient samples was improved (98 vs 92%, P = 10-7).

CONCLUSIONS

Correction of PIFs and application of an adaptive LoB increases specificity without a loss of sensitivity in ddPCR, leading to a higher accuracy in a real-life cohort of patients with stage IV NSCLC.

Simplified plasmid cloning with a universal MCS design and bacterial in vivo assembly

BACKGROUND

The ability to clone DNA sequences quickly and precisely into plasmids is essential for molecular biology studies. The recent development of seamless cloning technologies has made significant improvements in plasmid construction, but simple and reliable tools are always desirable for time- and labor-saving purposes.

RESULTS

We developed and standardized a plasmid cloning protocol based on a universal MCS (Multiple Cloning Site) design and bacterial in vivo assembly. With this method, the vector is linearized first by PCR (Polymerase Chain Reaction) or restriction digestion. Then a small amount (10 ~ 20 ng) of this linear vector can be mixed with a PCR-amplified insert (5× molar ratio against vector) and transformed directly into competent E. coli cells to obtain the desired clones through in vivo assembly. Since we used a 36-bp universal MCS as the homologous linker, any PCR-amplified insert with ~ 15 bp compatible termini can be cloned into the vector with high fidelity and efficiency. Thus, the need for redesigning insert-amplifying primers according to various vector sequences and the following PCR procedures was eliminated.

CONCLUSIONS

Our protocol significantly reduced hands-on time for preparing transformation reactions, had excellent reliability, and was confirmed to be a rapid and versatile plasmid cloning technique. The protocol contains mostly mixing steps, making it an extremely automation-friendly and promising tool in modern biology studies.

Validation of multiplex PCR sequencing assay of SIV

Background

The generation of accurate and reproducible viral sequence data is necessary to understand the diversity present in populations of RNA viruses isolated from clinical samples. While various sequencing methods are available, they often require high quality templates and high viral titer to ensure reliable data.

Methods

We modified a multiplex PCR and sequencing approach to characterize populations of simian immunodeficiency virus (SIV) isolated from nonhuman primates. We chose this approach with the aim of reducing the number of required input templates while maintaining fidelity and sensitivity. We conducted replicate sequencing experiments using different numbers of quantified viral RNA (vRNA) or viral cDNA as input material. We performed assays with clonal SIVmac239 to detect false positives, and we mixed SIVmac239 and a variant with 24 point mutations (SIVmac239-24X) to measure variant detection sensitivity.

Results

We found that utilizing a starting material of quantified viral cDNA templates had a lower rate of false positives and increased reproducibility when compared to that of quantified vRNA templates. This study identifies the importance of rigorously validating deep sequencing methods and including replicate samples when using a new method to characterize low frequency variants in a population with a small number of templates.

Conclusions

Because the need to generate reproducible and accurate sequencing data from diverse viruses from low titer samples, we modified a multiplex PCR and sequencing approach to characterize SIV from populations from non-human primates. We found that increasing starting template numbers increased the reproducibility and decreased the number of false positives identified, and this was further seen when cDNA was used as a starting material. Ultimately, we highlight the importance of vigorously validating methods to prevent overinterpretation of low frequency variants in a sample.

A SMRT approach for targeted amplicon sequencing of museum specimens (Lepidoptera)-patterns of nucleotide misincorporation

Natural history collections are a valuable resource for molecular taxonomic studies and for examining patterns of evolutionary diversification, particularly in the case of rare or extinct species. However, the recovery of sequence information is often complicated by DNA degradation. This article describes use of the Sequel platform (Pacific Biosciences) to recover the 658 bp barcode region of the mitochondrial cytochrome c oxidase I (COI) gene from 380 butterflies with an average age of 50 years. Nested multiplex PCR was employed for library preparation to facilitate sequence recovery from extracts with low concentrations of highly degraded DNA. By employing circular consensus sequencing (CCS) of short amplicons (circa 150 bp), full-length barcodes could be assembled without a reference sequence, an important advance from earlier protocols which required reference sequences to guide contig assembly. The Sequel protocol recovered COI sequences (499 bp on average) from 318 of 380 specimens (84%), much higher than for Sanger sequencing (26%). Because each read derives from a single molecule, it was also possible to quantify the incidence of substitutions arising from DNA damage. In agreement with past work on sequence changes induced by DNA degradation, the transition C/G → T/A was the most prevalent category of change, but its rate of occurrence (4.58E−4) was so low that it did not impede the recovery of reliable sequences. Because the current protocol recovers COI sequence from most museum specimens, and because sequence fidelity is unaffected by nucleotide misincorporations, large-scale sequence characterization of museum specimens is feasible.

Construction of a highly error-prone DNA polymerase for developing organelle mutation systems

A novel family of DNA polymerases replicates organelle genomes in a wide distribution of taxa encompassing plants and protozoans. Making error-prone mutator versions of gamma DNA polymerases revolutionised our understanding of animal mitochondrial genomes but similar advances have not been made for the organelle DNA polymerases present in plant mitochondria and chloroplasts. We tested the fidelities of error prone tobacco organelle DNA polymerases using a novel positive selection method involving replication of the phage lambda cI repressor gene. Unlike gamma DNA polymerases, ablation of 3'-5' exonuclease function resulted in a modest 5-8-fold error rate increase. Combining exonuclease deficiency with a polymerisation domain substitution raised the organelle DNA polymerase error rate by 140-fold relative to the wild type enzyme. This high error rate compares favourably with error-rates of mutator versions of animal gamma DNA polymerases. The error prone organelle DNA polymerase introduced mutations at multiple locations ranging from two to seven sites in half of the mutant cI genes studied. Single base substitutions predominated including frequent A:A (template: dNMP) mispairings. High error rate and semi-dominance to the wild type enzyme in vitro make the error prone organelle DNA polymerase suitable for elevating mutation rates in chloroplasts and mitochondria.

Guidelines for Sanger sequencing and molecular assay monitoring

Genetic sequencing, or DNA sequencing, using the Sanger technique has become widely used in the veterinary diagnostic community. This technology plays a role in verification of PCR results and is used to provide the genetic sequence data needed for phylogenetic analysis, epidemiologic studies, and forensic investigations. The Laboratory Technology Committee of the American Association of Veterinary Laboratory Diagnosticians has prepared guidelines for sample preparation, submission to sequencing facilities or instrumentation, quality assessment of nucleic acid sequence data performed, and for generating basic sequencing data and phylogenetic analysis for diagnostic applications. This guidance is aimed at assisting laboratories in providing consistent, high-quality, and reliable sequence data when using Sanger-based genetic sequencing as a component of their laboratory services.

MESA: automated assessment of synthetic DNA fragments and simulation of DNA synthesis, storage, sequencing and PCR errors

Summary

The development of de novo DNA synthesis, polymerase chain reaction (PCR), DNA sequencing and molecular cloning gave researchers unprecedented control over DNA and DNA-mediated processes. To reduce the error probabilities of these techniques, DNA composition has to adhere to method-dependent restrictions. To comply with such restrictions, a synthetic DNA fragment is often adjusted manually or by using custom-made scripts. In this article, we present MESA (Mosla Error Simulator), a web application for the assessment of DNA fragments based on limitations of DNA synthesis, amplification, cloning, sequencing methods and biological restrictions of host organisms. Furthermore, MESA can be used to simulate errors during synthesis, PCR, storage and sequencing processes.

Availability and implementation

MESA is available at mesa.mosla.de, with the source code available at github.com/umr-ds/mesa_dna_sim.

Contact

dominik.heider@uni-marburg.de

Supplementary information

Supplementary data are available at Bioinformatics online.

A Systematic Evaluation of High-Throughput Sequencing Approaches to Identify Low-Frequency Single Nucleotide Variants in Viral Populations

High-throughput sequencing such as those provided by Illumina are an efficient way to understand sequence variation within viral populations. However, challenges exist in distinguishing process-introduced error from biological variance, which significantly impacts our ability to identify sub-consensus single-nucleotide variants (SNVs). Here we have taken a systematic approach to evaluate laboratory and bioinformatic pipelines to accurately identify low-frequency SNVs in viral populations. Artificial DNA and RNA "populations" were created by introducing known SNVs at predetermined frequencies into template nucleic acid before being sequenced on an Illumina MiSeq platform. These were used to assess the effects of abundance and starting input material type, technical replicates, read length and quality, short-read aligner, and percentage frequency thresholds on the ability to accurately call variants. Analyses revealed that the abundance and type of input nucleic acid had the greatest impact on the accuracy of SNV calling as measured by a micro-averaged Matthews correlation coefficient score, with DNA and high RNA inputs (107 copies) allowing for variants to be called at a 0.2% frequency. Reduced input RNA (105 copies) required more technical replicates to maintain accuracy, while low RNA inputs (103 copies) suffered from consensus-level errors. Base errors identified at specific motifs identified in all technical replicates were also identified which can be excluded to further increase SNV calling accuracy. These findings indicate that samples with low RNA inputs should be excluded for SNV calling and reinforce the importance of optimising the technical and bioinformatics steps in pipelines that are used to accurately identify sequence variants.

Optimized high-fidelity 3DPCR to assess potential mitochondrial targeting by activation-induced cytidine deaminase

Activation‐induced cytidine deaminase (AID) initiates somatic hypermutation and class switch recombination of immunoglobulin genes in B cells, whereas off‐targeted AID activity contributes to oncogenic mutations and chromosomal translocations associated with B cell malignancies. Paradoxically, only a minority of AID is allowed to access the nuclear genome, but the majority of AID is retained in the cytoplasm. It is unknown whether cytoplasmic AID can access and target the mitochondrial genome [mitochondrial DNA (mtDNA)]. To address this issue, we developed high‐fidelity differential DNA denaturation PCR, which allowed the enrichment of genuine mtDNA mutations and therefore the identification of endogenous mtDNA mutation signatures in vitro. With this approach, we showed that AID targeting to mtDNA is a rare event in AID‐expressing lymphoma lines. Further biochemical and microscopic analysis revealed that a fraction of cytosol AID is associated with the outer membrane of mitochondria but unable to access the mitochondrial matrix. Together, our data suggested that the mitochondrial genome is protected from AID‐mediated mutagenesis by physical segregation of AID from accessing mtDNA within the mitochondrial matrix.

Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding?

True fungi (Fungi) and fungus-like organisms (e.g. Mycetozoa, Oomycota) constitute the second largest group of organisms based on global richness estimates, with around 3 million predicted species. Compared to plants and animals, fungi have simple body plans with often morphologically and ecologically obscure structures. This poses challenges for accurate and precise identifications. Here we provide a conceptual framework for the identification of fungi, encouraging the approach of integrative (polyphasic) taxonomy for species delimitation, i.e. the combination of genealogy (phylogeny), phenotype (including autecology), and reproductive biology (when feasible). This allows objective evaluation of diagnostic characters, either phenotypic or molecular or both. Verification of identifications is crucial but often neglected. Because of clade-specific evolutionary histories, there is currently no single tool for the identification of fungi, although DNA barcoding using the internal transcribed spacer (ITS) remains a first diagnosis, particularly in metabarcoding studies. Secondary DNA barcodes are increasingly implemented for groups where ITS does not provide sufficient precision. Issues of pairwise sequence similarity-based identifications and OTU clustering are discussed, and multiple sequence alignment-based phylogenetic approaches with subsequent verification are recommended as more accurate alternatives. In metabarcoding approaches, the trade-off between speed and accuracy and precision of molecular identifications must be carefully considered. Intragenomic variation of the ITS and other barcoding markers should be properly documented, as phylotype diversity is not necessarily a proxy of species richness. Important strategies to improve molecular identification of fungi are: (1) broadly document intraspecific and intragenomic variation of barcoding markers; (2) substantially expand sequence repositories, focusing on undersampled clades and missing taxa; (3) improve curation of sequence labels in primary repositories and substantially increase the number of sequences based on verified material; (4) link sequence data to digital information of voucher specimens including imagery. In parallel, technological improvements to genome sequencing offer promising alternatives to DNA barcoding in the future. Despite the prevalence of DNA-based fungal taxonomy, phenotype-based approaches remain an important strategy to catalog the global diversity of fungi and establish initial species hypotheses.

Analysis of p67 allelic sequences reveals a subtype of allele type 1 unique to buffalo-derived Theileria parva parasites from southern Africa

East Coast fever (ECF) and Corridor disease (CD) caused by cattle- and buffalo-derived T. parva respectively are the most economically important tick-borne diseases of cattle in the affected African countries. The p67 gene has been evaluated as a recombinant subunit vaccine against ECF, and for discrimination of T. parva parasites causing ECF and Corridor disease. The p67 allele type 1 was first identified in cattle-derived T. parva parasites from East Africa, where parasites possessing this allele type have been associated with ECF. Subsequent characterization of buffalo-derived T. parva parasites from South Africa where ECF was eradicated, revealed the presence of a similar allele type, raising concerns as to whether or not allele type 1 from parasites from the two regions is identical. A 900 bp central fragment of the gene encoding p67 was PCR amplified from T. parva DNA extracted from blood collected from cattle and buffalo in South Africa, Mozambique, Kenya, Tanzania and Uganda, followed by DNA sequence analysis. Four p67 allele types previously described were identified. A subtype of p67 allele type 1 was identified in parasites from clinical cases of CD and buffalo from southern Africa. Notably, p67 allele type 1 sequences from parasites associated with ECF in East Africa and CD in Kenya were identical. Analysis of two p67 B-cell epitopes (TpM12 and AR22.7) revealed amino acid substitutions in allele type 1 from buffalo-derived T. parva parasites from southern Africa. However, both epitopes were conserved in allele type 1 from cattle- and buffalo-derived T. parva parasites from East Africa. These findings reveal detection of a subtype of p67 allele type 1 associated with T. parva parasites transmissible from buffalo to cattle in southern Africa.

Population genetics and comparative mitogenomic analyses reveal cryptic diversity of Amphioctopus neglectus (Cephalopoda: Octopodidae)

This study presented 96 cox1 and 76 cox3 genes of Amphioctopus neglectus populations. Three distinct lineages were formed from phylogenetic trees and networks constructed using haplotypes. Mitogenomes of A. neglectus-a and A. neglectus-b as the representatives of two lineages separated from population genetics were sequenced to compare with A. neglectus at the genome-level. Amphioctopus neglectus-a showed significant differences with A. neglectus, mainly reflected in gene length, intergenic regions and the secondary structure of tandem repeat motifs. Notably, two sequence deletions in mitogenomes of the two representative species were detected in different positions of major non-coding regions, which were the most distinct differences with A. neglectus. Pairwise genetic distances and the phylogenetic analysis supported the relationship of (A. neglectus-a + (A. neglectus + A. neglectus-b)). This study suggested that A. neglectus-a should be considered as a potential cryptic species of this complex, while A. neglectus-b needed further verification to be defined.

CRISPR-sub: Analysis of DNA substitution mutations caused by CRISPR-Cas9 in human cells

CRISPR-Cas9 induces DNA cleavages at desired target sites in a guide RNA-dependent manner; DNA editing occurs through the resulting activity of DNA repair processes including non-homologous end joining (NHEJ), which is dominant in mammalian cells. NHEJ repair frequently causes small insertions and deletions (indels) near DNA cleavage sites but only rarely causes nucleotide substitutions. High-throughput sequencing is the primary means of assessing indel and substitution frequencies in bulk populations of cells in the gene editing field. However, it is difficult to detect bona fide substitutions, which are embedded among experimentally-induced substitution errors, in high-throughput sequencing data. Here, we developed a novel analysis method, named CRISPR-Sub, to statistically detect Cas9-mediated substitutions in high-throughput sequencing data by comparing Mock- and CRISPR-treated samples. We first pinpointed ‘hotspot positions’ in target sequences at which substitution mutations were quantitatively observed much more often (p > 0.001) in CRISPR- versus Mock-treated samples. We refer to the substitution mutations in defined hotspot positions as ‘apparent substitutions’ and ultimately calculated ‘apparent substitution frequencies’ for each target. By examining 51 endogenous target sites in HeLa cells, we found that the average apparent substitution frequency was 0.8% in all queries, that apparent substitutions frequently occur near CRISPR-Cas9 cleavage sites, and that nucleotide conversion showed no meaningful nucleotide preference patterns. Furthermore, we generated NHEJ-inhibited cell lines (LIG4^−/−) by knockout of the gene encoding ligase IV and found that the apparent substitution frequencies were significantly decreased in LIG4^−/− cells, strongly suggesting that DNA substitutions are generated by the NHEJ pathway.

Activation-induced cytidine deaminase can target multiple topologies of double-stranded DNA in a transcription-independent manner

Activation-induced cytidine deaminase (AID) mutates immunoglobulin genes and acts genome-wide. AID targets robustly transcribed genes, and purified AID acts on single-stranded (ss) but not double-stranded (ds) DNA oligonucleotides. Thus, it is believed that transcription is the generator of ssDNA for AID. Previous cell-free studies examining the relationship between transcription and AID targeting have employed a bacterial colony count assay wherein AID reverts an antibiotic resistance stop codon in plasmid substrates, leading to colony formation. Here, we established a novel assay where kb-long dsDNA of varying topologies is incubated with AID, with or without transcription, followed by direct sequencing. This assay allows for an unselected and in-depth comparison of mutation frequency and pattern of AID targeting in the absence of transcription or across a range of transcription dynamics. We found that without transcription, AID targets breathing ssDNA in supercoiled and, to a lesser extent, in relaxed dsDNA. The most optimal transcription only modestly enhanced AID action on supercoiled dsDNA in a manner dependent on RNA polymerase speed. These data suggest that the correlation between transcription and AID targeting may reflect transcription leading to AID-accessible breathing ssDNA patches naturally occurring in de-chromatinized dsDNA, as much as being due to transcription directly generating ssDNA.