Evolution and phylogenetic utility of CAD (rudimentary) among Mesozoic-aged Eremoneuran Diptera (Insecta)

Culicoides Latreille (Diptera: Ceratopogonidae) taxonomy: current challenges and future directions

Culicoides Latreille biting midges (Diptera: Ceratopogonidae) cause a significant biting nuisance to humans, livestock and equines, and are the biological vectors of a range of internationally important pathogens of both veterinary and medical importance. Despite their economic significance, the delimitation and identification of species and evolutionary relationships between species within this genus remains at best problematic. To date no phylogenetic study has attempted to validate the subgeneric classification of the genus and the monophyly of many of the subgenera remains doubtful. Many informal species groupings are also known to exist but few are adequately described, further complicating accurate identification. Recent contributions to Culicoides taxonomy at the species level have revealed a high correlation between morphological and molecular analyses although molecular analyses are revealing the existence of cryptic species. This review considers the methods for studying the systematics of Culicoides using both morphological and genetic techniques, with a view to understanding the factors limiting our current understanding of Culicoides biology and hence arbovirus epidemiology. In addition, we examine the global status of Culicoides identification, highlighting areas that are poorly addressed, including the potential implementation of emerging technologies.

Armored scale insect endosymbiont diversity at the species level: genealogical patterns of Uzinura diasipipdicola in the Chionaspis pinifoliae-Chionaspis heterophyllae species complex (Hemiptera: Coccoidea: Diaspididae)

Armored scale insects and their primary bacterial endosymbionts show nearly identical patterns of co-diversification when viewed at the family level, though the persistence of these patterns at the species level has not been explored in this group. Therefore we investigated genealogical patterns of co-diversification near the species level between the primary endosymbiont Uzinura diaspidicola and its hosts in the Chionaspis pinifoliae-Chionaspis heterophyllae species complex. To do this we generated DNA sequence data from three endosymbiont loci (rspB, GroEL, and 16S) and analyzed each locus independently using statistical parsimony network analyses and as a concatenated dataset using Bayesian phylogenetic reconstructions. We found that for two endosymbiont loci, 16S and GroEL, sequences from U. diaspidicola were broadly associated with host species designations, while for rspB this pattern was less clear as C. heterophyllae (species S1) shared haplotypes with several other Chionaspis species. We then compared the topological congruence of the phylogenetic reconstructions generated from a concatenated dataset of endosymbiont loci (including all three loci, above) to that from a concatenated dataset of armored scale hosts, using published data from two nuclear loci (28S and EF1α) and one mitochondrial locus (COI-COII) from the armored scale hosts. We calculated whether the two topologies were congruent using the Shimodaira-Hasegawa test. We found no significant differences (P = 0.4892) between the topologies suggesting that, at least at this level of resolution, co-diversification of U. diaspidicola with its armored scale hosts also occurs near the species level. This is the first such study of co-speciation at the species level between U. diaspidicola and a group of armored scale insects.

Diversification in Hawaiian long-legged flies (Diptera: Dolichopodidae: Campsicnemus): biogeographic isolation and ecological adaptation

Flies in the genus Campsicnemus have diversified into the second-largest adaptive radiation of Diptera in the Hawaiian Islands, with 179 Hawaiian endemic species currently described. Here we present the first phylogenetic analysis of Campsicnemus, with a focus on the Hawaiian fauna. We analyzed a combination of two nuclear (CAD, EF1α) and five mitochondrial (COI, COII, 12S, 16S, ND2) loci using Bayesian and maximum likelihood approaches to generate a phylogenetic hypothesis for the genus Campsicnemus. Our sampling included a total of 84 species (6 species from Europe, 1 from North America, 7 species from French Polynesia and 70 species from the Hawaiian Islands). The phylogenies were used to estimate divergence times, reconstruct biogeographic history, and infer ancestral ecological associations within this large genus. We found strong support for a South Pacific+Hawaiian clade, as well as for a monophyletic Hawaiian lineage. Divergence time estimates suggest that Hawaiian Islands were colonized approximately 4.6 million years ago, suggesting that most of the diversity within Campsicnemus evolved since the current high islands began forming ∼5 million years ago. We also observe a novel ecotype within the Pacific Campsicnemus; a widespread obligate water-skating form that has arisen multiple times across the Pacific Islands. Together, these analyses suggest that a combination of ecological, biogeographic and temporal factors have led to the impressive diversity of long-legged flies in Hawaii and elsewhere in the Pacific.

The utility of CAD in recovering Gondwanan vicariance events and the evolutionary history of Aciliini (Coleoptera: Dytiscidae)

BACKGROUND

Aciliini presently includes 69 species of medium-sized water beetles distributed on all continents except Antarctica. The pattern of distribution with several genera confined to different continents of the Southern Hemisphere raises the yet untested hypothesis of a Gondwana vicariance origin. The monophyly of Aciliini has been questioned with regard to Eretini, and there are competing hypotheses about the intergeneric relationship in the tribe. This study is the first comprehensive phylogenetic analysis focused on the tribe Aciliini and it is based on eight gene fragments. The aims of the present study are: 1) to test the monophyly of Aciliini and clarify the position of the tribe Eretini and to resolve the relationship among genera within Aciliini, 2) to calibrate the divergence times within Aciliini and test different biogeographical scenarios, and 3) to evaluate the utility of the gene CAD for phylogenetic analysis in Dytiscidae.

RESULTS

Our analyses confirm monophyly of Aciliini with Eretini as its sister group. Each of six genera which have multiple species are also supported as monophyletic. The origin of the tribe is firmly based in the Southern Hemisphere with the arrangement of Neotropical and Afrotropical taxa as the most basal clades suggesting a Gondwana vicariance origin. However, the uncertainty as to whether a fossil can be used as a stem-or crowngroup calibration point for Acilius influenced the result: as crowngroup calibration, the 95% HPD interval for the basal nodes included the geological age estimate for the Gondwana break-up, but as a stem group calibration the basal nodes were too young. Our study suggests CAD to be the most informative marker between 15 and 50 Ma. Notably, the 2000 bp CAD fragment analyzed alone fully resolved the tree with high support.

CONCLUSIONS

1) Molecular data confirmed Aciliini as a monophyletic group. 2) Bayesian optimizations of the biogeographical history are consistent with an influence of Gondwana break-up history, but were dependent on the calibration method. 3) The evaluation using a method of phylogenetic signal per base pair indicated Wnt and CAD as the most informative of our sampled genes.

Phylogeny, historical biogeography and divergence time estimates of the genus Colophon Gray (Coleoptera : Lucanidae)

The flightless Cape High-mountain stag beetle genus Colophon (Coleoptera: Lucanidae) is studied. All species are endemic to the Cape Floristic Region of the Western Cape, South Africa. The study aimed to determine the specific and phylogenetic status of the described species of Colophon and to determine the main factors driving their evolution by testing the hypothesis of a lowland origin. This was achieved by analysing the mitochondrial COI, 16S rRNA and nuclear CAD genes, using a maximum likelihood and Bayesian approach. Timing of key biogeographical events in the diversification of Colophon was estimated in BEAST. The combined molecular dataset supports the described species. High genetic divergence was found between the described taxa. The mean estimated divergence of the genus was mid-Cretaceous, with a split into two lineages during the early Paleocene to mid Eocene. Species divergence was shown to have occurred during mid to late Miocene. Most species of Colophon showed an allopatric distribution, although contact zones between geographically adjacent species are likely. Climate seems to be the main driving factor behind Colophon evolution and the hypothesis of a lowland origin appears to be supported.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Molecular phylogeny and biogeography of the Hawaiian craneflies Dicranomyia (Diptera: Limoniidae)

