Phylogenetic signal in the COI, 16S, and 28S genes for inferring relationships among genera of Microgastrinae (Hymenoptera; Braconidae): evidence of a high diversification rate in this group of parasitoids

Phylogenetic relationships among the Braconidae (Hymenoptera: Ichneumonoidea) inferred from partial 16S rDNA, 28S rDNA D2, 18S rDNA gene sequences and morphological characters

Phylogenetic relationships among the Braconidae were examined using homologous 16S rDNA, 28S rDNA D2 region, and 18S rDNA gene sequences and morphological data using both PAUP* 4.0 and MRBAYES 3.0B4 from 88 in-group taxa representing 35 subfamilies. The monophyletic nature of almost all subfamilies, of which multiple representatives are present in this study, is well-supported except for two subfamilies, Cenocoelinae and Neoneurinae that should probably be treated as tribal rank taxa in the subfamily Euphorinae. The topology of the trees generated in the present study supported the existence of three large generally accepted lineage or groupings of subfamilies: two main entirely endoparasitic lineages of this family, referred to as the "helconoid complex" and the "microgastroid complex," and the third "the cyclostome." The Aphidiinae was recovered as a member of the non-cyclostomes, probably a sister group of Euphorinae or Euphorinae-complex. The basal position of the microgastroid complex among the non-cyclostomes has been found in all our analyses. The cyclostomes were resolved as a monophyletic group in all analyses if two putatively misplaced groups (Mesostoa and Aspilodemon) were excluded from them. Certain well-supported relationships evident in this family from the previous analyses were recovered, such as a sister-group relationships of Alysiinae+Opiinae, of Braconinae+Doryctinae, and a close relationship between Macrocentrinae, Xiphozelinae, Homolobinae, and Charmontinae. The relationships of "Ichneutinae + ((Adeliinae + Cheloninae) + (Miracinae + (Cardiochilinae + Microgastrinae)))" was confirmed within the microgastroid complex. The position of Acampsohelconinae, Blacinae, and Trachypetinae is problematic.

A Secondary Structural Model of the 28S rRNA Expansion Segments D2 and D3 for Chalcidoid Wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.

Dracula ant phylogeny as inferred by nuclear 28S rDNA sequences and implications for ant systematics (Hymenoptera: Formicidae: Amblyoponinae)

Ants are one of the most ecologically and numerically dominant families of organisms in almost every terrestrial habitat throughout the world, though they include only about 1% of all described insect species. The development of eusociality is thought to have been a driving force in the striking diversification and dominance of this group, yet we know little about the evolution of the major lineages of ants and have been unable to clearly determine their primitive characteristics. Ants within the subfamily Amblyoponinae are specialized arthropod predators, possess many anatomically and behaviorally primitive characters and have been proposed as a possible basal lineage within the ants. We investigate the phylogenetic relationships among the members of the subfamily, using nuclear 28S rDNA sequence data. Outgroups for the analysis include members of the poneromorph and leptanillomorph (Apomyrma, Leptanilla) ant subfamilies, as well as three wasp families. Parsimony, maximum likelihood, and Bayesian analyses provide strong support for the monophyly of a clade containing the two genera Apomyrma+Mystrium (100% bpp; 97% ML bs; and 97% MP bs), and moderate support for the monophyly of the Amblyoponinae as long as Apomyrma (Apomyrminae) is included (87% bpp; 57% ML bs; and 76% MP bs). Analyses did not recover evidence of monophyly of the Amblyopone genus, while the monophyly of the other genera in the subfamily is supported. Based on these results we provide a morphological diagnosis of the Amblyoponinae that includes Apomyrma. Among the outgroup taxa, Typhlomyrmex grouped consistently with Ectatomma, supporting the recent placement of Typhlomyrmex in the Ectatomminae. The results of this present study place the included ant subfamilies into roughly two clades with the basal placement of Leptanilla unclear. One clade contains all the Amblyoponinae (including Apomyrma), Ponerinae, and Proceratiinae (Poneroid clade). The other clade contains members from subfamilies Cerapachyinae, Dolichoderinae, Ectatomminae, Formicinae, Myrmeciinae, and Myrmicinae (Formicoid clade).

A molecular phylogeny of eurytomid wasps inferred from DNA sequence data of 28S, 18S, 16S, and COI genes

