Plasmodium falciparum histidine rich protein 2 (pfhrp2): an additional genetic marker suitable for anti-malarial drug efficacy trials

Background Genotyping of the three Plasmodium falciparum polymorphic genes, msp1, msp2 and glurp, has been adopted as a standard strategy to distinguish recrudescence from new infection in drug efficacy clinical trials. However, the suitability of a particular gene is compromised in areas where its allelic variants distribution is significantly skewed, a phenomenon that might occur in isolated parasite populations or in areas of very low transmission. Moreover, observation of amplification bias has diminished the value of glurp as a marker. Methods The suitability of the polymorphic P. falciparum histidine-rich protein 2 (pfhrp2) gene was assessed to serve as an alternative marker using a PCR-sequencing or a PCR–RFLP protocol for genotyping of samples in drug efficacy clinical trials. The value of pfhrp2 was validated by side-by-side analyses of 5 admission-recrudescence sample pairs from Yemeni malaria patients. Results The outcome of the single pfhrp2 gene discrimination analysis has been found consistent with msp1, msp2 and glurp pool genotyping analysis for the differentiation of recrudescence from new infection. Conclusion The findings suggest that under the appropriate circumstances, pfhrp2 can serve as an additional molecular marker for monitoring anti-malarials efficacy. However, its use is restricted to endemic areas where only a minority of P. falciparum parasites lack the pfhrp2 gene.

Integration of mRNA and miRNA Analysis Reveals the Molecular Mechanism of Cotton Response to Salt Stress

To identify the regulatory network of known and novel microRNAs (miRNAs) and their targets responding to salt stress, a combined analysis of mRNA libraries, small RNA libraries, and degradome libraries were performed. In this study, we used unique molecular identifiers (UMIs), which are more sensitive, accurate, and reproducible than traditional methods of sequencing, to quantify the number of molecules and correct for amplification bias. We identified a total of 312 cotton miRNAs using seedlings at 0, 1, 3, and 6 h after NaCl treatment, including 80 known ghr-miRNAs and 232 novel miRNAs and found 155 miRNAs that displayed significant differential expression under salt stress. Among them, fifty-nine differentially expressed miRNAs were simultaneously induced in two or three tissues, while 66, 11, and 19 were specifically expressed in the roots, leaves, and stems, respectively. It is indicated there were different populations of miRNAs against salt stress in roots, leaves and stems. 399 candidate targets of salt-induced miRNAs showed significant differential expression before and after salt treatment, and 72 targets of 25 miRNAs were verified by degradome sequencing data. Furthermore, the regulatory relationship of miRNA-target gene was validated experimentally via 5′RLM-RACE, proving our data reliability. Gene ontology and KEGG pathway analysis found that salt-responsive miRNA targets among the differentially expressed genes were significantly enriched, and mainly involved in response to the stimulus process and the plant hormone signal transduction pathway. Furthermore, the expression levels of newly identified miRNA mir1 and known miRNAs miR390 and miR393 gradually decreased when subjected to continuous salt stress, while overexpression of these miRNAs both increased sensitivity to salt stress. Those newly identified miRNAs and mRNA pairs were conducive to genetic engineering and better understanding the mechanisms responding to salt stress in cotton.

Accurate and scalable variant calling from single cell DNA sequencing data with ProSolo

Accurate single cell mutational profiles can reveal genomic cell-to-cell heterogeneity. However, sequencing libraries suitable for genotyping require whole genome amplification, which introduces allelic bias and copy errors. The resulting data violates assumptions of variant callers developed for bulk sequencing. Thus, only dedicated models accounting for amplification bias and errors can provide accurate calls. We present ProSolo for calling single nucleotide variants from multiple displacement amplified (MDA) single cell DNA sequencing data. ProSolo probabilistically models a single cell jointly with a bulk sequencing sample and integrates all relevant MDA biases in a site-specific and scalable—because computationally efficient—manner. This achieves a higher accuracy in calling and genotyping single nucleotide variants in single cells in comparison to state-of-the-art tools and supports imputation of insufficiently covered genotypes, when downstream tools cannot handle missing data. Moreover, ProSolo implements the first approach to control the false discovery rate reliably and flexibly. ProSolo is implemented in an extendable framework, with code and usage at: https://github.com/prosolo/prosolo

Is amplification bias consequential in transposon sequencing (TnSeq) assays? A case study with a Staphylococcus aureus TnSeq library subjected to PCR-based and amplification-free enrichment methods

