Nuclear genes resolve mesozoic-aged divergences in the insect order Lepidoptera

BmPMFBP1 regulates the development of eupyrene sperm in the silkworm, Bombyx mori

Sperm deliver the male complement of DNA to the ovum, and thus play a key role in sexual reproduction. Accordingly, spermatogenesis has outstanding significance in fields as disparate as infertility treatments and pest-control, making it a broadly interesting and important focus for molecular genetics research in a wide range of species. Here we investigate spermatogenesis in the model lepidopteran insect Bombyx mori (silkworm moth), with particular focus on the gene PMFBP1 (polyamine modulated factor 1 binding protein 1). In humans and mouse, PMFBP1 is essential for spermatogenesis, and mutations of this gene are associated with acephalic spermatozoa, which cause infertility. We identified a B. mori gene labeled as “PMFBP1” in GenBank’s RefSeq database and sought to assess its role in spermatogenesis. Like in mammals, the silkworm version of this gene (BmPMFBP1) is specifically expressed in testes. We subsequently generated BmPMFBP1 mutants using a transgenic CRISPR/Cas9 system. Mutant males were sterile while the fertility of mutant females was comparable to wildtype females. In B. mori, spermatogenesis yields two types of sperm, the nucleated fertile eupyrene sperm, and anucleated unfertile apyrene sperm. Mutant males produced abnormal eupyrene sperm bundles but normal apyrene sperm bundles. For eupyrene sperm, nuclei were mislocated and disordered inside the bundles. We also found the BmPMFBP1 deficiency blocked the release of eupyrene sperm bundles from testes to ejaculatory seminalis. We found no obvious abnormalities in the production of apyrene sperm in mutant males, and double-matings with apyrene-deficient sex-lethal mutants rescued the ΔBmPMFBP1 infertility phenotype. These results indicate BmPMFBP1 functions only in eupyrene spermatogenesis, and highlight that distinct genes underlie the development of the two sperm morphs commonly found in Lepidoptera. Bioinformatic analyses suggest PMFBP1 may have evolved independently in lepidoptera and mammals, and that despite the shared name, are likely not homologous genes.

The presence of nucleated and anucleated dimorphic sperm produced by a single male is a notable characteristic of lepidopteran insects. Previously we identified the gene BmSxl and BmPnldc1 are required for apyrene and eupyrene sperm development, respectively. However, there remains very little known about the molecular mechanism of eupyrene and apyrene sperm development and function. In human and mouse, the protein PMFBP1 is related to acephalic spermatozoa syndrome. Here, we generate somatic mutants in the silkworm Bombyx mori for this gene. In the silkworm, BmPMFBP1 is essential for male fertility. Loss of BmPMFBP1 function results in defective eupyrene sperm bundles in the early elongation stage of spermatogenesis. The nuclei of the eupyrene sperm bundles were displaced from sperm heads while the apyrene sperm bundles were normal. This deficiency also results in the failure of the release of eupyrene sperm bundles from testes to the ejaculatory seminalis duct. Our study proves the important function of BmPMFBP1 in the development of the eupyrene sperm in the silkworm and identifies a potential target gene for lepidopteran pest control.

The Evolution of Endophagy in Herbivorous Insects

Herbivorous feeding inside plant tissues, or endophagy, is a common lifestyle across Insecta, and occurs in insect taxa that bore, roll, tie, mine, gall, or otherwise modify plant tissues so that the tissues surround the insects while they are feeding. Some researchers have developed hypotheses to explain the adaptive significance of certain endophytic lifestyles (e.g., miners or gallers), but we are unaware of previous efforts to broadly characterize the adaptive significance of endophagy more generally. To fill this knowledge gap, we characterized the limited set of evolutionary selection pressures that could have encouraged phytophagous insects to feed inside plants, and then consider how these factors align with evidence for endophagy in the evolutionary history of orders of herbivorous insects. Reviewing the occurrence of endophytic taxa of various feeding guilds reveals that the pattern of evolution of endophagy varies strongly among insect orders, in some cases being an ancestral trait (e.g., Coleoptera and Lepidoptera) while being more derived in others (e.g., Diptera). Despite the large diversity of endophagous lifestyles and evolutionary trajectories that have led to endophagy in insects, our consideration of selection pressures leads us to hypothesize that nutritionally based factors may have had a stronger influence on evolution of endophagy than other factors, but that competition, water conservation, and natural enemies may have played significant roles in the development of endophagy.

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths

Lepidoptera play key roles in many biological systems. Butterflies are hypothesized to have evolved contemporaneously with flowering plants, and moths are thought to have gained anti-bat defenses in response to echolocating predatory bats, but these hypotheses have largely gone untested. Using a transcriptomic, dated evolutionary tree of Lepidoptera, we demonstrate that the most recent common ancestor of Lepidoptera is considerably older than previously hypothesized. The oldest moths in crown Lepidoptera were present in the Carboniferous, some 300 million years ago, and began to diversify largely in synchrony with angiosperms. We show that multiple lineages of moths independently evolved hearing organs well before the origin of bats, rejecting the hypothesis that lepidopteran hearing organs arose in response to these predators.

Butterflies and moths (Lepidoptera) are one of the major superradiations of insects, comprising nearly 160,000 described extant species. As herbivores, pollinators, and prey, Lepidoptera play a fundamental role in almost every terrestrial ecosystem. Lepidoptera are also indicators of environmental change and serve as models for research on mimicry and genetics. They have been central to the development of coevolutionary hypotheses, such as butterflies with flowering plants and moths’ evolutionary arms race with echolocating bats. However, these hypotheses have not been rigorously tested, because a robust lepidopteran phylogeny and timing of evolutionary novelties are lacking. To address these issues, we inferred a comprehensive phylogeny of Lepidoptera, using the largest dataset assembled for the order (2,098 orthologous protein-coding genes from transcriptomes of 186 species, representing nearly all superfamilies), and dated it with carefully evaluated synapomorphy-based fossils. The oldest members of the Lepidoptera crown group appeared in the Late Carboniferous (∼300 Ma) and fed on nonvascular land plants. Lepidoptera evolved the tube-like proboscis in the Middle Triassic (∼241 Ma), which allowed them to acquire nectar from flowering plants. This morphological innovation, along with other traits, likely promoted the extraordinary diversification of superfamily-level lepidopteran crown groups. The ancestor of butterflies was likely nocturnal, and our results indicate that butterflies became day-flying in the Late Cretaceous (∼98 Ma). Moth hearing organs arose multiple times before the evolutionary arms race between moths and bats, perhaps initially detecting a wide range of sound frequencies before being co-opted to specifically detect bat sonar. Our study provides an essential framework for future comparative studies on butterfly and moth evolution.

