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.

Demographics of sources of HIV-1 transmission in Zambia: a molecular epidemiology analysis in the HPTN 071 PopART study.

BACKGROUND: In the last decade, universally available antiretroviral therapy (ART) has led to greatly improved health and survival of people living with HIV in sub-Saharan Africa, but new infections continue to appear. The design of effective prevention strategies requires the demographic characterisation of individuals acting as sources of infection, which is the aim of this study. METHODS: Between 2014 and 2018, the HPTN 071 PopART study was conducted to quantify the public health benefits of ART. Viral samples from 7124 study participants in Zambia were deep-sequenced as part of HPTN 071-02 PopART Phylogenetics, an ancillary study. We used these sequences to identify likely transmission pairs. After demographic weighting of the recipients in these pairs to match the overall HIV-positive population, we analysed the demographic characteristics of the sources to better understand transmission in the general population. FINDINGS: We identified a total of 300 likely transmission pairs. 178 (59·4%) were male to female, with 130 (95% CI 110-150; 43·3%) from males aged 25-40 years. Overall, men transmitted 2·09-fold (2·06-2·29) more infections per capita than women, a ratio peaking at 5·87 (2·78-15·8) in the 35-39 years source age group. 40 (26-57; 13·2%) transmissions linked individuals from different communities in the trial. Of 288 sources with recorded information on drug resistance mutations, 52 (38-69; 18·1%) carried viruses resistant to first-line ART. INTERPRETATION: HIV-1 transmission in the HPTN 071 study communities comes from a wide range of age and sex groups, and there is no outsized contribution to new infections from importation or drug resistance mutations. Men aged 25-39 years, underserved by current treatment and prevention services, should be prioritised for HIV testing and ART. FUNDING: National Institute of Allergy and Infectious Diseases, US President's Emergency Plan for AIDS Relief, International Initiative for Impact Evaluation, Bill & Melinda Gates Foundation, National Institute on Drug Abuse, and National Institute of Mental Health.

Targeted metagenomics reveals association between severity and pathogen co-detection in infants with respiratory syncytial virus.

Respiratory syncytial virus (RSV) is the leading cause of hospitalisation for respiratory infection in young children. RSV disease severity is known to be age-dependent and highest in young infants, but other correlates of severity, particularly the presence of additional respiratory pathogens, are less well understood. In this study, nasopharyngeal swabs were collected from two cohorts of RSV-positive infants <12 months in Spain, the UK, and the Netherlands during 2017-20. We show, using targeted metagenomic sequencing of >100 pathogens, including all common respiratory viruses and bacteria, from samples collected from 433 infants, that burden of additional viruses is common (111/433, 26%) but only modestly correlates with RSV disease severity. In contrast, there is strong evidence in both cohorts and across age groups that presence of Haemophilus bacteria (194/433, 45%) is associated with higher severity, including much higher rates of hospitalisation (odds ratio 4.25, 95% CI 2.03-9.31). There is no evidence for association between higher severity and other detected bacteria, and no difference in severity between RSV genotypes. Our findings reveal the genomic diversity of additional pathogens during RSV infection in infants, and provide an evidence base for future causal investigations of the impact of co-infection on RSV disease severity.

Phylogenetic Analysis of Hepatitis C Virus Infections in a Large Belgian Cohort Using Next-Generation Sequencing of Full-Length Genomes.

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium.

A highly virulent variant of HIV-1 circulating in the Netherlands.

We discovered a highly virulent variant of subtype-B HIV-1 in the Netherlands. One hundred nine individuals with this variant had a 0.54 to 0.74 log10 increase (i.e., a ~3.5-fold to 5.5-fold increase) in viral load compared with, and exhibited CD4 cell decline twice as fast as, 6604 individuals with other subtype-B strains. Without treatment, advanced HIV-CD4 cell counts below 350 cells per cubic millimeter, with long-term clinical consequences-is expected to be reached, on average, 9 months after diagnosis for individuals in their thirties with this variant. Age, sex, suspected mode of transmission, and place of birth for the aforementioned 109 individuals were typical for HIV-positive people in the Netherlands, which suggests that the increased virulence is attributable to the viral strain. Genetic sequence analysis suggests that this variant arose in the 1990s from de novo mutation, not recombination, with increased transmissibility and an unfamiliar molecular mechanism of virulence.