The Hawaiian Diptera offer an opportunity to compare patterns of diversification across large and small endemic radiations with varying species richness and levels of single island endemism. The craneflies (Limoniidae: Dicranomyia) represent a small radiation of 13 described species that have diversified within the Hawaiian Islands. We used Bayesian and maximum likelihood approaches to generate a molecular phylogeny of the Hawaiian Dicranomyia using a combination of nuclear and mitochondrial loci, estimated divergence times and reconstructed ancestral ranges. Divergence time estimation and ancestral range reconstruction suggest that the colonization that led to most of the diversity within the craneflies arrived prior to the formation of Kauai and demonstrates that the two major clades within that radiation contrast sharply in their patterns of diversification.

Environmental monitoring using next generation sequencing: rapid identification of macroinvertebrate bioindicator species

Introduction

Invertebrate communities are central to many environmental monitoring programs. In freshwater ecosystems, aquatic macroinvertebrates are collected, identified and then used to infer ecosystem condition. Yet the key step of species identification is often not taken, as it requires a high level of taxonomic expertise, which is lacking in most organizations, or species cannot be identified as they are morphologically cryptic or represent little known groups. Identifying species using DNA sequences can overcome many of these issues; with the power of next generation sequencing (NGS), using DNA sequences for routine monitoring becomes feasible.

Results

In this study, we test if NGS can be used to identify species from field-collected samples in an important bioindicator group, the Chironomidae. We show that Cytochrome oxidase I (COI) and Cytochrome B (CytB) sequences provide accurate DNA barcodes for chironomid species. We then develop a NGS analysis pipeline to identifying species using megablast searches of high quality sequences generated using 454 pyrosequencing against comprehensive reference libraries of Sanger-sequenced voucher specimens. We find that 454 generated COI sequences successfully identified up to 96% of species in samples, but this increased up to 99% when combined with CytB sequences. Accurate identification depends on having at least five sequences for a species; below this level species not expected in samples were detected. Incorrect incorporation of some multiplex identifiers (MID’s) used to tag samples was a likely cause, and most errors could be detected when using MID tags on forward and reverse primers. We also found a strong quantitative relationship between the number of 454 sequences and individuals showing that it may be possible to estimate the abundance of species from 454 pyrosequencing data.

Conclusions

Next generation sequencing using two genes was successful for identifying chironomid species. However, when detecting species from 454 pyrosequencing data sets it was critical to include known individuals for quality control and to establish thresholds for detecting species. The NGS approach developed here can lead to routine species-level diagnostic monitoring of aquatic ecosystems.

The hymenopteran tree of life: evidence from protein-coding genes and objectively aligned ribosomal data

Previous molecular analyses of higher hymenopteran relationships have largely been based on subjectively aligned ribosomal sequences (18S and 28S). Here, we reanalyze the 18S and 28S data (unaligned about 4.4 kb) using an objective and a semi-objective alignment approach, based on MAFFT and BAli-Phy, respectively. Furthermore, we present the first analyses of a substantial protein-coding data set (4.6 kb from one mitochondrial and four nuclear genes). Our results indicate that previous studies may have suffered from inflated support values due to subjective alignment of the ribosomal sequences, but apparently not from significant biases. The protein data provide independent confirmation of several earlier results, including the monophyly of non-xyelid hymenopterans, Pamphilioidea + Unicalcarida, Unicalcarida, Vespina, Apocrita, Proctotrupomorpha and core Proctotrupomorpha. The protein data confirm that Aculeata are nested within a paraphyletic Evaniomorpha, but cast doubt on the monophyly of Evanioidea. Combining the available morphological, ribosomal and protein-coding data, we examine the total-evidence signal as well as congruence and conflict among the three data sources. Despite an emerging consensus on many higher-level hymenopteran relationships, several problems remain unresolved or contentious, including rooting of the hymenopteran tree, relationships of the woodwasps, placement of Stephanoidea and Ceraphronoidea, and the sister group of Aculeata.

Molecular phylogeny of the burying beetles (Coleoptera: Silphidae: Nicrophorinae)

Burying beetles (Silphidae: Nicrophorus) are well-known for their monopolization of small vertebrate carcasses in subterranean crypts and complex biparental care behaviors. They have been the focus of intense behavioral, ecological, and conservation research since the 1980s yet no thorough phylogenetic estimate for the group exists. Herein, we infer relationships, test past hypotheses of relationships, and test biogeographic scenarios among 55 of the subfamily Nicrophorinae's currently valid and extant 72 species. Two mitochondrial genes, COI and COII, and two nuclear genes, the D2 region of 28S, and the protein coding gene CAD, provided 3,971 nucleotides for 58 nicrophorine and 5 outgroup specimens. Ten partitions, with each modeled by GTR+I+G, were used for a 100 M generation MrBayes analysis and maximum likelihood bootstrapping with Garli. The inferred Bayesian phylogeny was mostly well-resolved with only three weak branches of biogeographic relevance. The common ancestor of the subfamily and of the genus Nicrophorus was reconstructed as Old World with four separate transitions to the New World and four reverse colonizations of the Old World from the New. Divergence dating from analysis with BEAST indicate the genus Nicrophorus originated in the Cretaceous, 127-99 Ma. Most prior, pre-cladistic hypotheses of relationships were strongly rejected while most modern hypotheses were largely congruent with monophyletic groups in our estimated phylogeny. Our results reject a recent hypothesis that Nicrophorus morio Gebler, 1817 (NEW STATUS as valid species) is a subspecies of N. germanicus (L., 1758). Two subgenera of Nicrophorus are recognized: NecroxenusSemenov-Tian-Shanskij, 1933, and NicrophorusFabricius, 1775.

Molecular systematics of the Calliphoridae (Diptera: Oestroidea): evidence from one mitochondrial and three nuclear genes

Approximately 8% of calyptrate species diversity comes from the Calliphoridae, which includes flies of medical, veterinary, and forensic importance. The status of family Calliphoridae has for years been the central systematic problem of the superfamily Oestroidea, and phylogenetic relationships between the key groups of the Calliphoridae are unresolved and controversial. We reconstructed phylogenies of the Calliphoridae within the larger context of the other Oestroidea based on 5,189 bp of combined data from one mitochondrial (cytochrome oxidase subunit one) and three nuclear (carbamoylphosphate synthetase, elongation factor one alpha, and 28S ribosomal RNA) genes using maximum parsimony, maximum likelihood, and Bayesian methods. Trees obtained from the different phylogenetic methods were almost identical. Calliphoridae is polyphyletic, with the phylogenetic position of Mesembrinellinae still uncertain but clearly outside the lineage that includes other Calliphoridae and some noncalliphorids, and Polleniinae is the sister group of the family Tachinidae. Strong support for a sister group relationship between Rhiniinae and traditional calliphorid subfamilies conflicts with a recent proposal to give Rhiniinae family status. All calliphorid subfamilies (except Calliphorinae) for which we had more than one species were monophyletic. Melanomyinae was nested within Calliphorinae. Toxotarsinae was more closely related to Calliphorinae rather than, as indicated by morphology, to Chrysomyinae. Efforts to resolve the relationships of the Oestroid families were largely inconclusive, although the monophyly of the superfamily was strongly supported.

A molecular phylogeny for the leaf-roller moths (Lepidoptera: Tortricidae) and its implications for classification and life history evolution

BACKGROUND

Tortricidae, one of the largest families of microlepidopterans, comprise about 10,000 described species worldwide, including important pests, biological control agents and experimental models. Understanding of tortricid phylogeny, the basis for a predictive classification, is currently provisional. We present the first detailed molecular estimate of relationships across the tribes and subfamilies of Tortricidae, assess its concordance with previous morphological evidence, and re-examine postulated evolutionary trends in host plant use and biogeography.

METHODOLOGY/PRINCIPAL FINDINGS