As transposon sequencing (TnSeq) assays have become prolific in the microbiology field, it is of interest to scrutinize their potential drawbacks. TnSeq data consist of millions of nucleotide sequence reads that are generated by PCR amplification of transposon-genomic junctions. Reads mapping to the junctions are enumerated thus providing information on the number of transposon insertion mutations in each individual gene. Here we explore the possibility that PCR amplification of transposon insertions in a TnSeq library skews the results by introducing bias into the detection and/or enumeration of insertions. We compared the detection and frequency of mapped insertions when altering the number of PCR cycles, and when including a nested PCR, in the enrichment step. Additionally, we present nCATRAs - a novel, amplification-free TnSeq method where the insertions are enriched via CRISPR/Cas9-targeted transposon cleavage and subsequent Oxford Nanopore MinION sequencing. nCATRAs achieved 54 and 23% enrichment of the transposons and transposon-genomic junctions, respectively, over background genomic DNA. These PCR-based and PCR-free experiments demonstrate that, overall, PCR amplification does not significantly bias the results of TnSeq insofar as insertions in the majority of genes represented in our library were similarly detected regardless of PCR cycle number and whether or not PCR amplification was employed. However, the detection of a small subset of genes which had been previously described as essential is sensitive to the number of PCR cycles. We conclude that PCR-based enrichment of transposon insertions in a TnSeq assay is reliable, but researchers interested in profiling putative essential genes should carefully weigh the number of amplification cycles employed in their library preparation protocols. In addition, nCATRAs is comparable to traditional PCR-based methods (Kendall’s correlation=0.896–0.897) although the latter remain superior owing to their accessibility and high sequencing depth.

Amplification-free detection of viral RNA by super resolution imaging-based CRISPR/Cas13a System

RNA detection is crucial for biological research and clinical diagnosis. The current methods include both direct and amplification-based RNA detection. These methods require complicated procedures, suffering from low sensitivity, slow turnaround, and amplification bias. The CRISPR/Cas13a system is a direct RNA detection method via target RNA induced collateral cleavage activity. However, to detect low concentration RNA with CRISPR/Cas13a, target amplification is always required. Herein, we optimize the components of the CRISPR/Cas13a assay to enhance the sensitivity of viral RNA detection which improve the detection limit from 1 pM up to 100 fM. In addition, the integration of CRISPR/Cas13a biosensing and single molecule super resolution imaging is a novel strategy for direct and amplification-free RNA detection. After surface modification, fluorescent RNA reporters are immobilized on the glass coverslip surface and fluorescent signals are captured by total internal reflection fluorescence microscopy (TIRFM), shifting the measurement from spectroscopy to imaging. We quantify the fluorescence signal intensity before and after collateral cleavage of the CRISPR system when viral RNA is present and achieve a detection limit of 10 fM. Therefore, we provide a novel TIRFM-based system to visualize the CRISPR trans-cleavage for direct and robust RNA detection.

High-Throughput Sequencing of Phage Display Libraries Reveals Parasitic Enrichment of Indel Mutants Caused by Amplification Bias

The combination of phage display technology with high-throughput sequencing enables in-depth analysis of library diversity and selection-driven dynamics. We applied short-read sequencing of the mutagenized region on focused display libraries of two homologous nucleic acid modification eraser proteins—AlkB and FTO—biopanned against methylated DNA. This revealed enriched genotypes with small indels and concomitant doubtful amino acid motifs within the FTO library. Nanopore sequencing of the entire display vector showed additional enrichment of large deletions overlooked by region-specific sequencing, and further impacted the interpretation of the obtained amino acid motifs. We could attribute enrichment of these corrupted clones to amplification bias due to arduous FTO display slowing down host cell growth as well as phage production. This amplification bias appeared to be stronger than affinity-based target selection. Recommendations are provided for proper sequence analysis of phage display data, which can improve motive discovery in libraries of proteins that are difficult to display.

Long-read whole genome analysis of human single cells

With long-read sequencing we have entered an era where individual genomes are routinely assembled to near-completion and where complex genetic variation can efficiently be resolved. Here we demonstrate that long reads can be applied also to study the genomic architecture of individual human cells. Clonally expanded CD8+ T-cells from a human donor were used as starting material for a droplet-based multiple displacement amplification (dMDA) method designed to ensure long molecule lengths and minimal amplification bias. Sequencing of two single cells was performed on the PacBio Sequel II system, generating over 2.5 million reads and ~20Gb HiFi data (>QV20) per cell, achieving up to 40% genome coverage. This data allowed for single nucleotide variant (SNV) detection, including in genomic regions inaccessible by short reads. Over 1000 high-confidence structural variants (SVs) per cell were discovered in the PacBio data, which is four times more than the number of SVs detected in Illumina dMDA data from clonally related cells. In addition, several putative clone-specific somatic SV events could be identified. Single-cell de novo assembly resulted in 454-598 Mb assembly sizes and 35-42 kb contig N50 values. 1762 (12.8%) of expected gene models were found to be complete in the best single-cell assembly. The de novo constructed mitochondrial genomes were 100% identical for the two single cells subjected to PacBio sequencing, although mitochondrial heteroplasmy was also observed. In summary, the work presented here demonstrates the utility of long-read sequencing towards understanding the extent and distribution of complex genetic variation at the single cell level.

Computational Evaluation of DNA Metabarcoding for Universal Diagnostics of Invasive Insect Pests

