Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia

Molecular evidence for multiple origins and high genetic differentiation of non-native winter crane fly, Trichocera maculipennis (Diptera: Trichoceridae), in the maritime Antarctic

Native biodiversity and ecosystems of Antarctica safeguarded from biological invasion face recent threats from non-native species, accelerated by increasing human activities and climate changes. Over two decades ago, the winter crane fly, Trichocera maculipennis, was first detected on King George Island. It has now successfully colonized several research stations across King George Island. To understand the origin, genetic diversity, and population structure of this Holarctic species, we conducted mitochondrial DNA cytochrome c oxidase subunit I (COI) sequence analysis across both its native and invasive ranges. In parallel, we performed microsatellite loci analysis within the invasive ranges, utilizing 12 polymorphic microsatellite markers. Furthermore, we compared body sizes among adult males and females collected from three different locations of King George Island. Our COI sequence analysis exhibited two different lineages present on King George Island. Lineage I was linked to Arctic Svalbard and Polish cave populations and Lineage II was related to Canadian Terra Nova National Park populations, implying multiple origins. Microsatellite analysis further exhibited high levels of genetic diversity and significant levels of genetic differentiation among invasive populations. Body sizes of adult T. maculipennis were significantly different among invasive populations but were not attributed to genetics. This significant genetic diversity likely facilitated the rapid colonization and establishment of T. maculipennis on King George Island, contributing to their successful invasion. Molecular analysis results revealed a substantial amount of genetic variation within invasive populations, which can serve as management units for invasive species control. Furthermore, the genetic markers we developed in the study will be invaluable tools for tracking impending invasion events and the travel routes of new individuals. Taken together, these findings illustrate the highly invasive and adaptable characteristics of T. maculipennis. Therefore, immediate action is necessary to mitigate their ongoing invasion and facilitate their eradication.

Complex multiple introductions drive fall armyworm invasions into Asia and Australia

The fall armyworm (FAW) Spodoptera frugiperda is thought to have undergone a rapid 'west-to-east' spread since 2016 when it was first identified in western Africa. Between 2018 and 2020, it was recorded from South Asia (SA), Southeast Asia (SEA), East Asia (EA), and Pacific/Australia (PA). Population genomic analyses enabled the understanding of pathways, population sources, and gene flow in this notorious agricultural pest species. Using neutral single nucleotide polymorphic (SNP) DNA markers, we detected genome introgression that suggested most populations in this study were overwhelmingly C- and R-strain hybrids (n = 252/262). SNP and mitochondrial DNA markers identified multiple introductions that were most parsimoniously explained by anthropogenic-assisted spread, i.e., associated with international trade of live/fresh plants and plant products, and involved 'bridgehead populations' in countries to enable successful pest establishment in neighbouring countries. Distinct population genomic signatures between Myanmar and China do not support the 'African origin spread' nor the 'Myanmar source population to China' hypotheses. Significant genetic differentiation between populations from different Australian states supported multiple pathways involving distinct SEA populations. Our study identified Asia as a biosecurity hotspot and a FAW genetic melting pot, and demonstrated the use of genome analysis to disentangle preventable human-assisted pest introductions from unpreventable natural pest spread.

Resistance Bioassays and Allele Characterization Inform Analysis of Spodoptera frugiperda (Lepidoptera: Noctuidae) Introduction Pathways in Asia and Australia

The fall armyworm (FAW) Spodoptera frugiperda (Smith; Lepidoptera: Noctuidae) is present in over 70 countries in Africa, Asia, and Oceania. Its rapid dispersal since 2016 when it was first reported in western Africa, and associated devastation to agricultural productivity, highlight the challenges posed by this pest. Currently, its management largely relies on insecticide sprays and transgenic Bacillus thuringiensis toxins, therefore understanding their responses to these agents and characteristics of any resistance genes enables adaptive strategies. In Australia, S. frugiperda was reported at the end of January 2020 in northern Queensland and by March 2020, also in northern Western Australia. As an urgent first response we undertook bioassays on two Australian populations, one each from these initial points of establishment. To assist with preliminary sensitivity assessment, two endemic noctuid pest species, Helicoverpa armigera (Hübner; Lepidoptera, Noctuidae) and Spodoptera litura (Fabricius; Lepidoptera, Noctuidae), were concurrently screened to obtain larval LC50 estimates against various insecticides. We characterized known resistance alleles from the VGSC, ACE-1, RyR, and ABCC2 genes to compare with published allele frequencies and bioassay responses from native and invasive S. frugiperda populations. An approximately 10× LC50 difference for indoxacarb was detected between Australian populations, which was approximately 28× higher than that reported from an Indian population. Characterization of ACE-1 and VGSC alleles provided further evidence of multiple introductions in Asia, and multiple pathways involving genetically distinct individuals in Australia. The preliminary bioassay results and resistance allele patterns from invasive S. frugiperda populations suggest multiple introductions have contributed to the pest's spread and challenge the axiom of its rapid 'west-to-east' spread.

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.

