A decade of stability for wMel Wolbachia in natural Aedes aegypti populations

Mosquitoes carrying Wolbachia endosymbionts are being released in many countries for arbovirus control. The wMel strain of Wolbachia blocks Aedes-borne virus transmission and can spread throughout mosquito populations by inducing cytoplasmic incompatibility. Aedes aegypti mosquitoes carrying wMel were first released into the field in Cairns, Australia, over a decade ago, and with wider releases have resulted in the near elimination of local dengue transmission. The long-term stability of Wolbachia effects is critical for ongoing disease suppression, requiring tracking of phenotypic and genomic changes in Wolbachia infections following releases. We used a combination of field surveys, phenotypic assessments, and Wolbachia genome sequencing to show that wMel has remained stable in its effects for up to a decade in Australian Ae. aegypti populations. Phenotypic comparisons of wMel-infected and uninfected mosquitoes from near-field and long-term laboratory populations suggest limited changes in the effects of wMel on mosquito fitness. Treating mosquitoes with antibiotics used to cure the wMel infection had limited effects on fitness in the next generation, supporting the use of tetracycline for generating uninfected mosquitoes without off-target effects. wMel has a temporally stable within-host density and continues to induce complete cytoplasmic incompatibility. A comparison of wMel genomes from pre-release (2010) and nine years post-release (2020) populations show few genomic differences and little divergence between release locations, consistent with the lack of phenotypic changes. These results indicate that releases of Wolbachia-infected mosquitoes for population replacement are likely to be effective for many years, but ongoing monitoring remains important to track potential evolutionary changes.     

The endosymbiont Wolbachia in Eurasian populations of Drosophila melanogaster

The endosymbiotic [alpha]-proteobacteria Wolbachia is widely spread among arthropods and Filariidae nematodes. This bacterium is transmitted vertically via a transovarian route. Wolbachia is a cause of several reproductive abnormalities in the host species. We analyzed the isofemale lines created using flies collected from Drosophila melanogaster natural populations for infection with the endosymbiont Wolbachia. Wolbachia were genotyped according to five variable markers: the presence of insertion sequence IS5 in two loci, the copy number of two minisatellite repeats, and an inversion. Overall, 665 isofemale lines isolated from the populations of D. melanogaster from Ukraine, Belarus, Moldova, Caucasus, Central Asia, Ural, Udmurtia, Altai, West and East Siberia, and Far East in 1974 through 2005 were used in the work. The samples from Ukrainian, Altaian, and Middle Asian populations were largest. The infection rate of D. melanogaster populations from Middle Asia, Altaian, and Eastern Europe (Ukraine, Moldavia, and Belarus) with Wolbachia amounted to 64, 56, and 39%, respectively. The D. melanogaster population from the Caucasus displayed heterogeneity in the genotypes of this cytoplasmic infection. The Wolbachia genotype wMel, detected in all the populations studied, was the most abundant. The genotype wMelCS2 was always present in the populations from Middle Asia and Altai and was among the rare variants in the D. melanogaster populations from the Eastern Europe. Single instances of the Wolbachia genotype wMelCS occurred in a few flies from the Central Asian and Altai populations, but was not found this genotype in the other regions.     

A wMel Wolbachia variant in Aedes aegypti from field-collected Drosophila melanogaster with increased phenotypic stability under heat stress

Mosquito-borne diseases remain a major cause of morbidity and mortality. Population replacement strategies involving the wMel strain of Wolbachia are being used widely to control mosquito-borne diseases. However, these strategies may be influenced by temperature because wMel is vulnerable to heat. wMel infections in Drosophila melanogaster are genetically diverse, but few transinfections of wMel variants have been generated in Aedes aegypti. Here, we successfully transferred a wMel variant (termed wMelM) originating from a field-collected D. melanogaster into Ae. aegypti. The new wMelM variant (clade I) is genetically distinct from the original wMel transinfection (clade III), and there are no genomic differences between wMelM in its original and transinfected host. We compared wMelM with wMel in its effects on host fitness, temperature tolerance, Wolbachia density, vector competence, cytoplasmic incompatibility and maternal transmission under heat stress in a controlled background. wMelM showed a higher heat tolerance than wMel, likely due to higher overall densities within the mosquito. Both wMel variants had minimal host fitness costs, complete cytoplasmic incompatibility and maternal transmission, and dengue virus blocking under laboratory conditions. Our results highlight phenotypic differences between Wolbachia variants and wMelM shows potential as an alternative strain in areas with strong seasonal temperature fluctuations.     

Electrophoretic variability in natural populations of Drosophila melanogaster and Drosophila simulans

Genetic variation at 59 gene loci coding for enzymes (50) and larval proteins (9) has been studied in sympatric populations of Drosophila melanogaster and D. simulans from insular and continental origin. The average number of alleles per locus, the mean proportion of polymorphic loci and the mean heterozygosity are similar both within and between species. There are however some significant differences between D. simulans populations in the genotypic frequencies for four polymorphic loci.     

