Assessment of whole genome amplification-induced bias through high-throughput, massively parallel whole genome sequencing

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Rescue of Mycobacterium bovis DNA Obtained from Cultured Samples during Official Surveillance of Animal TB: Key Steps for Robust Whole Genome Sequence Data Generation

Epidemiological surveillance of animal tuberculosis (TB) based on whole genome sequencing (WGS) of Mycobacterium bovis has recently gained track due to its high resolution to identify infection sources, characterize the pathogen population structure, and facilitate contact tracing. However, the workflow from bacterial isolation to sequence data analysis has several technical challenges that may severely impact the power to understand the epidemiological scenario and inform outbreak response. While trying to use archived DNA from cultured samples obtained during routine official surveillance of animal TB in Portugal, we struggled against three major challenges: the low amount of M. bovis DNA obtained from routinely processed animal samples; the lack of purity of M. bovis DNA, i.e., high levels of contamination with DNA from other organisms; and the co-occurrence of more than one M. bovis strain per sample (within-host mixed infection). The loss of isolated genomes generates missed links in transmission chain reconstruction, hampering the biological and epidemiological interpretation of data as a whole. Upon identification of these challenges, we implemented an integrated solution framework based on whole genome amplification and a dedicated computational pipeline to minimize their effects and recover as many genomes as possible. With the approaches described herein, we were able to recover 62 out of 100 samples that would have otherwise been lost. Based on these results, we discuss adjustments that should be made in official and research laboratories to facilitate the sequential implementation of bacteriological culture, PCR, downstream genomics, and computational-based methods. All of this in a time frame supporting data-driven intervention.

Impact of template denaturation prior to whole genome amplification on gene detection in high GC-content species, Burkholderia mallei and B. pseudomallei

OBJECTIVE

In this study, we sought to determine the types and prevalence of antimicrobial resistance determinants (ARDs) in Burkholderia spp. strains using the Antimicrobial Resistance Determinant Microarray (ARDM).

RESULTS

Whole genome amplicons from 22 B. mallei (BM) and 37 B. pseudomallei (BP) isolates were tested for > 500 ARDs using ARDM v.3.1. ARDM detected the following Burkholderia spp.-derived genes, aac(6), blaBP/MBL-3, blaABPS, penA-BP, and qacE, in both BM and BP while blaBP/MBL-1, macB, blaOXA-42/43 and penA-BC were observed in BP only. The method of denaturing template for whole genome amplification greatly affected the numbers and types of genes detected by the ARDM. BlaTEM was detected in nearly a third of BM and BP amplicons derived from thermally, but not chemically denatured templates. BlaTEM results were confirmed by PCR, with 81% concordance between methods. Sequences from 414-nt PCR amplicons (13 preparations) were 100% identical to the Klebsiella pneumoniae reference gene. Although blaTEM sequences have been observed in B. glumae, B. cepacia, and other undefined Burkholderia strains, this is the first report of such sequences in BM/BP/B. thailandensis (BT) clade. These results highlight the importance of sample preparation in achieving adequate genome coverage in methods requiring untargeted amplification before analysis.

Development and proof-of-concept demonstration of a clinical metagenomics method for the rapid detection of bloodstream infection

BACKGROUND

The timely and accurate diagnosis of bloodstream infection (BSI) is critical for patient management. With longstanding challenges for routine blood culture, metagenomics is a promising approach to rapidly provide sequence-based detection and characterisation of bloodborne bacteria. Long-read sequencing technologies have successfully supported the use of clinical metagenomics for syndromes such as respiratory illness, and modified approaches may address two requisite factors for metagenomics to be used as a BSI diagnostic: depletion of the high level of host DNA to then detect the low abundance of microbes in blood.

METHODS

Blood samples from healthy donors were spiked with different concentrations of four prevalent causative species of BSI. All samples were then subjected to a modified saponin-based host DNA depletion protocol and optimised DNA extraction, whole genome amplification and debranching steps in preparation for sequencing, followed by bioinformatical analyses. Two related variants of the protocol are presented: 1mL of blood processed without bacterial enrichment, and 5mL of blood processed following a rapid bacterial enrichment protocol-SepsiPURE.

RESULTS

After first identifying that a large proportion of host mitochondrial DNA remained, the host depletion process was optimised by increasing saponin concentration to 3% and scaling the reaction to allow more sample volume. Compared to non-depleted controls, the 3% saponin-based depletion protocol reduced the presence of host chromosomal and mitochondrial DNA < 106 and < 103 fold respectively. When the modified depletion method was further combined with a rapid bacterial enrichment method (SepsiPURE; with 5mL blood samples) the depletion of mitochondrial DNA improved by a further > 10X while also increasing detectable bacteria by > 10X. Parameters during DNA extraction, whole genome amplification and long-read sequencing were also adjusted, and subsequently amplicons were detected for each input bacterial species at each of the spiked concentrations, ranging from 50-100 colony forming units (CFU)/mL to 1-5 CFU/mL.

CONCLUSION

In this proof-of-concept study, four prevalent BSI causative species were detected in under 12 h to species level (with antimicrobial resistance determinants) at concentrations relevant to clinical blood samples. The use of a rapid and precise metagenomic protocols has the potential to advance the diagnosis of BSI.

ddRAD Sequencing and DNA Barcoding

Double-digest restriction site-associated DNA sequencing is a library preparation protocol that enables capturing variable sites across the genome including single-nucleotide polymorphisms (SNPs). These SNPs can be utilized to gain evolutionary insights into patterns observed in DNA barcodes, to infer population structure and phylogenies, to detect gene flow and introgression, and to perform species delimitation analyses. The protocol includes chemically shearing genomic DNA with restriction enzymes, unique tagging, size selection, and amplification of the resulting DNA fragments. Here we provide a detailed description of each step of the protocol, as well as information on essential equipment and common issues encountered during laboratory work.

Evaluation of multiple displacement amplification for metagenomic analysis of low biomass samples