Dissecting the Species-Specific Virome in Culicoides of Thrace

Biting midges (Culicoides) are vectors of arboviruses of both veterinary and medical importance. The surge of emerging and reemerging vector-borne diseases and their expansion in geographical areas affected by climate change has increased the importance of understanding their capacity to contribute to novel and emerging infectious diseases. The study of Culicoides virome is the first step in the assessment of this potential. In this study, we analyzed the RNA virome of 10 Culicoides species within the geographical area of Thrace in the southeastern part of Europe, a crossing point between Asia and Europe and important path of various arboviruses, utilizing the Ion Torrent next-generation sequencing (NGS) platform and a custom bioinformatics pipeline based on TRINITY assembler and alignment algorithms. The analysis of the RNA virome of 10 Culicoides species resulted in the identification of the genomic signatures of 14 novel RNA viruses, including three fully assembled viruses and four segmented viruses with at least one segment fully assembled, most of which were significantly divergent from previously identified related viruses from the Solemoviridae, Phasmaviridae, Phenuiviridae, Reoviridae, Chuviridae, Partitiviridae, Orthomyxoviridae, Rhabdoviridae, and Flaviviridae families. Each Culicoides species carried a species-specific set of viruses, some of which are related to viruses from other insect vectors in the same area, contributing to the idea of a virus-carrier web within the ecosystem. The identified viruses not only expand our current knowledge on the virome of Culicoides but also set the basis of the genetic diversity of such viruses in the area of southeastern Europe. Furthermore, our study highlights that such metagenomic approaches should include as many species as possible of the local virus-carrier web that interact and share the virome of a geographical area.

Spatial and temporal distribution of Culicoides species in the New England region of New South Wales, Australia between 1990 and 2018

Culicoides are one of the smallest hematophagous flies measuring 1-5 mm in size with only females seeking blood for egg development. The present study investigated spatio-temporal distribution of Culicoides species trapped between 1990 and 2018 at 13 sites in the New England region of NSW, Australia using automated light traps. Trapping locations were divided into three subregions (tablelands, slopes and plains). Nineteen Culicoides species were identified. Culicoides marksi and C. austropalpalis were the most abundant and widespread species. Culicoides brevitarsis, the principal vector of livestock diseases in New South Wales comprised 2.9% of the total catch and was detected in 12 of the 13 locations in the study. Abundance as determined by Log10 Culicoides count per trapping event for the eight most abundant species did not vary significantly with season but trended towards higher counts in summer for C. marksi (P = 0.09) and C. austropalpalis (P = 0.05). Significant geographic variation in abundance was observed for C. marksi, C. austropalpalis and C. dycei with counts decreasing with increasing altitude from the plains to the slopes and tablelands. Culicoides victoriae exhibited the reverse trend in abundance (P = 0.08). Greater abundance during the warmer seasons and at lower altitudes for C. marksi and C. austropalpalis was indicative of temperature and rainfall dependence in this region with moderate summer dominance in rainfall. The Shannon-Wiener diversity index of species was higher on the tablelands (H = 1.59) than the slopes (H = 1.33) and plains (H = 1.08) with evenness indices of 0.62, 0.46 and 0.39 respectively. Culicoides species on the tablelands were more diverse than on the slopes and plains where C. marksi and C. austropalpalis dominated. The temporal and spatial variation in abundance, diversity and evenness of species reported in this diverse region of Australia provides additional insight into Culicoides as pests and disease vectors and may contribute to future modelling studies.

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex

Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Genomic analysis of bluetongue virus episystems in Australia and Indonesia

The distribution of bluetongue viruses (BTV) in Australia is represented by two distinct and interconnected epidemiological systems (episystems)—one distributed primarily in the north and one in the east. The northern episystem is characterised by substantially greater antigenic diversity than the eastern episystem; yet the forces that act to limit the diversity present in the east remain unclear. Previous work has indicated that the northern episystem is linked to that of island South East Asia and Melanesia, and that BTV present in Indonesia, Papua New Guinea and East Timor, may act as source populations for new serotypes and genotypes of BTV to enter Australia’s north. In this study, the genomes of 49 bluetongue viruses from the eastern episystem and 13 from Indonesia were sequenced and analysed along with 27 previously published genome sequences from the northern Australian episystem. The results of this analysis confirm that the Australian BTV population has its origins in the South East Asian/Melanesian episystem, and that incursions into northern Australia occur with some regularity. In addition, the presence of limited genetic diversity in the eastern episystem relative to that found in the north supports the presence of substantial, but not complete, barriers to gene flow between the northern and eastern Australian episystems. Genetic bottlenecks between each successive episystem are evident, and appear to be responsible for the reduction in BTV genetic diversity observed in the north to south–east direction.

Sellers' Revisited: A Big Data Reassessment of Historical Outbreaks of Bluetongue and African Horse Sickness due to the Long-Distance Wind Dispersion of Culicoides Midges