Using partial DNA sequence data from nuclear 28S and 18S genes and mitochondrial 16S and COI genes, we reconstructed a phylogeny of the family Eurytomidae. Both maximum parsimony and Bayesian methods were employed. The analysis revealed a significant incongruence between the mitochondrial genes and the nuclear genes, and we chose the results from the nuclear genes as our preferred hypothesis. Our phylogeny suggested that the family Eurytomidae is not a monophyletic group; neither are the genera Eurytoma and Bruchophagus. The monophyly of genera Sycophila and Plutarchia was well supported, as was the close association of the genera Aiolomorphus, Tenuipetiolus, Bephratelloides, and Phylloxeroxenus. Our phylogeny also revealed an anticipated pattern, in which species groups from the genera Eurytoma and Bruchophagus are often more closely related to other small genera than to other species groups of the same genus. Subsequent taxonomic revisions include elevating the subfamily Rileyinae to a family status and the divisions of the genera Eurytoma and Bruchophagus.

Tempo of speciation in a butterfly genus from the Southeast Asian tropics, inferred from mitochondrial and nuclear DNA sequence data

Molecular systematics is frequently beset with phylogenetic results that are not fully resolved. Researchers either state that the absence of resolution is due to character conflict, explosive speciation, or some combination of the two, but seldom do they carefully examine their data to distinguish between these causes. In this study, we exhaustively analyze a set of nuclear and mitochondrial nucleotide data for the Asian tropical butterfly genus Arhopala so as to highlight the causes of polytomies in the phylogenetic trees, and, as a result, to infer important biological events in the history of this genus. We began by using non-parametric statistical methods to determine whether the ambiguously resolved regions in these trees represent hard or soft polytomies. In addition we determined how this correlated to number of inferred changes on branches, using parametric maximum likelihood estimations. Based on congruent patterns in both mitochondrial and nuclear DNA sequences, we concluded that at two stages in the history of Arhopala there have been accelerated instances of speciation. One event, at the base of the phylogeny, generated many of the groups and subgroups currently recognized in this genus, while a later event generated another major clade consisting of both Oriental and Papuan species groups. Based on comparisons of closely related taxa, the ratio of instantaneous rate of evolution between mitochondrial and nuclear DNA evolution is established at approximately 3:1. The earliest radiation is dated between 7 and 11 Ma by a molecular clock analysis, setting the events generating much of the diversity of Arhopala at well before the Pleistocene. Periodical flooding of the Sunda plateau during interglacial periods was, therefore, not responsible for generating the major divisions in the genus Arhopala. Instead, we hypothesize that large-scale climatic changes taking place in the Miocene have induced the early acceleration in speciation.

Phylogenetic evidence for an ancient rapid radiation of Caribbean sponge-dwelling snapping shrimps (Synalpheus)

A common challenge in reconstructing phylogenies involves a high frequency of short internal branches, which makes basal relationships difficult to resolve. Often it is not clear whether this pattern results from insufficient or inappropriate data, versus from a rapid evolutionary radiation. The snapping shrimp genus Synalpheus, which contains in excess of 100 species and is a prominent component of coral-reef faunas worldwide, provides an example. Its taxonomy has long been problematic due to the subtlety of diagnostic characters and apparently widespread variability within species. Here we use partial mt COI and 16S rRNA sequences and morphological characters to reconstruct relationships among 31 species in the morphologically well-defined gambarelloides species group, a putative clade of obligate sponge associates that is mostly endemic to the Caribbean and contains the only known eusocial marine animals. Analysis of the combined data produced a single tree with good support for many terminal clades and for relationships with outgroups, but poor support for branches near the base of the gambarelloides group. Most basal branches are extremely short and terminal branches are long, suggesting a relatively ancient, but rapid radiation of the gambarelloides group. This hypothesis is supported by significant departure from a null model of temporally random cladogenesis. Calibration of divergence times among gambarelloides-group species using data from three geminate pairs of Synalpheus species separated by the isthmus of Panamá suggests a major radiation between approximately 5 and 7 Mya, a few My before final closure of the Panamanian seaway during a period of spreading carbonate environments in the Caribbean; a second, smaller radiation occurred approximately 4 Mya. This molecular evidence for a rapid radiation among Caribbean marine organisms in the late Miocene/early Pliocene is strikingly similar to patterns documented from fossil data for several other Caribbean reef-associated invertebrate taxa. The similar patterns and timing of cladogenesis evidenced by molecular and fossil data for different Caribbean and East Pacific taxa suggests that the radiation involved a wide range of organisms, and strengthens the case that poor basal resolution in the gambarelloides group of Synalpheus reflects a real evolutionary phenomenon. The rapid radiation also helps explain the historical difficulty of diagnosing species in Synalpheus.