Combining multiple displacement amplification (MDA) with metagenomics enables the analysis of samples with extremely low DNA concentrations, making them suitable for high-throughput sequencing. Although amplification bias and nonspecific amplification have been reported from MDA-amplified samples, the impact of MDA on metagenomic datasets is not well understood. We compared three MDA methods (i.e. bulk MDA, emulsion MDA, and primase MDA) for metagenomic analysis of two DNA template concentrations (approx. 1 and 100 pg) derived from a microbial community standard "mock community" and two low biomass environmental samples (i.e. borehole fluid and groundwater). We assessed the impact of MDA on metagenome-based community composition, assembly quality, functional profiles, and binning. We found amplification bias against high GC content genomes but relatively low nonspecific amplification such as chimeras, artifacts, or contamination for all MDA methods. We observed MDA-associated representational bias for microbial community profiles, especially for low-input DNA and with the primase MDA method. Nevertheless, similar taxa were represented in MDA-amplified libraries to those of unamplified samples. The MDA libraries were highly fragmented, but similar functional profiles to the unamplified libraries were obtained for bulk MDA and emulsion MDA at higher DNA input and across these MDA libraries for the groundwater sample. Medium to low-quality bins were possible for the high input bulk MDA metagenomes for the most simple microbial communities, borehole fluid, and mock community. Although MDA-based amplification should be avoided, it can still reveal meaningful taxonomic and functional information from samples with extremely low DNA concentration where direct metagenomics is otherwise impossible.

The Production of Recombinant African Swine Fever Virus Lv17/WB/Rie1 Strains and Their In Vitro and In Vivo Characterizations

Lv17/WB/Rie1-Δ24 was produced via illegitimate recombination mediated by low-dilution serial passage in the Cos7 cell line and isolated on PAM cell culture. The virus contains a huge ~26.4 Kb deletion in the left end of its genome. Lv17/WB/Rie1-ΔCD-ΔGL was generated via homologous recombination, crossing two ASFV strains (Lv17/WB/Rie1-ΔCD and Lv17/WB/Rie1-ΔGL containing eGFP and mCherry markers) during PAM co-infection. The presence of unique parental markers in the Lv17/WB/Rie1-ΔCD-ΔGL genome indicates at least two recombination events during the crossing, suggesting that homologous recombination is a relatively frequent event in the ASFV genome during replication in PAM. Pigs infected with Lv17/WB/Rie1-Δ24 and Lv17/WB/Rie1/ΔCD-ΔGL strains have shown mild clinical signs despite that ASFV could not be detected in their sera until a challenge infection with the Armenia/07 ASFV strain. The two viruses were not able to induce protective immunity in pigs against a virulent Armenia/07 challenge.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Flow Sorting, Whole Genome Amplification and Next-Generation Sequencing as Combined Tools to Study Heterogeneous Acute Lymphoblastic Leukemia

Next-generation sequencing (NGS) methods have been introduced for immunoglobulin (IG)/T-cell receptor (TR) gene rearrangement analysis in acute lymphoblastic leukemia (ALL) and lymphoma (LBL). These methods likely constitute faster and more sensitive approaches to analyze heterogenous cases of ALL/LBL, yet it is not known whether gene rearrangements constituting low percentages of the total sequence reads represent minor subpopulations of malignant cells or background IG/TR gene rearrangements in normal B-and T-cells. In a comparison of eight cases of B-cell precursor ALL (BCP-ALL) using both the EuroClonality NGS method and the IdentiClone multiplex-PCR/gene-scanning method, the NGS method identified between 29% and 139% more markers than the gene-scanning method, depending on whether the NGS data analysis used a threshold of 5% or 1%, respectively. As an alternative to using low thresholds, we show that IG/TR gene rearrangements in subpopulations of cancer cells can be discriminated from background IG/TR gene rearrangements in normal B-and T-cells through a combination of flow cytometry cell sorting and multiple displacement amplification (MDA)-based whole genome amplification (WGA) prior to the NGS. Using this approach to investigate the clonal evolution in a BCP-ALL patient with double relapse, clonal TR rearrangements were found in sorted leukemic cells at the time of second relapse that could be identified at the time of diagnosis, below 1% of the total sequence reads. These data emphasize that caution should be exerted when interpreting rare sequences in NGS experiments and show the advantage of employing the flow sorting of malignant cell populations in NGS clonality assessments.

A primer-independent DNA polymerase-based method for competent whole-genome amplification of intermediate to high GC sequences

Multiple displacement amplification (MDA) has proven to be a useful technique for obtaining large amounts of DNA from tiny samples in genomics and metagenomics. However, MDA has limitations, such as amplification artifacts and biases that can interfere with subsequent quantitative analysis. To overcome these challenges, alternative methods and engineered DNA polymerase variants have been developed. Here, we present new MDA protocols based on the primer-independent DNA polymerase (piPolB), a replicative-like DNA polymerase endowed with DNA priming and proofreading capacities. These new methods were tested on a genomes mixture containing diverse sequences with high-GC content, followed by deep sequencing. Protocols relying on piPolB as a single enzyme cannot achieve competent amplification due to its limited processivity and the presence of ab initio DNA synthesis. However, an alternative method called piMDA, which combines piPolB with Φ29 DNA polymerase, allows proficient and faithful amplification of the genomes. In addition, the prior denaturation step commonly performed in MDA protocols is dispensable, resulting in a more straightforward protocol. In summary, piMDA outperforms commercial methods in the amplification of genomes and metagenomes containing high GC sequences and exhibits similar profiling, error rate and variant determination as the non-amplified samples.

Y chromosome sequence and epigenomic reconstruction across human populations

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets.

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.

Histology-based and cytology-based needle sampling for targeted next-generation sequencing in the indeterminate thyroid tumors

PURPOSE

To establish the optimal and minimally invasive diagnostic approach for targeted next-generation sequencing (NGS) in the indeterminate thyroid tumors.

METHODS

The patients with indeterminate thyroid tumors were prospectively recruited and analyzed in a single tertiary medical center. We performed FNA and core needle biopsy (CNB) at the surgical specimens to confirm the quality of each sampling procedure. Cytological diagnosis by FNA, histological diagnosis by CNB and confirmed diagnosis by final surgery were compared to demonstrate the agreement among these approaches for the indeterminate thyroid tumors. The quality of the samples from FNA and CNB was evaluated, respectively to determine the optimal approach for targeted NGS. Finally, we performed ultrasound-guided CNB and FNA (US-CNB and US-FNA) on one case to confirm the clinical feasibility of being a pre-operative minimally invasive diagnostic approach.

