Epigenetic age estimation of wild mice using faecal samples.

Age is a key parameter in population ecology, with a myriad of biological processes changing with age as organisms develop in early life then later senesce. As age is often hard to accurately measure with non-lethal methods, epigenetic methods of age estimation (epigenetic clocks) have become a popular tool in animal ecology and are often developed or calibrated using captive animals of known age. However, studies typically rely on invasive blood or tissue samples, which limit their application in more sensitive or elusive species. Moreover, few studies have directly assessed how methylation patterns and epigenetic age estimates compare across environmental contexts (e.g. captive or laboratory-based vs. wild animals). Here, we built a targeted epigenetic clock from laboratory house mice (strain C57BL/6, Mus musculus) using DNA from non-invasive faecal samples, and then used it to estimate age in a population of wild mice (Mus musculus domesticus) of unknown age. This laboratory mouse-derived epigenetic clock accurately predicted adult wild mice to be older than juveniles and showed that wild mice typically increased in epigenetic age over time, but with wide variation in epigenetic ageing rate among individuals. Our results also suggested that, for a given body mass, wild mice had higher methylation across targeted CpG sites than laboratory mice (and consistently higher epigenetic age estimates as a result), even among the smallest, juvenile mice. This suggests wild and laboratory mice may display different CpG methylation levels from very early in life and indicates caution is needed when developing epigenetic clocks on laboratory animals and applying them in the wild.

Automated detection and classification of polioviruses from nanopore sequencing reads using piranha.

Widespread surveillance, rapid detection, and appropriate intervention will be critical for successful eradication of poliovirus. Using deployable next-generation sequencing (NGS) approaches, such as Oxford Nanopore Technologies' MinION, the time from sample to result can be significantly reduced compared to cell culture and Sanger sequencing. We developed piranha (poliovirus investigation resource automating nanopore haplotype analysis), a 'sequencing reads-to-report' solution to aid routine poliovirus testing of both stool and environmental samples and alleviate the bioinformatic bottleneck that often exists for laboratories adopting novel NGS approaches. Piranha can be used for efficient intratypic differentiation of poliovirus serotypes, for classification of Sabin-like polioviruses, and for detection of wild-type and vaccine-derived polioviruses. It produces interactive, distributable reports, as well as summary comma-separated values files and consensus poliovirus FASTA sequences. Piranha optionally provides phylogenetic analysis, with the ability to incorporate a local database, processing from raw sequencing reads to an interactive, annotated phylogeny in a single step. The reports describe each nanopore sequencing run with interpretable plots, enabling researchers to easily detect the presence of poliovirus in samples and quickly disseminate their results. Poliovirus eradication efforts are hindered by the lack of real-time detection and reporting, and piranha can be used to complement direct detection sequencing approaches.

SARS-CoV-2 lineage assignments using phylogenetic placement/UShER are superior to pangoLEARN machine-learning method.

With the rapid spread and evolution of SARS-CoV-2, the ability to monitor its transmission and distinguish among viral lineages is critical for pandemic response efforts. The most commonly used software for the lineage assignment of newly isolated SARS-CoV-2 genomes is pangolin, which offers two methods of assignment, pangoLEARN and pUShER. PangoLEARN rapidly assigns lineages using a machine-learning algorithm, while pUShER performs a phylogenetic placement to identify the lineage corresponding to a newly sequenced genome. In a preliminary study, we observed that pangoLEARN (decision tree model), while substantially faster than pUShER, offered less consistency across different versions of pangolin v3. Here, we expand upon this analysis to include v3 and v4 of pangolin, which moved the default algorithm for lineage assignment from pangoLEARN in v3 to pUShER in v4, and perform a thorough analysis confirming that pUShER is not only more stable across versions but also more accurate. Our findings suggest that future lineage assignment algorithms for various pathogens should consider the value of phylogenetic placement.

A framework for automated scalable designation of viral pathogen lineages from genomic data.

Pathogen lineage nomenclature systems are a key component of effective communication and collaboration for researchers and public health workers. Since February 2021, the Pango dynamic lineage nomenclature for SARS-CoV-2 has been sustained by crowdsourced lineage proposals as new isolates were sequenced. This approach is vulnerable to time-critical delays as well as regional and personal bias. Here we developed a simple heuristic approach for dividing phylogenetic trees into lineages, including the prioritization of key mutations or genes. Our implementation is efficient on extremely large phylogenetic trees consisting of millions of sequences and produces similar results to existing manually curated lineage designations when applied to SARS-CoV-2 and other viruses including chikungunya virus, Venezuelan equine encephalitis virus complex and Zika virus. This method offers a simple, automated and consistent approach to pathogen nomenclature that can assist researchers in developing and maintaining phylogeny-based classifications in the face of ever-increasing genomic datasets.

