High molecular weight DNA extraction strategies for long-read sequencing of complex metagenomes

Comparison of nine extraction methods for bacterial identification using the ONT MinION sequencer

The Anthrax mailings bioterrorism attack in 2001 revealed the need for universal and rapid microbial forensic analyses on unknown biological evidence. However, the gold standard for bacterial identification includes culturing isolates, which is laborious. Molecular approaches for bacterial identification revolve around 16S ribosomal gene sequencing using Sanger or next generation sequencing (NGS) platforms, but these techniques are laboratory-based and can also be time-consuming. The Oxford Nanopore Technologies (ONT) MinION sequencer can generate long read lengths that span the entire bacterial 16S rRNA gene and accurately identify the species level. This platform can be used in the field, allowing on-site evidence analysis. However, it requires higher quantities of pure DNA compared to other sequencing platforms; thus, the extraction method for bacterial DNA is critical for downstream analysis, which to date are tailored toward a priori knowledge of the species' taxonomic grouping. During an attack, the investigative team may not know what species they are handling; therefore, identifying an extraction method that can handle all bacterial groups and generate clean DNA for the MinION is useful for microbial forensic analysis. The purpose of this study was to identify a "universal" extraction method that can be coupled with ONT MinION sequencing for use in forensic situations for rapid identification. It also evaluated the cloud-based data analysis software provided by ONT, EPI2ME. No "universal" extraction method was identified as optimal for downstream MinION sequencing. However, the DNeasy PowerSoil Kit and Noda et al. Chelex-100 method gave comparable sequencing results and could be used as rapid extraction techniques. This study showed that the ONT 16S Barcoding Kit 1-24 coupled with the 16S FASTQ workflow might not be the best for use in forensic situations where species-level identification needs to be obtained, as most alignments were approximately 89% accurate. In all seven test organisms and nine extraction methods, accurate species identification was only obtained in 63% of the cases.

Diverse plasmid systems and their ecology across human gut metagenomes revealed by PlasX and MobMess

Plasmids alter microbial evolution and lifestyles by mobilizing genes that often confer fitness in changing environments across clades. Yet our ecological and evolutionary understanding of naturally occurring plasmids is far from complete. Here we developed a machine-learning model, PlasX, which identified 68,350 non-redundant plasmids across human gut metagenomes and organized them into 1,169 evolutionarily cohesive 'plasmid systems' using our sequence containment-aware network-partitioning algorithm, MobMess. Individual plasmids were often country specific, yet most plasmid systems spanned across geographically distinct human populations. Cargo genes in plasmid systems included well-known determinants of fitness, such as antibiotic resistance, but also many others including enzymes involved in the biosynthesis of essential nutrients and modification of transfer RNAs, revealing a wide repertoire of likely fitness determinants in complex environments. Our study introduces computational tools to recognize and organize plasmids, and uncovers the ecological and evolutionary patterns of diverse plasmids in naturally occurring habitats through plasmid systems.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

KombOver: Efficient k-core and K-truss based characterization of perturbations within the human gut microbiome

The microbes present in the human gastrointestinal tract are regularly linked to human health and disease outcomes. Thanks to technological and methodological advances in recent years, metagenomic sequencing data, and computational methods designed to analyze metagenomic data, have contributed to improved understanding of the link between the human gut microbiome and disease. However, while numerous methods have been recently developed to extract quantitative and qualitative results from host-associated microbiome data, improved computational tools are still needed to track microbiome dynamics with short-read sequencing data. Previously we have proposed KOMB as a de novo tool for identifying copy number variations in metagenomes for characterizing microbial genome dynamics in response to perturbations. In this work, we present KombOver (KO), which includes four key contributions with respect to our previous work: (i) it scales to large microbiome study cohorts, (ii) it includes both k-core and K-truss based analysis, (iii) we provide the foundation of a theoretical understanding of the relation between various graph-based metagenome representations, and (iv) we provide an improved user experience with easier-to-run code and more descriptive outputs/results. To highlight the aforementioned benefits, we applied KO to nearly 1000 human microbiome samples, requiring less than 10 minutes and 10 GB RAM per sample to process these data. Furthermore, we highlight how graph-based approaches such as k-core and K-truss can be informative for pinpointing microbial community dynamics within a myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) cohort. KO is open source and available for download/use at: https://github.com/treangenlab/komb.