Phylogeny and Evolution of Lepidoptera

Until recently, deep-level phylogeny in Lepidoptera, the largest single radiation of plant-feeding insects, was very poorly understood. Over the past two decades, building on a preceding era of morphological cladistic studies, molecular data have yielded robust initial estimates of relationships both within and among the ∼43 superfamilies, with unsolved problems now yielding to much larger data sets from high-throughput sequencing. Here we summarize progress on lepidopteran phylogeny since 1975, emphasizing the superfamily level, and discuss some resulting advances in our understanding of lepidopteran evolution.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

Evolution of larval life mode of Oecophoridae (Lepidoptera: Gelechioidea) inferred from molecular phylogeny

Phylogenetic relationships within family Oecophoridae have been poorly understood. Consequently the subfamily and genus level classifications with this family problematic. A comprehensive phylogenetic analysis of Oecophoridae, the concealer moths, was performed based on analysis of 4444 base pairs of mitochondrial COI, nuclear ribosomal RNA genes (18S and 28S) and nuclear protein coding genes (IDH, MDH, Rps5, EF1a and wingless) for 82 taxa. Data were analyzed using maximum likelihood (ML), parsimony (MP) and Bayesian (BP) phylogenetic frameworks. Phylogenetic analyses indicated that (i) genera Casmara, Tyrolimnas and Pseudodoxia did not belong to Oecophoridae, suggesting that Oecophoridae s. authors was not monophyletic; (ii) other oecophorids comprising two subfamilies, Pleurotinae and Oecophorinae, were nested within the same clade, and (iii) Martyringa, Acryptolechia and Periacmini were clustered with core Xyloryctidae. They appeared to be sister lineage with core Oecophoridae. BayesTraits were implemented to explore the ancestral character states to infer historical microhabitat patterns and sheltering strategy of larvae. Reconstruction of ancestral microhabitat of oecophorids indicated that oecophorids might have evolved from dried plant feeders and further convergently specialized. The ancestral larva sheltering strategy of oecophorids might have used a silk tube by making itself, shifting from mining leaves.

Intragenomic variations of multicopy ITS2 marker in Agrodiaetus blue butterflies (Lepidoptera, Lycaenidae)

The eukaryotic ribosomal DNA cluster consists of multiple copies of three genes, 18S, 5. 8S and 28S rRNAs, separated by multiple copies of two internal transcribed spacers, ITS1 and ITS2. It is an important, frequently used marker in both molecular cytogenetic and molecular phylogenetic studies. Despite this, little is known about intragenomic variations within the copies of eukaryotic ribosomal DNA genes and spacers. Here we present data on intraindividual variations of ITS2 spacer in three species of Agrodiaetus Hübner, 1822 blue butterflies revealed by cloning technique. We demonstrate that a distinctly different intragenomic ITS2 pattern exists for every individual analysed. ITS2 sequences of these species show significant intragenomic variation (up to 3.68% divergence), setting them apart from each other on inferred phylogenetic tree. This variation is enough to obscure phylogenetic relationships at the species level.

Say goodbye to tribes in the new house fly classification: A new molecular phylogenetic analysis and an updated biogeographical narrative for the Muscidae (Diptera)

House flies are one of the best known groups of flies and comprise about 5000 species worldwide. Despite over a century of intensive taxonomic research on these flies, classification of the Muscidae is still poorly resolved. Here we brought together the most diverse molecular dataset ever examined for the Muscidae, with 142 species in 67 genera representing all tribes and all biogeographic regions. Four protein coding genes were analyzed: mitochondrial CO1 and nuclear AATS, CAD (region 4) and EF1-α. Maximum likelihood and Bayesian approaches were used to analyze five different partitioning schemes for the alignment. We also used Bayes factors to test monophyly of the traditionally accepted tribes and subfamilies. Most subfamilial taxa were not recovered in our analyses, and accordingly monophyly was rejected by Bayes factor tests. Our analysis consistently found three main clades of Muscidae and so we propose a new classification with only three subfamilies without tribes. Additionally, we provide the first timeframe for the diversification of all major lineages of house flies and examine contemporary biogeographic hypotheses in light of this timeframe. We conclude that the muscid radiation began in the Paleocene to Eocene and is congruent with the final stages of the breakup of Gondwana, which resulted in the complete separation of Antarctica, Australia, and South America. With this newly proposed classification and better understanding of the timing of evolutionary events, we provide new perspectives for integrating morphological and ecological evolutionary understanding of house flies, their taxonomy, phylogeny, and biogeography.

Molecular characterization and phylogenetic analysis of the eukaryotic translation initiation factor 4A gene in Antheraea pernyi (Lepdoptera: Saturniidae)

Eukaryotic initiation factor 4A (eIF-4A) is an essential component for protein translation in eukaryotes. The eIF-4A gene (ApeIF-4A) was isolated and characterized from Antheraea pernyi (Guérin-Méneville) (Lepidoptera: Saturniidae). The obtained cDNA sequence was 1,435-bp long with an open reading frame of 1,266 bp encoding 421 amino acids. The predicted amino acid sequence shared several conserved features as found in known eIF-4As and revealed 74 and 78% identities with eIF-4As of Homo sapiens L. and Drosophila melanogaster (Meigen), respectively. Reverse transcription-polymerase chain reaction (RT-PCR) analysis showed that ApeIF-4A was transcribed at four developmental stages and in all tissues tested, suggesting that it plays an important role in development of A. pernyi. Homologous alignment suggested that eIF-4As are highly conserved throughout evolution of eukaryote organisms. Phylogenetic trees based on the amino acid and nucleotide sequences of eIF-4A demonstrated a similar topology with the classical systematics, suggesting that it has the potential value in phylogenetic inference of eukaryotes.

