Search

Abstract Background Wolbachia symbiosis in Aedes aegypti is an emerging biocontrol measure against dengue. However, assessing its real-world efficacy is challenging due to the non-randomised, field-based nature of most intervention studies. This research re-evaluates the spatial–temporal impact of Wolbachia interventions on dengue incidence using a large battery of quasi-experimental methods and assesses each method’s validity. Methods A systematic search for Wolbachia intervention data was conducted via PubMed. Efficacy was reassessed using commonly-used quasi-experimental approaches coupled with extensive robustness checks, including geospatial placebo tests and a simulation study. Intervention efficacies across multiple study sites were computed using high-resolution aggregations to examine heterogeneities across sites and study periods. We further designed a stochastic simulation framework to assess the methods’ ability to reproduce true intervention efficacies (IE). Results Wolbachia interventions in Singapore, Malaysia, and Brazil were associated to significant decreases in dengue incidence, with reductions ranging from 48.17% to 69.19%. IEs varied with location and duration. Malaysia showed increasing efficacy over time, while Brazil exhibited initial success with subsequent decline, hinting at operational challenges. Singapore's strategy was highly effective despite partial saturation. Simulations identified Synthetic Control Methods (SCM) and its variant, count Synthetic Control Method (cSCM), as superior in precision, with the smallest percentage errors in efficacy estimation. These methods also demonstrated robustness in placebo tests. Conclusions Wolbachia interventions exhibit consistent protective effects against dengue. SCM and cSCM provided the most precise and robust estimates of IEs, validated across simulated and real-world settings.

Abstract Wolbachia is a genus of maternally inherited endosymbionts that can affect reproduction of their hosts and influence metabolic processes. The pollinator, Valisia javana, is common in the male syconium of the dioecious fig Ficus hirta. Based on a high-quality chromosome-level V. javana genome with PacBio long-read and Illumina short-read sequencing, we discovered a sizeable proportion of Wolbachia sequences and used these to assemble two novel Wolbachia strains belonging to supergroup A. We explored its phylogenetic relationship with described Wolbachia strains based on MLST sequences and the possibility of induction of CI (cytoplasmic incompatibility) in this strain by examining the presence of cif genes known to be responsible for CI in other insects. We also identified mobile genetic elements including prophages and insertion sequences, genes related to biotin synthesis and metabolism. A total of two prophages and 256 insertion sequences were found. The prophage WOjav1 is cryptic (structure incomplete) and WOjav2 is relatively intact. IS5 is the dominant transposon family. At least three pairs of type I cif genes with three copies were found which may cause strong CI although this needs experimental verification; we also considered possible nutritional effects of the Wolbachia by identifying genes related to biotin production, absorption and metabolism. This study provides a resource for further studies of Wolbachia-pollinator-host plant interactions.

Abstract Background Wolbachia wMel is the most commonly used strain in rear and release strategies for Aedes aegypti mosquitoes that aim to inhibit the transmission of arboviruses such as dengue, Zika, Chikungunya and yellow fever. However, the long-term establishment of wMel in natural Ae. aegypti populations raises concerns that interactions between Wolbachia wMel and Ae. aegypti may lead to changes in the host genome, which could affect useful attributes of Wolbachia that allow it to invade and suppress disease transmission. Results We applied an evolve-and-resequence approach to study genome-wide genetic changes in Ae. aegypti from the Cairns region, Australia, where Wolbachia wMel was first introduced more than 10 years ago. Mosquito samples were collected at three different time points in Gordonvale, Australia, covering the phase before (2010) and after (2013 and 2018) Wolbachia releases. An additional three locations where Wolbachia replacement happened at different times across the last decade were also sampled in 2018. We found that the genomes of mosquito populations mostly remained stable after Wolbachia release, with population differences tending to reflect the geographic location of the populations rather than Wolbachia infection status. However, outlier analysis suggests that Wolbachia may have had an influence on some genes related to immune response, development, recognition and behavior. Conclusions Ae. aegypti populations remained geographically distinct after Wolbachia wMel releases in North Australia despite their Wolbachia infection status. At some specific genomic loci, we found signs of selection associated with Wolbachia, suggesting potential evolutionary impacts can happen in the future and further monitoring is warranted.

