Phylogenetic analysis of New World screwworm fly, Cochliomyia hominivorax, suggests genetic isolation of some Caribbean island populations following colonization from South America

Development of a diagnostic single nucleotide polymorphism (SNP) panel for identifying geographic origins of Cochliomyia hominivorax, the New World screwworm

The New World screwworm, Cochliomyia hominivorax, causes myiasis in livestock, humans, and other warm-blooded animals in much of South America and the Caribbean. It has been eradicated from North and Central America using the sterile insect technique and a biological barrier is currently maintained at the Panama - Colombian border. However, C. hominivorax is still a threat to eradicated areas as outbreaks can and do occur. In order to identify the origin of a fly involved in an outbreak scenario, diagnostic tools would be beneficial. Recently, the geographic population structure of this species was identified using single nucleotide polymorphisms (SNPs). Here we characterize the three major regional clusters: South America, the Inner Caribbean, and the Outer Caribbean. The objective of this study was to develop a SNP (single nucleotide polymorphism) panel to distinguish between these three clusters. A panel was developed using two unique SNPs per region for a total of six SNPs. This diagnostic SNP assay will allow for rapid source determination of flies from future incursions in order to intercept introductory pathways and aid in the control of New World screwworm.

Geographic Population Genetic Structure of the New World Screwworm, Cochliomyia hominivorax (Diptera: Calliphoridae), Using SNPs

The New World screwworm, Cochliomyia hominivorax (Coquerel 1858) (Diptera: Calliphoridae), is a serious parasite of livestock, humans, and other warm-blooded animals. It has been eradicated from the northern parts of its historical range down to the Panama-Colombian border where a permanent barrier zone is maintained. This eradication was accomplished through using the sterile insect technique (SIT). In 2016 there was an outbreak of C. hominivorax in the Florida Keys. In only six months, this pest was successfully re-eradicated using SIT, but the geographic origin of the invasion has yet to be resolved. It was previously determined that the Florida flies most likely represented a single invasion, and it was recommended that a finer-scale genetic assessment should be completed. Thus, this current proof-of-concept study aimed to develop a population genetic database using single nucleotide polymorphisms (SNPs) to reference outbreaks and potentially identify the origin of the Florida outbreak. This initial database consists of wild-caught samples from 4 geographic locations as well as laboratory colony samples that originated from 7 additional locations using a genotyping by sequencing (GBS) approach. Geographic population structuring was identified for twelve populations that clustered according to geographic location. The Florida outbreak samples appeared similar to samples from the outer Caribbean cluster which included samples from Dominican Republic and Trinidad and Tobago, however, these results will be further clarified with the replacement of laboratory colony samples with future wild-caught samples.

The toad fly Lucilia bufonivora: its evolutionary status and molecular identification

The blow fly genus Lucilia is composed largely of saprophages and facultative myasis agents, including the economically important species Lucilia cuprina (Wiedemann) (Diptera: Calliphoridae) and Lucilia sericata (Meigen). Only one species is generally recognized as an obligate agent of myiasis, Lucilia bufonivora Moniez, and this is an obligate parasite of toads. Lucilia silvarum (Meigen), a sister species, behaves mainly as a carrion breeder; however, it has also been reported as a facultative parasite of amphibians. Morphologically, these species are almost identical, and historically this has led to misidentification, taxonomic ambiguity and a paucity of studies of L. bufonivora. In this study, dipterous larvae were analysed from toad myiasis cases from the U.K., The Netherlands and Switzerland, together with adult specimens of fly species implicated in amphibian parasitism: L. bufonivora, L. silvarum and Lucilia elongata Shannon (from North America). Partial sequences of two genes, cox1 and ef1α, were amplified. Seven additional blow fly species were analysed as outgroups. Bayesian inference trees of cox1, ef1α and a combined-gene dataset were constructed. All larvae isolated from toads were identified as L. bufonivora and no specimens of L. silvarum were implicated in amphibian myiasis. This study confirms L. silvarum and L. bufonivora as distinct sister species and provides unambiguous molecular identification of L. bufonivora.