Biological identifications through DNA barcodes

Although much biological research depends upon species diagnoses, taxonomic expertise is collapsing. We are convinced that the sole prospect for a sustainable identification capability lies in the construction of systems that employ DNA sequences as taxon 'barcodes'. We establish that the mitochondrial gene cytochrome c oxidase I (COI) can serve as the core of a global bioidentification system for animals. First, we demonstrate that COI profiles, derived from the low-density sampling of higher taxonomic categories, ordinarily assign newly analysed taxa to the appropriate phylum or order. Second, we demonstrate that species-level assignments can be obtained by creating comprehensive COI profiles. A model COI profile, based upon the analysis of a single individual from each of 200 closely allied species of lepidopterans, was 100% successful in correctly identifying subsequent specimens. When fully developed, a COI identification system will provide a reliable, cost-effective and accessible solution to the current problem of species identification. Its assembly will also generate important new insights into the diversification of life and the rules of molecular evolution.

Estimating the age of the polydnavirus/braconid wasp symbiosis

Polydnaviruses are essential components mediating host-parasitoid relationships between some braconid wasps and their caterpillar hosts largely by suppressing or misdirecting the host immune systems. The polydnavirus-wasp relationship is an unusual apparent mutualism between viruses and eukaryotes and remarkably has evolved to the stage where the two entities no longer can be considered separate. Estimations of the age of the polydnavirus-bearing clade of braconid wasps based on separate calculations from the mitochondrial 16S rRNA and COI genes and the nuclear 28S rRNA gene, calibrated using fossil data, converge to indicate a date of origin of approximately 73.7 +/- 10 million years ago. This range provides an upper bound on the time during which these wasps and viruses have been functionally associated.

Effects of the polydnavirus of Cotesia congregata on the immune system and development of non-habitual hosts of the parasitoid

Polydnaviruses (PDV) are obligate mutualistic symbionts found in association with some groups of parasitic Hymenoptera. In these groups, they suppress the immune response of the parasitoid's host and are required for successful parasitoid reproduction. Several PDV effects have been described in different experimental systems, but no clear picture of PDV mode of immunosuppression has emerged. No study to date has directly tested if PDV modes of action are evolutionarily conserved or divergent among parasitoid taxa within the Ichneumonoidea. We hypothesize the divergence in PDV mode of immunosuppression can be detected by identifying points of divergence in the immune response of different host species to PDV from one parasitoid species. This study tests the effects of purified PDV from Cotesia congregata on the immune response of three larval lepidopteran species that naturally are hosts of parasitoid species that differ in taxonomic relatedness to C. congregata. Here we demonstrate that despite associations with distantly related parasitoids (Ichneumonidae and Braconidae), Manduca sexta and Heliothis virescens showed similar patterns of increased glucose dehydrogenase (GLD) activity, suppressed cellular encapsulation in vitro, and increased time to pupation. In contrast, Lymantria dispar showed no response to C. congregata PDV across any of the parameters measured, even though it has an evolutionary association with several parasitoids closely related to C. congregata and within the Microgastrinae. The PDV immunosuppression in H. virescens and M. sexta does not correlate with host molecular phylogeny either. The suborganismal effects shown in M. sexta and H. virescens translated into significantly reduced pupation success in M. sexta only. Results demonstrate that while some PDV modes of immunosuppression in hosts may be divergent, others may be conserved across broad host groups.

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.

Limitations of relative apparent synapomorphy analysis (RASA) for measuring phylogenetic signal

