Mosquitoes harbor a high diversity of RNA viruses, including many that impact human health. Despite a growing effort to describe the extent and nature of the mosquito virome, little is known about how these viruses persist, spread, and interact with both their hosts and other microbes. To address this issue we performed a metatranscriptomics analysis of 12 Western Australian mosquito populations structured by species and geographic location. Our results identified the complete genomes of 24 species of RNA viruses from a diverse range of viral families and orders, among which 19 are newly described. Comparisons of viromes revealed a striking difference between the two mosquito genera, with viromes of mosquitoes of the Aedes genus exhibiting substantially less diversity and lower abundances than those of mosquitoes of the Culex genus, within which the viral abundance reached 16.87% of the total non-rRNA. In addition, there was little overlap in viral diversity between the two genera, although the viromes were very similar among the three Culex species studied, suggesting that the host taxon plays a major role in structuring virus diversity. In contrast, we found no evidence that geographic location played a major role in shaping RNA virus diversity, and several viruses discovered here exhibited high similarity (95 to 98% nucleotide identity) to those from Indonesia and China. Finally, using abundance-level and phylogenetic relationships, we were able to distinguish potential mosquito viruses from those present in coinfecting bacteria, fungi, and protists. In sum, our metatranscriptomics approach provides important insights into the ecology of mosquito RNA viruses. IMPORTANCE Studies of virus ecology have generally focused on individual viral species. However, recent advances in bulk RNA sequencing make it possible to utilize metatranscriptomic approaches to reveal both complete virus diversity and the relative abundance of these viruses. We used such a metatranscriptomic approach to determine key aspects of the ecology of mosquito viruses in Western Australia. Our results show that RNA viruses are some of the most important components of the mosquito transcriptome, and we identified 19 new virus species from a diverse set of virus families. A key result was that host genetic background plays a more important role in shaping virus diversity than sampling location, with Culex species harboring more viruses at higher abundance than those from Aedes mosquitoes.KEYWORDS Australia, ecology, evolution, mosquito, phylogeny, transcriptomics, virome

Zika virus (ZIKV) has emerged since 2013 as a significant global human health threat following outbreaks in the Pacific Islands and rapid spread throughout South and Central America. Severe congenital and neurological sequelae have been linked to ZIKV infections. Assessing the ability of common mosquito species to transmit ZIKV and characterizing variation in mosquito transmission of different ZIKV strains is important for estimating regional outbreak potential and for prioritizing local mosquito control strategies for Aedes and Culex species. In this study, we evaluated the laboratory vector competence of Aedes aegypti, Culex quinquefasciatus, and Culex tarsalis that originated in areas of California where ZIKV cases in travelers since 2015 were frequent. We compared infection, dissemination, and transmission rates by measuring ZIKV RNA levels in cohorts of mosquitoes that ingested blood meals from type I interferon-deficient mice infected with either a Puerto Rican ZIKV strain from 2015 (PR15), a Brazilian ZIKV strain from 2015 (BR15), or an ancestral Asian-lineage Malaysian ZIKV strain from 1966 (MA66). With PR15, Cx. quinquefasciatus was refractory to infection (0%, N = 42) and Cx. tarsalis was infected at 4% (N = 46). No ZIKV RNA was detected in saliva from either Culex species 14 or 21 days post feeding (dpf). In contrast, Ae. aegypti developed infection rates of 85% (PR15; N = 46), 90% (BR15; N = 20), and 81% (MA66; N = 85) 14 or 15 dpf. Although MA66-infected Ae. aegypti showed higher levels of ZIKV RNA in mosquito bodies and legs, transmission rates were not significantly different across virus strains (P = 0.13, Fisher’s exact test). To confirm infectivity and measure the transmitted ZIKV dose, we enumerated infectious ZIKV in Ae. aegypti saliva using Vero cell plaque assays. The expectorated plaque forming units PFU varied by viral strain: MA66-infected expectorated 13±4 PFU (mean±SE, N = 13) compared to 29±6 PFU for PR15-infected (N = 13) and 35±8 PFU for BR15-infected (N = 6; ANOVA, df = 2, F = 3.8, P = 0.035). These laboratory vector competence results support an emerging consensus that Cx. tarsalis and Cx. quinquefasciatus are not vectors of ZIKV. These results also indicate that Ae. aegypti from California are efficient laboratory vectors of ancestral and contemporary Asian lineage ZIKV.