RESULTS

A total of 6 female patients (average age: 50.83 ± 15.18 years) with indeterminate thyroid tumors (average size: 1.79 ± 0.91 cm) were recruited for further analyses. The pathological diagnoses could be obtained by CNB in the first five cases, and the sample quality of CNB for targeted NGS was better than that of FNA, even after 10X dilution. The gene mutations associated with thyroid malignancy could be detected by NGS. In the case treated with US-CNB, the pathological and targeted NGS results were successfully obtained, which suggested the possibility of thyroid malignancy to facilitate immediate decision of subsequent treatment.

CONCLUSION

CNB could serve as a minimally invasive diagnostic approach in the indeterminate thyroid tumors by providing pathological diagnoses and qualified samples for detection of mutated genes, which facilitates appropriate and immediate management.

A fast machine-learning-guided primer design pipeline for selective whole genome amplification

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Evaluating the application value of NGS-based PGT-A by screening cryopreserved MDA products of embryos from PGT-M cycles with known transfer outcomes

PURPOSE

To determine the application value of next-generation sequencing (NGS)-based preimplantation genetic testing for aneuploidies (PGT-A).

METHODS

We conducted a retrospective case-control study on a cohort of frozen-thawed embryo transfer (FET) cycles following preimplantation genetic testing for monogenic disorders (PGT-M) between 2014 and 2017. Cycles that produced live births or early miscarriages were divided into live birth group (n = 76) or miscarriage group (n = 19), respectively. The NGS-based aneuploidy screening was performed on the multiple displacement amplification (MDA) products of the embryonic trophectoderm biopsy samples that were cryopreserved following PGT-M.

RESULTS

In the live birth group, 75% (57/76) embryos were euploid and 14.5% (11/76) were aneuploid. The remaining 10.5% (8/76) embryos were NGS-classified mosaic with the high- (≥ 50%) and low-level (< 50%) mosaicism rates at 7.9% (6/76) and 2.6% (2/76), respectively. In the miscarriage group, only 23.5% (4/17) embryos were aneuploid, while 58.8% (10/17) were euploid and 17.6% (3/17) were NGS-classified mosaic with the high- and low-level mosaicism rates at 11.8% (2/17) and 5.9% (1/17), respectively. For live birth and miscarriage groups, the transferable rate was 82.9% (63/76) and 70.6% (12/17), respectively, whereas the untransferable rate was 17.1% (13/76) and 29.4% (5/17), respectively.

CONCLUSION

The application of NGS-based PGT-A remains questionable, as it may cause at least one in six embryos with reproductive potential to be discarded and prevent miscarriage in less than one in three embryos in single-gene disease carriers.

Long-read sequencing on the SMRT platform enables efficient haplotype linkage analysis in preimplantation genetic testing for β-thalassemia

PURPOSE

This study aimed to evaluate the value of long-read sequencing for preimplantation haplotype linkage analysis.

METHODS

The genetic material of the three β-thalassemia mutation carrier couples was sequenced using single-molecule real-time sequencing in the 7.7-kb region of the HBB gene and a 7.4-kb region that partially overlapped with it to detect the presence of 17 common HBB gene mutations in the Chinese population and the haplotypes formed by the continuous array of single-nucleotide polymorphisms linked to these mutations. By using the same method to analyze multiple displacement amplification products of embryos from three families and comparing the results with those of the parents, it could be revealed whether the embryos carry disease-causing mutations without the need for a proband.

RESULTS

The HBB gene mutations of the three couples were accurately detected, and the haplotype linked to the pathogenic site was successfully obtained without the need for a proband. A total of 68.75% (22/32) of embryos from the three families successfully underwent haplotype linkage analysis, and the results were consistent with the results of NGS-based mutation site detection.

CONCLUSION

This study supports long-read sequencing as a potential tool for preimplantation haplotype linkage analysis.

New Approaches to Multi-Parametric HIV-1 Genetics Using Multiple Displacement Amplification: Determining the What, How, and Where of the HIV-1 Reservoir

Development of potential HIV-1 curative interventions requires accurate characterization of the proviral reservoir, defined as host-integrated viral DNA genomes that drive rebound of viremia upon halting ART (antiretroviral therapy). Evaluation of such interventions necessitates methods capable of pinpointing the rare, genetically intact, replication-competent proviruses within a background of defective proviruses. This evaluation can be achieved by identifying the distinct integration sites of intact proviruses within host genomes and monitoring the dynamics of these proviruses and host cell lineages over longitudinal sampling. Until recently, molecular genetic approaches at the single proviral level have been generally limited to one of a few metrics, such as proviral genome sequence/intactness, host-proviral integration site, or replication competency. New approaches, taking advantage of MDA (multiple displacement amplification) for WGA (whole genome amplification), have enabled multiparametric proviral characterization at the single-genome level, including proviral genome sequence, host-proviral integration site, and phenotypic characterization of the host cell lineage, such as CD4 memory subset and antigen specificity. In this review, we will examine the workflow of MDA-augmented molecular genetic approaches to study the HIV-1 reservoir, highlighting technical advantages and flexibility. We focus on a collection of recent studies in which investigators have used these approaches to comprehensively characterize intact and defective proviruses from donors on ART, investigate mechanisms of elite control, and define cell lineage identity and antigen specificity of infected CD4+ T cell clones. The highlighted studies exemplify how these approaches and their future iterations will be key in defining the targets and evaluating the impacts of HIV curative interventions.

Fifty Generations of Amitosis: Tracing Asymmetric Allele Segregation in Polyploid Cells with Single-Cell DNA Sequencing

