Application of different DNA extraction procedures, library preparation protocols and sequencing platforms: impact on sequencing results

Transcriptomic and genetic profiling in a spontaneous non-human primate model of hypertrophic cardiomyopathy and sudden cardiac death

Hypertrophic cardiomyopathy (HCM) afflicts humans, cats, pigs, and rhesus macaques. Disease sequelae include congestive heart failure, thromboembolism, and sudden cardiac death (SCD). Sarcomeric mutations explain some human and cat cases, however, the molecular basis in rhesus macaques remains unknown. RNA-Seq of the LV tissues of five HCM-affected and seven healthy control rhesus macaques was employed for differential transcriptomic analyses. DNA from 15 severely HCM-affected and 21 healthy geriatric rhesus macaques were selected for whole-genome sequencing. A genome-wide association study (GWAS) of disease status and SCD outcome was performed. 614 down- and 1,065 upregulated differentially expressed genes (DEGs) were identified between groups. The top DEG (MAFF) was overexpressed in affected animals (log2FoldChange = 4.71; PAdjusted-value = 1.14E-133). Channelopathy-associated enriched terms were identified in ~ 57% of downregulated DEGs providing transcriptomic evidence of hypertrophic and arrhythmic disease processes. For GWAS, no putative variant withstood segregation. Polygenic modeling analysis resulted in poor prediction power and burden testing could not explain HCM by an association of multiple variants in any gene. Neither single nor compound genetic variant(s), or identified polygenic profile, suggest complex genotype-phenotype interactions in rhesus macaques. Brought forth is an established dataset of robustly phenotyped rhesus macaques as an open-access resource for future cardiovascular disease genetic studies.

Effect of alkaline leaching of phosphogypsum on sulfate reduction activity and bacterial community composition using different sources of anaerobic microbial inoculum

Phosphogypsum (PG), a by-product of the phosphate industry, is high in sulfate, (SO42-), which makes it an excellent substrate for sulfate-reducing bacteria (SRB) to produce hydrogen sulfide. This work aimed to optimize SO42- leaching from PG to achieve a high biological reduction of SO42- and generate high sulfide concentrations for subsequent use in the biological recovery of elemental sulfur. Five SRB consortia were isolated and enriched from: IS (Industrial sludges), MS (Marine sediments), WC (Winogradsky column), SNV (petroleum industry sediments) and PG (stored Phosphogypsum). The five consortia showed reduction activity when using PG leachate (with water) as source of SO42- and lactate, acetate, or glucose as the electron donor. The highest reduction rate (81.5 %) was registered using lactate and the IS consortium (81.5 %) followed by MS (79 %) and PG (71 %). To enhance the concentration of leached SO42- from PG for future utilization with the isolated consortia, PG was treated with NaOH solutions (2 % and 5 %). SO42- release of 97 % was achieved with a 5 % concentration and the resulting leachate was further diluted to target a SO42- concentration of 12.4 g·L-1 for utilization with the isolated consortia. Compared to water leachate, a significantly higher reduction rate was registered (2 g·L-1 of SO42) using the IS consortium, demonstrating limited inhibition effect of sulfide- concentration on SRB functionalities. Moreover, metagenomic analysis of the consortia revealed that using PG as a source of SO42- increased the abundance of Deltaproteobacteria, including known SRB like Desulfovibrio, Desulfomicrobium, and Desulfosporosinus, as well as novel SRB genera (Cupidesulfovibrio, Desulfocurvus, Desulfococcus) that showed, for the first time, significant potential as novel sulfate-reducers using PG as a SO42- source.

Results of the 2020 Genomic Proficiency Test for the network of European Union Reference Laboratory for Antimicrobial Resistance assessing whole-genome-sequencing capacities

The global surveillance and outbreak investigation of antimicrobial resistance (AMR) is amidst a paradigm shift from traditional biology to bioinformatics. This is due to developments in whole-genome-sequencing (WGS) technologies, bioinformatics tools, and reduced costs. The increased use of WGS is accompanied by challenges such as standardization, quality control (QC), and data sharing. Thus, there is global need for inter-laboratory WGS proficiency test (PT) schemes to evaluate laboratories' capacity to produce reliable genomic data. Here, we present the results of the first iteration of the Genomic PT (GPT) organized by the Global Capacity Building Group at the Technical University of Denmark in 2020. Participating laboratories sequenced two isolates and corresponding DNA of Salmonella enterica, Escherichia coli and Campylobacter coli, using WGS methodologies routinely employed at their laboratories. The participants' ability to obtain consistently good-quality WGS data was assessed based on several QC WGS metrics. A total of 21 laboratories from 21 European countries submitted WGS and meta-data. Most delivered high-quality sequence data with only two laboratories identified as overall underperforming. The QC metrics, N50 and number of contigs, were identified as good indicators for high-sequencing quality. We propose QC thresholds for N50 greater than 20 000 and 25 000 for Campylobacter coli and Escherichia coli, respectively, and number of contigs >200 bp greater than 225, 265 and 100 for Salmonella enterica, Escherichia coli and Campylobacter coli, respectively. The GPT2020 results confirm the importance of systematic QC procedures, ensuring the submission of reliable WGS data for surveillance and outbreak investigation to meet the requirements of the paradigm shift in methodology.

