Skip to main content

Anchored enrichment dataset for true flies (order Diptera) reveals insights into the phylogeny of flower flies (family Syrphidae)

Abstract

Background

Anchored hybrid enrichment is a form of next-generation sequencing that uses oligonucleotide probes to target conserved regions of the genome flanked by less conserved regions in order to acquire data useful for phylogenetic inference from a broad range of taxa. Once a probe kit is developed, anchored hybrid enrichment is superior to traditional PCR-based Sanger sequencing in terms of both the amount of genomic data that can be recovered and effective cost. Due to their incredibly diverse nature, importance as pollinators, and historical instability with regard to subfamilial and tribal classification, Syrphidae (flower flies or hoverflies) are an ideal candidate for anchored hybrid enrichment-based phylogenetics, especially since recent molecular phylogenies of the syrphids using only a few markers have resulted in highly unresolved topologies. Over 6200 syrphids are currently known and uncovering their phylogeny will help us to understand how these species have diversified, providing insight into an array of ecological processes, from the development of adult mimicry, the origin of adult migration, to pollination patterns and the evolution of larval resource utilization.

Results

We present the first use of anchored hybrid enrichment in insect phylogenetics on a dataset containing 30 flower fly species from across all four subfamilies and 11 tribes out of 15. To produce a phylogenetic hypothesis, 559 loci were sampled to produce a final dataset containing 217,702 sites. We recovered a well resolved topology with bootstrap support values that were almost universally >95 %. The subfamily Eristalinae is recovered as paraphyletic, with the strongest support for this hypothesis to date. The ant predators in the Microdontinae are sister to all other syrphids. Syrphinae and Pipizinae are monophyletic and sister to each other. Larval predation on soft-bodied hemipterans evolved only once in this family.

Conclusions

Anchored hybrid enrichment was successful in producing a robustly supported phylogenetic hypothesis for the syrphids. Subfamilial reconstruction is concordant with recent phylogenetic hypotheses, but with much higher support values. With the newly designed probe kit this analysis could be rapidly expanded with further sampling, opening the door to more comprehensive analyses targeting problem areas in syrphid phylogenetics and ecology.

Background

Thanks in part to modern molecular techniques, the field of biological systematics has made great advances in assembling the Tree of Life. Well-supported phylogenetic hypotheses, based partly or entirely on phylogenomic datasets, now exist for many major animal groups, including (holometabolous) insects [1, 2], birds [35], mammals [6], and squamates [7]. Phylogenomic analyses have been made possible by the dramatically decreasing costs of genome/transcriptome sequencing of non-model organisms [8]. However, for many phylogenetic questions, a dense, comprehensive sampling of genomes/transcriptomes is a still prohibitively expensive enterprise. In order to generate these comprehensive phylogenomic data sets, several cost-effective alternatives to whole genome or transcriptome sequencing have been proposed.

One such method is hybrid enrichment [9, 10], which uses oligonucleotide probes or “baits” targeting specific areas of the genome in question. These probes hybridize to genomic fragments containing the loci of interest, allowing them to be amplified and sequenced using high-throughput sequencing. Originally developed for medical research on human diseases [10, 11], hybrid enrichment is a flexible technique for which applications in phylogenomic research are just beginning to be realized [4, 1214]. Unlike traditional polymerase chain reaction (PCR), hybrid enrichment techniques can be used to isolate and amplify many loci in a single reaction, and thus greatly improve the representation of single species in terms of gene coverage in phylogenomic analyses. Furthermore, once a probe kit is developed the cost of a project increases primarily by the number of taxa added (unlike Sanger sequencing which increases by the number of taxa and loci added) [13].

Two major hybrid enrichment methods are currently used for phylogenetic studies: the ultraconserved element (UCE) approach [12] and anchored hybrid enrichment (AHE) [13]. The UCE approach targets highly conserved noncoding regions of the genome [12] while AHE targets highly conserved regions primarily in the coding portion of the genome; specifically, it targets these regions flanked by less conserved regions in an attempt to acquire more data useful for phylogenetic inference [13]. AHE probe kits are also designed to target a wide range of taxonomic groups: the initial probe kit was designed for use across all vertebrate taxa. This was accomplished by comparing the complete genome of five model organisms [15]. While recent studies have used the UCE approach to study ants [16] and a related exon-capture method to study brittle stars [17], the present study is the first invertebrate project conducted using the AHE technique, utilizing the first iteration of insect-specific probes to construct a phylogenetic hypothesis of the dipteran family Syrphidae.

Syrphidae is a large and relatively well-known family of Diptera with over 6200 described species worldwide [18]. The family has traditionally been divided into three subfamilies: Syrphinae, Microdontinae, and Eristalinae [19]. However, Pipizini, a tribe of historically uncertain placement, has recently been elevated to subfamilial level (i.e. Pipizinae) [20]. In addition, latest phylogenetic studies using molecular sequence data [21] and combined molecular and morphological data [22] recover Eristalinae as paraphyletic. Finally, the Microdontinae have been alternately placed within what would now be considered Eristalinae [23, 24], within Syrphinae [25], or as a separate family [2628] but are currently considered a subfamily [22, 2935]. In summary, there is no phylogenetic consensus of subfamilial relationships.

The current tribal division of the family is based mostly on adult morphological characters and larval biology [36]. A total of 15 tribes are recognized: Microdontini and Spheginobacchini, in Microdontinae; Brachyopini, Callicerini, Cerioidini, Eristalini, Merodontini, Milesiini, Rhingiini, Sericomyiini, and Volucellini, in Eristalinae; and Bacchini, Paragini, Syrphini, and Toxomerini, in Syrphinae [20]. The subfamily Pipizinae has no tribal subdivision. The classification into tribes has not been generally accepted, and the relationships among them have never been studied in detail for the entire family [27, 3739]. Some of the genera have been placed in different tribes and some tribes have even been placed in different subfamilies. For instance, Spheginobacchini has been placed within eristalines, syrphines and microdontines [22, 40, 41] as well as “Pipizini” [20]. Moreover, some tribes are not supported by the last molecular phylogenetic studies, such as Brachyopini, Bacchini or Toxomerini, or their placement within a subfamily is uncertain or unresolved as there is no agreement among different works, e.g. Paragini, Volucellini, Merodontini, and Callicerini [20, 39, 42, 43].

Adults of most species of flower flies are conspicuous flower visitors, where they feed on both pollen and nectar [44]. This behaviour has earned the family the common name “flower flies” (also known as “hoverflies”), and has also generated a large amount of interest in the family as pollinators in both natural ecosystems and agricultural crops [4550]. The only exception are the microdontines, whose adults are rarely seen on flowers, and in some species they do not feed at all [51]. In contrast to the relatively uniform behaviour of the adults, syrphid larvae display an extraordinary diversity of life histories for a single family, including terrestrial and aquatic predators, inquilines in ant, wasp and bumblebee nests, saprophages, mycophages, root borers, stem miners, leaf miners, and wood borers in decaying logs [40, 52, 53]. Larvae of Microdontinae are inquilines in ants’ nests feeding on eggs, larvae and pupae [54], but also may parasite ant pupae [55]. Immature stages of Eristalinae include saprophages in a wide range of decaying organic media from dung to dead wood, some phytophages in various plants, and some predaceous species, i.e. species of the genus Volucella Geoffrey, 1762 are wasp- and bee-brood predators, and larvae of Nepenthosyrphus Meijere, 1932 are sit-and-wait aquatic predators in the phytotelmata of pitcher plants in SE Asia [40, 53, 5659]. Larvae of Pipizinae and Syrphinae share a similar feeding mode, but while known pipizine larvae are predatory mostly on woolly or root aphids with waxy secretions and gall-forming hemipterans, the majority of syrphine larvae prey on a broader range of soft-bodied arthropods such as aphids, coccids and psyllids, but also on Thysanoptera, immature Coleoptera, and Lepidoptera caterpillars [60]. The larvae of some Neotropical syrphines develop as stem borers and leaf miners in plants or as pollen feeders [6164]. This high diversity of natural histories makes syrphid immatures interesting and economically important as they can be biological control agents of plant pests and invasive weeds, re-cyclers of dead plant and animal matter, and pests of some ornamental plants [40, 53, 65, 66].

