Castanet: A pipeline for rapid analysis of targeted multi-pathogen genomic data.

MOTIVATION: Target enrichment strategies generate genomic data from multiple pathogens in a single process, greatly improving sensitivity over metagenomic sequencing and enabling cost-effective, high throughput surveillance and clinical applications. However, uptake by research and clinical laboratories is constrained by an absence of computational tools that are specifically designed for the analysis of multi-pathogen enrichment sequence data. Here we present an analysis pipeline, Castanet, for use with multi-pathogen enrichment sequencing data. Castanet is designed to work with short-read data produced by existing targeted enrichment strategies, but can be readily deployed on any BAM file generated by another methodology. Also included are an optional graphical interface and installer script. RESULTS: In addition to genome reconstruction, Castanet reports method-specific metrics that enable quantification of capture efficiency, estimation of pathogen load, differentiation of low-level positives from contamination, and assessment of sequencing quality. Castanet can be used as a traditional end-to-end pipeline for consensus generation, but its strength lies in the ability to process a flexible, pre-defined set of pathogens of interest directly from multi-pathogen enrichment experiments. In our tests, Castanet consensus sequences were accurate reconstructions of reference sequences, including in instances where multiple strains of the same pathogen were present. Castanet performs effectively on standard computers and can process the entire output of a 96-sample enrichment sequencing run (50M reads) using a single batch process command, in $<$2 h. AVAILABILITY AND IMPLEMENTATION: Source code freely available under GPL-3 license at https://github.com/MultipathogenGenomics/castanet, implemented in Python 3.10 and supported in Ubuntu Linux 22.04. SUPPLEMENTARY INFORMATION: Supplementary data available at the journal's web site. Supporting sequencing data available at https://www.ebi.ac.uk/ena/browser/view/PRJEB77004.

The mutational spectrum of SARS-CoV-2 genomic and antigenomic RNA.

The raw material for viral evolution is provided by intra-host mutations occurring during replication, transcription or post-transcription. Replication and transcription of Coronaviridae proceed through the synthesis of negative-sense 'antigenomes' acting as templates for positive-sense genomic and subgenomic RNA. Hence, mutations in the genomes of SARS-CoV-2 and other coronaviruses can occur during (and after) the synthesis of either negative-sense or positive-sense RNA, with potentially distinct patterns and consequences. We explored for the first time the mutational spectrum of SARS-CoV-2 (sub)genomic and anti(sub)genomic RNA. We use a high-quality deep sequencing dataset produced using a quantitative strand-aware sequencing method, controlled for artefacts and sequencing errors, and scrutinized for accurate detection of within-host diversity. The nucleotide differences between negative- and positive-sense strand consensus vary between patients and do not show dependence on age or sex. Similarities and differences in mutational patterns between within-host minor variants on the two RNA strands suggested strand-specific mutations or editing by host deaminases and oxidative damage. We observe generally neutral and slight negative selection on the negative strand, contrasting with purifying selection in ORF1a, ORF1b and S genes of the positive strand of the genome.

Endemic HBV among hospital in-patients in Bangladesh, including evidence of occult infection.

Bangladesh is one of the top-ten most heavily burdened countries for viral hepatitis, with hepatitis B (HBV) infections responsible for the majority of cases. Recombinant and occult HBV infections (OBI) have been reported previously in the region. We investigated an adult fever cohort (n=201) recruited in Dhaka, to determine the prevalence of HBV and OBI. A target-enrichment deep sequencing pipeline was applied to samples with HBV DNA >3.0 log10 IU ml-1. HBV infection was present in 16/201 (8 %), among whom 3/16 (19 %) were defined as OBI (HBsAg-negative but detectable HBV DNA). Whole genome deep sequences (WGS) were obtained for four cases, identifying genotypes A, C and D. One OBI case had sufficient DNA for sequencing, revealing multiple polymorphisms in the surface gene that may contribute to the occult phenotype. We identified mutations associated with nucleos(t)ide analogue resistance in 3/4 samples sequenced, although the clinical significance in this cohort is unknown. The high prevalence of HBV in this setting illustrates the importance of opportunistic clinical screening and DNA testing of transfusion products to minimise OBI transmission. WGS can inform understanding of diverse disease phenotypes, supporting progress towards international targets for HBV elimination.

Simultaneous Viral Whole-Genome Sequencing and Differential Expression Profiling in Respiratory Syncytial Virus Infection of Infants.

Targeted metagenomics using strand-specific libraries with target enrichment is a sensitive, generalized approach to pathogen sequencing and transcriptome profiling. Using this method, we recovered 13 (76%) complete human respiratory syncytial virus (RSV) genomes from 17 clinical respiratory samples, reconstructed the phylogeny of the infecting viruses, and detected differential gene expression between 2 RSV subgroups, specifically, a lower expression of the P gene and a higher expression of the M2 gene in RSV-A than in RSV-B. This methodology can help to relate viral genetics to clinical phenotype and facilitate ongoing population-level RSV surveillance and vaccine development. Clinical Trials Registration. NCT03627572 and NCT03756766.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings.

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Performance of a high-throughput next-generation sequencing method for analysis of HIV drug resistance and viral load.

