Mutational signature dynamics indicate SARS-CoV-2's evolutionary capacity is driven by host antiviral molecules.

The COVID-19 pandemic has been characterised by sequential variant-specific waves shaped by viral, individual human and population factors. SARS-CoV-2 variants are defined by their unique combinations of mutations and there has been a clear adaptation to more efficient human infection since the emergence of this new human coronavirus in late 2019. Here, we use machine learning models to identify shared signatures, i.e., common underlying mutational processes and link these to the subset of mutations that define the variants of concern (VOCs). First, we examined the global SARS-CoV-2 genomes and associated metadata to determine how viral properties and public health measures have influenced the magnitude of waves, as measured by the number of infection cases, in different geographic locations using regression models. This analysis showed that, as expected, both public health measures and virus properties were associated with the waves of regional SARS-CoV-2 reported infection numbers and this impact varies geographically. We attribute this to intrinsic differences such as vaccine coverage, testing and sequencing capacity and the effectiveness of government stringency. To assess underlying evolutionary change, we used non-negative matrix factorisation and observed three distinct mutational signatures, unique in their substitution patterns and exposures from the SARS-CoV-2 genomes. Signatures 1, 2 and 3 were biased to C→T, T→C/A→G and G→T point mutations. We hypothesise assignments of these mutational signatures to the host antiviral molecules APOBEC, ADAR and ROS respectively. We observe a shift amidst the pandemic in relative mutational signature activity from predominantly Signature 1 changes to an increasingly high proportion of changes consistent with Signature 2. This could represent changes in how the virus and the host immune response interact and indicates how SARS-CoV-2 may continue to generate variation in the future. Linkage of the detected mutational signatures to the VOC-defining amino acids substitutions indicates the majority of SARS-CoV-2's evolutionary capacity is likely to be associated with the action of host antiviral molecules rather than virus replication errors.

SARS-CoV-2 Omicron BA.5 Infections in Vaccinated Persons, Rural Uganda.

We describe a cluster of COVID-19 breakthrough infections after vaccination in Kyamulibwa, Kalungu District, Uganda. All but 1 infection were from SARS-CoV-2 Omicron strain BA.5.2.1. We identified 6 distinct genotypes by genome sequencing. Infections were mild, suggesting vaccination is not protective against infection but may limit disease severity.

Virus genomes reveal factors that spread and sustained the Ebola epidemic.

The 2013-2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic 'gravity' model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.

Automated download and clean-up of family-specific databases for kmer-based virus identification.