Hence, a robust phylogeny of syrphids is crucial to tackle the evolution of mimicry [67], to test the coevolution of microdontines and their ant hosts [54], to infer the evolution of larval life histories and the biology of the common ancestor, and to study the evolution of migratory behaviour.

The aim of the current study was to develop a set of AHE probes for use in Diptera, and to use the newly developed probe set to address the systematic position of the more problematic (e.g. unstable placements, unique morphology) taxa within Syrphidae, especially at the subfamilial and tribal level. Due to their high level of diversity, myriad of larval life histories, historical intractability of a robust subfamilial phylogenetic hypothesis, and economic and ecological significance, Syrphidae are an attractive model organism to test the utility of AHE. The project was accomplished by utilizing AHE to obtain genomic data from 559 nuclear gene regions (374 used in the final analyses). Although the main goal of this study was to elucidate phylogenetic relationships within the family Syrphidae, sequence data from a total of 12 cyclorrhaphan Diptera families were captured, illustrating the flexibility of the technique.

Although the current study includes all major clades of Syrphidae, the phylogeny proposed here will eventually form the basis for a much larger and more thoroughly sampled phylogenetic study (http://www.canacoll.org/Diptera/Staff/Skevington/Syrphidae/Syrphidae_World_Phylogeny.htm). This initiative is being conducted by a large group of entomologists and promises to be the largest phylogenetic collaboration attempted on a single family of insects.

Methods

Anchored hybrid enrichment laboratory data collection

Data were collected following the general methods of Lemmon et al. [13] through the Center for Anchored Phylogenomics at Florida State University (www.anchoredphylogeny.com). Briefly, 50ul of each genomic DNA sample, with quantity ranging from 11.5 to 985.3 ng) was sonicated to a fragment size of ~150-350 base pairs (bp) using a Covaris E220 Focused-ultrasonicator with Covaris microTUBES. Subsequently, library preparation and indexing were performed on a Beckman-Coulter Biomek FXp liquid-handling robot following a protocol modified from Meyer and Kirschner [68]. One important modification is a size-selection step after blunt-end repair using SPRIselect beads (Beckman-Coulter Inc.; 0.9× ratio of bead to sample volume). Indexed samples were then pooled at equal quantities (typically 12–16 samples per pool), and enrichments were performed on each multi-sample pool using an Agilent Custom SureSelect kit (Agilent Technologies), designed as specified above. After enrichment, the three enrichment pools were pooled in equal quantities for sequencing in one PE150 Illumina HiSeq2000 lane. Sequencing was performed in the Translational Science Laboratory in the College of Medicine at Florida State University.

Probe development

We began with nucleotide alignments of 4485 protein coding genes for 13 insect species identified by Niehuis et al. [69]. Each alignment contained up to 11 members of Holometabola from five orders (Diptera, Hymenoptera, Lepidoptera, Strepsiptera, and Coleoptera) and two non-holometabolous insects (used as outgroup) from two orders (Anoplura and Hemiptera). A full list of the species and their higher taxonomy is given in Table 1. We then selected a preliminary set of loci containing > =6 taxa and at least one consecutive 120 bp region with >50 % pairwise sequence identity. Sequences for each species were extracted, and exon boundaries were then identified using published genomes (see Table 1 for details) and custom scripts that identified matches between the transcript sequences (Table 2) and the genomes using 40-mers.

Table 1 Voucher specimens used to determine exon boundaries for initial probe site selection
Table 2 Diptera genomes and transcriptomes used to develop probe kit