Three Distinct Annotation Platforms Differ in Detection of Antimicrobial Resistance Genes in Long-Read, Short-Read, and Hybrid Sequences Derived from Total Genomic DNA or from Purified Plasmid DNA

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

In vitro and in silico parameters for precise cgMLST typing of Listeria monocytogenes

Background

Whole genome sequencing analyzed by core genome multi-locus sequence typing (cgMLST) is widely used in surveillance of the pathogenic bacteria Listeria monocytogenes. Given the heterogeneity of available bioinformatics tools to define cgMLST alleles, our aim was to identify parameters influencing the precision of cgMLST profiles.

Methods

We used three L. monocytogenes reference genomes from different phylogenetic lineages and assessed the impact of in vitro (i.e. tested genomes, successive platings, replicates of DNA extraction and sequencing) and in silico parameters (i.e. targeted depth of coverage, depth of coverage, breadth of coverage, assembly metrics, cgMLST workflows, cgMLST completeness) on cgMLST precision made of 1748 core loci. Six cgMLST workflows were tested, comprising assembly-based (BIGSdb, INNUENDO, GENPAT, SeqSphere and BioNumerics) and assembly-free (i.e. kmer-based MentaLiST) allele callers. Principal component analyses and generalized linear models were used to identify the most impactful parameters on cgMLST precision.

Results

The isolate’s genetic background, cgMLST workflows, cgMLST completeness, as well as depth and breadth of coverage were the parameters that impacted most on cgMLST precision (i.e. identical alleles against reference circular genomes). All workflows performed well at ≥40X of depth of coverage, with high loci detection (> 99.54% for all, except for BioNumerics with 97.78%) and showed consistent cluster definitions using the reference cut-off of ≤7 allele differences.

Conclusions

This highlights that bioinformatics workflows dedicated to cgMLST allele calling are largely robust when paired-end reads are of high quality and when the sequencing depth is ≥40X.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08437-4.

Whole Genome Sequencing Provides an Added Value to the Investigation of Staphylococcal Food Poisoning Outbreaks

Through staphylococcal enterotoxin (SE) production, Staphylococcus aureus is a common cause of food poisoning. Detection of staphylococcal food poisoning (SFP) is mostly performed using immunoassays, which, however, only detect five of 27 SEs described to date. Polymerase chain reactions are, therefore, frequently used in complement to identify a bigger arsenal of SE at the gene level (se) but are labor-intensive. Complete se profiling of isolates from different sources, i.e., food and human cases, is, however, important to provide an indication of their potential link within foodborne outbreak investigation. In addition to complete se gene profiling, relatedness between isolates is determined with more certainty using pulsed-field gel electrophoresis, Staphylococcus protein A gene typing and other methods, but these are shown to lack resolution. We evaluated how whole genome sequencing (WGS) can offer a solution to these shortcomings. By WGS analysis of a selection of S. aureus isolates, including some belonging to a confirmed foodborne outbreak, its added value as the ultimate multiplexing method was demonstrated. In contrast to PCR-based se gene detection for which primers are sometimes shown to be non-specific, WGS enabled complete se gene profiling with high performance, provided that a database containing reference sequences for all se genes was constructed and employed. The custom compiled database and applied parameters were made publicly available in an online user-friendly interface. As an all-in-one approach with high resolution, WGS additionally allowed inferring correct isolate relationships. The different DNA extraction kits that were tested affected neither se gene profiling nor relatedness determination, which is interesting for data sharing during SFP outbreak investigation. Although confirming the production of enterotoxins remains important for SFP investigation, we delivered a proof-of-concept that WGS is a valid alternative and/or complementary tool for outbreak investigation.

Evaluating coverage bias in next-generation sequencing of Escherichia coli

Whole-genome sequencing is essential to many facets of infectious disease research. However, technical limitations such as bias in coverage and tagmentation, and difficulties characterising genomic regions with extreme GC content have created significant obstacles in its use. Illumina has claimed that the recently released DNA Prep library preparation kit, formerly known as Nextera Flex, overcomes some of these limitations. This study aimed to assess bias in coverage, tagmentation, GC content, average fragment size distribution, and de novo assembly quality using both the Nextera XT and DNA Prep kits from Illumina. When performing whole-genome sequencing on Escherichia coli and where coverage bias is the main concern, the DNA Prep kit may provide higher quality results; though de novo assembly quality, tagmentation bias and GC content related bias are unlikely to improve. Based on these results, laboratories with existing workflows based on Nextera XT would see minor benefits in transitioning to the DNA Prep kit if they were primarily studying organisms with neutral GC content.