Appropriate design and selection of PCR primers plays a critical role in determining the sensitivity and specificity of a metabarcoding assay. Despite several studies applying metabarcoding to insect pest surveillance, the diagnostic performance of the short "mini-barcodes" required by high-throughput sequencing platforms has not been established across the broader taxonomic diversity of invasive insects. We address this by computationally evaluating the diagnostic sensitivity and predicted amplification bias for 68 published and novel cytochrome c oxidase subunit 1 (COI) primers on a curated database of 110,676 insect species, including 2,625 registered on global invasive species lists. We find that mini-barcodes between 125-257 bp can provide comparable resolution to the full-length barcode for both invasive insect pests and the broader Insecta, conditional upon the subregion of COI targeted and the genetic similarity threshold used to identify species. Taxa that could not be identified by any barcode lengths were phylogenetically clustered within "problem groups", many arising through taxonomic inconsistencies rather than insufficient diagnostic information within the barcode itself. Substantial variation in predicted PCR bias was seen across published primers, with those including 4-5 degenerate nucleotide bases showing almost no mismatch to major insect orders. While not completely universal, a single COI mini-barcode can successfully differentiate the majority of pest and non-pest insects from their congenerics, even at the small amplicon size imposed by 2 x 150 bp sequencing. We provide a ranked summary of high-performing primers and discuss the bioinformatic steps required to curate reliable reference databases for metabarcoding studies.

Applicability of DNA-based identifications for the WFD-guided monitoring using macroinvertebrates: a large-scale DNA metabarcoding study for implementing routine ecological status assessments in Iberian rivers

Over the last decade, remarkable improvements have been made in the field of metabarcoding-based tools for routine ecological status assessments. However, important issues are yet to be solved to fulfil the European Water Framework Directive (WFD) requirements and standards. These limitations, which include problems related to e.g. the lack of a complete COI macroinvertebrate barcode database available for the Iberian Peninsula Murria 2020, or the scarce recovery of specific taxa due to DNA extraction and/or PCR amplification bias, are especially difficult to overcome for routine freshwater macroinvertebrate monitoring. For that purpose, a large-scale study is on going to test how metabarcoding data can infer existing macroinvertebrate morphotaxonomy-based biotic indexes and ecological status of Iberian rivers. Freshwater macroinvertebrates were selected as a Biological Quality Element and identified by using both morphological and metabarcoding approaches. The mitochondrial gene for cytochrome c oxidase subunit I (COI) was used as a DNA Barcode. Taxonomic coverage, taxonomic composition metrics and ecological status obtained from both approaches were analysed. Physical and chemical variables obtained during the routine biomonitoring, as well as other ecological parameters including biodiversity indexes, were also assessed. Multivariate data analysis of these environmental and biotic data obtained from both approaches were compared. Results seem to support the hypothesis Kuntke 2019 that the DNA-metabarcoding approach might deliver similar quality assessments results to the morphological approach, though some refinement must be done at the different steps of the process prior to establish a reliable procedure allowing the alternative use of both methods giving similar results for the ecological status classes marked by the WFD.

MIPP-Seq: ultra-sensitive rapid detection and validation of low-frequency mosaic mutations

Background Mosaic mutations contribute to numerous human disorders. As such, the identification and precise quantification of mosaic mutations is essential for a wide range of research applications, clinical diagnoses, and early detection of cancers. Currently, the low-throughput nature of single allele assays (e.g., allele-specific ddPCR) commonly used for genotyping known mutations at very low alternate allelic fractions (AAFs) have limited the integration of low-level mosaic analyses into clinical and research applications. The growing importance of mosaic mutations requires a more rapid, low-cost solution for mutation detection and validation. Methods To overcome these limitations, we developed Multiple Independent Primer PCR Sequencing (MIPP-Seq) which combines the power of ultra-deep sequencing and truly independent assays. The accuracy of MIPP-seq to quantifiable detect and measure extremely low allelic fractions was assessed using a combination of SNVs, insertions, and deletions at known allelic fractions in blood and brain derived DNA samples. Results The Independent amplicon analyses of MIPP-Seq markedly reduce the impact of allelic dropout, amplification bias, PCR-induced, and sequencing artifacts. Using low DNA inputs of either 25 ng or 50 ng of DNA, MIPP-Seq provides sensitive and quantitative assessments of AAFs as low as 0.025% for SNVs, insertion, and deletions. Conclusions MIPP-Seq provides an ultra-sensitive, low-cost approach for detecting and validating known and novel mutations in a highly scalable system with broad utility spanning both research and clinical diagnostic testing applications. The scalability of MIPP-Seq allows for multiplexing mutations and samples, which dramatically reduce costs of variant validation when compared to methods like ddPCR. By leveraging the power of individual analyses of multiple unique and independent reactions, MIPP-Seq can validate and precisely quantitate extremely low AAFs across multiple tissues and mutational categories including both indels and SNVs. Furthermore, using Illumina sequencing technology, MIPP-seq provides a robust method for accurate detection of novel mutations at an extremely low AAF.