Index hopping on the Illumina HiseqX platform and its consequences for ancient DNA studies

Scalable genotyping of microbial colonies

The sequence of the 16S region is taxonomically informative and widely used for genotyping microbes. While it is easy and inexpensive to genotype several isolates by Sanger sequencing the 16S region, this method becomes quite costly if scaled to many isolates. High-throughput sequencing provides one potential avenue for obtaining 16S sequences at scale but presents additional challenges. First, DNA purification workflows for high-throughput sample preparation are labour-intensive and expensive. Second, cost-effective multiplexing and library preparation schemes are difficult to implement for many libraries on a single sequencing run. Therefore, we implemented a scalable protocol for isolate genotyping involving colony polymerase chain reaction (PCR) with simple cell lysis as well as a four-barcode indexing scheme that enables scalable multiplexing and streamlined library preparation by amplifying with four primers simultaneously in a single reaction. We tested this protocol on 93 colonies cultured from environmental samples, and we were able to ascertain the identity of ~90% of microbial isolates.

Sequencing RNA from old, dried specimens reveals past viromes and properties of long-surviving RNA

Recovery of virus sequences from old samples provides an opportunity to study virus evolution and reconstruct historic virus-host interactions. Studies of old virus sequences have mainly relied on DNA or on RNA from fixed or frozen samples. The millions of specimens in natural history museums represent a potential treasure trove of old virus sequences, but it is not clear how well RNA survives in old samples. We experimentally assessed the stability of RNA in insects stored dry at room temperature over 72 weeks. Although RNA molecules grew fragmented, RNA yields remained surprisingly constant. RT-qPCR of host and virus RNA showed minimal differences between dried and frozen specimens. To assess RNA survival in much older samples we acquired Drosophila specimens from North American entomological collections. We recovered sequences from known and novel viruses including several coding complete virus genomes from a fly collected in 1908. We found that the virome of D. melanogaster has changed little over the past century. Galbut virus, the most prevalent virus infection in contemporary D. melanogaster, was also the most common in historic samples. Finally, we investigated the genomic and physical features of surviving RNA. RNA that survived was fragmented, chemically damaged, and preferentially double stranded or contained in ribonucleoprotein complexes. This showed that RNA - especially certain types of RNA - can survive in biological specimens over extended periods in the absence of fixation or freezing and confirms the utility of dried specimens to provide a clearer understanding of virus evolution.

Contributions of Hanseniaspora species to Pinot Noir microbial terroir in Oregon's Willamette Valley wine region

The apiculate yeast genus Hanseniaspora has appeared frequently in enological research for more than 100 years, mostly focused upon the species H. uvarum due to its notable capacity to cause spoilage. Recently, there has been increased research into the potential benefits of other Hanseniaspora species, such as H. vineae, in producing more complex wines. Furthermore, large-scale DNA sequencing-based (metabarcoding) vineyard ecology studies have suggested that Hanseniaspora species may not be evenly distributed. To address potential differences across geographical areas in Oregon, we sampled extensively from 12 vineyards within the Willamette Valley American Viticultural Area (AVA), across 2 sub-AVAs (Eola-Amity Hills and Yamhill-Carlton). Metabarcoding was then used to assess the contribution of Hanseniaspora to the grape berry fungal community and the impact of wine processing on diversity. While 6 of the 23 recognized Hanseniaspora species were present on Pinot Noir grapes in the Willamette Valley AVA, differences between vineyards were driven by the abundance of H. uvarum. Significant positive correlations between the amount of H. uvarum present in must and at cold soak, and then cold soak to early ferment were observed. While intuitive, it is worth noting that no prior studies have observed this across such a large number of grape samples from different vineyards. Our results provide clear evidence that the abundance of H. uvarum on grapes may be an important predictor of potential impacts on wine quality, particularly if performing cold soak, which acts as an enrichment step.

IMPORTANCE

Hanseniaspora yeasts are frequently found in uninoculated wine fermentations, and depending upon the species present, their contributions to the wine may be positive or negative. We found that in Oregon's Willamette Valley, the most common species of Hanseniaspora in Pinot Noir vineyards was the known spoilage organism, H. uvarum. This species was one of the strongest contributors to differences in fungal communities between different vineyards and was enriched during typical Pinot Noir processing. These results support Hanseniaspora as an integral and functional component of vineyard "microbial terroir" within Oregon.

Comparative Performance Evaluation of Double-Stranded RNA High-Throughput Sequencing for the Detection of Viral Infection in Temperate Fruit Crops

There is limited information on the compared performances of biological, serological. and molecular assays with high-throughput sequencing (HTS) for viral indexing in temperate fruit crops. Here, using a range of samples of predetermined virological status, we compared two performance criteria (inclusivity and analytical sensitivity) of enzyme-linked immunosorbent assay (ELISA), molecular hybridization, reverse transcription (RT)-PCR, and double-stranded RNA (dsRNA) HTS for the detection of a total of 14 viruses (10 genera) and four viroids (three genera). When undiluted samples from individual plants were used, ELISA had the lowest performance, with an overall detection rate of 68.7%, followed by RT-PCR (82.5%) and HTS (90.7%; 100% if considering only viruses). The lower performance of RT-PCR reflected the inability to amplify some isolates as a consequence of point mutations affecting primer-binding sites. In addition, HTS identified viruses that had not been identified by other assays in nearly two-thirds of the samples. Analysis of serial dilutions of fruit tree samples allowed comparison of analytical sensitivities for various viruses. ELISA showed the lowest analytical sensitivity, but RT-PCR showed higher analytical sensitivity than HTS for most of the samples. Overall, these results confirm the superiority of HTS over biological indexing in terms of speed and inclusivity and show that while the absolute analytical sensitivity of RT-PCR tends to be higher than that of HTS, PCR inclusivity is affected by viral genetic diversity. Taken together, these results make a strong case for the implementation of HTS-based approaches in fruit tree viral testing protocols supporting quarantine and certification programs.

Pinpointing the microbiota of tardigrades: What is really there?

Microbiota are considered significant in the biology of tardigrades, yet their diversity and distribution remain largely unexplored. This is partly due to the methodological challenges associated with studying the microbiota of small organisms that inhabit microbe-rich environments. In our study, we characterized the microbiota of 31 species of cultured tardigrades using 16S rRNA amplicon sequencing. We employed various sample preparation strategies and multiple types of controls and estimated the number of microbes in samples using synthetic DNA spike-ins. We also reanalysed data from previous tardigrade microbiome studies. Our findings suggest that the microbial communities of cultured tardigrades are predominantly composed of bacterial genotypes originating from food, medium, or reagents. Despite numerous experiments, we found it challenging to identify strains that were enriched in certain tardigrades, which would have indicated likely symbiotic associations. Putative tardigrade-associated microbes rarely constituted more than 20% of the datasets, although some matched symbionts identified in other studies. We also uncovered serious contamination issues in previous tardigrade microbiome studies, casting doubt on some of their conclusions. We concluded that tardigrades are not universally dependent on specialized microbes. Our work underscores the need for rigorous safeguards in studies of the microbiota of microscopic organisms and serves as a cautionary tale for studies involving samples with low microbiome abundance.

