Genomics of phage-bacterium-host interaction of mushroom pathogenic Pseudomonads

The most commonly cultivated mushroom in Europe and North America is the Agaricus bisporus, also known as the button mushroom or Portobello mushroom. Bacterial diseases of Agaricus bisporus caused by Pseudomonas species are a cause of significant crop loss and downgrading of produce, resulting in considerable economic cost. Bacteriophage have long been an attractive option for biocontrol of bacterial contamination of food products, however the precise genetic interactions between phage, bacterium and host are often inadequately explored. This project aims to explore the genetic interactions between the mushroom pathogenic bacterium Pseudomonas tolaasii and Pseudomonas agarici with the newly identified bacteriophage Pseudomonas phage NV1 and Pseudomonas phage ϕNV3. Full genome sequencing has been performed on the P. tolaasii strain 2192T and P. agarici NCPPB 2472, and the genomes of both mined for potential biosynthetic clusters involved in virulence as well as genes involved in phage resistance. Within the genome of P. tolaasii 2192T, putative non-ribosomal peptide synthases have been identified which are hypothesised to be involved in the production of the tolaasin toxin involved in disease symptom appearance on mushroom surfaces. Within the genome of P. agarici NCPPB 2472, a biosynthetic cluster was identified that is hypothesised to produce the siderophore achromobactin, an important virulence factor. P. agarici NCPPB 2472 was identified as possessing a single Type I-F CRISPR/Cas system, predicted to be involved in the development of phage resistance, as well as complete operon predicted to be involved in the production of the exopolysaccharide alginate. A third Pseudomonas species was identified on the surface of disease free mushrooms which was identified as a potentially new species of Pseudomonas, named Pseudomonas sp. NS1. The genome of P. sp. NS1 was likewise sequenced and mined for potential biosynthetic gene clusters which identified a cluster demonstrated to be involved in the production of White Line Inducing Principle. The Pseudomonas phages NV1 and ϕNV3 were isolated from environmental samples and identified to be narrow host range phage specific for P. tolaasii 2192T and P. agarici NCPPB 2472 respectively. Both phage NV1 and ϕNV3 were identified as new species of the Luz24likeviruses and phiKMVlikeviruses respectively. The genomes of both phages were isolated and sequenced, with phage ϕNV3 identified as containing a conserved Signal-Arrest-Release endolysin system, which was confirmed by in vitro protein expression. Likewise, the lysis cassette proteins of ϕNV3 were identified and investigated via protein complementation assay in vitro. The full growth characteristics and life cycle of phage ϕNV3 has been investigated and reported in this study and a broad-host range mutant of ϕNV3 identified which has allowed the T7-like tail protein of ϕNV3 to be identified as the host specificity determinant. Transcriptome analyses of non-infected P. agarici NCPPB 2472 and P. agarici NCPPB 2472 infected with phage ϕNV3 at a multiplicity of infection (MOI) of 1 at 40 min post infection, were performed in triplicate using RNA-seq. A reliable method has been established that will be useful in future studies, although comparative gene expression analysis revealed no significant differences in expression between the two treatments at the multiplicity of infection and time point chosen in this case a significant quantity of phage transcripts were detected, demonstrating active phage infection

A comparison of phenotypic and WGS drug susceptibility testing in Mycobacterium tuberculosis isolates from the Republic of Korea