We sequenced up to five nuclear genes (6,633 bp) in each of 52 tortricids spanning all three subfamilies and 19 of the 22 tribes, plus up to 14 additional genes, for a total of 14,826 bp, in 29 of those taxa plus all 14 outgroup taxa. Maximum likelihood analyses yield trees that, within Tortricidae, differ little among data sets and character treatments and are nearly always strongly supported at all levels of divergence. Support for several nodes was greatly increased by the additional 14 genes sequenced in just 29 of 52 tortricids, with no evidence of phylogenetic artifacts from deliberately incomplete gene sampling. There is strong support for the monophyly of Tortricinae and of Olethreutinae, and for grouping of these to the exclusion of Chlidanotinae. Relationships among tribes are robustly resolved in Tortricinae and mostly so in Olethreutinae. Feeding habit (internal versus external) is strongly conserved on the phylogeny. Within Tortricinae, a clade characterized by eggs being deposited in large clusters, in contrast to singly or in small batches, has markedly elevated incidence of polyphagous species. The five earliest-branching tortricid lineages are all species-poor tribes with mainly southern/tropical distributions, consistent with a hypothesized Gondwanan origin for the family.

CONCLUSIONS/SIGNIFICANCE

We present the first robustly supported phylogeny for Tortricidae, and a revised classification in which all of the sampled tribes are now monophyletic.

Phylogeny of Bembidion and related ground beetles (Coleoptera: Carabidae: Trechinae: Bembidiini: Bembidiina)

The phylogeny of the large genus Bembidion and related genera is inferred from four nuclear protein-coding genes (CAD, wingless, arginine kinase, and topoisomerase I), ribosomal DNA (28S and 18S), and the mitochondrial gene cytochrome oxidase I (COI). 230 of the more than 1200 species of Bembidion are sampled, as well as 26 species of five related genera, and 14 outgroups. Nuclear copies (numts) of COI were found sparsely scattered through sampled species. The resulting phylogeny, based upon individual gene analyses and combined analyses using maximum likelihood and parsimony, is very well supported at most nodes. Additional analyses explored the evidence, and corroborate the phylogeny. Seven analyses, each with one of the seven genes removed from the combined matrix, were also conducted, and yielded maximum likelihood bootstrap trees sharing over 92% of their nodes with the original, well-resolved bootstrap trees based on the complete set of seven genes. All key nodes were present in all seven analyses missing a single gene, indicating that support for these nodes comes from at least two genes. In addition, the inferred maximum likelihood tree based on the combined matrix is well-behaved and self-predicting, in that simulated evolution of sequences on the inferred tree under the inferred model of evolution yields a matrix from which all but one of the model tree's clades are recovered with bootstrap value >50, suggesting that internal branches in the tree may be of a length to yield sequences sufficient to allow their inference. All likelihood analyses were conducted under both a proportion-invariable plus gamma site-to-site rate variation model, as well as a simpler gamma model. The choice of model did not have a major effect on inferred phylogenies or their bootstrap values. The inferred phylogeny shows that Bembidarenas is not closely related to Bembidiina, and Phrypeus is likely distant as well; the remaining genera of Bembidiina form a monophyletic group. Lionepha, formerly considered a subgenus of Bembidion, is shown to be outside of the clade of Asaphidion+Bembidion, and is separated as its own genus. B. (Phyla) obtusum is quite isolated within Bembidion, and there is some evidence that the remaining Bembidion form a clade. Within Bembidion, there are three large clades that are well-supported, the Bembidion, Odontium, and Ocydromus Series. The Bembidion Series contains Bembidion (s. str.), Notaphus, Furcacampa, Emphanes, Trepanedoris, Diplocampa, and related Holarctic species; all species from South America, Australia, New Zealand; and most species from southern Africa and Madagascar. All species in South America, except for members of Notaphus and Nothocys, form a clade, the Antiperyphanes Complex, which has independently radiated into body forms and niches occupied by multiple, independent Northern-Hemisphere forms. All species from New Zealand, including Zecillenus, and Australian species formerly placed in Ananotaphus together form a clade. Bembidion (s. str.) and Cyclolopha are in a clade with the Old World, Southern Hemisphere lineages Notaphocampa, Sloanephila, and Omotaphus. The large subgenus Notaphus appears to have originated in South America, with all Northern Hemisphere Notaphus arising from within a south-temperate grade. All major variation in frontal furrows on the head is contained within the Bembidion Series. The Odontium Series contains subgenera Hirmoplataphus and Hydriomicrus, which together are the sister clade of Odontium, Bracteon, Ochthedromus, Pseudoperyphus, and Microserrullula. The very large Ocydromus Series, dominant in the Holarctic region, includes the Ocydromus Complex, with many subgenera, including Hypsipezum and Leuchydrium; the phylogeny within this group is notably at odds with the current classification. Also included in the Ocydromus Series are Nepha and Bembidionetolitzkya, as well as the Princidium Complex, in which the intertidal B. (Cillenus) laterale falls. Outside these three series are a number of smaller groups, including the Plataphus Complex (containing Blepharoplataphus, Plataphus, the latter including Plataphodes); the Hydrium Complex (Metallina, Chlorodium, and Hydrium, which contains Eurytrachelus), whose sister group might be subgenus Andrewesa; Trechonepha and Liocosmius, which might be sisters; and B. (Melomalus) planatum, which is not close to Plataphus. There is some evidence that these groups plus the Ocydromus and Odontium Series form a clade. A few enigmatic groups were harder to place. The sister group of the pair Philochthus plus Philochthemphanes might be B. wickhami; Eupetedromus is well outside the three major series and not related to Notaphus; the high-elevation Asian group Hoquedela is a very isolated lineage. Notaphiellus is removed from synonymy with Nothocys, and placed in synonymy with Notaphus; Plataphodes is synonymized with Plataphus, as Plataphus is paraphyletic otherwise; Eurytrachelus is synonymized with Hydrium. A new subgenus, Lindrochthus, is described to house the distinctive B. wickhami. The implications of the inferred phylogeny for some morphological characters used in Bembidiina systematics are explored, and some of the most widely used (e.g., location of discal seta ed3 on the elytron, and shape of the shoulder) are shown to be notably homoplastic. For example, the location of elytral seta ed3 has undergone at least nine transitions between two states.

Can deliberately incomplete gene sample augmentation improve a phylogeny estimate for the advanced moths and butterflies (Hexapoda: Lepidoptera)?

This paper addresses the question of whether one can economically improve the robustness of a molecular phylogeny estimate by increasing gene sampling in only a subset of taxa, without having the analysis invalidated by artifacts arising from large blocks of missing data. Our case study stems from an ongoing effort to resolve poorly understood deeper relationships in the large clade Ditrysia ( > 150,000 species) of the insect order Lepidoptera (butterflies and moths). Seeking to remedy the overall weak support for deeper divergences in an initial study based on five nuclear genes (6.6 kb) in 123 exemplars, we nearly tripled the total gene sample (to 26 genes, 18.4 kb) but only in a third (41) of the taxa. The resulting partially augmented data matrix (45% intentionally missing data) consistently increased bootstrap support for groupings previously identified in the five-gene (nearly) complete matrix, while introducing no contradictory groupings of the kind that missing data have been predicted to produce. Our results add to growing evidence that data sets differing substantially in gene and taxon sampling can often be safely and profitably combined. The strongest overall support for nodes above the family level came from including all nucleotide changes, while partitioning sites into sets undergoing mostly nonsynonymous versus mostly synonymous change. In contrast, support for the deepest node for which any persuasive molecular evidence has yet emerged (78-85% bootstrap) was weak or nonexistent unless synonymous change was entirely excluded, a result plausibly attributed to compositional heterogeneity. This node (Gelechioidea + Apoditrysia), tentatively proposed by previous authors on the basis of four morphological synapomorphies, is the first major subset of ditrysian superfamilies to receive strong statistical support in any phylogenetic study. A "more-genes-only" data set (41 taxa×26 genes) also gave strong signal for a second deep grouping (Macrolepidoptera) that was obscured, but not strongly contradicted, in more taxon-rich analyses.

A combination of molecular and morphological approaches resolves species in the taxonomically difficult genus Procladius Skuse (Diptera: Chironomidae) despite high intra-specific morphological variation