Since Zika virus (ZIKV) emerged as a global human health threat, numerous studies have pointed to Aedes aegypti as the primary vector due to its high competence and propensity to feed on humans. The majority of vector competence studies have been conducted between 26-28˚C, but arboviral extrinsic incubation periods (EIPs), and therefore transmission efficiency, are known to be affected strongly by temperature. To better understand the relationship between ZIKV EIPs and temperature, we evaluated the effect of adult mosquito exposure temperature on ZIKV infection, dissemination, and transmission in Ae. aegypti at four temperatures: 18˚C, 21˚C, 26˚C, and 30˚C. Mosquitoes were exposed to viremic mice infected with a 2015 Puerto Rican ZIKV strain, and engorged mosquitoes were sorted into the four temperatures with 80% RH and constant access to 10% sucrose. ZIKV infection, dissemination, and transmission rates were assessed via RT-qPCR from individual mosquito bodies, legs and wings, and saliva, respectively, at three to five time points per temperature from three to 31 days, based on expectations from other flavivirus EIPs. The median time from ZIKV ingestion to transmission (median EIP, EIP 50) at each temperature was estimated by fitting a generalized linear mixed model for each temperature. EIP 50 ranged from 5.1 days at 30˚C to 24.2 days at 21˚C. At 26˚C, EIP 50 was 9.6 days. At 18˚C, only 15% transmitted by day 31 so EIP 50 could not be estimated. This is among the first studies to characterize the effects of temperature on ZIKV EIP in Ae. aegypti, and the first to do so based on feeding of mosquitoes on a live, viremic host. This information is critical for modeling ZIKV transmission dynamics to understand geographic and seasonal limits of ZIKV risk; it is especially relevant for determining risk in subtropical regions with established Ae. aegypti populations and relatively high rates of return travel from the tropics (e.g. California or Florida), as these regions typically experience cooler temperature ranges than tropical regions.

Ross River virus (RRV), an alphavirus of the Togaviridae family, is the most medically significant mosquito-borne virus of Australia. Past RRV phylogenetic and evolutionary analyses have been based on partial genome analyses only. Three geographically distinct RRV lineages, the Eastern, the Western, and the supposedly extinct North-Eastern lineage, were classified previously. We sought to expand on past phylogenies through robust genome-scale phylogeny to better understand RRV genetic diversity and evolutionary dynamics. We analyzed 106 RRV complete coding sequences, which included 13 genomes available on NCBI and 94 novel sequences derived for this study, sampled throughout Western Australia (1977–2014) and during the substantial Pacific Islands RRV epidemic (1979–1980). Our final data set comprised isolates sampled over 59 years (1959–2018) from a range of locations. Four distinct genotypes were defined, with the newly described genotype 4 (G4) found to be the contemporary lineage circulating in Western Australia. The prior geographical classification of RRV lineages was not supported by our findings, with evidence of geographical and temporal cocirculation of distinct genetic groups. Bayesian Markov chain Monte Carlo (MCMC) analysis revealed that RRV lineages diverged from a common ancestor approximately 94 years ago, with distinct lineages emerging roughly every 10 years over the past 50 years in periodic bursts of genetic diversity. Our study has enabled a more robust analysis of RRV evolutionary history and resolved greater genetic diversity that had been previously defined by partial E2 gene analysis. IMPORTANCE Ross River virus (RRV) causes the most common mosquito-borne infection in Australia and causes a significant burden of suffering to infected individuals as well as being a large burden to the Australian economy. The genetic diversity of RRV and its evolutionary history have so far only been studied using partial E2 gene analysis with a limited number of isolates. Robust whole-genome analysis has not yet been conducted. This study generated 94 novel near-whole-genome sequences to investigate the evolutionary history of RRV to better understand its genetic diversity through comprehensive whole-genome phylogeny. A better understanding of RRV genetic diversity will enable better diagnostics, surveillance, and potential future vaccine design.