Amitosis is a widespread form of unbalanced nuclear division whose biomedical and evolutionary significance remain unclear. Traditionally, insights into the genetics of amitosis have been gleaned by assessing the rate of phenotypic assortment. Though powerful, this experimental approach relies on the availability of phenotypic markers. Leveraging Paramecium tetraurelia, a unicellular eukaryote with nuclear dualism and a highly polyploid somatic nucleus, we probe the limits of single-cell whole-genome sequencing to study the consequences of amitosis. To this end, we first evaluate the suitability of single-cell sequencing to study the AT-rich genome of P. tetraurelia, focusing on common sources of genome representation bias. We then asked: can alternative rearrangements of a given locus eventually assort after a number of amitotic divisions? To address this question, we track somatic assortment of developmentally acquired Internal Eliminated Sequences (IESs) up to 50 amitotic divisions post self-fertilization. To further strengthen our observations, we contrast empirical estimates of IES retention levels with in silico predictions obtained through mathematical modeling. In agreement with theoretical expectations, our empirical findings are consistent with a mild increase in variation of IES retention levels across successive amitotic divisions of the macronucleus. The modest levels of somatic assortment in P. tetraurelia suggest that IESs retention levels are largely sculpted at the time of macronuclear development, and remain fairly stable during vegetative growth. In forgoing the requirement for phenotypic assortment, our approach can be applied to a wide variety of amitotic species and could facilitate the identification of environmental and genetic factors affecting amitosis.

The Opportunities and Challenges of Molecular Tagging Next-Generation Sequencing in Liquid Biopsy

Liquid biopsy (LB) is a promising tool that is rapidly evolving as a standard of care in early and advanced stages of cancer settings. Next-generation sequencing (NGS) methods have become essential in molecular diagnostics and clinical laboratories dealing with LB analytes, i.e., cell-free DNA and RNA. The sensitivity and high-throughput capacity of NGS enable us to overcome technical issues that are mainly attributable to low-abundance (below 1% mutated allelic frequency) tumour genetic material circulating within biological fluids. In this context, the introduction of unique molecular identifiers (UMIs), also known as molecular barcodes, applied to various NGS platforms greatly improved the characterization of rare genetic alterations, as they resulted in a drastic reduction in background noise while maintaining high levels of positive predictive value and sensitivity. Different UMI strategies have been developed, such as single (e.g., safe-sequencing system, Safe-SeqS) or double (duplex-sequencing system, Duplex-Seq) strand-based labelling, and, currently, considerable results corroborate their potential implementation in a routine laboratory. Recently, the US Food and Drug Administration approved the clinical use of two comprehensive UMI-based NGS assays (FoundationOne Liquid CDx and Guardant360 CDx) in cfDNA mutational assessment. However, to definitively translate LB into clinical practice, UMI-based NGS protocols should meet certain feasibility requirements in terms of cost-effectiveness, wet laboratory performance and easy access to web-source and bioinformatic tools for downstream molecular data.

Relevance of ddRADseq method for species and population delimitation of closely related and widely distributed wolf spiders (Araneae, Lycosidae)

Although species delimitation is often controversial, emerging DNA-based and classical morphology-based methods are rarely compared using large-scale samplings, even less in the case of widely distributed species that have distant, allopatric populations. In the current study, we examined species boundaries within two wolf spider species of the genus Pardosa (Araneae, Lycosidae), P. riparia and P. palustris. Wolf spiders constitute an excellent model for testing the relevance of traditional vs. modern methods in species and population delimitation because several closely related species are distributed over cross-continental geographic ranges. Allopatric populations of the two Pardosa species were sampled across Europe to Far East Russia (latitudinal range > 150°) and several dozen individuals were studied using morphological characters (morphometry of three measures for both sexes, plus five in males only and two in females only), DNA barcoding (COI sequencing) and double-digest restriction site associated DNA sequencing (ddRADseq). The results obtained allow for changing the taxonomic status of two Far East Russian populations to subspecies and ddRADseq proved to be a powerful tool for taxonomic research despite scarce sampling and inherent subjectivity of species delimitation in allopatry. Overall, this study pleads for both multi-criteria and more population-based studies in taxonomy.

Template length, concentration and guanidine and cytosine content influence on multiple displacement amplification efficiency

Detection of low abundance human health pathogens in environmental samples is a challenge for water monitoring. This limitation can be overcome by the introduction of multiple displacement amplification (MDA) where a minute amount of genetic material can be amplified using a phi-29 DNA polymerase. However, the genetic makeup and the concentration of the polynucleotides might influence the amplification process due to inherent assay bias. Herein, a series of experiments were designed to demonstrate the effect of genome length, guanidine and cytosine content, and template concentration on the efficiency of MDA. Quantitative polymerase chain reaction (qPCR) was performed to quantify pre- and post-MDA concentrations of selected genes. Linear regression between pre- and post-MDA log gene copies L^-1 of both environmental and lab-grown samples showed a positive correlation (F = 77.59, P < 0.001, R² = 0.7, slope = 1.01). Correlation between relative polynucleotide increase after MDA and target organism length and gene target guanidine and cytosine (G + C) content (F = 4.3, P = 0.02) shows that lower G + C and higher genome length is favored in the MDA process. The MDA process was shown to favor a longer genome over a shorter genome (1.19 and 1.04 change in log gene copy L^-1, respectively) and a lower G + C content over a higher G + C content (1.11 and 0.61 change in log gene copy L^-1, respectively). There was no MDA bias observed when polynucleotides had the same G + C and genome length but different initial concentrations. This study highlights the need for increased caution when interpreting relative abundance of organisms amplified by MDA such as in next generation sequencing.

Adsorption Sequencing as a Rapid Method to Link Environmental Bacteriophages to Hosts

An important viromics challenge is associating bacteriophages to hosts. To address this, we developed adsorption sequencing (AdsorpSeq), a readily implementable method to measure phages that are preferentially adsorbed to specific host cell envelopes. AdsorpSeq thus captures the key initial infection cycle step. Phages are added to cell envelopes, adsorbed phages are isolated through gel electrophoresis, after which adsorbed phage DNA is sequenced and compared with the full virome. Here, we show that AdsorpSeq allows for separation of phages based on receptor-adsorbing capabilities. Next, we applied AdsorpSeq to identify phages in a wastewater virome that adsorb to cell envelopes of nine bacteria, including important pathogens. We detected 26 adsorbed phages including common and rare members of the virome, a minority being related to previously characterized phages. We conclude that AdsorpSeq is an effective new tool for rapid characterization of environmental phage adsorption, with a proof-of-principle application to Gram-negative host cell envelopes.

Development of a DNA-based biosensor for the fast and sensitive detection of ochratoxin A in urine