Defining a metagenomic threshold for detecting low abundances of Providencia alcalifaciens in canine faecal samples

Introduction

Acute haemorrhagic diarrhoea syndrome (AHDS) in dogs is a condition of unknown aetiology. Providencia alcalifaciens is suspected to play a role in the disease as it was commonly found in dogs suffering from AHDS during a Norwegian outbreak in 2019. The role of this bacterium as a constituent of the canine gut microbiota is unknown, hence this study set out to investigate its occurrence in healthy dogs using metagenomics.

Materials and methods

To decrease the likelihood of false detection, we established a metagenomic threshold for P. alcalifaciens by spiking culture-negative stool samples with a range of bacterial dilutions and analysing these by qPCR and shotgun metagenomics. The detection limit for P. alcalifaciens was determined and used to establish a metagenomic threshold. The threshold was validated on naturally contaminated faecal samples with known cultivation status for P. alcalifaciens. Finally, the metagenomic threshold was used to determine the occurrence of P. alcalifaciens in shotgun metagenomic datasets from canine faecal samples (n=362) collected in the HUNT One Health project.

Results

The metagenomic assay and qPCR had a detection limit of 1.1x103 CFU P. alcalifaciens per faecal sample, which corresponded to a Cq value of 31.4 and 569 unique k-mer counts by shotgun metagenomics. Applying this metagenomic threshold to 362 faecal metagenomic datasets from healthy dogs, P. alcalifaciens was found in only 1.1% (95% CI [0.0, 6.8]) of the samples, and then in low relative abundances (median: 0.04%; range: 0.00 to 0.81%). The sensitivity of the qPCR and shotgun metagenomics assay was low, as only 40% of culture-positive samples were also positive by qPCR and metagenomics.

Discussion

Using our detection limit, the occurrence of P. alcalifaciens in faecal samples from healthy dogs was low. Given the low sensitivity of the metagenomic assay, these results do not rule out a significantly higher occurrence of this bacterium at a lower abundance.

Rapid and sensitive detection of genome contamination at scale with FCS-GX

Assembled genome sequences are being generated at an exponential rate. Here we present FCS-GX, part of NCBI's Foreign Contamination Screen (FCS) tool suite, optimized to identify and remove contaminant sequences in new genomes. FCS-GX screens most genomes in 0.1-10 min. Testing FCS-GX on artificially fragmented genomes demonstrates high sensitivity and specificity for diverse contaminant species. We used FCS-GX to screen 1.6 million GenBank assemblies and identified 36.8 Gbp of contamination, comprising 0.16% of total bases, with half from 161 assemblies. We updated assemblies in NCBI RefSeq to reduce detected contamination to 0.01% of bases. FCS-GX is available at https://github.com/ncbi/fcs/ or https://doi.org/10.5281/zenodo.10651084 .

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Evaluation of next-generation sequencing performance for in vitro detection of viruses in biological products

Next-Generation Sequencing (NGS) can detect nucleic acid sequences in a massively parallel sequencing. This technology is expected to be widely applied for the detection of viral contamination in biologics. The recently published ICH-Q5A (R2) draft indicates that NGS could be an alternative or supplement to in vitro viral tests. To examine the performance of NGS for the in vitro detection of viruses, adenovirus type 5 (Ad5), a model virus, was inoculated into Vero cells, which are the most popular indicator cells for the detection of adventitious viruses in the in vitro test. Total RNA extracted from the Vero cells infected with Ad5 was serially diluted with that from non-infected Vero cells, and each sample was analyzed using short- or long-read NGSs. The limits of detection of both NGS methods were almost the same and both methods were sensitive enough to detect viral sequences as long as there was at least one copy in one assay. Although the multiplexing in NGS carries the risk of cross-contamination among the samples, which could lead to false positives, this technology has the potential to become a rapid and sensitive method for detecting adventitious agents in biologics.

Clinical evaluation of cell-free and cellular metagenomic next-generation sequencing of infected body fluids

INTRODUCTION

Previous studies have evaluated metagenomic next-generation sequencing (mNGS) of cell-free DNA (cfDNA) for pathogen detection in blood and body fluid samples. However, no study has assessed the diagnostic efficacy of mNGS using cellular DNA.

OBJECTIVES

This is the first study to systematically evaluate the efficacy of cfDNA and cellular DNA mNGS for pathogen detection.

METHODS

A panel of seven microorganisms was used to compare cfDNA and cellular DNA mNGS assays concerning limits of detection (LoD), linearity, robustness to interference, and precision. In total, 248 specimens were collected between December 2020 and December 2021. The medical records of all the patients were reviewed. These specimens were analysed using cfDNA and cellular DNA mNGS assays, and the mNGS results were confirmed using viral qPCR, 16S rRNA, and internal transcribed spacer (ITS) amplicon next-generation sequencing.

RESULTS

The LoD of cfDNA and cellular DNA mNGS was 9.3 to 149 genome equivalents (GE)/mL and 27 to 466 colony-forming units (CFU)/mL, respectively. The intra- and inter-assay reproducibility of cfDNA and cellular DNA mNGS was 100%. Clinical evaluation revealed that cfDNA mNGS was good at detecting the virus in blood samples (receiver operating characteristic (ROC) area under the curve (AUC), 0.9814). In contrast, the performance of cellular DNA mNGS was better than that of cfDNA mNGS in high host background samples. Overall, the diagnostic efficacy of cfDNA combined with cellular DNA mNGS (ROC AUC, 0.8583) was higher than that of cfDNA (ROC AUC, 0.8041) or cellular DNA alone (ROC AUC, 0.7545).

CONCLUSION

Overall, cfDNA mNGS is good for detecting viruses, and cellular DNA mNGS is suitable for high host background samples. The diagnostic efficacy was higher when cfDNA and cellular DNA mNGS were combined.

Benchmarking of virome metagenomic analysis approaches using a large, 60+ members, viral synthetic community

IMPORTANCE