The “Asian tiger mosquito”, Aedes albopictus, is highly invasive, an aggressive biter and a major arbovirus vector. It is not currently present on mainland Australia despite being intercepted on numerous occasions at international ports and infesting the Torres Strait of Australia since at least 2004. In the current paper, we describe the invasion and current status of Ae. albopictus in the Torres Strait, as well as research conducted to assess the threat of this species becoming established in arbovirus transmission cycles on the Australian mainland. Genetic analysis of the invading population demonstrated that the Indonesian region was the likely origin of the invasion and not Papua New Guinea (PNG) as initially suspected. There was also intermixing between Torres Strait, PNG and Indonesian populations, indicating that the species could be re-introduced into the Torres Strait compromising any successful eradication programme. Vector competence experiments with endemic and exotic viruses revealed that Ae. albopictus from the Torres Strait are efficient alphavirus vectors, but less efficient flavivirus vectors. Ae.albopictus obtains blood meals from a range of vertebrate hosts (including humans), indicating that it could play a role in both zoonotic and human-mosquito arbovirus transmission cycles in Australia. Predictive models coupled with climate tolerance experiments suggest that a Torres Strait strain of Ae. albopictus could colonise southern Australia by overwintering in the egg stage before proliferating in the warmer months. Cohabitation experiments demonstrated that the presence of Aedes notoscriptus larvae in containers would not prevent the establishment of Ae. albopictus. Evidence from these studies, coupled with global experience suggests that we need to be prepared for the imminent invasion of Australia by Ae. albopictus by thoroughly understanding its biology and being willing to embrace emerging control technologies.

In 2005, established populations of Aedes albopictus (Skuse) were discovered in the Torres Strait, the region that separates Papua New Guinea from northern Australia. This increased the potential for this species to be introduced to mainland Australia. Because it is an arbovirus vector elsewhere, we undertook laboratory-based infection and transmission experiments to determine the potential for Ae. albopictus from the Torres Strait to become infected with and transmit the four major Australian endemic arboviruses--Murray Valley encephalitis virus, West Nile virus Kunjin strain (WNV(KUN)), Ross River virus (RRV), and Barmah Forest virus--as well as the exotic Japanese encephalitis virus. Ae. albopictus is susceptible to infection with all viruses, with infection rates ranging between 8% for WNV(KUN) and 71% for RRV. Transmission rates of approximately 25% were observed for RRV and Barmah Forest virus, but these were < 17% for Murray Valley encephalitis virus, WNV(KUN), and Japanese encephalitis virus. Given its relative vector competence for alphaviruses, we also examined the replication kinetics and extrinsic incubation periods required for transmission of RRV and chikungunya virus. Despite lower body titers, more mosquitoes reared and maintained at 28 degrees C became infected with and transmitted the virus than those reared and maintained at 22 degrees C. The minimum time between Ae. albopictus consuming an infected bloodmeal and transmitting chikungunya virus was 2 d at 28 degrees C and 4 d at 22 degrees C, and for RRV, it was 4 d, irrespective of the temperature. Given its opportunistic feeding habits and aggressive biting behavior, the establishment of Ae. albopictus on the Australian mainland could have a considerable impact on alphavirus transmission.

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.

The discovery of hepaciviruses in non-human hosts has accelerated following the advancement of high-throughput sequencing technology. Hepaciviruses have now been described in reptiles, fish, birds, and an extensive array of mammals. Using metagenomic sequencing on pooled samples of field-collected Culex annulirostris mosquitoes, we discovered a divergent hepacivirus-like sequence, named Jogalong virus, from the Kimberley region in northern Western Australia. Using PCR, we screened the same 300 individual mosquitoes and found just a single positive sample (1/300, 0.33%). Phylogenetic analysis of the hepacivirus NS5B protein places Jogalong virus within the genus Hepacivirus but on a distinct and deeply rooted monophyletic branch shared with duck hepacivirus, suggesting a notably different evolutionary history. Vertebrate barcoding PCR targeting two mitochondrial genes, cytochrome c oxidase subunit I and cytochrome b, indicated that the Jogalong virus-positive mosquito had recently fed on the tawny frogmouth (Podargus strigoides), although it is currently unknown whether this bird species contributes to the natural ecology of this virus.