Genome-Scale Phylogeny and Evolutionary Analysis of Ross River Virus Reveals Periodic Sweeps of Lineage Dominance in Western Australia, 1977-2014

Genomic Analysis of Sindbis Virus Reveals Uncharacterized Diversity within the Australasian Region, and Support for Revised SINV Taxonomy

Sindbis virus (SINV) is a widely dispersed mosquito-borne alphavirus. Reports of Sindbis disease are largely restricted to northern Europe and South Africa. SINV is frequently sampled in Australian mosquito-based arbovirus surveillance programs, but human disease has rarely been reported. Molecular epidemiological studies have characterized six SINV genotypes (G1-G6) based on E2 gene phylogenies, mostly comprising viruses derived from the African-European zoogeographical region and with limited representation of Australasian SINV. In this study, we conducted whole genome sequencing of 66 SINV isolates sampled between 1960 and 2014 from countries of the Australasian region: Australia, Malaysia, and Papua New Guinea. G2 viruses were the most frequently and widely sampled, with three distinct sub-lineages defined. No new G6 SINV were identified, confirming geographic restriction of these viruses to south-western Australia. Comparison with global SINV characterized large-scale nucleotide and amino acid sequence divergence between African-European G1 viruses and viruses that circulate in Australasia (G2 and G3) of up to 26.83% and 14.55%, respectively, divergence that is sufficient for G2/G3 species demarcation. We propose G2 and G3 are collectively a single distinct alphavirus species that we name Argyle virus, supported by the inapparent or mild disease phenotype and the higher evolutionary rate compared with G1. Similarly, we propose G6, with 24.7% and 12.61% nucleotide and amino acid sequence divergence, is a distinct alphavirus species that we name Thomson's Lake virus.

Susceptibility of Aedes albopictus, Ae. aegypti and human populations to Ross River virus in Kuala Lumpur, Malaysia

BACKGROUND

Emerging arboviruses such as chikungunya and Zika viruses have unexpectedly caused widespread outbreaks in tropical and subtropical regions recently. Ross River virus (RRV) is endemic in Australia and has epidemic potential. In Malaysia, Aedes mosquitoes are abundant and drive dengue and chikungunya outbreaks. We assessed risk of an RRV outbreak in Kuala Lumpur, Malaysia by determining vector competence of local Aedes mosquitoes and local seroprevalence as a proxy of human population susceptibility.

METHODOLOGY/PRINCIPAL FINDINGS

We assessed oral susceptibility of Malaysian Ae. aegypti and Ae. albopictus by real-time PCR to an Australian RRV strain SW2089. Replication kinetics in midgut, head and saliva were determined at 3 and 10 days post-infection (dpi). With a 3 log10 PFU/ml blood meal, infection rate was higher in Ae. albopictus (60%) than Ae. aegypti (15%; p<0.05). Despite similar infection rates at 5 and 7 log10 PFU/ml blood meals, Ae. albopictus had significantly higher viral loads and required a significantly lower median oral infectious dose (2.7 log10 PFU/ml) than Ae. aegypti (4.2 log10 PFU/ml). Ae. albopictus showed higher vector competence, with higher viral loads in heads and saliva, and higher transmission rate (RRV present in saliva) of 100% at 10 dpi, than Ae. aegypti (41%). Ae. aegypti demonstrated greater barriers at either midgut escape or salivary gland infection, and salivary gland escape. We then assessed seropositivity against RRV among 240 Kuala Lumpur inpatients using plaque reduction neutralization, and found a low rate of 0.8%.

CONCLUSIONS/SIGNIFICANCE

Both Ae. aegypti and Ae. albopictus are susceptible to RRV, but Ae. albopictus displays greater vector competence. Extensive travel links with Australia, abundant Aedes vectors, and low population immunity places Kuala Lumpur, Malaysia at risk of an imported RRV outbreak. Surveillance and increased diagnostic awareness and capacity are imperative to prevent establishment of new arboviruses in Malaysia.

Fine-scale genomic tracking of Ross River virus using nanopore sequencing

BACKGROUND

Ross River virus (RRV) is Australia's most common and widespread mosquito-transmitted arbovirus and is of significant public health concern. With increasing anthropogenic impacts on wildlife and mosquito populations, it is important that we understand how RRV circulates in its endemic hotspots to determine where public health efforts should be directed. Current surveillance methods are effective in locating the virus but do not provide data on the circulation of the virus and its strains within the environment. This study examined the ability to identify single nucleotide polymorphisms (SNPs) within the variable E2/E3 region by generating full-length haplotypes from a range of mosquito trap-derived samples.

METHODS

A novel tiled primer amplification workflow for amplifying RRV was developed with analysis using Oxford Nanopore Technology's MinION and a custom ARTIC/InterARTIC bioinformatic protocol. By creating a range of amplicons across the whole genome, fine-scale SNP analysis was enabled by specifically targeting the variable region that was amplified as a single fragment and established haplotypes that informed spatial-temporal variation of RRV in the study site in Victoria.

RESULTS