Dancing the Nanopore limbo - Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BACKGROUND

While genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.

RESULTS

We generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.

CONCLUSION

Our results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

Evaluation of an Adapted Semi-Automated DNA Extraction for Human Salivary Shotgun Metagenomics

Recent attention has highlighted the importance of oral microbiota in human health and disease, e.g., in Parkinson's disease, notably using shotgun metagenomics. One key aspect for efficient shotgun metagenomic analysis relies on optimal microbial sampling and DNA extraction, generally implementing commercial solutions developed to improve sample collection and preservation, and provide high DNA quality and quantity for downstream analysis. As metagenomic studies are today performed on a large number of samples, the next evolution to increase study throughput is with DNA extraction automation. In this study, we proposed a semi-automated DNA extraction protocol for human salivary samples collected with a commercial kit, and compared the outcomes with the DNA extraction recommended by the manufacturer. While similar DNA yields were observed between the protocols, our semi-automated DNA protocol generated significantly higher DNA fragment sizes. Moreover, we showed that the oral microbiome composition was equivalent between DNA extraction methods, even at the species level. This study demonstrates that our semi-automated protocol is suitable for shotgun metagenomic analysis, while allowing for improved sample treatment logistics with reduced technical variability and without compromising the structure of the oral microbiome.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Dynamic genetic adaptation of Bacteroides thetaiotaomicron during murine gut colonization

To understand how a bacterium ultimately succeeds or fails in adapting to a new host, it is essential to assess the temporal dynamics of its fitness over the course of colonization. Here, we introduce a human-derived commensal organism, Bacteroides thetaiotaomicron (Bt), into the guts of germ-free mice to determine whether and how the genetic requirements for colonization shift over time. Combining a high-throughput functional genetics assay and transcriptomics, we find that gene usage changes drastically during the first days of colonization, shifting from high expression of amino acid biosynthesis genes to broad upregulation of diverse polysaccharide utilization loci. Within the first week, metabolism becomes centered around utilization of a predominant dietary oligosaccharide, and these changes are largely sustained through 6 weeks of colonization. Spontaneous mutations in wild-type Bt also evolve around this locus. These findings highlight the importance of considering temporal colonization dynamics in developing more effective microbiome-based therapies.

16S rRNA gene-based microbiota profiles from diverse avian faeces are largely independent of DNA preservation and extraction method

The avian gut microbiota has been the subject of considerable recent attention, with potential implications for diverse fields such as the poultry industry, microbial ecology, and conservation. Faecal microbiotas are frequently used as a non-invasive proxy for the gut microbiota, however the extraction of high-quality microbial DNA from avian faeces has often proven challenging. Here we aimed to evaluate the performance of two DNA preservation methods (95% ethanol and RNAlater) and five extraction approaches (IndiSpin Pathogen Kit, QIAamp PowerFecal Pro DNA Kit, MicroGEM PrepGEM Bacteria Kit, ZymoBIOMICS DNA Miniprep Kit, and an in-house phase separation-based method) for studying the avian gut microbiota. Systematic testing of the efficacy of these approaches on faecal samples from an initial three avian species (chicken, ostrich, and the flightless parrot kākāpō) revealed substantial differences in the quality, quantity and integrity of extracted DNA, but negligible influence of applied method on 16S rRNA gene-based microbiota profiles. Subsequent testing with a selected combination of preservation and extraction method on 10 further phylogenetically and ecologically diverse avian species reiterated the efficacy of the chosen approach, with bacterial community structure clustering strongly by technical replicates for a given avian species. Our finding that marked differences in extraction efficacy do not appear to influence 16S rRNA gene-based bacterial community profiles provides an important foundation for ongoing research on the avian gut microbiota.

Comparison of 6 DNA extraction methods for isolation of high yield of high molecular weight DNA suitable for shotgun metagenomics Nanopore sequencing to detect bacteria

BACKGROUND

Oxford Nanopore Technologies (ONT) offers an accessible platform for long-read sequencing, which improves the reconstruction of genomes and helps to resolve complex genomic contexts, especially in the case of metagenome analysis. To take the best advantage of long-read sequencing, DNA extraction methods must be able to isolate pure high molecular weight (HMW) DNA from complex metagenomics samples, without introducing any bias. New methods released on the market, and protocols developed at the research level, were specifically designed for this application and need to be assessed.

RESULTS