Molecular approaches for identifying aquatic macroinvertebrate species are increasingly being used but there is ongoing debate about the number of DNA markers needed to differentiate species accurately. Here, we use two mitochondrial genes (cytochrome oxidase I, cytochrome b) and a nuclear gene (carbamoylphosphate synthetase) to differentiate species variation within the taxonomically challenging chironomid genus Procladius from southern Australia, a genus which is important for pollution monitoring. The mitochondrial genes indicated cryptic species that were subsequently linked to morphological variation at the larval and pupal stage. Two species previously described based on morphological criteria were linked to molecular markers, and there was evidence for additional cryptic species. Each genetic marker provided different information, highlighting the importance of considering multiple genes when dissecting taxonomically difficult groups, particularly those used in pollution monitoring.

Increased gene sampling strengthens support for higher-level groups within leaf-mining moths and relatives (Lepidoptera: Gracillariidae)

BACKGROUND

Researchers conducting molecular phylogenetic studies are frequently faced with the decision of what to do when weak branch support is obtained for key nodes of importance. As one solution, the researcher may choose to sequence additional orthologous genes of appropriate evolutionary rate for the taxa in the study. However, generating large, complete data matrices can become increasingly difficult as the number of characters increases. A few empirical studies have shown that augmenting genes even for a subset of taxa can improve branch support. However, because each study differs in the number of characters and taxa, there is still a need for additional studies that examine whether incomplete sampling designs are likely to aid at increasing deep node resolution. We target Gracillariidae, a Cretaceous-age (~100 Ma) group of leaf-mining moths to test whether the strategy of adding genes for a subset of taxa can improve branch support for deep nodes. We initially sequenced ten genes (8,418 bp) for 57 taxa that represent the major lineages of Gracillariidae plus outgroups. After finding that many deep divergences remained weakly supported, we sequenced eleven additional genes (6,375 bp) for a 27-taxon subset. We then compared results from different data sets to assess whether one sampling design can be favored over another. The concatenated data set comprising all genes and all taxa and three other data sets of different taxon and gene sub-sampling design were analyzed with maximum likelihood. Each data set was subject to five different models and partitioning schemes of non-synonymous and synonymous changes. Statistical significance of non-monophyly was examined with the Approximately Unbiased (AU) test.

RESULTS

Partial augmentation of genes led to high support for deep divergences, especially when non-synonymous changes were analyzed alone. Increasing the number of taxa without an increase in number of characters led to lower bootstrap support; increasing the number of characters without increasing the number of taxa generally increased bootstrap support. More than three-quarters of nodes were supported with bootstrap values greater than 80% when all taxa and genes were combined. Gracillariidae, Lithocolletinae + Leucanthiza, and Acrocercops and Parectopa groups were strongly supported in nearly every analysis. Gracillaria group was well supported in some analyses, but less so in others. We find strong evidence for the exclusion of Douglasiidae from Gracillarioidea sensu Davis and Robinson (1998). Our results strongly support the monophyly of a G.B.R.Y. clade, a group comprised of Gracillariidae + Bucculatricidae + Roeslerstammiidae + Yponomeutidae, when analyzed with non-synonymous changes only, but this group was frequently split when synonymous and non-synonymous substitutions were analyzed together.

CONCLUSIONS

1) Partially or fully augmenting a data set with more characters increased bootstrap support for particular deep nodes, and this increase was dramatic when non-synonymous changes were analyzed alone. Thus, the addition of sites that have low levels of saturation and compositional heterogeneity can greatly improve results. 2) Gracillarioidea, as defined by Davis and Robinson (1998), clearly do not include Douglasiidae, and changes to current classification will be required. 3) Gracillariidae were monophyletic in all analyses conducted, and nearly all species can be placed into one of six strongly supported clades though relationships among these remain unclear. 4) The difficulty in determining the phylogenetic placement of Bucculatricidae is probably attributable to compositional heterogeneity at the third codon position. From our tests for compositional heterogeneity and strong bootstrap values obtained when synonymous changes are excluded, we tentatively conclude that Bucculatricidae is closely related to Gracillariidae + Roeslerstammiidae + Yponomeutidae.

Molecular phylogeny of the blowfly genus Chrysomya

Chrysomya Robineau-Desvoidy (Diptera: Calliphoridae) is a genus of blowfly commonly observed in tropical and subtropical countries of the Old World. Species in this genus are vectors of bacteria, protozoans and helminths, cause myiasis, are predators of other carrion insects, and are important forensic indicators. Hypotheses concerning the evolution of sex determination, larval anatomy and genome size in Chrysomya have been difficult to evaluate because a robust phylogeny of the genus was lacking. Similarly, the monophyly of subgenera was uncertain. The phylogeny of Chrysomya spp. was reconstructed based on 2386 bp of combined mitochondrial cytochrome oxidase subunit I (COI) and nuclear carbamoylphosphate synthetase (CPS) genes. Maximum parsimony (MP), maximum likelihood (ML) and Bayesian analysis (BA) differed only slightly in the resulting tree topology. Chrysomya was monophyletic. Monogenic reproduction is almost certainly derived rather than, as has been suggested, primitive within the genus, and tuberculate larvae probably evolved twice. Genome size is more likely to have decreased over evolutionary time rather than, as has been suggested, increased within the genus, but its correlation with developmental time was not observed. The subgenera Microcalliphora, Eucompsomyia and Achoetandrus were recovered as monophyletic.

Overcoming the effects of rogue taxa: Evolutionary relationships of the bee flies

Bombyliidae (5000 sp.), or bee flies, are a lower brachyceran family of flower-visiting flies that, as larvae, act as parasitoids of other insects. The evolutionary relationships are known from a morphological analysis that yielded minimal support for higher-level groupings. We use the protein-coding gene CAD and 28S rDNA to determine phylogeny and to test the monophyly of existing subfamilies, the divisions Tomophtalmae, and 'the sand chamber subfamilies'. Additionally, we demonstrate that consensus networks can be used to identify rogue taxa in a Bayesian framework. Pruning rogue taxa post-analysis from the final tree distribution results in increased posterior probabilities. We find 8 subfamilies to be monophyletic and the subfamilies Heterotropinae and Mythicomyiinae to be the earliest diverging lineages. The large subfamily Bombyliinae is found to be polyphyletic and our data does not provide evidence for the monophyly of Tomophthalmae or the 'sand chamber subfamilies'.

Systematics and biogeography of the Gondwanan Orthocladiinae (Diptera: Chironomidae)

Restrictions to effective dispersal and gene flow caused by the fragmentation of ancient supercontinents are considered to have driven diversification and speciation on disjunct landmasses globally. Investigating the role that these processes have played in the development of diversity within and among taxa is crucial to understanding the origins and evolution of regional biotas. Within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae), a group of genera that are distributed across the austral continents (Australia, New Zealand, South America) have been proposed to represent a relict Gondwanan clade. We used a molecular approach to resolve relationships among taxa with the aim to determine the relative roles that vicariance and dispersal may have played in the evolution of this group. Continental biotas did not form monophyletic groups, in accordance with expectations given existing morphological evidence. Patterns of phylogenetic relationships among taxa did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent. Likewise, divergence time estimates, particularly for New Zealand taxa, largely post-dated continental fragmentation and implied instead that several transoceanic dispersal events may have occurred post-vicariance. Passive dispersal of gravid female chironomid adults is the most likely mechanism for transoceanic movement, potentially facilitated by West Wind Drift or anti-cyclone fronts. Estimated timings of divergence among Australian and South American Botryocladius, on the other hand, were congruent with the proposed ages of separation of the two continents from Antarctica. Taken together, these data suggest that a complex relationship between both vicariance and dispersal may explain the evolution of this group. The sampling regime we implemented here was the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly revealed several novel taxa that will require formal description.

Episodic radiations in the fly tree of life