Do longer sequences improve the accuracy of identification of forensically important Calliphoridae species?

Species identification is a crucial step in forensic entomology. In several cases the calculation of the larval age allows the estimation of the minimum Post-Mortem Interval (mPMI). A correct identification of the species is the first step for a correct mPMI estimation. To overcome the difficulties due to the morphological identification especially of the immature stages, a molecular approach can be applied. However, difficulties in separation of closely related species are still an unsolved problem. Sequences of 4 different genes (COI, ND5, EF-1α, PER) of 13 different fly species collected during forensic experiments (Calliphora vicina, Calliphora vomitoria, Lucilia sericata, Lucilia illustris, Lucilia caesar, Chrysomya albiceps, Phormia regina, Cynomya mortuorum, Sarcophaga sp., Hydrotaea sp., Fannia scalaris, Piophila sp., Megaselia scalaris) were evaluated for their capability to identify correctly the species. Three concatenated sequences were obtained combining the four genes in order to verify if longer sequences increase the probability of a correct identification. The obtained results showed that this rule does not work for the species L. caesar and L. illustris. Future works on other DNA regions are suggested to solve this taxonomic issue.

Molecular Characterization of the 2016 New World Screwworm (Diptera: Calliphoridae) Outbreak in the Florida Keys

New World screwworm (NWS), Cochliomyia hominivorax (Coquerel 1858) (Diptera: Calliphoridae), is a myiasis-causing fly that can be a serious threat to the health of livestock, wildlife, and humans. Its progressive eradication from the southern United States, Mexico, and Central America from the 1950s to 2000s is an excellent example of successful pest management using sterile insect technique (SIT). In late 2016, autochthonous NWS were detected in the Florida Keys, representing this species' first invasion in the United States in >30 yr. Rapid use of quarantine and SIT was successful in eliminating the infestation by early 2017; however, the geographic source of this infestation remains unknown. Here, we use amplicon sequencing to generate mitochondrial and nuclear sequence data representing all confirmed cases of NWS from this infestation, and compare these sequences to preexisting data sets sampling the native distribution of NWS. We ask two questions regarding the FL Keys outbreak. First, is this infestation the result of a single invasion from one source, or multiple invasions from different sources? And second, what is the geographic origin of this invasion? We found virtually no sequence variation between specimens collected from the FL Keys outbreak, which is consistent with a single source of introduction. However, we also found very little geographic resolution in any of the data sets, which precludes identification of the source of this outbreak. Our lack of success in answering our second question speaks to the need for finer-scale genetic or genomic assessments of NWS population structure, which would facilitate source determination of potential future outbreaks.

DNA-barcoding of forensically important blow flies (Diptera: Calliphoridae) in the Caribbean Region

Correct identification of forensically important insects, such as flies in the family Calliphoridae, is a crucial step for them to be used as evidence in legal investigations. Traditional identification based on morphology has been effective, but has some limitations when it comes to identifying immature stages of certain species. DNA-barcoding, using COI, has demonstrated potential for rapid and accurate identification of Calliphoridae, however, this gene does not reliably distinguish among some recently diverged species, raising questions about its use for delimitation of species of forensic importance. To facilitate DNA based identification of Calliphoridae in the Caribbean we developed a vouchered reference collection from across the region, and a DNA sequence database, and further added the nuclear ITS2 as a second marker to increase accuracy of identification through barcoding. We morphologically identified freshly collected specimens, did phylogenetic analyses and employed several species delimitation methods for a total of 468 individuals representing 19 described species. Our results show that combination of COI + ITS2 genes yields more accurate identification and diagnoses, and better agreement with morphological data, than the mitochondrial barcodes alone. All of our results from independent and concatenated trees and most of the species delimitation methods yield considerably higher diversity estimates than the distance based approach and morphology. Molecular data support at least 24 distinct clades within Calliphoridae in this study, recovering substantial geographic variation for Lucilia eximia, Lucilia retroversa, Lucilia rica and Chloroprocta idioidea, probably indicating several cryptic species. In sum, our study demonstrates the importance of employing a second nuclear marker for barcoding analyses and species delimitation of calliphorids, and the power of molecular data in combination with a complete reference database to enable identification of taxonomically and geographically diverse insects of forensic importance.