Allopatric distribution and diversification without niche shift in a bryophyte-feeding basal moth lineage (Lepidoptera: Micropterigidae)

The Lepidoptera represent one of the most successful radiations of plant-feeding insects, which predominantly took place within angiosperms beginning in the Cretaceous period. Angiosperm colonization is thought to underlie the evolutionary success of the Lepidoptera because angiosperms provide an enormous range of niches for ecological speciation to take place. By contrast, the basal lepidopteran lineage, Micropterigidae, remained unassociated with angiosperms since Jurassic times but nevertheless achieved a modest diversity in the Japanese Archipelago. We explored the causes and processes of diversification of the Japanese micropterigid moths by performing molecular phylogenetic analysis and extensive ecological surveying. Phylogenetic analysis recovered a monophyletic group of approximately 25 East Asian endemic species that feed exclusively on the liverwort Conocephalum conicum, suggesting that niche shifts hardly played a role in their diversification. Consistent with the low flying ability of micropterigid moths, the distributions of the Conocephalum specialists are each localized and allopatric, indicating that speciation by geographical isolation has been the major process shaping the diversity of Japanese Micropterigidae. To our knowledge, this is the largest radiation of herbivorous insects that does not accompany any apparent niche differentiation. We suggest that the significance of non-ecological speciation during the diversification of the Lepidoptera is commonly underestimated.

Molecular data sheds light on the classification of long-legged flies (Diptera:Dolichopodidae)

Dolichopodidae (long-legged flies) is the world’s fourth largest dipteran family, but a phylogeny based on a broad global taxon sample is still lacking. We present here a first molecular phylogenetic hypothesis for Dolichopodidae, based on 157 dolichopodid species in 68 genera and 15 subfamilies from the Old and New World, and seven empidoid species (Empididae, Hybotidae) as outgroups. Both relatively fast-evolving mitochondrial markers (COI, 12S, 16S) and a more conserved nuclear marker (18S) were used, the latter being widely employed to study the phylogeny at higher taxonomic levels. We present strong evidence for Microphorinae as sister group to Dolichopodidae sensu stricto, and for the monophyletic Parathalassiinae as part of Dolichopodidae sensu stricto. Monophyly of Achalcinae, Dolichopodinae, and Sciapodinae is supported and Stolidosomatinae are placed within Sympycninae. Diaphorinae, Medeterinae, Neurigoninae, Rhaphiinae, and Sympycninae are paraphyletic, and Hydrophorinae and Peloropeodinae polyphyletic. Our broad taxon sample allows us to gain new insights into the complex systematics of Dolichopodidae. Our results highlight several problems with the traditional classification, which have considerable consequences for the systematic status of some taxa. The poor resolution observed in deep divergences supports previous hypotheses suggesting a rapid early radiation of Dolichopodidae.

Molecular evolution of lepidopteran silk proteins: insights from the ghost moth, Hepialus californicus

Silk production has independently evolved in numerous arthropod lineages, such as Lepidoptera, the moths and butterflies. Lepidopteran larvae (caterpillars) synthesize silk proteins in modified salivary glands and spin silk fibers into protective tunnels, escape lines, and pupation cocoons. Molecular sequence data for these proteins are necessary to determine critical features of their function and evolution. To this end, we constructed an expression library from the silk glands of the ghost moth, Hepialus californicus, and characterized light chain fibroin and heavy chain fibroin gene transcripts. The predicted H. californicus silk fibroins share many elements with other lepidopteran and trichopteran fibroins, such as conserved placements of cysteine, aromatic, and polar amino acid residues. Further comparative analyses were performed to determine site-specific signatures of selection and to assess whether fibroin genes are informative as phylogenetic markers. We found that purifying selection has constrained mutation within the fibroins and that light chain fibroin is a promising molecular marker. Thus, by characterizing the H. californicus fibroins, we identified key functional amino acids and gained insight into the evolutionary processes that have shaped these adaptive molecules.

Assessing the value of DNA barcodes and other priority gene regions for molecular phylogenetics of Lepidoptera

BACKGROUND

Despite apparently abundant amounts of observable variation and species diversity, the order Lepidoptera exhibits a morphological homogeneity that has provided only a limited number of taxonomic characters and led to widespread use of nucleotides for inferring relationships. This study aims to characterize and develop methods to quantify the value of priority gene regions designated for Lepidoptera molecular systematics. In particular, I assess how the DNA barcode segment of the mitochondrial COI gene performs across a broad temporal range given its number one position of priority, most sequenced status, and the conflicting opinions on its phylogenetic performance.

METHODOLOGY/PRINCIPAL FINDINGS

Gene regions commonly sequenced for lepidoptera phylogenetics were scored using multiple measures across three categories: practicality, which includes universality of primers and sequence quality; phylogenetic utility; and phylogenetic signal. I found that alternative measures within a category often appeared correlated, but high scores in one category did not necessarily translate into high scores in another. The DNA barcode was easier to sequence than other genes, and had high scores for utility but low signal above the genus level.

CONCLUSIONS/SIGNIFICANCE

Given limited financial resources and time constraints, careful selection of gene regions for molecular phylogenetics is crucial to avoid wasted effort producing partially informative data. This study introduces an approach to assessing the value of gene regions prior to the initiation of new studies and presents empirical results to help guide future selections.

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.

Repeated climate-linked host shifts have promoted diversification in a temperate clade of leaf-mining flies

A central but little-tested prediction of "escape and radiation" coevolution is that colonization of novel, chemically defended host plant clades accelerates insect herbivore diversification. That theory, in turn, exemplifies one side of a broader debate about the relative influence on clade dynamics of intrinsic (biotic) vs. extrinsic (physical-environmental) forces. Here, we use a fossil-calibrated molecular chronogram to compare the effects of a major biotic factor (repeated shift to a chemically divergent host plant clade) and a major abiotic factor (global climate change) on the macroevolutionary dynamics of a large Cenozoic radiation of phytophagous insects, the leaf-mining fly genus Phytomyza (Diptera: Agromyzidae). We find one of the first statistically supported examples of consistently elevated net diversification accompanying shift to new plant clades. In contrast, we detect no significant direct effect on diversification of major global climate events in the early and late Oligocene. The broader paleoclimatic context strongly suggests, however, that climate change has at times had a strong indirect influence through its effect on the biotic environment. Repeated rapid Miocene radiation of these flies on temperate herbaceous asterids closely corresponds to the dramatic, climate-driven expansion of seasonal, open habitats.