APOBEC3 deaminase editing in mpox virus as evidence for sustained human transmission since at least 2016.

Historically, mpox has been characterized as an endemic zoonotic disease that transmits through contact with the reservoir rodent host in West and Central Africa. However, in May 2022, human cases of mpox were detected spreading internationally beyond countries with known endemic reservoirs. When the first cases from 2022 were sequenced, they shared 42 nucleotide differences from the closest mpox virus (MPXV) previously sampled. Nearly all these mutations are characteristic of the action of APOBEC3 deaminases, host enzymes with antiviral function. Assuming APOBEC3 editing is characteristic of human MPXV infection, we developed a dual-process phylogenetic molecular clock that-inferring a rate of ~6 APOBEC3 mutations per year-estimates that MPXV has been circulating in humans since 2016. These observations of sustained MPXV transmission present a fundamental shift to the perceived paradigm of MPXV epidemiology as a zoonosis and highlight the need for revising public health messaging around MPXV as well as outbreak management and control.

SARS-CoV-2 evolution in the Omicron era.

Since SARS-CoV-2 BA.5 (Omicron) emerged and spread in 2022, Omicron lineages have markedly diversified. Here we review the evolutionary trajectories and processes that underpin the emergence of these lineages, and identify the most prevalent sublineages. We discuss the potential origins of second-generation BA.2 lineages. Simple and complex recombination, antigenic drift and convergent evolution have enabled SARS-CoV-2 to accumulate mutations that alter its antigenicity. We also discuss the potential evolutionary trajectories of SARS-CoV-2 in the future.

SCORPIO: a utility for defining and classifying mutation constellations of virus genomes.

SUMMARY: Scorpio provides a set of command line utilities for classifying, haplotyping, and defining constellations of mutations for an aligned set of genome sequences. It was developed to enable exploration and classification of variants of concern within the SARS-CoV-2 pandemic, but can be applied more generally to other species. AVAILABILITY AND IMPLEMENTATION: Scorpio is an open-source project distributed under the GNU GPL version 3 license. Source code and binaries are available at https://github.com/cov-lineages/scorpio, and binaries are also available from Bioconda. SARS-CoV-2 specific definitions can be installed as a separate dependency from https://github.com/cov-lineages/constellations.

Correction: Implementation of Genomic Surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings.

[This corrects the article DOI: 10.1371/journal.pgph.0001455.].

Sensitive poliovirus detection using nested PCR and nanopore sequencing: a prospective validation study.

Timely detection of outbreaks is needed for poliovirus eradication, but gold standard detection in the Democratic Republic of the Congo takes 30 days (median). Direct molecular detection and nanopore sequencing (DDNS) of poliovirus in stool samples is a promising fast method. Here we report prospective testing of stool samples from suspected polio cases, and their contacts, in the Democratic Republic of the Congo between 10 August 2021 and 4 February 2022. DDNS detected polioviruses in 62/2,339 (2.7%) of samples, while gold standard combination of cell culture, quantitative PCR and Sanger sequencing detected polioviruses in 51/2,339 (2.2%) of the same samples. DDNS provided case confirmation in 7 days (median) in routine surveillance conditions. DDNS enabled confirmation of three serotype 2 circulating vaccine-derived poliovirus outbreaks 23 days (mean) earlier (range 6-30 days) than the gold standard method. The mean sequence similarity between sequences obtained by the two methods was 99.98%. Our data confirm the feasibility of implementing DDNS in a national poliovirus laboratory.

Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Implementation of genomic surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings.