In this paper, a novel DNA-based biosensor is proposed, which is based on paramagnetic microbeads carrying an ochratoxin A (OTA) capture aptamer. A sandwich-like detection complex is linked to the capture aptamer and is able to trigger, in presence of OTA, an isothermal rolling circle amplification (RCA) reaction. This latter generated autocatalytic units with a peroxidase activity (DNAzyme) that, in presence of a proper substrate, gave a blue-coloured product visible by the naked eye. The capture aptamer, blocked onto magnetic beads, allowed the specific capture of OTA in liquid samples. The modified detection aptamer, annealed to a circularized probe, was then used to detect the toxin capture event. Indeed, in the presence of OTA and an isothermal enzyme, the circular DNA was amplified, producing a single-stranded and tandem repeated long homologous copy of its sequence. In the DNA strand, a self-catalytic structure was formed with hemin as the catalytic core, inducing the development of blue colour in the presence of ABTS and hydrogen peroxide. The results showed that the biosensor has high sensitivity and selectivity for the detection of OTA, as low as 1.09 × 10^-12 ng/mL. Moreover, the proposed biosensor was successfully used for the detection of OTA in naturally contaminated rat urine. Accuracy and repeatability data obtained in recovery experiments were satisfying, being recoveries >95% with relative standard deviations in the range 3.6-15%. For the first time, an aptasensor was successfully applied to detect OTA in biological fluids. It can be used for mycotoxin biomonitoring and assessment of individual exposure.

De novo sequence assembly requires bioinformatic checking of chimeric sequences

De novo assembly of sequence reads from next generation sequencing platforms is a common strategy for detecting presence and sequencing of viruses in biospecimens. Amplification artifacts and presence of several related viruses in the same specimen can lead to assembly of erroneous, chimeric sequences. We now report that such chimeras can also occur between viral and non-viral biological sequences incorrectly joined together which may cause erroneous detection of viruses, highlighting the importance of performing a chimera checking step in bioinformatics pipelines. Using Illumina NextSeq and metagenomic sequencing, we analyzed 80 consecutive non-melanoma skin cancers (NMSCs) from 11 immunosuppressed patients together with 11 NMSCs from patients who had only developed 1 NMSC. We aligned high-quality reads against a Human Papillomavirus (HPV) database and found HPV sequences in 9/91 specimens. A previous bioinformatic analysis of the same crude sequencing data from some of these samples had found an additional 3 specimens to be HPV-positive after performing de novo assembly. The reason for the discrepancy was investigated and found to be mostly caused by chimeric sequences containing both viral and non-viral sequences. Non-viral sequences were present in these 3 samples. To avoid erroneous detection of HPV when performing sequencing, we thus developed a novel script to identify HPV chimeric sequences.

Encoding multiple digital DNA signals in a single analog channel

For many analytic and biomedical applications, the presence of an analyte above or below a critical concentration is more informative for decision making than the actual concentration value. Straightforward analog-to-digital signal conversion does not take full advantage of the precision and dynamic range of modern sensors. Here, we present and experimentally demonstrate an analog-to-multiple-digital signal conversion, reporting digital signals that indicate whether the concentrations of specific DNA sequences exceed respective threshold values. These threshold values can be individually programmed for each target sequence. Experimentally, we showed representation of four DNA targets’ information in a single fluorescence channel.

Optimization of a WGA-Free Molecular Tagging-Based NGS Protocol for CTCs Mutational Profiling

Molecular characterization of Circulating Tumor Cells (CTCs) is still challenging, despite attempts to minimize the drawbacks of Whole Genome Amplification (WGA). In this paper, we propose a Next-Generation Sequencing (NGS) optimized protocol based on molecular tagging technology, in order to detect CTCs mutations while skipping the WGA step. MDA-MB-231 and MCF-7 cell lines, as well as leukocytes, were sorted into pools (2–5 cells) using a DEPArray™ system and were employed to set up the overall NGS procedure. A substantial reduction of reagent volume for the preparation of libraries was performed, in order to fit the limited DNA templates directly derived from cell lysates. Known variants in TP53, KRAS, and PIK3CA genes were detected in almost all the cell line pools (35/37 pools, 94.6%). No additional alterations, other than those which were expected, were found in all tested pools and no mutations were detected in leukocytes. The translational value of the optimized NGS workflow is confirmed by sequencing CTCs pools isolated from eight breast cancer patients and through the successful detection of variants. In conclusion, this study shows that the proposed NGS molecular tagging approach is technically feasible and, compared to traditional NGS approaches, has the advantage of filtering out the artifacts generated during library amplification, allowing for the reliable detection of mutations and, thus, making it highly promising for clinical use.

SCGid: a consensus approach to contig filtering and genome prediction from single-cell sequencing libraries of uncultured eukaryotes

MOTIVATION

Whole-genome sequencing of uncultured eukaryotic genomes is complicated by difficulties in acquiring sufficient amounts of tissue. Single-cell genomics (SCG) by multiple displacement amplification provides a technical workaround, yielding whole-genome libraries which can be assembled de novo. Downsides of multiple displacement amplification include coverage biases and exacerbation of contamination. These factors affect assembly continuity and fidelity, complicating discrimination of genomes from contamination and noise by available tools. Uncultured eukaryotes and their relatives are often underrepresented in large sequence data repositories, further impairing identification and separation.

RESULTS

We compare the ability of filtering approaches to remove contamination and resolve eukaryotic draft genomes from SCG metagenomes, finding significant variation in outcomes. To address these inconsistencies, we introduce a consensus approach that is codified in the SCGid software package. SCGid parallelly filters assemblies using different approaches, yielding three intermediate drafts from which consensus is drawn. Using genuine and mock SCG metagenomes, we show that our approach corrects for variation among draft genomes predicted by individual approaches and outperforms them in recapitulating published drafts in a fast and repeatable way, providing a useful alternative to available methods and manual curation.

AVAILABILITY AND IMPLEMENTATION

The SCGid package is implemented in python and R. Source code is available at http://www.github.com/amsesk/SCGid under the GNU GPL 3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Technical pitfalls that bias comparative microbial community analyses of aquatic disease Ian Hewson

