Contamination-controlled high-throughput whole genome sequencing for influenza A viruses using the MiSeq sequencer

Recovering high-quality bacterial genomes from cross-contaminated cultures: a case study of marine Vibrio campbellii

BACKGROUND

Environmental monitoring of bacterial pathogens is critical for disease control in coastal marine ecosystems to maintain animal welfare and ecosystem function and to prevent significant economic losses. This requires accurate taxonomic identification of environmental bacterial pathogens, which often cannot be achieved by commonly used genetic markers (e.g., 16S rRNA gene), and an understanding of their pathogenic potential based on the information encoded in their genomes. The decreasing costs of whole genome sequencing (WGS), combined with newly developed bioinformatics tools, now make it possible to unravel the full potential of environmental pathogens, beyond traditional microbiological approaches. However, obtaining a high-quality bacterial genome, requires initial cultivation in an axenic culture, which is a bottleneck in environmental microbiology due to cross-contamination in the laboratory or isolation of non-axenic strains.

RESULTS

We applied WGS to determine the pathogenic potential of two Vibrio isolates from coastal seawater. During the analysis, we identified cross-contamination of one of the isolates and decided to use this dataset to evaluate the possibility of bioinformatic contaminant removal and recovery of bacterial genomes from a contaminated culture. Despite the contamination, using an appropriate bioinformatics workflow, we were able to obtain high quality and highly identical genomes (Average Nucleotide Identity value 99.98%) of one of the Vibrio isolates from both the axenic and the contaminated culture. Using the assembled genome, we were able to determine that this isolate belongs to a sub-lineage of Vibrio campbellii associated with several diseases in marine organisms. We also found that the genome of the isolate contains a novel Vibrio plasmid associated with bacterial defense mechanisms and horizontal gene transfer, which may offer a competitive advantage to this putative pathogen.

CONCLUSIONS

Our study shows that, using state-of-the-art bioinformatics tools and a sufficient sequencing effort, it is possible to obtain high quality genomes of the bacteria of interest and perform in-depth genomic analyses even in the case of a contaminated culture. With the new isolate and its complete genome, we are providing new insights into the genomic characteristics and functional potential of this sub-lineage of V. campbellii. The approach described here also highlights the possibility of recovering complete bacterial genomes in the case of non-axenic cultures or obligatory co-cultures.

Carryover Contamination-Controlled Amplicon Sequencing Workflow for Accurate Qualitative and Quantitative Detection of Pathogens: a Case Study on SARS-CoV-2

Carryover contamination during amplicon sequencing workflow (AMP-Seq) put the accuracy of the high-throughput detection for pathogens at risk. The purpose of this study is to develop a carryover contaminations-controlled AMP-Seq (ccAMP-Seq) workflow to enable accurate qualitative and quantitative detection for pathogens. By using the AMP-Seq workflow to detect SARS-CoV-2, Aerosols, reagents and pipettes were identified as potential sources of contaminations and ccAMP-Seq was then developed. ccAMP-Seq used filter tips and physically isolation of experimental steps to avoid cross contamination, synthetic DNA spike-ins to compete with contaminations and quantify SARS-CoV-2, dUTP/uracil DNA glycosylase system to digest the carryover contaminations, and a new data analysis procedure to remove the sequencing reads from contaminations. Compared to AMP-Seq, the contamination level of ccAMP-Seq was at least 22-folds lower and the detection limit was also about an order of magnitude lower-as low as one copy/reaction. By testing the dilution series of SARS-CoV-2 nucleic acid standard, ccAMP-Seq showed 100% sensitivity and specificity. The high sensitivity of ccAMP-Seq was further confirmed by the detection of SARS-CoV-2 from 62 clinical samples. The consistency between qPCR and ccAMP-Seq was 100% for all the 53 qPCR-positive clinical samples. Seven qPCR-negative clinical samples were found to be positive by ccAMP-Seq, which was confirmed by extra qPCR tests on subsequent samples from the same patients. This study presents a carryover contamination-controlled, accurate qualitative and quantitative amplicon sequencing workflow that addresses the critical problem of pathogen detection for infectious diseases. IMPORTANCE Accuracy, a key indicator of pathogen detection technology, is compromised by carryover contamination in the amplicon sequencing workflow. Taking the detection of SARS-CoV-2 as case, this study presents a new carryover contamination-controlled amplicon sequencing workflow. The new workflow significantly reduces the degree of contamination in the workflow, thereby significantly improving the accuracy and sensitivity of the SARS-CoV-2 detection and empowering the ability of quantitative detection. More importantly, the use of the new workflow is simple and economical. Therefore, the results of this study can be easily applied to other microorganism, which has great significance for improving the detection level of microorganism.