Targeted capture and sequencing of Orientia tsutsugamushi genomes from chiggers and humans.

Scrub typhus is a febrile disease caused by Orientia tsutsugamushi, transmitted by larval stage Trombiculid mites (chiggers), whose primary hosts are small mammals. The phylogenomics of O. tsutsugamushi in chiggers, small mammals and humans remains poorly understood. To combat the limitations imposed by the low relative quantities of pathogen DNA in typical O. tsutsugamushi clinical and ecological samples, along with the technical, safety and cost limitations of cell culture, a novel probe-based target enrichment sequencing protocol was developed. The method was designed to capture variation among conserved genes and facilitate phylogenomic analysis at the scale of population samples. A whole-genome amplification step was incorporated to enhance the efficiency of sequencing by reducing duplication rates. This resulted in on-target capture rates of up to 93% for a diverse set of human, chigger, and rodent samples, with the greatest success rate in samples with real-time PCR Ct values below 35. Analysis of the best-performing samples revealed phylogeographic clustering at local, provincial and international scales. Applying the methodology to a comprehensive set of samples could yield a more complete understanding of the ecology, genomic evolution and population structure of O. tsutsugamushi and other similarly challenging organisms, with potential benefits in the development of diagnostic tests and vaccines.

Distinct patterns of within-host virus populations between two subgroups of human respiratory syncytial virus.

Human respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infection in young children globally, but little is known about within-host RSV diversity. Here, we characterised within-host RSV populations using deep-sequencing data from 319 nasopharyngeal swabs collected during 2017-2020. RSV-B had lower consensus diversity than RSV-A at the population level, while exhibiting greater within-host diversity. Two RSV-B consensus sequences had an amino acid alteration (K68N) in the fusion (F) protein, which has been associated with reduced susceptibility to nirsevimab (MEDI8897), a novel RSV monoclonal antibody under development. In addition, several minor variants were identified in the antigenic sites of the F protein, one of which may confer resistance to palivizumab, the only licensed RSV monoclonal antibody. The differences in within-host virus populations emphasise the importance of monitoring for vaccine efficacy and may help to explain the different prevalences of monoclonal antibody-escape mutants between the two subgroups.

SARS-CoV-2 within-host diversity and transmission.

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

Oxford Nanopore MinION Sequencing Enables Rapid Whole Genome Assembly of Rickettsia typhi in a Resource-Limited Setting.

The infrastructure challenges and costs of next-generation sequencing have been largely overcome, for many sequencing applications, by Oxford Nanopore Technologies' portable MinION sequencer. However, the question remains open whether MinION-based bacterial whole genome sequencing is by itself sufficient for the accurate assessment of phylogenetic and epidemiological relationships between isolates and whether such tasks can be undertaken in resource-limited settings. To investigate this question, we sequenced the genome of an isolate of Rickettsia typhi, an important and neglected cause of fever across much of the tropics and subtropics, for which only three genomic sequences previously existed. We prepared and sequenced libraries on a MinION in Vientiane, Lao PDR, using v9.5 chemistry, and in parallel, we sequenced the same isolate on the Illumina platform in a genomics laboratory in the United Kingdom. The MinION sequence reads yielded a single contiguous assembly, in which the addition of Illumina data revealed 226 base-substitution and 5,856 indel errors. The combined assembly represents the first complete genome sequence of a human R. typhi isolate collected in the last 50 years and differed from the genomes of existing strains collected over a 90-year time period at very few sites, with no rearrangements. Filtering based on the known error profile of MinION data improved the accuracy of the nanopore-only assembly. However, the frequency of false-positive errors remained greater than true sequence divergence from recorded sequences. Although nanopore-only sequencing cannot yet recover phylogenetic signals in R. typhi, such an approach may be applicable for more diverse organisms.

Case Report: Application of hepatitis B virus (HBV) deep sequencing to distinguish between acute and chronic infection.