Molecular phylogeny of the forensically important genus Cochliomyia (Diptera: Calliphoridae)

Cochliomyia Townsend includes several abundant and one of the most broadly distributed, blow flies in the Americas, and is of significant economic and forensic importance. For decades, Cochliomyia hominivorax (Coquerel) and Cochliomyia macellaria (Fabricius) have received attention as livestock parasites and primary indicator species in forensic entomology. However, Cochliomyia minima Shannon and Cochliomyia aldrichi Del Ponte have only been subject to basic taxonomy and faunistic studies. Here we present the first complete phylogeny of Cochliomyia including numerous specimens per species, collected from 13 localities in the Caribbean. Four genes, the mitochondrial COI and the nuclear EF-1α, 28S rRNA, and ITS2, were analyzed. While we found some differences among gene trees, a concatenated gene matrix recovered a robustly supported monophyletic Cochliomyia with Compsomyiops Townsend as its sister group and recovered the monophyly of Cochliomyia hominivorax, Cochliomyia macellaria and Cochliomyia minima. Our results support a close relationship between Cochliomyia minima and Cochliomyia aldrichi. However, we found Cochliomyia aldrichi containing Cochliomyia minima, indicating recent speciation, or issues with the taxonomy of the group. We provide basic information on habitat preference, distribution and feeding habits of Cochliomyia minima and Cochliomyia aldrichi that will be useful for future forensic studies in the Caribbean.

Traumatic Myiasis: A Neglected Disease in a Changing World

Traumatic myiasis, the parasitic infestation by fly larvae in traumatic lesions of the tissues of living vertebrates, is a serious medical condition in humans and a welfare and economic issue in domestic animals. New molecular studies are providing insights into its evolution and epidemiology. Nevertheless, its incidence in humans is generally underreported, particularly in tropical and subtropical regions. Myiasis in domestic animals has been studied more extensively, but continuous management is difficult and expensive. A key concern is the inadvertent introduction and global spread of agents of myiasis into nonendemic areas, facilitated by climate change and global transport. The incursion of the New World screwworm fly (Cochliomyia hominivorax) into Libya is the most notable of many such range shifts and demonstrates the potential risks of these parasites and the costs of removing them once established in a geographic area. Nevertheless, the insect agents of myiasis can be of societal benefit to forensic science and in medicine as an aid to wound treatment (larval therapy).

Gene flow and genetic structure of Bactrocera carambolae (Diptera, Tephritidae) among geographical differences and sister species, B. dorsalis, inferred from microsatellite DNA data

The Carambola fruit fly, Bactroceracarambolae, is an invasive pest in Southeast Asia. It has been introduced into areas in South America such as Suriname and Brazil. Bactroceracarambolae belongs to the Bactroceradorsalis species complex, and seems to be separated from Bactroceradorsalis based on morphological and multilocus phylogenetic studies. Even though the Carambola fruit fly is an important quarantine species and has an impact on international trade, knowledge of the molecular ecology of Bactroceracarambolae, concerning species status and pest management aspects, is lacking. Seven populations sampled from the known geographical areas of Bactroceracarambolae including Southeast Asia (i.e., Indonesia, Malaysia, Thailand) and South America (i.e., Suriname), were genotyped using eight microsatellite DNA markers. Genetic variation, genetic structure, and genetic network among populations illustrated that the Suriname samples were genetically differentiated from Southeast Asian populations. The genetic network revealed that samples from West Sumatra (Pekanbaru, PK) and Java (Jakarta, JK) were presumably the source populations of Bactroceracarambolae in Suriname, which was congruent with human migration records between the two continents. Additionally, three populations of Bactroceradorsalis were included to better understand the species boundary. The genetic structure between the two species was significantly separated and approximately 11% of total individuals were detected as admixed (0.100 ≤ Q ≤ 0.900). The genetic network showed connections between Bactroceracarambolae and Bactroceradorsalis groups throughout Depok (DP), JK, and Nakhon Sri Thammarat (NT) populations. These data supported the hypothesis that the reproductive isolation between the two species may be leaky. Although the morphology and monophyly of nuclear and mitochondrial DNA sequences in previous studies showed discrete entities, the hypothesis of semipermeable boundaries may not be rejected. Alleles at microsatellite loci could be introgressed rather than other nuclear and mitochondrial DNA. Bactroceracarambolae may be an incipient rather than a distinct species of Bactroceradorsalis. Regarding the pest management aspect, the genetic sexing Salaya5 strain (SY5) was included for comparison with wild populations. The SY5 strain was genetically assigned to the Bactroceracarambolae cluster. Likewise, the genetic network showed that the strain shared greatest genetic similarity to JK, suggesting that SY5 did not divert away from its original genetic makeup. Under laboratory conditions, at least 12 generations apart, selection did not strongly affect genetic compatibility between the strain and wild populations. This knowledge further confirms the potential utilization of the Salaya5 strain in regional programs of area-wide integrated pest management using SIT.