BACKGROUND: WGS has significant potential to help tackle the major public health problem of TB. The Republic of Korea has the third highest rates of TB of all Organisation for Economic Cooperation and Development countries but there has been very limited use of WGS in TB to date. OBJECTIVES: A retrospective comparison of Mycobacterium tuberculosis (MTB) clinical isolates from 2015 to 2017 from two centres in the Republic of Korea using WGS to compare phenotypic drug susceptibility testing (pDST) and WGS drug susceptibility predictions (WGS-DSP). METHODS: Fifty-seven MTB isolates had DNA extracted and were sequenced using the Illumina HiSeq platform. The WGS analysis was performed using bwa mem, bcftools and IQ-Tree; resistance markers were identified using TB profiler. Phenotypic susceptibilities were carried out at the Supranational TB reference laboratory (Korean Institute of Tuberculosis). RESULTS: For first-line antituberculous drugs concordance for rifampicin, isoniazid, pyrazinamide and ethambutol was 98.25%, 92.98%, 87.72% and 85.96%, respectively. The sensitivity of WGS-DSP compared with pDST for rifampicin, isoniazid, pyrazinamide and ethambutol was 97.30%, 92.11%, 78.95% and 95.65%, respectively. The specificity for these first-line antituberculous drugs was 100%, 94.74%, 92.11% and 79.41%, respectively. The sensitivity and specificity for second-line drugs ranged from 66.67% to 100%, and from 82.98% to 100%, respectively. CONCLUSIONS: This study confirms the potential role for WGS in drug susceptibility prediction, which would reduce turnaround times. However, further larger studies are needed to ensure current databases of drug resistance mutations are reflective of the TB present in the Republic of Korea

Single base mutations in the nucleocapsid gene of SARS-CoV-2 affects amplification efficiency of sequence variants and may lead to assay failure

Reverse transcriptase quantitative PCR (RT-qPCR) is the main diagnostic assay used to detect SARS-CoV-2 RNA in respiratory samples. RT-qPCR is performed by specifically targeting the viral genome using complementary oligonucleotides called primers and probes. This approach relies on prior knowledge of the genetic sequence of the target. Viral genetic variants with changes to the primer/probe binding region may reduce the performance of PCR assays and have the potential to cause assay failure. In this work we demonstrate how two single nucleotide variants (SNVs) altered the amplification curve of a diagnostic PCR targeting the Nucleocapsid (N) gene and illustrate how threshold setting can lead to false-negative results even where the variant sequence is amplified. We also describe how in silico analysis of SARS-CoV-2 genome sequences available in the COVID-19 Genomics UK Consortium (COG-UK) and GISAID databases was performed to predict the impact of sequence variation on the performance of 22 published PCR assays. The vast majority of published primer and probe sequences contain sequence mismatches with at least one SARS-CoV-2 lineage. We recommend that visual observation of amplification curves is included as part of laboratory quality procedures, even in high throughput settings where thresholds are set automatically and that in silico analysis is used to monitor the potential impact of new variants on established assays. Ideally comprehensive in silico analysis should be applied to guide selection of highly conserved genomic regions to target with future SARS-CoV-2 PCR assays

Untargeted metagenomics protocol for the diagnosis of infection from CSF and tissue from sterile sites

Metagenomic next-generation sequencing (mNGS) is an untargeted technique capable of detecting all microbial nucleic acid within a sample. This protocol outlines our wet laboratory method for mNGS of cerebrospinal fluid (CSF) specimens and tissues from sterile sites. We use this method routinely in our clinical service, processing 178 specimens over the past 2.5 years in a laboratory that adheres to ISO:15189 standards. We have successfully used this protocol to diagnose multiple cases of encephalitis and hepatitis

Genomic investigations of unexplained acute hepatitis in children

Since its first identification in Scotland, over 1,000 cases of unexplained paediatric hepatitis in children have been reported worldwide, including 278 cases in the UK1. Here we report an investigation of 38 cases, 66 age-matched immunocompetent controls and 21 immunocompromised comparator participants, using a combination of genomic, transcriptomic, proteomic and immunohistochemical methods. We detected high levels of adeno-associated virus 2 (AAV2) DNA in the liver, blood, plasma or stool from 27 of 28 cases. We found low levels of adenovirus (HAdV) and human herpesvirus 6B (HHV-6B) in 23 of 31 and 16 of 23, respectively, of the cases tested. By contrast, AAV2 was infrequently detected and at low titre in the blood or the liver from control children with HAdV, even when profoundly immunosuppressed. AAV2, HAdV and HHV-6 phylogeny excluded the emergence of novel strains in cases. Histological analyses of explanted livers showed enrichment for T cells and B lineage cells. Proteomic comparison of liver tissue from cases and healthy controls identified increased expression of HLA class 2, immunoglobulin variable regions and complement proteins. HAdV and AAV2 proteins were not detected in the livers. Instead, we identified AAV2 DNA complexes reflecting both HAdV-mediated and HHV-6B-mediated replication. We hypothesize that high levels of abnormal AAV2 replication products aided by HAdV and, in severe cases, HHV-6B may have triggered immune-mediated hepatic disease in genetically and immunologically predisposed children

