Modelling the Wolbachia incompatible insect technique: strategies for effective mosquito population elimination.

BACKGROUND: The Wolbachia incompatible insect technique (IIT) shows promise as a method for eliminating populations of invasive mosquitoes such as Aedes aegypti (Linnaeus) (Diptera: Culicidae) and reducing the incidence of vector-borne diseases such as dengue, chikungunya and Zika. Successful implementation of this biological control strategy relies on high-fidelity separation of male from female insects in mass production systems for inundative release into landscapes. Processes for sex-separating mosquitoes are typically error-prone and laborious, and IIT programmes run the risk of releasing Wolbachia-infected females and replacing wild mosquito populations. RESULTS: We introduce a simple Markov population process model for studying mosquito populations subjected to a Wolbachia-IIT programme which exhibit an unstable equilibrium threshold. The model is used to study, in silico, scenarios that are likely to yield a successful elimination result. Our results suggest that elimination is best achieved by releasing males at rates that adapt to the ever-decreasing wild population, thus reducing the risk of releasing Wolbachia-infected females while reducing costs. CONCLUSIONS: While very high-fidelity sex separation is required to avoid establishment, release programmes tend to be robust to the release of a small number of Wolbachia-infected females. These findings will inform and enhance the next generation of Wolbachia-IIT population control strategies that are already showing great promise in field trials.

Wolbachia-induced apoptosis associated with increased fecundity in Laodelphax striatellus (Hemiptera: Delphacidae).

Wolbachia influence the fitness of their invertebrate hosts. They have effects on reproductive incompatibility and egg production. Although the former are well characterized, the mechanistic basis of the latter is unclear. Here, we investigate whether apoptosis, which has been implicated in fecundity in model insects, influences the interaction between fecundity and Wolbachia in the planthopper Laodelphax striatellus. Wolbachia-infected females produced about 30% more eggs than uninfected females. We used the terminal deoxyribonucleotidyl transferase (TDT)-mediated dUTP-digoxigenin nick end labeling staining to visualize apoptosis. Microscopic observations indicated that the Wolbachia strain wStri increased the number of ovarioles that contained apoptotic nurse cells in both young and aged adult females. The frequency of apoptosis was much higher in the infected females. The increased fecundity appeared to be a result of apoptosis of nurse cells, which provide nutrients to the growing oocytes. In addition, cell apoptosis inhibition by caspase messenger RNA interference in Wolbachia-infected L. striatellus markedly decreased egg numbers. Together, these data suggest that wStri might enhance fecundity by increasing the number of apoptotic cells in the ovaries in a caspase-dependent manner. Our findings establish a link between Wolbachia-induced apoptosis and egg production effects mediated by Wolbachia, although the way in which the endosymbiont influences caspase levels remains to be determined.

Evolutionary divergence of mitochondrial genomes in two Tetranychus species distributed across different climates.

There is increasing evidence that mitochondrial genomes (mitogenomes) can be under selection, whereas the selective regimes shaping mitogenome evolution remain largely unclear. To test for mitogenome evolution in relation to the climate adaptation, we explored mtDNA variation in two spider mite (Tetranychus) species that distribute across different climates. We sequenced 26 complete mitogenomes of Tetranychus truncates, which occurs in both warm and cold regions, and nine complete mitogenomes of Tetranychus pueraricola, which is restricted to warm regions. Patterns of evolution in the two species' mitogenomes were compared through a series of dN /dS methods and physicochemical profiles of amino acid replacements. We found that: (1) the mitogenomes of both species were under widespread purifying selection; (2) elevated directional adaptive selection was observed in the T. truncatus mitogenome, perhaps linked to the cold climates adaptation of T. truncatus; and (3) the strength of selection varied across genes, and diversifying positive selection detected on ND4 and ATP6 pointed to their crucial roles during adaptation to different climatic conditions. This study gained insight into the mitogenome evolution in relation to the climate adaptation.