We report here efforts to benchmark performance of two widespread approaches for virome analysis, which target either virion-associated nucleic acids (VANA) or highly purified double-stranded RNAs (dsRNAs). This was achieved using synthetic communities of varying complexity levels, up to a highly complex community of 72 viral agents (115 viral molecules) comprising isolates from 21 families and 61 genera of plant viruses. The results obtained confirm that the dsRNA-based approach provides a more complete representation of the RNA virome, in particular, for high complexity ones. However, for viromes of low to medium complexity, VANA appears a reasonable alternative and would be the preferred choice if analysis of DNA viruses is of importance. Several parameters impacting performance were identified as well as a direct relationship between the completeness of virome description and sample sequencing depth. The strategy, results, and tools used here should prove useful in a range of virome analysis efforts.

An integrative protocol for one-step PCR amplicon library construction and accurate demultiplexing of pooled sequencing data

High-throughput sequencing of amplicons has been widely used to precisely and efficiently identify species compositions and analyze community structures, greatly promoting biological studies involving large amounts of complex samples, especially those involving environmental and pathogen-monitoring ones. Commercial library preparation kits for amplicon sequencing, which generally require multiple steps, including adapter ligation and indexing, are expensive and time-consuming, especially for applications at a large scale. To overcome these limitations, a "one-step PCR approach" has been previously proposed for constructions of amplicon libraries using long fusion primers. However, efficient amplifications of target genes and accurate demultiplexing of pooled sequencing data remain to be addressed. To tackle these, we present an integrative protocol for one-step PCR amplicon library construction (OSPALC). High-quality reads have been generated by this approach to reliably identify species compositions of mock bacterial communities and environmental samples. With this protocol, the amplicon library is constructed through one regular PCR with long primers, and the total cost per DNA/cDNA sample decreases to just 7% of the typical cost via the multi-step PCR approach. Empirically tested primers and optimized PCR conditions to construct OSPALC libraries for 16S rDNA V4 regions are demonstrated as a case study. Tools to design primers targeting at any genomic regions are also presented. In principle, OSPALC can be readily applied to construct amplicon libraries of any target genes using DNA or RNA samples, and will facilitate research in numerous fields.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00182-1.

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

Holocene deglaciation drove rapid genetic diversification of Atlantic walrus

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

Rapid and easy construction of a simplified amplicon sequencing (simplified AmpSeq) library for marker-assisted selection

Marker-assisted selection (MAS) is fundamental for plant breeding programs, as it can identify desirable seedlings at a young stage and reduce the cost, time and space needed for plant maintenance, especially for perennial crops. To facilitate the process of genotyping, which is time consuming and laborious, we developed a simplified amplicon sequencing (simplified AmpSeq) library construction method for next-generation sequencing that can be applied to MAS in breeding programs. The method is based on one-step PCR with a mixture of two primer sets: the first consisting of tailed target primers, the second of primers that contain flow-cell binding sites, indexes and tail sequences complementary to those in the first set. To demonstrate the process of MAS using s implified AmpSeq, we created databases of genotypes for important traits by using cultivar collections including triploid cultivars and segregating seedlings of Japanese pear (Pyrus pyrifolia Nakai), Japanese chestnut (Castanea crenata Sieb. et Zucc.) and apple (Malus domestica Borkh.). Simplified AmpSeq has the advantages of high repeatability, ability to estimate allele number in polyploid species and semi-automatic evaluation using target allele frequencies. Because this method provides high flexibility for designing primer sets and targeting any variant, it will be useful for plant breeding programs.

Rapid and sensitive detection of genome contamination at scale with FCS-GX

Assembled genome sequences are being generated at an exponential rate. Here we present FCS-GX, part of NCBI's Foreign Contamination Screen (FCS) tool suite, optimized to identify and remove contaminant sequences in new genomes. FCS-GX screens most genomes in 0.1-10 minutes. Testing FCS-GX on artificially fragmented genomes demonstrates sensitivity >95% for diverse contaminant species and specificity >99.93%. We used FCS-GX to screen 1.6 million GenBank assemblies and identified 36.8 Gbp of contamination (0.16% of total bases), with half from 161 assemblies. We updated assemblies in NCBI RefSeq to reduce detected contamination to 0.01% of bases. FCS-GX is available at https://github.com/ncbi/fcs/.

Multi-factorial examination of amplicon sequencing workflows from sample preparation to bioinformatic analysis

BACKGROUND

The development of sequencing technologies to evaluate bacterial microbiota composition has allowed new insights into the importance of microbial ecology. However, the variety of methodologies used among amplicon sequencing workflows leads to uncertainty about best practices as well as reproducibility and replicability among microbiome studies. Using a bacterial mock community composed of 37 soil isolates, we performed a comprehensive methodological evaluation of workflows, each with a different combination of methodological factors spanning sample preparation to bioinformatic analysis to define sources of artifacts that affect coverage, accuracy, and biases in the resulting compositional profiles.

RESULTS

Of the workflows examined, those using the V4-V4 primer set enabled the highest level of concordance between the original mock community and resulting microbiome sequence composition. Use of a high-fidelity polymerase, or a lower-fidelity polymerase with an increased PCR elongation time, limited chimera formation. Bioinformatic pipelines presented a trade-off between the fraction of distinct community members identified (coverage) and fraction of correct sequences (accuracy). DADA2 and QIIME2 assembled V4-V4 reads amplified by Taq polymerase resulted in the highest accuracy (100%) but had a coverage of only 52%. Using mothur to assemble and denoise V4-V4 reads resulted in a coverage of 75%, albeit with marginally lower accuracy (99.5%).

CONCLUSIONS

Optimization of microbiome workflows is critical for accuracy and to support reproducibility and replicability among microbiome studies. These considerations will help reveal the guiding principles of microbial ecology and impact the translation of microbiome research to human and environmental health.

DNA High-Throughput Sequencing for Arthropod Gut Content Analysis to Evaluate Effectiveness and Safety of Biological Control Agents

The search for effective biological control agents without harmful non-target effects has been constrained by the use of impractical (field direct observation) or imprecise (cage experiments) methods. While advances in the DNA sequencing methods, more specifically the development of high-throughput sequencing (HTS), have been quickly incorporated in biodiversity surveys, they have been slow to be adopted to determine arthropod prey range, predation rate and food web structure, and critical information to evaluate the effectiveness and safety of a biological control agent candidate. The lack of knowledge on how HTS methods could be applied by ecological entomologists constitutes part of the problem, although the lack of expertise and the high cost of the analysis also are important limiting factors. In this review, we describe how the latest HTS methods of metabarcoding and Lazaro, a method to identify prey by mapping unassembled shotgun reads, can serve biological control research, showing both their power and limitations. We explain how they work to determine prey range and also how their data can be used to estimate predation rates and subsequently be translated into food webs of natural enemy and prey populations helping to elucidate their role in the community. We present a brief history of prey detection through molecular gut content analysis and also the attempts to develop a more precise formula to estimate predation rates, a problem that still remains. We focused on arthropods in agricultural ecosystems, but most of what is covered here can be applied to natural systems and non-arthropod biological control candidates as well.