Genomic outposts serve the phylogenomic pioneers: designing novel nuclear markers for genomic DNA extractions of lepidoptera

Increasing the number of characters used in phylogenetic studies is the next crucial step towards generating robust and stable phylogenetic hypotheses - i.e., strongly supported and consistent across reconstruction method. Here we describe a genomic approach to finding new protein-coding genes for systematics in nonmodel taxa, which can be PCR amplified from standard, slightly degraded genomic DNA extracts. We test this approach on Lepidoptera, searching the draft genomic sequence of the silk moth Bombyx mori, for exons > 500 bp in length, removing annotated gene families, and compared remaining exons with butterfly EST databases to identify conserved regions for primer design. These primers were tested on a set of 65 taxa primarily in the butterfly family Nymphalidae. We were able to identify and amplify six previously unused gene regions (Arginine Kinase, GAPDH, IDH, MDH, RpS2, and RpS5) and two rarely used gene regions (CAD and DDC) that when added to the three traditional gene regions (COI, EF-1alpha and wingless) gave a data set of 8114 bp. Phylogenetic robustness and stability increased with increasing numbers of genes. Smaller taxanomic subsets were also robust when using the full gene data set. The full 11-gene data set was robust and stable across reconstruction methods, recovering the major lineages and strongly supporting relationships within them. Our methods and insights should be applicable to taxonomic groups having a single genomic reference species and several EST databases from taxa that diverged less than 100 million years ago.

Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes

Colletidae comprise approximately 2500 species of bees primarily distributed in the southern continents (only two colletid genera are widely distributed: Colletes and Hylaeus). Previously published studies have failed to resolve phylogenetic relationships on a worldwide basis and this has been a major barrier to the progress of research regarding systematics and evolution of colletid bees. For this study, data from four nuclear gene loci: elongation factor-1alpha (F2 copy), opsin, wingless, and 28S rRNA were analyzed for 122 species of colletid bees, representing all subfamilies and tribes currently recognized; 22 species belonging to three other bee families were used as outgroups. Bayesian, maximum likelihood, and parsimony methods were employed to investigate the phylogenetic relationships within Colletidae and resulted in highly congruent and well-resolved trees. The phylogenetic results show that Colletidae are monophyletic and that all traditionally recognized subfamilies (except Paracolletinae) are also strongly supported as monophyletic. Our phylogenetic hypothesis provides a framework within which broad questions related to the taxonomy, biogeography, morphology, evolution, and ecology of colletid bees can be addressed.

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.

Why do Manduca sexta feed from white flowers? Innate and learnt colour preferences in a hawkmoth

Flower colour is an important signal used by flowering plants to attract pollinators. Many anthophilous insects have an innate colour preference that is displayed during their first foraging bouts and which could help them locate their first nectar reward. Nevertheless, learning capabilities allow insects to switch their colour preferences with experience and thus, to track variation in floral nectar availability. Manduca sexta, a crepuscular hawkmoth widely studied as a model system for sensory physiology and behaviour, visits mostly white, night-blooming flowers lacking UV reflectance throughout its range in the Americas. Nevertheless, the spectral sensitivity of the feeding behaviour of naïve moths shows a narrow peak around 450 nm wavelengths, suggesting an innate preference for the colour blue. Under more natural conditions (i.e. broader wavelength reflectance) than in previous studies, we used dual choice experiments with blue- and white-coloured feeders to investigate the innate preference of naïve moths and trained different groups to each colour to evaluate their learning capabilities. We confirmed the innate preference of M. sexta for blue and found that these moths were able to switch colour preferences after training experience. These results unequivocally demonstrate that M. sexta moths innately prefer blue when presented against white flower models and offer novel experimental evidence supporting the hypothesis that learning capabilities could be involved in their foraging preferences, including their widely observed attraction to white flowers in nature.

Ancient rapid radiations of insects: challenges for phylogenetic analysis

Phylogenies of major groups of insects based on both morphological and molecular data have sometimes been contentious, often lacking the data to distinguish between alternative views of relationships. This paucity of data is often due to real biological and historical causes, such as shortness of time spans between divergences for evolution to occur and long time spans after divergences for subsequent evolutionary changes to obscure the earlier ones. Another reason for difficulty in resolving some of the relationships using molecular data is the limited spectrum of genes so far developed for phylogeny estimation. For this latter issue, there is cause for current optimism owing to rapid increases in our knowledge of comparative genomics. At least some historical patterns of divergence may, however, continue to defy our attempts to completely reconstruct them with confidence, at least using current strategies.

The phylogenetic relationships of flies in the superfamily Empidoidea (Insecta: Diptera)

We conducted a molecular phylogenetic study of the Empidoidea, a diverse group of 10,000 species of true flies, with two major goals: to reconstruct a taxonomically complete and robustly supported phylogeny for the group and to use this information to assess several competing classifications for the clade. We amassed 3900+ nucleotides of coding data from the carbamoylphosphate synthase domain of the rudimentary locus (CAD) and 1200+ nucleotides from the large nuclear ribosomal subunit (28S) from 72 and 71 species, respectively, representing several orthorrhaphan and cyclorrhaphan families and all previously recognized empidoidean subfamilies. Independent and combined phylogenetic analyses of these data were conducted using parsimony, maximum likelihood, and Bayesian criteria. The combined matrix included 61 taxa for which both CAD and 28S sequences were obtained. Analyses of CAD first and second codon positions alone and when concatenated with 28S sequences yielded trees with similar and largely stable topologies. Analyses of 28S data alone supported many clades although resolution is limited by low sequence divergence. The following major empidoid clades were recovered with convincing support in a majority of analyses: Atelestidae, Empidoidea exclusive of Atelestidae, Hybotidae sensu lato, Dolichopodidae+Microphorinae (including Parathallassius), and Empididae sensu lato (including Brachystomatinae, Ceratomerinae, Clinocerinae, Empidinae, Hemerodromiinae, Oreogetoninae, and Trichopezinae). The branching arrangement among these four major clades was Atelestidae, Hybotidae, Dolichopodidae/Microphorinae, Empididae. Previously recognized subclades recovered with robust support included Hybotinae, Brachystomatinae, Tachydromiinae, Clinocerinae (in part), Hemerodromiinae, Empidinae, and Empidiini.