Evidence of weak genetic structure and recent gene flow between Bactrocera dorsalis s.s. and B. papayae, across Southern Thailand and West Malaysia, supporting a single target pest for SIT applications

Background

Bactrocera dorsalis s.s. (Hendel) and B. papayae Drew & Hancock, are invasive pests belonging to the B. dorsalis complex. Their species status, based on morphology, is sometimes arguable. Consequently, the existence of cryptic species and/or population isolation may decrease the effectiveness of the sterile insect technique (SIT) due to an unknown degree of sexual isolation between released sterile flies and wild counterparts. To evaluate the genetic relationship and current demography in wild populations for guiding the application of area-wide integrated pest management using SIT, seven microsatellite-derived markers from B. dorsalis s.s. and another five from B. papayae were used for surveying intra- and inter-specific variation, population structure, and recent migration among sympatric and allopatric populations of the two morphological forms across Southern Thailand and West Malaysia.

Results

Basic genetic variations were not significantly different among forms, populations, and geographical areas (P > 0.05). Nonetheless, two sets of microsatellite markers showed significantly different levels of polymorphisms. Genetic differentiation between intra- and inter-specific differences was significant, but low. Seventeen populations revealed three hypothetical genetic clusters (K = 3) regardless of forms and geographical areas. The genetic structure of sympatric populations slightly changed during the different years of collection. Recent gene flow (m ≥ 0.10) was frequently detected whether samples were sympatric or allopatric. Ninety-five of 379 individuals distributed across the given area were designated as recent migrants or of admixed ancestry. As a consequence of substantial migration, no significant correlation between genetic and geographic distances was detected (R² = 0.056, P = 0.650).

Conclusions

According to the 12 microsatellite variations, weak population structure and recent gene flow suggest that there is no status for cryptic species between B. dorsalis s.s. and B. papayae forms in Southern Thailand and West Malaysia. Both forms can be treated as a single target pest for the SIT program in an area-wide sense. Additionally, the result of species identification based on molecular data and morphological character are not congruent. The use of independent, multiple approaches in the characterization of the target population may ensure the effectiveness and feasibility of SIT-based control in the target area.

Parasite epidemiology in a changing world: can molecular phylogeography help us tell the wood from the trees?

SUMMARY Molecular phylogeography has revolutionised our ability to infer past biogeographic events from cross-sectional data on current parasite populations. In ecological parasitology, this approach has been used to address fundamental questions concerning host-parasite co-evolution and geographic patterns of spread, and has raised many technical issues and problems of interpretation. For applied parasitologists, the added complexity inherent in adding population genetic structure to perceived parasite distributions can sometimes seem to cloud rather than clarify approaches to control. In this paper, we use case studies firstly to illustrate the potential extent of cryptic diversity in parasite and parasitoid populations, secondly to consider how anthropogenic influences including movement of domestic animals affect the geographic distribution and host associations of parasite genotypes, and thirdly to explore the applied relevance of these processes to parasites of socio-economic importance. The contribution of phylogeographic approaches to deeper understanding of parasite biology in these cases is assessed. Thus, molecular data on the emerging parasites Angiostrongylus vasorum in dogs and wild canids, and the myiasis-causing flies Lucilia spp. in sheep and Cochliomyia hominovorax in humans, lead to clear implications for control efforts to limit global spread. Broader applications of molecular phylogeography to understanding parasite distributions in an era of rapid global change are also discussed.