Does increased heat resistance result in higher susceptibility to predation? A test using Drosophila melanogaster selection and hardening.

Heat resistance of ectotherms can be increased both by plasticity and evolution, but these effects may have trade-offs resulting from biotic interactions. Here, we test for predation costs in Drosophila melanogaster populations with altered heat resistance produced by adult hardening and directional selection for increased heat resistance. In addition, we also tested for genetic trade-offs by testing heat resistance in lines that have evolved under increased predation risk. We show that while 35/37 °C hardening increases heat resistance as expected, it does not increase predation risk from jumping spiders or mantids; in fact, there was an indication that survival may have increased under predation following a triple 37 °C compared to a single 35 °C hardening treatment. Flies that survived a 39 °C selection cycle showed lower survival under predation, suggesting a predation cost of exposure to a more severe heat stress. There was, however, no correlated response to selection because survival did not differ between control and selected lines after selection was relaxed for one or two generations. In addition, lines selected for increased predation risk did not differ in heat resistance. Our findings suggest independent evolutionary responses to predation and heat as measured in laboratory assays, and no costs of heat hardening on susceptibility to predation.

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics, genomics and marine LiDAR data.

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long-term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.

Spatial and Temporal Variation in Aedes aegypti and Aedes albopictus (Diptera: Culicidae) Numbers in the Yogyakarta Area of Java, Indonesia, With Implications for Wolbachia Releases.

of mosquito vector populations, particularly through Wolbachia endosymbionts. The success of these strategies depends on understanding the dynamics of vector populations. In preparation for Wolbachia releases around Yogyakarta, we have studied Aedes populations in five hamlets. Adult monitoring with BioGent- Sentinel (BG-S) traps indicated that hamlet populations had different dynamics across the year; while there was an increase in Aedes aegypti (L.) and Aedes albopictus (Skuse) numbers in the wet season, species abundance remained relatively stable in some hamlets but changed markedly (>2 fold) in others. Local rainfall a month prior to monitoring partly predicted numbers of Ae. aegypti but not Ae. albopictus. Site differences in population size indicated by BG-S traps were also evident in ovitrap data. Egg or larval collections with ovitraps repeated at the same location suggested spatial autocorrelation (<250 m) in the areas of the hamlets where Ae. aegypti numbers were high. Overall, there was a weak negative association (r<0.43) between Ae. aegypti and Ae. albopictus numbers in ovitraps when averaged across collections. Ae. albopictus numbers in ovitraps and BG-S traps were positively correlated with vegetation around areas where traps were placed, while Ae. aegypti were negatively correlated with this feature. These data inform intervention strategies by defining periods when mosquito densities are high, highlighting the importance of local site characteristics on populations, and suggesting relatively weak interactions between Ae. aegypti and Ae. albopictus. They also indicate local areas within hamlets where consistently high mosquito densities may influence Wolbachia invasions and other interventions.

Mitochondrial DNA variants help monitor the dynamics of Wolbachia invasion into host populations.

Wolbachia is the most widespread endosymbiotic bacterium of insects and other arthropods that can rapidly invade host populations. Deliberate releases of Wolbachia into natural populations of the dengue fever mosquito, Aedes aegypti, are used as a novel biocontrol strategy for dengue suppression. Invasion of Wolbachia through the host population relies on factors such as high fidelity of the endosymbiont transmission and limited immigration of uninfected individuals, but these factors can be difficult to measure. One way of acquiring relevant information is to consider mitochondrial DNA (mtDNA) variation alongside Wolbachia in field-caught mosquitoes. Here we used diagnostic mtDNA markers to differentiate infection-associated mtDNA haplotypes from those of the uninfected mosquitoes at release sites. Unique haplotypes associated with Wolbachia were found at locations outside Australia. We also performed mathematical and qualitative analyses including modelling the expected dynamics of the Wolbachia and mtDNA variants during and after a release. Our analyses identified key features in haplotype frequency patterns to infer the presence of imperfect maternal transmission of Wolbachia, presence of immigration and possibly incomplete cytoplasmic incompatibility. We demonstrate that ongoing screening of the mtDNA variants should provide information on maternal leakage and immigration, particularly in releases outside Australia. As we demonstrate in a case study, our models to track the Wolbachia dynamics can be successfully applied to temporal studies in natural populations or Wolbachia release programs, as long as there is co-occurring mtDNA variation that differentiates infected and uninfected populations.