Dissecting the ancient rapid radiation of microgastrine wasp genera using additional nuclear genes

Previous estimates of a generic level phylogeny for the ubiquitous parasitoid wasp subfamily Microgastrinae (Hymenoptera) have been problematic due to short internal branches deep in the phylogeny. These short branches might be attributed to a rapid radiation among the taxa, the use of genes that are unsuitable for the levels of divergence being examined, or insufficient quantity of data. We added over 1200 nucleotides from four nuclear genes to a dataset derived from three genes to produce a dataset of over 3000 nucleotides per taxon. While the number of well-supported short branches in the phylogeny increased, we still did not obtain strong bootstrap support for every node. Parametric and nonparametric bootstrap simulations projected that an enormous, and likely unobtainable, amount of data would be required to get bootstrap support greater than 50% for every node. However, a marked increase in the number of well-supported nodes was seen when we conducted a Bayesian analysis of a combined dataset generated from morphological characters added to the seven gene dataset. Our results suggest that, in some cases, combining morphological and genetic characters may be the most practical way to increase support for short branches deep in a phylogeny.

Gelechioidea (Insecta: Lepidoptera) systematics: A reexamination using combined morphology and mitochondrial DNA data

Phylogenies of Gelechioidea (Insecta: Lepidoptera) historically have been in disagreement and definitions vary at the family and subfamily levels. Addition of new taxa or new characters drastically changes relationships indicating that current phylogenetic schemes require more investigation. This study is the first phylogenetic analysis of Gelechioidea to include molecular data. Here we present a combined analysis using Maximum Parsimony to investigate sister-group relationships within Gelechioidea. The addition of Cytochrome oxidase I and II to revised published morphological matrices gives 453 parsimony informative characters for the 42 taxa for which we have sequence data. The combined analysis resulted in two trees with mostly novel sister-group relationships. These results challenge current concepts of Gelechioidea, suggesting that traditional morphological characters that have united taxa may not be homologous structures and are in need of further investigation. A combination of morphological data with new molecular data will be the most robust method of study for Gelechioidea phylogenetics.

When and where did troidine butterflies (Lepidoptera : Papilionidae) evolve? Phylogenetic and biogeographic evidence suggests an origin in remnant Gondwana in the Late Cretaceous

The age, geographic origin and time of major radiation of the butterflies (Hesperioidea + Papilionoidea + Hedyloidea) are largely unknown. The general modern view is that butterflies arose during the Late Jurassic/Cretaceous in the southern hemisphere (southern Pangea/Gondwana before continental breakup), but this is not universally accepted, and is a best guess based largely on circumstantial evidence. The extreme paucity of fossils and lack of modern, higher-level phylogenies of extant monophyletic groups have been major impediments towards determining reliable estimates of either their age or geographic origin. Here we present a phylogenetic and historical biogeographic analysis of a higher butterfly taxon, the swallowtail tribe Troidini. We analysed molecular data for three protein-encoding genes, mitochondrial ND5 and COI–COII, and nuclear EF–1α, both separately and in combination using maximum parsimony (with and without character weighting and transition/transversion weighting), maximum likelihood and Bayesian methods. Our sample included representatives of all 10 genera of Troidini and distant ingroup taxa (Baroniinae, Parnassiinae, Graphiini, Papilionini), with Pieridae as outgroup. Analysis of the combined dataset (4326 bp; 1012 parsimony informative characters) recovered the Troidini as a well supported monophyletic group and the monophyly of its two subtribes, Battina and Troidina. The most parsimonious biogeographic hypothesis suggests a southern origin of the tribe in remnant Gondwana (Madagascar–Greater India–Australia–Antarctica–South America) sometime after the rifting and final separation of Africa in the Late Cretaceous (<90 Mya). Although an ancient vicariance pattern is proposed, at least four relatively recent dispersal/extinction events are needed to reconcile anomalies in distribution, most of which can be explained by geological and climatic events in South-east Asia and Australia during the late Tertiary. Application of a molecular clock based on a rate smoothing programme to estimate various divergence times based on vicariance events, revealed two peculiarities in our biogeographic vicariance model that do not strictly accord with current understanding of the temporal breakup of Gondwana: (1) the troidine fauna of Greater India did not become isolated from Gondwana (Antarctica) until the end of the Cretaceous (c. 65 Mya), well after Madagascar separated from Greater India (84 Mya); and (2) the faunas of Greater India, Australia and South America diverged simultaneously, also at the K/T boundary. A recent published estimate of the time (31 Mya) of divergence between Cressida Swainson (Australia) and Euryades Felder & Felder (South America) is shown to be in error.

Single-copy nuclear genes recover cretaceous-age divergences in bees

We analyzed the higher level phylogeny of the bee family Halictidae based on the coding regions of three single-copy nuclear genes (long-wavelength [LW] opsin, wingless, and elongation factor 1-alpha [EF-1 alpha]). Our combined data set consisted of 2,234 aligned nucleotide sites (702 base pairs [bp] for LW opsin, 405 bp for wingless, and 1,127 bp for EF-1 alpha) and 779 parsimony-informative sites. We included 58 species of halictid bees from 33 genera, representing all subfamilies and tribes, and rooted the trees using seven outgroups from other bee families: Colletidae, Andrenidae, Melittidae, and Apidae. We analyzed the separate and combined data sets by a variety of methods, including equal weights parsimony, maximum likelihood, and Bayesian methods. Analysis of the combined data set produced a strong phylogenetic signal with high bootstrap and Bremer support and high posterior probability well into the base of the tree. The phylogeny recovered the monophyly of the Halictidae and of all four subfamilies and both tribes, recovered relationships among the subfamilies and tribes congruent with morphology, and provided robust support for the relationships among the numerous genera in the tribe Halictini, sensu Michener (2000). Using our combined nucleotide data set, several recently described halictid fossils from the Oligocene and Eocene, and recently developed Bayesian methods, we estimated the antiquity of major clades within the family. Our results indicate that each of the four subfamilies arose well before the Cretaceous-Tertiary boundary and suggest that the early radiation of halictid bees involved substantial African-South American interchange roughly coincident with the separation of these two continents in the late Cretaceous. This combination of single-copy nuclear genes is capable of recovering Cretaceous-age divergences in bees with high levels of support. We propose that LW opsin, wingless, and EF-1 alpha(F2 copy) may be useful in resolving relationships among bee families and other Cretaceous-age insect lineages.