Novel Zoonotic Avian Influenza A(H3N8) Virus in Chicken, Hong Kong, China

Zoonotic and pandemic influenza continue to pose threats to global public health. Pandemics arise when novel influenza A viruses, derived in whole or in part from animal or avian influenza viruses, adapt to transmit efficiently in a human population that has little population immunity to contain its onward transmission. Viruses of previous pandemic concern, such as influenza A(H7N9), arose from influenza A(H9N2) viruses established in domestic poultry acquiring a hemagglutinin and neuraminidase from influenza A viruses of aquatic waterfowl. We report a novel influenza A(H3N8) virus in chicken that has emerged in a similar manner and that has been recently reported to cause zoonotic disease. Although they are H3 subtype, these avian viruses are antigenically distant from contemporary human influenza A(H3N2) viruses, and there is little cross-reactive immunity in the human population. It is essential to heighten surveillance for these avian A(H3N8) viruses in poultry and in humans.

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.

Next Generation Sequencing and Bioinformatics Methodologies for Infectious Disease Research and Public Health: Approaches, Applications, and Considerations for Development of Laboratory Capacity

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics.

Novel Multiplexed Amplicon-Based Sequencing to Quantify SARS-CoV-2 RNA from Wastewater

The application of wastewater-based epidemiology (WBE) to support the global response to the COVID-19 pandemic has shown encouraging outcomes. The accurate, sensitive, and high-throughput detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in municipal wastewater is critical for WBE. Here, we present a novel approach based on multiplexed amplicon-based sequencing, namely the ATOPlex platform, for detecting SARS-CoV-2. The ATOPlex platform is capable of quantifying SARS-CoV-2 RNA at concentrations that are at least 1 order of magnitude lower than the detection limit of reverse transcription quantitative polymerase chain reaction (RT-qPCR). Robust and accurate phylogenetic placement can be done at viral concentrations 4 times lower than the detection limit of RT-qPCR. We further found that the solid fraction in wastewater harbors a considerable amount of viral RNA, highlighting the need to extract viral RNA from the solid and liquid fractions of wastewater. This study delivers a highly sensitive, phylogenetically informative, and high-throughput analytical workflow that facilitates the application of WBE.

Method comparison of targeted influenza A virus typing and whole-genome sequencing from respiratory specimens of companion animals

Epidemics of H3N8 and H3N2 influenza A viruses (IAVs) in dogs, along with recognition of spillover infections from IAV strains typically found in humans or other animals, have emphasized the importance of efficient laboratory testing. Given the lack of active IAV surveillance or immunization requirements for dogs, cats, or horses imported into the United States, serotype prediction and whole-genome sequencing of positive specimens detected at veterinary diagnostic laboratories are also needed. The conserved sequences at the ends of the viral genome segments facilitate universal amplification of all segments of viral genomes directly from respiratory specimens. Although several methods for genomic analysis have been reported, no optimization focusing on companion animal strains has been described, to our knowledge. We compared 2 sets of published universal amplification primers using 26 IAV-positive specimens from dogs, horses, and a cat. Libraries prepared from the resulting amplicons were sequenced using Illumina chemistry, and reference-based assemblies were generated from the data produced by both methods. Although both methods produced high-quality data, coverage profiles and base calling differed between the 2 methods. The sequence data were also used to identify the subtype of the IAV strains sequenced and then compared to standard PCR assays for neuraminidase types N2 and N8.

Technical Validation of a Hepatitis C Virus Whole Genome Sequencing Assay for Detection of Genotype and Antiviral Resistance in the Clinical Pathway

