- Open Access
Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae)
BMC Ecology and Evolution volume 21, Article number: 70 (2021)
The common name of the Flesh flies (Sarcophagidae) usually relates them with organisms feeding on decomposing organic matter, although the biology of one of the largest radiations among insects also includes predation, coprophagy, and even kleptoparasitism. The question of whether the ancestor of all sarcophagids was a predator or a decomposer, or in association to which host have sarcophagids evolved, has thus always piqued the curiosity of flesh fly specialists. Such curiosity has often been hindered by both the impossibility of having a well-supported phylogeny of Sarcophagidae and its sister group to trace live habits and the scarcity of information on the biology of the group. Using a phylogenomic dataset of protein-encoding ultraconserved elements from representatives of all three subfamilies of Sarcophagidae as ingroup and a large Calyptratae outgroup, a robust phylogenetic framework and timescale are generated to understand flesh fly systematics and the evolution of their life histories.
The evolutionary history for Sarcophagidae reconstructed here differs considerably from previous hypotheses. Within subfamily Sarcophaginae, a group of predatory flies, including genera Lepidodexia and Boettcheria, emerged as sister-group to the rest of Sarcophaginae. The genera Oxysarcodexia, Ravinia, and Tricharaea, long considered archaic and early-branching coprophagous and sarcosaprophagous lineages, were found nested well within the Sarcophaginae as sister-group to the sarcosaprophagous Microcerella. Predation on invertebrates is suggested as the ancestral and dominant strategy throughout the early evolution of flesh flies. Several transitions from predation to sarcosaprophagy and coprophagy occur across the sarcophagid phylogenetic tree, in contrast with almost no transitions from sarcosaprophagy or coprophagy to predatory habits. Regarding the morphological evolution of flesh flies, there might be a concerted evolution of male genitalia traits, such as the phallotrema position and the juxta, or the vesica and the folding of the phallotrema. One diversification rate shift was inferred in the evolution of sarcophagids, which is related to the origin of genus Sarcophaga.
This study has a significant impact on understanding sarcophagid evolution and highlights the importance of having a robust phylogenetic framework to reconstruct the ancestral character state of biological and morphological characters. I discuss the evolution of life histories of the family in relation to their hosts or substrates and outline how sarcosaprophagy, coprophagy, and kleptoparasitism behavior on various hosts may have evolved from predation on invertebrates. This study provides a phylogenetic framework for further physiological and comparative genomic work between predatory, sarcosaprophagous, coprophagous, and kleptoparasitic lineages, which could also have significant implications for the evolution of diverse life histories in other Diptera.
The flesh flies (Sarcophagidae, ~ 3000 species) have been recognized as one of the largest insect radiations  with an age estimated at 23 Ma  and a myriad of different life habits, whose evolutionary patterns are not yet explained. The common name of flesh flies usually relates them directly with organisms feeding on decomposing organic matter, although their biology is one of the most diverse among insects. In fact, they constitute a diverse assemblage of flies that vary considerably in their biological requirements for larval food resources and feeding habits. Published data of the sarcophagid larval biology document the complexities of their nutritional options and biological relationships [3, 4], which include parasites, parasitoids on various hosts, kleptoparasites (inquilines), predators, sarcosaprophages, and coprophages . Some species are habitat- or host-specific and accurately placed in one of these categories. Others show a more opportunistic or generalist approach, utilizing instead a variety of hosts or several different food sources. However, for many of these species, precise ecological requirements remain essentially unknown. Thus, limited biological information on the group and little understanding of the phylogenetic relationships of sarcophagids limit inferences on the evolutionary pattern of this diverse ecology and its impacts on flesh-fly diversification.
Sarcophagidae species are organized into three subfamilies, Miltogramminae, Paramacronychiinae, and Sarcophaginae . The most diverse subfamily is Sarcophaginae with around 2,000 species in 46 genera  that have colonized almost all terrestrial ecosystems. Only few Sarcophaginae species occur in subarctic regions but most are found in tropical areas, showing a peak of diversity in the Neotropical region . Only the three genera of Sarcophaginae Blaesoxipha Loew, Ravinia Robineau-Desvoidy, and Sarcophaga Meigen have been able to extend their distribution out of the Americas, although the latter one also colonized Afro-Eurasia and Australasia/Oceania and diversified becoming the largest radiation within this subfamily with almost 900 species worldwide [9,10,11].
Flesh flies of Sarcophaginae seem to be a versatile group able to extend their distributions and dominate different feeding substrates. In urban environments, some flesh flies indistinctly feed on garbage and corpses of animals and humans becoming mechanical vectors of important diseases of great importance for public health [12, 13]. Their affinity for decomposing matter of human corpses makes them useful as indicators of time and place of death in forensic investigations [14, 15]. This preference for decomposing organic matter is defined as sarcosaprophagy, which is different from coprophagy, as carrion and dung produce different profiles of volatile organic compounds. Despite of these differences in volatile profiles, some sarcophagid species are able to feed both on carrion and dung. Other sarcophagids seem to have evolved more specialized relationships with their feeding substrates, which sometimes are not as ‘passive’ as carrion and dung, but consist of a broad range of living organisms targeted by flesh flies as preys or hosts. Sarcophagids are able to develop on terrestrial gastropods and arthropods (mostly insects) including millipedes, scorpions, beetles, grasshoppers, cockroaches, and mantids either as predators or parasites [5, 16, 17]. Even more specialized sarcophagids behave as kleptoparasitoids on soil-nesting wasps, bees and ants , a strategy where some species have adapted to usurp other species’ resources instead of collecting their own. Carrion and dung-flies provide ecosystem services such as nutrient recycling [19,20,21], which are essential for the sustainability and well being of urban, rural, and wild ecosystems, while the ecosystem impact and services of predator and kleptoparasitic flies is entirely unclear.
In contrast to other widespread fly groups, flesh flies have been largely neglected in phylogenetic studies. Still, the most representative phylogenetic hypotheses included only a small taxon sampling representing less than one third of the flesh fly species [9, 18, 22,23,24] or they were based on datasets of maximum nine genes reporting very low statistical support [18, 22, 23, 25]. As a result, basic knowledge on phylogenetic relationships is lacking for most flesh fly lineages and their evolutionary history remains unraveled. Specifically, the homology characterization of male terminalia structures, which are essential for species recognition in diversity surveys and species delimitation analyses, and the evolution of larval feeding substrates and habits are still little studied.
Large-scale evolutionary studies are currently dramatically benefitting from modern genomic-based approaches, yet no attempts have been made to apply these methods to elucidate the evolution of sarcophagids. The most popular approaches for understanding biological diversification in time and space include two different strategies of reduced-representation of genome using targeted enrichment, i.e., anchored hybrid enrichment (AHE) [9, 26, 27] and ultraconserved elements (UCE) [28,29,30,31,32,33]. Targeted enrichment is especially popular as it produces large amounts of molecular data from highly fragmented DNA and/or sub-optimally preserved samples such as those coming from museum specimens [34, 35]. Particularly, studies using UCEs have capitalized on museum specimens [30, 31, 33,34,35,36,37] and set the stage for the development of an entire field known as “museomics”. Using protein-encoding UCEs brings the best of two major approaches of reduced representation sequencing strategies, i.e., transcriptomes and UCE. Thus, UCEs can be applied to degraded DNA while transcriptomes ensure targeting protein-encoding genomic regions. The use of protein-encoding UCEs to resolve phylogenetic relationships within dipterans has only started to be explored , and as Sarcophagidae is one the largest radiations of flies that account for the majority of fly life on Earth , they constitute an interesting taxon to be studied with the combined use of the UCE targeted-enrichment method and massively parallel DNA sequencing technologies.
Thus, I aimed to reconstruct the phylogenetic relationships of flesh flies and provide a discussion of the implications for morphological and biological character evolution. In an effort to resolve possible incongruences, I assembled a novel UCE-based dataset of Sarcophagidae and a large representation of related outgroup species of Oestroidea and other Calyptratae, and assessed the robustness of phylogenetic estimates using concatenation and gene tree-based approaches. I additionally conducted an analysis of diversification rates and a reconstruction of ancestral character states (i.e., male terminalia characters and larval natural history) of the major lineages of flesh flies and Oestoidea outgroups.
UCE probes and capture results
All DNA extractions from pinned museum specimens (= 24), specimens preserved in 96% ethanol (= 66) and liquid nitrogen (= 19) succeeded. Similarly, all existing DNA aliquots (= 32) had ample well-preserved DNA (above 1 ng/µL) for the present study. Thus, 141 DNA extractions were enriched using the UCE probes and sequenced on two Illumina Hi-Seq lanes. 17% of the sampled taxa consisted of pinned museum specimens (3–54 years old), of which 100% resulted in successful UCE enrichment. I recovered more UCE loci than the average for 11 pinned, dried specimens with age range of 16–48 years. See Additional file 1 to compare specimen age versus total DNA extracted and UCEs captured coded for preservation method.
Sequencing of libraries produced an average of 2,374,209 raw paired-end reads per sample. Trinity assembled reads into 1070–75,722 contigs with average of 9,608 contigs assembled per sample. These contigs had average lengths of 242.9–516.3 bp. From the total assembled contigs, a total of 2018 UCE loci out of 2,581 UCE targets were recovered across all taxa with an average of 1214 UCE loci per sample and average lengths ranging from 230.5 to 625 bp. Summary results of empirically generated UCE data processed are presented in Additional file 1.
I analyzed datasets containing 64–2018 UCE loci from 141 representatives, including 110 ingroup species of Sarcophagidae and 31 outgroup species from 10 Calyptratae families, with a total concatenated aligned length ranging between 20 and 551 kB.
Phylogenetic relationships for Sarcophagidae and 10 other Calyptratae fly families were inferred from 17 datasets having varying loci occupancy and coded as nucleotides and amino acids using a concatenated ML approach as well as by reconstructing a species tree estimated from UCE gene trees. Both methods returned two topologies largely congruent at genus- and species-level relationships, but topology A received higher branch support (BS, Bootstrap support = 96–100) than topology B (Fig. 1a, b). Most datasets, i.e., 12 out of 17, returned the topology A (Table 1). Datasets producing topology A included 1271–2018 UCE loci (Fig. 1a), while only one dataset, which included 59 UCE loci coded as nucleotides, produced topology B (Fig. 1b). An alternative topology, which is an intermediate between topologies A and B, was produced by four datasets having 288–936 UCE loci coded as nucleotides or as amino acids (Fig. 1a, b). For illustration purposes, the most recurrently recovered topology (topology A), i.e., the coalescent-based tree reconstructed from the dataset containing 2018 UCE loci coded as nucleotides, which received strong statistical support, is depicted in Fig. 2.
Overall, datasets including more UCE loci (1271–2018) produced topologies that were more congruent among them, while datasets containing fewer UCE loci (59–936) (Table 1) produced topologies showing alternative relationships at genus- and subgenus-level. Similarly, UCE loci coded as amino acids, which by principle produce datasets with fewer characters, produced alternative topologies that were incongruent with the phylogenetic relationships recovered using datasets that included more characters. Resulting topologies from all of the analyses are included in Additional file 2.
Within Sarcophagidae, the subfamilies were recovered as monophyletic. Among the Sarcophagidae genera represented by two or more species, only the genera Peckia Robineau-Desvoidy and Titanogrypa Townsend were not recovered as monophyletic. Within the genus Sarcophaga, among subgenera represented by two or more species, only subgenera Helicophagella Enderlein and Liopygia Enderlein were recovered as paraphyletic in all analyses.
Subfamily-level relationships within Sarcophagidae received full support (BS = 100, local posterior probability (LPP) = 100) and were identical across all phylogenetic analyses, with Sarcophaginae as sister to (Miltogramminae + Paramacronychiinae). Regarding genus-level relationships, all analyses produced topologies that were largely congruent and well supported (BS = 89–100, LPP = 100) with two exceptions, as follows: (1) genus Emdenimyia Lopes as sister to (Tripanurga Brauer & Bergenstamm + Boettcheria Parker) in 12 out of 17 analyses (Fig. 1a) or to (Lepidodexia Brauer & Bergenstamm (Tripanurga + Boettcheria)) in 5 out of 17 analyses (Fig. 1b); (2) genus Chrysagria Townsend as sister to (Peckiamyia Dodge + Retrocitomyia Lopes) in 16 out of 17 analyses (Fig. 1a) or to (Helicobia Coquillett ((Chrysagria (Peckiamyia + Retrocitomyia)) + (remaining Sarcophaginae))) in 1 out of 17 analyses (Fig. 1b). Regarding subgenus- and species-level relationships, all analyses produced topologies that were largely congruent, with few exceptions, all of them within the genus Sarcophaga. Seven different topologies are reconstructed for this genus. One of these topologies is overall well supported and recurrent among analyses, which is depicted in Fig. 2. Differences among the seven topologies for Sarcophaga include (1) variable phylogenetic positions of clade H containing species of subgenera Neosarcophaga Shewell and Myorhina Robineau-Desvoidy (among others) and clade J containing subgenera Bellieriomima Rohdendorf and Helicophagella (among others) that are recovered as sister groups or as laddered sister-groups to the remaining Sarcophaga, and (2) the position of subgenus Parasarcophaga Johnston & Tiegs as sister to (Robineauella Enderlein (Liosarcophaga Enderlein + Rosellea Rohdendorf)) or to (Liosarcophaga + Rosellea) in clade M (Fig. 2, Additional file 2).
Using the function makeChronosCalib in the R package ape, I estimated the timescale of the evolution of Sarcophagidae. The estimation of ages used the topology A from the coalescent-based tree reconstructed from the dataset containing 2018 UCE loci and coded as nucleotides. The estimated age for the most recent common ancestor (MRCA) of the clade (Miltogramminae + Paramacronychiinae) is 19.6 Ma and for Sarcophaginae is 21.8 Ma, while the estimated age for the MRCA of the subfamily Miltogramminae is 7.8 Ma and for Paramacronychiinae is 14.9 Ma. The estimated age for the MRCA of genus Sarcophaga, the richest of the whole family, is 14.1 Ma. The Miocene was an important setting in generating the presently recognized generic diversity of Sarcophagidae, as the majority of genera appear to have evolved between 1 and 15 Ma. However, these age estimates include considerable uncertainty and should be read with caution, given the lack of fossils for sarcophagids and closely related fly families.
Ancestral character reconstruction of life history and male terminalia traits
The ancestral character state reconstruction (ACR) for life-history traits using Mesquite suggests that phytophagy and predation are the most parsimonious states for the MRCA of Oestroidea, while the most parsimonious state for the MRCA of the clade containing the Oestroidea families Calliphoridae, Polleniidae, Rhiniidae, Rhinophoridae, and Tachinidae is predation on invertebrates (Fig. 3). The most parsimonious state for the MRCA of (Oestridae + Sarcophagidae) is predation on invertebrates or vertebrates. Similarly, the most parsimonious state for the MRCA of Sarcophagidae, and for each of its subfamilies, is predation on invertebrates. Within subfamily Sarcophaginae, several MRCAs have made use of other resources than invertebrates. For example, the most parsimonious state for the MRCA of the clade B is a diverse set of phenotypes including feces, invertebrates, and vertebrates as larval food resources and coprophagy, sarcosaprophagy and, predation as larval feeding habits. More specialized phenotypes are estimated within this clade, with coprophagy as the most parsimonious state for the MRCA of the clade (Oxysarcodexia Townsend + Oxyvinia Dodge). Similarly, a specialized phenotype is estimated for the clade containing genera Engelimyia Lopes, Lipoptilocnema Townsend, and Peckia, whose MRCA has sarcosaprophagy on vertebrates as the most parsimonious state. The most parsimonious state for the MRCA of the genus Sarcophaga includes predation and sarcosaprophagy on invertebrates.
In the corHMM-based ACR, the two evaluated models (ARD = all transition rates are allowed and are independently estimated; ER = all transitions among the specified number of rate classes are the same) provided a differential best fit to the life history and morphological characters (Additional file 3). I found stronger support for a MRCA of Sarcophagidae using invertebrates than vertebrates as larval food resource (Fig. 4a), and overwhelming support for a MRCA of Sarcophagidae with predation as larval feeding habit (Fig. 4b). There is further strong support for invertebrates as larval food resource of the MRCA of each of the Sarcophagidae subfamilies, and very strong support for predation as the larval feeding habit of the MRCA of Sarcophaginae while slightly weaker support for predation as the larval feeding habit of the MRCA of Miltogramminae and Paramacronychiinae. Within the Sarcophaginae, transitions in larval food resources are observed for the clade B (Fig. 4a). For example, the MRCA of this clade is possibly a predator on invertebrates, while the MRCA of genus Microcerella Macquart is a sarcosaprophage on vertebrates and the MRCA of the remaining taxa in clade C is possibly a predator or coprophagous. Similarly, the MRCA of the clade containing genera Engelimyia, Lipoptilocnema, Peckia, and Sarcophaga is most likely a predator on invertebrates, while the MRCA of the subclade containing genera Engelimyia, Lipoptilocnema, and Peckia is a sarcosaprophage on vertebrates. There is further strong support for the MRCA of genus Sarcophaga as a predator on invertebrates, with several subsequent transitions towards sarcosaprophagy on vertebrates and coprophagy. Thus, predominant sarcosaprophagy and coprophagy have evolved independently only a few times within the subfamily Sarcophaginae. Similarly, few taxa have adapted to a combination of life habits including coprophagy, sarcosaprophagy, and predation from a mostly predator ancestor. Such is the case of some species of subgenera of Sarcophaga, like Bercaea Robineau-Desvoidy, Helicophagella, Liopygia, Neobellieria Blanchard, Parasarcophaga, and Thyrsocnema Enderlein, which are subgenera usually having one widespread species with a variety of records as coprophage, sarcosaprophage and/or predator, while the remaining species of these subgenera have species with more specialized habits, but in many cases, species with less available information on their habits.
Regarding the corHMM-based ACR for male terminalia traits, I found strong support for the MRCA of Miltogramminae and Paramacronychiinae having the posterior margin of the abdominal ST5 straight or with a shallow concavity, while it forms a cleft in the MRCA of Sarcophaginae (Additional file 4). Within the later, there is strong support for a transition from a posterior margin of the abdominal ST5 forming a cleft towards a straight or a shallow concavity shape in the MRCA of the clade containing Cistudinomyia Townsend, Dexosarcophaga Townsend, Nephochaetopteryx Townsend, Oxysarcodexia, Oxyvinia, Sarcofahrtiopsis Hall, and Tricharaea Thomson. Regarding the outline of dorsal surface of cercal prong in lateral view, the MRCA of the three subfamilies have a straight or almost straight outline (Additional file 5), while the MRCA of the clade E has the outline of dorsal surface of cercal prong in lateral view with a proximal hump. There is support for a transition from a straight or almost straight cercal outline to an outline with a subapical saddle-shaped concavity followed by a hump in the MRCA of the clade containing genera Engelimyia, Lipoptilocnema, Peckia, Sarcophaga, and Villegasia Dodge. Regarding the connection between basi- and distiphallus, there is strong support for the MRCA of Miltogramminae and Paramacronychiinae having a continuous connection, while this connection is non-continuous in the MRCA of Sarcophaginae and continuous in the MRCA of the clade C (Additional file 6). This connection between basi- and distiphallus is studied in more detail with the inclusion of the character shape of connection between basi- and distiphallus, for which the ACR recovers a very strong support for a serial transition between the presence of a distinct hinge between basi- and distiphallus, a partially sclerotized connection and a fully sclerotized connection. More specifically, strong support is obtained for the MRCA of the clade B having a distinct hinge between basi- and distiphallus, while the MRCA of genus Microcerella has a partially sclerotized connection and the MRCA of the clade C has a fully sclerotized connection (Fig. 5). The ACR for the harpes gives numerous independent origins for these structures within Sarcophaginae (Additional file 7), while the ACR for the vesica provides a strong support for a single origin of this structure in the MRCA of Sarcophaginae, and a reversal to absence of vesica in the MRCA of genus Blaesoxipha (Additional file 8). Similarly, a folded phallotrema, forming three openings, seems to have evolved only once in the MRCA of Sarcophagidae (Additional file 9). There is support for a phallotrema placed in a ventral position with regard to the phallic tube in the MRCA of all sarcophagids, with the independent evolution of an apical phallotrema in the MRCA of Miltogramminae (Additional file 10). The ACR for the acrophallic levers supports only one origin in the MRCA of the clade C (Additional file 11). The number of styli of the phallus is estimated to be one in the MRCA of (Miltogramminae + Paramacronychiinae) while three in the MRCA of Sarcophaginae, with a reduction to two styli in the MRCA of the large clade F (Fig. 6a). There is support for a single origin of the capitis in the MRCA of Sarcophagidae, with two independent reductions in the MRCA of Miltogramminae and in four subgenera of the genus Peckia (Additional file 12). Regarding the evolution of the median process, the ARC provides support for a single reduction of this structure in the MRCA of the large clade F (Additional file 13). Finally, the juxta is estimated to have evolved only once in the ancestor of all sarcophagids and have reduced only once in the MRCA of the subfamily Miltogramminae (Additional file 14).
Genus-level diversification of Sarcophagidae
Rates of diversification were estimated in Sarcophagidae to investigate whether significant rate shifts have occurred over time, using (1) the time calibrated topology from the coalescent-based tree reconstructed from the dataset containing 2018 UCE loci and coded as nucleotides and (2) a set of sampling fractions based on species estimates (Additional file 15). I plotted the mean marginal posterior density of speciation rates in Fig. 6b. This phylorate graph shows an increase in estimated rate along the branch leading to genus Sarcophaga (semi-transparent black circle in Fig. 6b) excluding subgenus Neobellieria (light blue branch in Fig. 6b). The highest rate increase occurs along the branch leading to a clade containing several subgenera of Sarcophaga (indicated with an asterisk on Fig. 6b). All of the top nine most credible rate shift configurations support an increase in estimated rate along the branch leading to or within the genus Sarcophaga (Fig. 7). Some of the top nine most credible rate shift configurations support additional rate increases in other branches but with very low frequency (= posterior probability support) (f = 0.06–0.0098) (Fig. 7c–e, g–i). Moreover, the frequency for the best shift configuration was high (f = 0.57), and the alternative eight best shift configurations were generally consistent with this (Fig. 7).
The family Sarcophagidae is a highly diverse group, being one of the largest insect radiations among the living organisms . The present study corroborates previous findings regarding the monophyly of sarcophagids [18, 23,24,25, 33, 39,40,41,42,43] and its three subfamilies [6, 9, 10, 18, 25, 33]. Several studies have challenged the traditional classification and questioned subfamily-level relationships [25, 40, 44], but more evidence is accumulating in support of Paramacronychiinae as sister to Miltogramminae [8, 9, 18, 23, 33]. Regarding the phylogenetic relationships within Miltogramminae and Paramacronychiinae, the relationships recovered here are consistent with those estimated using other molecular datasets. For example, Eumacronychia Townsend is sister to the remaining Miltogramminae  and the genera Dexagria Rohdendorf, Paramacronychia Brauer & Bergenstamm, and Brachicoma Rondani are early-branching lineages within Paramacronychiinae .
The relationships within the subfamily Sarcophaginae have historically been a challenge for phylogeneticists. Previous morphology-based studies suggested phylogenetic hypotheses for the genera of this subfamily [7, 45], which are not fully supported by Sanger-based phylogenies [22, 23, 25]. Most of these studies have important differences in taxon sampling and molecular markers, and many of them received weak statistical support. More recent next generation sequencing (NGS)-based phylogenies using transcriptomes , anchored hybrid enrichment , and UCEs  received stronger statistical support, but had limited taxon sampling for genera. In the present study, the use of UCE data augmented by improved taxon sampling results in a much better resolved phylogeny, with most nodes receiving maximum support (Fig. 2).
The new phylogeny does not generally contradict well-supported clades in the previous NGS-based studies but does uncover many novel relationships that are revealed by the inclusion of a wide sampling of genera. These novel relationships include (1) the placement of the clade A or (Lepidodexia (Emdenimyia (Boettcheria + Tripanurga))) as sister to the remaining sarcophagids, (2) the sister-group relationship between Microcerella and the clade C (Fig. 2), (3) the sister-group relationship between Cistudinomyia and Tricharaea, (4) the sister-group relationship between Dexosarcophaga and (Nephochaetopteryx + Sarcofahrtiopsis), (5) the sister-group relationship between Argoravinia Townsend and (Malacophagomyia Lopes (Titanogrypa + Rafaelia Townsend + Tulaeopoda Townsend)) in clade D, (6) the sister-group relationship between (Comasarcophaga (Mecynocorpus + Fletcherimyia)) and (Emblemasoma Aldrich (Spirobolomyia Townsend + Blaesoxipha)), (7) the sister-group relationship between Helicobia and a large clade containing genera Chrysagria, Peckiamyia, Retrocitomyia, Villegasia, Peckia, Engelimyia, Lipoptilocnema, and Sarcophaga, and (8) the sister-group relationship between (Peckia (Pattonella) Enderlein (Peckia (remaining subgenera) (Engelimyia + Lipoptilocnema))) and Sarcophaga. These results partially agree with phylogenetic relationships for Sarcophaginae based on morphological and molecular characters. For example, the morphology-based ‘Blaesoxipha clade’ of Buenaventura & Pape  is represented here by Blaesoxipha, Comasarcophaga, Fletcherimyia, Mecynocorpus, and Spirobolomyia and is recovered as a well-supported clade, although the clade recovered here also includes Emblemasoma. The internal relationships of this clade support those previously recovered by morphology  and UCEs  with Blaesoxipha as sister to Spirobolomyia and Comasarcophaga closely related to Fletcherimyia. The results partially support the morphology-based ‘Sarcophaga clade’ of Buenaventura & Pape , which includes genera Chrysagria, Engelimyia, Helicobia, Lipoptilocnema, Peckia, and Sarcophaga, although the clade recovered here also includes Peckiamyia, Retrocitomyia, and Villegasia. The present study also supports the existence of clades that previously had only weak or ambiguous support such as a Boettcheria + Tripanurga  and Peckiamyia + Retrocitomyia [7, 46].
Within the largest radiation of Sarcophaginae, phylogenetic relationships in the hyperdiverse genus Sarcophaga closely match those recovered in previous molecular phylogenies, with the Nearctic subgenus Neobellieria as sister-group to the remaining Sarcophaga species [10, 11, 33]. An early-branching clade including subgenera Heteronychia Brauer & Bergenstamm and Brasia Strand is also supported by data from an anchored hybrid enrichment analysis . The clade (Sarcophaga + Thyrsocnema) has also been supported in previous molecular phylogenies [9,10,11, 33], as well as its close relationship to Myorhina [9, 11]. For the first time the Nearctic endemic subgenus Neosarcophaga is included in a phylogenetic study and all analyses recovered it in clade H as the sister to the Palaearctic clade (Myorhina (Sarcophaga + Thyrsocnema)). The relationships within clade J with subgenera Bellieriomima, Helicophagella, Mauritiella Verves, and Asceloctella Enderlein (among others) and those of clade K with subgenera Pandelleana Rohdendorf and Sarcorohdendorfia Baranov (among others) highly resemble those recovered in previous NGS-based phylogenies [9, 33]. Within clade K, the sister-group relationship between Stackelbergeola Rohdendorf and Rohdendorfisca Grunin is also consistent with previous Sanger- and NGS-based phylogenies [9, 33, 47]. A close relationship between Bercaea and Liopygia has been suggested in previous phylogenies , which is supported here with Liopygia rendered paraphyletic by Bercaea. Similarly, a close relationship between subgenera Liosarcophaga, Parasarcophaga, Robineauella, and Rosellea is supported by a previous molecular analysis  and it is here confirmed with clade M that has (Robineauella (Parasarcophaga (Liosarcophaga + Rosellea))). Most of the similar phylogenetic results of the present study and those of previous NGS-based phylogenies are due to taxon sampling compatibility and not to data duplication or data similarity, as the loci analyzed in these studies are not compatible but they give a significantly consistent phylogenetic signal.
Evolution of life habits
Flesh flies successfully feed on a breadth of live or dead hosts (vertebrates, invertebrates) and substrates (feces), yet these vary tremendously in accessibility, volatile profile, attractiveness to potential competitors, as well as the set of behavioral, chemical, and morphological specializations needed for the gravid females and the first instar larvae to survive interactions with their hosts. A significant distinction can usually be made between a community of predatory species and a community of decomposer species living as sarcosaprophages, coprophages, or kleptoparasitic flies. There is a general fidelity in the choice of larval feeding habits, with some flesh-fly genera adapted to behave as predators on invertebrates (Blaesoxipha, Boettcheria, Chrysagria, Cistudinomyia, Emblemasoma, Emdenimyia, Lepidodexia, Malacophagomyia, Rafaelia, Spirobolomyia, Tripanurga, and several subgenera of Sarcophaga) and a group of genera adapted as decomposers of organic matter (sarcosaprograges and coprophages) either on invertebrates or vertebrates (Engelimyia, Lipoptilocnema, Microcerella, Nephochaetopteryx, Peckia, Peckiamyia, Oxysarcodexia, Oxyvinia, Ravinia, Retrocitomyia, Sarcofahrtiopsis, Tricharaea, and Villegasia). Very few Sarcophagidae taxa can be considered generalists, as usually only one or two species in particular genera are able to feed on different trophic substrates according to their availability. However, it has historically been unclear how this manifold ecology of predators, sarcosaprograges, coprophages, and kleptoparasites evolved within the flesh flies.
The evolutionary patterns giving origin to the diverse ecology of sarcophagids and hosts to which these flies have been associated in their evolution have not been fully addressed. The lack of a well-supported phylogeny of Sarcophagidae and its sister family (Oestridae) has been the most limiting factor. However, the scarcity of information on the biology of the group has also been a relevant limitation. The possibility of having a robust phylogeny for Sarcophagidae has improved over the years with the advance of phylogenetic methods and the efficiency in NGS-based molecular techniques and morphological methods for accessing and collecting different sources of data for phylogenetic reconstruction. The knowledge on the biology of the group, especially of some genera, continues to be a limitation, which can be mitigated with the use of algorithms for ancestral character state reconstruction that allow for missing data. However, this comes with costs in the uncertainty of the estimations.
Earlier hypotheses suggested that at least the subfamily Sarcophaginae had a coprophage or saprophage ancestor, and was undergoing the change from coprophagous and saprophagous nutrition to pure parasitic nutrition . Regarding subfamily Miltogramminae, a recent hypothesis suggests that larvae of ancient miltogrammines were sarcosaprophagous . An ancestral sarcosaprophagous Sarcophagidae fly feeding on invertebrates as primary food source was recently supported by the first formal analysis of this question, and the ancestors of the subfamilies Miltogramminae and Paramacronychiinae were further reconstructed as sarcosaprophagous feeding on both vertebrates and invertebrates . These hypotheses are only partially supported by the data presented here. The results presented here support kleptoparasitism (or inquilinism) of Miltogramminae arising from a predator (i.e., parasitic) ancestor [49, 50]. In contrast to previous hypotheses, the MRCA of all sarcophagids is here estimated to have been a predator on invertebrates, which is consistently supported by both the Mesquite-based (including the sister-family Oestridae, among other outgroups) and the corHMM-based (including only Sarcophagidae) ancestral character state reconstructions. This same combination of larval feeding habit and larval food resource is also supported for the MRCA of each of the three subfamilies.
Differences in conclusions between previous hypotheses and the present analyses are probably due to (a) differences in taxon sampling and (b) differences in the topologies used for ACR. The MRCA of Miltogramminae subfamily is estimated as sarcosaprophagous when species of only this subfamily are analyzed (e.g., ). However, when the scope of the analysis includes other subfamilies and outgroups of other Oestroidea families, such as Oestridae, i.e., the sister-group of Sarcophagidae, the MRCA of Miltogramminae is reconstructed as a predator on invertebrates. Differences in topologies, especially regarding the earliest divergences, have an important impact on the ACR. For example, the inclusion of the genus Galopagomyia Bischof as sister to the remaining Paramacronychiinae in the supertree of Yan et al.  dramatically affects the estimated ancestral life history for the subfamilies Miltogramminae and Paramacronychiinae, giving sarcosaprophagy on invertebrates or vertebrates as the ancestral character state of the MRCA of these subfamilies. The phylogenetic position of Galopagomyia has only once been evaluated in a phylogenetic analysis including only Paramacronychiinae species and two Sarcophaginae species as outgroup . In that analysis, the sister-group relationship between Galopagomyia and the remaining Paramacronychiinae is supported by two character states, i.e., the color of the tegula with respect to the basicosta and the shape of the posterior margin of the ST5 , both of which can vary greatly across taxa of the family and do not constitute synapomorphies for the clade Galopagomyia + remaining Paramacronychiinae. For example, a black tegula contrasting with yellowish or light brown basicosta as observed in Galopagomyia is also found in many non-related genera of Sarcophaginae . Similarly, many non-related genera of Sarcophaginae share with Galopagomyia a distinctly emarginated and either broadly U-shaped or distinctly V-shaped posterior margin of the ST5. The evaluation of these characters in the broad context of the family would most likely not support the position of Galopagomyia as sister to the remaining Paramacronychiinae. Therefore, the phylogenetic position of Galopagomyia is considered uncertain. This taxon was manually added to the supertree used for ACR in Yan et al. , and therefore its phylogenetic position continues to be uncertain. Without Galopagomyia, the supertree used for ACR in Yan et al.  would have genus Agria Robineau-Desvoidy, a predator on invertebrates, as sister to the remaining Paramacronychiinae. Thus, as in Yan et al.  the MRCA of Miltogramminae is estimated as a predator or sarcosaprophage on vertebrates or invertebrates, then an ACR with an adjusted topology (without Galopagomyia) would most likely estimate the MRCA of both Miltogramminae and Paramacronychiinae as a predator on invertebrates, which would be consistent with results presented here. Another topology effect in the estimations of Yan et al.  is related to the taxa populating the early phylogenetic divergences within Sarcophaginae. The MRCA of Sarcophaginae is estimated as a sarcosaprophage on invertebrates in Yan et al.’s  supertree, as it has Tricharaea and Sarcofahrtiopsis, which are sarcosaprohages/coprophages and sarcosaprohages, respectively, as laddered sister-groups to the remaining Sarcophaginae. In contrast, the present ACR estimates the MRCA of Sarcophaginae as a predator on invertebrates based on an UCE-based phylogeny providing robust support for (Lepidodexia (Emdenimyia (Boettcheria + Tripanurga))) or clade A as sister to the remaining Sarcophaginae (Figs. 3 and 4). Not only the phylogenetic position of this clade of invertebrate predators determines the ACR estimation for the MRCA of Sarcophaginae, but the subsequent divergence of another taxon, Udamopyga Hall, which is coded as predator or sarcosaprophage on invertebrates.
An interesting transition from the predatory habits on invertebrates towards non-predatory habits in the early divergences of Sarcophaginae is observed within the clade B (Figs. 3 and 4). The MRCA of clade B is estimated as a predator on invertebrates, while the MRCA of Microcerella is estimated as a sarcosaprophage on vertebrates and the MRCA of its sister clade C is estimated as a coprophagous. Within clade C, other interesting results are related to the habitats occupied by closely related taxa and their use of available resources. Thus, the sister-group relationship between Tricharaea and Cistudinomyia tells the story of taxa living on beaches (although Tricharaea can also be found in other environments), which are adapted to predate on a specific host like turtles in the case of Cistudinomyia, and a more plastic taxon like Tricharaea using both dead invertebrates and vertebrates as well as feces to feed on. Another transition from predatory habits to sarcosaprophagy occurs in the MRCA of the clade containing genera Engelimyia, Lipoptilocnema, and Peckia, and such transition is also observed in a few species of genera Blaesoxipha and Argoravinia, and in the MRCA of (Sarcosolomonia Baranov + Sarcorohdendorfia) within the genus Sarcophaga. In contrast, there are no transitions from sarcosaprophagy or coprophagy to predatory habits, but clearly overall sarcosaprophagy and coprophagy have evolved in a nonrandom fashion in Sarcophaginae. No transitions from an exclusively sarcosaprophagous habit back to predation are supported (although exceptions could occur in Peckia and some uncertainty remains regarding Helicobia). This may indicate benefits of an adaptation to sarcosaprophagy that prevent the reversal to a predatory lifestyle. These results also indicate that the origin of sarcosaprophagy would have taken place in the context of an existing ecology of an ancestor with a predator-host established relationship. This would also support a hypothesis of a gradual transition from predation to sarcosaprophagy or coprophagy, with predator flesh flies gradually attacking not only healthy hosts, but also injured or weakened hosts, or even dead hosts. A transition from predation to sarcosaprophagy or coprophagy would imply less risk for a predator flesh fly, while the opposite could mean a sarcosaprophagous flesh fly attacking a healthy host that could eventually counterattack. More detailed species-level analyses estimating ancestral larval feeding habits within genera including ‘generalist species’ (in genera such as Argoravinia, Peckia, Sarcophaga) are necessary, which in turn could lead to slightly different conclusions at the genus level.
The observation of non-reversals from predation to sarcosaprophagy or coprophagy could also be used in future research using a different approach to model the parameter process (transitions between the different rate classes) in a hidden Markov model to assume that sarcosaprophagy and coprophagy are not lost once they evolve, which contrasts with the assumption for the models used here where all transitions among the specified number of rate classes are the same (ER) or all transition rates are allowed and are independently estimated (ARD). Such a different approach would allow the inference of a biologically relevant, but unmeasured ‘hidden’ character that could have influenced the evolution of the observed characters here.
Evolution of male terminalia traits
Male terminalia traits evolved almost equally convergently and non-convergently in multiple lineages across all three main clades of Sarcophagidae, as it has been found in other studies , although there is a slight dominance of non-convergent traits here. The six convergent traits were the shape of posterior margin of the abdominal ST5, outline of dorsal surface of cercal prong, connection between basi- and distiphallus, shape of connection between basi- and distiphallus, harpes, and capitis, while the seven non-convergent traits were the vesica, phallotrema configuration, phallotrema position with regard to phallic tube, acrophallic levers, number of styli, median process, and juxta.
Regarding convergent traits, there are degrees in convergence with character states evolving twice in the phylogeny, while other characters have multiple independent origins across the tree. For example, a straight posterior margin of the abdominal ST5 evolves twice independently in the MRCA of (Miltogramminae + Paramacronychiinae) and in the MRCA of the clade containing genera Cistudinomyia, Dexosarcophaga, Nephochaetopteryx, Oxysarcodexia, Oxyvinia, Sarcofahrtiopsis, and Tricharaea. An absent or reduced capitis also evolves twice in subfamily Miltogramminae and some subgenera of genus Peckia. Whereas characters like the outline of dorsal surface of cercal prong and the harpes have multiple independent origins and few reversals across the sarcophagid tree. The multiple origins of the harpes and possibly also the loss of capitis could be explained by a homology definition problem, as at least the harpes are difficult to delimitate and could be confused with other accessory appendages of the phallus.
An interesting case of convergence is observed for the traits related to the connection and shape of connection between basi- and distiphallus. In general, a continuous connection between basi- and distiphallus evolves twice in the family, once in the MRCA of clade (Miltogramminae + Paramacronychiinae) and another time in the MRCA of clade C (Additional file 6). A second character, which looks into the details of this connection, shows that there is a transition in the degree of sclerotization originating the continuous connection between basi- and distiphallus in clade C (Fig. 5). Thus, the MRCA of clade B has a distinct hinge between basi- and distiphallus (Fig. 5a, b), while the MRCA of genus Microcerella has an intermediate character state between a distinct hinge and a fully sclerotized connection (Fig. 5c), and the MRCA of clade C has a fully sclerotized connection (Fig. 5d–f). The intermediate character state of the genus Microcerella consists of a hinge on the dorsal side of the phallus and a sclerotized, paler, rigid and tubular area on the ventral side of the phallus between basi- and distiphallus. This sclerotized, paler, rigid and tubular ventral area between basi- and distiphallus had been described before  but not analyzed in a broader phylogenetic context of the family. Such transitions in the degree of sclerotization in the morphological evolution of Sarcophagidae or for lineages within Miltogramminae or Paramacronychiinae have not been reported before.
The diversification of the subfamily Sarcophaginae, which includes 2/3 of the diversity of the family, is marked by the concerted evolution of a set of phallic traits that were found to be non-convergent. Some of these phallic traits involve complex structures like the juxta, vesica, and a complex acrophallus (with various styli), which have only isolated reversals or losses. The loss of complex structures as irreversible over time is a concept known as Dollo’s law . Although this evolutionary principle is still commonly accepted, a number of cases where it is apparently violated have been proposed. Here I found that most of the complex phallic structures (e.g., juxta, vesica, a complex acrophallus with more than one styli) are rarely lost once they have evolved, and only the harpes seem to be the exception. The juxta originates most probably in the ancestor of all sarcophagids, and became a more complex structure separate by a hinge from the rest of the phallus in the clade A (Lepidodexia (Emdenimyia (Boettcheria + Tripanurga))) (Additional file 14), which constitutes the first branching within Sarcophaginae. The phallotrema placed in a ventral position with regard to phallic tube follows the same evolutionary pattern as the juxta (Additional file 10). The vesica and folding of the phallotrema evolve in the MRCA of Sarcophagidae (Additional files 8, 9), while the acrophallic levers evolve in the MRCA of clade C (Additional file 11). There might be some correlation between morphological characters, such as the concerted evolution of a phallotrema placed in a ventral position with regard to phallic tube and the origin of the juxta. The observed pattern in Miltogramminae suggests certain dependency between these characters, as species in this subfamily have a phallotrema placed in an apical position and absence of juxta. Similarly, concerted evolution is observed between the origin of the vesica and the folding of the phallotrema. Similar to the previous couple of characters, Miltogramminae lacks both vesica and the folding of the phallotrema. Interestingly, the vesica is particularly ornamented and complex in the genera Cistudinomyia, Dexosarcophaga, Nephochaetopteryx, Oxysarcodexia, Oxyvinia, Sarcofahrtiopsis, and Tricharaea, and it seems to be functionally related to an extrusion of the styli during mating, which is mediated by the acrophallic levers , another trait showing a non-convergent evolution that matches the evolutionary pattern of an ornamented vesica.
Regarding the relation between traits and diversification rate, the only trait having an evolutionary pattern close to the increased rate of diversification identified along the branch leading to genus Sarcophaga (Fig. 6b) is the number of styli (Fig. 6a). The MRCA of Sarcophaginae was estimated to have had three styli, which transitions to two styli in the MRCA of the large clade F that contains Sarcophaga and other genera (Fig. 2). There is no evident explanation on how the reduction in number of styli could have been related to the massive radiation within the genus Sarcophaga.
Diversification of Sarcophagidae
The majority of the flesh fly diversity is represented by the subfamily Sarcophaginae, which contains three of the largest, most species-rich genera within the Sarcophagidae: Blaesoxipha, Lepidodexia, and Sarcophaga. In the present study, one diversification rate shift was inferred, which is associated unsurprisingly with the genus Sarcophaga that is the most species-rich genus and among the geographically most widespread taxa of flesh flies. The genera Blaesoxipha and Lepidodexia are not associated with any diversification rate shift. Incomplete taxon sampling could have influenced the age estimates for Blaesoxipha and Lepidodexia in the dating analyses. In addition, Lepidodexia is the least studied among the most species-rich genera, which could suffer from an underestimated number of species. The sensitivity of BAMM analyses to the selected rate shift prior in estimations of diversification rates and rate shifts (under certain circumstances) has been questioned [53, 54]. Similarly, a tendency to overestimate diversification rates in smaller clades, which may result in a potential underestimation of rate shifts overall, has also been identified and criticized . The results presented here are most likely not biased by these sensitivity and overestimation issues affecting the diversification rate estimates, given that I recovered only one statistically significant shift for the most species-rich genus within Sacophagidae, which essentially confirms observations based on taxonomic species diversity. These results are concordant with studies showing lineages within Sarcophagidae [9, 10, 18, 23] as the dominant fast-evolving groups of Oestroidea. Furthermore, these results support a super-radiation within the genus Sarcophaga, as recent studies suggest [9,10,11, 23, 33].
The dynamics of evolutionary diversification are usually linked to ecological opportunity and the evolution of a key innovation. This means that increases in diversification rates could be the result of a lineage that evolved and diversified into previously inaccessible environmental niche space because of an ecological opportunity event (e.g., dispersal to and colonization of a new environment or extinction of a previously dominant group) or with the evolution of a key innovation, such as a novel trait (be it morphological, physiological, or genetic), resulting in rapid speciation as niche space is partitioned unencumbered by biological interactions, such as competition and predation . Evidence for morphological diversification under increased diversification rates has been recovered in some clades of insects  and mammals [58, 59] using phylogenetic methods. Similarly, there is evidence for positive shifts in diversification rate of ants [57, 60, 61] and beetles  possibly associated with increased ecological opportunity in the form of biogeographic dispersal. In the case of the increased diversification rate inferred here for the genus Sarcophaga, this could be the result of an ecological opportunity event. Previous research showed this genus evolved in the Nearctic region and experienced a rapid radiation occurring in the Nearctic region with a subsequent dispersal into the Palaearctic region . This is supported by the present results, although the diversification rate inferred here shows that the rate increase occurred once some Sarcophaga lineages had dispersed into the Palaearctic region and not before in the Nearctic region. Thus, the increased diversification rate inferred here for the large part of the genus Sarcophaga could be the result of a lineage that evolved and diversified into a previously inaccessible environmental niche space. Thus, the observed resemblance between the evolutionary pattern of traits such as the number of styli and pattern of diversification rate might not explain the rapid radiation of the non-Nearctic lineages of Sarcophaga. Even if the evolutionary pattern of this trait closely resemble the diversification pattern of Sarcophaga, the reduction in the number of styli is observed for the large clade F, which includes Sarcophaga but also many other genera (i.e., Chrysagria, Helicobia, Peckiamyia, Retrocitomyia, Villegasia, Peckia, Engelimyia, and Lipoptilocnema) that did not show increases in diversification rates and did not expand their biogeographic distribution beyond the Neotropical region. This hypothesis of diversification associated with biogeographic dispersal and ecological opportunity could be further tested with extended species-level phylogenetic, biogeographic, and biological data, while hypotheses regarding the evolutionary significance and influence of innovative features such as morphological traits in the diversification of sarcophagids will possibly present more difficult challenges.
The phylogenomic approach combining taxon-specific, protein encoding, UCE probes with a large ingroup and outgroup sampling obtained a well-supported phylogeny for Sarcophagidae at the subfamily, genus, and species level. I was able to firmly place a group of predatory flies including the genera Lepidodexia and Boettcheria as sister-group to the rest of Sarcophaginae, which contrasts with previous phylogenetic hypotheses for sarcophagids. Similarly, genera Oxysarcodexia, Ravinia, and Tricharaea, long considered archaic and early-branching coprophagous and sarcosaprophagous lineages, were found nested well within the Sarcophaginae as sister-group to the sarcosaprophagous Microcerella. The ACR estimated predation on invertebrates as the ancestral and dominant strategy throughout the early evolution of flesh flies independent of the type of ancestral character reconstruction analysis performed (Mesquite- or corHMM-based). This strategy is also supported for the most recent common ancestor of each of the three subfamilies. These estimations also suggest that sarcosaprophagy, coprophagy, and kleptoparasitism evolved from predation. Thus, several transitions from predation to sarcosaprophagy and coprophagy were estimated across the phylogenetic tree of sarcophagids, in contrast with almost no transitions from sarcosaprophagy or coprophagy to predatory habits. The evolution of morphological traits seems not related with the change in feeding habit, although there might be a correlation between morphological characters only, such as the concerted evolution of the position of the phallotrema and the origin of the juxta, or the origin of the vesica and the folding of the phallotrema. The only diversification rate shift inferred in the evolution of Sarcophagidae is associated unsurprisingly with the richest and geographically most widespread taxon, the genus Sarcophaga, which could be the result of an increased ecological opportunity in the form of biogeographic dispersal and colonization of the Palaearctic region. As suggested in early studies on Sarcophaginae, a correlation could be present between feeding habits and physiological rather than morphological change , which is a hypothesis that should be tested in future research. Future studies using UCEs should include and phylogenetically place enigmatic taxa such as genera Carinoclypeus Dodge, Galopagomyia, Sarcodexiopsis Townsend, Sarothromyiops Townsend, Sinopiella Lopes & Tibana, Tapacura Tibana & Lopes, and Thomazomyia Lopes. By reconstructing a robust phylogeny and highlighting patterns of life histories and morphological evolution, this study has established the framework for further physiological and comparative genomic work between predatory, sarcosaprophage, coprophage, and kleptoparasitic lineages, which could also have significant implications for the evolution of diverse life histories in other Diptera.
Taxon sampling included 110 ingroup Sarcophagidae species of the subfamilies Miltogramminae (4 spp.), Paramacronychiinae (4 spp.), Sarcophaginae (102 spp.), and 31 outgroup species representing the families Anthomyiidae (2 spp.), Calliphoridae (7 spp.), Fanniidae (2 spp.), Mesembrinellidae (2 spp.), Oestridae (4 spp.), Polleniidae (1 sp.), Rhiniidae (4 spp.), Rhinophoridae (3 spp.), Scathophagidae (1 sp.), and Tachinidae (5 spp.). The selected outgroups represent fly families that could question the monophyly of Sarcophagidae and that are relevant representatives to evaluate the evolution of life history within the superfamily Oestroidea. Even tough the monophyly of Sarcophagidae has been studied and tested before [24, 25, 33], the present study includes a large outgroup in order to produce a robust phylogeny for the ancestral character state reconstruction. Thus, the present study includes the closest relatives of Sarcophagidae, i.e., fly families of the superfamily Oestroidea (Calliphoridae, Mesembrinellidae, Oestridae, Polleniidae, Rhiniidae, Rhinophoridae, and Tachinidae, except Mystacinobiidae and Ulurumyiidae) and more distant relatives belonging to the Muscoidea grade (i.e., Anthomyiidae, Fanniidae, Scathophagidae). All specimens included in this study were collected in accordance with local regulations and all necessary permits were obtained. Voucher specimens have been deposited at the Entomological Collection of Alicante University (CEUA), Museum für Naturkunde (MFN), Musée royal de l’Afrique centrale (RMCA), National Museum of Natural History (USNM), Natural History Museum of Denmark (NHMD), and Entomological Collection of Wright State University (JOSC) (see Additional file 1 for specimen identifiers and collection data).
UCE data collection
Genomic DNA was obtained from DNA aliquots, pinned museum specimens, specimens preserved in 96% Ethanol, and in Liquid Nitrogen. I used 32 existing DNA aliquots from previous molecular studies [10, 11, 33, 47], which were stored in a − 20 °C freezer. Genomic DNA was extracted from 24 pinned specimens, 66 specimens preserved in 96% Ethanol, and 19 specimens collected and placed directly in empty vials stored in Liquid Nitrogen in the National Museum of Natural History (USNM) Biorepository, as indicated in Additional file 1, where specimen identity, preservation method, targeted tissue for extraction (i.e., thorax, abdomen, legs or the whole body excluding terminalia), collection data and corresponding repositories at natural history museums of all specimens is provided. Dust, pollen, and other forms of accumulated debris on pinned specimens were removed using sterilized forceps and a soft paintbrush. DNA was non-destructively extracted from the thorax of pinned specimens, while it was destructively extracted from specimens preserved in 96% ethanol and liquid nitrogen by grinding the tissue with a sterile pestle. DNA extractions used a DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA, USA) and followed the manufacturer’s protocol, but to maximize DNA yield the Proteinase K digestion ran for 48 h at 56 °C and DNA was eluted twice in 50 µL (total volume 100 µL). To estimate size of the genomic DNA, 10 µL of each extract were run for 40 min at 100 volts on 1.5% agarose SB (sodium borate) gels.
Library preparation, target enrichment, and sequencing of UCEs
Extracted genomic DNA was quantified using a Qubit fluorometer (High sensitivity kit, Life Technologies, Inc.). DNA (0.3–590 ng, 107.1 ng mean) was sheared to a target size of approximately 500–600 bp by sonication (Q800, Qsonica LLC.), depending on prior degradation and fragmentation of DNA. This sheared DNA was used as input for library preparation following a protocol for UCEs by Faircloth et al.  and detailed in Blaimer et al. . For adapter ligation, I used Tru-Seq-style adapters  and PCR amplified 50% of the resulting library volume (15 µL) with a reaction mix of 25 µL HiFi HotStart polymerase (Kapa Biosystems), 2.5 µL each of Illumina TruSeq-style i5 and i7 primers (5 µM each), and 5 µL double-distilled water (ddH20). I used the following thermal protocol: 98 °C for 45 s; 13 cycles of 98 °C for 15 s, 65 °C for 30 s, 72 °C for 60 s, and final extension at 72 °C for 5 m. After rehydrating (in 23 µL pH 8 Elution Buffer) and purifying reactions using 1.0× speedbeads, 8 libraries were combined at equimolar ratios into enrichment pools with final concentrations of 85.4–189 ng/µL. Each pool was enriched using a set of 5,117 custom-designed probes targeting 2,581 UCE loci in Calyptratae flies  using the MYcroarray MYBaits kit , except I used a 0.1x concentration of the standard MYBaits concentration, and added 0.7 µL of 500 µL custom blocking oligos designed against the custom sequence tags. The pool hybridization reaction ran for 24 h at 65 °C. Subsequently, I bound all pools to streptavidin beads (MyOne C1, Life Technologies) and washed bound libraries according to a standard target enrichment protocol . Post-enrichment amplification was performed on beads with the KAPA Hifi HotStart ReadyMix using the following thermal profile: 98 °C for 45 s; 18 cycles of 98 °C for 15 s, 60 °C for 30 s, 72 °C for 60 s; and a final extension of 72 °C for 5 m. Post-enrichment libraries were purified using 1.0x speedbeads (Sera-mag, GE Healthcare) and rehydrated the enriched pools in 22 µL TLE. Post-enrichment library concentration was quantified via qPCR using a SYBR® FAST qPCR kit (Kapa Biosystems) on a ViiA™ 7 (Life Technologies). Based on the size-adjusted concentrations estimated by qPCR, I pooled libraries at equimolar concentrations and size-selected for 250–800 with a BluePippin (SageScience) (1.5% agarose, 250 bp–1.5 kb), and the pool-of-pools was quality checked on an Agilent 2200 TapeStation. The pooled libraries were sequenced using two lanes of a 125-bp paired-end Illumina HiSeq 2500 run (University of Utah Genomics Core Facility).
Processing and alignment of UCE data
Illumiprocessor , based on the package Trimmomatic , was used to trim the demultiplexed FASTQ data output for adapter contamination and low-quality bases. Cleaned reads were assembled using Trinity . All further data processing relied on the PHYLUCE package [69, 70] with Python scripts designed by the Smithsonian Institution Bioinformatics Group (available at www.github.com/SmithsonianWorkshops/Targeted_Enrichment/blob/master/phyluce.md). Summary statistics were computed on the data using the phyluce_assembly_get_fastq_stats.py script. Average sequencing coverage and contig length across assembled contigs were calculated using the phyluce_assembly_get_trinity_coverage.py script. To identify contigs representing enriched UCE loci from each species, species-specific contig assemblies were aligned to a FASTA file of all enrichment baits (min_coverage = 70, min_identity = 80), and sequence coverage statistics (avg, min, max) for contigs containing UCE loci were calculated. I created FASTA files for each UCE locus containing sequence data for taxa present at that particular locus and aligned these using MAFFT  (min-length = 20, no-trim). Alignments were trimmed using Gblocks  with relaxed settings (b1 = 0.5, b2 = 0.5, b3 = 12, b4 = 7). The alignment design used two different configurations. First, in an exploratory stage, concatenated UCE alignments having varying loci occupancy (0.1–0.9, nine datasets) were submitted to phylogenetic reconstruction and analyzed as nucleotides. Second, UCE alignments having varying loci occupancy (0.1, 0.3, 0.6, 0.75, four datasets) were submitted to phylogenetic reconstruction and analyzed both as nucleotides and amino acids. All of these datasets (Table 1) were designed to evaluate the relative contribution of varying amounts of UCE loci as nucleotides or amino acids to the construction of the phylogenetic tree.
First, datasets having varying loci occupancy and coded as nucleotides (analyses 1–9 in Table 1) were analyzed as concatenated datasets under Maximum Likelihood (ML). ML best tree (model GTRGAMMAI) and bootstrap searches (N = 100) of the nine concatenated datasets were conducted in RAxML v8.2.7 . As the parameters Gamma and the proportion of invariable sites cannot be optimized independently from each other, I also analyzed the nine concatenated datasets using the GTRCAT model, which showed a reduction in conflict between resulting topologies, therefore these results are preferred. Subsequently, each of the eight datasets (analyses 10–17 in Table 1) having varying loci occupancy and coded as nucleotides and amino acids was analyzed using gene trees in a multi-coalescent species tree analysis with ASTRAL-III . Data were partitioned by individual UCE loci using the Sliding-Window Site Characteristics approach and site characteristics such as entropy implemented in the SWSC-EN algorithm, which generates partitions that account for heterogeneity in rates and patterns of molecular evolution within each UCE . A partitioning scheme from the by-locus character sets was selected with PartitionFinder2 . Then, I sequentially ran an ML analysis for the best tree and 1000 replicates of ultrafast bootstrap on each locus for the gene tree estimations using IQ-TREE . A multi-coalescent species tree analysis was carried out in ASTRAL-III using gene trees (one tree search per gene) estimated by 100 ML searches conducted in RAxML. Statistical supports by ASTRAL-III are local posterior probabilities (LPP), which are branch support values that measure the support for a quadripartition, not a bipartition.
All of the above phylogenomic analyses were performed on the Smithsonian Institution High Performance Cluster Hydra (SI/HPC) using Python scripts (designed by the Smithsonian Institution Bioinformatics Group, Michael Lloyd, and some modified by Bonnie Blaimer and myself) (available at www.github.com/SmithsonianWorkshops/Targeted_Enrichment/blob/master/phyluce.md). Tips of final trees were renamed using a Perl script (designed by Michael Lloyd and available at www.github.com/MikeWLloyd/Tree-Tip-Replacer).
Diversification rates of Sarcophagidae were estimated using the coalescent species-tree obtained from the dataset having loci occupancy 0.10 (includes 2018 UCE loci) and analyzed as nucleotides, which produced a strongly supported topology consistent with most of trees obtained from the remaining analyses. All of the outgroups were trimmed using the R package phytools v0.7-70 . The coalescent species tree was transformed into a chronogram using correlated rates with makeChronosCalib in the R package ape v5.4-1 , setting the age of the most recent common ancestor of Sarcophagidae at 23 Ma (11.4–38.0) . This age is uncertain due to the lack of fossils for sarcophagids and other groups of Oestroidea, but it is conservative with respect to other age estimates.
Diversification rate estimation
An evaluation on whether shifts in diversification rates occurred over time in the evolution of Sarcophagidae was conducted using BAMM v2.5 [80, 81] and the associated R package BAMMtools v2.1.5 . Incomplete sampling was accounted by using clade-specific sampling probabilities. To calculate sampling probabilities, a species-richness matrix was assembled. This matrix lists each tip (i.e., a species) of the phylogeny, which is assigned to a clade (i.e., its respective genus). Thus, the sampling fraction for each clade is calculated as the number of sampled species divided by number of described species in that clade (Additional file 15). The chronogram obtained above plus the set of sampling fractions were used to perform BAMM analyses, which followed the guidelines in the BAMM documentation (http://bamm-project.org/). The function setBAMMpriors within BAMMtools was used to obtain appropriate priors for speciation-extinction analyses, and the expected number of shifts was left at the default value (= 1). The Bayesian analysis included four MCMC chains with a length of 25 million generations, sampling every 10,000 generations, with discarded burn-in of 10%. Convergence was visualized using plots of log likelihoods of all sampled generations and all model parameters over time were examined focusing especially on effective sample size (ESS) values with the R package coda v0.19–4 . Diversification rates were analyzed and visualized with various functions in BAMMtools, as follows: computeBayesFactors compared the evidence for models with at least one diversification shift to the evidence for the null model (zero diversification shifts) and identified the best-supported model of rate shifts; credibleShiftSet computed 95% credible set of distinct shift configurations that account for 95% of the probability of the data; getBestShiftConfiguration extracted the rate shift configuration with the maximum posterior probability (MAP); plot.bammdata plotted a ‘phylorate’ graph showing mean, model-averaged diversification rates along branches of the phylogeny.
Evolution of sarcophagid life histories and male terminalia traits
Since I was interested in the evolutionary trajectory of life histories and male terminalia in sarcophagids, I inferred ancestral character states using the dated phylogeny. Biological data were compiled mostly from the literature [3, 4, 51, 172,173,174,175,176, 83,84,85,86,87,88,89,90,91, 5, 92,93,94,95,96,97,98,99,100,101, 6, 102,103,104,105,106,107,108,109,110,111, 13, 112,113,114,115,116,117,118,119,120,121, 14, 122,123,124,125,126,127,128,129,130,131, 18, 132,133,134,135,136,137,138,139,140,141, 20, 142,143,144,145,146,147,148,149,150,151, 47, 152,153,154,155,156,157,158,159,160,161, 49, 162,163,164,165,166,167,168,169,170,171] and personal collecting experience. Definition of life-history characters followed Yan et al. , where the characters larval food resource and larval feeding habits were selected to represent the breeding and feeding strategies of sarcophagid and Calyptratae larvae. The character larval food resource had five character states as follows: 0) invertebrates, (1) vertebrates, (2) feces, (3) plants, and (4) unknown. In the present study, the character larval feeding habits had six character states as follows: 0) kleptoparasitism, (1) sarcosaprophagy, (2) predation, (3) coprophagy, (4) phytophagy, (5) unknown. Yan et al.  did not include phytophagy as a character state of larval feeding habits. Following Yan et al. , predation was defined to include parasitism and parasitoidism; sarcosaprophagy (feeding on carrion) and coprophagy (feeding on faeces) were defined as distinct character states because their associated volatile profiles are different ; and kleptoparasitism was defined as a distinct character state due to the combined behavioural, chemical, and morphological especializations needed for the gravid females and the first instar larvae to survive interactions with their hosts .
Thirteen morphological characters of the male terminalia were selected from the dataset of Buenaventura & Pape . The selected characters are: (1) Shape of posterior margin of the abdominal sternite 5 (ST5), (2) outline of dorsal surface of cercal prong in lateral view, (3) connection between basi- and distiphallus, (4) shape of connection between basi- and distiphallus, (5) harpes, (6) presence of vesica, (7) phallotrema configuration, (8) phallotrema position with regard to phallic tube, (9) presence of acrophallic levers, (10) number of styli, (11) presence of capitis, (12) median process, and (13) presence of juxta.
Using the biological and morphological characters described above, two datasets were produced. The first dataset includes all taxa sampled (ingroup and outgroup) for which only biological characters are scored (Additional file 16). The second dataset includes only ingroup species for which all of the 15 biological and morphological characters are scored (Additional file 17). The character state phytophagy of larval feeding habits was only included in the first dataset as this character state was sampled in the outgroup species and it is excluded from the second dataset as it was not sampled in the ingroup species. Morphological characters of the male terminalia were scored for the ingroup only and included only in the second dataset, as the homology of these characters needs to be further studied in comparative analyses across the whole Calyptratae clade. Species of flesh flies generally have marked preferences for carrion or faeces, but for some species, there is no certainty about their preferences. Adults of some species have been observed as visitors of both carrion and feces, but there is no information about their larval feeding strategies. In such cases, species were scored as polymorphic when feeding data supporting the polymorphism was available. When data was not available (unknown), species were scored with ‘?’ for the first dataset and with ‘n/a’ for the second dataset.
For ACR, I used the trace character option with Maximum parsimony under unordered states assumption in the software Mesquite version 3.04 (build 725) [179, 180] to trace the life-history characters both in outgroup and ingroup taxa using the first dataset. In parallel, also for ACR, I used the rayDISC function in the package corHMM in R v4.0.0 (https://www.R-project.org/) which can analyze multivariate traits and the time-calibrated phylogeny estimated while employing the median range calibration on the root node. The corHMM-based ACR was performed on each biological and morphological character independently using the second dataset. Thus, each ACR analysis used two models available in corHMM, ‘equal rates’ (ER) and ‘all rates different’ (ARD). All analyses were carried out twice with the same settings to ascertain robustness of reconstructions. The fit of these models (i.e., ER and ARD) for the ACR analysis on each character was compared using the resulting lnL scores and corrected Akaike information criterion (AICc, test corrected for small sample sizes) values (Additional file 3). Results are presented for ACR with models significantly better considering AICc estimates. When either ER or ARD are not significantly better, then results are presented for the model with the lowest AICc estimate.
Availability of data and materials
The UCE data supporting the conclusions of this article are available in the Museum für Naturkunde repository under accession https://doi.org/10.7479/t80n-ta02 . The datasets (biological and morphological data) supporting the conclusions of this article are included within the article and its additional files.
Ancestral character state reconstruction
Anchored hybrid enrichment
Corrected Akaike information criterion
All rates different
Entomological Collection of Alicante University
Effective sample size
Entomological Collection of Wright State University
Local posterior probability
Million years ago
Maximum posterior probability
Markov Chain Monte Carlo
Most recent common ancestor
Next generation sequencing
Natural History Museum of Denmark
National Museum of Natural History
Musée royal de l’Afrique central
Smithsonian Institution High Performance Cluster Hydra
Abdominal sternite 5
Sliding-Window Site Characteristics Entropy
United States National Museum of Natural History
Museum für Naturkunde
Scholl JP, Wiens JJ. Diversification rates and species richness across the Tree of Life. Proc R Soc B Biol Sci. 2016;283:1838.
Cerretti P, Stireman JO, Pape T, O’Hara JE, Marinho MAT, Rognes K, et al. First fossil of an oestroid fly (Diptera: Calyptratae: Oestroidea) and the dating of oestroid divergences. PLoS ONE. 2017;12(8):1–24. https://doi.org/10.1371/journal.pone.0182101.
Lopes HS. Collecting and rearing Sarcophagid flies (Diptera) in Brazil during forty years. An Acad Bras Cienc. 1973;45:279–91.
Davis E, Turner W. Biology, distribution and abundance of flesh flies (Diptera: Sarcophagidae) of the Wallowa- Whitman National Forest in northeastern Oregon. Melanderia. 1978;30:111–59.
Pape T, Dahlem G. Sarcophagidae. In: Zumbado M, Brown BV, Borkent A, Cumming JM, Wood DM, Woodley NE, editors. Manual of Central American Diptera, Vol 2. Ottawa: NRC Research Press; 2010. p. 1313–35.
Pape T. Catalogue of the Sarcophagidae of the World (Insecta: Diptera). Mem Entomol Int. 1996;8:1–558.
Buenaventura E, Pape T. Phylogeny, evolution and male terminalia functionality of Sarcophaginae (Diptera: Sarcophagidae). Zool J Linn Soc. 2017;183(4):808–906. https://doi.org/10.1093/zoolinnean/zlx070.
Yan L, Buenaventura E, Pape T, Kutty SN, Bayless KM, Zhang D. A phylotranscriptomic framework for flesh fly evolution (Diptera, Calyptratae, Sarcophagidae). Cladistics. 2020;1:1–33.
Buenaventura E, Szpila K, Cassel BK, Wiegmann BM, Pape T. Anchored hybrid enrichment challenges the traditional classification of flesh flies (Diptera: Sarcophagidae). Syst Entomol. 2019;45(2):281–301.
Buenaventura E, Pape T. Multilocus and multiregional phylogeny reconstruction of the genus Sarcophaga (Diptera, Sarcophagidae). Mol Phylogenet Evol. 2017;107:619–29.
Buenaventura E, Whitmore D, Pape T. Molecular phylogeny of the hyperdiverse genus Sarcophaga (Diptera: Sarcophagidae), and comparison between algorithms for identification of rogue taxa. Cladistics. 2017;33(2):109–33.
Buenaventura E, Camacho G, Garcia A, Wolff M. Sarcophagidae (Diptera) de importancia forense en Colombia: claves taxonómicas, notas sobre su biología y distribución. Rev Colomb Entomol. 2009;35(2):189–96.
Byrd JH, Castner JL. Forensic Entomology: The Utility of Arthropods in Legal Investigations. New York: CRC Press; 2011. 705 p.
Wells JD, Pape T, Sperling FA. DNA-based identification and molecular systematics of forensically important Sarcophagidae (Diptera). J Forensic Sci. 2001;46(5):1098–102.
Meiklejohn KA, Wallman JF, Dowton M. DNA barcoding identifies all immature life stages of a forensically important flesh fly (Diptera: Sarcophagidae). J Forensic Sci. 2013 Jan;58(1):1–7.
Eggleton P, Belshaw R. Insect parasitoids: and evolutionary overview. Philos Trans R Soc B Biol Sci. 1992;337:1–20.
Feener Jr. DH, Brown BV. Diptera as parasitoids. Annu Rev Entomol. 1997;42:73–97.
Piwczyński M, Pape T, Deja-Sikora E, Sikora M, Akbarzadeh K, Szpila K. Molecular phylogeny of Miltogramminae (Diptera: Sarcophagidae): implications for classification, systematics and evolution of larval feeding strategies. Mol Phylogenet Evol. 2017;116:49–60. http://linkinghub.elsevier.com/retrieve/pii/S1055790317302634.
Denno RF, Cothran WR. Competitive interactions and ecological strategies of sarcophagid and calliphorid flies inhabiting rabbit carrion. Ann Entomol Soc Am. 1976;69(1):5. http://www.ingentaconnect.com/content/esa/aesa/1976/00000069/00000001/art00026.
Bänziger H, Pape T. Flowers, faeces and cadavers: Natural feeding and laying habits of flesh flies in Thailand (Diptera: Sarcophagidae, Sarcophaga spp.). J Nat Hist. 2004;38:1677–94. http://www.tandfonline.com/doi/abs/10.1080/0022293031000156303.
Raguso RA. Don’t forget the flies: dipteran diversity and its consequences for floral ecology and evolution. Appl Entomol Zool. 2020;55(1):1–7. https://doi.org/10.1007/s13355-020-00668-9.
Stamper T, Dahlem G, Cookman C, Debry R. Phylogenetic relationships of flesh flies in the subfamily Sarcophaginae based on three mtDNA fragments (Diptera: Sarcophagidae). Syst Entomol. 2012;38(1):35–44.
Piwczyński M, Szpila K, Grzywacz A, Pape T. A large-scale molecular phylogeny of flesh flies (Diptera: Sarcophagidae). Syst Entomol. 2014 Oct 17;39(4):783–99.
Kutty SN, Meusemann K, Bayless KM, Marinho MAT, Pont AC, Zhou X, et al. Phylogenomic analysis of Calyptratae: resolving the phylogenetic relationships within a major radiation of Diptera. Cladistics. 2019;35(6):605–22.
Kutty SN, Pape T, Wiegmann BM, Meier R. Molecular phylogeny of the Calyptratae (Diptera: Cyclorrhapha) with an emphasis on the superfamily Oestroidea and the position of Mystacinobiidae and McAlpine’s fly. Vol. 35, Systematic Entomology. 2010. p. 614–35.
Young AD, Lemmon AR, Skevington JH, Mengual X, Ståhls G, Reemer M, et al. Anchored enrichment dataset for true flies (order Diptera) reveals insights into the phylogeny of flower flies (family Syrphidae). BMC Evol Biol [Internet]. 2016;16:143. https://doi.org/10.1186/s12862-016-0714-0.
Lemmon AR, Emme S a., Lemmon EM. Anchored hybrid enrichment for massively high-throughput phylogenomics. Syst Biol. 2012;61(5):727–44.
Crawford NG, Parham JF, Sellas AB, Faircloth BC, Glenn TC, Papenfuss TJ, et al. A phylogenomic analysis of turtles. Mol Phylogenet Evol [Internet]. 2015;83:250–7. https://doi.org/10.1016/j.ympev.2014.10.021.
Sun K, Meiklejohn KA, Faircloth BC, Glenn TC, Braun EL, Kimball RT. The evolution of peafowl and other taxa with ocelli (eyespots): a phylogenomic approach. Proc Biol Sci. 2014;281(1790):20140823–20140823.
Faircloth BC, Sorenson L, Santini F, Alfaro ME. A phylogenomic perspective on the radiation of ray-finned fishes based upon targeted sequencing of ultraconserved elements (UCEs). PLoS ONE. 2013;8:6.
Wood HM, González VL, Lloyd M, Coddington J, Scharff N. Next-generation museum genomics: Phylogenetic relationships among palpimanoid spiders using sequence capture techniques (Araneae: Palpimanoidea). Mol Phylogenet Evol. 2018;127:907–18.
Blaimer BB, LaPolla JS, Branstetter MG, Lloyd MW, Brady SG. Phylogenomics, biogeography and diversification of obligate mealybug-tending ants in the genus Acropyga. Mol Phylogenet Evol. 2016;102:20–9.
Buenaventura E, Lloyd MW, Perilla López JM, González VL, Thomas-Cabianca A, Dikow T. Protein-encoding ultraconserved elements provide a new phylogenomic perspective of Oestroidea flies (Diptera: Calyptratae). Syst Entomol. 2021;46(1):5–27. https://doi.org/10.1111/syen.12443.
Blaimer BB, Lloyd MW, Guillory WX, Brady SG. Sequence capture and phylogenetic utility of genomic ultraconserved elements obtained from pinned insect specimens. PLoS ONE. 2016;11(8):1–20.
McCormack JE, Tsai WLE, Faircloth BC. Sequence capture of ultraconserved elements from bird museum specimens. Mol Ecol Resour. 2016;16(5):1189–203.
Baca SM, Alexander A, Gustafson GT, Short AEZ. Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga.’ Syst Entomol. 2017;42(4):786–95.
Derkarabetian S, Benavides LR, Giribet G. Sequence capture phylogenomics of historical ethanol-preserved museum specimens: Unlocking the rest of the vault. Mol Ecol Resour. 2019;19(6):1531–44. https://doi.org/10.1111/1755-0998.13072.
Wiegmann BM, Trautwein MD, Winkler IS, Barr NB, Kim J-W, Lambkin CL, et al. Episodic radiations in the fly tree of life. Proc Natl Acad Sci USA. 2011;108(14):5690–5.
Marinho MAT, Junqueira ACM, Paulo DF, Esposito MC, Villet MH, Azeredo-Espin AML. Molecular phylogenetics of Oestroidea (Diptera: Calyptratae) with emphasis on Calliphoridae: Insights into the inter-familial relationships and additional evidence for paraphyly among blowflies. Vol. 65, Molecular Phylogenetics and Evolution. 2012. p. 840–54.
Pape T. Phylogeny of the Tachinidae family-group (Diptera: Calyptratae). Tijdschr voor Entomol. 1992;135:43–86.
Zhang D, Yan L, Zhang M, Chu H, Cao J, Li K, et al. Phylogenetic inference of Calyptrates, with the first mitogenomes for Gasterophilinae (Diptera: Oestridae) and Paramacronychiinae (Diptera: Sarcophagidae). Int J Biol Sci. 2016;12(5):489–504.
Yan L, Zhang M, Gao Y, Pape T, Zhang D. First mitogenome for the subfamily Miltogramminae (Diptera: Sarcophagidae) and its phylogenetic implications. Eur J Entomol. 2017;114:422–9.
Song Z, Wang X, Liang G. Molecular evolution and phylogenetic utility of the internal transcribed spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera). J Mol Evol. 2008 Nov;67(5):448–64.
Pape T. A new genus of Paramacronychiinae (Diptera: Sarcophagidae), argued from a genus-level cladistic analysis. Syst Entomol. 1998;23:187–200.
Giroux M, Pape T, Wheeler TA. Towards a phylogeny of the flesh flies (Diptera: Sarcophagidae): morphology and phylogenetic implications of the acrophallus in the subfamily Sarcophaginae. Zool J Linn Soc. 2010;158:740–78.
Buenaventura E, Pape T. Phylogeny of the Peckia-genus group: evolution of male genitalia in the major necrophagous guild of Neotropical flesh flies (Diptera: Sarcophagidae). Org Divers Evol. 2015;15:301–31.
Zhang M, Buenaventura E, Pape T, Zhang D. Molecules and morphology unite Sarcophaga (Stackelbergeola) Rohdendorf and S. (Rohdendorfisca) Grunin within megadiverse Sarcophaga Meigen (sensu lato) (Diptera, Sarcophagidae). Invertebr Syst. 2016;30:463–78.
Roback SS. The evolution and taxonomy of the Sarcophaginae (Diptera, Sarcophagidae). B III, Illinois Biol Monogr. 1954;23:1–181.
Povolný D, Verves YG. The Flesh-Flies of Central Europe (Insecta, Diptera, Sarcophagidae). Spixiana, Zeitschrift für Zool. 1997;24:1–263. http://www.pfeil-verlag.de/04biol/e2_24.html.
Verves YG. The origin of inquilinism of miltogrammatin flies (Miltogrammatinae Sarcophagidae, Diptera). Visnyk kyyivs Univ (Ser Biol). 1976;18:106–108.
Mulieri PR, Mariluis JC, Patitucci LD, Olea MS. The Sarcophaginae (Diptera: Sarcophagidae) of Southern South America. I. The species of Microcerella Macquart from the Patagonian Region. Zootaxa. 2015;3933(1):1–88.
Dollo L. The laws of evolution. Bull Soc Bel Geol Paleontol. 1893;7:164–166.
Moore BR, Höhna S, May MR, Rannala B, Huelsenbeck JP. Critically evaluating the theory and performance of Bayesian analysis of macroevolutionary mixtures. Proc Natl Acad Sci USA. 2016;113:9569–74.
Rabosky DL, Mitchell JS, Chang J. Is BAMM flawed? Theoretical and practical concerns in the analysis of multi-rate diversification models. Syst Biol. 2017;66:477–498.
Meyer ALS, Wiens JJ. Estimating diversification rates for higher taxa: BAMM can give problematic estimates of rates and rate shifts. Evolution. 2018;72:39–53.
Yoder JB, Clancey E, Des Roches S, Eastman JM, Gentry L, Godsoe W, et al. Ecological opportunity and the origin of adaptive radiations. J Evol Biol. 2010;23(8):1581–96.
Price SL, Etienne RS, Powell S. Tightly congruent bursts of lineage and phenotypic diversification identified in a continental ant radiation. Evolution. 2016;70(4):903–12.
Cooper N, Purvis A. Body size evolution in mammals: complexity in tempo and mode. Am Nat. 2010;175(6):727–38.
Slater GJ, Price SA, Santini F, Alfaro ME. Diversity versus disparity and the radiation of modern cetaceans. Proc R Soc B Biol Sci. 2010;277(1697):3097–3104.
Blaimer BB, Ward PS, Schultz TR, Fisher BL, Brady SG. Paleotropical diversification dominates the evolution of the hyperdiverse ant tribe Crematogastrini (Hymenoptera: Formicidae). Insect Syst Divers. 2018;2:5.
Price SL, Powell S, Kronauer DJC, Tran LAP, Pierce NE, Wayne RK. Renewed diversification is associated with new ecological opportunity in the Neotropical turtle ants. J Evol Biol. 2014;27(2):242–58.
Toussaint EFA, Hendrich L, Shaverdo H, Balke M. Mosaic patterns of diversification dynamics following the colonization of Melanesian islands. Sci Rep. 2015;5:1–11.
Faircloth BC, Branstetter MG, White ND, Brady SG. Target enrichment of ultraconserved elements from arthropods provides a genomic perspective on relationships among Hymenoptera. Mol Ecol Resour. 2015;15:489–501.
Faircloth BC, Glenn TC. Not all sequence tags are created equal: designing and validating sequence identification tags robust to indels. PLoS ONE. 2012;7:8.
Blumenstiel B, Cibulskis K, Fisher S, DeFelice M, Barry A, Fennell T, et al. Targeted exon sequencing by in-solution hybrid selection. Curr Protoc Hum Genet. 2010;66(1):18.4.1–18.4.24.
Faircloth BC. Illumiprocessor: a trimmomatic wrapper for parallel adapter and quality trimming. 2013;25:78
Bolger AM, Usadel B, Lohse M. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–652.
Faircloth BC. PHYLUCE is a software package for the analysis of conserved genomic loci. Bioinformatics. 2016;32:786–788.
Faircloth BC, McCormack JE, Crawford NG, Harvey MG, Brumfield RT, Glenn TC. Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales. Syst Biol. 2012;61(5):717–26.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17(4):540–52.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:15312–23. http://bioinformatics.oxfordjournals.org/content/early/2014/01/21/bioinformatics.btu033.abstract.
Zhang C, Rabiee M, Sayyari E, Mirarab S. ASTRAL-III: Polynomial time species tree reconstruction from partially resolved gene trees. BMC Bioinformatics. 2018;19(Suppl 6):15–30. https://doi.org/10.1186/s12859-018-2129-y.
Tagliacollo VA, Lanfear R. Estimating improved partitioning schemes for ultraconserved elements. Mol Biol Evol. 2018;35:1798–1811.
Lanfear R, Frandsen P, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2017;34:772–773.
Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74.
Revell L. phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Paradis E, Schliep K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35:526–8.
Rabosky DL, Santini F, Eastman J, Smith SA, Sidlauskas B, Chang J, et al. Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation. Nat Commun. 2013;4:1958.
Rabosky D, Grundler M, Anderson C, Title P, Shi J, Brown J, et al. BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees. Methods Ecol Evol. 2014;5:701–7.
Plummer M, Best N, Cowles K, Vines K. CODA: convergence diagnosis and output analysis for MCMC. R News. 2006;6:7–11.
Spofford MG, Kurczewski FE. Nearctic species of Miltogrammini (Diptera: Sarcophagidae) associated with species of Aculeata (Hymenoptera: Vespoidea, Pompiloidea, Sphecoidea, Apoidea). J Kansas Entomol Soc. 1989;62:254–267.
Nandi BC. The fauna of India and the adjacent countries. Calcutta: Zoological Survey of India; 2002. 1–608 p.
Szpila K, Akbarzadeh K, Pape T. First description of the first instar larva of Sphecapatoclea and Sphecapatodes (Diptera: Sarcophagidae). Zool Anz. 2017;266:129–35.
Broderick A, Hancock E. Insect infestation of Mediterranean marine turtle eggs. Herpetol Rev. 1997;28(4):190–1.
Szpila K. The first instar of European Miltogramminae (Diptera, Sarcophagidae). Torun: Wydawnictwo Naukowe UMK; 2010.
Xu WT, Zhang D, Pape T. Biology of Eumacronychia Townsend, with a redescription of E. persolla Reinhard, 1965 (Diptera, Sarcophagidae). Zookeys. 2018;2018(783):55–65.
Tereshkin A. Entomophagous parasites of nun moth (Lymantria monacha L.) in Byelorussia (in the period of outbreak of mass reproduction in 1976–1984). Tereshkin, A.M., Minsk-Nowosibirsk Minsk-Novosibirsk. 1988. p. 1–177.
Carvalho-Filho FDS, Esposito MC. Revision of Argoravinia Townsend (Diptera: Sarcophagidae) of Brazil with the description of two new species. Zootaxa. 2012;3256:1–26.
Evans HE, Matthews RW, Callan EM. Observations of the nesting biology of Rubrica surinamensis (DeGeer) (Hymenoptera, Sphecidae). Psyche (Stuttg). 1974;81:334–352.
Jirón LF, Vargas LG, Vargas-Alvarado E. Four muscoid flies (Sarcophagidae and Muscidae) associated with human cadavers in Costa Rica. Brenesia. 1983;21:3–5.
Pape T. The world Blaesoxipha Loew, 1861 (Diptera: Sarcophagidae). Entomol Scand Suppl. 1994;45:1–247.
Dahlem G, Downes W. Revision of the genus Boettcheria in America North of Mexico (Diptera: Sarcophagidae). Insecta Mundi 1996 p. 77–103. http://digitalcommons.unl.edu/insectamundi/8.
Verves YG, Khrokalo LA. Family Sarcophagidae. In: Gerlach J, editor. Diptera of the Seychelles Islands. Pensoft Se. Sofia, Bulgaria: Pensoft Publishers; 2009. p. 270–303.
Lopes HS, Achoy RM. On Chrysagria (Diptera, Sarcophagidae), with descriptions of two new species, one of them living on larvae of Altina sp. (Coleoptera, Chrysomelidae). Rev Bras Biol. 1986;46:273–6.
Eberhard WG. Use of host aggregations as mating sites by Chrysagria alticophaga (Diptera, Sarcophagidae). Psyche (Stuttg). 1989;96(1–2):45–8.
Greene C. The puparia and larvae of sarcophagid flies. Proc United States Natl Museum. 1925;66(29):1–26.
Peters J.A. The box turtle as a host for dipterous parasites. Am Midl Nat. 1948;40:472–4.
Dodge HR. Sarcophagid flies parasitic on reptiles. (Diptera, Sarcophagidae). Proc Entomol Soc Washingt. 1955;57:183–7.
Boyer T, Boyer D. Chapter 7 – Turtles, tortoises and terrapins. In: Mader D, editor. Reptile Medicine and Surgery. 2nd editio. Saint Louis: W.B. Saunders; 2006. p. 78–99.
Neck R, Lopes HS. On some North American Sarcophagidae and Calliphoridae (Diptera) reared from Gastropoda. Rev Bras Biol. 1973;33:183–92.
Lopes HS. Sôbre uma espécie nova de Dexosarcophaga Townsend, 1917 (Dipt. Sarcophagidae) cujas larvas vivem in ninho de Camponotus (Hymenoptera, Formicidae). Rev Bras Biol. 1968;28:521–3.
Stucky BJ. Infection behavior, life history, and host parasitism rates of Emblemasoma erro (Diptera: Sarcophagidae), an acoustically hunting parasitoid of the cicada Tibicen dorsatus (Hemiptera: Cicadidae). Zool Stud. 2015;54:1–17.
Lopes HS. Contribuição ao conhecimento do gênero Emdenimyia (Diptera, Sarcophagidae). Rev Bras Biol. 1969;29:429 436.
Dahlem GA, Naczi RFC. Flesh flies (Diptera: Sarcophagidae) associated with north american Pitcher plants (Sarraceniaceae), with descriptions of three new species. Ann Entomol Soc Am. 2006 Mar 1;99(2):218–40.
Plank H. Natural enemies of the sugar cane moth stalk borer in Cuba. Ann Entomol Soc Am. 1929;22:621–40.
Stegmaier C. Notes on some Sarcophagidae (Diptera) reared from snails (Mollusca) in Florida. Florida Entomol. 1972;55(4):237–42.
Reeves WK, Pape T, Adler PH. Biological notes on New World Sarcophagidae (Diptera). Stud dipterologica. 2000;7(2):497–500.
Lopes HS. On Engelimyia, new genus for Sarcophaga cassidifera Engel (Diptera, Sarcophagidae). Rev Bras Biol. 1975;34:425–30.
Mulieri PR, Schaefer EF, Duré MI, González CE. A new flesh fly species (Diptera: Sarcophagidae) parasitic on leptodactylid frogs. Parasitol Res. 2018;117(3):809–18.
Gómez-hoyos DA, Suárez-joaqui T, Marín-gómez OH. Flesh fly myiasis (Diptera: Sarcophagidae) in Pristimantis thectopternus (Anura: Strabomantidae) from Colombia. Herpetolohy Notes. 2012;5:27–9.
Mello-Patiu CA, Luna-Dias C. Myiasis in the neotropical amphibian Hypsiboas beckeri (Anura: Hylidae) by a new species of Lepidodexia (Diptera: Sarcophagidae). J Parasitol. 2010;96(4):685–8. http://www.ncbi.nlm.nih.gov/pubmed/20486737.
Mulieri PR, Mello-Patiu CA, Aballay FH. Taxonomic revision of Lipoptilocnema (Diptera: Sarcophagidae), with notes on natural history and forensic importance of its species. J Med Entomol. 2016;7:1–26. http://jme.oxfordjournals.org/lookup/doi/https://doi.org/10.1093/jme/tjw099.
Lopes HS. Sôbre Malacophagomyia g.n. (Diptera, Sarcophagidae) cujas larvas vivem em cadávares de “Gastropoda” (Mollusca). Rev Bras Biol. 1966;26:315–21.
Curran CH, Walley G. Sarcophagidae. In: Curran CH, editor. The Diptera of Kartabo, Bartica District, British Guiana, with descriptions of new species from other British Guiana localities. New York: Bulletin American Museum of Natural History; 1934. p. 474–91.
Lopes HS. Contribuição ao conhecimento do gênero Udamopyga Hall e de outro sarphagideos que vivem em moluscos no Brasil (Diptera). Rev Entomol. 1940;11:924–54.
Souza CM, Pape T, Thyssen PJ. Oxysarcodexia Townsend, 1917 (Diptera: Sarcophagidae)-a centennial conspectus. Vol. 4841, Zootaxa. 2020. 1–126 p.
Knipling E. A comparative study of the first-instar larvae of the genus Sarcophaga (Calliphoridae, Diptera), with notes on the biology. J Parasitol. 1936;22(5):417–54.
Sanders D, Dobson R. The insect complex associated with bovine manure in Indiana. Ann Entomol Soc Am. 1966;59:955–9.
d’Almeida JM, Mello RP de. Comportamento de dípteros muscóides frente a substratos de oviposição, em laboratório, no Rio de Janeiro, RJ, Brasil. Mem Inst Oswaldo Cruz. 1996;91(1):131–6.
Buenaventura E, Pape T. Revision of the New World genus Peckia Robineau-Desvoidy (Diptera: Sarcophagidae). Zootaxa. 2013 Mar 11;3622(1):87.
Méndez J, Pape T. Biology and immature stages of Peckia gulo (Fabricius, 1805) (Diptera: Sarcophagidae). Stud dipterologica. 2003;9(2):371–4.
D’Almeida JM. Ovipositional substrates used by Calyptrate Diptera in Tijuca Forest, Rio de Janeiro. Mem Inst Oswaldo Cruz. 1994;89(2):261–4.
Lopes HS. Algumas espécies de Sarcophaga do Brasil. Rev Entomol. 1935;5:38–46.
Coffey M. Studies on the association of flies (Diptera) with dung in southeastern Washington. Ann Entomol Soc Am. 1966;59(1):207–18.
Pickens LG. The life history and predatory efficiency of Ravinia lherminieri (Diptera: Sarcophagidae) on the face fly (Diptera: Muscidae). Can Entomol. 1981;113:523–6.
Wong E, Dahlem G, Stamper T, Debry R. Discordance between morphological species identification and mtDNA phylogeny in the flesh fly genus Ravinia (Diptera: Sarcophagidae). Invertebr Syst. 2014;29:1–11.
Lopes HS. Sarcophagid flies (Diptera) from Pacatuba, State of Ceará. Rev Bras Biol. 1975;34:271–94.
Dias E, Neves D, Lopes HS. Estudos sobre Fauna de Sarcophagidae (Diptera) de Belo Horizonte - Minas Gerais I - Levantamento Taxonômico e Sinantrópico. Mem Inst Oswaldo Cruz. 1984;79:83–91.
Tibana R, Xerez R. Uma nova espécie de Retrocitomyia Lopes, 1982 (Diptera, Sarcophagidae). Rev Bras Biol. 1985;45:485–488.
Firmino Alves CA, dos Santos WE, Creão-duarte AJ. Diptera (Insecta) de importância forense da região Neotropical. Entomotropica. 2014;29(2):77–94.
De Sousa JRP, Da Silva Carvalho-Filho F, Esposito MC, De Meyer M. Distribution and abundance of Necrophagous Flies (Diptera: Calliphoridae and Sarcophagidae) in Maranhão, Northeastern Brazil. J Insect Sci. 2015;15:1.
Pereira De Sousa JR, Da Silva Carvalho-Filho F, Juen L, Esposito MC. Evaluating the effects of different vegetation types on necrophagous fly communities (Diptera: Calliphoridae; Sarcophagidae): Implications for conservation. PLoS ONE. 2016;11(10):1–23.
De Mello-Patiu CA, Salazar-Souza M. Retrocitomyia Lopes, 1982 (Diptera: Sarcophagidae): New species, new records, key to males, and an updated catalog. Zootaxa. 2016;4171(3):534–48.
Pape T, Méndez J. Two new species of Sarcofahrtiopsis (Diptera: Sarcophagidae). Zootaxa. 2004;485(1):1.
Pape T, Méndez J. A new species of Sarcofahrtiopsis Hall, 1933 from Panama (Diptera: Sarcophagidae). Ann Zool. 2002;52:339–42.
Curran CH. Sarcophaginae of the American Museum Congo Expedition (Diptera). Am Museum Novit. 1934;727:1–31.
Johnston TH, Tiegs CW. New and little-known sarcophagid flies from south-eastern Queensland. Proc R Soc Queensl. 1921;33:46–90.
Fuller M. Notes on the biology and early stages of Helicobia australis J. & T. (Sarcophagidae), a dipterous insect associated with grasshoppers. Proc Linnaean Soc New South Wales. 1938;63:133–8.
Hall SCB, Parmenter CJ. Necrotic egg and hatchling remains are key factors attracting dipterans to sea turtle (Caretta caretta, Chelonia mydas, Natator depressus) nests in Central Queensland, Australia. Copeia. 2008;(1):75–81.
Verves YG. Review of subtribes Harpagophallina and Heteronychiina (Diptera, Sarcophagidae). Zool Zhurnal. 1989;68(7):89–96.
Zumpt F. Myiasis in man and animals in the Old World. London: Butterworths; 1965. 267 p.
Coupland J, Barker G. Chapter 3 – Diptera as predators and parasitoids of terrestrial gastropods, with emphasis on Phoridae, Calliphoridae, Sarcophagidae, Muscidae and Fanniidae. In: Barker G, editor. Natural enemies of terrestrial molluscs. Wallingford: CABI Publishers; 2004. p. 85–158.
Pape T. The Sarcophagidae (Diptera) of Fennoscandia and Denmark. Fauna Entomol Scand. 1987;19:1–203.
Povolný D. Thyrsocnema platariae sp. n., a new species of flesh-fly from Greek Macedonia (Diptera, Sarcophagidae). Acta Entomol Bohemoslov. 1992;89:383–6.
Pape T. Family: SARCOPHAGIDAE. A taxonomic database to all flesh flies. Natural History Museum of Denmark. 2014. https://diptera.dk/sarco/index.php.
Eberhardt A. Untersuchungen über das Schmarotzen von Sarcophaga carnaria an Regenwürmern und Vergleich der Biologie einiger Sarcophaga Arten. Zeitschrift für Morphol und Ökologie der Tiere. 1955;43(616–647).
Enderlein G. Klassifikation der Sarcophagiden. Sarcophagiden-Studien I. Arch für Klassif und Phylogenetische Entomol. 1928;1:1–56.
Povolný D. A tentative list of fleshflies inhabiting alpine zone of the protected area Velká Fatra (Diptera, Sarcophagidae). Acta Univ Agric (Ser A). 1993;41:197–206.
Berner L. Sur le parasitisme des hélicidés par des mouches du genre Sarcophaga. Bull du Muséum d’Histoire Nat Marseille. 1973;33:87–94.
Verves YG. The study of sarcophagids (Diptera, Sarcophagidae) - parasites of terrestrial gastropods. Vestn Zool. 1976;3:28–50.
Povolný D, Verves YG. A preliminary list of Bulgarian Sarcophaginae (Diptera). Acta Entomol Musei Natl Pragae. 1990;43:283–328.
Sanjean J. Taxonomic studies of Sarcophaga larvae of New York, with notes on the adults. Mem New York Agric Exp Stn. 1957;349:1–115.
Pérez-Moreno S, Marcos-García MA, Rojo S. Comparative morphology of early stages of two Mediterranean Sarcophaga Meigen, 1826 (Diptera; Sarcophagidae) and a review of the feeding habits of Palaearctic species. Micron. 2006 Jan;37(2):169–79.
Lukin L. Human cutaneous myiasis in Brisbane: a prospective study. Med J Aust. 1989;150:237–40.
Martínez-Sánchez A, Marcos-Garcia MA, Rojo S. Biodiversidad de la comunidad de dipteros sarcosaprofagos en ambientes insulares del sudeste Iberico (Diptera: Calliphoridae, Muscidae, Sarcophagidae). Nouv Rev d’Entomologie. 2005;22:251–65.
Shiota T, Yoshida Y, Hirai S, Torii S. Intestinal myiasis caused by Parasarcophaga crassipalpis (Diptera, Sarcophagidae). Pediatrics. 1990;85:215–7.
Pape T. Two new species of Sarcophaga Meigen from Madeira and mainland Portugal (Diptera: Sarcophagidae). Tijdschr voor Entomol. 1990;133:39–42.
Blackith R, Blackith R, Pape T. Taxonomy and systematics of Helicophagella Enderlein, 1928 (Diptera, Sarcophagidae) with the description of a new species and a revised catalogue. Stud dipterologica. 1998;4(2):383–434.
Blackith RM, Pape T. Redescription of Sarcophaga discifera of authors, not Pandellé, and a new species of Sarcophaga Meigen from Oman (Diptera: Sarcophagidae). Stud dipterologica. 1999;6(2):239–50.
Giroux M, Wheeler TA. Systematics and phylogeny of the subgenus Sarcophaga (Neobellieria) (Diptera: Sarcophagidae). Ann Entomol Soc Am. 2009;102:567–87.
Prado e Castro C, García MD, Arnaldos MI, González-Mora D. Sarcophagidae (Diptera) atraídos a cadáveres de cochinillo, con nuevas citas para la fauna portuguesa. Graellsia. 2010;66(2):285.
Shazia, Tayyib M, Suhail A, Yousuf MJ. Systematics and Population of Sarcophagid Flies in Faisalabad (Pakistan). Int J Agric Biol. 2006;8(6):809–11.
Singh D, Bharti M. Some notes on the nocturnal larviposition by two species of Sarcophaga (Diptera: Sarcophagidae). Forensic Sci Int. 2008;177(1):19–20.
Blackith R, Blackith R, O’Connor J. A check-list of Irish flesh-flies (Diptera: Sarcophagidae: Sarcophagini) and their known distribution. Irish Nat J. 1998;24(11):427–34.
Xue WQ, Verves YG, Du J. A review of subtribe Boettcheriscina Verves 1990 (Diptera: Sarcophagidae), with descriptions of a new species and genus from China. Ann la Soc Entomol Fr. 2011;47(3–4):303–29.
Verves Y, Barták M, Kubík Š, Civelek HS. New records of Sarcophagidae from Turkey (Diptera). Zookeys. 2017;2017(703):129–58.
Eberhardt A, Steiner G. Untersuchungen über das Schmarotzen von Sarcophaga spp. in Regenwürmern. Zeitschrift für Morphol und Ökologie der Tiere. 1952;41:147–60.
Kirchberg E. Zur Larvennahrung einiger heimischer Sarcophaga-Arten, insbesondere zur Frage, ob S. carnaria L. als obligatorischer Regenwurmparasit anzusehen sei. (Diptera, Tachinidae). Zeitschrift für Morphol und Ökologie der Tiere. 1954;43:99–112.
Kirchberg E. Zucht von Sarcophaga carnaria L. (Diptera, Tachinidae) aus einer Freilandpopulation von Regenwürmern des Genus Allolobophora Eisen (Oligoch., Lumbricidae). (Zur Kenntnis der Gattung Sarcophaga Mg. III.). Anzeiger für Schädlingskunde, Pflanzenschutz und Umweltschutz. 1961;34(1):6–7.
Pape T, Bänziger H. Three new species of Sarcophaga Meigen found during ecological studies on flesh flies (Diptera: Sarcophagidae) in Thailand. Entomol Sci. 2003;6:49–56. https://doi.org/10.1046/j.1343-8786.2003.00004.x/full.
Aldrich JM. Sarcophaga and allies in North America, Vol. 1. La Fayette: Entomological Society of America, Thomas Say Foundation; 1916.
Hall DG. Biology of Sarothromyia femoralis var. simplex Aldrich. (Diptera, Calliphoridae). Ann Entomol Soc Am. 1932;25(3):641–647.
Vogt R. Turtle egg (Graptemys: Emydidae) infestations by fly larvae. Copeia. 1981;1981:457–9.
d’Almeida JM. Substratos para criação de dipteros caliptratos em área rural do Estado do Rio de Janeiro. Arq Univ Fed Rural do Rio Janeiro. 1987;9:13–22.
Stavert JR, Drayton BA, Beggs JR, Gaskett AC. The volatile organic compounds of introduced and native dung and carrion and their role in dung beetle foraging behaviour. Ecol Entomol. 2014;39:556–565.
Iyengar E V. Kleptoparasitic interactions throughout the animal kingdom and a re-evaluation, based on participant mobility, of the conditions promoting the evolution of kleptoparasitism. Biol J Linn Soc. 2008;93:745–762.
Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Evolution. 2008;62:1103–18.
Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 3.04 http://mesquiteproject.org. 2015.
Buenaventura E. Data from: Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae) [Dataset] [Internet]. Museum für Naturkunde Berlin (MfN)-Leibniz Institute for Evolution and Biodiversity Science. 2021. p. 1. https://doi.org/10.7479/t80n-ta02.
EB thanks Rebecca Dikow and Vanessa González (Smithsonian Data Science Lab), Matthew Kweskin (Laboratory of Analytical Biology, NMNH), and Michael Lloyd (The Jackson Laboratory) for their support during the UCE data analysis on the Smithsonian Institution High Performance Cluster Hydra. EB thanks Bonnie Blaimer (Museum für Naturkunde) for sharing scripts for coalescent-based analyses. EB extends special thanks to Torsten Dikow (NMNH) for his constant support during her work at the Smithsonian Institution and for teaming up to secure funding from the Global Genome Initiative. EB thanks Mauren Turcatel (Field Museum) for her companionship and support during the collecting trip to Florida and Wendy Porras for her assistance and support during the collecting trip to Costa Rica. EB thanks Arianna Thomas-Cabianca (CEUA, RMCA), John Stireman and Juan Manuel Perilla-López (JOSC), Marta Wolff (Universidad de Antioquia), Thomas Pape (NHMD), and Torsten Dikow (USNM), for facilitating accession to Museum specimens. EB thanks all of the anonymous reviewers for delivering a constructive and thorough review that greatly improved the manuscript at all levels.
Open Access funding enabled and organized by Projekt DEAL. EB was supported by a Peter Buck Postdoctoral Fellowship from the Smithsonian Institution and the National Museum of Natural History (NMNH) and received research funds from the Curtis W. Sabrosky Endowment Fund (with contributions by Holly Williams) from the Diptera unit in the Department of Entomology of the NMNH. EB received funding from the Global Genome Initiative of the Smithsonian Institution (GGI-Exploratory-2016-044 and GGI-Rolling-2016/2017). EB was supported by the Museum für Naturkunde during the analyses and writing stages of this project. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Specimens for this study were collected at the following locations, and with respective institutions providing authorizations for the collection and exportation (see Additional file 1 for collecting location and repository of spcimens). Costa Rica: Ministerio de Ambiente y Energía (Scientific Research Authorization Number R-045-2017-OT-CONAGEBIO to EB), Sistema Nacional de Áreas de Conservación (Permit Number 064-2017-ACC-PI to EB). USA: National Park Service; Florida Department of Environmental Protection (Permit Number 04111710 to EB); in addition to public land and roadside collecting sites in Washington and California do not require permitted collecting. Colombia: Entomological Collection Universidad de Antioquia, Autoridad Nacional de Licencias Ambientales (Authorization No. 0524 to Universidad de Antioquia). In addition, specimens (or only their body parts or DNA aliquots, see Additional file 1 for details) from the following locations were borrowed from respective museums, which guaranteed accession. South Africa, Kenya: CEUA and RMCA. Brazil, Guatemala, Iran, Italy, Peru, Turkey, Ukraine, United Arab Emirates, and Venezuela: MFN. Bahamas, Bolivia, Denmark, Ecuador, Israel, Malaysia, New Zealand, South Africa, and Venezuela: USNM. Canada, USA: JOSC. Australia, Burundi, Croatia, France, Kenya, Tanzania, Turkey, and Uganda: NHMD.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Specimen data. Specimen identity, preservation method, targeted tissue for extraction, collection data, specimen identifiers and repositories at Natural History Museums, DNA input data, contig data, UCE capture data. Abbreviations: CEUA, Entomological Collection Alicante University; JOSC, Entomological Collection Wright State University; MFN, Museum für Naturkunde Berlin; NHMD, Natural History Museum of Denmark; RMCA, Musée royal de l'Afrique centrale; USNM, National Museum of Natural History.
Phylogenetic relationships for Sarcophagidae and 10 other Calyptratae fly families inferred from 17 datasets having varying loci occupancy and coded as nucleotides and amino acids using a concatenated ML approach as well as by reconstructing a species tree estimated from UCE gene trees. Datasets are described in Table 1.
Comparison of models used for ancestral character state reconstructions. Models compared are ‘equal rates’ (ER, one transition rate) and ‘all rates different’ (ARD). Significantly better values are in bold. Characters analysed are: (1) larval food resource, (2) larval feeding habits, (3) abdominal ST5, shape of posterior margin, (4) cercal prong, outline of dorsal surface in lateral view, (5) phallus, connection between basi- and distiphallus, (6) phallus, connection between basi- and distiphallus, shape, (7) harpes, (8) vesica, (9) phallotrema, configuration, (10) phallotrema, position with regard to phallic tube, (11) acrophallic levers, (12) styli, number, (13) capitis, (14) median process, and (15) juxta.
Evolution of shape of posterior margin of the abdominal sternite 5 in Sarcophagidae. Ancestral character state reconstruction for shape of posterior margin of the abdominal sternite 5 using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of shape of posterior margin of the outline of dorsal surface of cercal prong in lateral view. Ancestral character state reconstruction for shape of posterior margin of the abdominal sternite 5 using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ARD) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the connection between basi- and distiphallus in Sarcophagidae. Ancestral character state reconstruction for connection between basi- and distiphallus using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the harpes in Sarcophagidae. Ancestral character state reconstruction for harpes using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ARD) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the vesica in Sarcophagidae. Ancestral character state reconstruction for vesica using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the configuration of the phallotrema in Sarcophagidae. Ancestral character state reconstruction for configuration of the phallotrema using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the phallotrema position with regard to phallic tube in Sarcophagidae. Ancestral character state reconstruction for phallotrema position with regard to phallic tube using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the acrophallic levers in Sarcophagidae. Ancestral character state reconstruction for acrophallic levers using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the capitis in Sarcophagidae. Ancestral character state reconstruction for capitis using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the median process in Sarcophagidae. Ancestral character state reconstruction for median process using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Evolution of the juxta in Sarcophagidae. Ancestral character state reconstruction for juxta using maximum likelihood and the rayDISC function in the R package corHMM. Only reconstruction of the best fitting model (ER) is shown. Pie proportions represent state probabilities estimated for each internal node. Character states are indicated in insets at the bottom.
Species richness matrix used for diversification analyses. This table lists the assigned clade for each species, the number of species sampled, and the number of nominal species described. From these cladespecific values, sampling fractions were calculated and used for BAMM analyses.
Dataset of biological characters including ingroup and outgroup. Characters: 1) Larval food resource: (0) invertebrates; (1) vertebrates; (2) feces; (3) plants; (4) ? (unknown); 2) larval feeding habits: (0) kleptoparasitism; (1) sarcosaprophagy; (2) predation; (3) coprophagy; (4) phytophagy; (5) ? (unknown).
Dataset of biological and morphological characters including ingroup only. Characters: 1) Larval food resource: (0) invertebrates; (1) vertebrates; (2) feces; (3) n/a (unknown); 2) larval feeding habits: (0) kleptoparasitism; (1) sarcosaprophagy; (2) predation; (3) coprophagy; (4) n/a (unknown); 3) abdominal ST5, shape of posterior margin: (0) straight or with a shallow concavity; (1) forming a cleft; 4) cercal prong, outline of dorsal surface in lateral view: (0) straight or almost straight; (1) swollen or curved uniformly; (2) with a proximal hump; (3) with a subapical saddle-shaped concavity followed by a hump; 5) phallus, connection between basi- and distiphallus: (0) continuous; (1) non-continuous; 6) phallus, connection between basi- and distiphallus, shape: (0) desclerotized band; (1) distinct hinge; (2) partially sclerotized; (3) fully sclerotized; 7) harpes: (0) absent; (1) present; 8) vesica: (0) absent; (1) present; 9) phallotrema, configuration: (0) not folded, forming one opening; (1) folded, forming three openings; 10) phallotrema, position with regard to phallic tube: (0) apical; (1) ventral; 11) acrophallic levers: (0) present; (1) absent; 12) styli, number: (0) one; (1) three; (2) two; 13) capitis: (0) present; (1) absent; 14) median process: (0) present; (1) absent or reduced; 15) juxta: (0) absent; (1) present.
About this article
Cite this article
Buenaventura, E. Museomics and phylogenomics with protein-encoding ultraconserved elements illuminate the evolution of life history and phallic morphology of flesh flies (Diptera: Sarcophagidae). BMC Ecol Evo 21, 70 (2021). https://doi.org/10.1186/s12862-021-01797-7
- Targeted enrichment
- Ultraconserved elements
- Rapid radiation