In this study, with different bacterial cocktail mixes, analyzed as pure or spiked in a synthetic fecal matrix, we evaluated the performances of 6 DNA extraction methods using various cells lysis and purification techniques, from quick and easy, to more time-consuming and gentle protocols, including a portable method for on-site application. In addition to the comparison of the quality, quantity and purity of the extracted DNA, the performance obtained when doing Nanopore sequencing on a MinION flow cell was also tested. From the obtained results, the Quick-DNA HMW MagBead Kit (Zymo Research) was selected as producing the best yield of pure HMW DNA. Furthermore, this kit allowed an accurate detection, by Nanopore sequencing, of almost all the bacterial species present in a complex mock community.

CONCLUSION

Amongst the 6 tested methods, the Quick-DNA HMW MagBead Kit (Zymo Research) was considered as the most suitable for Nanopore sequencing and would be recommended for bacterial metagenomics studies using this technology.

Application of nanopore adaptive sequencing in pathogen detection of a patient with Chlamydia psittaci infection

Introduction

Nanopore sequencing has been widely used in clinical metagenomic sequencing for pathogen detection with high portability and real-time sequencing. Oxford Nanopore Technologies has recently launched an adaptive sequencing function, which can enrich on-target reads through real-time alignment and eject uninteresting reads by reversing the voltage across the nanopore. Here we evaluated the utility of adaptive sequencing in clinical pathogen detection.

Methods

Nanopore adaptive sequencing and standard sequencing was performed on a same flow cell with a bronchoalveolar lavage fluid sample from a patient with Chlamydia psittacosis infection, and was compared with the previous mNGS results.

Results

Nanopore adaptive sequencing identified 648 on-target stop receiving reads with the longest median read length(688bp), which account for 72.4% of all Chlamydia psittaci reads and 0.03% of total reads in enriched group. The read proportion matched to C. psittaci in the stop receiving group was 99.85%, which was much higher than that of the unblock (<0.01%) and fail to adapt (0.02%) groups. Nanopore adaptive sequencing generated similar data yield of C. psittaci compared with standard nanopore sequencing. The proportion of C. psittaci reads in adaptive sequencing is close to that of standard nanopore sequencing and mNGS, but generated lower genome coverage than mNGS.

Discussion

Nanopore adaptive sequencing can effectively identify target C. psittaci reads in real-time, but how to increase the targeted data of pathogens still needs to be further evaluated.

Capturing marine microbiomes and environmental DNA: A field sampling guide

The expanding interest in marine microbiome and eDNA sequence data has led to a demand for sample collection and preservation standard practices to enable comparative assessments of results across studies and facilitate meta-analyses. We support this effort by providing guidelines based on a review of published methods and field sampling experiences. The major components considered here are environmental and resource considerations, sample processing strategies, sample storage options, and eDNA extraction protocols. It is impossible to provide universal recommendations considering the wide range of eDNA applications; rather, we provide information to design fit-for-purpose protocols. To manage scope, the focus here is on sampling collection and preservation of prokaryotic and microeukaryotic eDNA. Even with a focused view, the practical utility of any approach depends on multiple factors, including habitat type, available resources, and experimental goals. We broadly recommend enacting rigorous decontamination protocols, pilot studies to guide the filtration volume needed to characterize the target(s) of interest and minimize PCR inhibitor collection, and prioritizing sample freezing over (only) the addition of preservation buffer. An annotated list of studies that test these parameters is included for more detailed investigation on specific steps. To illustrate an approach that demonstrates fit-for-purpose methodologies, we provide a protocol for eDNA sampling aboard an oceanographic vessel. These guidelines can aid the decision-making process for scientists interested in sampling and sequencing marine microbiomes and/or eDNA.

A Fast and Easy Method to Co-extract DNA and RNA from an Environmental Microbial Sample

We herein propose a fast and easy DNA and RNA co-extraction method for environmental microbial samples. It combines bead beating and phenol-chloroform phase separation followed by the separation and purification of DNA and RNA using the Qiagen AllPrep DNA/RNA mini kit. With a handling time of ~3 h, our method simultaneously extracted high-quality DNA (peak size >10-15‍ ‍kb) and RNA (RNA integrity number >6) from lake bacterioplankton filtered samples. The method is also applicable to low-biomass samples (expected DNA or RNA yield <50‍ ‍ng) and eukaryotic microbial samples, providing an easy option for more versatile eco-genomic applications.