Invasion of Wolbachia at the residential block level is associated with local abundance of Stegomyia aegypti, yellow fever mosquito, populations and property attributes.

Wolbachia can suppress dengue and control mosquito populations and this depends on the successful invasion of Wolbachia-infected mosquitoes into local populations. Ovitrap data collected during the recent invasion of wMel-infected Stegomyia aegypti (Diptera: Culicidae) (Linnaeus) into Gordonvale near Cairns, Australia, were used to identify variables that help predict the success of localized invasion. Based on the variance in Wolbachia frequencies across Gordonvale as well as at another release site at Yorkeys Knob in comparison to simulations, it was estimated that on average 2-4 females contributed eggs to an ovitrap. By collating ovitrap data from two collection periods at the start of the release from residential blocks, it was found that uninfected mosquitoes had a patchy distribution across the release site. Residential blocks with relatively high uninfected mosquito numbers were less easily invaded by Wolbachia than blocks with low numbers. The numbers of uninfected mosquitoes in ovitraps were negatively correlated with the proportion of brick houses in a residential block, whereas local Wolbachia frequencies were correlated positively with this variable as well as negatively with the amount of shading in a yard and availability of breeding sites. These findings point to proxy measures for predicting the ease of localized invasion of Wolbachia.

Evidence for adaptive divergence of thermal responses among Bemisia tabaci populations from tropical Colombia following a recent invasion.

There is an increasing evidence that populations of ectotherms can diverge genetically in response to different climatic conditions, both within their native range and (in the case of invasive species) in their new range. Here, we test for such divergence in invasive whitefly Bemisia tabaci populations in tropical Colombia, by considering heritable variation within and between populations in survival and fecundity under temperature stress, and by comparing population differences with patterns established from putatively neutral microsatellite markers. We detected significant differences among populations linked to mean temperature (for survival) and temperature variation (for fecundity) in local environments. A QST  - FST analysis indicated that phenotypic divergence was often larger than neutral expectations (QST  > FST ). Particularly, for survival after a sublethal heat shock, this divergence remained linked to the local mean temperature after controlling for neutral divergence. These findings point to rapid adaptation in invasive whitefly likely to contribute to its success as a pest species. Ongoing evolutionary divergence also provides challenges in predicting the likely impact of Bemisia in invaded regions.

Balancing genetic uniqueness and genetic variation in determining conservation and translocation strategies: a comprehensive case study of threatened dwarf galaxias, Galaxiella pusilla (Mack) (Pisces: Galaxiidae).

Genetic markers are widely used to define and manage populations of threatened species based on the notion that populations with unique lineages of mtDNA and well-differentiated nuclear marker frequencies should be treated separately. However, a danger of this approach is that genetic uniqueness might be emphasized at the cost of genetic diversity, which is essential for adaptation and is potentially boosted by mixing geographically separate populations. Here, we re-explore the issue of defining management units, focussing on a detailed study of Galaxiella pusilla, a small freshwater fish of national conservation significance in Australia. Using a combination of microsatellite and mitochondrial markers, 51 populations across the species range were surveyed for genetic structure and diversity. We found an inverse relationship between genetic differentiation and genetic diversity, highlighting a long-term risk of deliberate isolation of G. pusilla populations based on protection of unique lineages. Instead, we adopt a method for identifying genetic management units that takes into consideration both uniqueness and genetic variation. This produced a management framework to guide future translocation and re-introduction efforts for G. pusilla, which contrasted to the framework based on a more traditional approach that may overlook important genetic variation in populations.

