Microsatellite and Wolbachia analysis in Rhagoletis cerasi natural populations: population structuring and multiple infections

Rhagoletis cerasi (Diptera: Tephritidae) is a major pest of sweet and sour cherries in Europe and parts of Asia. Despite its economic significance, there is a lack of studies on the genetic structure of R. cerasi populations. Elucidating the genetic structure of insects of economic importance is crucial for developing phenological‐predictive models and environmental friendly control methods. All natural populations of R. cerasi have been found to harbor the endosymbiont Wolbachia pipientis, which widely affects multiple biological traits contributing to the evolution of its hosts, and has been suggested as a tool for the biological control of insect pests and disease vectors. In the current study, the analysis of 18 R. cerasi populations collected in Greece, Germany, and Russia using 13 microsatellite markers revealed structuring of R. cerasi natural populations, even at close geographic range. We also analyzed the Wolbachia infection status of these populations using 16S rRNA‐, MLST‐ and wsp‐based approaches. All 244 individuals screened were positive for Wolbachia. Our results suggest the fixation of the wCer1 strain in Greece while wCer2, wCer4, wCer5, and probably other uncharacterized strains were also detected in multiply infected individuals. The role of Wolbachia and its potential extended phenotypes needs a thorough investigation in R. cerasi. Our data suggest an involvement of this symbiont in the observed restriction in the gene flow in addition to a number of different ecological factors.     

Tracing horizontal Wolbachia movements among bees (Anthophila): a combined approach using multilocus sequence typing data and host phylogeny

The endosymbiotic bacterium Wolbachia enhances its spread via vertical transmission by generating reproductive effects in its hosts, most notably cytoplasmic incompatibility (CI). Additionally, frequent interspecific horizontal transfer is evident from a lack of phylogenetic congruence between Wolbachia and its hosts. The mechanisms of this lateral transfer are largely unclear. To identify potential pathways of Wolbachia movements, we performed multilocus sequence typing of Wolbachia strains from bees (Anthophila). Using a host phylogeny and ecological data, we tested various models of horizontal endosymbiont transmission. In general, Wolbachia strains seem to be randomly distributed among bee hosts. Kleptoparasite‐host associations among bees as well as other ecological links could not be supported as sole basis for the spread of Wolbachia. However, cophylogenetic analyses and divergence time estimations suggest that Wolbachia may persist within a host lineage over considerable timescales and that strictly vertical transmission and subsequent random loss of infections across lineages may have had a greater impact on Wolbachia strain distribution than previously estimated. Although general conclusions about Wolbachia movements among arthropod hosts cannot be made, we present a framework by which precise assumptions about shared evolutionary histories of Wolbachia and a host taxon can be modelled and tested.     

Diversity and evolution of the Wolbachia endosymbionts of Bemisia (Hemiptera: Aleyrodidae) whiteflies

Wolbachia is the most prevalent symbiont described in arthropods to date. Wolbachia can manipulate host reproduction, provide nutrition to insect hosts and protect insect hosts from pathogenic viruses. So far, 13 supergroups of Wolbachia have been identified. The whitefly Bemisia tabaci is a complex containing more than 28 morphologically indistinguishable cryptic species. Some cryptic species of this complex are invasive. In this study, we report a comprehensive survey of Wolbachia in B. tabaci and its relative B. afer from 1658 insects representing 54 populations across 13 provinces of China and one state of Australia. Based on the results of PCR or sequencing of the 16S rRNA gene, the overall rates of Wolbachia infection were 79.6% and 0.96% in the indigenous and invasive Bemisia whiteflies, respectively. We detected a new Wolbachia supergroup by sequencing five molecular marker genes including 16S rRNA, groEL, gltA, hcpA, and fbpA genes. Data showed that many protein‐coding genes have limitations in detecting and classifying newly identified Wolbachia supergroups and thus raise a challenge to the known Wolbachia MLST standard analysis system. Besides, the other Wolbachia strains detected from whiteflies were clustered into supergroup B. Phylogenetic trees of whitefly mitochondrial cytochrome oxidase subunit I and Wolbachia multiple sequencing typing genes were not congruent. In addition, Wolbachia was also detected outside the special bacteriocytes in two cryptic species by fluorescence in situ hybridization, indicating the horizontal transmission of Wolbachia. Our results indicate that members of Wolbachia are far from well explored.     

Wolbachia infections and superinfections in cytoplasmically incompatible populations of the European cherry fruit fly Rhagoletis cerasi (Diptera,Tephritidae)