Abstract Background Dengue suppression often relies on control of the mosquito vector, Aedes aegypti, through applications of insecticides of which the pyrethroid group has played a dominant role. Insecticide resistance is prevalent in Ae. aegypti around the world, and the resulting reduction of insecticide efficacy is likely to exacerbate the impact of dengue. Dengue has been a public health problem in Saudi Arabia, particularly in Jeddah, since its discovery there in the 1990s, and insecticide use for vector control is widespread throughout the city. An alternative approach to insecticide use, based on blocking dengue transmission in mosquitoes by the endosymbiont Wolbachia, is being trialed in Jeddah following the success of this approach in Australia and Malaysia. Knowledge of insecticide resistance status of mosquito populations in Jeddah is a prerequisite for establishing a Wolbachia-based dengue control program as releases of Wolbachia mosquitoes succeed when resistance status of the release population is similar to that of the wild population. Methods WHO resistance bioassays of mosquitoes with deltamethrin, permethrin and DDT were used in conjunction with TaqMan® SNP Genotyping Assays to characterize mutation profiles of Ae. aegypti. Results Screening of the voltage-sensitive sodium channel (Vssc), the pyrethroid target site, revealed mutations at codons 989, 1016 and 1534 in Ae. aegypti from two districts of Jeddah. The triple mutant homozygote (1016G/1534C/989P) was confirmed from Al Safa and Al Rawabi. Bioassays with pyrethroids (Type I and II) and DDT showed that mosquitoes were resistant to each of these compounds based on WHO definitions. An association between Vssc mutations and resistance was established for the Type II pyrethroid, deltamethrin, with one genotype (989P/1016G/1534F) conferring a survival advantage over two others (989S/1016V/1534C and the triple heterozygote). An indication of synergism of Type I pyrethroid activity with piperonyl butoxide suggests that detoxification by cytochrome P450s accounts for some of the pyrethroid resistance response in Ae. aegypti populations from Jeddah. Conclusions The results provide a baseline for monitoring and management of resistance as well as knowledge of Vssc genotype frequencies required in Wolbachia release populations to ensure homogeneity with the target field population. Vssc mutation haplotypes observed show some similarity with those from Ae. aegypti in southeast Asia and the Indo-Pacific, but the presence of the triple mutant haplotype in three genotypes indicates that the species in this region may have a unique population history. Graphical Abstract

Abstract Background The two-spotted spider mite, Tetranychus urticae, is a major agricultural pest with a cosmopolitan distribution, and its polyphagous habits provide a model for investigating herbivore-plant interactions. There are two body color forms of T. urticae with a different host preference. Comparative genomics and transcriptomics are used here to investigate differences in responses of the forms to host plants at the molecular level. Biological responses of the two forms sourced from multiple populations are also presented. Results We carried out principal component analysis of transcription changes in three red and three green T. urticae populations feeding on their original host (common bean), and three hosts to which they were transferred: cotton, cucumber and eggplant. There were differences among the forms in gene expression regardless of their host plant. In addition, different changes in gene expression were evident in the two forms when responding to the same host transfer. We further compared biological performance among populations of the two forms after feeding on each of the four hosts. Fecundity of 2-day-old adult females showed a consistent difference between the forms after feeding on bean. We produced a 90.1-Mb genome of the red form of T. urticae with scaffold N50 of 12.78 Mb. Transcriptional profiles of genes associated with saliva, digestion and detoxification showed form-dependent responses to the same host and these genes also showed host-specific expression effects. Conclusions Our research revealed that forms of T. urticae differ in host-determined transcription responses and that there is form-dependent plasticity in the transcriptomic responses. These differences may facilitate the extreme polyphagy shown by spider mites, although fitness differences on hosts are also influenced by population differences unrelated to color form.