Effect of Wolbachia on insecticide susceptibility in lines of Aedes aegypti.

Two stable infections of Wolbachia pipientis, wMelPop and wMel, now established in Aedes aegypti, are being used in a biocontrol program to suppress the transmission of dengue. Any effects of Wolbachia infection on insecticide resistance of mosquitoes may undermine the success of this program. Bioassays of Ae. aegypti were conducted to test for differences in response to insecticides between Wolbachia infected (wMelPop, wMel) and uninfected lines. Insecticides screened were bifenthrin, the pyrethroid commonly used for adult knockdown, as well as larvicides: Bacillus thuringiensis var. israelensis, the organophosphate, temephos and the insect growth regulator, s-methoprene. While differences in response between lines were detected for some insecticides, no obvious or consistent effects related to presence of Wolbachia infection were observed. Spreading Wolbachia infections are, therefore, unlikely to affect the efficacy of traditional chemical control of mosquito outbreaks.

A proline repeat polymorphism of the Frost gene of Drosophila melanogaster showing clinal variation but not associated with cold resistance.

Genetic polymorphisms underlying adaptive shifts in thermal responses are poorly known even though studies are providing a detailed understanding of these responses at the cellular and physiological levels. The Frost gene of Drosophila melanogaster is a prime candidate for thermal adaptation; it is up-regulated under cold stress and knockdown of this gene influences cold resistance. Here we describe an amino-acid INDEL polymorphism in proline repeat number in the structural component of this gene. The two main repeats, accounting for more than 90% of alleles in eastern Australia, show a strong clinal pattern; the 6P allele was at a high frequency in tropical locations, and the 10P allele was common in temperate populations. However, the frequency of these alleles was not associated with three different assays of cold resistance. Adult transcription level of Frost was also unrelated to cold resistance as measured through post chill coma mobility. The functional significance of the proline repeat polymorphism therefore remains unclear despite its clinal pattern. The data also demonstrate the feasibility of using Roche/454 sequencing for establishing clinal patterns.

Population genetics of the wheat curl mite (Aceria tosichella Keifer) in Australia: implications for the management of wheat pathogens.

The wheat curl mite (WCM), Aceria tosichella Keifer, is a polyphagous eriophyoid mite and the primary vector of wheat streak mosaic virus (WSMV) and five other viral pathogens in cereals. Previous research using molecular markers and a series of laboratory experiments found A. tosichella in Australia to consist of two genetically distinct lineages, which have broad overlapping distributions and differ in their ability to transmit WSMV under controlled conditions. This pattern of transmission also appears to be apparent in the field, whereby a strong association between WSMV detection and a single WCM lineage has been detected. In this study, we conduct a population genetic analysis and provide information on the genetic structure of the Australian viruliferous WCM lineage. We assessed genetic differentiation of 16 WCM populations using nine microsatellite markers. Strong evidence for extensive gene flow and low genetic structuring throughout the Australian wheatbelt was evident, with an exception for Western Australian and far north Queensland populations that appear to be genetically isolated. The data also indicate genetic patterns consistent with an arrhenotokous parthenogenetic mode of reproduction. Implications of these findings are discussed with reference to the management of WCM and associated cereal pathogens in Australia and overseas.

Changes in the genetic structure of Aedes aegypti (Diptera: Culicidae) populations in Queensland, Australia, across two seasons: implications for potential mosquito releases.