In this paper we use hypothetical and empirical data matrices to evaluate the ability of relative apparent synapomorphy analysis (RASA) to measure phylogenetic signal, select outgroups, and identify terminals subject to long-branch attraction. In all cases, except for equal character-state frequencies, RASA indicated extraordinarily high levels of phylogenetic information for hypothetical data matrices that are uninformative regarding relationships among the terminals. Yet, regardless of the number of characters or character-state frequencies, RASA failed to detect phylogenetic signal for hypothetical matrices with strong phylogenetic signal. In our empirical example, RASA indicated increasing phylogenetic signal for matrices for which the strict consensus of the most parsimonious trees is increasingly poorly resolved, clades are increasingly poorly supported, and for which many relationships are in conflict with more widely sampled analyses. RASA is an ineffective approach to identify outgroup terminal(s) with the most plesiomorphic character states for the ingroup. Our hypothetical example demonstrated that RASA preferred outgroup terminals with increasing numbers of convergent character states with ingroup terminals, and rejected the outgroup terminal with all plesiomorphic character states. Our empirical example demonstrated that RASA, in all three cases examined, selected an ingroup terminal, rather than an outgroup terminal, as the best outgroup. In no case was one of the two outgroup terminals even close to being considered the optimal outgroup by RASA. RASA is an ineffective means of identifying problematic long-branch terminals. In our hypothetical example, RASA indicated a terminal as being a problematic long-branch terminal in spite of the terminal being on a zero-length branch and having no possibility of undergoing long-branch attraction with another terminal. RASA also failed to identify actual problematic long-branch terminals that did undergo long-branch attraction, but only after following Lyons-Weiler and Hoelzer's (1997) three-step process to identify and remove terminals subject to long-branch attraction. We conclude that RASA should not be used for any of these purposes.

Species-specific identification of human hookworms by PCR of the mitochondrial cytochrome oxidase I gene

Significant differences in the life histories of the human hookworms Ancylostoma duodenale and Necator americanus necessitate their differentiation for epidemiological studies and the design of control programs. Current methods of identification require time-consuming, labor-intensive techniques. A polymerase chain reaction (PCR)-based method that enables rapid species identification is described. The mitochondrial cytochrome oxidase I genes of both species were sequenced, and species-specific primer sets were designed. The primers were used in PCR to amplify 585-bp fragments of the cytochrome oxidase gene from individual hookworm eggs, larvae, and adults. The technique was also able to identify mixed infections containing equal amounts of eggs from each species. The technique is rapid, technically simple, and sensitive and will permit the accurate identification of human hookworms in epidemiological field studies.

Large amounts of genetic divergence among Italian species of the genus Orchesella (Insecta, collembola) and the relationships of two new species

The phylogenetic position of two putative new species of the collembolan genus Orchesella was investigated by comparison with four other Italian species of the genus using a fragment of the mitochondrial gene encoding for subunit I of cytochrome c oxidase (COI). The gene showed the well-known A + T bias, typical of insect mitochondrial DNA, although A + T content was not as high as that observed in species belonging to more derived insect orders. The large number of variable sites in 3rd codon positions (85.2% variable) suggested that these sites contain significant homoplasy due to multiple hits. Despite the lack of morphological differentiation, the COI portion examined shows remarkable levels of genetic divergence between the putative species and their closest relatives. Phylogenetic analysis suggests that one of the putative new species is related to O. villosa, whereas the other is included in a clade with O. cincta and O. ranzii. The species O. chiantica appears to be related to O. villosa, agreeing with previous allozyme data.

Resolving phylogeny at the family level by mitochondrial cytochrome oxidase sequences: phylogeny of carrion beetles (Coleoptera, Silphidae)

We investigated the phylogenetic relationships of carrion beetles (Coleoptera, Silphidae) using 2094 bp of their mitochondrial cytochrome oxidase subunit I and II and tRNA leucine gene sequences. Shorter fragments of this gene region previously have been used to establish generic relationships in insects. In this study, they provided more than sufficient resolution, although the third positions of the protein-coding sequences reached saturation for the deeper divergences. This first published phylogeny for the Silphidae comprises 23 species from 13 genera sampled across the geographic range of the family. In addition, we included species from three related families as outgroups. One of these families, the Agyrtidae, was, until recently, included in the Silphidae, but its resolution here justifies its current position as a separate family. The silphid subfamilies Nicrophorinae and Silphinae are monophyletic in all analyses. All genera for which several species were sampled are supported as monophyletic groups, with the exception of the genus Silpha. European and North American representatives of two Nicrophorus species described from both continents are supported as each others' closest relatives. The lineage that colonized Gondwanaland and that most likely originated in the Palearctic is the most basal within the Silphinae.