Ecology, Not Host Phylogeny, Shapes the Oral Microbiome in Closely Related Species

Host-associated microbiomes are essential for a multitude of biological processes. Placed at the contact zone between external and internal environments, the little-studied oral microbiome has important roles in host physiology and health. Here, we investigate the roles of host evolutionary relationships and ecology in shaping the oral microbiome in three closely related gorilla subspecies (mountain, Grauer's, and western lowland gorillas) using shotgun metagenomics of 46 museum-preserved dental calculus samples. We find that the oral microbiomes of mountain gorillas are functionally and taxonomically distinct from the other two subspecies, despite close evolutionary relationships and geographic proximity with Grauer's gorillas. Grauer's gorillas show intermediate bacterial taxonomic and functional, and dietary profiles. Altitudinal differences in gorilla subspecies ranges appear to explain these patterns, suggesting a close connection between dental calculus microbiomes and the environment, likely mediated through diet. This is further supported by the presence of gorilla subspecies-specific phyllosphere/rhizosphere taxa in the oral microbiome. Mountain gorillas show a high abundance of nitrate-reducing oral taxa, which may promote adaptation to a high-altitude lifestyle by modulating blood pressure. Our results suggest that ecology, rather than evolutionary relationships and geographic distribution, shape the oral microbiome in these closely related species.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol

Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

Assessing the performance of historical skins and bones for museomics using wolf specimens as a case study

Advances in the field of museomics have promoted a high sampling demand for natural history collections (NHCs), eventually resulting in damage to invaluable resources to understand historical biodiversity. It is thus essential to achieve a consensus about which historical tissues present the best sources of DNA. In this study, we evaluated the performance of different historical tissues from Iberian wolf NHCs in genome-wide assessments. We targeted three tissues—bone (jaw and femur), maxilloturbinal bone, and skin—that have been favored by traditional taxidermy practices for mammalian carnivores. Specifically, we performed shotgun sequencing and target capture enrichment for 100,000 single nucleotide polymorphisms (SNPs) selected from the commercial Canine HD BeadChip across 103 specimens from 1912 to 2005. The performance of the different tissues was assessed using metrics based on endogenous DNA content, uniquely high-quality mapped reads after capture, and enrichment proportions. All samples succeeded as DNA sources, regardless of their collection year or sample type. Skin samples yielded significantly higher amounts of endogenous DNA compared to both bone types, which yielded equivalent amounts. There was no evidence for a direct effect of tissue type on capture efficiency; however, the number of genotyped SNPs was strictly associated with the starting amount of endogenous DNA. Evaluation of genotyping accuracy for distinct minimum read depths across tissue types showed a consistent overall low genotyping error rate (<7%), even at low (3x) coverage. We recommend the use of skins as reliable and minimally destructive sources of endogenous DNA for whole-genome and target enrichment approaches in mammalian carnivores. In addition, we provide a new 100,000 SNP capture array validated for historical DNA (hDNA) compatible to the Canine HD BeadChip for high-quality DNA. The increasing demand for NHCs as DNA sources should encourage the generation of genomic datasets comparable among studies.

A 14,000-year-old genome sheds light on the evolution and extinction of a Pleistocene vulture

The New World Vulture [Coragyps] occidentalis (L. Miller, 1909) is one of many species that were extinct by the end of the Pleistocene. To understand its evolutionary history we sequenced the genome of a 14,000 year old [Coragyps] occidentalis found associated with megaherbivores in the Peruvian Andes. occidentalis has been viewed as the ancestor, or possibly sister, to the extant Black Vulture Coragyps atratus, but genomic data shows occidentalis to be deeply nested within the South American clade of atratus. Coragyps atratus inhabits lowlands, but the fossil record indicates that occidentalis mostly occupied high elevations. Our results suggest that occidentalis evolved from a population of atratus in southwestern South America that colonized the High Andes 300 to 400 kya. The morphological and morphometric differences between occidentalis and atratus may thus be explained by ecological diversification following from the natural selection imposed by this new and extreme, high elevation environment. The sudden evolution of a population with significantly larger body size and different anatomical proportions than atratus thus constitutes an example of punctuated evolution.

14,000 year old DNA reveals the evolutionary dynamics and adaptations of South American vultures.

Detecting and removing sample contamination in phylogenomic data: an example and its implications for Cicadidae phylogeny (Insecta: Hemiptera)

Contamination of a genetic sample with DNA from one or more non-target species is a continuing concern of molecular phylogenetic studies, both Sanger sequencing studies and Next-Generation Sequencing (NGS) studies. We developed an automated pipeline for identifying and excluding likely cross-contaminated loci based on detection of bimodal distributions of patristic distances across gene trees. When the contamination occurs between samples within a dataset, comparisons between a contaminated sample and its contaminant taxon will yield bimodal distributions with one peak close to zero patristic distance. This new method does not rely on a priori knowledge of taxon relatedness nor does it determine the causes(s) of the contamination. Exclusion of putatively contaminated loci from a dataset generated for the insect family Cicadidae showed that these sequences were affecting some topological patterns and branch supports, although the effects were sometimes subtle, with some contamination-influenced relationships exhibiting strong bootstrap support. Long tip branches and outlier values for one anchored phylogenomic pipeline statistic (AvgNHomologs) were correlated with the presence of contamination. While the AHE markers used here, which target hemipteroid taxa, proved effective in resolving deep and shallow level Cicadidae relationships in aggregate, individual markers contained inadequate phylogenetic signal, in part probably due to short length. The cleaned dataset, consisting of 429 loci, from 90 genera representing 44 of 56 current Cicadidae tribes, supported three of the four sampled Cicadidae subfamilies in concatenated-matrix maximum likelihood (ML) and multispecies coalescent-based species tree analyses, with the fourth subfamily weakly supported in the ML trees. No well-supported patterns from previous family-level Sanger sequencing studies of Cicadidae phylogeny were contradicted. One taxon (Aragualna plenalinea) did not fall with its current subfamily in the genetic tree, and this genus and its tribe Aragualnini is reclassified to Tibicininae following morphological re-examination. Only subtle differences were observed in trees after removal of loci for which divergent base frequencies were detected. Greater success may be achieved by increased taxon sampling and developing a probe set targeting a more recent common ancestor and longer loci. Searches for contamination are an essential step in phylogenomic analyses of all kinds and our pipeline is an effective solution.

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Library Preparation and Sequencing Platform Introduce Bias in Metagenomic-Based Characterizations of Microbiomes

Metagenomics is increasingly used to describe microbial communities in biological specimens. Ideally, the steps involved in the processing of the biological specimens should not change the microbiome composition in a way that it could lead to false interpretations of inferred microbial community composition. Common steps in sample preparation include sample collection, storage, DNA isolation, library preparation, and DNA sequencing. Here, we assess the effect of three library preparation kits and two DNA sequencing platforms. Of the library preparation kits, one involved a PCR step (Nextera), and two were PCR free (NEXTflex and KAPA). We sequenced the libraries on Illumina HiSeq and NextSeq platforms. As example microbiomes, two pig fecal samples and two sewage samples of which aliquots were stored at different storage conditions (immediate processing and storage at -80°C) were assessed. All DNA isolations were performed in duplicate, totaling 80 samples, excluding controls. We found that both library preparation and sequencing platform had systematic effects on the inferred microbial community composition. The different sequencing platforms introduced more variation than library preparation and freezing the samples. The results highlight that all sample processing steps need to be considered when comparing studies. Standardization of sample processing is key to generating comparable data within a study, and comparisons of differently generated data, such as in a meta-analysis, should be performed cautiously. IMPORTANCE Previous research has reported effects of sample storage conditions and DNA isolation procedures on metagenomics-based microbiome composition; however, the effect of library preparation and DNA sequencing in metagenomics has not been thoroughly assessed. Here, we provide evidence that library preparation and sequencing platform introduce systematic biases in the metagenomic-based characterization of microbial communities. These findings suggest that library preparation and sequencing are important parameters to keep consistent when aiming to detect small changes in microbiome community structure. Overall, we recommend that all samples in a microbiome study are processed in the same way to limit unwanted variations that could lead to false conclusions. Furthermore, if we are to obtain a more holistic insight from microbiome data generated around the world, we will need to provide more detailed sample metadata, including information about the different sample processing procedures, together with the DNA sequencing data at the public repositories.

Dual UMIs and Dual Barcodes With Minimal PCR Amplification Removes Artifacts and Acquires Accurate Antibody Repertoire

Antibody repertoire sequencing (Rep-seq) has been widely used to reveal repertoire dynamics and to interrogate antibodies of interest at single nucleotide-level resolution. However, polymerase chain reaction (PCR) amplification introduces extensive artifacts including chimeras and nucleotide errors, leading to false discovery of antibodies and incorrect assessment of somatic hypermutations (SHMs) which subsequently mislead downstream investigations. Here, a novel approach named DUMPArts, which improves the accuracy of antibody repertoires by labeling each sample with dual barcodes and each molecule with dual unique molecular identifiers (UMIs) via minimal PCR amplification to remove artifacts, is developed. Tested by ultra-deep Rep-seq data, DUMPArts removed inter-sample chimeras, which cause artifactual shared clones and constitute approximately 15% of reads in the library, as well as intra-sample chimeras with erroneous SHMs and constituting approximately 20% of the reads, and corrected base errors and amplification biases by consensus building. The removal of these artifacts will provide an accurate assessment of antibody repertoires and benefit related studies, especially mAb discovery and antibody-guided vaccine design.

How challenging RADseq data turned out to favor coalescent-based species tree inference. A case study in Aichryson (Crassulaceae)

Analysing multiple genomic regions while incorporating detection and qualification of discordance among regions has become standard for understanding phylogenetic relationships. In plants, which usually have comparatively large genomes, this is feasible by the combination of reduced-representation library (RRL) methods and high-throughput sequencing enabling the cost effective acquisition of genomic data for thousands of loci from hundreds of samples. One popular RRL method is RADseq. A major disadvantage of established RADseq approaches is the rather short fragment and sequencing range, leading to loci of little individual phylogenetic information. This issue hampers the application of coalescent-based species tree inference. The modified RADseq protocol presented here targets ca. 5,000 loci of 300-600nt length, sequenced with the latest short-read-sequencing (SRS) technology, has the potential to overcome this drawback. To illustrate the advantages of this approach we use the study group Aichryson Webb & Berthelott (Crassulaceae), a plant genus that diversified on the Canary Islands. The data analysis approach used here aims at a careful quality control of the long loci dataset. It involves an informed selection of thresholds for accurate clustering, a thorough exploration of locus properties, such as locus length, coverage and variability, to identify potential biased data and a comparative phylogenetic inference of filtered datasets, accompanied by an evaluation of resulting BS support, gene and site concordance factor values, to improve overall resolution of the resulting phylogenetic trees. The final dataset contains variable loci with an average length of 373nt and facilitates species tree estimation using a coalescent-based summary approach. Additional improvements brought by the approach are critically discussed.

Contamination Issue in Viral Metagenomics: Problems, Solutions, and Clinical Perspectives

We describe the most common internal and external sources and types of contamination encountered in viral metagenomic studies and discuss their negative impact on sequencing results, particularly for low-biomass samples and clinical applications. We also propose some basic recommendations for reducing the background noise in viral shotgun metagenomic (SM) studies, which would limit the bias introduced by various classes of contaminants. Regardless of the specific viral SM protocol, contamination cannot be totally avoided; in particular, the issue of reagent contamination should always be addressed with high priority. There is an urgent need for the development and validation of standards for viral metagenomic studies especially if viral SM protocols will be more widely applied in diagnostics.

Revisiting the effect of PCR replication and sequencing depth on biodiversity metrics in environmental DNA metabarcoding

Environmental DNA (eDNA) metabarcoding is an increasingly popular tool for measuring and cataloguing biodiversity. Because the environments and substrates in which DNA is preserved differ considerably, eDNA research often requires bespoke approaches to generating eDNA data. Here, we explore how two experimental choices in eDNA study design-the number of PCR replicates and the depth of sequencing of PCR replicates-influence the composition and consistency of taxa recovered from eDNA extracts. We perform 24 PCR replicates from each of six soil samples using two of the most common metabarcodes for Fungi and Viridiplantae (ITS1 and ITS2), and sequence each replicate to an average depth of ~84,000 reads. We find that PCR replicates are broadly consistent in composition and relative abundance of dominant taxa, but that low abundance taxa are often unique to one or a few PCR replicates. Taxa observed in one out of 24 PCR replicates make up 21-29% of the total taxa detected. We also observe that sequencing depth or rarefaction influences alpha diversity and beta diversity estimates. Read sampling depth influences local contribution to beta diversity, placement in ordinations, and beta dispersion in ordinations. Our results suggest that, because common taxa drive some alpha diversity estimates, few PCR replicates and low read sampling depths may be sufficient for many biological applications of eDNA metabarcoding. However, because rare taxa are recovered stochastically, eDNA metabarcoding may never fully recover the true amplifiable alpha diversity in an eDNA extract. Rare taxa drive PCR replicate outliers of alpha and beta diversity and lead to dispersion differences at different read sampling depths. We conclude that researchers should consider the complexity and unevenness of a community when choosing analytical approaches, read sampling depths, and filtering thresholds to arrive at stable estimates.

The oral microbiota of wild bears in Sweden reflects the history of antibiotic use by humans

Following the advent of industrial-scale antibiotic production in the 1940s,¹ antimicrobial resistance (AMR) has been on the rise and now poses a major global health threat in terms of mortality, morbidity, and economic burden.^2,3 Because AMR can be exchanged between humans, livestock, and wildlife, wild animals can be used as indicators of human-associated AMR contamination of the environment.⁴ However, AMR is a normal function of natural environments and is present in host-associated microbiomes, which makes it challenging to distinguish between anthropogenic and natural sources.^4,5 One way to overcome this difficulty is to use historical samples that span the period from before the mass production of antibiotics to today. We used shotgun metagenomic sequencing of dental calculus, the calcified form of the oral microbial biofilm, to determine the abundance and repertoire of AMR genes in the oral microbiome of Swedish brown bears collected over the last 180 years. Our temporal metagenomics approach allowed us to establish a baseline of natural AMR in the pre-antibiotics era and to quantify a significant increase in total AMR load and diversity of AMR genes that is consistent with patterns of national human antibiotic use. We also demonstrated a significant decrease in total AMR load in bears in the last two decades, which coincides with Swedish strategies to mitigate AMR. Our study suggests that public health policies can be effective in limiting human-associated AMR contamination of the environment and wildlife.

Concatenation of paired-end reads improves taxonomic classification of amplicons for profiling microbial communities

Background

Taxonomic classification of genetic markers for microbiome analysis is affected by the numerous choices made from sample preparation to bioinformatics analysis. Paired-end read merging is routinely used to capture the entire amplicon sequence when the read ends overlap. However, the exclusion of unmerged reads from further analysis can result in underestimating the diversity in the sequenced microbial community and is influenced by bioinformatic processes such as read trimming and the choice of reference database. A potential solution to overcome this is to concatenate (join) reads that do not overlap and keep them for taxonomic classification. The use of concatenated reads can outperform taxonomic recovery from single-end reads, but it remains unclear how their performance compares to merged reads. Using various sequenced mock communities with different amplicons, read length, read depth, taxonomic composition, and sequence quality, we tested how merging and concatenating reads performed for genus recall and precision in bioinformatic pipelines combining different parameters for read trimming and taxonomic classification using different reference databases.

Results

The addition of concatenated reads to merged reads always increased pipeline performance. The top two performing pipelines both included read concatenation, with variable strengths depending on the mock community. The pipeline that combined merged and concatenated reads that were quality-trimmed performed best for mock communities with larger amplicons and higher average quality sequences. The pipeline that used length-trimmed concatenated reads outperformed quality trimming in mock communities with lower quality sequences but lost a significant amount of input sequences for taxonomic classification during processing. Genus level classification was more accurate using the SILVA reference database compared to Greengenes.

Conclusions

Merged sequences with the addition of concatenated sequences that were unable to be merged increased performance of taxonomic classifications. This was especially beneficial in mock communities with larger amplicons. We have shown for the first time, using an in-depth comparison of pipelines containing merged vs concatenated reads combined with different trimming parameters and reference databases, the potential advantages of concatenating sequences in improving resolution in microbiome investigations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04410-2.

Handling of targeted amplicon sequencing data focusing on index hopping and demultiplexing using a nested metabarcoding approach in ecology

High-throughput sequencing platforms are increasingly being used for targeted amplicon sequencing because they enable cost-effective sequencing of large sample sets. For meaningful interpretation of targeted amplicon sequencing data and comparison between studies, it is critical that bioinformatic analyses do not introduce artefacts and rely on detailed protocols to ensure that all methods are properly performed and documented. The analysis of large sample sets and the use of predefined indexes create challenges, such as adjusting the sequencing depth across samples and taking sequencing errors or index hopping into account. However, the potential biases these factors introduce to high-throughput amplicon sequencing data sets and how they may be overcome have rarely been addressed. On the example of a nested metabarcoding analysis of 1920 carabid beetle regurgitates to assess plant feeding, we investigated: (i) the variation in sequencing depth of individually tagged samples and the effect of library preparation on the data output; (ii) the influence of sequencing errors within index regions and its consequences for demultiplexing; and (iii) the effect of index hopping. Our results demonstrate that despite library quantification, large variation in read counts and sequencing depth occurred among samples and that the sequencing error rate in bioinformatic software is essential for accurate adapter/primer trimming and demultiplexing. Moreover, setting an index hopping threshold to avoid incorrect assignment of samples is highly recommended.

Isothermal amplifications - a comprehensive review on current methods

The introduction of nucleic acid amplification techniques has revolutionized the field of medical diagnostics in the last decade. The advent of PCR catalyzed the increasing application of DNA, not just for molecular cloning but also for molecular based diagnostics. Since the introduction of PCR, a deeper understanding of molecular mechanisms and enzymes involved in DNA/RNA replication has spurred the development of novel methods devoid of temperature cycling. Isothermal amplification methods have since been introduced utilizing different mechanisms, enzymes, and conditions. The ease with which isothermal amplification methods have allowed nucleic acid amplification to be carried out has had a profound impact on the way molecular diagnostics are being designed after the turn of the millennium. With all the advantages isothermal amplification brings, the issues or complications surrounding each method are heterogeneous making it difficult to identify the best approach for an end-user. This review pays special attention to the various isothermal amplification methods by classifying them based on the mechanistic characteristics which include reaction formats, amplification information, promoter, strand break, and refolding mechanisms. We would also compare the efficiencies and usefulness of each method while highlighting the potential applications and detection methods involved. This review will serve as an overall outlook on the journey and development of isothermal amplification methods as a whole.

Characterizing Microbiomes via Sequencing of Marker Loci: Techniques To Improve Throughput, Account for Cross-Contamination, and Reduce Cost

New approaches to characterizing microbiomes via high-throughput sequencing provide impressive gains in efficiency and cost reduction compared to approaches that were standard just a few years ago. However, the speed of method development has been such that staying abreast of the latest technological advances is challenging. Moreover, shifting laboratory protocols to include new methods can be expensive and time consuming. To facilitate adoption of new techniques, we provide a guide and review of recent advances that are relevant for single-locus sequence-based study of microbiomes—from extraction to library preparation—including a primer regarding the use of liquid-handling automation in small-scale academic settings. Additionally, we describe several amendments to published techniques to improve throughput, track contamination, and reduce cost. Notably, we suggest adding synthetic DNA molecules to each sample during nucleic acid extraction, thus providing a method of documenting incidences of cross-contamination. We also describe a dual-indexing scheme for Illumina sequencers that allows multiplexing of many thousands of samples with minimal PhiX input. Collectively, the techniques that we describe demonstrate that laboratory technology need not impose strict limitations on the scale of molecular microbial ecology studies.

IMPORTANCE New methods to characterize microbiomes reduce technology-imposed limitations to study design, but many new approaches have not been widely adopted. Here, we present techniques to increase throughput and reduce contamination alongside a thorough review of current best practices.

The Assembled and Annotated Genome of the Fairy-Ring Fungus Marasmius oreades

Marasmius oreades is a basidiomycete fungus that grows in so called "fairy rings," which are circular, underground mycelia common in lawns across temperate areas of the world. Fairy rings can be thought of as natural, long-term evolutionary experiments. As each ring has a common origin and expands radially outwards over many years, different sectors will independently accumulate mutations during growth. The genotype can be followed to the next generation, as mushrooms producing the sexual spores are formed seasonally at the edge of the ring. Here, we present new genomic data from 95 single-spore isolates of the species, which we used to construct a genetic linkage map and an updated version of the genome assembly. The 44-Mb assembly was anchored to 11 linkage groups, producing chromosome-length scaffolds. Gene annotation revealed 13,891 genes, 55% of which contained a pfam domain. The repetitive fraction of the genome was 22%, and dominated by retrotransposons and DNA elements of the KDZ and Plavaka groups. The level of assembly contiguity we present is so far rare in mushroom-forming fungi, and we expect studies of genomics, transposons, phylogenetics, and evolution to be facilitated by the data we present here of the iconic fairy-ring mushroom.

Microbiome Structure of the Aphid Myzus persicae (Sulzer) Is Shaped by Different Solanaceae Plant Diets

Myzus persicae (Sulzer) is an important insect pest in agriculture that has a very broad host range. Previous research has shown that the microbiota of insects has implications for their growth, development, and environmental adaptation. So far, there is little detailed knowledge about the factors that influence and shape the microbiota of aphids. In the present study, we aimed to investigate diet-induced changes in the microbiome of M. persicae using high-throughput sequencing of bacterial 16S ribosomal RNA gene fragments in combination with molecular and microbiological experiments. The transfer of aphids to different plants from the Solanaceae family resulted in a substantial decrease in the abundance of the primary symbiont Buchnera. In parallel, a substantial increase in the abundance of Pseudomonas was observed; it accounted for up to 69.4% of the bacterial community in M. persicae guts and the attached bacteriocytes. In addition, we observed negative effects on aphid population dynamics when they were transferred to pepper plants (Capsicum annuum L.). The microbiome of this treatment group showed a significantly lower increase in the abundance of Pseudomonas when compared with the other Solanaceae plant diets, which might be related to the adaptability of the host to this diet. Molecular quantifications of bacterial genera that were substantially affected by the different diets were implemented as an additional verification of the microbiome-based observations. Complementary experiments with bacteria isolated from aphids that were fed with different plants indicated that nicotine-tolerant strains occur in Solanaceae-fed specimens, but they were not restricted to them. Overall, our mechanistic approach conducted under controlled conditions provided strong indications that the aphid microbiome shows responses to different plant diets. This knowledge could be used in the future to develop environmentally friendly methods for the control of insect pests in agriculture.

Handling of spurious sequences affects the outcome of high-throughput 16S rRNA gene amplicon profiling

16S rRNA gene amplicon sequencing is a popular approach for studying microbiomes. However, some basic concepts have still not been investigated comprehensively. We studied the occurrence of spurious sequences using defined microbial communities based on data either from the literature or generated in three sequencing facilities and analyzed via both operational taxonomic units (OTUs) and amplicon sequence variants (ASVs) approaches. OTU clustering and singleton removal, a commonly used approach, delivered approximately 50% (mock communities) to 80% (gnotobiotic mice) spurious taxa. The fraction of spurious taxa was generally lower based on ASV analysis, but varied depending on the gene region targeted and the barcoding system used. A relative abundance of 0.25% was found as an effective threshold below which the analysis of spurious taxa can be prevented to a large extent in both OTU- and ASV-based analysis approaches. Using this cutoff improved the reproducibility of analysis, i.e., variation in richness estimates was reduced by 38% compared with singleton filtering using six human fecal samples across seven sequencing runs. Beta-diversity analysis of human fecal communities was markedly affected by both the filtering strategy and the type of phylogenetic distances used for comparison, highlighting the importance of carefully analyzing data before drawing conclusions on microbiome changes. In summary, handling of artifact sequences during bioinformatic processing of 16S rRNA gene amplicon data requires careful attention to avoid the generation of misleading findings. We propose the concept of effective richness to facilitate the comparison of alpha-diversity across studies.

Diagnostic Validation of a Clinical Laboratory-Oriented Targeted RNA Sequencing System for Detecting Gene Fusions in Hematologic Malignancies

Targeted RNA sequencing (RNA-seq) is a highly accurate method for sequencing transcripts of interest with a high resolution and throughput. However, RNA-seq has not been widely performed in clinical molecular laboratories because of the complexity of data processing and interpretation. We developed and validated a customized RNA-seq panel and data processing protocol for fusion detection using 4 analytical validation samples and 51 clinical samples, covering seven types of hematologic malignancies. Analytical validation showed that the results for target gene coverage and between- and within-run precision and linearity tests were reliable. Using clinical samples, RNA-seq based on filtering and prioritization strategies detected all 25 known fusions previously found by multiplex reverse transcriptase-PCR and fluorescence in situ hybridization. It also detected nine novel fusions. Known fusions detected by RNA-seq included two IGH rearrangements supported by expression analysis. Novel fusions included six that targeted just one partner gene. In addition, 18 disease- and drug resistance-associated transcript variants in ABL1, GATA2, IKZF1, JAK2, RUNX1, and WT1 were designated simultaneously. Expression analysis showed distinct clustering according to subtype and lineage. In conclusion, this study showed that our customized RNA-seq system had a reliable and stable performance for fusion detection, with enhanced diagnostic yield for hematologic malignancies in a clinical diagnostic setting.

Multiplexed detection of SARS-CoV-2 and other respiratory infections in high throughput by SARSeq

The COVID-19 pandemic has demonstrated the need for massively-parallel, cost-effective tests monitoring viral spread. Here we present SARSeq, saliva analysis by RNA sequencing, a method to detect SARS-CoV-2 and other respiratory viruses on tens of thousands of samples in parallel. SARSeq relies on next generation sequencing of multiple amplicons generated in a multiplexed RT-PCR reaction. Two-dimensional, unique dual indexing, using four indices per sample, enables unambiguous and scalable assignment of reads to individual samples. We calibrate SARSeq on SARS-CoV-2 synthetic RNA, virions, and hundreds of human samples of various types. Robustness and sensitivity were virtually identical to quantitative RT-PCR. Double-blinded benchmarking to gold standard quantitative-RT-PCR performed by human diagnostics laboratories confirms this high sensitivity. SARSeq can be used to detect Influenza A and B viruses and human rhinovirus in parallel, and can be expanded for detection of other pathogens. Thus, SARSeq is ideally suited for differential diagnostic of infections during a pandemic.

Massively parallel but cost-effective testing is essential to monitor the spread of pathogenic agents. Here the authors present SARSseq, which uses a dual indexing strategy in a multiplexed RT-PCR reaction to diagnose SARS-CoV-2 at scale.

Predicting sample success for large-scale ancient DNA studies on marine mammals

In recent years, nonhuman ancient DNA studies have begun to focus on larger sample sizes and whole genomes, offering the potential to reveal exciting and hitherto unknown answers to ongoing biological and archaeological questions. However, one major limitation to such studies is the substantial financial and time investments still required during sample screening, due to uncertainty regarding successful sample selection. This study investigates the effect of a wide range of sample properties including latitude, sample age, skeletal element, collagen preservation, and context on endogenous content and DNA damage profiles for 317 ancient and historic pinniped samples collected from across the North Atlantic and surrounding regions. Using generalised linear and mixed-effect models, we found that a range of factors affected DNA preservation within each of the species under consideration. The most important findings were that endogenous content varied significantly within species according to context, the type of skeletal element, the collagen content and collection year. There also appears to be an effect of the sample's geographic origin, with samples from the Arctic generally showing higher endogenous content and lower damage rates. Both latitude and sample age were found to have significant relationships with damage levels, but only for walrus samples. Sex, ontogenetic age and extraction material preparation were not found to have any significant relationship with DNA preservation. Overall, skeletal element and sample context were found to be the most influential factors and should therefore be considered when selecting samples for large-scale ancient genome studies.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

Leaf-associated fungal and viral communities of wild plant populations differ between cultivated and natural ecosystems

Plants are colonized by diverse fungal and viral communities that influence their growth and survival as well as ecosystem functioning. Viruses interact with both plants and the fungi they host. Our understanding of plant-fungi-virus interactions is very limited, especially in wild plants. Combining metagenomic and culturomic approaches, we assessed the richness, diversity, and composition of leaf-associated fungal and viral communities from pools of herbaceous wild plants representative of four sites corresponding to cultivated or natural ecosystems. We identified 161 fungal families and 18 viral families comprising 249 RNA-dependent RNA polymerase-based operational taxonomic units (RdRp OTUs) from leaves. Fungal culturomics captured 12.3% of the fungal diversity recovered with metagenomic approaches and, unexpectedly, retrieved viral OTUs that were almost entirely different from those recovered by leaf metagenomics. Ecosystem management had a significant influence on both leaf mycobiome and virome, with a higher fungal community richness in natural ecosystems and a higher viral family richness in cultivated ecosystems, suggesting that leaf-associated fungal and viral communities are under the influence of different ecological drivers. Both the leaf-associated fungal and viral community compositions showed a strong site-specificity. Further research is needed to confirm these trends and unravel the factors structuring plant-fungi-virus interactions in wild plant populations.

Bioinformatics resources for SARS-CoV-2 discovery and surveillance

In early January 2020, the novel coronavirus (SARS-CoV-2) responsible for a pneumonia outbreak in Wuhan, China, was identified using next-generation sequencing (NGS) and readily available bioinformatics pipelines. In addition to virus discovery, these NGS technologies and bioinformatics resources are currently being employed for ongoing genomic surveillance of SARS-CoV-2 worldwide, tracking its spread, evolution and patterns of variation on a global scale. In this review, we summarize the bioinformatics resources used for the discovery and surveillance of SARS-CoV-2. We also discuss the advantages and disadvantages of these bioinformatics resources and highlight areas where additional technical developments are urgently needed. Solutions to these problems will be beneficial not only to the prevention and control of the current COVID-19 pandemic but also to infectious disease outbreaks of the future.

Biogeography and Independent Diversification in the Protist Symbiont Community of Heterotermes tenuis

The eukaryotic microbiome of “lower” termites is highly stable and host-specific. This is due to the mutually obligate nature of the symbiosis and the direct inheritance of protists by proctodeal trophallaxis. However, vertical transmission is occasionally imperfect, resulting in daughter colonies that lack one or more of the expected protist species. This phenomenon could conceivably lead to regional differences in protist community composition within a host species. Here, we have characterized the protist symbiont community of Heterotermes tenuis (Hagen) (Blattodea: Rhinotermitidae) from samples spanning South and Central America. Using light microscopy, single cell isolation, and amplicon sequencing, we report eight species-level protist phylotypes belonging to four genera in the phylum Parabasalia. The diversity and distribution of each phylotype’s 18S rRNA amplicon sequence variants (ASVs) mostly did not correlate with geographical or host genetic distances according to Mantel tests, consistent with the lack of correlation we observed between host genetic and geographical distances. However, the ASV distances of Holomastigotoides Ht3 were significantly correlated with geography while those of Holomastigotoides Ht1 were significantly correlated with host phylogeny. These results suggest mechanisms by which termite-associated protist species may diversify independently of each other and of their hosts, shedding light on the coevolutionary dynamics of this important symbiosis.

Swab-Seq: A high-throughput platform for massively scaled up SARS-CoV-2 testing

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is due to the high rates of transmission by individuals who are asymptomatic at the time of transmission^1,2. Frequent, widespread testing of the asymptomatic population for SARS-CoV-2 is essential to suppress viral transmission. Despite increases in testing capacity, multiple challenges remain in deploying traditional reverse transcription and quantitative PCR (RT-qPCR) tests at the scale required for population screening of asymptomatic individuals. We have developed SwabSeq, a high-throughput testing platform for SARS-CoV-2 that uses next-generation sequencing as a readout. SwabSeq employs sample-specific molecular barcodes to enable thousands of samples to be combined and simultaneously analyzed for the presence or absence of SARS-CoV-2 in a single run. Importantly, SwabSeq incorporates an in vitro RNA standard that mimics the viral amplicon, but can be distinguished by sequencing. This standard allows for end-point rather than quantitative PCR, improves quantitation, reduces requirements for automation and sample-to-sample normalization, enables purification-free detection, and gives better ability to call true negatives. After setting up SwabSeq in a high-complexity CLIA laboratory, we performed more than 80,000 tests for COVID-19 in less than two months, confirming in a real world setting that SwabSeq inexpensively delivers highly sensitive and specific results at scale, with a turn-around of less than 24 hours. Our clinical laboratory uses SwabSeq to test both nasal and saliva samples without RNA extraction, while maintaining analytical sensitivity comparable to or better than traditional RT-qPCR tests. Moving forward, SwabSeq can rapidly scale up testing to mitigate devastating spread of novel pathogens.