The COVID-19 pandemic highlighted the importance of global genomic surveillance to monitor the emergence and spread of SARS-CoV-2 variants and inform public health decision-making. Until December 2020 there was minimal capacity for viral genomic surveillance in most Caribbean countries. To overcome this constraint, the COVID-19: Infectious disease Molecular epidemiology for PAthogen Control & Tracking (COVID-19 IMPACT) project was implemented to establish rapid SARS-CoV-2 whole genome nanopore sequencing at The University of the West Indies (UWI) in Trinidad and Tobago (T&T) and provide needed SARS-CoV-2 sequencing services for T&T and other Caribbean Public Health Agency Member States (CMS). Using the Oxford Nanopore Technologies MinION sequencing platform and ARTIC network sequencing protocols and bioinformatics pipeline, a total of 3610 SARS-CoV-2 positive RNA samples, received from 17 CMS, were sequenced in-situ during the period December 5th 2020 to December 31st 2021. Ninety-one Pango lineages, including those of five variants of concern (VOC), were identified. Genetic analysis revealed at least 260 introductions to the CMS from other global regions. For each of the 17 CMS, the percentage of reported COVID-19 cases sequenced by the COVID-19 IMPACT laboratory ranged from 0·02% to 3·80% (median = 1·12%). Sequences submitted to GISAID by our study represented 73·3% of all SARS-CoV-2 sequences from the 17 CMS available on the database up to December 31st 2021. Increased staffing, process and infrastructural improvement over the course of the project helped reduce turnaround times for reporting to originating institutions and sequence uploads to GISAID. Insights from our genomic surveillance network in the Caribbean region directly influenced non-pharmaceutical countermeasures in the CMS countries. However, limited availability of associated surveillance and clinical data made it challenging to contextualise the observed SARS-CoV-2 diversity and evolution, highlighting the need for development of infrastructure for collecting and integrating genomic sequencing data and sample-associated metadata.

Comparison of Eleven RNA Extraction Methods for Poliovirus Direct Molecular Detection in Stool Samples.

Direct detection by PCR of poliovirus RNA in stool samples provides a rapid diagnostic and surveillance tool that can replace virus isolation by cell culture in global polio surveillance. The sensitivity of direct detection methods is likely to depend on the choice of RNA extraction method and sample volume. We report a comparative analysis of 11 nucleic acid extraction methods (7 manual and 4 semiautomated) for poliovirus molecular detection using stool samples (n = 59) that had been previously identified as poliovirus positive by cell culture. To assess the effect of RNA recovery methods, extracted RNA using each of the 11 methods was tested with a poliovirus-specific reverse transcription-quantitative PCR (RT-qPCR), a pan-poliovirus RT-PCR (near-whole-genome amplification), a pan-enterovirus RT-PCR (entire capsid region), and a nested VP1 PCR that is the basis of a direct detection method based on nanopore sequencing. We also assessed extracted RNA integrity and quantity. The overall effect of extraction method on poliovirus PCR amplification assays tested in this study was found to be statistically significant (P < 0.001), thus indicating that the choice of RNA extraction method is an important component that needs to be carefully considered for any diagnostic based on nucleic acid amplification. Performance of the methods was generally consistent across the different assays used. Of the 11 extraction methods tested, the MagMAX viral RNA isolation kit used manually or automatically was found to be the preferable method for poliovirus molecular direct detection considering performance, cost, and processing time. IMPORTANCE Poliovirus, the causative agent of poliomyelitis, is a target of global eradication led by the World Health Organization since 1988. Direct molecular detection and genomic sequencing without virus propagation in cell culture is arguably a critical tool in the final stages of polio eradication. Efficient recovery of good-quality viral RNA from stool samples is a prerequisite for direct detection by nucleic acid amplification. We tested 11 nucleic acid extraction methods to identify those facilitating sensitive, fast, simple, and cost-effective extraction, with flexibility for manual and automated protocols considered. Several different PCR assays were used to compare the recovered viral RNA to test suitability for poliovirus direct molecular detection. Our findings highlight the importance of choosing a suitable RNA extraction protocol and provide useful information to diagnostic laboratories and researchers facing the choice of RNA extraction method for direct molecular virus detection from stool.

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK.

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organization as Alpha. Originating in early autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK and the imposition of new restrictions, in particular, the English national lockdown in November 2020. While these interventions succeeded in reducing the absolute number of cases, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of the SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically infected individual. We conclude that the latter provides the best explanation of the observed behaviour and dynamics of the variant, although the individual need not be immunocompromised, as persistently infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs and find that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations and a lack of the rapid evolutionary rate on its ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms), it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

Urgent need for a non-discriminatory and non-stigmatizing nomenclature for monkeypox virus.

We propose a novel, non-discriminatory classification of monkeypox virus diversity. Together with the World Health Organization, we named three clades (I, IIa and IIb) in order of detection. Within IIb, the cause of the current global outbreak, we identified multiple lineages (A.1, A.2, A.1.1 and B.1) to support real-time genomic surveillance.