A genomic BAC library and a new BAC-GFP vector to study the holocentric pest Spodoptera frugiperda

Two genomic tools for the study of Lepidoptera and the holocentric structure of their chromosomes are presented in this paper. A bacterial artificial chromosome (BAC) library was constructed using nuclear DNA partially digested with HindIII from eggs of Spodoptera frugiperda. The library contains a total of 36,864 clones with an average insert size of 125 kb, which corresponds to approximately 11.5 genome equivalents. Hybridization screening of the library was performed with eight single-copy genes, giving an average hit of 10 clones per marker gene. Colinearity between the genome and BACs was demonstrated at the triose phosphate isomerase (tpi) locus. Probing of the library with a PCR fragment internal to the 18S ribosomal gene allowed an estimation of the rDNA locus size close to 115 repeats per haploid genome. A new vector (pBAC3.6eGFP) for transient transfection into S. frugiperda cell lines has been constructed. It is based on the BAC vector, pBAC3.6e, in which a gene encoding GFP was inserted under the control of the densovirus P9 promoter. This vector has the advantage to accommodate large genomic inserts and to be transfected in a large lepidopteran host range. It was used to construct a second BAC library from Sf9 cell nuclear DNA in order to allow a comparison between somatic and cell line genome organization.

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

We sequenced nearly the entire carbomoylphosphate synthase (CPS) domain of CAD, or rudimentary, (ca. 4 kb) from 29 species of flies representing all major clades within Eremoneura, or higher flies, and several orthorrhaphous brachyceran outgroups. We compared these sequences with orthologs from Anopheles gambiae and Drosophila melanogaster to assess structure, compositional bias, and phylogenetic utility. CAD is large (6.6+ kb), complex (comprised of three major and myriad minor functional domains) and relatively free of introns, extreme nucleotide bias (except third codon positions), and large hypervariable regions. The CPS domain possesses moderate levels of nonsynonymous divergence among taxa of intermediate evolutionary age and conveys considerable phylogenetic signal. Phylogenetic analysis of CPS sequences under varying methods and assumptions resulted in well-resolved, strongly supported trees concordant with many traditional ideas about higher dipteran phylogeny and with prior inferences from 28S rDNA. The most robustly supported major eremoneuran clades were Cyclorrhapha, Platypezoidea, Eumuscomorpha, Empidoidea, Atelestidae, Empidoidea exclusive of Atelestidae, Hybotidae s.l., Microphoridae+Dolichopodidae, and Empididae s. str. Because CAD is ubiquitous, apparently single copy (at least within holometabolous insects), readily obtained from several insect orders using primers described herein, and exhibits considerable phylogenetic utility, it should have wide applicability in insect molecular systematics.

How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets

We analyzed 12 combined mitochondrial and nuclear gene datasets in seven orders of insects using both equal weights parsimony (to evaluate phylogenetic utility) and Bayesian methods (to investigate substitution patterns). For the Bayesian analyses we used relatively complex models (e.g., general time reversible models with rate variation) that allowed us to quantitatively compare relative rates among genes and codon positions, patterns of rate variation among genes, and substitution patterns within genes. Our analyses indicate that nuclear and mitochondrial genes differ in a number of important ways, some of which are correlated with phylogenetic utility. First and most obviously, nuclear genes generally evolve more slowly than mitochondrial genes (except in one case), making them better markers for deep divergences. Second, nuclear genes showed universally high values of CI and (generally) contribute more to overall tree resolution than mitochondrial genes (as measured by partitioned Bremer support). Third, nuclear genes show more homogeneous patterns of among-site rate variation (higher values of alpha than mitochondrial genes). Finally, nuclear genes show more symmetrical transformation rate matrices than mitochondrial genes. The combination of low values of alpha and highly asymmetrical transformation rate matrices may explain the overall poor performance of mitochondrial genes when compared to nuclear genes in the same analysis. Our analyses indicate that some parameters are highly correlated. For example, A/T bias was positively and significantly associated with relative rate and CI was positively and significantly associated with alpha (the shape of the gamma distribution). These results provide important insights into the substitution patterns that might characterized high quality genes for phylogenetic analysis: high values of alpha, unbiased base composition, and symmetrical transformation rate matrices. We argue that insect molecular systematists should increasingly focus on nuclear rather than mitochondrial gene datasets because nuclear genes do not suffer from the same substitutional biases that characterize mitochondrial genes.

Sequence and phylogenetic analysis of the complete mitochondrial genome of the flour beetle Tribolium castanaeum

We describe the first complete mitochondrial genome sequence from a representative of the insect order Coleoptera, the flour beetle Tribolium castaneum. The 15,881 bp long Tribolium mitochondrial genome encodes 13 putative proteins, two ribosomal RNAs and 22 tRNAs canonical for animal mitochondrial genomes. Their arrangement is identical to that in Drosophila melanogaster, which is considered ancestral for insects and crustaceans (Boore et al., 1998; Hwang, et al., 2001a). Nucleotide composition, amino acid composition, and codon usage fall within the range of values observed in other insect mitochondrial genomes. Most notable features are the use of TCT as tRNA(Ser(AGN)) anticodon instead of GCT, which is used in most other arthropod species, and the relative scarcity of special sequence motifs in the 1431 bp long control region. Phylogenetic analysis confirmed resolving power in the conserved regions of the mitochondrial proteome regarding diversification events, which predate the emergence of pterygote insects, while little resolution was obtained at the level of basal perygote diversification. The partition of faster evolving amino acid sites harbored strong support for joining Lepidoptera with Diptera, which is consistent with a monophyletic Mecopterida.