Choice of direct acting antiviral (DAA) therapy for Hepatitis C Virus (HCV) in the United Kingdom and similar settings usually requires knowledge of the genotype and, in some cases, antiviral resistance (AVR) profile of the infecting virus. To determine these, most laboratories currently use Sanger technology, but next-generation sequencing (NGS) offers potential advantages in throughput and accuracy. However, NGS poses unique technical challenges, which require idiosyncratic development and technical validation approaches. This applies particularly to virology, where sequence diversity is high and the amount of starting genetic material is low, making it difficult to distinguish real data from artifacts. We describe the development and technical validation of a sequence capture-based HCV whole genome sequencing (WGS) assay to determine viral genotype and AVR profile. We use clinical samples of known subtypes and viral loads, and simulated FASTQ datasets to validate the analytical performances of both the wet laboratory and bioinformatic pipeline procedures. We show high concordance of the WGS assay compared to current "gold standard" Sanger assays. Specificity was 92.3 and 96.1% for AVR and genotyping, respectively. Discordances were due to the inability of Sanger assays to assign the correct subtype or accurately call mixed drug-resistant variants. We show high repeatability and reproducibility with >99.8% sequence similarity between sequence runs as well as high precision for variant frequency detection at >98.8% in the 95th percentile. Post-sequencing bioinformatics quality control workflows allow the accurate distinction between mixed infections, cross-contaminants and recombinant viruses at a threshold of >5% for the minority population. The sequence capture-based HCV WGS assay is more accurate than legacy AVR and genotyping assays. The assay has now been implemented in the clinical pathway of England's National Health Service HCV treatment programs, representing the first validated HCV WGS pipeline in clinical service. The data generated will additionally provide granular national-level genomic information for public health policy making and support the WHO HCV elimination strategy.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

The difference of detection rate of avian influenza virus in the wild bird surveillance using various methods

Korea is located within the East Asian-Australian flyway of wild migratory birds during the fall and winter seasons. Consequently, the likelihood of introduction of numerous subtypes and pathotypes of the Avian influenza (AI) virus to Korea has been thought to be very high. In the current study, we surveyed wild bird feces for the presence of AI virus that had been introduced to Korea between September 2017 and February 2018. To identify and characterize the AI virus, we employed commonly used methods, namely, virus isolation (VI) via egg inoculation, real-time reverse transcription-polymerase chain reaction (rRT-PCR), conventional RT-PCR (cRT-PCR) and a newly developed next generation sequencing (NGS) approach. In this study, 124 out of 11,145 fresh samples of wild migratory birds tested were rRT-PCR positive; only 52.0% of VI positive samples were determined as positive by rRT-PCR from fecal supernatant. Fifty AI virus specimens were isolated from fresh fecal samples and typed. The cRT-PCR subtyping results mostly coincided with the NGS results, although NGS detected the presence of 11 HA genes and four NA genes that were not detected by cRT-PCR. NGS analysis confirmed that 12% of the identified viruses were mixed-subtypes which were not detected by cRT-PCR. Prevention of the occurrence of AI virus requires a workflow for rapid and accurate virus detection and verification. However, conventional methods of detection have some limitations. Therefore, different methods should be combined for optimal surveillance, and further studies are needed in aspect of the introduction and application of new methods such as NGS.

Molecular insights into evolution, mutations and receptor-binding specificity of influenza A and B viruses from outpatients and hospitalized patients in Singapore

BACKGROUND

This study compared the genomes of influenza viruses that caused mild infections among outpatients and severe infections among hospitalized patients in Singapore, and characterized their molecular evolution and receptor-binding specificity.

METHODS

The complete genomes of influenza A/H1N1, A/H3N2 and B viruses that caused mild infections among outpatients and severe infections among inpatients in Singapore during 2012-2015 were sequenced and characterized. Using various bioinformatics approaches, we elucidated their evolutionary, mutational and structural patterns against the background of global and vaccine strains.

RESULTS