OBJECTIVES: To evaluate the performance of a high-throughput research assay for HIV drug resistance testing based on whole genome next-generation sequencing (NGS) that also quantifies HIV viral load. METHODS: Plasma samples (n = 145) were obtained from HIV-positive MSM (HPTN 078). Samples were analysed using clinical assays (the ViroSeq HIV-1 Genotyping System and the Abbott RealTime HIV-1 Viral Load assay) and a research assay based on whole-genome NGS (veSEQ-HIV). RESULTS: HIV protease and reverse transcriptase sequences (n = 142) and integrase sequences (n = 138) were obtained using ViroSeq. Sequences from all three regions were obtained for 100 (70.4%) of the 142 samples using veSEQ-HIV; results were obtained more frequently for samples with higher viral loads (93.5% for 93 samples with >5000 copies/mL; 50.0% for 26 samples with 1000-5000 copies/mL; 0% for 23 samples with <1000 copies/mL). For samples with results from both methods, drug resistance mutations (DRMs) were detected in 33 samples using ViroSeq and 42 samples using veSEQ-HIV (detection threshold: 5.0%). Overall, 146 major DRMs were detected; 107 were detected by both methods, 37 were detected by veSEQ-HIV only (frequency range: 5.0%-30.6%) and two were detected by ViroSeq only. HIV viral loads estimated by veSEQ-HIV strongly correlated with results from the Abbott RealTime Viral Load assay (R2 = 0.85; n = 142). CONCLUSIONS: The NGS-based veSEQ-HIV method provided results for most samples with higher viral loads, was accurate for detecting major DRMs, and detected mutations at lower levels compared with a method based on population sequencing. The veSEQ-HIV method also provided HIV viral load data.

Transfusion-transmitted hepatitis C: A cluster of cases in transfusion-dependent thalassaemia patients in Sri Lanka.

OBJECTIVES: To report the clinical and virologic epidemiology of a recent epidemic of hepatitis C in thalassaemia patients in Sri Lanka. BACKGROUND: Transfusion-dependent thalassaemia patients remain at risk for hepatitis C virus (HCV). Here, we report a cluster of recent HCV infections in Sri Lankan thalassaemia patients and examine the phylogenetic relationship of viral sequences. METHODS: We conducted two prospective cross-sectional surveys of 513 patients in four Sri Lankan thalassaemia centres in 2014/2015 and re-surveyed one centre in 2016. We screened for anti-HCV antibodies using the CTK Biotech enzyme-linked immunosorbent assay (ELISA) kits and confirmed active infection by reverse transcription-polymerase chain reaction (RT-PCR) for HCV-RNA. HCV genomes were sequenced by unbiased target enrichment. RESULTS: Anti-HCV antibodies were found in 116/513 (22.6%) of patients initially tested. Active hepatitis C infection was found in 26 patients with no cases of active hepatitis B infection. Of 26 patients with HCV, two were infected with genotype 1(a), and the rest had 3(a). In a single centre (Ragama), 122 patients (120 new cases and two previously tested, but negative) were retested for anti-HCV antibodies. 32/122 (26.2%) patients were seropositive. Twenty-three (23/122; 18.8%) of these new cases were confirmed by HCV PCR (all genotype 3[a]). CONCLUSION: There is a significant cluster of recent HCV cases in multiply transfused thalassaemia patients in several centres in Sri Lanka. Most of the viruses shared a close phylogenetic relationship. The results are consistent with recent continuing transfusion-transmitted HCV infection. Routine surveillance for HCV of chronically transfused patients is required irrespective of screening of blood products.

Prolonged Zika Virus RNA Detection in Semen of Immunosuppressed Patient.

Zika virus RNA has been detected in semen samples collected <370 days after symptom onset. We report unusual persistence of Zika virus RNA in semen, confirmed by sequencing at 515 days after symptom onset and detectable for >900 days, in a patient with immunosuppression.

PHYLOSCANNER: Inferring Transmission from Within- and Between-Host Pathogen Genetic Diversity.

A central feature of pathogen genomics is that different infectious particles (virions and bacterial cells) within an infected individual may be genetically distinct, with patterns of relatedness among infectious particles being the result of both within-host evolution and transmission from one host to the next. Here, we present a new software tool, phyloscanner, which analyses pathogen diversity from multiple infected hosts. phyloscanner provides unprecedented resolution into the transmission process, allowing inference of the direction of transmission from sequence data alone. Multiply infected individuals are also identified, as they harbor subpopulations of infectious particles that are not connected by within-host evolution, except where recombinant types emerge. Low-level contamination is flagged and removed. We illustrate phyloscanner on both viral and bacterial pathogens, namely HIV-1 sequenced on Illumina and Roche 454 platforms, HCV sequenced with the Oxford Nanopore MinION platform, and Streptococcus pneumoniae with sequences from multiple colonies per individual. phyloscanner is available from https://github.com/BDI-pathogens/phyloscanner.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing.

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.