Wolbachia is an obligately intracellular, maternally inherited bacterium which has been detected in many arthropods. Wolbachia infections disperse in host populations by mechanisms such as cytoplasmic incompatibility (CI). CI leads to embryonic mortality which occurs when infected males mate with uninfected females or females with a different Wolbachia strain. Populations of the European cherry fruit fly Rhagoletis cerasi (Diptera, Tephritidae) were found to be infected by two different Wolbachia strains, wCer1 and wCer2. Superinfections with both strains occurred throughout southern and central Europe and infections with wCer1 were found in northern, western and eastern Europe. Strong unidirectional CI between European populations of R. cerasi were first reported in the 1970s. From the conformity in the recent geographical distribution of the Wolbachia infections and the CI expression patterns found 25 years ago it was deduced that wCer2 potentially causes CI in R. cerasi. The comparison of the geographical distributions indicated that wCer1 + 2 must have spread into wCer1-infected populations in some areas. In other regions, a spread of wCer1 + 2 was probably prevented by dispersal barriers. There, a sharp transition from infected to superinfected populations suggested regional isolation between wCer1 and wCer1 + 2-infected populations.     

Epidemiology in evolutionary time: the case of Wolbachia horizontal transmission between arthropod host species

Wolbachia are bacterial endosymbionts that manipulate the reproduction of their arthropod hosts. Although theory suggests that infections are frequently lost within host species due to the evolution of resistance, Wolbachia infect a huge number of species worldwide. This apparent paradox suggests that horizontal transmission between host species has been a key factor in shaping the global Wolbachia pandemic. Because Wolbachia infections are thus acquired and lost like any other infection, we use a standard epidemiological model to analyse Wolbachia horizontal transmission dynamics over evolutionary time. Conceptually modifying the model, we apply it not to transmission between individuals but between species. Because, on evolutionary timescales, infections spread frequently between closely related species and occasionally over large phylogenetic distances, we represent the set of host species as a small‐world network that satisfies both requirements. Our model reproduces the effect of basic epidemiological parameters, which demonstrates the validity of our approach. We find that the ratio between transmission rate and recovery rate is crucial for determining the proportion of infected species (incidence) and that, in a given host network, the incidence may still be increasing over evolutionary time. Our results also point to the importance of occasional transmission over long phylogenetic distances for the observed high incidence levels of Wolbachia. In conclusion, we are able to explain why Wolbachia are so abundant among arthropods, although selection for resistance within hosts often leads to infection loss. Furthermore, our unorthodox approach of using epidemiology in evolutionary time can be applied to all symbionts that use horizontal transmission to infect new hosts.     

Invasion genetics of American cherry fruit fly in Europe and signals of hybridization with the European cherry fruit fly

The American cherry fruit fly is an invasive pest species in Europe, of serious concern in tart cherry production as well as for the potential to hybridize with the European cherry fruit fly, Rhagoletis cerasi L. (Diptera: Tephritidae), which might induce new pest dynamics. In the first European reports, the question arose whether only the eastern American cherry fruit fly, Rhagoletis cingulata (Loew) (Diptera: Tephritidae), is present, or also the closely related western American cherry fruit fly, Rhagoletis indifferens Curran. In this study, we investigate the species status of European populations by comparing these with populations of both American species from their native ranges, the invasion dynamics in German (first report in 1993) and Hungarian (first report in 2006) populations, and we test for signals of hybridization with the European cherry fruit fly. Although mtDNA sequence genealogy could not separate the two American species, cross‐species amplification of 14 microsatellite loci separated them with high probabilities (0.99–1.0) and provided evidence for R. cingulata in Europe. German and Hungarian R. cingulata populations differed significantly in microsatellite allele frequencies, mtDNA haplotype and wing pattern distributions, and both were genetically depauperate relative to North American populations. The diversity suggests independent founding events in Germany and Hungary. Within each country, R. cingulata displayed little or no structure in any trait, which agrees with rapid local range expansions. In cross‐species amplifications, signals of hybridization between R. cerasi and R. cingulata were found in 2% of R. cingulata individuals and in 3% of R. cerasi. All putative hybrids had R. cerasi mtDNA indicating that the original between‐species mating involved R. cerasi females and R. cingulata males.     

Evidence for spatial clines and mixed geographic modes of speciation for North American cherry‐infesting Rhagoletis (Diptera: Tephritidae) flies