The accessibility of high-throughput DNA sequencing technologies has attracted the application of comparative microbial analyses to study diseases. These studies present a window into host microbiome diversity and composition that can be used to address ecological theory in the context of host biology and behavior. Recently, comparative microbiome studies have been used to study non-vertebrate aquatic diseases to elucidate microorganisms potentially involved in disease processes or in disease prevention. These investigations suffer from many well-described biases, especially prior to sequence analyses, that could lead to misleading conclusions. Microbiome-focused studies of aquatic metazoan diseases provide valuable documentation of microbial ecology, although, they are only a starting point for establishing disease etiology, which demands quantitative validation through targeted approaches. The microbiome approach to understanding disease is most useful after laboratory diagnostics guided by pathology have failed to identify a causative agent. This opinion piece presents several technical pitfalls which may affect wider interpretation of microbe-host interactions through comparative microbial community analyses and provides recommendations, based on studies in non-aquatic systems, for incorporation into future aquatic disease research.

Testing the impact of Whole Genome Amplification on genome comparison using the polyploid flagellated Giardia duodenalis as a model

The availability of high quality genomic DNA in sufficient amounts to perform Next Generation Sequencing (NGS) experiments is challenging for pathogens that cannot be cultivated in vitro, as is the case for many parasites. Therefore, Whole Genome Amplification (WGA) of genomic DNA is used to overcome this limitation. In this study, we evaluated the effect of WGA using the intestinal flagellated protozoan Giardia duodenalis as a model, due to its genome compactness (12 Mb), the presence of two diploid nuclei with variable levels of allelic sequence heterogeneity (ASH), and the availability of reference genomes. We selected one isolate (ZX15) belonging to the same genetic group of the reference isolate WB, namely Assemblage A, sub-Assemblage AI. Genomic DNA from the ZX15 isolate (GEN dataset) and that obtained by WGA of 1 ng of the same genomic DNA (WGA dataset) were sequenced on a HiSeq Illumina platform. Trimmed reads from the GEN and WGA experiments were mapped against the WB reference genome, showing the presence of a very small number of mutations (846 and 752, respectively). The difference in the number of mutations is largely accounted by local variation in coverage and not by bias introduced by WGA. No significant difference were observed in the distribution of mutations in coding and non-coding regions, in the proportion of heterozygous mutations (ASH), or in the transition/transversion ratio of Single Nucleotide Variants within coding sequences. We conclude that the quantitative and qualitative impact of WGA on the identification of mutations is limited, and that this technique can be used to conduct comparative genomics studies.

Obtaining Genome Sequences of Mutualistic Bacteria in Single Microcystis Colonies

Cells of Microcystis are associated with heterotrophic bacteria and organized in colonies in natural environment, which are basic elements in the mass occurrence of cyanobacterial species. Analyzing these colonies by using metagenomics is helpful to understand species composition and relationship. Meanwhile, the difference in population abundance among Microcystis colonies could be used to recover genome bins from metagenome assemblies. Herein, we designed a pipeline to obtain high-quality genomes of mutualistic bacteria from single natural Microcystis colonies. Single colonies were lysed, and then amplified by using multiple displacement amplification to overcome the DNA quantity limit. A two-step assembly was performed after sequencing and scaffolds were grouped into putative bins based on their differential-coverage among species. We analyzed six natural colonies of three prevailing Microcystis species from Lake Taihu. Clustering results proved that colonies of the same species were similar in the microbial community composition. Eight putative population genome bins with wide bacterial diversity and different GC content were identified based on coverage difference among colonies. At the phylum level, proteobacteria was the most abundant besides cyanobacteria. Six of the population bins were further refined into nearly complete genomes (completeness > 90%).

Bacterial Single Cell Whole Transcriptome Amplification in Microfluidic Platform Shows Putative Gene Expression Heterogeneity

Single cell RNA sequencing is a technology that provides the capability of analyzing the transcriptome of a single cell from a population. So far, single cell RNA sequencing has been focused mostly on human cells due to the larger starting amount of RNA template for subsequent amplification. One of the major challenges of applying single cell RNA sequencing to microbial cells is to amplify the femtograms of the RNA template to obtain sufficient material for downstream sequencing with minimal contamination. To achieve this goal, efforts have been focused on multiround RNA amplification, but would introduce additional contamination and bias. In this work, we for the first time coupled a microfluidic platform with multiple displacement amplification technology to perform single cell whole transcriptome amplification and sequencing of Porphyromonas somerae, a microbe of interest in endometrial cancer, as a proof-of-concept demonstration of using single cell RNA sequencing tool to unveil gene expression heterogeneity in single microbial cells. Our results show that the bacterial single-cell gene expression regulation is distinct across different cells, supporting widespread heterogeneity.

Population Structure and Dynamics of Helminthic Infection: Schistosomiasis