A bioinformatic and laboratory pipeline was successfully designed and implemented on mosquito whole trap homogenates. Resulting data showed that genotyping could be conducted in real time and that whole trap consensus of the viruses (with major SNPs) could be determined in a timely manner. Minor variants were successfully detected from the variable E2/E3 region of RRV, which allowed haplotype determination within complex mosquito homogenate samples.

CONCLUSIONS

The novel bioinformatic and wet laboratory methods developed here will enable fast detection and characterisation of RRV isolates. The concepts presented in this body of work are transferable to other viruses that exist as quasispecies in samples. The ability to detect minor SNPs, and thus haplotype strains, is critically important for understanding the epidemiology of viruses their natural environment.

Early Genomic Surveillance and Phylogeographic Analysis of Getah Virus, a Reemerging Arbovirus, in Livestock in China

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

Enhanced Arbovirus Surveillance with High-Throughput Metatranscriptomic Processing of Field-Collected Mosquitoes

Surveillance programs are essential for the prevention and control of mosquito-borne arboviruses that cause serious human and animal diseases. Viral metatranscriptomic sequencing can enhance surveillance by enabling untargeted, high-throughput arbovirus detection. We used metatranscriptomic sequencing to screen field-collected mosquitoes for arboviruses to better understand how metatranscriptomics can be utilised in routine surveillance. Following a significant flood event in 2016, more than 56,000 mosquitoes were collected over seven weeks from field traps set up in Victoria, Australia. The traps were split into samples of 1000 mosquitoes or less and sequenced on the Illumina HiSeq. Five arboviruses relevant to public health (Ross River virus, Sindbis virus, Trubanaman virus, Umatilla virus, and Wongorr virus) were detected a total of 33 times in the metatranscriptomic data, with 94% confirmed using reverse transcription quantitative PCR (RT-qPCR). Analysis of metatranscriptomic cytochrome oxidase I (COI) sequences enabled the detection of 12 mosquito and two biting midge species. Screening of the same traps by an established public health arbovirus surveillance program corroborated the metatranscriptomic arbovirus and mosquito species detections. Assembly of genome sequences from the metatranscriptomic data also led to the detection of 51 insect-specific viruses, both known and previously undescribed, and allowed phylogenetic comparison to past strains. We have demonstrated how metatranscriptomics can enhance surveillance by enabling untargeted arbovirus detection, providing genomic epidemiological data, and simultaneously identifying vector species from large, unsorted mosquito traps.

The Putative Roles and Functions of Indel, Repetition and Duplication Events in Alphavirus Non-Structural Protein 3 Hypervariable Domain (nsP3 HVD) in Evolution, Viability and Re-Emergence

Alphavirus non-structural proteins 1–4 (nsP1, nsP2, nsP3, and nsP4) are known to be crucial for alphavirus RNA replication and translation. To date, nsP3 has been demonstrated to mediate many virus–host protein–protein interactions in several fundamental alphavirus mechanisms, particularly during the early stages of replication. However, the molecular pathways and proteins networks underlying these mechanisms remain poorly described. This is due to the low genetic sequence homology of the nsP3 protein among the alphavirus species, especially at its 3′ C-terminal domain, the hypervariable domain (HVD). Moreover, the nsP3 HVD is almost or completely intrinsically disordered and has a poor ability to form secondary structures. Evolution in the nsP3 HVD region allows the alphavirus to adapt to vertebrate and insect hosts. This review focuses on the putative roles and functions of indel, repetition, and duplication events that have occurred in the alphavirus nsP3 HVD, including characterization of the differences and their implications for specificity in the context of virus–host interactions in fundamental alphavirus mechanisms, which have thus directly facilitated the evolution, adaptation, viability, and re-emergence of these viruses.

Genome Sequence Analysis of First Ross River Virus Isolate from Papua New Guinea Indicates Long-Term, Local Evolution

Ross River virus (RRV) is the most medically significant mosquito-borne virus of Australia, in terms of human morbidity. RRV cases, characterised by febrile illness and potentially persistent arthralgia, have been reported from all Australian states and territories. RRV was the cause of a large-scale epidemic of multiple Pacific Island countries and territories (PICTs) from 1979 to 1980, involving at least 50,000 cases. Historical evidence of RRV seropositivity beyond Australia, in populations of Papua New Guinea (PNG), Indonesia and the Solomon Islands, has been documented. We describe the genomic characterisation and timescale analysis of the first isolate of RRV to be sampled from PNG to date. Our analysis indicates that RRV has evolved locally within PNG, independent of Australian lineages, over an approximate 40 year period. The mean time to most recent common ancestor (tMRCA) of the unique PNG clade coincides with the initiation of the PICTs epidemic in mid-1979. This may indicate that an ancestral variant of the PNG clade was seeded into the region during the epidemic, a period of high RRV transmission. Further epidemiological and molecular-based surveillance is required in PNG to better understand the molecular epidemiology of RRV in the general Australasian region.

V. R, Jha RM, Kumari A, Gandotra S, Rao V, Faruq M, Kumar S, Reshma G. B, Varma G. N, Roy SS, Sengupta A, Chattopadhyay S, Singhal K, Pradhan S, Jha D, Naushin S, Wadhwa S, Tyagi N, Poojary M, Scaria V, Sivasubbu S. High throughput detection and genetic epidemiology of SARS-CoV-2 using COVIDSeq next-generation sequencing