Abstract Background In the inner city of Yogyakarta, Indonesia, insecticide resistance is expected in the main dengue vector, Aedes aegypti, because of the intensive local application of pyrethroid insecticides. However, detailed information about the nature of resistance in this species is required to assist the release of Wolbachia mosquitoes in a dengue control program, so that we can ensure that insecticide resistance in the strain of Ae. aegypti being released matches that of the background population. Methods High-resolution melt genotyping was used to screen for kdr mutations associated with pyrethroid resistance in the voltage-sensitive sodium channel (V SSC ) gene in Ae. aegypti of some areas in the inner city of Yogyakarta. Results The results show that the V1016G mutation predominated, with individuals homozygous for the 1016G allele at a frequency of 82.1% and the mutant allele G at a frequency of 92%. Two patterns of co-occurrence of mutations were detected in this study, homozygous individuals V1016G/S989P; and heterozygous individuals V1016G/F1534C/S989P. We found the simultaneous occurrence of kdr mutations V1016G and F1534C at all collection sites, but not within individual mosquitoes. Homozygous mutants at locus 1016 were homozygous wild-type at locus 1534 and vice versa, and heterozygous V1016G were also heterozygous for F1534C. The most common tri-locus genotype co-occurrences were homozygous mutant 1016GG and homozygous wild-type FF1534, combined with homozygous mutant 989PP (GG/FF/PP) at a frequency of 38.28%. Conclusions Given the relatively small differences in frequency of resistance alleles across the city area, locality variations in resistance should have minor implications for the success of Wolbachia mosquito trials being undertaken in the Yogyakarta area.

Abstract Background Host-associated microbial communities play an important role in the fitness of insect hosts. However, the factors shaping microbial communities in wild populations, including genetic background, ecological factors, and interactions among microbial species, remain largely unknown. Results Here, we surveyed microbial communities of the small brown planthopper (SBPH, Laodelphax striatellus) across 17 geographical populations in China and Japan by using 16S rRNA amplicon sequencing. Using structural equation models (SEM) and Mantel analyses, we show that variation in microbial community structure is likely associated with longitude, annual mean precipitation (Bio12), and mitochondrial DNA variation. However, a Wolbachia infection, which is spreading to northern populations of SBPH, seems to have a relatively greater role than abiotic factors in shaping microbial community structure, leading to sharp decreases in bacterial taxon diversity and abundance in host-associated microbial communities. Comparative RNA-Seq analyses between Wolbachia-infected and -uninfected strains indicate that the Wolbachia do not seem to alter the immune reaction of SBPH, although Wolbachia affected expression of metabolism genes. Conclusion Together, our results identify potential factors and interactions among different microbial species in the microbial communities of SBPH, which can have effects on insect physiology, ecology, and evolution. Video Abstract

Abstract Background Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species. Results We find disproportionately high rates of gene gains in internal branches in the species’ phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species’ differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one. Conclusion Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change.

Abstract Background Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China. Results Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage. Conclusions Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

Abstract Background With Aedes aegypti mosquitoes now being released in field programmes aimed at disease suppression, there is interest in identifying factors influencing the mating and invasion success of released mosquitoes. One factor that can increase release success is size: released males may benefit competitively from being larger than their field counterparts. However, there could be a risk in releasing only large males if small field females avoid these males and instead prefer small males. Here we investigate this risk by evaluating mating success for mosquitoes differing in size. Results We measured mating success indirectly by coupling size with Wolbachia-infected or uninfected mosquitoes and scoring cytoplasmic incompatibility. Large females showed no evidence of a mating preference, whereas small males were relatively more successful than large males when mating with small females, exhibiting an advantage of around 20–25%. Conclusions Because field females typically encompass a wide range of sizes while laboratory reared (and released) males typically fall into a narrow size range of large mosquitoes, these patterns can influence the success of release programmes which rely on cytoplasmic incompatibility to suppress populations and initiate replacement invasions. Releases could include some small males generated under low food or crowded conditions to counter this issue, although this would need to be weighed against issues associated with costs of producing males of various size classes.