Flies are one of four superradiations of insects (along with beetles, wasps, and moths) that account for the majority of animal life on Earth. Diptera includes species known for their ubiquity (Musca domestica house fly), their role as pests (Anopheles gambiae malaria mosquito), and their value as model organisms across the biological sciences (Drosophila melanogaster). A resolved phylogeny for flies provides a framework for genomic, developmental, and evolutionary studies by facilitating comparisons across model organisms, yet recent research has suggested that fly relationships have been obscured by multiple episodes of rapid diversification. We provide a phylogenomic estimate of fly relationships based on molecules and morphology from 149 of 157 families, including 30 kb from 14 nuclear loci and complete mitochondrial genomes combined with 371 morphological characters. Multiple analyses show support for traditional groups (Brachycera, Cyclorrhapha, and Schizophora) and corroborate contentious findings, such as the anomalous Deuterophlebiidae as the sister group to all remaining Diptera. Our findings reveal that the closest relatives of the Drosophilidae are highly modified parasites (including the wingless Braulidae) of bees and other insects. Furthermore, we use micro-RNAs to resolve a node with implications for the evolution of embryonic development in Diptera. We demonstrate that flies experienced three episodes of rapid radiation--lower Diptera (220 Ma), lower Brachycera (180 Ma), and Schizophora (65 Ma)--and a number of life history transitions to hematophagy, phytophagy, and parasitism in the history of fly evolution over 260 million y.

Phylogeny of the African ball-rolling dung beetle genus Epirinus Reiche (Coleoptera:Scarabaeidae:Scarabaeinae)

Recent phylogenetic analyses have pointed to the dung beetle genus Epirinus as the putative African ancestral roller group. Consequently, we tested the roller status of species in the genus with observational studies and constructed a molecular phylogeny based on partial sequences of two mitochondrial and two nuclear genes for 16 of 29 species. Tested species were confirmed to be dung rollers. Monophyly of the genus was confirmed, lending support to the synonymy of the wingless genus Endroedyantus with Epirinus. Moreover, the phylogenetic hypothesis was found to have a similar topology with a previously published one based on morphological data. A combined molecular and morphology analysis showed congruence between the molecular and morphological datasets. The loss of flight in some species and estimated divergence dates within Epirinus are discussed.

Evolutionary relationships among primary endosymbionts of the mealybug subfamily phenacoccinae (hemiptera: Coccoidea: Pseudococcidae)

Mealybugs (Coccoidea: Pseudococcidae) are sap-sucking plant parasites that harbor bacterial endosymbionts within specialized organs. Previous studies have identified two subfamilies, Pseudococcinae and Phenacoccinae, within mealybugs and determined the primary endosymbionts (P-endosymbionts) of the Pseudococcinae to be Betaproteobacteria ("Candidatus Tremblaya princeps") containing Gammaproteobacteria secondary symbionts. Here, the P-endosymbionts of phenacoccine mealybugs are characterized based on 16S rRNA from the bacteria of 20 species of phenacoccine mealybugs and four outgroup Puto species (Coccoidea: Putoidae) and aligned to more than 100 published 16S rRNA sequences from symbiotic and free-living bacteria. Phylogenetic analyses recovered three separate lineages of bacteria from the Phenacoccinae, and these are considered to be the P-endosymbionts of their respective mealybug hosts, with those from (i) the mealybug genus Rastrococcus belonging to the Bacteroidetes, (ii) the subterranean mealybugs, tribe Rhizoecini, also within Bacteroidetes, in a clade sister to cockroach endosymbionts (Blattabacterium), and (iii) the remaining Phenacoccinae within the Betaproteobacteria, forming a well-supported sister group to "Candidatus Tremblaya princeps." Names are proposed for two strongly supported lineages: "Candidatus Brownia rhizoecola" for P-endosymbionts of Rhizoecini and "Candidatus Tremblaya phenacola" for P-endosymbionts of Phenacoccinae excluding Rastrococcus and Rhizoecini. Rates of nucleotide substitution among lineages of Tremblaya were inferred to be significantly faster than those of free-living Betaproteobacteria. Analyses also recovered a clade of Gammaproteobacteria, sister to the P-endosymbiont lineage of aphids ("Candidatus Buchnera aphidicola"), containing the endosymbionts of Putoidae, the secondary endosymbionts of pseudococcine mealybugs, and the endosymbionts of several other insect groups.

A multigene phylogeny of the fly superfamily Asiloidea (Insecta): Taxon sampling and additional genes reveal the sister-group to all higher flies (Cyclorrhapha)

Asiloidea are a group of 9 lower brachyceran fly families, considered to be the closest relative to the large Metazoan radiation Eremoneura (Cyclorrhapha+Empidoidea). The evidence for asiloid monophyly is limited, and few characters define the relationships between the families of Asiloidea and Eremoneura. Additionally, enigmatic genera, Hilarimorpha and Apystomyia, retain morphological characters of both asiloids and higher flies. We use the nuclear protein-coding gene CAD and 28S rDNA to test the monophyly of Asiloidea and to resolve its relationship to Eremoneura. We explore the effects of taxon sampling on support values and topological stability, the resolving power of additional genes, and hypothesis testing using four-cluster likelihood mapping. We find that: (1) the 'asiloid' genus Apystomyia is sister to Cyclorrhapha, (2) the remaining asiloids are monophyletic at the exclusion of the family Bombyliidae, and (3) our best estimate of relationships places the asiloid flies excluding Bombyliidae as the sister-group to Eremoneura, though high support is lacking.

Did dung beetles arise in Africa? A phylogenetic hypothesis based on five gene regions

Scarabaeinae dung beetle phylogenetics are poorly understood, although recent phylogenetic reconstruction based on morphology and molecular analyses are congruent on the hypothesis that the oldest Gondwana tribes are polyphyletic and that representatives of the ancestral groups are found in Africa. We present a molecular phylogeny of the African representatives of the two oldest tribes, Canthonini (the putative ancestor of all "rollers") and Dichotomiini (thought to represent the ancestral "tunnellers"), based on partial sequence data from two mitochondrial and three nuclear genes, extracted from eight of the nine dichotomiine genera and 17 of the 23 genera of Canthonini. Three well-supported lineages were consistently obtained. Divergence times estimated the origin of the two tribes at around 56 million years ago (MYA) with the splits of the three dung beetle lineages being estimated to have taken place between 40 and 34MYA. The ages of these splits and subsequent radiation of the modern dung beetle groups concur with those predicted by the fossil record and coincide with the proposed age of radiation of the large mammal groups with whose dung most African dung beetles are still associated. Dispersal of dung beetle groups from Africa is proposed as a biogeographic model, and evidence is presented that dung beetles disperse quickly and widely across continents, and even oceans.

Placement of Conopidae (Diptera) within Schizophora based on mtDNA and nrDNA gene regions

The first attempt to phylogenetically place Conopidae using molecular characters, as well as the largest molecular analysis of relationships within Schizophora (Diptera) to date, is presented. Twenty-eight taxa from 11 acalyptrate families and seven acalyptrate superfamilies are represented. Nearly 12,800 bp of sequence data from 10 genes representing both mitochondrial (cytochrome oxidase I (COI), cytochrome b (cytB), and 12S) and nuclear genes (28S, the carbamoyl phosphate synthetase region of CAD (CAD), elongation factor-1alpha (EF-1alpha), white, alanyl-tRNA synthetase (AATS), triose phosphate isomerase (TPI), and phosphogluconate dehydrogenase (PGD)) are analysed. Parsimony and Bayesian analyses strongly support the monophyly of both Conopidae and Schizophora. While in the parsimony analysis, Conopidae are placed as sister to the remaining Schizophora, the Bayesian analysis recovers a Conopidae+Lauxaniidae clade. The value of nuclear, mitochondrial, ribosomal, and protein-coding gene sequence data for answering phylogenetic questions at different levels of divergence is evaluated.

Hennig's orphans revisited: testing morphological hypotheses in the "Opomyzoidea" (Diptera: Schizophora)