Shotgun metagenomic sequencing of the first case of monkeypox virus in Brazil, 2022.

Monkeypox virus (MPXV), a zoonotic virus endemic to the African continent, has been reported in 33 non-endemic countries since May 2022. We report an almost complete genome of the first confirmed case of MPXV in Brazil. Shotgun metagenomic sequencing was completed in 18 hours, from DNA extraction to consensus sequence generation.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Making genomic surveillance deliver: A lineage classification and nomenclature system to inform rabies elimination.

The availability of pathogen sequence data and use of genomic surveillance is rapidly increasing. Genomic tools and classification systems need updating to reflect this. Here, rabies virus is used as an example to showcase the potential value of updated genomic tools to enhance surveillance to better understand epidemiological dynamics and improve disease control. Previous studies have described the evolutionary history of rabies virus, however the resulting taxonomy lacks the definition necessary to identify incursions, lineage turnover and transmission routes at high resolution. Here we propose a lineage classification system based on the dynamic nomenclature used for SARS-CoV-2, defining a lineage by phylogenetic methods for tracking virus spread and comparing sequences across geographic areas. We demonstrate this system through application to the globally distributed Cosmopolitan clade of rabies virus, defining 96 total lineages within the clade, beyond the 22 previously reported. We further show how integration of this tool with a new rabies virus sequence data resource (RABV-GLUE) enables rapid application, for example, highlighting lineage dynamics relevant to control and elimination programmes, such as identifying importations and their sources, as well as areas of persistence and routes of virus movement, including transboundary incursions. This system and the tools developed should be useful for coordinating and targeting control programmes and monitoring progress as countries work towards eliminating dog-mediated rabies, as well as having potential for broader application to the surveillance of other viruses.

Pango lineage designation and assignment using SARS-CoV-2 spike gene nucleotide sequences.

BACKGROUND: More than 2 million SARS-CoV-2 genome sequences have been generated and shared since the start of the COVID-19 pandemic and constitute a vital information source that informs outbreak control, disease surveillance, and public health policy. The Pango dynamic nomenclature is a popular system for classifying and naming genetically-distinct lineages of SARS-CoV-2, including variants of concern, and is based on the analysis of complete or near-complete virus genomes. However, for several reasons, nucleotide sequences may be generated that cover only the spike gene of SARS-CoV-2. It is therefore important to understand how much information about Pango lineage status is contained in spike-only nucleotide sequences. Here we explore how Pango lineages might be reliably designated and assigned to spike-only nucleotide sequences. We survey the genetic diversity of such sequences, and investigate the information they contain about Pango lineage status. RESULTS: Although many lineages, including the main variants of concern, can be identified clearly using spike-only sequences, some spike-only sequences are shared among tens or hundreds of Pango lineages. To facilitate the classification of SARS-CoV-2 lineages using subgenomic sequences we introduce the notion of designating such sequences to a "lineage set", which represents the range of Pango lineages that are consistent with the observed mutations in a given spike sequence. CONCLUSIONS: We find that many lineages, including the main variants-of-concern, can be reliably identified by spike alone and we define lineage-sets to represent the lineage precision that can be achieved using spike-only nucleotide sequences. These data provide a foundation for the development of software tools that can assign newly-generated spike nucleotide sequences to Pango lineage sets.

Genomics-informed outbreak investigations of SARS-CoV-2 using civet.

The scale of data produced during the SARS-CoV-2 pandemic has been unprecedented, with more than 13 million sequences shared publicly at the time of writing. This wealth of sequence data provides important context for interpreting local outbreaks. However, placing sequences of interest into national and international context is difficult given the size of the global dataset. Often outbreak investigations and genomic surveillance efforts require running similar analyses again and again on the latest dataset and producing reports. We developed civet (cluster investigation and virus epidemiology tool) to aid these routine analyses and facilitate virus outbreak investigation and surveillance. Civet can place sequences of interest in the local context of background diversity, resolving the query into different 'catchments' and presenting the phylogenetic results alongside metadata in an interactive, distributable report. Civet can be used on a fine scale for clinical outbreak investigation, for local surveillance and cluster discovery, and to routinely summarise the virus diversity circulating on a national level. Civet reports have helped researchers and public health bodies feedback genomic information in the appropriate context within a timeframe that is useful for public health.