Abstract Background Populations of herbivorous insects may become genetically differentiated because of local adaptation to different hosts and climates as well as historical processes, and further genetic divergence may occur following the development of reproductive isolation among populations. Here we investigate the population genetic structure of the orchard pest peach fruit moth (PFM) Carposina sasakii (Lepidoptera: Carposinidae) in China, which shows distinct biological differences when characterized from different host plants. Genetic diversity and genetic structure were assessed among populations from seven plant hosts and nine regions using 19 microsatellite loci and a mitochondrial sequence. Results Strong genetic differentiation was found among geographical populations representing distinct geographical regions, but not in host-associated populations collected from the same area. Mantel tests based on microsatellite loci indicated an association between genetic differentiation and geographical distance, and to a lesser extent environmental differentiation. Approximate Bayesian Computation analyses supported the scenario that PFM likely originated from a southern area and dispersed northwards before the last glacial maximum during the Quaternary. Conclusions Our analyses suggested a strong impact of geographical barriers and historical events rather than host plants on the genetic structure of the PFM; however, uncharacterized environmental factors and host plants may also play a role. Studies on adaptive shifts in this moth should take into account geographical and historical factors.

Abstract Background Distinguishing orthologous and paralogous relationships between genes across multiple species is essential for comparative genomic analyses. Various computational approaches have been developed to resolve these evolutionary relationships, but strong trade-offs between precision and recall of orthologue prediction remains an ongoing challenge. Results Here we present Orthonome, an orthologue prediction pipeline, designed to reduce the trade-off between orthologue capture rates (recall) and accuracy of multi-species orthologue prediction. The pipeline compares sequence domains and then forms sequence-similar clusters before using phylogenetic comparisons to identify inparalogues. It then corrects sequence similarity metrics for fragment and gene length bias using a novel scoring metric capturing relationships between full length as well as fragmented genes. The remaining genes are then brought together for the identification of orthologues within a phylogenetic framework. The orthologue predictions are further calibrated along with inparalogues and gene births, using synteny, to identify novel orthologous relationships. We use 12 high quality Drosophila genomes to show that, compared to other orthologue prediction pipelines, Orthonome provides orthogroups with minimal error but high recall. Furthermore, Orthonome is resilient to suboptimal assembly/annotation quality, with the inclusion of draft genomes from eight additional Drosophila species still providing >6500 1:1 orthologues across all twenty species while retaining a better combination of accuracy and recall than other pipelines. Orthonome is implemented as a searchable database and query tool along with multiple-sequence alignment browsers for all sets of orthologues. The underlying documentation and database are accessible at http://www.orthonome.com . Conclusion We demonstrate that Orthonome provides a superior combination of orthologue capture rates and accuracy on complete and draft drosophilid genomes when tested alongside previously published pipelines. The study also highlights a greater degree of evolutionary conservation across drosophilid species than earlier thought.

Abstract Background Dengue fever, the most prevalent global arboviral disease, represents an important public health problem in Indonesia. Control of dengue relies on the control of its main vector, the mosquito Aedes aegypti, yet nothing is known about the population history and genetic structure of this insect in Indonesia. Our aim was to assess the spatio-temporal population genetic structure of Ae. aegypti in Yogyakarta, a densely populated region on Java with common dengue outbreaks. Methods We used multiple marker systems (microsatellites, nuclear and mitochondrial genome-wide single nucleotide polymorphisms generated via Restriction-site Associated DNA sequencing) to analyze 979 Ae. aegypti individuals collected from the Yogyakarta city and the surrounding hamlets during the wet season in 2011 and the following dry season in 2012. We employed individual- and group-based approaches for inferring genetic structure. Results We found that Ae. aegypti in Yogyakarta has spatially structured and seasonally stable populations. The spatial structuring was significant for the nuclear and mitochondrial markers, while the temporal structuring was non-significant. Nuclear markers identified three main genetic clusters, showing that hamlets have greater genetic isolation from each other and from the inner city sites. However, one hamlet experienced unrestricted mosquito interbreeding with the inner city, forming a single genetic cluster. Genetic distance was poorly correlated with the spatial distance among mosquito samples, suggesting stronger influence of human-assisted gene flow than active mosquito movement on spatial genetic structure. A star-shaped mitochondrial haplotype network and a significant R 2 test statistic (R 2 = 0.0187, P = 0.001) support the hypothesis that Ae. aegypti in Yogyakarta originated from a small or homogeneous source and has undergone a relatively recent demographic expansion. Conclusion We report the first insights into the spatio-temporal genetic structure and the underlying processes in the dengue fever mosquito from Yogyakarta, Indonesia. Our results provide valuable information on the effectiveness of local control measures as well as guidelines for the implementation of novel biocontrol strategies such as release of Wolbachia-infected mosquitoes.