The acalyptrate fly superfamily Opomyzoidea, as currently recognized, is a poorly-known group of 14 families. The composition of this group and relationships among included families have been controversial. Furthermore, the delimitation of two opomyzoid families, Aulacigastridae and Periscelididae, has been unstable with respect to placement of the genera Stenomicra, Cyamops, and Planinasus. To test the monophyly of Opomyzoidea, previously proposed relationships between families, and the position of the three problematic genera, we sequenced over 3300bp of nucleotide sequence data from the 28S ribosomal DNA and CAD (rudimentary) genes from 29 taxa representing all opomyzoid families, as well as 13 outgroup taxa. Relationships recovered differed between analyses, and only branches supporting well-established monophyletic families were recovered with high support, with a few exceptions. Opomyzoidea and its included subgroup, Asteioinea, were found to be non-monophyletic. Stenomicra, Cyamops, and Planinasus group consistently with Aulacigastridae, contrary to recent classifications. Xenasteiidae and Australimyzidae, two small, monogeneric families placed in separate superfamilies, were strongly supported as sister groups.

Molecular evolution of the pDo500 satellite DNA family in Dolichopoda cave crickets (Rhaphidophoridae)

BACKGROUND

Non-coding satellite DNA (satDNA) usually has a high turn-over rate frequently leading to species specific patterns. However, some satDNA families evolve more slowly and can be found in several related species. Here, we analyzed the mode of evolution of the pDo500 satDNA family of Dolichopoda cave crickets. In addition, we discuss the potential of slowly evolving satDNAs as phylogenetic markers.

RESULTS

We sequenced 199 genomic or PCR amplified satDNA repeats of the pDo500 family from 12 Dolichopoda species. For the 38 populations under study, 39 pDo500 consensus sequences were deduced. Phylogenetic analyses using Bayesian, Maximum Parsimony, and Maximum Likelihood approaches yielded largely congruent tree topologies. The vast majority of pDo500 sequences grouped according to species designation. Scatter plots and statistical tests revealed a significant correlation between genetic distances for satDNA and mitochondrial DNA. Sliding window analyses showed species specific patterns of variable and conserved regions. The evolutionary rate of the pDo500 satDNA was estimated to be 1.63-1.78% per lineage per million years.

CONCLUSIONS

The pDo500 satDNA evolves gradually at a rate that is only slightly faster than previously published rates of insect mitochondrial COI sequences. The pDo500 phylogeny was basically congruent with the previously published mtDNA phylogenies. Accordingly, the slowly evolving pDo500 satDNA family is indeed informative as a phylogenetic marker.

Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology

Background

Phylogenetic analyses provide a framework for examining the evolution of morphological and molecular diversity, interpreting patterns in biogeography, and achieving a stable classification. The generic and suprageneric relationships within mosquitoes (Diptera: Culicidae) are poorly resolved, making these subjects difficult to address.

Results

We carried out maximum parsimony and maximum likelihood, including Bayesian, analyses on a data set consisting of six nuclear genes and 80 morphological characters to assess their ability to resolve relationships among 25 genera. We also estimated divergence times based on sequence data and fossil calibration points, using Bayesian relaxed clock methods. Strong support was recovered for the basal position and monophyly of the subfamily Anophelinae and the tribes Aedini and Sabethini of subfamily Culicinae. Divergence times for major culicid lineages date to the early Cretaceous.

Conclusions

Deeper relationships within the family remain poorly resolved, suggesting the need for additional taxonomic sampling. Our results support the notion of rapid radiations early in the diversification of mosquitoes.

Toward reconstructing the evolution of advanced moths and butterflies (Lepidoptera: Ditrysia): an initial molecular study

BACKGROUND

In the mega-diverse insect order Lepidoptera (butterflies and moths; 165,000 described species), deeper relationships are little understood within the clade Ditrysia, to which 98% of the species belong. To begin addressing this problem, we tested the ability of five protein-coding nuclear genes (6.7 kb total), and character subsets therein, to resolve relationships among 123 species representing 27 (of 33) superfamilies and 55 (of 100) families of Ditrysia under maximum likelihood analysis.

RESULTS

Our trees show broad concordance with previous morphological hypotheses of ditrysian phylogeny, although most relationships among superfamilies are weakly supported. There are also notable surprises, such as a consistently closer relationship of Pyraloidea than of butterflies to most Macrolepidoptera. Monophyly is significantly rejected by one or more character sets for the putative clades Macrolepidoptera as currently defined (P < 0.05) and Macrolepidoptera excluding Noctuoidea and Bombycoidea sensu lato (P < or = 0.005), and nearly so for the superfamily Drepanoidea as currently defined (P < 0.08). Superfamilies are typically recovered or nearly so, but usually without strong support. Relationships within superfamilies and families, however, are often robustly resolved. We provide some of the first strong molecular evidence on deeper splits within Pyraloidea, Tortricoidea, Geometroidea, Noctuoidea and others.Separate analyses of mostly synonymous versus non-synonymous character sets revealed notable differences (though not strong conflict), including a marked influence of compositional heterogeneity on apparent signal in the third codon position (nt3). As available model partitioning methods cannot correct for this variation, we assessed overall phylogeny resolution through separate examination of trees from each character set. Exploration of "tree space" with GARLI, using grid computing, showed that hundreds of searches are typically needed to find the best-feasible phylogeny estimate for these data.

CONCLUSION

Our results (a) corroborate the broad outlines of the current working phylogenetic hypothesis for Ditrysia, (b) demonstrate that some prominent features of that hypothesis, including the position of the butterflies, need revision, and (c) resolve the majority of family and subfamily relationships within superfamilies as thus far sampled. Much further gene and taxon sampling will be needed, however, to strongly resolve individual deeper nodes.

Testing the monophyly of Calocidae (Insecta: Trichoptera) based on multiple molecular data

Calocidae constitute a hypothesised monophyletic group of caddisflies (Trichoptera) being geographically restricted to New Zealand (one genus) and Australia (five genera). This analysis tests the monophyly of the family based on sequences from five different molecular genes. The complete data set includes 29 species and covers a complete genus representation of the Calocidae as well as representatives of other families in which one or more calocid genera have been classified. Sequences from two mitochondrial (cytochrome oxidase I and 16S) and three nuclear (elongation factor 1-alpha, RNA polymerase-II, and Cadherin) genes were used, resulting in a 3958bp data set and 37.1% parsimony informative characters. The Cadherin (CAD) and RNA polymerase-II (POL-II) genes are used for the first time for revealing Trichoptera phylogenies. The character matrix was analyzed by using maximum parsimony (MP) and Bayesian criteria, the latter by applying three different partition strategies for comparison. Two most parsimonious trees were found, differing in the position of one clade within the sister-group to a monophyletic Calocidae. The Bayesian tree based on the maximum number of partitions differs from trees based on a reduced partition analysis with respect to taxa outside the current circumscription of Calocidae. Both the MP and Bayesian analyses left Calocidae monophyletic, with a monophyletic clade of all Australian genera being sister-group to the New Zealand genus. The results from the agreement subtree analysis demonstrates that CAD performs well both separately and in combination with other genes and adds substantial resolution to the calocid phylogeny in a combined MP analysis.

Single-copy nuclear genes resolve the phylogeny of the holometabolous insects

BACKGROUND

Evolutionary relationships among the 11 extant orders of insects that undergo complete metamorphosis, called Holometabola, remain either unresolved or contentious, but are extremely important as a context for accurate comparative biology of insect model organisms. The most phylogenetically enigmatic holometabolan insects are Strepsiptera or twisted wing parasites, whose evolutionary relationship to any other insect order is unconfirmed. They have been controversially proposed as the closest relatives of the flies, based on rDNA, and a possible homeotic transformation in the common ancestor of both groups that would make the reduced forewings of Strepsiptera homologous to the reduced hindwings of Diptera. Here we present evidence from nucleotide sequences of six single-copy nuclear protein coding genes used to reconstruct phylogenetic relationships and estimate evolutionary divergence times for all holometabolan orders.

RESULTS