An important criterion for understanding speciation is the geographic context of population divergence. Three major modes of allopatric, parapatric, and sympatric speciation define the extent of spatial overlap and gene flow between diverging populations. However, mixed modes of speciation are also possible, whereby populations experience periods of allopatry, parapatry, and/or sympatry at different times as they diverge. Here, we report clinal patterns of variation for 21 nuclear‐encoded microsatellites and a wing spot phenotype for cherry‐infesting Rhagoletis (Diptera: Tephritidae) across North America consistent with these flies having initially diverged in parapatry followed by a period of allopatric differentiation in the early Holocene. However, mitochondrial DNA (mtDNA) displays a different pattern; cherry flies at the ends of the clines in the eastern USA and Pacific Northwest share identical haplotypes, while centrally located populations in the southwestern USA and Mexico possess a different haplotype. We hypothesize that the mitochondrial difference could be due to lineage sorting but more likely reflects a selective sweep of a favorable mtDNA variant or the spread of an endosymbiont. The estimated divergence time for mtDNA suggests possible past allopatry, secondary contact, and subsequent isolation between USA and Mexican fly populations initiated before the Wisconsin glaciation. Thus, the current genetics of cherry flies may involve different mixed modes of divergence occurring in different portions of the fly''s range. We discuss the need for additional DNA sequencing and quantification of prezygotic and postzygotic reproductive isolation to verify the multiple mixed‐mode hypothesis for cherry flies and draw parallels from other systems to assess the generality that speciation may commonly involve complex biogeographies of varying combinations of allopatric, parapatric, and sympatric divergence.     

Isolation and characterization of microsatellite loci from the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae)

We report the isolation and development of 81 novel primers for amplifying microsatellite loci in the Rhagoletis pomonella sibling species complex, and the sequencing, characterization and analysis of basic population genetic parameters for nine of these genes. We also report the successful cross‐species amplification of several of these loci. The R. pomonella sibling species complex is a textbook example of genetic differentiation in sympatry via host‐plant shifting. Microsatellite markers can be useful for mapping host‐plant‐associated adaptations in Rhagoletis that generate reproductive isolation and facilitate speciation, as well as for resolving the genetic structure and evolutionary history of fly populations.     

Occurrence and Genetic Characteristics of Cryptosporidium hominis and Cryptosporidium andersoni in Horses from Southwestern China

A total of 333 fecal specimens from horses in southwestern China were genotyped based on analysis of the small subunit rRNA (SSU rRNA) gene. Cryptosporidium hominis and Cryptosporidium andersoni were identified in 2 and 4 stool specimens, respectively. The identification of C. hominis was confirmed by sequence analysis of the 70‐kDa heat shock protein (HSP70) and oocyst wall protein (COWP) genes. Subtyping analysis of the 60‐kDa glycoprotein (GP60) gene sequence of C. hominis revealed a new rare subtype Id, named IdA15; only three Id isolates have been reported in humans to date. Multilocus sequence typing (MLST) analysis indicated that the C. andersoni subtype was A6, A5, A2, and A1 at the four minisatellite loci (MS1, MS2, MS3, and MS16, respectively). This is the first report to identify the presence of C. andersoni and C. hominis in horses in southwestern China and the first to identify a rare zoonotic subtype Id of C. hominis in horses. These findings suggest that infected horses may act as potential reservoirs of Cryptosporidium to transmit infections to humans.     

The hitchhiker's guide to Europe: the infection dynamics of an ongoing Wolbachia invasion and mitochondrial selective sweep in Rhagoletis cerasi

Wolbachia is a maternally inherited and ubiquitous endosymbiont of insects. It can hijack host reproduction by manipulations such as cytoplasmic incompatibility (CI) to enhance vertical transmission. Horizontal transmission of Wolbachia can also result in the colonization of new mitochondrial lineages. In this study, we present a 15‐year‐long survey of Wolbachia in the cherry fruit fly Rhagoletis cerasi across Europe and the spatiotemporal distribution of two prevalent strains, wCer1 and wCer2, and associated mitochondrial haplotypes in Germany. Across most of Europe, populations consisted of either 100% singly (wCer1) infected individuals with haplotype HT1, or 100% doubly (wCer1&2) infected individuals with haplotype HT2, differentiated only by a single nucleotide polymorphism. In central Germany, singly infected populations were surrounded by transitional populations, consisting of both singly and doubly infected individuals, sandwiched between populations fixed for wCer1&2. Populations with fixed infection status showed perfect association of infection and mitochondria, suggesting a recent CI‐driven selective sweep of wCer2 linked with HT2. Spatial analysis revealed a range expansion for wCer2 and a large transition zone in which wCer2 splashes appeared to coalesce into doubly infected populations. Unexpectedly, the transition zone contained a large proportion (22%) of wCer1&2 individuals with HT1, suggesting frequent intraspecific horizontal transmission. However, this horizontal transmission did not break the strict association between infection types and haplotypes in populations outside the transition zone, suggesting that this horizontally acquired Wolbachia infection may be transient. Our study provides new insights into the rarely studied Wolbachia invasion dynamics in field populations.     