Inversion polymorphism in Indian natural populations of Drosophila melanogaster

Six natural populations (three urban and three rural) of Drosophila melanogaster from India were analysed for chromosome inversions, revealing the presence of 19 different paracentric autosomal inversions. One new inversion has also been detected in a laboratory stock established from flies collected from Kerala. In total 20 different paracentric inversions in Indian D. melanogaster have been detected during the present study, and of these, 4 are common cosmopolitans; 2 are rare cosmopolitans; 7 are recurrent endemics; and 7 are unique endemics. The quantitative data clearly show that the urban populations are different from the rural ones with respect to inversion polymorphism.     

From parasite to mutualist: rapid evolution of Wolbachia in natural populations of Drosophila

Wolbachia are maternally inherited bacteria that commonly spread through host populations by causing cytoplasmic incompatibility, often expressed as reduced egg hatch when uninfected females mate with infected males. Infected females are frequently less fecund as a consequence of Wolbachia infection. However, theory predicts that because of maternal transmission, these "parasites" will tend to evolve towards a more mutualistic association with their hosts. Drosophila simulans in California provided the classic case of a Wolbachia infection spreading in nature. Cytoplasmic incompatibility allowed the infection to spread through individual populations within a few years and from southern to northern California (more than 700 km) within a decade, despite reducing the fecundity of infected females by 15%-20% under laboratory conditions. Here we show that the Wolbachia in California D. simulans have changed over the last 20 y so that infected females now exhibit an average 10% fecundity advantage over uninfected females in the laboratory. Our data suggest smaller but qualitatively similar changes in relative fecundity in nature and demonstrate that fecundity-increasing Wolbachia variants are currently polymorphic in natural populations.     

Extensive paternal mtDNA leakage in natural populations of Drosophila melanogaster

Strict maternal inheritance is considered a hallmark of animal mtDNA. Although recent reports suggest that paternal leakage occurs in a broad range of species, it is still considered an exceptionally rare event. To evaluate the impact of paternal leakage on the evolution of mtDNA, it is essential to reliably estimate the frequency of paternal leakage in natural populations. Using allele‐specific real‐time quantitative PCR (RT‐qPCR), we show that heteroplasmy is common in natural populations with at least 14% of the individuals carrying multiple mitochondrial haplotypes. However, the average frequency of the minor mtDNA haplotype is low (0.8%), which suggests that this pervasive heteroplasmy has not been noticed before due to a lack of power in sequencing surveys. Based on the distribution of mtDNA haplotypes in the offspring of heteroplasmic mothers, we found no evidence for strong selection against one of the haplotypes. We estimated that the rate of paternal leakage is 6% and that at least 100 generations are required for complete sorting of mtDNA haplotypes. Despite the high proportion of heteroplasmic individuals in natural populations, we found no evidence for recombination between mtDNA molecules, suggesting that either recombination is rare or recombinant haplotypes are counter‐selected. Our results indicate that evolutionary studies using mtDNA as a marker might be biased by paternal leakage in this species.     

Contrasting patterns of natural variation in global Drosophila melanogaster populations

Despite the popularity of Drosophila melanogaster in functional and evolutionary genetics, the global pattern of natural variation has not yet been comprehensively described in this species. For the first time, we report a combined survey using neutral microsatellites and mitochondrial sequence variation jointly. Thirty-five populations originating from five continents were compared. In agreement with previous microsatellite studies, sub-Saharan African populations were the most variable ones. Consistent with previous reports of a single 'out of Africa' habitat expansion, we found that non-African populations contained a subset of the African alleles. The pattern of variation detected for the mitochondrial sequences differed substantially. The most divergent haplotypes were detected in the Mediterranean region while Africa harboured most haplotypes, which were all closely related. In the light of the well-established African origin of D. melanogaster, our results cast severe doubts about the suitability of mtDNA for biogeographic inference in this model organism.     

Polymorphism of yellow locus in Drosophila melanogaster from natural populations

We studied molecular characteristics of yellow (y; 1-0.0) locus, which determines the body coloration of phenotypically wild-type and mutant alleles isolated from geographically distant populations of Drosophila melanogaster in different years. According to Southern data, restrictions map of yellow locus of all studied strains differ from each other as well as from that of Oregon stock. FISH analysis shows that in the neighborhood of yellow locus in X chromosome neither P nor hobo elements are found in y1-775 stock, while only hobo is found there in y1-859 and y1-866 stocks, only P element in y+sn849 stock, and both elements in y1-719 stock. Thus, all studied mutant variants of yellow are of independent origin. Yellow locus residing at the very end of X chromosome (region 1A5-8 of cytologic map) carries significantly more transposon than retrotransposon-induced mutations, as compared to white locus (regions 3C2). It is possible that transposons are more active than retrotransposons at the chromosomal ends of D. melanogaster.     

Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster

Restriction map variation in 64 X chromosome lines extracted from three different populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing restriction enzymes for a 106-kb region encompassing the yellow gene and the achaete-scute complex that is located in the region of reduced crossing over close to the telomere. Nine restriction site polymorphisms (out of 176 sites scored) and 19 length polymorphisms (15 insertions and 4 deletions) were detected. The estimated level of heterozygosity per nucleotide, H = 0.0003, is much lower than that reported for autosomal and sex-linked loci located in regions with normal levels of crossing over. The overall frequency of polymorphic restriction sites is reduced. Six out of nine restriction site polymorphisms are unique and the other three have frequencies less than 0.17. Some large insertions have reached relatively high frequencies, 0.08 to 0.17. Consistent with the theoretically predicted negative relationship between crossing over and the magnitude of linkage disequilibrium, an increase in the relative number of nonrandom associations was observed in the y-ac-sc region.     

Abdominal pigmentation in Drosophila melanogaster females from natural Indian populations

Females of Drosophila melanogaster collected from five geographically distant populations in India were analysed for the intensity of pigmentation in the 5th, 6th and 7th segments of the abdomen. In all three segments, this intensity was found to vary among individuals of any given population, and, furthermore, different populations differ with respect to this phenotypic trait. Statistical analysis revealed significant intra- and interpopulational variation. A clinical pattern was detected: females from populations closer to the equator tended to have lighter cuticle, in which case differences between the three segments could not be detected and all three segments responded both independently and jointly to latitudinal variation, as indicated by a statistically significant positive correlation between latitude and pigmentation score. This is the first report on abdominal pigmentation analysis in natural populations of D. melanogaster that provides evidence that phenotypic flexibility reflects temperature differences, as a result of which abdominal pigmentation shows geographic differentiation.     

Genetic heterogeneity in natural populations of Drosophila melanogaster for ability to withstand dessication

Summary Strains set up from single inseminated females of D. melanogaster derived from two wild populations have been shown to differ in their ability to withstand dessication, as measured by mortalities after 16 hours in a dry environment, thus there are genes segregating in wild populations for ability to withstand dessication. A more detailed study on strains from one of the wild populations, showed that strains with high wet and dry weights lose water by dessication relatively less rapidly and have lower mortalities, than strains with lower wet and dry weights.Variability within and between five inbred strains was studied with results as above. Heritabilities for wet weight, dry weight, and mortality were 0.40, 0.41 and 0.60 respectively, showing the likelihood that the traits would be amenable to further genetic analysis.The relevance of the results are discussed in relation to stress to high temperatures, and the ecology of the species in general.The competent technical assistance of Miss Clare Escott is gratefully acknowledged. Financial support from the Australian Institute of Nuclear Science and Engineering and the Australian Research Grants Committee, is gratefully acknowledged.     

Structure of variation in enzyme activity in natural Drosophila melanogaster populations

Enzyme activity variation was assessed in several isofemale lines originating from two Hungarian Drosophila melanogaster populations. Samples from each population were taken from from two villages; 8-9 isofemale lines were established from each village. The activities of ADH, alphaGPDH, IDH and 6PGDH were determined in the adults (in the F1 generation) and in the larvae (in the F3 generation) as well. Enzyme activities were measured on starch gel after electrophoresis. The activity of the enzyme was detected in a single individual and it was also possible to determine its genotype. The results showed that most of the variation occurred within sites for all four enzymes. This within site variation was more or less equally partitioned into within and between isofemale line (family) components. A smaller portion of variation was attributable to the differences between the populations. Nevertheless, adult alphaGPDH, and larval IDH and 6PGDH activities exhibited significant differences between the two populations. Variation in larval activities of all enzymes was higher than that of the adults, but 6PGDH had considerably higher variation in the adults. The greater variation in larval activities probably reflected the greater environmental variation in the microhabitat of the larvae compared to that of the adults. Larval activities of the investigated enzymes showed much stronger correlation than adult activities. The correlation pattern in the adults differed greatly between the two populations.     

Restriction-map variation at the zeste-tko region in natural populations of Drosophila melanogaster

Restriction-map variation in 64 X chromosome lines extracted from three different natural populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing enzymes for a 20-kb region including the zeste and tko genes. Ten restriction-site and four length polymorphisms (two insertions and two deletions) were detected. Contrary to the predicted lower level of variation for genes on the X chromosome, the level of variation attributable to nucleotide substitution (estimated heterozygosity/nucleotide = 0.004) was similar to that previously reported for autosomal loci. The amount of insertion/deletion variation in the studied region was within the range observed in autosomal regions and thus not explainable by a simple selection model against the effects of insertional mutations. A general lack of linkage disequilibrium between polymorphic sites was observed.