Our results strongly support Hymenoptera as the earliest branching holometabolan lineage, the monophyly of the extant orders, including the fleas, and traditionally recognized groupings of Neuropteroidea and Mecopterida. Most significantly, we find strong support for a close relationship between Coleoptera (beetles) and Strepsiptera, a previously proposed, but analytically controversial relationship. Exploratory analyses reveal that this relationship cannot be explained by long-branch attraction or other systematic biases. Bayesian divergence times analysis, with reference to specific fossil constraints, places the origin of Holometabola in the Carboniferous (355 Ma), a date significantly older than previous paleontological and morphological phylogenetic reconstructions. The origin and diversification of most extant insect orders began in the Triassic, but flourished in the Jurassic, with multiple adaptive radiations producing the astounding diversity of insect species for which these groups are so well known.

CONCLUSION

These findings provide the most complete evolutionary framework for future comparative studies on holometabolous model organisms and contribute strong evidence for the resolution of the 'Strepsiptera problem', a long-standing and hotly debated issue in insect phylogenetics.

Phylogeny and biogeography of hawkmoths (Lepidoptera: Sphingidae): evidence from five nuclear genes

BACKGROUND

The 1400 species of hawkmoths (Lepidoptera: Sphingidae) comprise one of most conspicuous and well-studied groups of insects, and provide model systems for diverse biological disciplines. However, a robust phylogenetic framework for the family is currently lacking. Morphology is unable to confidently determine relationships among most groups. As a major step toward understanding relationships of this model group, we have undertaken the first large-scale molecular phylogenetic analysis of hawkmoths representing all subfamilies, tribes and subtribes.

METHODOLOGY/PRINCIPAL FINDINGS

The data set consisted of 131 sphingid species and 6793 bp of sequence from five protein-coding nuclear genes. Maximum likelihood and parsimony analyses provided strong support for more than two-thirds of all nodes, including strong signal for or against nearly all of the fifteen current subfamily, tribal and sub-tribal groupings. Monophyly was strongly supported for some of these, including Macroglossinae, Sphinginae, Acherontiini, Ambulycini, Philampelini, Choerocampina, and Hemarina. Other groupings proved para- or polyphyletic, and will need significant redefinition; these include Smerinthinae, Smerinthini, Sphingini, Sphingulini, Dilophonotini, Dilophonotina, Macroglossini, and Macroglossina. The basal divergence, strongly supported, is between Macroglossinae and Smerinthinae+Sphinginae. All genes contribute significantly to the signal from the combined data set, and there is little conflict between genes. Ancestral state reconstruction reveals multiple separate origins of New World and Old World radiations.

CONCLUSIONS/SIGNIFICANCE

Our study provides the first comprehensive phylogeny of one of the most conspicuous and well-studied insects. The molecular phylogeny challenges current concepts of Sphingidae based on morphology, and provides a foundation for a new classification. While there are multiple independent origins of New World and Old World radiations, we conclude that broad-scale geographic distribution in hawkmoths is more phylogenetically conserved than previously postulated.

Evidence for a composite anterior determinant in the hover fly Episyrphus balteatus (Syrphidae), a cyclorrhaphan fly with an anterodorsal serosa anlage

Most insect embryos develop from a monolayer of cells around the yolk, but only part of this blastoderm forms the embryonic rudiment. Another part forms extra-embryonic serosa. Size and position of the serosa anlage vary between species, and previous work raises the issue of whether such differences co-evolve with the mechanisms that establish anteroposterior (AP) polarity of the embryo. AP polarity of the Drosophila embryo depends on bicoid, which is necessary and sufficient to determine the anterior body plan. Orthologs of bicoid have been identified in various cyclorrhaphan flies and their occurrence seems to correlate with a mid-dorsal serosa or amnioserosa anlage. Here, we introduce with Episyrphus balteatus (Syrphidae) a cyclorrhaphan model for embryonic AP axis specification that features an anterodorsal serosa anlage. Current phylogenies place Episyrphus within the clade that uses bicoid mRNA as anterior determinant, but no bicoid-like sequence could be identified in this species. Using RNA interference (RNAi) and ectopic mRNA injection, we obtained evidence that pattern formation along the entire AP axis of the Episyrphus embryo relies heavily on the precise regulation of caudal, and that anterior pattern formation in particular depends on two localized factors rather than one. Early zygotic activation of orthodenticle is separated from anterior repression of caudal, two distinct functions which in Drosophila are performed jointly by bicoid, whereas hunchback appears to be regulated by both factors. Furthermore, we found that overexpression of orthodenticle is sufficient to confine the serosa anlage of Episyrphus to dorsal blastoderm. We discuss our findings in a phylogenetic context and propose that Episyrphus employs a primitive cyclorrhaphan mechanism of AP axis specification.

Monophyly of terrestrial adephagan beetles as indicated by three nuclear genes (Coleoptera: Carabidae and Trachypachidae)

The beetle suborder Adephaga is traditionally divided into two sections on the basis of habitat, terrestrial Geadephaga and aquatic Hydradephaga. Monophyly of both groups is uncertain, and the relationship of the two groups has implications for inferring habitat transitions within Adephaga. Here we examine phylogenetic relationships of these groups using evidence provided by DNA sequences from all four suborders of beetles, including 60 species of Adephaga, four Archostemata, three Myxophaga, and ten Polyphaga. We studied 18S ribosomal DNA and 28S ribosomal DNA, aligned with consideration of secondary structure, as well as the nuclear protein-coding gene wingless. Independent and combined Bayesian, likelihood, and parsimony analyses of all three genes supported placement of Trachypachidae in a monophyletic Geadephaga, although for analyses of 28S rDNA and some parsimony analyses only if Coleoptera is constrained to be monophyletic. Most analyses showed limited support for the monophyly of Hydradephaga. Outside of Adephaga, there is support from the ribosomal genes for a sister group relationship between Adephaga and Polyphaga. Within the small number of sampled Polyphaga, analyses of 18S rDNA, wingless, and the combined matrix supports monophyly of Polyphaga exclusive of Scirtoidea. Unconstrained analyses of the evolution of habitat suggest that Adephaga was ancestrally aquatic with one transition to terrestrial. However, in analyses constrained to disallow changes from aquatic to terrestrial habitat, the phylogenies imply two origins of aquatic habit within Adephaga.

Evaluating nuclear protein-coding genes for phylogenetic utility in beetles

Although nuclear protein-coding genes have proven broadly useful for phylogenetic inference, relatively few such genes are regularly employed in studies of Coleoptera, the most diverse insect order. We increase the number of loci available for beetle systematics by developing protocols for three genes previously unused in beetles (alpha-spectrin, RNA polymerase II and topoisomerase I) and by refining protocols for five genes already in use (arginine kinase, CAD, enolase, PEPCK and wingless). We evaluate the phylogenetic performance of each gene in a Bayesian framework against a presumably known test phylogeny. The test phylogeny covers 31 beetle specimens and two outgroup taxa of varying age, including three of the four extant beetle suborders and a denser sampling in Adephaga and in the carabid genus Bembidion. All eight genes perform well for Cenozoic divergences and accurately separate closely related species within Bembidion, but individual genes differ markedly in accuracy over the older Mesozoic and Permian divergences. The concatenated data reconstruct the test phylogeny with high support in both Bayesian and parsimony analyses, indicating that combining data from multiple nuclear loci will be a fruitful approach for assembling the beetle tree of life.

A phylogeny of the oil bee tribe Ctenoplectrini (Hymenoptera: Anthophila) based on mitochondrial and nuclear data: evidence for early Eocene divergence and repeated out-of-Africa dispersal