The possibility that outbreaks of bluetongue (BT) and African horse sickness (AHS) might occur via long-distance wind dispersion (LDWD) of their insect vector (Culicoides spp.) was proposed by R. F. Sellers in a series of papers published between 1977 and 1991. These investigated the role of LDWD by means of visual examination of the wind direction of synoptic weather charts. Based on the hypothesis that simple wind direction analysis, which does not allow for wind speed, might have led to spurious conclusions, we reanalyzed six of the outbreak scenarios described in Sellers' papers. For this reanalysis, we used a custom-built Big Data application ("TAPPAS") which couples a user-friendly web-interface with an established atmospheric dispersal model ("HYSPLIT"), thus enabling more sophisticated modeling than was possible when Sellers undertook his analyzes. For the two AHS outbreaks, there was strong support from our reanalysis of the role of LDWD for that in Spain (1966), and to a lesser degree, for the outbreak in Cyprus (1960). However, for the BT outbreaks, the reassessments were more complex, and for one of these (western Turkey, 1977) we could discount LDWD as the means of direct introduction of the virus. By contrast, while the outbreak in Cyprus (1977) showed LDWD was a possible means of introduction, there is an apparent inconsistency in that the outbreaks were localized while the dispersion events covered much of the island. For Portugal (1956), LDWD from Morocco on the dates suggested by Sellers is very unlikely to have been the pathway for introduction, and for the detection of serotype 2 in Florida (1982), LDWD from Cuba would require an assumption of a lengthy survival time of the midges in the air column. Except for western Turkey, the BT reanalyses show the limitation of LDWD modeling when used by itself, and indicates the need to integrate susceptible host population distribution (and other covariate) data into the modeling process. A further refinement, which will become increasingly important to assess LDWD, will be the use of virus and vector genome sequence data collected from potential source and the incursion sites.

Mitochondrial DNA and trade data support multiple origins of Helicoverpa armigera (Lepidoptera, Noctuidae) in Brazil

The Old World bollworm Helicoverpa armigera is now established in Brazil but efforts to identify incursion origin(s) and pathway(s) have met with limited success due to the patchiness of available data. Using international agricultural/horticultural commodity trade data and mitochondrial DNA (mtDNA) cytochrome oxidase I (COI) and cytochrome b (Cyt b) gene markers, we inferred the origins and incursion pathways into Brazil. We detected 20 mtDNA haplotypes from six Brazilian states, eight of which were new to our 97 global COI-Cyt b haplotype database. Direct sequence matches indicated five Brazilian haplotypes had Asian, African, and European origins. We identified 45 parsimoniously informative sites and multiple substitutions per site within the concatenated (945 bp) nucleotide dataset, implying that probabilistic phylogenetic analysis methods are needed. High diversity and signatures of uniquely shared haplotypes with diverse localities combined with the trade data suggested multiple incursions and introduction origins in Brazil. Increasing agricultural/horticultural trade activities between the Old and New Worlds represents a significant biosecurity risk factor. Identifying pest origins will enable resistance profiling that reflects countries of origin to be included when developing a resistance management strategy, while identifying incursion pathways will improve biosecurity protocols and risk analysis at biosecurity hotspots including national ports.

Genotyping of whole genome amplified reduced representation libraries reveals a cryptic population of Culicoides brevitarsis in the Northern Territory, Australia

BACKGROUND

The advent of genotyping by Next Generation Sequencing has enabled rapid discovery of thousands of single nucleotide polymorphism (SNP) markers and high throughput genotyping of large populations at an affordable cost. Genotyping by sequencing (GBS), a reduced representation library sequencing method, allows highly multiplexed sequencing of genomic subsets. This method has limitations for small organisms with low amounts of genomic DNA, such as the bluetongue virus (BTV) vectors, Culicoides midges.

RESULTS

This study employed the GBS method to isolate SNP markers de novo from whole genome amplified Culicoides brevitarsis genomic DNA. The individuals were collected from regions representing two different Australian patterns of BTV strain distribution: the Northern Territory (NT) and the east coast. We isolated 8145 SNPs using GBS. Phylogenetic analysis conducted using the filtered 3263 SNPs revealed the presence of a distinct C. brevitarsis sub-population in the NT and this was confirmed by analysis of mitochondrial DNA. Two loci showed a very strong signal for selection and were unique to the NT population. Bayesian analysis with STRUCTURE indicated a possible two-population cluster.

CONCLUSIONS

The results suggest that genotyping vectors with high density markers in combination with biological and environmental data is useful. However, more extensive sampling over a wider spatial and temporal range is needed. The presence of sub-structure in populations and loci under natural selection indicates the need for further investigation of the role of vectors in shaping the two Australian systems of BTV transmission. The described workflow is transferable to genotyping of small, non-model organisms, including arthropod vectors of pathogens of economic and medical importance.