The rapid emergence of coronavirus disease 2019 (COVID-19) as a global pandemic affecting millions of individuals globally has necessitated sensitive and high-throughput approaches for the diagnosis, surveillance, and determining the genetic epidemiology of SARS-CoV-2. In the present study, we used the COVIDSeq protocol, which involves multiplex-PCR, barcoding, and sequencing of samples for high-throughput detection and deciphering the genetic epidemiology of SARS-CoV-2. We used the approach on 752 clinical samples in duplicates, amounting to a total of 1536 samples which could be sequenced on a single S4 sequencing flow cell on NovaSeq 6000. Our analysis suggests a high concordance between technical duplicates and a high concordance of detection of SARS-CoV-2 between the COVIDSeq as well as RT-PCR approaches. An in-depth analysis revealed a total of six samples in which COVIDSeq detected SARS-CoV-2 in high confidence which were negative in RT-PCR. Additionally, the assay could detect SARS-CoV-2 in 21 samples and 16 samples which were classified inconclusive and pan-sarbeco positive respectively suggesting that COVIDSeq could be used as a confirmatory test. The sequencing approach also enabled insights into the evolution and genetic epidemiology of the SARS-CoV-2 samples. The samples were classified into a total of 3 clades. This study reports two lineages B.1.112 and B.1.99 for the first time in India. This study also revealed 1,143 unique single nucleotide variants and added a total of 73 novel variants identified for the first time. To the best of our knowledge, this is the first report of the COVIDSeq approach for detection and genetic epidemiology of SARS-CoV-2. Our analysis suggests that COVIDSeq could be a potential high sensitivity assay for the detection of SARS-CoV-2, with an additional advantage of enabling the genetic epidemiology of SARS-CoV-2.

A cross-reactive antibody protects against Ross River virus musculoskeletal disease despite rapid neutralization escape in mice

Arthritogenic alphaviruses cause debilitating musculoskeletal disease and historically have circulated in distinct regions. With the global spread of chikungunya virus (CHIKV), there now is more geographic overlap, which could result in heterologous immunity affecting natural infection or vaccination. Here, we evaluated the capacity of a cross-reactive anti-CHIKV monoclonal antibody (CHK-265) to protect against disease caused by the distantly related alphavirus, Ross River virus (RRV). Although CHK-265 only moderately neutralizes RRV infection in cell culture, it limited clinical disease in mice independently of Fc effector function activity. Despite this protective phenotype, RRV escaped from CHK-265 neutralization in vivo, with resistant variants retaining pathogenic potential. Near the inoculation site, CHK-265 reduced viral burden in a type I interferon signaling-dependent manner and limited immune cell infiltration into musculoskeletal tissue. In a parallel set of experiments, purified human CHIKV immune IgG also weakly neutralized RRV, yet when transferred to mice, resulted in improved clinical outcome during RRV infection despite the emergence of resistant viruses. Overall, this study suggests that weakly cross-neutralizing antibodies can protect against heterologous alphavirus disease, even if neutralization escape occurs, through an early viral control program that tempers inflammation.

The induction of broadly neutralizing antibodies is a goal of many antiviral vaccine programs. In this study, we show that cross-reactive monoclonal and polyclonal antibodies developed after CHIKV infection or immunization with relatively weak cross-neutralizing activity can protect against RRV-induced musculoskeletal disease in mice. Even though RRV rapidly escaped from neutralization, antibody therapy reduced inflammation in musculoskeletal tissues and decreased viral burden near the site of infection in a manner that required type I interferon signaling. These studies in mice show that broadly reactive antibodies with limited neutralizing activity still can confer protection against heterologous alphaviruses.

Phylogenetic and Timescale Analysis of Barmah Forest Virus as Inferred from Genome Sequence Analysis

Barmah Forest virus (BFV) is a medically important mosquito-borne alphavirus endemic to Australia. Symptomatic disease can be a major cause of morbidity, associated with fever, rash, and debilitating arthralgia. BFV disease is similar to that caused by Ross River virus (RRV), the other major Australian alphavirus. Currently, just four BFV whole-genome sequences are available with no genome-scale phylogeny in existence to robustly characterise genetic diversity. Thirty novel genome sequences were derived for this study, for a final 34-taxon dataset sampled over a 44 year period. Three distinct BFV genotypes were characterised (G1–3) that have circulated in Australia and Papua New Guinea (PNG). Evidence of spatio-temporal co-circulation of G2 and G3 within regions of Australia was noted, including in the South West region of Western Australia (WA) during the first reported disease outbreaks in the state’s history. Compared with RRV, the BFV population appeared more stable with less frequent emergence of novel lineages. Preliminary in vitro assessment of RRV and BFV replication kinetics found that RRV replicates at a significantly faster rate and to a higher, more persistent titre compared with BFV, perhaps indicating mosquitoes may be infectious with RRV for longer than with BFV. This investigation resolved a greater diversity of BFV, and a greater understanding of the evolutionary dynamics and history was attained.