The bee tribe Ctenoplectrini, with two genera, comprises nine species in tropical Africa and ten in Asia and Australia. Most of them collect floral oil, pollen, and nectar from Cucurbitaceae, but three species are thought to be cleptoparasites. The unusual morphology of Ctenoplectrini has made it difficult to infer their closest relatives, in turn preventing an understanding of these bees' geographic and temporal origin. We used two mitochondrial and two nuclear markers (4741 nucleotides) generated for most of the species to test the monophyly of the tribe, its relationships to other Apidae, and its biogeographic history. Ctenoplectrini are strongly supported as monophyletic and closest to the Long-horned bees, Eucerini. The presumably cleptoparasitic species form a clade (Ctenoplectrina) that is sister to the remaining species (Ctenoplectra), confirming the independent evolution of cleptoparasitism in this tribe. Tree topology and molecular dating together suggest that Ctenoplectrini originated in Africa in the Early Eocene and that Ctenoplectra dispersed twice from Africa to Asia, sometime in the Late Eocene, 30-40 my ago, from where one species reached the Australian continent via Indonesia and New Guinea in the mid-Miocene, c. 13 my ago. Dry and cool mid-Miocene climates also coincide with the divergence between Ctenoplectra bequaerti from West Africa and Ctenoplectra terminalis from East and South Africa, perhaps related to fragmentation of the equatorial African rainforest belt.

A phylogeny of long-tongued horse flies (Diptera:Tabanidae:Philoliche) with the first cladistic review of higher relationships within the family

Long-tongued horse flies (Diptera : Tabanidae : Pangoniinae) have proboscis lengths at least as long as their heads, the longest belonging to the Old World genus Philoliche (Wiedemann, 1820). These long proboscides are used to probe for nectar in deep-throated flowers. For some flower species, these flies are the only known pollinators. Although horse flies are both vectors of disease and important pollinators, there has been no previous study of phylogenetic relationships below the subfamily level. The first comprehensive phylogenetic hypothesis for Philoliche (and consequently, Tabanidae) based on molecular data from one nuclear (CAD) and one mitochondrial (COI) gene is presented here. Using an exemplar approach, species from eight of nine tribes in all three subfamilies in Tabanidae were sampled, giving some of the first insights into relationships among the family as a whole. All nine subgenera of Philoliche, and multiple geographic representatives of the subgenus P. (Philoliche) (Wiedemann, 1820) in southern Africa were also sampled. Within the subgenus Philoliche, molecular support was found for a previously synonymised species. In general, these analyses sustain earlier intuitive classifications, but do not support the monophyly of all currently recognised subfamilies.

Genetic Analysis of Fungicide-Resistant Sclerotinia homoeocarpa Isolates from Tennessee and Northern Mississippi

Sclerotinia homoeocarpa is the causal agent of dollar spot disease that reduces the uniformity and aesthetic value of golf putting greens. Fungicide-resistant isolates of S. homoeocarpa were collected from putting greens at 10 locations across Tennessee and northern Mississippi. Genetic diversity among the 60 isolates was investigated using vegetative compatibility, conserved gene sequences, and amplified fragment length polymorphism (AFLP). Six tester strains were paired with Tennessee and northern Mississippi isolates on potato dextrose agar. Some of the 60 isolates were delineated into vegetative compatibility groups, but fungicide resistance could not be associated with a particular vegetative compatibility group. Genetic similarities of isolates at the vegetative compatibility level could be attributed to founder effects. Sequencing the regions of CAD, EF1-α, β-tubulin, and internal transcribed spacers revealed 100% homology among isolates. Capillary gel electrophoresis and analysis of AFLP fragments indicated 86 to 100% similarity between the isolates. Vegetative compatibility and molecular data indicate that the populations of the pathogen are clonal. Isolates did not cluster according to fungicide resistance during unweighted pair group with arithmetic means analysis, but did appear to cluster according to vegetative compatibility group and location. Although associations could not be made between molecular markers and fungicide resistance, links between vegetative compatibility and AFLP markers may provide a foundation from which other studies could be performed.

Phylogeny of pteromalid parasitic wasps (Hymenoptera: Pteromalidae): Initial evidence from four protein-coding nuclear genes

Chalcidoidea (approximately 22,000 described species) is the most ecologically diverse superfamily of parasitic Hymenoptera and plays a major role in the biological control of insect pests. However, phylogenetic relationships both within and between chalcidoid families have been poorly understood, particularly for the large family Pteromalidae and relatives. Forty-two taxa, broadly representing Chalcidoidea but concentrated in the 'pteromalid lineage,' were sequenced for 4620 bp of protein-coding sequence from four nuclear genes for which we present new primers. These are: CAD (1719 bp) DDC (708 bp), enolase (1149 bp), and PEPCK (1044 bp). The combined data set was analyzed using parsimony, maximum likelihood, and Bayesian methods. Statistical significance of the apparent non-monophyly of some taxonomic groups on our trees was evaluated using the approximately unbiased test of Shimodaira [Shimodaira, H. 2002. An approximately unbiased test of phylogenetic tree selection. Syst. Biol. 51(3), 492-508]. In accord with previous studies, we find moderate to strong support for monophyly of Chalcidoidea, a sister-group relationship of Mymaridae to the remainder of Chalcidoidea, and a relatively basal placement of Encarsia (Aphelinidae) within the latter. The 'pteromalid lineage' of families is generally recovered as monophyletic, but the hypothesis of monophyly for Pteromalidae, which appear paraphyletic with respect to all other families sampled in that lineage, is decisively rejected (P < 10(-14)). Within Pteromalidae, monophyly was strongly supported for nearly all tribes represented by multiple exemplars, and for two subfamilies. All other multiply-represented subfamilies appeared para- or polyphyletic in our trees, although monophyly was significantly rejected only for Miscogasterinae, Ormocerinae, and Colotrechninae. The limited resolution obtained in the analyses presented here reinforces the idea that reconstruction of pteromalid phylogeny is a difficult problem, possibly due to rapid radiation of many chalcidoid taxa. Initial phylogenetic comparisons of life history traits suggest that the ancestral chalcidoid was small-bodied and parasitized insect eggs.

Phylogeny and Bayesian divergence time estimations of small-headed flies (Diptera: Acroceridae) using multiple molecular markers

The first formal analysis of phylogenetic relationships among small-headed flies (Acroceridae) is presented based on DNA sequence data from two ribosomal (16S and 28S) and two protein-encoding genes: carbomoylphosphate synthase (CPS) domain of CAD (i.e., rudimentary locus) and cytochrome oxidase I (COI). DNA sequences from 40 species in 22 genera of Acroceridae (representing all three subfamilies) were compared with outgroup exemplars from Nemestrinidae, Stratiomyidae, Tabanidae, and Xylophagidae. Parsimony and Bayesian simultaneous analyses of the full data set recover a well-resolved and strongly supported hypothesis of phylogenetic relationships for major lineages within the family. Molecular evidence supports the monophyly of traditionally recognised subfamilies Philopotinae and Panopinae, but Acrocerinae are polyphyletic. Panopinae, sometimes considered "primitive" based on morphology and host-use, are always placed in a more derived position in the current study. Furthermore, these data support emerging morphological evidence that the type genus Acrocera Meigen, and its sister genus Sphaerops, are atypical acrocerids, comprising a sister lineage to all other Acroceridae. Based on the phylogeny generated in the simultaneous analysis, historical divergence times were estimated using Bayesian methodology constrained with fossil data. These estimates indicate Acroceridae likely evolved during the late Triassic but did not diversify greatly until the Cretaceous.

Low nuclear DNA variation supports a recent origin of Hawaiian Hylaeus bees (Hymenoptera: Colletidae)

Previous phylogenetic work on the Hawaiian bees of the genus Hylaeus, based on mitochondrial DNA and morphology, appeared to support a recent origin for the group, but support for the resulting tree was weak. Four nuclear genes with varying evolutionary rates -- arginine kinase, EF-1alpha, opsin, and wingless -- were sequenced for a reduced taxon set in an attempt to find one or more data set that would provide better support. All showed very low variation (<2%) in the ingroup. Comparison among genes revealed a much higher than expected rate of evolution in mtDNA, especially at first and second positions. While the data from the nuclear genes showed insufficient variation for phylogenetic analysis, the strong sequence similarity among the Hawaiian species supports the previous hypothesis of a recent origin for the group.