A morphological and molecular analysis was undertaken with the objective of identifying markers for geographical populations of Old World screwworm flies, Chrysomya bezziana Villeneuve (Diptera: Calliphoridae). The morphological analysis involved 192 adult flies from 14 countries, and the molecular analysis involved 45 larvae or adults from 14 populations in 11 countries. Principal components and cluster analysis of 10 morphological characters indicated that flies from Papua New Guinea (PNG) were a distinct group and most similar to flies from nearby Asian islands (Java, Sabah). There was poor resolution of other geographical regions, but some support for clustering of flies from Africa or India. Cladistic analysis of mitochondrial DNA sequences gave strong support for recognizing two races of Old World screwworm, one from sub-Saharan Africa and the other from the Gulf region and Asia. This latter race could be further divided into two lineages, i.e. one from mainland Asia (from Iraq to the Malay Peninsula) and the other from two islands of PNG.
Phylogenetic, genealogical and population relationships of Chrysomya bezziana, the Old World screwworm fly (OWSF), were inferred from DNA sequences of mitochondrial cytochrome b (cyt b), nuclear elongation factor-1α (EF-1α) and nuclear white eye colour (white), using sequences of Chrysomya megacephala and Chrysomya rufifacies as outgroups. Cyt b (717bp, 754 specimens), EF-1α (361bp, 256 specimens) and white (577bp, 242 specimens) were analysed from up to two African and nine Asian countries, including 10 Indonesian islands. We show that OWSF occurs as distinctive African and Asian lineages based on cyt b and white, and that there is a marked differentiation between Sumatran and Javan populations in Indonesia, supported by the genealogy and analysis of molecular variance of cyt b alone. Four cyt b sub-lineages are recognised in Asia: only 2.1 occurs on the Asian mainland, from Yemen to Peninsular Malaysia; only 2.2, 2.3 and 2.4 occur in central Indonesia; 2.4 predominates on New Guinea; and 2.1 co-occurs with others only on Sumatra in western Indonesia. This phylogeography and the genetic distances between cyt b haplotypes indicate pre-historic, natural dispersal of OWSF eastwards into Indonesia and other Malesian islands, followed by vicariant evolution in New Guinea and central Indonesia. OWSF is absent from Australia, where there is surveillance for importation or natural invasion. Judged by cyt b haplotype markers, there is currently little spread of OWSF across sea barriers, despite frequent shipments of Australian livestock through Indonesian seas to the Middle East Gulf region. These findings will inform plans for integrated pest management, which could be applied progressively, for example starting in East Nusa Tenggara (central Indonesia) where OWSF has regional cyt b markers, and progressing westwards to Java where any invasion from Sumatra is unlikely. Cyt b markers would help identify the source of any re-emergence in treated areas.