Molecular characterization,virulence and horizontal transmission of Beauveria pseudobassiana from Dendroctonus micans (Kug.) (Coleoptera: Curculionidae)

The great spruce bark beetle, Dendroctonus micans (Kugelann) (Coleoptera: Curculionidae), has been a potential threat for Turkey and the entire Eurasian spruce forests for many years. Control strategies which have been applied so far are still insufficient to prevent its damage. A previous study has shown that a Beauveria isolate (ARSEF 9271) proved to be an efficient microbial control agent against the great spruce bark beetle. In this study, this isolate was identified as B. pseudobassiana based on the partial sequence of EF1‐α and ITS sequence. A conidial suspension (1 × 10⁸/ml) of this fungus caused 100% mortality on both larvae and adults of D. micans within 5 and 6 days, respectively. Also, it caused 100% mycosis value on both larvae and adults. Mortality values of horizontal transmission experiments between larvae and adults which were contaminated with 1 × 10⁶/ml spore suspension at 25%, 50%, 75% and 100% rates were determined as 100% after 15 days at 20°C under the laboratory conditions. We also determined the decrease of the damage in spruce wood block (15 × 25 cm) when the contamination rate of the larvae was increased. Our results indicate that B. pseudobassiana ARSEF 9271 seems to be a very promising biocontrol agent against D. micans.     

Global genetic diversity,lineage distribution,and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well‐known example of a worldwide‐distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.     

Obligate annual and successive facultative diapause establish a bet‐hedging strategy of Rhagoletis cerasi (Diptera: Tephritidae) in seasonally unpredictable environments

To cope with temporal and spatial heterogeneity of habitats, herbivorous insects in the temperate zone usually enter diapause that facilitates synchronization of their life cycle with specific stages of host plants, such as fruit ripening. In the present study, we address those factors regulating dormancy responses as part of a ‘longer strategy’ to persist and thrive in temperate environments, focusing on Rhagoletis cerasi, a univoltine, oligophagous species, which overwinters as pupae and emerges when host fruits are available for oviposition at local scale. To ensure population survival and reproduction at habitats with ecological heterogeneity, R. cerasi has evolved a sophisticated diapause strategy based on a combination of local adaptation and diversified bet‐hedging strategies. Diapause duration is determined both by (i) the adaptive response to local host fruit phenology patterns (annual diapause) and (ii) the plastic responses to unpredictable inter‐annual (temporal) climatic variability that drives a proportion of the populations to extend dormancy by entering a second, successive, facultative cycle of prolonged diapause as part of a bet‐hedging strategy. Besides the dormant periods, post‐diapause development (which varies among populations) exerts ‘fine tune’ adjustments that assure synchronization and may correct possible errors. Adults emerging from pupae with prolonged diapause are larger in body size compared with counterparts emerging during the first year of diapause. However, female fecundity rates are reduced, followed by an extended post‐oviposition period, whereas adult longevity remains unaffected. Overall, it appears that R. cerasi populations are adapted to ecological conditions of local habitats and respond plastically to unpredictable environmental (climatic) conditions.     

Identification of Belgian mosquito species (Diptera: Culicidae) by DNA barcoding

Since its introduction in 2003, DNA barcoding has proven to be a promising method for the identification of many taxa, including mosquitoes (Diptera: Culicidae). Many mosquito species are potential vectors of pathogens, and correct identification in all life stages is essential for effective mosquito monitoring and control. To use DNA barcoding for species identification, a reliable and comprehensive reference database of verified DNA sequences is required. Hence, DNA sequence diversity of mosquitoes in Belgium was assessed using a 658 bp fragment of the mitochondrial cytochrome oxidase I (COI) gene, and a reference data set was established. Most species appeared as well‐supported clusters. Intraspecific Kimura 2‐parameter (K2P) distances averaged 0.7%, and the maximum observed K2P distance was 6.2% for Aedes koreicus. A small overlap between intra‐ and interspecific K2P distances for congeneric sequences was observed. Overall, the identification success using best match and the best close match criteria were high, that is above 98%. No clear genetic division was found between the closely related species Aedes annulipes and Aedes cantans, which can be confused using morphological identification only. The members of the Anopheles maculipennis complex, that is Anopheles maculipennis s.s. and An. messeae, were weakly supported as monophyletic taxa. This study showed that DNA barcoding offers a reliable framework for mosquito species identification in Belgium except for some closely related species.     