Responses of bacterial communities and organic matter degradation in surface sediment to Macrobrachium nipponense bioturbation

The excessive accumulation of organic matter (OM) in sediments in aquaculture ponds is a potential environmental threat due to the risk of endogenous water pollution and eutrophication. From the perspective of inhibiting OM accumulation to prevent endogenous water pollution, the present study investigated the OM degradation states, variations of bacterial communities and basic environmental factors in sediments with/without Macrobrachium nipponense treatment/control groups in triplicate for effects of bioturbation on OM degradation in 90-day incubation. The total organic carbon (TOC) and total nitrogen (TN) in the M. nipponense treatment were higher than in the control at the 30th and 60th days, while no significant differences between treatment and control were found at the end of the experiment. Significantly higher oxidation-reduction potential (ORP) and more extensively degraded OM were observed in the M. nipponense treatment. Eleven significantly differential bacterial taxa were enriched in the sediments of M. nipponense treatment, of which eight (Actinobacteria, Chitinophagales, Chitinophagaceae, Flavihumibacter, Marinilabiliaceae, Cytophaga xylanolytica group, Christensenellaceae, and Christensenellaceae R-7 group) were significantly correlated with at least two OM degradation indicators. The functional groups chemoheterotrophy, aerobic chemoheterotrophy, xylanolysis, ureolysis, and intracellular parasites were enhanced by M. nipponense and were negatively correlated with OM degradation indictors. Overall, the M. nipponense bioturbation effectively increased the ORP to provide better conditions for OM degradation, altered the taxonomic composition and functional groups to enhance the bacterial ability for OM degradation, and finally promoted the OM degradation of the surface sediment in an artificial aquaculture system.

Impact of DNA extraction on whole genome sequencing analysis for characterization and relatedness of Shiga toxin-producing Escherichia coli isolates

Whole genome sequencing (WGS) has proven to be the ultimate tool for bacterial isolate characterization and relatedness determination. However, standardized and harmonized workflows, e.g. for DNA extraction, are required to ensure robust and exchangeable WGS data. Data sharing between (inter)national laboratories is essential to support foodborne pathogen control, including outbreak investigation. This study evaluated eight commercial DNA preparation kits for their potential influence on: (i) DNA quality for Nextera XT library preparation; (ii) MiSeq sequencing (data quality, read mapping against plasmid and chromosome references); and (iii) WGS data analysis, i.e. isolate characterization (serotyping, virulence and antimicrobial resistance genotyping) and phylogenetic relatedness (core genome multilocus sequence typing and single nucleotide polymorphism analysis). Shiga toxin-producing Escherichia coli (STEC) was selected as a case study. Overall, data quality and inferred phylogenetic relationships between isolates were not affected by the DNA extraction kit choice, irrespective of the presence of confounding factors such as EDTA in DNA solution buffers. Nevertheless, completeness of STEC characterization was, although not substantially, influenced by the plasmid extraction performance of the kits, especially when using Nextera XT library preparation. This study contributes to addressing the WGS challenges of standardizing protocols to support data portability and to enable full exploitation of its potential.

The Benefits of Whole Genome Sequencing for Foodborne Outbreak Investigation from the Perspective of a National Reference Laboratory in a Smaller Country

Gradually, conventional methods for foodborne pathogen typing are replaced by whole genome sequencing (WGS). Despite studies describing the overall benefits, National Reference Laboratories of smaller countries often show slower uptake of WGS, mainly because of significant investments required to generate and analyze data of a limited amount of samples. To facilitate this process and incite policy makers to support its implementation, a Shiga toxin-producing Escherichia coli (STEC) O157:H7 (stx1+, stx2+, eae+) outbreak (2012) and a STEC O157:H7 (stx2+, eae+) outbreak (2013) were retrospectively analyzed using WGS and compared with their conventional investigations. The corresponding results were obtained, with WGS delivering even more information, e.g., on virulence and antimicrobial resistance genotypes. Besides a universal, all-in-one workflow with less hands-on-time (five versus seven actual working days for WGS versus conventional), WGS-based cgMLST-typing demonstrated increased resolution. This enabled an accurate cluster definition, which remained unsolved for the 2013 outbreak, partly due to scarce epidemiological linking with the suspect source. Moreover, it allowed detecting two and one earlier circulating STEC O157:H7 (stx1+, stx2+, eae+) and STEC O157:H7 (stx2+, eae+) strains as closely related to the 2012 and 2013 outbreaks, respectively, which might have further directed epidemiological investigation initially. Although some bottlenecks concerning centralized data-sharing, sampling strategies, and perceived costs should be considered, we delivered a proof-of-concept that even in smaller countries, WGS offers benefits for outbreak investigation, if a sufficient budget is available to ensure its implementation in surveillance. Indeed, applying a database with background isolates is critical in interpreting isolate relationships to outbreaks, and leveraging the true benefit of WGS in outbreak investigation and/or prevention.