While disease and outbreaks are mainly clonal for bacteria and other asexually reproducing organisms, sexual reproduction in schistosomes and other helminths usually results in unique individuals. For sexually reproducing organisms, the traits conserved in clones will instead be conserved in the group of organisms that tends to breed together, the population. While the same tools are applied to characterize DNA, how results are interpreted can be quite different at times (see another article in this collection, http://www.asmscience.org/content/journal/microbiolspec/10.1128/microbiolspec.AME-0002-2018). It is difficult to know what the real effect any control program has on the parasite population without assessing the health of this population, how they respond to the control measure, and how they recover, if they do. This review, part of the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology of Infectious Diseases, concentrates on one approach using pooled samples to study schistosome populations and shows how this and other approaches have contributed to our understanding of this parasite family's biology and epidemiology. *This article is part of a curated collection.

Genome-scale phylogenetics reveals a monophyletic Zoopagales (Zoopagomycota, Fungi)

Previous genome-scale phylogenetic analyses of Fungi have under sampled taxa from Zoopagales; this order contains many predacious or parasitic genera, and most have never been grown in pure culture. We sequenced the genomes of 4 zoopagalean taxa that are predators of amoebae, nematodes, or rotifers and the genome of one taxon that is a parasite of amoebae using single cell sequencing methods with whole genome amplification. Each genome was a metagenome, which was assembled and binned using multiple techniques to identify the target genomes. We inferred phylogenies with both super matrix and coalescent approaches using 192 conserved proteins mined from the target genomes and performed ancestral state reconstructions to determine the ancestral trophic lifestyle of the clade. Our results indicate that Zoopagales is monophyletic. Ancestral state reconstructions provide moderate support for mycoparasitism being the ancestral state of the clade.

Evaluating methods for Avian avulavirus-1 whole genome sequencing

Background

Avian avulavirus-1 (AAvV-1, previously Newcastle Disease Virus) is responsible for poultry and wild birds' disease outbreaks. Numerous whole genome sequencing methods were reported for this virus. These methods included cloning, specific primers amplification, shotgun PCR approaches, Sequence Independent Single Primer Amplification and next generation sequencing platform kits.

Methods

Three methods were used to sequence 173 Israeli Avian avulavirus-1 field isolates and one vaccine strain (VH). The sequencing was performed on Proton and Ion Torrent Personal Genome Machine and to a lesser extent, Illumina MiSeq and NextSeq sequencers. Target specific primers (SP) and Sequence Independent Single Primer Amplification (SISPA) products sequenced via the Ion torrent sequencer had a high error rate and truncated genomes. All the next generation sequencing platform sequencing kits generated high sequence accuracy and near-complete genomic size.

Results

A high level of mutations was observed in the intergenic regions between the avian avulavirus-1 genes. Within genes, multiple regions are more mutated than the Fusion region currently used for typing.

Conclusions

Our findings suggest that the whole genome sequencing by the Ion torrent sequencing kit is sufficient. However, when higher fidelity is desired, the Illumina NextSeq and Proton torrent sequencing kits were found to be preferable.

Contributions and Challenges of High Throughput qPCR for Determining Antimicrobial Resistance in the Environment: A Critical Review

Expansion in whole genome sequencing and subsequent increase in antibiotic resistance targets have paved the way of high throughput qPCR (HT-qPCR) for analyzing hundreds of antimicrobial resistance genes (ARGs) in a single run. A meta-analysis of 51 selected studies is performed to evaluate ARGs abundance trends over the last 7 years. WaferGen^TM SmartChip is found to be the most widely used HT-qPCR platform among others for evaluating ARGs. Up till now around 1000 environmental samples (excluding biological replicates) from different parts of the world have been analyzed on HT-qPCR. Calculated detection frequency and normalized ARGs abundance (ARGs/16S rRNA gene) reported in gut microbiome studies have shown a trend of low ARGs as compared to other environmental matrices. Disparities in the HT-qPCR data analysis which are causing difficulties to researchers in precise interpretation of results have been highlighted and a possible way forward for resolving them is also suggested. The potential of other amplification technologies and point of care or field deployable devices for analyzing ARGs have also been discussed in the review. Our review has focused on updated information regarding the role, current status and future perspectives of HT-qPCR in the field of antimicrobial resistance.

Pushing the limits of whole genome amplification: successful sequencing of RADseq library from a single microhymenopteran (Chalcidoidea, Trichogramma)

A major obstacle to high-throughput genotyping of microhymenoptera is their small size. As species are difficult to discriminate, and because complexes may exist, the sequencing of a pool of specimens is hazardous. Thus, one should be able to sequence pangenomic markers (e.g., RADtags) from a single specimen. To date, whole genome amplification (WGA) prior to library construction is still a necessity as at most 10 ng of DNA can be obtained from single specimens (sometimes less). However, this amount of DNA is not compatible with manufacturer's requirements for commercial kits. Here we test the accuracy of the GenomiPhi kit V2 on Trichogramma wasps by comparing RAD libraries obtained from the WGA of single specimens (F0 and F1 generation, about1 ng input DNA for the WGA (0.17-2.9 ng)) and a biological amplification of genomic material (the pool of the progeny of the F1 generation). Globally, we found that 99% of the examined loci (up to 48,189 for one of the crosses, 109 bp each) were compatible with the mode of reproduction of the studied model (haplodiploidy) and Mendelian inheritance of alleles. The remaining 1% (0.01% of the analysed nucleotides) could represent WGA bias or other experimental/analytical bias. This study shows that the multiple displacement amplification method on which the GenomiPhi kit relies, could also be of great help for the high-throughput genotyping of microhymenoptera used for biological control, or other organisms from which only a very small amount of DNA can be extracted, such as human disease vectors (e.g.,  sandflies, fleas, ticks etc.).

Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses

BACKGROUND

Viruses are key players regulating microbial ecosystems. Exploration of viral assemblages is now possible thanks to the development of metagenomics, the most powerful tool available for studying viral ecology and discovering new viruses. Unfortunately, several sources of bias lead to the misrepresentation of certain viruses within metagenomics workflows, hindering the shift from merely descriptive studies towards quantitative comparisons of communities. Therefore, benchmark studies on virus enrichment and random amplification protocols are required to better understand the sources of bias.

RESULTS

We assessed the bias introduced by viral enrichment on mock assemblages composed of seven DNA viruses, and the bias from random amplification methods on human saliva DNA viromes, using qPCR and deep sequencing, respectively. While iodixanol cushions and 0.45 μm filtration preserved the original composition of nuclease-protected viral genomes, low-force centrifugation and 0.22 μm filtration removed large viruses. Comparison of unamplified and randomly amplified saliva viromes revealed that multiple displacement amplification (MDA) induced stochastic bias from picograms of DNA template. However, the type of bias shifted to systematic using 1 ng, with only a marginal influence by amplification time. Systematic bias consisted of over-amplification of small circular genomes, and under-amplification of those with extreme GC content, a negative bias that was shared with the PCR-based sequence-independent, single-primer amplification (SISPA) method. MDA based on random priming provided by a DNA primase activity slightly outperformed those based on random hexamers and SISPA, which may reflect differences in ability to handle sequences with extreme GC content. SISPA viromes showed uneven coverage profiles, with high coverage peaks in regions with low linguistic sequence complexity. Despite misrepresentation of certain viruses after random amplification, ordination plots based on dissimilarities among contig profiles showed perfect overlapping of related amplified and unamplified saliva viromes and strong separation from unrelated saliva viromes. This result suggests that random amplification bias has a minor impact on beta diversity studies.

CONCLUSIONS

Benchmark analyses of mock and natural communities of viruses improve understanding and mitigate bias in metagenomics surveys. Bias induced by random amplification methods has only a minor impact on beta diversity studies of human saliva viromes.

Using metagenomics to investigate human and environmental resistomes

Antibiotic resistance is a global health concern declared by the WHO as one of the largest threats to modern healthcare. In recent years, metagenomic DNA sequencing has started to be applied as a tool to study antibiotic resistance in different environments, including the human microbiota. However, a multitude of methods exist for metagenomic data analysis, and not all methods are suitable for the investigation of resistance genes, particularly if the desired outcome is an assessment of risks to human health. In this review, we outline the current state of methods for sequence handling, mapping to databases of resistance genes, statistical analysis and metagenomic assembly. In addition, we provide an overview of important considerations related to the analysis of resistance genes, and recommend some of the currently used tools and methods that are best equipped to inform research and clinical practice related to antibiotic resistance.

Effective BAC clone anchoring with genotyping-by-sequencing and Diversity Arrays Technology in a large genome cereal rye

Identification of bacterial artificial chromosome (BAC) clones containing specific sequences is a prerequisite for many applications, such as physical map anchoring or gene cloning. Existing BAC library screening strategies are either low-throughput or require a considerable initial input of resources for platform establishment. We describe a high-throughput, reliable, and cost-effective BAC library screening approach deploying genotyping platforms which are independent from the availability of sequence information: a genotyping-by-sequencing (GBS) method DArTSeq and the microarray-based Diversity Arrays Technology (DArT). The performance of these methods was tested in a very large and complex rye genome. The DArTseq approach delivered superior results: a several fold higher efficiency of addressing genetic markers to BAC clones and anchoring of BAC clones to genetic map and also a higher reliability. Considering the sequence independence of the platform, the DArTseq-based library screening can be proposed as an attractive method to speed up genomics research in resource poor species.

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build-up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double-digest restriction site-associated DNA sequencing (ddRADseq) to obtain and analyse genomic-scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.

Functional metagenomics reveals a novel carbapenem-hydrolyzing mobile beta-lactamase from Indian river sediments contaminated with antibiotic production waste

Evolution has provided environmental bacteria with a plethora of genes that give resistance to antibiotic compounds. Under anthropogenic selection pressures, some of these genes are believed to be recruited over time into pathogens by horizontal gene transfer. River sediment polluted with fluoroquinolones and other drugs discharged from bulk drug production in India constitute an environment with unprecedented, long-term antibiotic selection pressures. It is therefore plausible that previously unknown resistance genes have evolved and/or are promoted here. In order to search for novel resistance genes, we therefore analyzed such river sediments by a functional metagenomics approach. DNA fragments providing resistance to different antibiotics in E. coli were sequenced using Sanger and PacBio RSII platforms. We recaptured the majority of known antibiotic resistance genes previously identified by open shot-gun metagenomics sequencing of the same samples. In addition, seven novel resistance gene candidates (six beta-lactamases and one amikacin resistance gene) were identified. Two class A beta-lactamases, bla_RSA1 and bla_RSA2, were phylogenetically close to clinically important ESBLs like bla_GES, bla_BEL and bla_L2, and were further characterized for their substrate spectra. The blaRSA1 protein, encoded as an integron gene cassette, efficiently hydrolysed penicillins, first generation cephalosporins and cefotaxime, while blaRSA2 was an inducible class A beta-lactamase, capable of hydrolyzing carbapenems albeit with limited efficiency, similar to the L2 beta-lactamase from Stenotrophomonas maltophilia. All detected novel genes were associated with plasmid mobilization proteins, integrons, and/or other resistance genes, suggesting a potential for mobility. This study provides insight into a resistome shaped by an exceptionally strong and long-term antibiotic selection pressure. An improved knowledge of mobilized resistance factors in the external environment may make us better prepared for the resistance challenges that we may face in clinics in the future.

Microfluidic Device for Aptamer-Based Cancer Cell Capture and Genetic Mutation Detection

We present a microfluidic device for specifically capturing cancer cells and isolating their genomic DNA (gDNA) for specific amplification and sequence analysis. To capture cancer cells within the device, nucleic acid aptamers that specifically bind to cancer cells were immobilized within a channel containing micropillars designed to increase capture efficiency. The captured cells were lysed in situ, and their gDNA was isolated by physical entanglement within a second smaller-dimensioned micropillar array. This type of isolation allows the gDNA to be retained and purified within the channel and enables amplification and analysis to be performed on the gDNA without the loss of the original template. We developed a technique for selectively amplifying genes from whole gDNA using multiple displacement amplification. The amplified gene samples were sequenced, and the resulting sequence information was compared against the known wild-type gene to identify any mutations. We have tested cervical and ovarian cancer cells for mutations in the TP53 gene using this technology. This approach offers a way to monitor multiple genetic mutations in the same small population of cells, which is beneficial given the wide diversity in cancer cells, and therefore it requires very few cells to be extracted from a patient sample.

Evaluation of single-cell genomics to address evolutionary questions using three SAGs of the choanoflagellate Monosiga brevicollis

Single-cell genomics (SCG) appeared as a powerful technique to get genomic information from uncultured organisms. However, SCG techniques suffer from biases at the whole genome amplification step that can lead to extremely variable numbers of genome recovery (5-100%). Thus, it is unclear how useful can SCG be to address evolutionary questions on uncultured microbial eukaryotes. To provide some insights into this, we here analysed 3 single-cell amplified genomes (SAGs) of the choanoflagellate Monosiga brevicollis, whose genome is known. Our results show that each SAG has a different, independent bias, yielding different levels of genome recovery for each cell (6-36%). Genes often appear fragmented and are split into more genes during annotation. Thus, analyses of gene gain and losses, gene architectures, synteny and other genomic features can not be addressed with a single SAG. However, the recovery of phylogenetically-informative protein domains can be up to 55%. This means SAG data can be used to perform accurate phylogenomic analyses. Finally, we also confirm that the co-assembly of several SAGs improves the general genomic recovery. Overall, our data show that, besides important current limitations, SAGs can still provide interesting and novel insights from poorly-known, uncultured organisms.