Abstract Background Introduced Wolbachia bacteria can influence the susceptibility of Aedes aegypti mosquitoes to arboviral infections as well as having detrimental effects on host fitness. Previous field trials demonstrated that the wMel strain of Wolbachia effectively and durably invades Ae. aegypti populations. Here we report on trials of a second strain, wMelPop-PGYP Wolbachia, in field sites in northern Australia (Machans Beach and Babinda) and central Vietnam (Tri Nguyen, Hon Mieu Island), each with contrasting natural Ae. aegypti densities. Methods Mosquitoes were released at the adult or pupal stages for different lengths of time at the sites depending on changes in Wolbachia frequency as assessed through PCR assays of material collected through Biogents-Sentinel (BG-S) traps and ovitraps. Adult numbers were also monitored through BG-S traps. Changes in Wolbachia frequency were compared across hamlets or house blocks. Results Releases of adult wMelPop-Ae. aegypti resulted in the transient invasion of wMelPop in all three field sites. Invasion at the Australian sites was heterogeneous, reflecting a slower rate of invasion in locations where background mosquito numbers were high. In contrast, invasion across Tri Nguyen was relatively uniform. After cessation of releases, the frequency of wMelPop declined in all sites, most rapidly in Babinda and Tri Nguyen. Within Machans Beach the rate of decrease varied among areas, and wMelPop was detected for several months in an area with a relatively low mosquito density. Conclusions These findings highlight challenges associated with releasing Wolbachia-Ae. aegypti combinations with low fitness, albeit strong virus interference properties, as a means of sustainable control of dengue virus transmission.

Abstract Background Recent releases have been carried out with Aedes aegypti mosquitoes infected with the w MelPop mosquito cell-line adapted (w MelPop-CLA) strain of Wolbachia. This infection introduced from Drosophila provides strong blockage of dengue and other arboviruses but also has large fitness costs in laboratory tests. The releases were used to evaluate the fitness of released infected mosquitoes, and (following termination of releases) to test for any effects of w MelPop-CLA on wing size and shape when mosquitoes were reared under field conditions. Methods We monitored gravid females via double sticky traps to assess the reproductive success of w MelPop-CLA-infected females and also sampled the overall mosquito population post-release using Biogent Sentinel traps. Morphometric analyses were used to evaluate infection effects on wing shape as well as size. Results Oviposition success as assessed through double sticky traps was unrelated to size of released mosquitoes. However, released mosquitoes with lower wing loading were more successful. Furthermore, w MelPop-CLA-infected mosquitoes had 38.3% of the oviposition success of uninfected mosquitoes based on the predicted infection frequency after release. Environmental conditions affected wing shape and particularly size across time in uninfected mosquitoes, but not in naturally-reared w MelPop-CLA-infected mosquitoes. Although the overall size and shape do not differ between naturally-reared w MelPop-CLA-infected and uninfected mosquitoes, the infected mosquitoes tended to have smaller wings than uninfected mosquitoes during the cooler November in comparison to December. Conclusion These results confirm the lower fitness of w MelPop-CLA infection under field conditions, helping to explain challenges associated with a successful invasion by this strain. In the long run, invasion may depend on releasing strains carrying insecticide resistance or egg desiccation resistance, combined with an active pre-release population suppression program.