Evaluating metagenomics and targeted approaches for diagnosis and surveillance of viruses

: Background : Metagenomics is a powerful approach for the detection of unknown and novel pathogens. Workflows based on Illumina short-read sequencing are becoming established in diagnostic laboratories. However, high sequencing depth requirements, long turnaround times, and limited sensitivity hinder broader adoption. We investigated whether we could overcome these limitations using protocols based on untargeted sequencing with Oxford Nanopore Technologies (ONT), which offers real-time data acquisition and analysis, or a targeted panel approach, which allows the selective sequencing of known pathogens and could improve sensitivity. Methods: We evaluated detection of viruses with readily available untargeted metagenomic workflows using Illumina and ONT, and an Illumina-based enrichment approach using the Twist Bioscience Comprehensive Viral Research Panel (CVRP), which targets 3153 viruses. We tested samples consisting of a dilution series of a six-virus mock community in a human DNA/RNA background, designed to resemble clinical specimens with low microbial abundance and high host content. Protocols were designed to retain the host transcriptome, since this could help confirm the absence of infectious agents. We further compared the performance of commonly used taxonomic classifiers. Results: Capture with the Twist CVRP increased sensitivity by at least 10–100-fold over untargeted sequencing, making it suitable for the detection of low viral loads (60 genome copies per ml (gc/ml)), but additional methods may be needed in a diagnostic setting to detect untargeted organisms. While untargeted ONT had good sensitivity at high viral loads (60,000 gc/ml), at lower viral loads (600–6000 gc/ml), longer and more costly sequencing runs would be required to achieve sensitivities comparable to the untargeted Illumina protocol. Untargeted ONT provided better specificity than untargeted Illumina sequencing. However, the application of robust thresholds standardized results between taxonomic classifiers. Host gene expression analysis is optimal with untargeted Illumina sequencing but possible with both the CVRP and ONT. Conclusions: Metagenomics has the potential to become standard-of-care in diagnostics and is a powerful tool for the discovery of emerging pathogens. Untargeted Illumina and ONT metagenomics and capture with the Twist CVRP have different advantages with respect to sensitivity, specificity, turnaround time and cost, and the optimal method will depend on the clinical context

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control

Multicenter benchmarking of short and long read wet lab protocols for clinical viral metagenomics

Metagenomics is gradually being implemented for diagnosing infectious diseases. However, in-depth protocol comparisons for viral detection have been limited to individual sets of experimental workflows and laboratories. In this study, we present a benchmark of metagenomics protocols used in clinical diagnostic laboratories initiated by the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS). A mock viral reference panel was designed to mimic low biomass clinical specimens. The panel was used to assess the performance of twelve metagenomic wet lab protocols currently in use in the diagnostic laboratories of participating ENNGS member institutions. Both Illumina and Nanopore, shotgun and targeted capture probe protocols were included. Performance metrics sensitivity, specificity, and quantitative potential were assessed using a central bioinformatics pipeline. Overall, viral pathogens with loads down to 104 copies/ml (corresponding to CT values of 31 in our PCR assays) were detected by all the evaluated metagenomic wet lab protocols. In contrast, lower abundant mixed viruses of CT values of 35 and higher were detected only by a minority of the protocols. Considering the reference panel as the gold standard, optimal thresholds to define a positive result were determined per protocol, based on the horizontal genome coverage. Implementing these thresholds, sensitivity and specificity of the protocols ranged from 67 to 100 % and 87 to 100 %, respectively. A variety of metagenomic protocols are currently in use in clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implying the need for standardization of metagenomic analysis for use in clinical settings.</p