The phylogenetic trees of the 8 gene segments revealed that the outpatient and inpatient strains overlapped with representative global and vaccine strains. We observed a cluster of inpatients with A/H3N2 strains that were closely related to vaccine strain A/Texas/50/2012(H3N2). Several protein sites could accurately discriminate between outpatient versus inpatient strains, with site 221 in neuraminidase (NA) achieving the highest accuracy for A/H3N2. Interestingly, amino acid residues of inpatient but not outpatient isolates at those sites generally matched the corresponding residues in vaccine strains, except at site 145 of hemagglutinin (HA). This would be especially relevant for future surveillance of A/H3N2 strains in relation to their antigenicity and virulence. Furthermore, we observed a trend in which the HA proteins of influenza A/H3N2 and A/H1N1 exhibited enhanced ability to bind both avian and human host cell receptors. In contrast, the binding ability to each receptor was relatively stable for the HA of influenza B.

CONCLUSIONS

Overall, our findings extend our understanding of the molecular and structural evolution of influenza virus strains in Singapore within the global context of these dynamic viruses.

Divergent evolutionary trajectories of influenza B viruses underlie their contemporaneous epidemic activity

Two influenza B viruses (Victoria and Yamagata) cocirculate in humans and contribute to the estimated 290,000–650,000 annual influenza-attributed deaths. Here, we analysed influenza B genomic data to understand the causes of a recent surge in human influenza B infections. We found that evolution is acting differently on Yamagata and Victoria viruses and that this has led to the cocirculation of a diverse group of influenza B viruses. If this phenomenon continues, we could potentially witness the emergence of 3 or more distinct influenza B viruses that could require their own vaccine component, thereby complicating influenza vaccine formulation and highlighting the urgency of developing universal influenza vaccines.

Influenza B viruses have circulated in humans for over 80 y, causing a significant disease burden. Two antigenically distinct lineages (“B/Victoria/2/87-like” and “B/Yamagata/16/88-like,” termed Victoria and Yamagata) emerged in the 1970s and have cocirculated since 2001. Since 2015 both lineages have shown unusually high levels of epidemic activity, the reasons for which are unclear. By analyzing over 12,000 influenza B virus genomes, we describe the processes enabling the long-term success and recent resurgence of epidemics due to influenza B virus. We show that following prolonged diversification, both lineages underwent selective sweeps across the genome and have subsequently taken alternate evolutionary trajectories to exhibit epidemic dominance, with no reassortment between lineages. Hemagglutinin deletion variants emerged concomitantly in multiple Victoria virus clades and persisted through epistatic mutations and interclade reassortment—a phenomenon previously only observed in the 1970s when Victoria and Yamagata lineages emerged. For Yamagata viruses, antigenic drift of neuraminidase was a major driver of epidemic activity, indicating that neuraminidase-based vaccines and cross-reactivity assays should be employed to monitor and develop robust protection against influenza B morbidity and mortality. Overall, we show that long-term diversification and infrequent selective sweeps, coupled with the reemergence of hemagglutinin deletion variants and antigenic drift of neuraminidase, are factors that contributed to successful circulation of diverse influenza B clades. Further divergence of hemagglutinin variants with poor cross-reactivity could potentially lead to circulation of 3 or more distinct influenza B viruses, further complicating influenza vaccine formulation and highlighting the urgent need for universal influenza vaccines.

Improved detection of influenza A virus from blue-winged teals by sequencing directly from swab material

Abstract

The greatest diversity of influenza A virus (IAV) is found in wild aquatic birds of the orders Anseriformes and Charadriiformes. In these birds, IAV replication occurs mostly in the intestinal tract. Fecal, cloacal, and/or tracheal swabs are typically collected and tested by real-time RT-PCR (rRT-PCR) and/or by virus isolation in embryonated chicken eggs in order to determine the presence of IAV. Virus isolation may impose bottlenecks that select variant populations that are different from those circulating in nature, and such bottlenecks may result in artifactual representation of subtype diversity and/or underrepresented mixed infections. The advent of next-generation sequencing (NGS) technologies provides an opportunity to explore to what extent IAV subtype diversity is affected by virus isolation in eggs. In the present work, we evaluated the advantage of sequencing by NGS directly from swab material of IAV rRT-PCR-positive swabs collected during the 2013-14 surveillance season in Guatemala and compared to results from NGS after virus isolation. The results highlight the benefit of sequencing IAV genomes directly from swabs to better understand subtype diversity and detection of alternative amino acid motifs that could otherwise escape detection using traditional methods of virus isolation. In addition, NGS sequencing data from swabs revealed reduced presence of defective interfering particles compared to virus isolates. We propose an alternative workflow in which original swab samples positive for IAV by rRT-PCR are first subjected to NGS before attempting viral isolation. This approach should speed the processing of samples and better capture natural IAV diversity.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.3h2n106.