The molecular systematics of blowflies and screwworm flies (Diptera: Calliphoridae) using 28S rRNA, COX1 and EF-1α: insights into the evolution of dipteran parasitism

The Calliphoridae include some of the most economically significant myiasis-causing flies in the world - blowflies and screwworm flies - with many being notorious for their parasitism of livestock. However, despite more than 50 years of research, key taxonomic relationships within the family remain unresolved. This study utilizes nucleotide sequence data from the protein-coding genes COX1 (mitochondrial) and EF1α (nuclear), and the 28S rRNA (nuclear) gene, from 57 blowfly taxa to improve resolution of key evolutionary relationships within the family Calliphoridae. Bayesian phylogenetic inference was carried out for each single-gene data set, demonstrating significant topological difference between the three gene trees. Nevertheless, all gene trees supported a Calliphorinae-Luciliinae subfamily sister-lineage, with respect to Chrysomyinae. In addition, this study also elucidates the taxonomic and evolutionary status of several less well-studied groups, including the genus Bengalia (either within Calliphoridae or as a separate sister-family), genus Onesia (as a sister-genera to, or sub-genera within, Calliphora), genus Dyscritomyia and Lucilia bufonivora, a specialised parasite of frogs and toads. The occurrence of cross-species hybridisation within Calliphoridae is also further explored, focusing on the two economically significant species Lucilia cuprina and Lucilia sericata. In summary, this study represents the most comprehensive molecular phylogenetic analysis of family Calliphoridae undertaken to date.

Enabling technologies to improve area-wide integrated pest management programmes for the control of screwworms

The economic devastation caused in the past by the New World screwworm fly Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae) to the livestock industry in the U.S.A., Mexico and the rest of Central America was staggering. The eradication of this major livestock pest from North and Central America using the sterile insect technique (SIT) as part of an area-wide integrated pest management (AW-IPM) programme was a phenomenal technical and managerial accomplishment with enormous economic implications. The area is maintained screwworm-free by the weekly release of 40 million sterile flies in the Darien Gap in Panama, which prevents migration from screwworm-infested areas in Columbia. However, the species is still a major pest in many areas of the Caribbean and South America and there is considerable interest in extending the eradication programme to these countries. Understanding New World screwworm fly populations in the Caribbean and South America, which represent a continuous threat to the screwworm-free areas of Central America and the U.S.A., is a prerequisite to any future eradication campaigns. The Old World screwworm fly Chrysomya bezziana Villeneuve (Diptera: Calliphoridae) has a very wide distribution ranging from Southern Africa to Papua New Guinea and, although its economic importance is assumed to be less than that of its New World counterpart, it is a serious pest in extensive livestock production and a constant threat to pest-free areas such as Australia. In the 1980s repeated introductions and an expansion of Old World screwworm populations were reported in the Middle East; in the 1990s it invaded Iraq and since late 2007 it has been reported in Yemen, where a severe outbreak of myiasis occurred in 2008. Small-scale field trials have shown the potential of integrating the SIT in the control of this pest and various international organizations are considering using the release of sterile insects as part of an AW-IPM approach on a much wider scale. Wohlfahrtia magnifica (Schiner) (Diptera: Sarcophagidae) is a screwworm of temperate regions, which, although of limited agricultural importance, has invaded several new locations in the past few years. This special issue reports on the results of a 6-year project funded by the Joint Food and Agriculture Organization of the United Nations/International Atomic Energy Agency (FAO/IAEA) Programme of Nuclear Techniques in Food and Agriculture entitled 'Enabling Technologies for the Expansion of the SIT for Old and New World Screwworm'. A major goal of the project was to better understand population genetic variation in screwworms as an aid to the identification of isolated populations. The project also addressed issues related to genetic sexing, cuticular hydrocarbons, population dynamics, genetic transformation and chromosome analysis.