Diseases transmitted by mosquitoes could be controlled if vector populations were replaced with strains that have reduced vector competency. Such a strategy is being developed for control of dengue virus which is transmitted by Aedes aegypti (L.) (Diptera: Culicidae). Mosquitoes artificially infected with the bacterium, Wolbachia pipientis Hertig, are being assessed as candidates for release at the adult stage with the aim of replacement of the wild population. Wolbachia can reduce the capacity of Ae. aegypti to transmit dengue virus and has potential to be driven through the natural population via a system of cytoplasmic incompatibility. Deployment of benign mosquito strains will be influenced by population size and structure of wild-type Ae. aegypti in proposed release areas, as well as rates of gene flow among populations in the wet and dry tropical seasons. Mosquitoes from northern Queensland were screened with genetic markers to find an optimal locality for release of a benign strain of Ae. aegypti. The inland towns of Chillagoe and Charters Towers and the coastal town of Ingham had mosquito populations that were partly genetically isolated from mosquitoes in other areas across both seasons. These locations may be suitable release sites if it is important for the released strain to be restricted during initial phases of implementation. Smaller genetic differences were also evident among other regions and were consistent over two seasons (wet and dry).

The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations.

Dengue fever is the most important mosquito-borne viral disease of humans with more than 50 million cases estimated annually in more than 100 countries. Disturbingly, the geographic range of dengue is currently expanding and the severity of outbreaks is increasing. Control options for dengue are very limited and currently focus on reducing population abundance of the major mosquito vector, Aedes aegypti. These strategies are failing to reduce dengue incidence in tropical communities and there is an urgent need for effective alternatives. It has been proposed that endosymbiotic bacterial Wolbachia infections of insects might be used in novel strategies for dengue control. For example, the wMelPop-CLA Wolbachia strain reduces the lifespan of adult A. aegypti mosquitoes in stably transinfected lines. This life-shortening phenotype was predicted to reduce the potential for dengue transmission. The recent discovery that several Wolbachia infections, including wMelPop-CLA, can also directly influence the susceptibility of insects to infection with a range of insect and human pathogens has markedly changed the potential for Wolbachia infections to control human diseases. Here we describe the successful transinfection of A. aegypti with the avirulent wMel strain of Wolbachia, which induces the reproductive phenotype cytoplasmic incompatibility with minimal apparent fitness costs and high maternal transmission, providing optimal phenotypic effects for invasion. Under semi-field conditions, the wMel strain increased from an initial starting frequency of 0.65 to near fixation within a few generations, invading A. aegypti populations at an accelerated rate relative to trials with the wMelPop-CLA strain. We also show that wMel and wMelPop-CLA strains block transmission of dengue serotype 2 (DENV-2) in A. aegypti, forming the basis of a practical approach to dengue suppression.

Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission.

Genetic manipulations of insect populations for pest control have been advocated for some time, but there are few cases where manipulated individuals have been released in the field and no cases where they have successfully invaded target populations. Population transformation using the intracellular bacterium Wolbachia is particularly attractive because this maternally-inherited agent provides a powerful mechanism to invade natural populations through cytoplasmic incompatibility. When Wolbachia are introduced into mosquitoes, they interfere with pathogen transmission and influence key life history traits such as lifespan. Here we describe how the wMel Wolbachia infection, introduced into the dengue vector Aedes aegypti from Drosophila melanogaster, successfully invaded two natural A. aegypti populations in Australia, reaching near-fixation in a few months following releases of wMel-infected A. aegypti adults. Models with plausible parameter values indicate that Wolbachia-infected mosquitoes suffered relatively small fitness costs, leading to an unstable equilibrium frequency <30% that must be exceeded for invasion. These findings demonstrate that Wolbachia-based strategies can be deployed as a practical approach to dengue suppression with potential for area-wide implementation.

A combination of molecular and morphological approaches resolves species in the taxonomically difficult genus Procladius Skuse (Diptera: Chironomidae) despite high intra-specific morphological variation.