Seasonal Genetic Drift of Human Influenza A Virus Quasispecies Revealed by Deep Sequencing

After a pandemic wave in 2009 following their introduction in the human population, the H1N1pdm09 viruses replaced the previously circulating, pre-pandemic H1N1 virus and, along with H3N2 viruses, are now responsible for the seasonal influenza type A epidemics. So far, the evolutionary potential of influenza viruses has been mainly documented by consensus sequencing data. However, like other RNA viruses, influenza A viruses exist as a population of diverse, albeit related, viruses, or quasispecies. Interest in this quasispecies nature has increased with the development of next generation sequencing (NGS) technologies that allow a more in-depth study of the genetic variability. NGS deep sequencing methodologies were applied to determine the whole genome genetic heterogeneity of the three categories of influenza A viruses that circulated in humans between 2007 and 2012 in France, directly from clinical respiratory specimens. Mutation frequencies and single nucleotide polymorphisms were used for comparisons to address the level of natural intrinsic heterogeneity of influenza A viruses. Clear differences in single nucleotide polymorphism profiles between seasons for a given subtype also revealed the constant genetic drift that human influenza A virus quasispecies undergo.

Universal Target Capture of HIV Sequences From NGS Libraries

Background: Global surveillance of viral sequence diversity is needed to keep pace with the constant evolution of HIV. Recent next generation sequencing (NGS) methods have realized the goal of sequencing circulating virus directly from patient specimens. Yet, a simple, universal approach that maximizes sensitivity and sequencing capacity remains elusive. Here we present a novel HIV enrichment strategy to yield near complete genomes from low viral load specimens. Methodology: A non-redundant biotin-labeled probe set (HIV-xGen; n = 652) was synthesized to tile all HIV-1 (groups M, N, O, and P) and HIV-2 (A and B) strains. Illumina Nextera barcoded libraries of either gene-specific or randomly primed cDNA derived from infected plasma were hybridized to probes in a single pool and unbound sequences were washed away. Captured viral cDNA was amplified by Illumina adaptor primers, sequenced on a MiSeq, and NGS reads were demultiplexed for alignment with CLC Bio software. Results: HIV-xGen probes selectively captured and amplified reads spanning the entirety of the HIV phylogenetic tree. HIV sequences clearly present in unenriched libraries of specimens but previously not observed due to high host background levels, insufficient sequencing depth or the extent of multiplexing, were now enriched by >1,000-fold. Thus, xGen selection not only substantially increased the depth of existing sequence, but also extended overall genome coverage by an average of 40%. We characterized 50 new, diverse HIV strains from clinical specimens and demonstrated a viral load cutoff of approximately log 3.5 copies/ml for full length coverage. Genome coverage was <20% for 5/10 samples with viral loads <log 3.5 copies/ml and >90% for 35/40 samples with higher viral loads. Conclusions: Characterization of >20 complete genomes at a time is now possible from a single probe hybridization and MiSeq run. With the versatility to capture all HIV strains and the sensitivity to detect low titer specimens, HIV-xGen will serve as an important tool for monitoring HIV sequence diversity.

INSaFLU: an automated open web-based bioinformatics suite "from-reads" for influenza whole-genome-sequencing-based surveillance

BACKGROUND

A new era of flu surveillance has already started based on the genetic characterization and exploration of influenza virus evolution at whole-genome scale. Although this has been prioritized by national and international health authorities, the demanded technological transition to whole-genome sequencing (WGS)-based flu surveillance has been particularly delayed by the lack of bioinformatics infrastructures and/or expertise to deal with primary next-generation sequencing (NGS) data.

RESULTS