Deep sequencing of the full-length hepatitis B virus (HBV) genome provides the opportunity to determine the extent to which viral diversity, genotype, polymorphisms, insertions and deletions may influence presentation and outcomes of disease. Increasing experience with analysis of HBV genomic data opens up the potential for using these data to inform insights into pathophysiology of infection and to underpin decision making in clinical practice. We here set out to undertake whole genome HBV sequencing from an adult who presented acutely unwell with a new diagnosis of HBV infection, and tested positive for both HBV anti-core IgM and IgG, possibly representing either acute hepatitis B infection (AHB) or chronic hepatitis B with an acute reactivation (CHB-AR). The distinction between these two scenarios may be important in predicting prognosis and underpinning treatment decisions, but can be challenging based on routine laboratory tests. Through application of deep whole-genome sequencing we typed the isolate as genotype-D1, and identified several minority variants including G1764A and G1986A substitutions in the pre-core promoter and pre-core regions, which support CHB-AR rather than AHB. In the longer term, enhanced deep sequencing data for HBV may provide improved evidence to distinguish between acute and chronic infection, to predict outcomes and to stratify treatment.

Epstein-Barr virus reactivation in sepsis due to community-acquired pneumonia is associated with increased morbidity and an immunosuppressed host transcriptomic endotype.

Epstein-Barr virus (EBV) reactivation is common in sepsis patients but the extent and nature of this remains unresolved. We sought to determine the incidence and correlates of EBV-positivity in a large sepsis cohort. We also hypothesised that EBV reactivation would be increased in patients in whom relative immunosuppression was the major feature of their sepsis response. To identify such patients we aimed to use knowledge of sepsis response subphenotypes based on transcriptomic studies of circulating leukocytes, specifically patients with a Sepsis Response Signature endotype (SRS1) that we have previously shown to be associated with increased mortality and features of immunosuppression. We assayed EBV from the plasma of intensive care unit (ICU) patients with sepsis due to community-acquired pneumonia. In total 730 patients were evaluated by targeted metagenomics (n = 573 patients), digital droplet PCR (n = 565), or both (n = 408). We had previously analysed gene expression in peripheral blood leukocytes for a subset of individuals (n = 390). We observed a 37% incidence of EBV-positivity. EBV reactivation was associated with longer ICU stay (12.9 vs 9.2 days; p = 0.004) and increased organ failure (day 1 SOFA score 6.9 vs 5.9; p = 0.00011). EBV reactivation was associated with the relatively immunosuppressed SRS1 endotype (p = 0.014) and differential expression of a small number of biologically relevant genes. These findings are consistent with the hypothesis that viral reactivation in sepsis is a consequence of immune compromise and is associated with increasing severity of illness although further mechanistic studies are required to definitively illustrate cause and effect.

Evidence of tenofovir resistance in chronic hepatitis B virus (HBV) infection: An observational case series of South African adults.

INTRODUCTION: Tenofovir disoproxil fumarate (TDF) is widely recommended for treatment of chronic hepatitis B virus (HBV) infection because it is safe, affordable and has a high genetic barrier to resistance. TDF resistance associated mutations (RAMs) have been reported, but data are limited, particularly for Africa. We set out to identify potential RAMs in individuals with detectable HBV viraemia on TDF treatment. METHODS: We recruited adults with chronic HBV infection from Cape Town, South Africa, identifying individuals with a TDF resistance phenotype, defined as persistent HBV vireamia despite >12 months of TDF treatment. We sequenced HBV DNA using MiSeq Illumina with whole genome target enrichment, and sought potential TDF RAMs, based on a pre-defined list of polymorphisms. RESULTS: Among 66 individuals with chronic HBV (genotypes A and D), three met our clinical definition for TDF resistance, of whom two were coinfected with HIV. In one participant, the consensus HBV sequence contained nine polymorphisms that have been described in association with TDF resistance. Significant treatment non-adherence in this individual was unlikely, as HIV RNA was suppressed. TDF RAMs were also present in HBV sequences from the other two participants, but other factors including treatment non-adherence may also have had a role in failure of HBV DNA suppression in these cases. DISCUSSION: Our findings add to the evidence that RAMs in HBV reverse transcriptase may underpin a TDF resistant phenotype. This is the first time these RAMs have been reported from Africa in association with clinical evidence of TDF resistance.