Reconstructing species phylogeny of the carabid beetles Ohomopterus using multiple nuclear DNA sequences: heterogeneous information content and the performance of simultaneous analyses

We attempted a phylogenetic reconstruction for the carabid subgenus Ohomopterus (genus Carabus), a notable case of radiation with mitochondrial introgression across species. Sequence data from five nuclear single copy loci were used, including wingless (Wg), phosphoenolpyruvate carboxykinase (PepCK), cytochrome c (Cytc), elongation factor-1alpha (EF-1alpha), and an anonymous single copy locus (Carab1). Sequences of Cytc, EF-1alpha, and Carab1 included intron or intron-like parts with length variation. The analysis of individual loci resulted in low resolution of the phylogenetic relationships, and the monophyly of several morphologically recognized species for which multiple specimens were analyzed was not revealed. Several specimens were heterozygous, with non-monophyletic alleles observed in three of the five loci at which alleles in heterozygotes were separated. In a simultaneous analysis of the five loci with ambiguously aligned parts eliminated and heterozygotic sites treated as missing, the resulting tree was well resolved, but the branch support was generally weak because of conflicting phylogenetic signals from different loci. We also attempted to incorporate allelic sequence data plus the ambiguously aligned parts in the analysis, by using all possible combinations of alleles from different loci in heterozygotic individuals, but the resultant tree was not supported more strongly. Nonetheless, these simultaneous analyses provided support for the monophyly of several species and species groups, and revealed the basic evolutionary trend of OHOMOPTERUS: initial widespread groups with simpler genitalia and the origination of exaggerated genitalia in a derived clade. This study exemplifies problems inherent in the phylogenetic reconstruction of closely related organisms where low levels of variation limit the information content from each locus, while heterozygosity, different phylogenetic history of multiple loci, and alignment ambiguity further hamper phylogenetic reconstruction unless several loci converge on a uniform signal.

Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

The organization, structure, and nucleotide variability of the ribosomal repeat unit was compared among families, genera, and species of cockroaches (Insecta:Blattodea). Sequence comparisons and molecular phylogenetic analyses were used to describe rDNA repeat unit variation at differing taxonomic levels. A reverse similar 1200 bp fragment of the 28S rDNA sequence was assessed for its potential utility in reconstructing higher-level phylogenetic relationships in cockroaches. Parsimony and maximum likelihood analyses of these data strongly support the expected pattern of relationships among cockroach groups. The examined 5' end of the 28S rDNA is shown to be an informative marker for larger studies of cockroach phylogeny. Comparative analysis of the nucleotide sequences of the rDNA internal transcribed spacers (ITS1 and ITS2) among closely related species of Blattella and Periplaneta reveals that ITS sequences can vary widely in primary sequence, length, and folding pattern. Secondary structure estimates for the ITS region of Blattella species indicate that variation in this spacer region can also influence the folding pattern of the 5.8S subunit. These results support the idea that ITS sequences play an important role in the stability and function of the rRNA cluster.

The long and short of sperm polymorphisms in insects

Production of more than one morphological type of sperm in a common testis has been documented for a variety of invertebrates, including gastropods, spiders, centipedes, and insects. This unusual phenomenon is difficult to explain by current theory, particularly since available evidence indicates that one sperm type is often incapable of effecting fertilization. In this review we critically examine evidence on the distribution and development of sperm heteromorphisms among insects in light of competing hypotheses for the evolutionary origin, maintenance, and function of a non-fertilizing class of sperm. To date, no single hypothesis, including alternatives which assume non-fertilizing sperm are non-adaptive, or that they provision, facilitate, or compete with fertilizing sperm, has received strong empirical support by any group of insects. The diversity of sperm heteromorphisms suggests that non-fertilizing sperm may have different functions in different clades or even serve multiple functions within a clade. We suggest that insight could be gained from (1) new models for the evolution of sperm polymorphism, (2) comparative studies that focus on multiple traits simultaneously (e.g. sperm number, proportion, length, and remating rate) and utilize clades in which more than one gain or loss of sperm heteromorphism has been documented (e.g. Pentatomidae, Carabidae, or Diopsidae), and (3) experimental studies that exploit individual variation or directly manipulate the composition of the male ejaculate.

A comparative study of odorant binding protein genes: differential expression of the PBP1-GOBP2 gene cluster inManduca sexta(Lepidoptera) and the organization of OBP genes inDrosophila melanogaster(Diptera)

Insects discriminate odors using sensory organs called olfactory sensilla, which display a wide range of phenotypes. Sensilla express ensembles of proteins, including odorant binding proteins (OBPs), olfactory receptors (ORs) and odor degrading enzymes (ODEs); odors are thought to be transported to ORs by OBPs and subsequently degraded by ODEs. These proteins belong to multigene families. The unique combinatorial expression of specific members of each of these gene families determines, in part, the phenotype of a sensillum and what odors it can detect. Furthermore, OBPs, ORs and ODEs are expressed in different cell types, suggesting the need for cell–cell communication to coordinate their expression. This report examines the OBP gene family. In Manduca sexta, the genes encoding PBP1Msex and GOBP2Msex are sequenced, shown to be adjacent to one another, and characterized together with OBP gene structures of other lepidoptera and Drosophila melanogaster. Expression of PBP1Msex, GOBP1Msex and GOBP2Msex is characterized in adult male and female antenna and in larval antenna and maxilla. The genomic organization of 25 D. melanogaster OBPs are characterized with respect to gene locus, gene cluster, amino acid sequence similarity, exon conservation and proximity to OR loci, and their sequences are compared with 14 M. sexta OBPs. Sensilla serve as portals of important behavioral information, and genes supporting sensilla function are presumably under significant evolutionary selective pressures. This study provides a basis for studying the evolution of the OBP gene family, the regulatory mechanisms governing the coordinated expression of OBPs, ORs and ODEs, and the processes that determine specific sensillum phenotypes.