We developed and implemented INSaFLU ("INSide the FLU"), which is the first influenza-oriented bioinformatics free web-based suite that deals with primary NGS data (reads) towards the automatic generation of the output data that are actually the core first-line "genetic requests" for effective and timely influenza laboratory surveillance (e.g., type and sub-type, gene and whole-genome consensus sequences, variants' annotation, alignments and phylogenetic trees). By handling NGS data collected from any amplicon-based schema, the implemented pipeline enables any laboratory to perform multi-step software intensive analyses in a user-friendly manner without previous advanced training in bioinformatics. INSaFLU gives access to user-restricted sample databases and projects management, being a transparent and flexible tool specifically designed to automatically update project outputs as more samples are uploaded. Data integration is thus cumulative and scalable, fitting the need for a continuous epidemiological surveillance during the flu epidemics. Multiple outputs are provided in nomenclature-stable and standardized formats that can be explored in situ or through multiple compatible downstream applications for fine-tuned data analysis. This platform additionally flags samples as "putative mixed infections" if the population admixture enrolls influenza viruses with clearly distinct genetic backgrounds, and enriches the traditional "consensus-based" influenza genetic characterization with relevant data on influenza sub-population diversification through a depth analysis of intra-patient minor variants. This dual approach is expected to strengthen our ability not only to detect the emergence of antigenic and drug resistance variants but also to decode alternative pathways of influenza evolution and to unveil intricate routes of transmission.

CONCLUSIONS

In summary, INSaFLU supplies public health laboratories and influenza researchers with an open "one size fits all" framework, potentiating the operationalization of a harmonized multi-country WGS-based surveillance for influenza virus. INSaFLU can be accessed through https://insaflu.insa.pt .

Sample tracking in microbiome community profiling assays using synthetic 16S rRNA gene spike-in controls

Workflows for microbiome community profiling by high-throughput sequencing are prone to sample mix-ups and cross-contamination due to the complexity of the procedures and large number of samples typically analyzed in parallel. We employed synthetic 16S rRNA gene spike-in controls to establish a method for tracking of sample identity and detection of cross-contamination in microbiome community profiling assays based on 16S rRNA gene amplicon sequencing (16S-seq). Results demonstrated that combinatorial sample tracking mixes (STMs) can be reliably resolved by Illumina sequencing and faithfully represent their sample of origin. In a single-blinded experiment, addition of STMs at low levels was shown to be sufficient to unambiguously identify and resolve swapped samples. Using artificial admixtures of individually SMT-tagged samples, we further established the ability to detect and quantify cross-contamination down to a level of approximately 1%. The utility of our technique was underscored through detection of an unplanned case of cross-contamination that occurred during this study. By enabling detection of sample mix-ups and cross-contamination throughout 16S-seq workflows, the present technique thus assures provenance of sequence data on a per-sample basis. The method can be readily implemented in standard 16S-seq workflows and its routine application is expected to enhance the reliability of 16S-seq data.

Increased sequencing depth does not increase captured diversity of arbuscular mycorrhizal fungi

The arrival of 454 sequencing represented a major breakthrough by allowing deeper sequencing of environmental samples than was possible with existing Sanger approaches. Illumina MiSeq provides a further increase in sequencing depth but shorter read length compared with 454 sequencing. We explored whether Illumina sequencing improves estimates of arbuscular mycorrhizal (AM) fungal richness in plant root samples, compared with 454 sequencing. We identified AM fungi in root samples by sequencing amplicons of the SSU rRNA gene with 454 and Illumina MiSeq paired-end sequencing. In addition, we sequenced metagenomic DNA without prior PCR amplification. Amplicon-based Illumina sequencing yielded two orders of magnitude higher sequencing depth per sample than 454 sequencing. Initial analysis with minimal quality control recorded five times higher AM fungal richness per sample with Illumina sequencing. Additional quality control of Illumina samples, including restriction of the marker region to the most variable amplicon fragment, revealed AM fungal richness values close to those produced by 454 sequencing. Furthermore, AM fungal richness estimates were not correlated with sequencing depth between 300 and 30,000 reads per sample, suggesting that the lower end of this range is sufficient for adequate description of AM fungal communities. By contrast, metagenomic Illumina sequencing yielded very few AM fungal reads and taxa and was dominated by plant DNA, suggesting that AM fungal DNA is present at prohibitively low abundance in colonised root samples. In conclusion, Illumina MiSeq sequencing yielded higher sequencing depth, but similar richness of AM fungi in root samples, compared with 454 sequencing.