Variation of within‐crop microhabitat use by Drosophila suzukii (Diptera: Drosophilidae) in blackberry

Drosophila suzukii Matsumara (Diptera: Drosophilidae) is an invasive vinegar fly that infests ripe and ripening soft skinned fruits. In the south‐eastern United States, blackberry (Rubus spp.) crops are heavily impacted by D. suzukii, and current management tactics rely on the use of broad‐spectrum insecticides targeted to adult populations. An improved understanding of D. suzukii biology and ecology are necessary to create sustainable management options. Knowledge of how D. suzukii interacts with resources will enable targeted management actions in the future. In this present study, we monitored larval infestation throughout the blackberry canopy and found that infestation was highest in the inner portion of the canopy and lower in more exposed locations. We also documented higher humidity within the cane canopy relative to the edge of the field. A difference in abiotic conditions may create within‐crop microhabitats that D. suzukii is able to exploit. Future research will explore how to take advantage of these microhabitats in pest management programs.     

Behavioural tactics used by invasive cane toads (Rhinella marina) to exploit apiaries in Australia

Behavioural flexibility plays a key role in facilitating the ability of invasive species to exploit anthropogenically‐created resources. In Australia, invasive cane toads (Rhinella marina) often gather around commercial beehives (apiaries), whereas native frogs do not. To document how toads use this resource, we spool‐tracked cane toads in areas containing beehives and in adjacent natural habitat without beehives, conducted standardized observations of toad feeding behaviour, and ran prey‐manipulation trials to compare the responses of cane toads versus native frogs to honeybees as potential prey. Toads feeding around beehives travelled shorter distances per night, and hence used different microhabitats, than did toads from nearby control sites without beehives. The toads consumed live bees from the hive entrance (rather than dead bees from the ground), often climbing on top of one another to gain access to the hive entrance. Prey manipulation trials confirm that bee movement is the critical stimulus that elicits the toads’ feeding response; and in standardized trials, native frogs consumed bees less frequently than did toads. In summary, cane toads flexibly modify their movements, foraging behaviour and dietary composition to exploit the nutritional opportunities created by commercial beehives, whereas native anurans do not.     

Field evaluation of the relative attractiveness of enriched ginger root oil (EGO)lure and trimedlure for African Ceratitis species (Diptera: Tephritidae)

The males of some fruit flies (Diptera: Tephritidae) are known to be attracted to specific parapheromones. The trapping results between trimedlure (TML) and enriched ginger root oil (EGOlure) were compared at two experimental sites in Morogoro (Central Tanzania) for a period of 12 weeks co‐inciding with the main citrus season. Both attractants captured a comparable diversity of fruit flies, except that EGOlure also attracted fruit flies, such as Ceratitis cosyra, not normally found in TML‐baited traps. Both EGOlure and TML attracted mainly or exclusively male fruit flies, but the catches with EGOlure were equal or superior to those with TML. It is concluded that EGOlure should be considered as a suitable alternative for TML in detection, monitoring and control programs for African fruit flies of the genus Ceratitis. It has the added advantage that it combines the attractiveness with regard to species spectrum of both TML and terpinyl acetate.     

Existence of species complex largely reduced barcoding success for invasive species of Tephritidae: a case study in Bactrocera spp.

Fruit flies in the family Tephritidae are the economically important pests that have many species complexes. DNA barcoding has gradually been verified as an effective tool for identifying species in a wide range of taxonomic groups, and there are several publications on rapid and accurate identification of fruit flies based on this technique; however, comprehensive analyses of large and new taxa for the effectiveness of DNA barcoding for fruit flies identification have been rare. In this study, we evaluated the COI barcode sequences for the diagnosis of fruit flies using 1426 sequences for 73 species of Bactrocera distributed worldwide. Tree‐based [neighbour‐joining (NJ)]; distance‐based, such as Best Match (BM), Best Close Match (BCM) and Minimum Distance (MD); and character‐based methods were used to evaluate the barcoding success rates obtained with maintaining the species complex in the data set, treating a species complex as a single taxon unit, and removing the species complex. Our results indicate that the average divergence between species was 14.04% (0.00–25.16%), whereas within a species this was 0.81% (0.00–9.71%); the existence of species complexes largely reduced the barcoding success for Tephritidae, for example relatively low success rates (74.4% based on BM and BCM and 84.8% based on MD) were obtained when the sequences from species complexes were included in the analysis, whereas significantly higher success rates were achieved if the species complexes were treated as a single taxon or removed from the data set – BM (98.9%), BCM (98.5%) and MD (97.5%), or BM (98.1%), BCM (97.4%) and MD (98.2%).