18S rDNA phylogeny of the Tubificidae (Clitellata) and its constituent taxa: dismissal of the Naididae

The phylogeny of the Tubificidae, and of most of its subfamilies and some of its genera, is revisited, on the basis of sequences of 18S ribosomal DNA in a selection of species. Forty-six new 18S sequences of Naididae (6), Tubificidae (37), Phreodrilidae (1), Lumbriculidae (1), and Enchytraeidae (1) are reported and aligned together with corresponding sequences of 21 previously studied taxa. The 18S gene of Insulodrilus bifidus provides the first molecular evidence that phreodrilids are closely related to tubificids, corroborating previous conclusions based on morphology. The data further support the monophyletic status of Tubificidae, provided that the "Naididae" is regarded a part of this family; "naidids" may not even constitute a monophyletic group. It is thus suggested that the family name Naididae is formally suppressed as a junior synonym of the Tubificidae. The 18S gene also resolves a number of relationships within the tubificids. Among the subfamilies, Tubificinae is supported, Rhyacodrilinae and Phallodrilinae are revealed as nonmonophyletic, and Limnodriloidinae remains unresolved. Most tubificid genera tested for monophyly are corroborated by the data, only one (Tubifex) is refuted, and two (Tubificoides and Limnodriloides) are unresolved from other taxa. It is concluded that it will be valuable to expand the taxonomic sampling for 18S rDNA in clitellates, and in annelids in general, as this is likely to improve the resolution at many levels. However, it will be equally important to combine the annelid 18S data with other gene sequences and nonmolecular characters, to estimate the phylogeny of these common and diverse worms with greater precision.

A maximum-likelihood analysis of eight phylogenetic markers in gallwasps (Hymenoptera: Cynipidae): implications for insect phylogenetic studies

We assessed the utility of eight DNA sequence markers (5.8S rDNA, 18S rDNA, 28S rDNA, ITS regions, long-wavelength opsin, elongation factor 1-alpha, cytochrome b, and cytochrome oxidase I) in reconstructing phylogenetic relationships at various levels of divergence in gallwasps (Hymenoptera: Cynipidae), using a set of eight exemplar taxa. We report sequence divergence values and saturation levels and compare phylogenetic results of these sequences analyzed both separately and combined to a well-corroborated morphological phylogeny. Likelihood ratio tests were used to find the best evolutionary model fitting each of the markers. The likelihood model best explaining the data is, for most loci, parameter rich, with strong A-T bias for mitochondrial loci and strong rate heterogeneity for the majority of loci. Our data suggest that 28S rDNA, elongation factor 1-alpha, and long-wavelength opsin may be potentially useful markers for the resolution of cynipid and other insect within-family-level divergences (circa 50-100 mya old), whereas mitochondrial loci and ITS regions are most useful for lower-level phylogenetics. In contrast, the 18S rDNA marker is likely to be useful for the resolution of above-family-level relationships.

Phylogeny of feather mite subfamily Avenzoariinae (Acari: Analgoidea: Avenzoariidae) inferred from combined analyses of molecular and morphological data

Phylogenetic relationships among feather mites of the subfamily Avenzoariinae (Acari: Analgoidea: Avenzoariidae) were reconstructed by parsimony analysis of a combined data matrix. We analyzed 41 morphological characters and 246 molecular characters from a fragment of the 16S rDNA. Morphological trees were well supported at deep branches (genera and above), but showed much less support and resolution within genera. Molecular analyses produced trees with better resolution and support on terminal branches and worse support on basal branches. I(MF) index for the combined matrix pointed to the significant congruence of both data subsets with the whole of the data. The topology of the combined tree was close to the morphological tree in the deep branches and had well-resolved terminal branches as in the molecular tree. This suggests a considerable level of complimentarity between the two data sets. An analysis of association patterns of the mites and their hosts was conducted based on the results of the combined analyses for the Avenzoariinae and a phylogeny of their charadriiform hosts (compiled from various bird phylogeny hypotheses). The trees could be reconciled by the invoking of 12-13 cospeciation events, 6-7 duplications, 2 host shifts, and 26-29 sorting events. This suggests a high degree of cospeciation.

Elongation factor-2: a useful gene for arthropod phylogenetics

Robust resolution of controversial higher-level groupings within Arthropoda requires additional sources of characters. Toward this end, elongation factor-2 sequences (1899 nucleotides) were generated from 17 arthropod taxa (5 chelicerates, 6 crustaceans, 3 hexapods, 3 myriapods) plus an onychophoran and a tardigrade as outgroups. Likelihood and parsimony analyses of nucleotide and amino acid data sets consistently recovered Myriapoda and major chelicerate groups with high bootstrap support. Crustacea + Hexapoda (= Pancrustacea) was recovered with moderate support, whereas the conflicting group Myriapoda + Hexapoda (= Atelocerata) was never recovered and bootstrap values were always <5%. With additional nonarthropod sequences included, one indel supports monophyly of Tardigrada, Onychophora, and Arthropoda relative to molluscan, annelidan, and mammalian outgroups. New and previously published sequences from RNA polymerase II (1038 nucleotides) and elongation factor-1alpha (1092 nucleotides) were analyzed for the same taxa. A comparison of bootstrap values from the three genes analyzed separately revealed widely varying values for some clades, although there was never strong support for conflicting groups. In combined analyses, there was strong bootstrap support for the generally accepted clades Arachnida, Arthropoda, Euchelicerata, Hexapoda, and Pycnogonida, and for Chelicerata, Myriapoda, and Pancrustacea, whose monophyly is more controversial. Recovery of some additional groups was fairly robust to method of analysis but bootstrap values were not high; these included Pancrustacea + Chelicerata, Hexapoda + Cephalocarida + Remipedia, Cephalocarida + Remipedia, and Malaocostraca + Cirripedia. Atelocerata (= Myriapoda + Hexapoda) was never recovered. Elongation factor-2 is now the second protein-encoding, nuclear gene (in addition to RNA polymerase II) to support Pancrustacea over Atelocerata. Atelocerata is widely cited in morphology-based analyses, and the discrepancy between results derived from molecular and morphological data deserves greater attention.