Together with the alignments, the exon boundaries were used to identify suitable candidate regions (exons) to target using an Anchored Phylogenomics approach, as described by Lemmon et al. [13]. The following requirements were used to select 962 insect-wide targets: 1) the region was at least 150 bp in length, 2) the region contained no exon boundaries, and 3) the region contained no indels. Details of these targets are given in Additional file 1: Table S1. Concatenated alignments have been uploaded to the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra), with accession numbers (Biosample #) available in Table 3. The lengths of these targets ranged from 150 to 863 bp (mean = 187 bp) whereas the pairwise sequences similarity ranged from 45 to 84 % (mean = 66 %).

Table 3 Voucher specimens used in phylogenetic analysis. JSS = Jeff Skevington Specimen. All vouchers deposited in CNC

In order to develop an enrichment kit efficient for Diptera, we developed a reference database based on the Drosophila melanogaster sequences contained within the 962 target locus alignments, plus 13 established loci provided by Brian Wiegmann [70]. The database contained spaced k-mers derived from conserved sites within each locus. These were used to scan for homologous loci in seven Diptera genomes and 14 Diptera transcriptomes (see Table 2 for complete list). After the sequence best matching to the references was identified for each species x locus combination, alignments were estimated for each locus using MAFFT (Katoh and Standley, 2013; v7.023b with -genafpair and -maxiterate 1000 flags) [71]. Geneious v5.6.4 (Biomatters, available from http://www.geneious.com) was then used to select well-aligned regions that overlapped with the core insect regions, contained high taxon representation (>10 of 21 lineages), and contained low gaps. The 546 chosen anchor locus alignments contained 121–1497 sites (average of 588 sites) and 48 %-84 % pairwise sequence similarity (average = 69 %). The 13 functional locus alignments contained 185–3035 sites (average of 1758 sites) and 50 %-79 % pairwise sequence similarity (average = 66 %).

Finally, in order to ensure efficient enrichment, we checked for high-copy regions (e.g. microsatellites and transposable elements) in each of the seven genome-derived references as follows. First, a database was constructed for each species using all 15-mers found in the trimmed alignments for that species. We also added to the database all 15-mers that were 1 bp removed from the observed 15-mers. The genome for the species was then exhaustively scanned for the presence of these 15-mers and matches were tallied at the alignment positions at which the 15-mer was found. Alignment regions containing > 100,000 counts in any of the seven species were masked to prevent probe tiling across these regions. Probes of 120 bp were tiled uniformly at 1.72× tiling density (57,681 probes total). Final probe regions and the final probe sequences are available as Additional file 2: Table S2 and Additional file 3: Table S3. Scripts used for locus selection and design and alignments are available upon request from ARL.

In essence, the process for choosing probes for the Diptera kit was fundamentally the same as for choosing probes for the vertebrate kit (V1, Lemmon et al. 2012 [13]). The only difference was that alignments containing only genomes formed the basis of the vertebrate kit, whereas alignments containing both genomes and transcriptomes formed the basis of the Diptera kit.

Anchored hybrid enrichment bioinformatic data analysis

Paired-read merging

Typically, between 50 and 75 % of sequenced library fragments had an insert size between 150 and 300 bp. Since 150 bp paired-end sequencing was performed, this means that the majority of the paired reads overlap and thus should be merged prior to assembly. The overlapping reads were identified and merged following Rokyta [72]. In short, for each degree of overlap for each read we computed the probability of obtaining the observed number of matches by chance, and selected degree of overlap that produced the lowest probability, with a p-value less than 10−10 required to merge reads. When reads are merged, mismatches are reconciled using base-specific quality scores, which were combined to form the new quality scores for the merged read (see [72] for details). Reads failing to meet the probability criterion were kept separate in the assembly. The merging process produces three files one containing merged reads and two containing the unmerged reads.

Assembly

The reads were assembled into contigs using an assembler that makes use of both a divergent reference assembly approach to map reads to the probe regions and a de-novo assembly approach to extend the assembly into the flanks. The reference assembler uses a library of spaced 20-mers derived from the conserved sites of the alignments used during probe design. A preliminary match was called if at least 17 of 20 matches exist between a spaced k-mer and the corresponding positions in a read. Reads obtaining a preliminary match were then compared to an appropriate reference sequence used for probe design to determine the maximum number of matches out of 100 consecutive bases (all possible gap-free alignments between the read and the reference were considered). The read was considered mapped to the given locus if at least 55 matches were found. Once a read was mapped, an approximate alignment position was estimated using the position of the spaced 20-mer, and all 60-mers existing in the read were stored in a hash table used by the de-novo assembler. The de-novo assembler identified exact matches between a read and one of the 60-mers found in the hash table. Simultaneously using the two levels of assembly described above, the three read files were traversed repeatedly until an entire pass through the reads produced no additional mapped reads.

A list of all 60-mers found in the mapped reads was compiled, the 60-mers were clustered if found together in at least two reads. The 60-mer clusters were then used to separate the reads into clusters for contig estimation. Relative alignment positions of reads within each cluster were then refined in order to increase the agreement across the reads. Up to one gap was also inserted per read if needed to improve the alignment. Note that given sufficient coverage and an absence of contamination, each single-copy locus should produce a single assembly cluster. Low coverage (leading to a break in the assembly), contamination, and gene duplication, can all lead to an increased number of assembly clusters. A whole genome duplication, for example, would increase the number of clusters to two per locus.

Consensus bases were called from assembly clusters as follows. For each site an unambiguous base was called if the bases present were identical or if the polymorphism of that site could be explained as sequencing error, assuming a binomial probability model with the probability of error equal to 0.1 and alpha equal to 0.05. If the polymorphism could not be explained as sequencing error, the ambiguous base was called that corresponded to the IUPAC code. Called bases were soft-masked (made lowercase) for sites with coverage lower than five. A summary of the assembly results is presented in Additional file 4: Table S4.

Contamination filtering

In order to filter out possible low-level contaminants, consensus sequences derived from very low coverage assembly clusters (<10 reads) were removed from further analysis. After filtering, consensus sequences were grouped by locus (across individuals) in order to produce sets of homologs.

Orthology

Orthology was determined for each locus as follows. First, a pairwise distance measure was computed for pairs of homologs. To compute the pairwise distance between two sequences, we computed the percent of 20-mers observed in the two sequences that were found in both sequences. Note that the list of 20-mers was constructed from consecutive 20-mers as well as spaced 20-mers (every third base), in order to allow increased levels of sequence divergence. Using the distance matrix, we clustered the sequences using a Neighbor-Joining algorithm as follows: Pairwise distances were ranked from smallest to largest. Starting with the smallest value, pairs of sequences from the set of homologs (representing the next distance in the list) were joined into the same cluster. If one of the two sequences was already in a cluster, the clusters were merged. Clusters containing homologs originating from the same individual were not joined, such that when clustering was complete, each cluster contained at most one homolog per species. Sequence clusters containing fewer than 50 % of the species were removed from downstream processing.

Alignment (MAFFT)

Sequences in each orthologous set were aligned using MAFFT v7.023b [71], with --genafpair and --maxiterate 1000 flags.

Alignment trimming

In order to reduce the error in the data, the alignment for each locus was then trimmed/masked using the following procedure. First, each alignment site was identified as "conserved" if the most common character observed was present in > 40 % of the sequences. This step identified sites for which we were confident were aligned correctly for a sufficient portion of the taxa (typically third codon potions would not be included here). Second, 20 bp regions of each sequence that contained < 10 stable sites were masked. This step identified regions of each sequence that were not well aligned to the majority of the sequences and thus should be masked. Third, sites with fewer than 12 unmasked bases were removed from the alignment. This step identified large regions of the alignments that should be removed entirely from the alignment because they contain large quantities of missing data [73].

Taxon sampling

Representatives of all four Syrphidae subfamilies and 11 tribes were analysed. We also included taxa of another four dipteran families, i.e. Platypezidae [Platypeza sp.], Pipunculidae [Chalarus spurius (Fallén, 1816) and Pipunculus sp.ON12], Sepsidae [Themira nigricornis (Meigen, 1826)], and Tachinidae [Epalpus signifer (Walker, 1849)]. A total of 30 flower fly species were sampled (Table 3). Syrphid taxa come from four different Biogeographical Regions, but the majority are Nearctic specimens. Morphological identification of syrphids and pipunculids were provided by A.D.Y and J.H.S., other outgroup taxa were morphologically identified by colleagues at the Canadian National Collection of Insects, Arachnids, and Nematodes (CNC).

DNA extraction

Genomic DNA extractions were obtained with the QIAGEN DNeasy kit (Qiagen Inc., Santa Clara, CA, USA). Full specimens were extracted overnight at 56 °C, and total DNA was purified the following day following the manufacturer’s protocol. Following extraction, specimens were critical-point dried with the EM CPD300 (Leica Microsystems, Vienna, Austria) and deposited at CNC.

Vouchers

Specimens for the study were collected by Malaise trap or hand-collecting, preserved in 95-100 % ethanol, and placed in a −80 °C freezer until extraction. The voucher data and unique identifiers for the specimens used for the molecular study are presented in Table 3. Specimens have since been critical point dried, mounted, labeled and deposited in the Canadian National Collection of Insects, Arachnids and Nematodes.

The 5' region of the mitochondrial Cytochrome c Oxidase Subunit I (COI) gene was sequenced for each specimen in order to act as a surrogate voucher and allow linkage of the exemplars to a large molecular dataset being assembled. Amplification, purification, sequencing and contig assembly were carried out as described in Gibson et al. [74].

COI sequence alignment was straightforward as no indels (insertions or deletions) were found. The alignment was made by hand using Mesquite v2.74 [75] and translated into amino acids to ensure that there were no stop codons. Sequences were submitted to BOLD and uploaded from there to GenBank (Table 3).

Phylogeny estimation

A maximum likelihood (ML) tree (with 100 bootstrap replicates) for a single concatenated matrix was estimated using RAxML v7.2.6 [76], with the GTR + G substitution model partitioned by locus under default parameters. Platypeza was used to root the tree.

Results

Trimmed alignments contained 35 taxa and 217,702 sites (across 343 chosen loci), of which 89,534 sites were informative. The concatenated dataset was largely complete, with only 6 % missing data. Maximum Likelihood estimation (Fig. 1) of the present concatenated dataset produced a fully resolved tree, with 31/32 nodes (97 %) supported by >95 % bootstrap support (BS) values. As expected, Syrphidae was recovered as a monophyletic group with Microdon Meigen, 1803 as the sister to other lineages (BS = 100 %). The sister clade to the Syrphidae included Pipunculidae + Schizophora. The subfamilies Pipizinae and Syrphinae were resolved as clades. The potential monophyly of the subfamily Microdontinae could not be established (only one taxon included) and Eristalinae was resolved as non-monophyletic. A paraphyletic Eristalinae was placed sister to Syrphinae + Pipizinae. Within the eristalines, several tribes were resolved monophyletic based on the studied taxa. Merodontini (Eumerus Meigen, 1822 + Merodon Meigen, 1803) was recovered as a clade sister to the remainder of the Eristalinae. Volucellini (Graptomyza Weidemann, 1820 + Copestylum Macquart, 1846), Rhingiini (Rhingia Scopoli, 1763 + Cheilosia Meigen, 1822 + Ferdinandea Rondani, 1844) and Milesiini (Brachypalpus Macquart, 1834 + Xylota Meigen, 1822) were also found to be monophyletic. The three remaining tribes that were included in the analysis (Eristalini, Brachyopini, and Callicerini) only had a single member included, so potential monophyly could not be established. Within Syrphinae, three of the four tribes were included, i.e. Syrphini, Bacchini, and Paragini, but not Toxomerini. Bacchini was recovered as paraphyletic, with Melanostoma Schiner, 1860 placed as sister to the remainder of the Syrphinae, and a clade consisting of Baccha Fabricius, 1805 + Platycheirus Lepeletier & Serville, 1828 sister to Syrphinae excluding Melanostoma. Syrphini is a large tribe comprised of the majority of the syrphine genera, and formed a single clade with Paragus (the sole member of Paragini) resolved within it.

Fig. 1
figure 1

The ML phylogenetic tree based on the sequenced taxa using RAxML under the model GTR + G. Bootstrap support values are depicted above the nodes. Legend: black: outgroups; green: Microdontinae; orange: Eristalinae; red: Pipizinae; and blue: Syrphinae

Discussion

This analysis represents the first iteration of newly developed Diptera probes for AHE. While the probes were developed by analysing the genome of only 21 insect species, they were successfully used to extract sequence data from 18 additional Dipteran families (data not shown). Furthermore, while 559 loci were targeted designed, only 343 loci were included in the final analysis in order to minimize missing data. As more invertebrate genomes become available and probe kits are refined, ever larger datasets will be attainable from a broad spectrum of invertebrate taxa for a fraction of the cost of traditional Sanger sequencing methods [13].

The ML analysis produced a fully-resolved phylogram, with only one node with low bootstrap support (BS = 68 %) (see Fig. 1). While previous analyses have recovered similar phylogenies [2022], no previous works have recovered a fully-resolved tree with high support. A possible explanation for this surprising result is the high number of loci and bp included in our analysis bases on the newly-designed probes, which might allow fully resolved phylogenies for other dipteran families to be produced. The present analysis includes the largest genomic dataset ever created for the phylogenetic analysis of an insect/Diptera family, with 343 loci and 217,702 bp.

The two Pipunculidae taxa were recovered as sister to Epalpus signifer (Tachinidae) and Themira nigricornis (Sepsidae), both schizophoran flies. Although traditional morphological analyses [7780] have supported a sister group relationship between Syrphoidea (Pipunculidae + Syrphidae) and Schizophora, more recent morphological [81] and molecular [70] analyses suggest a sister group relationship between Pipunculidae and Schizophora, rendering Syrphoidea paraphyletic.

Placement of Microdontinae has a chequered history as pointed out in the introduction. The “ant flies” are morphologically very distinct from the remaining Syrphidae and all species with known larval histories are associated with ants. Larvae are either predatory or parasitoids in ant nests and have developed elaborate pheromone mimicry to carry out this feat [40, 55, 82, 83]. Strong morphological and ecological specializations within the group have made microdontines very difficult to place into phylogenetic context. Thompson [26] was the first to provide quantitative evidence that they are sister to all other Syrphidae species (based on adult morphology). Despite this, other contradictory hypotheses have continued to be proposed. For example, in their study of larval characters and evolution, Rotheray and Gilbert [84] presented a hypothesis supporting a sister-group relationship between Microdontinae and pipizines and syrphines. This hypothesis assumed a single predatory larval lineage within Syrphidae. Our study refutes this and supports Thompson [26] and several recent molecular studies using Sanger sequence data [2022, 35]. Proposals as per Thompson [27] and Speight (1987, 2014) [28, 85] have been made to elevate the ant flies to family status and although our present results do not refute this, it remains an argument largely based on the perceived level of morphological and ecological difference of ant flies from other syrphids. Microdontinae is a highly diverse clade and still understudied taxonomically and biologically [35, 41, 54, 55]. Only one species was available for the present study, but the inclusion of members of the Spheginobacchini as well as other taxa not closely related to the genus Microdon [41] will allow testing the relationships among the taxa of this subfamily and having a larger support on its placement among flower flies.

Eristalinae was recovered as paraphyletic in the present study. The monophyly of Eristalinae is supported by several studies and the currently followed classification follows this line of reasoning [23, 36, 86, 87]. In contrast, evidence from more recent surveys using adult morphological and/or molecular characters, with a very limited number of loci, resolve Eristalinae as paraphyletic [2022, 42]. Our analysis is the first to use AHE data from hundreds of loci, and ML analysis of the data provides support for a non-monophyletic Eristalinae (Fig. 1). In the present study, Merodontini was resolved as sister group of the other eristalines + syrphines + pipizines, and Volucellini and Helophilus Meigen, 1822 (Eristalini) were recovered in different nodes, with the other included eristaline tribes forming a clade, i.e. Rhingiini, Sphegina Meigen, 1822 (Brachyopini), Callicerini and Milesiini. Our taxon sampling is not enough to make conclusions about the tribal relationships within this subfamily. Consequently, a larger and broader taxon sampling is still required, including tribes that were not available for the present study such as Cerioidini and Sericomyiini, to understand how eristaline tribes are related. The only weakly supported node on the maximum likelihood tree is within the Eristalinae. Eristalinae is the subfamily with the highest number of species and larval biology diversity, and it is reflected in the classification with the recognition of nine tribes and several subtribes. Addition of more taxa in future studies will address the question of the monophyly of the subfamily and the tribes, and will also help to better understand larval evolution within this incredibly diverse group of flies.

Syrphinae and Pipizinae were reciprocally monophyletic and sister groups to each other. The placement of Pipizinae as sister to Syrphinae is a phylogenetic hypothesis that has gained increasing support in recent years, and last phylogenetic works have recovered Pipizinae either within Syrphinae [39, 88], or sister to it [2022, 84]. The frequent placement of Pipizinae within Eristalinae owes much to the fact that many early classification schemes were based largely or entirely on adult morphological characters. The present results strongly suggest a common origin of these two groups, which implies that predatory larvae feeding on soft-bodied arthropods have evolved only once in the evolution of the Syrphidae, and they corroborate previous surveys and the recent elevation of Pipizinae to subfamilial level [20]. Future studies will explore the interrelationships of the members of this subfamily and will test the hypothesis exposed by Vujić et al. [89].

Finally, the resolution of Syrphinae as a monophyletic group was not unexpected as virtually all existing flower fly phylogenies hypothesize that Syrphinae is a clade. In contrast, the current tribal classification within Syrphinae is not supported in our analyses in concordance with the last phylogenetic studies [20, 22, 39, 84, 90]. Bacchini was found to be paraphyletic, and its members (Melanostoma, Platycheirus and Baccha) were resolved in two groups, partly in agreement with previous studies [20, 22, 39]. Paragini, a syrphine tribe of historically uncertain placement, was resolved as sister to Scaeva Fabricius, 1805 + Betasyrphus Matsumura, 1917, making the current tribe Syrphini paraphyletic. Our results corroborate the hypothesis by Rotheray and Gilbert [38], using larval morphological characters, and by Mengual [91] and Mengual et al. [20], using molecular data alone or in combination with adult morphological characters respectively. Addition of more taxa and the inclusion of the tribe Toxomerini will help to understand the tribal classification of Syrphinae, to define new tribal groups, and, the most important, to study the evolution of predation within this group to answer why and how some taxa became phytophagous secondarily.

The scenario recovered in the present analysis using AHE data shows that predation evolved at least three times in different groups with distinct feeding strategies, viz. Pipizinae + Syrphinae, Microdontinae and Volucellini (although the genus Volucella was not studied). A key piece into this puzzle is the unknown biology of the immatures of Spheginobacchini, which would help to understand the relation between microdontines and the rest of flower flies. Excellent mimics of wasps and bumblebees appear in several groups, especially within Eristalinae in genera like Temnostoma Lepeletier and Serville, 1828, Spilomyia Meigen, 1803 or Volucella. The existence of a broad spectrum from non-mimics, through partial or imperfect mimics, to perfect mimics might indicate a multiple origin for mimicry. The same scenario is found when migratory species are taken into consideration based in our results. Species like Episyrphus balteatus (De Geer, 1776), Sphaerophoria scripta (Linnaeus, 1758) or members of Scaeva, Platycheirus and Helophilus are well-known migrants but little has been studied about the characteristics, origin and mechanisms of these migrations. A fully resolved exhaustively sampled phylogeny based on AHE has the potential to resolve these questions.

Conclusions

This is the first time that AHE technique is used on an extended and very diverse group of insects and represents the largest dataset assembled to bear on the phylogeny of a dipteran group. The price and repeatability using the present probe kit makes this technique a reliable methodology for future research using large output sequence datasets. Present results corroborate a number of earlier findings and hypotheses, although this dataset should be considered preliminary due to the small taxon sample.

The next step, that is building upon a framework with more thorough taxon sampling of the many morphologically highly diverse groups, will create the most comprehensive hypothesis ever made for a large lineage of flies. With such a high level of ecological and morphological diversity, a detailed phylogeny of Syrphidae will support future work in fields such as pollination biology and biological control, and will help to answer major challenging questions that remain open, such as the evolution of inquiline-host associations in myrmecophilic flies, the evolution of larval feeding behaviour, the development of perfect and imperfect mimicry, the origin and biogeography of the different taxon groups, as well as the patterns of migratory behaviour. As it stands, this study provides a test for previous phylogenetic work on syrphids and illustrates that anchored hybrid enrichment is a useful technique for rapidly assembling comprehensive, large datasets for phylogenetic hypothesis testing. Current anchored data collection and analysis pipelines allow 96 samples to be processed in as little as 3 weeks, from DNA extracts to trimmed alignments and preliminary phylogeny estimates (www.anchoredphylogeny.com).

Abbreviations

AHE, anchored hybrid enrichment; bp, base pairs; BS, bootstrap support; CNC, Canadian National Collection of Insects, Arachnids, and Nematodes; COI, Cytochrome c Oxidase Subunit I; indel, insertion or deletion; ML, maximum likelihood; PCR, polymerase chain reaction; UCE, ultraconserved element.

References

  1. Beutel RG, Friedrich F, Hörnschemeyer T, Pohl H, Hünefeld F, Beckmann F, Meier R, Misof B, Whiting MF, Vilhelmsen L. Morphological and molecular evidence converge upon a robust phylogeny of the megadiverse Holometabola. Cladistics. 2010;26:1–15.

    Article  Google Scholar 

  2. Misof B, Liu S, Meusemann K, Peters RS, Donath A, Mayer C, Frandsen PB, Ware J, Flouri T, Beutel RG. Phylogenomics resolves the timing and pattern of insect evolution. Science. 2014;346(6210):763–7.

    Article  CAS  PubMed  Google Scholar 

  3. Hackett SJ, Kimball RT, Reddy S, Bowie RC, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han K-L, Harshman J. A phylogenomic study of birds reveals their evolutionary history. Science. 2008;320(5884):1763–8.

    Article  CAS  PubMed  Google Scholar 

  4. McCormack JE, Harvey MG, Faircloth BC, Crawford NG, Glenn TC, Brumfield RT. A phylogeny of birds based on over 1,500 loci collected by target enrichment and high-throughput sequencing. PLoS One. 2013;8(1), e54848.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, Lemmon AR. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature. 2015;526:569–73.

    Article  CAS  PubMed  Google Scholar 

  6. Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simão TL, Stadler T. Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammal diversification. Science. 2011;334(6055):521–4.

    Article  CAS  PubMed  Google Scholar 

  7. Pyron RA, Burbrink FT, Wiens JJ. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evol Biol. 2013;13(1):93.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Glenn TC. Field guide to next‐generation DNA sequencers. Mol Ecol Resour. 2011;11(5):759–69.

    Article  CAS  PubMed  Google Scholar 

  9. Albert TJ, Molla MN, Muzny DM, Nazareth L, Wheeler D, Song X, Richmond TA, Middle CM, Rodesch MJ, Packard CJ. Direct selection of human genomic loci by microarray hybridization. Nat Methods. 2007;4(11):903–5.

    Article  CAS  PubMed  Google Scholar 

  10. Gnirke A, Melnikov A, Maguire J, Rogov P, LeProust EM, Brockman W, Fennell T, Giannoukos G, Fisher S, Russ C. Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing. Nat Biotechnol. 2009;27(2):182–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Bamshad MJ, Ng SB, Bigham AW, Tabor HK, Emond MJ, Nickerson DA, Shendure J. Exome sequencing as a tool for Mendelian disease gene discovery. Nat Rev Genet. 2011;12(11):745–55.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

  13. Lemmon AR, Emme SA, Lemmon EM. Anchored hybrid enrichment for massively high-throughput phylogenomics. Syst Biol. 2012;61(5):727–44.

    Article  CAS  PubMed  Google Scholar 

  14. Eytan RI, Evans BR, Dornburg A, Lemmon AR, Lemmon EM, Wainwright PC, Near TJ. Are 100 enough? Inferring acanthomorph teleost phylogeny using Anchored Hybrid Enrichment. BMC Evol Biol. 2015;15(1):1.

    Article  Google Scholar 

  15. Lemmon EM, Lemmon AR. High-throughput genomic data in systematics and phylogenetics. Annu Rev Ecol Evol Syst. 2013;44:99–121.

    Article  Google Scholar 

  16. Blaimer BB, Brady SG, Schultz TR, Lloyd MW, Fisher BL, Ward PS. Phylogenomic methods outperform traditional multi-locus approaches in resolving deep evolutionary history: a case study of formicine ants. BMC Evol Biol. 2015;15(1):1.

    Article  Google Scholar 

  17. Hugall AF, O’Hara TD, Hunjan S, Nilsen R, Moussalli A. An exon-capture system for the entire class Ophiuroidea. Mol Biol Evol. 2015;33(1):281–94.

  18. Pape T, Thompson FC. Systema Dipterorum. The Biosystematic Database of World Diptera. Version 1.5 2013. http://www.diptera.org/. Accessed 30 Apr 2015.

  19. Vockeroth JR, Thompson FC. Syrphidae. In: McAlpine JF, Peterson BV, Shewell GE, Teskey HJ, Vockeroth JR, Wood DM, editors. Manual of Nearctic Diptera, vol. 2. Ottawa: Canadian Government Publishing Centre; 1987. p. 713–43.

    Google Scholar 

  20. Mengual X, Ståhls G, Rojo S. Phylogenetic relationships and taxonomic ranking of pipizine flower flies (Diptera: Syrphidae) with implications for the evolution of aphidophagy. Cladistics. 2015;31(5):491–08.

  21. Skevington JH, Yeates DK. Phylogeny of the Syrphoidea (Diptera) inferred from mtDNA sequences and morphology with particular reference to classification of the Pipunculidae (Diptera). Mol Phylogen Evol. 2000;16(2):212–24.

    Article  CAS  Google Scholar 

  22. Ståhls G, Hippa H, Rotheray G, Muona J, Gilbert F. Phylogeny of Syrphidae (Diptera) inferred from combined analysis of molecular and morphological characters. Syst Entomol. 2003;28(4):433–50.

    Article  Google Scholar 

  23. Goffe CER. An outline of a revised classification of the Syrphidae (Diptera) on phylogenetic lines. Trans Soc Br Entomol. 1952;11(4):97–124.

    Google Scholar 

  24. Wirth WW, Sedman YS, Weems HVJ. Family Syrphidae. In: Stone A, Sabrosky CW, Wirth WW, Foote RH, Coulson JR, editors. A catalog of the Diptera of America north of Mexico. Vol. Agriculture Handbook no. 276. Washington: United States Department of Agriculture; 1965. p. 557–625.

    Google Scholar 

  25. Williston SW. Synopsis of the North American Syrphidae. Bulletin of the United States National Museum. 1886;31:i–xxx. 1–335.

    Google Scholar 

  26. Thompson FC. A new genus of Microdontine flies (Diptera: Syrphidae) with notes on the placement of the subfamily. Psyche. 1969;76:74–85.

    Article  Google Scholar 

  27. Thompson FC. A contribution to a generic revision of the Neotropical Milesinae (Diptera: Syrphidae). Arquivos de Zoologica. 1972;23(2):73–215.

    CAS  Google Scholar 

  28. Speight MCD. External morphology of adult Syrphidae (Diptera). Tijdschr Entomol. 1987;130:141–75.

    Google Scholar 

  29. Hull FM. The morphology and inter-relationship of the genera of syrphid flies, recent and fossil. Tran Zool Soc Lond. 1949;26:257–408.

    Article  Google Scholar 

  30. Knutson LV, Thompson FC, Vockeroth JR. Family Syrphidae. In: Delfinado MD, Hardy DE, editors. A Catalog of the Diptera of the Oriental Region Volume 2 Suborder Brachycera Through Division Aschiza, Suborder Cyclorrhapha, vol. 2. Honolulu: The University Press of Hawaii; 1975. p. 307–74.

    Google Scholar 

  31. Smith KGV, Vockeroth JR. 38. Family Syrphidae. In: Crosskey RW, editor. Catalogue of the Diptera of the Afrotropical Region. London: Publications of the British Museum (Natural History); 1980. p. 488–510.

    Google Scholar 

  32. Thompson FC, Vockeroth JR. 51. Family Syrphidae. BMSP. 1989;86:437–58.

    Google Scholar 

  33. Thompson FC, Vockeroth JR, Sedman YS. 46. Family Syrphidae. In: A Catalogue of the Diptera of the Americas South of the United States. Vol. Part 3. Cyclorrhapha. São Paulo: Museu de Zoologia, Universidade de São Paulo; 1976. p. 1–195.

    Google Scholar 

  34. Cheng X-Y, Thompson FC. A generic conspectus of the Microdontinae (Diptera: Syrphidae) with the description of two new genera from Africa and China. Zootaxa. 2008;1879:21–48.

    Google Scholar 

  35. Reemer M, Ståhls G. Generic revision and species classification of the Microdontinae (Diptera, Syrphidae). ZooKeys. 2013;288:1–213.

    Article  PubMed  Google Scholar 

  36. Thompson FC, Rotheray G. Family Syrphidae. In: Papp L, Darvas B, editors. Contributions to a manual of Palaearctic Diptera (with species reference to flies of economic importance). Vol. 3, Higher Brachycera. Budapest: Science Herald; 1998. p. 81–139.

    Google Scholar 

  37. Thompson FC. Notes on the status and relationships of some genera in the tribe Milesiini (Diptera: Syrphidae). Proc Entomol Soc Wash. 1975;77(3):291–305.

    Google Scholar 

  38. Rotheray GE, Gilbert FS. The phylogeny and systematics of European predacious Syrphidae (Diptera) based on larval and puparial stages. Zool J Linn Soc. 1989;95(1):29–70.

    Article  Google Scholar 

  39. Mengual X, Stahls G, Rojo S. First phylogeny of predatory flower flies (Diptera, Syrphidae, Syrphinae) using mitochondrial COI and nuclear 28S rRNA genes: conflict and congruence with the concurrent tribal classification. Cladistics. 2008;24:543–62.

    Article  Google Scholar 

  40. Rotheray G, Gilbert F. The Natural History of Hoverflies. Forrest Text: Cardigan, UK; 2011.

    Google Scholar 

  41. Reemer M, Ståhls G. Phylogenetic relationships of Microdontinae (Diptera: Syrphidae) based on molecular and morphological characters. Syst Entomol. 2013;38:661–88.

    Article  Google Scholar 

  42. Hippa H, Stahls G. Morphological characters of adult Syrphidae: descriptions and phylogenetic utility. Acta Zool Fenn. 2005;215:1–72.

    Google Scholar 

  43. Mengual X, Ståhls G, Rojo S. Is the mega-diverse genus Ocyptamus (Diptera, Syrphidae) monophyletic? Evidence from molecular characters including the secondary structure of 28S rRNA. Mol Phylogen Evol. 2012;62:191–205.

    Article  Google Scholar 

  44. Larson BM, Kevan PG, Inouye DW. Flies and flowers: taxonomic diversity of anthophiles and pollinators. Can Entomol. 2001;133(4):439–66.

    Article  Google Scholar 

  45. Barrett SC, Helenurm K. The reproductive biology of boreal forest herbs. I. Breeding systems and pollination. Can J Bot. 1987;65(10):2036–46.

    Article  Google Scholar 

  46. Zych M. Pollination biology of Heracleum sphondylium L. (Apiaceae): the advantages of being white and compact. Acta Soc Bot Pol. 2014;71(2):163–70.

    Article  Google Scholar 

  47. Jauker F, Wolters V. Hover flies are efficient pollinators of oilseed rape. Oecologia. 2008;156(4):819–23.

    Article  PubMed  Google Scholar 

  48. Rader R, Howlett BG, Cunningham SA, Westcott DA, Newstrom‐Lloyd LE, Walker MK, Teulon DA, Edwards W. Alternative pollinator taxa are equally efficient but not as effective as the honeybee in a mass flowering crop. J Appl Ecol. 2009;46(5):1080–7.

    Article  Google Scholar 

  49. Ssymank A, Kearns CA, Pape T, Thompson FC. Pollinating flies (Diptera): a major contribution to plant diversity and agricultureal production. Biodiversity. 2008;9(1–2):86–9.

    Article  Google Scholar 

  50. Rader R, Bartomeus I, Garibaldi LA, Garratt MP, Howlett BG, Winfree R, Cunningham SA, Mayfield MM, Arthur AD, Andersson GK. Non-bee insects are important contributors to global crop pollination. Proc Natl Acad Sci U S A. 2016;113(1):146–51.

  51. Reemer M. A review of Microdontinae (Diptera: Syrphidae) of Surinam, with a key to the Neotropical genera. Tijdschr Entomol. 2014;157:27–57.

    Article  Google Scholar 

  52. Heiss EM. A Classification of the larvae and puparia of the Syrphidae of Illinois exclusive of aquatic forms, vol. 187. Urbana: The University of Illinois Press; 1938.

    Google Scholar 

  53. Rotheray GE. Colour Guide to Hoverfly Larvae (Diptera, Syrphidae), vol. 9. Sheffield: Derek Whiteley; 1993.

    Google Scholar 

  54. Reemer M. Review and phylogenetic evaluation of associations between Microdontinae (Diptera: Syrphidae) and ants (Hymenoptera: Formicidae). Psyche. 2013;2013:1–9.

    Article  Google Scholar 

  55. Pérez-Lachaud GJ, Reemer M, Lachaud J-P. An unusual, but not unexpected, evolutionary step taken by syrphid flies: the first record of true primary parasitoidism of ants by Microdontinae. Biol J Linn Soc. 2014;111:462–72.

    Article  Google Scholar 

  56. Mogi M, Chan KL. Predatory habits of dipteran larvae inhabiting Nepenthes pitchers. Raffles Bull Zool. 1996;44(1):233–45.

    Google Scholar 

  57. Rotheray G. The predatory larvae of two Nepenthosyrphus species living in pitcher plants (Diptera, Syrphidae). Stud Dipterol. 2003;10:219–26.

    Google Scholar 

  58. Stuke J-H. Phylogenetische Rekonstruktion der Verwandtschaftsbeziehungen innerhalb der Gattung Cheilosia Meigen, 1822 anhand der Larvenstadien (Diptera: Syrphidae). [Phylogenetic relationships within the genus Cheilosia Meigen, 1822, as evidenced by the larval stages (Diptera: Syrphidae)]. Stud Dipterol. 2000;Supplement 8:1–118.

    Google Scholar 

  59. Morales GE, Wolff M. Insects associated with the composting process of solid urban waste separated at the source. Rev Bras Entomol. 2010;54(4):645–53.

    Article  Google Scholar 

  60. Rojo S, Gilbert FS, Marcos-Garcia MA, Nieto JM, Mier MP. A World Review of Predatory Hoverflies (Diptera, Syrhidae: Syrphinae) and Their Prey. Alicante: CIBIO Ediciones; 2003.

    Google Scholar 

  61. Nishida K, Rotheray G, Thompson FC. First non-predaceous syrphine flower fly (Diptera: Syrphidae): A new leaf-mining Allograpta from Costa Rica. Stud Dipterol. 2002;9(2):421–36.

    Google Scholar 

  62. Weng J, Rotheray G. Another non-predaceous syrphine flower fly (Diptera: Syrphidae): pollen feeding in the larva of Allograpta micrura. Stud Dipterol. 2009;15(2008):245–58.

    Google Scholar 

  63. Marin A, Comerma J. Nota sobre las larvas de Mesograpta polita (Say)(Syrphidae, Diptera) en espigas de maíz (Zea mays L.) en El Limón, Aragua, Venezuela [Note on the larvae of Mesograpta polita (Say)(Syrphidae, Diptera) on the staminate flowers of corn (Zea mays L.) in El Limón, Aragua, Venezuela]. Agronomía tropical. 1969;19:335–9.

    Google Scholar 

  64. Reemer M, Rotheray GE. Pollen feeding larvae in the presumed predatory syrphine genus Toxomerus Macquart (Diptera, Syrphidae). J Nat Hist. 2009;43(15):939–49.

    Article  Google Scholar 

  65. Grosskopf G. Biology and life history of Cheilosia urbana (Meigen) and Cheilosia psilophthalma (Becker), two sympatric hoverflies approved for the biological control of hawkweeds (Hieracium spp.) in New Zealand. Biol Control. 2005;35:142–54.

    Article  Google Scholar 

  66. Grosskopf G, Wilson L, Littlefield J. Host-range investigations of potential biological control agents of alien invasive hawkweeds (Hieracium spp.) in the USA and Canada: an overview. In: Julien MH, Sforza R, Bon MC, Evans CH, Hatcher PE, Hinz HL, Rector BG, CAB International, editors. Proceedings of the XII International Symposium on Biological Control of Weeds. Wallingford: Citeseer; 2008. p. 552–7.

    Google Scholar 

  67. Penney HD, Hassall C, Skevington JH, Abbott KR, Sherratt TN. A comparative analysis of the evolution of imperfect mimicry. Nature. 2012;483:461–6.

    Article  CAS  PubMed  Google Scholar 

  68. Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010;2010(6):t5448.

    Article  Google Scholar 

  69. Niehuis O, Hartig G, Grath S, Pohl H, Lehmann J, Tafer H, Donath A, Krauss V, Eisenhardt C, Hertel J, Petersen M, Mayer C, Meusemann K, Peters Ralph S, Stadler Peter F, Beutel Rolf G, Bornberg-Bauer E, McKenna Duane D, Misof B. Genomic and morphological evidence converge to resolve the enigma of Strepsiptera. Curr Biol. 2012;22(14):1309–13.

    Article  CAS  PubMed  Google Scholar 

  70. Wiegmann BM, Trautwein MD, Winkler IS, Barr NB, Kim J-W, Lambkin C, Bertone MA, Cassel BK, Bayless KM, Heimberg AM, Wheeler BM, Peterson KJ, Pape T, Sinclair BJ, Skevington JH, Blagoderov V, Caravask J, Kutty SN, Schmidt-Ott U, Kampmeier GE, Thompson FC, Grimaldi DA, Beckenbach AT, Courtney GW, Friedrich M, Meier R, Yeates DK. Episodic radiations in the fly tree of life. Proc Natl Acad Sci U S A. 2011;108(14):5690–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Rokyta DR, Lemmon AR, Margres MJ, Aronow K. The venom-gland transcriptome of the eastern diamondback rattlesnake (Crotalus adamanteus). BMC Genomics. 2012;13(1):1.

    Article  Google Scholar 

  73. Lemmon AR, Lemmon EM. A likelihood framework for estimating phylogeographic history on a continuous landscape. Syst Biol. 2008;57(4):544–61.

    Article  PubMed  Google Scholar 

  74. Gibson JF, Skevington JH, Kelso S. Placement of the Conopidae (Diptera) within the Schizophora based on ten mtDNA and nrDNA gene regions. Mol Phylogen Evol. 2010;56:91–103.

    Article  CAS  Google Scholar 

  75. Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 2.73 2010. http://mesquiteproject.org/. Accessed 30 May 2016.

  76. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.

    Article  CAS  PubMed  Google Scholar 

  77. Hennig W. Phylogenetic Systematics. Champaign: University of Illinois Press; 1966.

    Google Scholar 

  78. Griffiths GCD. The phylogenetic classification of Diptera Cyclorrhapha with special reference to the structure of the male postabdomen. The Hague: Dr. W. Junk N. V; 1972.

    Google Scholar 

  79. McAlpine JF. Phylogeny and classification of the Muscomorpha. In: McAlpine JF, editor. Manual of Nearctic Diptera. Vol. 3. Ottawa: Agriculture Canada Monograph No. 32, Research Branch, Agriculture Canada; 1989. p. 1397–505.

    Google Scholar 

  80. Rotheray GE, Gilbert F. Phylogenetic relationships and the larval head of the lower Cyclorrhapha (Diptera). Zool J Linn Soc. 2008;153:287–323.

    Article  Google Scholar 

  81. Tachi T. Homology of the metapleuron of Cyclorrhapha, with discussion of the paraphyly of Syrphoidea (Diptera: Aschiza). Insect Syst Evol. 2014;1–20.

  82. Howard RW, Akre RD, Garnett WB. Chemical mimicry in an obligate predator of carpenter ants (Hymenoptera: Formicidae). Ann Entomol Soc Am. 1990;83(3):607–16.

    Article  CAS  Google Scholar 

  83. Howard RW, Stanley Samuelson DW, Akre RD. Biosynthesis and chemical mimicry of cuticular hydrocarbons from the obligate predator, Microdon albicomatus Novak (Diptera: Syrphidae) and its ant prey, Myrmica incompleta Provancher (Hymenoptera: Formicidae). J Kans Entomol Soc. 1990;63(3):437–43.

    Google Scholar 

  84. Rotheray G, Gilbert F. Phylogeny of Palaearctic Syrphidae (Diptera): evidence from larval stages. Zool J Linn Soc. 1999;127(1):1–112.

    Article  Google Scholar 

  85. Speight MCD. StN key for the identification of the genera of European Syrphidae (Diptera). Syrph the Net, the database of European Syrphidae, vol. 79. Dublin: Syrph the Net publications; 2014.

    Google Scholar 

  86. Hartley JC. A Taxonomic account of the larvae of some British Syrphidae. Proc Zool Soc Lond (Reprinted). 1961;136(Part 4):505–73.

    Google Scholar 

  87. Vockeroth JR. The flower flies of the subfamily Syrphinae of Canada, Alaska and Greenland, vol. 18. Ottawa: Canada Communications Group - Publishing; 1992.

    Google Scholar 

  88. Cheng XY, Lu J, Huang CM, Zhou HZ, Dai ZH, Zhang GX. Determination of phylogenetic position of Pipizini (Diptera : Syrphidae): based on molecular biological and morphological data. Sci China Ser C. 2000;43(2):146–56.

    Article  CAS  Google Scholar 

  89. Vujić AS, Ståhls G, Ačanski J, Bartsch H, Bygebjerg R, Stefanović A. Systematics of Pipizini and taxonomy of European Pipiza Fallén: molecular and morphological evidence (Diptera, Syrphidae). Zool Scr. 2013;42(3):288–305.

    Article  Google Scholar 

  90. Miranda GFG, Skevington JH, Marshall SA, Kelso S. The genus Ocyptamus Macquart (Diptera: Syrphidae): a molecular phylogenetic analysis (in review). Arthropod Syst Phylo. 2016.

  91. Mengual X. The systematic position and phylogenetic relationships of Asiobaccha Violovitsh (Diptera, Syrphidae). J Asia-Pacif Entomol. 2015;18:397–408.

    Article  Google Scholar 

Download references

Acknowledgements

We are grateful for lab assistance provided by Brian K. Cassel, Dept. of Entomology, NC State University. The project could not have been completed without the support of transcriptomes supplied by 1KITE.

Funding

This study was supported by funding to JHS from Agriculture and Agri-Food Canada, and the Natural Sciences and Engineering Research Council of Canada, to ARL and EML by NSF IIP1313554 and to BMW by NSF DEB1257960.

Availability of data and materials

All data used in this project is available as tables in the main text in combination with supplemental tables provided. Concatenated alignments used for phylogenetic analysis have been uploaded to the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra), with accession numbers (Biosample #) available in Table 3.

Authors’ contributions

ADY: manuscript assembly, manuscript revisions, voucher identification. ARL: probe design, lab work oversight, analysis, manuscript revisions. JHS: project conceptualization, manuscript coordination, voucher collection, voucher identification, manuscript revisions, phylogenetic and taxonomic insights. XM: project conceptualization, manuscript revisions, voucher collection, phylogenetic and taxonomic insights. GS: project conceptualization, manuscript revisions, voucher collection, phylogenetic and taxonomic insights. MR: project conceptualization, manuscript revisions, voucher collection, taxonomic insights. KJ: project conceptualization, manuscript revisions, phylogenetic insights. SK: manuscript revisions, analysis, molecular lab work & COI sequencing. EML: manuscript revisions, lab work. MH: project conceptualization, manuscript revisions, taxonomic insights. MDM: project conceptualization, manuscript revisions. BM: manuscript revisions, probe development. BMW: manuscript revisions, probe coordination and development. All authors have read and approved the final version of this manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethical approval and consent to participate

Not applicable.

Animal ethics and client-owner consent

All specimens collected for research purposes comply with the Convention on Trade in Endangered Species of Wild Fauna and Flora and the IUCN Policy on Research Involving Species at Risk.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andrew Donovan Young.

Additional files

Additional file 1: Table S1.

Preliminary loci for enrichment across Insecta. Complete listing of all insect loci used to determine initial probe site selection. (DOCX 165 kb)

Additional file 2: Table S2.

Final Probe Regions. Complete listing of all regions used in alignments to generate probes. (XLSX 1446 kb)

Additional file 3: Table S3.

Final Probe Sequences. Complete listing of all final probes used for Diptera genome amplification. (XLSX 3101 kb)

Additional file 4: Table S4.

Assembly Results. Complete listing of assemblies produced from the raw reads. (XLSX 18 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Young, A.D., Lemmon, A.R., Skevington, J.H. et al. Anchored enrichment dataset for true flies (order Diptera) reveals insights into the phylogeny of flower flies (family Syrphidae). BMC Evol Biol 16, 143 (2016). https://doi.org/10.1186/s12862-016-0714-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-016-0714-0

Keywords