Molecular approaches for identifying aquatic macroinvertebrate species are increasingly being used but there is ongoing debate about the number of DNA markers needed to differentiate species accurately. Here, we use two mitochondrial genes (cytochrome oxidase I, cytochrome b) and a nuclear gene (carbamoylphosphate synthetase) to differentiate species variation within the taxonomically challenging chironomid genus Procladius from southern Australia, a genus which is important for pollution monitoring. The mitochondrial genes indicated cryptic species that were subsequently linked to morphological variation at the larval and pupal stage. Two species previously described based on morphological criteria were linked to molecular markers, and there was evidence for additional cryptic species. Each genetic marker provided different information, highlighting the importance of considering multiple genes when dissecting taxonomically difficult groups, particularly those used in pollution monitoring.

Dynamics of the "popcorn" Wolbachia infection in outbred Aedes aegypti informs prospects for mosquito vector control.

Forty percent of the world's population is at risk of contracting dengue virus, which produces dengue fever with a potentially fatal hemorrhagic form. The wMelPop Wolbachia infection of Drosophila melanogaster reduces life span and interferes with viral transmission when introduced into the mosquito Aedes aegypti, the primary vector of dengue virus. Wolbachia has been proposed as an agent for preventing transmission of dengue virus. Population invasion by Wolbachia depends on levels of cytoplasmic incompatibility, fitness effects, and maternal transmission. Here we characterized these traits in an outbred genetic background of a potential target population of Ae. aegypti using two crossing schemes. Cytoplasmic incompatibility was strong in this background, and the maternal transmission rate of Wolbachia was high. The infection substantially reduced longevity of infected adult females, regardless of whether adults came from larvae cultured under high or low levels of nutrition or density. The infection reduced the viability of diapausing and nondiapausing eggs. Viability was particularly low when eggs were laid by older females and when diapausing eggs had been stored for a few weeks. The infection affected mosquito larval development time and adult body size under different larval nutrition levels and densities. The results were used to assess the potential for wMelPop-CLA to invade natural populations of Ae. aegypti and to develop recommendations for the maintenance of fitness in infected mosquitoes that need to compete against field insects.

Genetic mapping of adaptive wing size variation in Drosophila simulans.

Many ecologically important traits exhibit latitudinal variation. Body size clines have been described repeatedly in insects across multiple continents, suggesting that similar selective forces are shaping these geographical gradients. It is unknown whether these parallel clinal patterns are controlled by the same or different genetic mechanism(s). We present here, quantitative trait loci (QTL) analysis of wing size variation in Drosophila simulans. Our results show that much of the wing size variation is controlled by a QTL on Chr 3L with relatively minor contribution from other chromosome arms. Comparative analysis of the genomic positions of the QTL indicates that the major QTL on Chr 3 are distinct in D. simulans and D. melanogaster, whereas the QTL on Chr 2R might overlap between species. Our results suggest that parallel evolution of wing size clines could be driven by non-identical genetic mechanisms but in both cases involve a major QTL as well as smaller effects of other genomic regions.

A comprehensive assessment of geographic variation in heat tolerance and hardening capacity in populations of Drosophila melanogaster from eastern Australia.

We examined latitudinal variation in adult and larval heat tolerance in Drosophila melanogaster from eastern Australia. Adults were assessed using static and ramping assays. Basal and hardened static heat knockdown time showed significant linear clines; heat tolerance increased towards the tropics, particularly for hardened flies, suggesting that tropical populations have a greater hardening response. A similar pattern was evident for ramping heat knockdown time at 0.06°C min(-1) increase. There was no cline for ramping heat knockdown temperature (CT(max) ) at 0.1°C min(-1) increase. Acute (static) heat knockdown temperature increased towards temperate latitudes, probably reflecting a greater capacity of temperate flies to withstand sudden temperature increases during summer in temperate Australia. Larval viability showed a quadratic association with latitude under heat stress. Thus, patterns of heat resistance depend on assay methods. Genetic correlations in thermotolerance across life stages and evolutionary potential for critical thermal limits should be the focus of future studies.