SUMMARY: Here, we present an automated pipeline for Download Of NCBI Entries (DONE) and continuous updating of a local sequence database based on user-specified queries. The database can be created with either protein or nucleotide sequences containing all entries or complete genomes only. The pipeline can automatically clean the database by removing entries with matches to a database of user-specified sequence contaminants. The default contamination entries include sequences from the UniVec database of plasmids, marker genes and sequencing adapters from NCBI, an E.coli genome, rRNA sequences, vectors and satellite sequences. Furthermore, duplicates are removed and the database is automatically screened for sequences from green fluorescent protein, luciferase and antibiotic resistance genes that might be present in some GenBank viral entries, and could lead to false positives in virus identification. For utilizing the database, we present a useful opportunity for dealing with possible human contamination. We show the applicability of DONE by downloading a virus database comprising 37 virus families. We observed an average increase of 16 776 new entries downloaded per month for the 37 families. In addition, we demonstrate the utility of a custom database compared to a standard reference database for classifying both simulated and real sequence data. AVAILABILITYAND IMPLEMENTATION: The DONE pipeline for downloading and cleaning is deposited in a publicly available repository (https://bitbucket.org/genomicepidemiology/done/src/master/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SARS-CoV-2 Variants, South Sudan, January-March 2021.

As the coronavirus pandemic continues, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequence data are required to inform vaccine efforts. We provide SARS-CoV-2 sequence data from South Sudan and document the dominance of SARS-CoV-2 lineage B.1.525 (Eta variant) during the country's second wave of infection.

Correction to: Whole genome sequencing and phylogenetic analysis of human metapneumovirus strains from Kenya and Zambia.

Following the publication of this article [1], it was noted that due to a typesetting error the figure legends were paired incorrectly.

Whole genome sequencing and phylogenetic analysis of human metapneumovirus strains from Kenya and Zambia.

BACKGROUND: Human metapneumovirus (HMPV) is an important cause of acute respiratory illness in young children. Whole genome sequencing enables better identification of transmission events and outbreaks, which is not always possible with sub-genomic sequences. RESULTS: We report a 2-reaction amplicon-based next generation sequencing method to determine the complete genome sequences of five HMPV strains, representing three subgroups (A2, B1 and B2), directly from clinical samples. In addition to reporting five novel HMPV genomes from Africa we examined genetic diversity and sequence patterns of publicly available HMPV genomes. We found that the overall nucleotide sequence identity was 71.3 and 80% for HMPV group A and B, respectively, the diversity between HMPV groups was greater at amino acid level for SH and G surface protein genes, and multiple subgroups co-circulated in various countries. Comparison of sequences between HMPV groups revealed variability in G protein length (219 to 241 amino acids) due to changes in the stop codon position. Genome-wide phylogenetic analysis showed congruence with the individual gene sequence sets except for F and M2 genes. CONCLUSION: This is the first genomic characterization of HMPV genomes from African patients.

Main Routes of Entry and Genomic Diversity of SARS-CoV-2, Uganda.

We established rapid local viral sequencing to document the genomic diversity of severe acute respiratory syndrome coronavirus 2 entering Uganda. Virus lineages closely followed the travel origins of infected persons. Our sequence data provide an important baseline for tracking any further transmission of the virus throughout the country and region.

Genomic sequence of yellow fever virus from a Dutch traveller returning from the Gambia-Senegal region, the Netherlands, November 2018.

In November 2018, yellow fever was diagnosed in a Dutch traveller returning from a bicycle tour in the Gambia-Senegal region. A complete genome sequence of yellow fever virus (YFV) from the case was generated and clustered phylogenetically with YFV from the Gambia and Senegal, ruling out importation into the Netherlands from recent outbreaks in Brazil or Angola. We emphasise the need for increased public awareness of YFV vaccination before travelling to endemic countries.

Whole-Genome Next-Generation Sequencing to Study Within-Host Evolution of Norovirus (NoV) Among Immunocompromised Patients With Chronic NoV Infection.

Background: The genus Norovirus comprises large genetic diversity, and new GII.4 variants emerge every 2-3 years. It is unknown in which host these new variants originate. Here we study whether prolonged shedders within the immunocompromised population could be a reservoir for newly emerging strains. Methods: Sixty-five fecal samples from 16 immunocompromised patients were retrospectively selected. Isolated viral RNA was enriched by hybridization with a custom norovirus whole-genome RNA bait set and deep sequenced on the Illumina MiSeq platform. Results: Patients shed virus for average 352 days (range, 76-716 days). Phylogenetic analysis showed distinct GII.4 variants in 3 of 13 patients (23%). The viral mutation rates were variable between patients but did not differ between various immune status groups. All within-host GII.4 viral populations showed amino acid changes at blocking epitopes over time, and the majority of VP1 amino acid mutations were located at the capsid surface. Conclusions: This study found viruses in immunocompromised hosts that are genetically distinct from viruses circulating in the general population, and these patients therefore may contain a reservoir for newly emerging strains. Future studies need to determine whether these new strains are of risk to other immunocompromised patients and the general population.

Severe acute respiratory infection caused by swine influenza virus in a child necessitating extracorporeal membrane oxygenation (ECMO), the Netherlands, October 2016.

In October 2016, a severe infection with swine influenza A(H1N1) virus of the Eurasian avian lineage occurred in a child with a previous history of eczema in the Netherlands, following contact to pigs. The patient's condition deteriorated rapidly and required life support through extracorporeal membrane oxygenation. After start of oseltamivir treatment and removal of mucus plugs, the patient fully recovered. Monitoring of more than 80 close unprotected contacts revealed no secondary cases.

Deep sequencing of norovirus genomes defines evolutionary patterns in an urban tropical setting.

UNLABELLED: Norovirus is a highly transmissible infectious agent that causes epidemic gastroenteritis in susceptible children and adults. Norovirus infections can be severe and can be initiated from an exceptionally small number of viral particles. Detailed genome sequence data are useful for tracking norovirus transmission and evolution. To address this need, we have developed a whole-genome deep-sequencing method that generates entire genome sequences from small amounts of clinical specimens. This novel approach employs an algorithm for reverse transcription and PCR amplification primer design using all of the publically available norovirus sequence data. Deep sequencing and de novo assembly were used to generate norovirus genomes from a large set of diarrheal patients attending three hospitals in Ho Chi Minh City, Vietnam, over a 2.5-year period. Positive-selection analysis and direct examination of protein changes in the virus over time identified codons in the regions encoding proteins VP1, p48 (NS1-2), and p22 (NS4) under positive selection and expands the known targets of norovirus evolutionary pressure. IMPORTANCE: The high transmissibility and rapid evolutionary rate of norovirus, combined with a short-lived host immune responses, are thought to be the reasons why the virus causes the majority of pediatric viral diarrhea cases. The evolutionary patterns of this RNA virus have been described in detail for only a portion of the virus genome and never for a virus from a detailed urban tropical setting. We provide a detailed sequence description of the noroviruses circulating in three Ho Chi Minh City hospitals over a 2.5-year period. This study identified patterns of virus change in known sites of host immune response and identified three additional regions of the virus genome under selection that were not previously recognized. In addition, the method described here provides a robust full-genome sequencing platform for community-based virus surveillance.

Evolution of increased positive charge on the SARS-CoV-2 spike protein may be adaptation to human transmission.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve and infect individuals. The exterior surface of the SARS-CoV-2 virion is dominated by the spike protein, and the current work examined spike protein biochemical features that have changed during the 3 years in which SARS-CoV-2 has infected humans. Our analysis identified a striking change in spike protein charge, from -8.3 in the original Lineage A and B viruses to -1.26 in most of the current Omicron viruses. We conclude that in addition to immune selection pressure, the evolution of SARS-CoV-2 has also altered viral spike protein biochemical properties, which may influence virion survival and promote transmission. Future vaccine and therapeutic development should also exploit and target these biochemical properties.

Dehalococcoides mccartyi strain CBDB1 takes up protons from the cytoplasm to reductively dehalogenate organohalides indicating a new modus of proton motive force generation.

Proton translocation across the cytoplasmic membrane is a vital process for all organisms. Dehalococcoides strains are strictly anaerobic organohalide respiring bacteria that lack quinones and cytochromes but express a large membrane-bound protein complex (OHR complex) proposed to generate a proton gradient. However, its functioning is unclear. By using a dehalogenase-based enzyme activity assay with deuterium-labelled water in various experimental designs, we obtained evidence that the halogen atom of the halogenated electron acceptor is substituted with a proton from the cytoplasm. This suggests that the protein complex couples exergonic electron flux through the periplasmic subunits of the OHR complex to the endergonic transport of protons from the cytoplasm across the cytoplasmic membrane against the proton gradient to the halogenated electron acceptor. Using computational tools, we located two proton-conducting half-channels in the AlphaFold2-predicted structure of the OmeB subunit of the OHR complex, converging in a highly conserved arginine residue that could play a proton gatekeeper role. The cytoplasmic proton half-channel in OmeB is connected to a putative proton-conducting path within the reductive dehalogenase subunit. Our results indicate that the reductive dehalogenase and its halogenated substrate serve as both electron and proton acceptors, providing insights into the proton translocation mechanism within the OHR complex and contributing to a better understanding of energy conservation in D. mccartyi strains. Our results reveal a very simple mode of energy conservation in anaerobic bacteria, showing that proton translocation coupled to periplasmic electron flow might have importance also in other microbial processes and biotechnological applications.

Viral epidemic preparedness: a perspective from five clinical microbiology laboratories in Europe.

BACKGROUND: Pandemic preparedness is critical to respond effectively to existing and emerging/new viral pathogens. Important lessons have been learned during the last pandemic at various levels. This revision discusses some of the major challenges and potential ways to address them in the likely event of future pandemics. OBJECTIVES: To identify critical points of readiness that may help us accelerate the response to future pandemics from a clinical microbiology laboratory perspective with a focus on viral diagnostics and genomic sequencing. The potential areas of improvement identified are discussed from the sample collection to information reporting. SOURCES: Microbiologists and researchers from five countries reflect on challenges encountered during the COVID-19 pandemic, review published literature on prior and current pandemics, and suggest potential solutions in preparation for future outbreaks. CONTENT: Major challenges identified in the pre-analytic and post-analytic phases from sample collection to result reporting are discussed. From the perspective of clinical microbiology laboratories, the preparedness for a new pandemic should focus on zoonotic viruses. Laboratory readiness for scalability is critical and should include elements related to material procurement, training personnel, specific funding programmes, and regulatory issues to rapidly implement "in-house" tests. Laboratories across various countries should establish (or re-use) operational networks to communicate to respond effectively, ensuring the presence of agile circuits with full traceability of samples. IMPLICATIONS: Laboratory preparedness is paramount to respond effectively to emerging and re-emerging viral infections and to limit the clinical and societal impact of new potential pandemics. Agile and fully traceable methods for sample collection to report are the cornerstone of a successful response. Expert group communication and early involvement of information technology personnel are critical for preparedness. A specific budget for pandemic preparedness should be ring-fenced and added to the national health budgets.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8.

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

Near-Complete Genome Sequences of Measles Virus Strains from 10 Years of Uganda Country-wide Surveillance.

Measles remains a global health challenge despite the availability of a safe and effective vaccine. Sporadic outbreaks of measles virus infections continue in Uganda. We report eight near-complete genome sequences of measles virus strains from Uganda cases from 2011 to 2020, providing useful data for assessing vaccine escape and local/regional transmission.

Identification of missed viruses by metagenomic sequencing of clinical respiratory samples from Kenya.

Pneumonia remains a major cause of mortality and morbidity. Most molecular diagnoses of viruses rely on polymerase chain reaction (PCR) assays that however can fail due to primer mismatch. We investigated the performance of routine virus diagnostics in Kilifi, Kenya, using random-primed viral next generation sequencing (viral NGS) on respiratory samples which tested negative for the common viral respiratory pathogens by a local standard diagnostic panel. Among 95 hospitalised pneumonia patients and 95 household-cohort individuals, analysis of viral NGS identified at least one respiratory-associated virus in 35 (37%) and 23 (24%) samples, respectively. The majority (66%; 42/64) belonged to the Picornaviridae family. The NGS data analysis identified a number of viruses that were missed by the diagnostic panel (rhinovirus, human metapneumovirus, respiratory syncytial virus and parainfluenza virus), and these failures could be attributed to PCR primer/probe binding site mismatches. Unexpected viruses identified included parvovirus B19, enterovirus D68, coxsackievirus A16 and A24 and rubella virus. The regular application of such viral NGS could help evaluate assay performance, identify molecular causes of missed diagnoses and reveal gaps in the respiratory virus set used for local screening assays. The results can provide actionable information to improve the local pneumonia diagnostics and reveal locally important viral pathogens.

Spike Protein Cleavage-Activation in the Context of the SARS-CoV-2 P681R Mutation: an Analysis from Its First Appearance in Lineage A.23.1 Identified in Uganda.

Based on its predicted ability to affect transmissibility and pathogenesis, surveillance studies have highlighted the role of a specific mutation (P681R) in the S1/S2 furin cleavage site of the SARS-CoV-2 spike protein. Here we analyzed A.23.1, first identified in Uganda, as a P681R-containing virus several months prior to the emergence of B.1.617.2 (Delta variant). We performed assays using peptides mimicking the S1/S2 from A.23.1 and B.1.617 and observed significantly increased cleavability with furin compared to both an original B lineage (Wuhan-Hu1) and B.1.1.7 (Alpha variant). We also performed cell-cell fusion and functional infectivity assays using pseudotyped particles and observed an increase in activity for A.23.1 compared to an original B lineage spike. However, these changes in activity were not reproduced in the B lineage spike bearing only the P681R substitution. Our findings suggest that while A.23.1 has increased furin-mediated cleavage linked to the P681R substitution, this substitution needs to occur on the background of other spike protein changes to enable its functional consequences. IMPORTANCE During the course of the SARS-CoV-2 pandemic, viral variants have emerged that often contain notable mutations in the spike gene. Mutations that encode changes in the spike S1/S2 (furin) activation site have been considered especially impactful. The S1/S2 change from proline to arginine at position 681 (P681R) first emerged in the A.23.1 variant in Uganda, and subsequently occurred in the more widely transmitted Delta variant. We show that the A.23.1 spike is more readily activated by the host cell protease furin, but that this is not reproduced in an original SARS-CoV-2 spike containing the P681R mutation. Changes to the S1/S2 (furin) activation site play a role in SARS-CoV-2 infection and spread, but successful viruses combine these mutations with other less well identified changes, occurring as part of natural selection.

Association of envelope-specific B-cell differentiation and viral selective pressure signatures in HIV-1 CRF01_AE infection.

OBJECTIVE: In HIV type 1 (HIV-1) infection, virus-specific B-cell and neutralizing antibody (NAb) responses are impaired but exert selective pressure on target viral Envelope (Env) resulting in prominent sequence diversification among geographical areas. The basal induction patterns of HIV Env-specific B cells and their interaction with HIV Env awaits clarification. DESIGN: We investigated the relationship of Env polymorphisms and Env-specific B-cell responses in treatment-naive HIV-1 CRF01_AE-infected Vietnamese. METHODS: Samples of 43 HIV-1 CRF01_AE infection-identified individuals were divided into acute-phase ( n  = 12) and chronic-phase ( n  = 31) by combined criteria of serological recent-infection assay and clinical parameters. We quantified subcloning-based polymorphic residue site numbers in plasma-derived Env variable region 1-5 (V1-V5)-coding regions within each individual, designating their summation within each region as variant index. Peripheral blood Env gp 140-specific B-cell responses and plasma neutralizing activity of Env pseudoviruses were examined to analyze their relationship with variant index. RESULTS: HIV-1 CRF01_AE Env gp140-specific total B-cell and plasma cell (CD19 + IgD - CD27 + CD38 + CD138 + ) responses were determined. In chronic-phase samples, significant correlation of variant index in all Env V1-V5 regions with Env-specific plasma cell responses was shown, and V1-V5 total variant index correlated stronger with Env-specific plasma cell as compared with total Env-specific B-cell responses. Env V5 variant index was significantly higher in chronic-phase cross-neutralizers of V5-polymorphic/VRC01-insensitive CRF01_AE Env. CONCLUSION: Results revealed the association between circulating Env-specific plasma cell responses and Env polymorphisms, implicating selective pressure on Env by plasma cell-derived antibodies and conversely suggests that Env-specific B-cell induction alone is insufficient for exerting Env selective pressure in HIV infection.