Target-enriched long-read sequencing (TELSeq) contextualizes antimicrobial resistance genes in metagenomes

BACKGROUND

Metagenomic data can be used to profile high-importance genes within microbiomes. However, current metagenomic workflows produce data that suffer from low sensitivity and an inability to accurately reconstruct partial or full genomes, particularly those in low abundance. These limitations preclude colocalization analysis, i.e., characterizing the genomic context of genes and functions within a metagenomic sample. Genomic context is especially crucial for functions associated with horizontal gene transfer (HGT) via mobile genetic elements (MGEs), for example antimicrobial resistance (AMR). To overcome this current limitation of metagenomics, we present a method for comprehensive and accurate reconstruction of antimicrobial resistance genes (ARGs) and MGEs from metagenomic DNA, termed target-enriched long-read sequencing (TELSeq).

RESULTS

Using technical replicates of diverse sample types, we compared TELSeq performance to that of non-enriched PacBio and short-read Illumina sequencing. TELSeq achieved much higher ARG recovery (>1,000-fold) and sensitivity than the other methods across diverse metagenomes, revealing an extensive resistome profile comprising many low-abundance ARGs, including some with public health importance. Using the long reads generated by TELSeq, we identified numerous MGEs and cargo genes flanking the low-abundance ARGs, indicating that these ARGs could be transferred across bacterial taxa via HGT.

CONCLUSIONS

TELSeq can provide a nuanced view of the genomic context of microbial resistomes and thus has wide-ranging applications in public, animal, and human health, as well as environmental surveillance and monitoring of AMR. Thus, this technique represents a fundamental advancement for microbiome research and application. Video abstract.

Identification of a CCG-Enriched Expanded Allele in Patients with Myotonic Dystrophy Type 1 Using Amplification-Free Long-Read Sequencing

Myotonic dystrophy type 1 (DM1) exhibits highly heterogeneous clinical manifestations caused by an unstable CTG repeat expansion reaching up to 4000 CTG. The clinical variability depends on CTG repeat number, CNG repeat interruptions, and somatic mosaicism. Currently, none of these factors are simultaneously and accurately determined due to the limitations of gold standard methods used in clinical and research laboratories. An amplicon method for targeting the DMPK locus using single-molecule real-time sequencing was recently developed to accurately analyze expanded alleles. However, amplicon-based sequencing still depends on PCR, and the inherent bias toward preferential amplification of smaller repeats can be problematic in DM1. Thus, an amplification-free long-read sequencing method was developed by using CRISPR/Cas9 technology in DM1. This method was used to sequence the DMPK locus in patients with CTG repeat expansion ranging from 130 to >1000 CTG. We showed that elimination of PCR amplification improves the accuracy of measurement of inherited repeat number and somatic repeat variations, two key factors in DM1 severity and age at onset. For the first time, an expansion composed of >85% CCG repeats was identified by using this innovative method in a DM1 family with an atypical clinical profile. No-amplification targeted sequencing represents a promising method that can overcome research and diagnosis shortcomings, with translational implications for clinical and genetic counseling in DM1.

Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes

Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.

Functional characterization of prokaryotic dark matter: the road so far and what lies ahead

Eight-hundred thousand to one trillion prokaryotic species may inhabit our planet. Yet, fewer than two-hundred thousand prokaryotic species have been described. This uncharted fraction of microbial diversity, and its undisclosed coding potential, is known as the "microbial dark matter" (MDM). Next-generation sequencing has allowed to collect a massive amount of genome sequence data, leading to unprecedented advances in the field of genomics. Still, harnessing new functional information from the genomes of uncultured prokaryotes is often limited by standard classification methods. These methods often rely on sequence similarity searches against reference genomes from cultured species. This hinders the discovery of unique genetic elements that are missing from the cultivated realm. It also contributes to the accumulation of prokaryotic gene products of unknown function among public sequence data repositories, highlighting the need for new approaches for sequencing data analysis and classification. Increasing evidence indicates that these proteins of unknown function might be a treasure trove of biotechnological potential. Here, we outline the challenges, opportunities, and the potential hidden within the functional dark matter (FDM) of prokaryotes. We also discuss the pitfalls surrounding molecular and computational approaches currently used to probe these uncharted waters, and discuss future opportunities for research and applications.