Priors and Posteriors in Bayesian Timing of Divergence Analyses: The Age of Butterflies Revisited.

The need for robust estimates of times of divergence is essential for downstream analyses, yet assessing this robustness is still rare. We generated a time-calibrated genus-level phylogeny of butterflies (Papilionoidea), including 994 taxa, up to 10 gene fragments and an unprecedented set of 12 fossils and 10 host-plant node calibration points. We compared marginal priors and posterior distributions to assess the relative importance of the former on the latter. This approach revealed a strong influence of the set of priors on the root age but for most calibrated nodes posterior distributions shifted from the marginal prior, indicating significant information in the molecular data set. Using a very conservative approach we estimated an origin of butterflies at 107.6 Ma, approximately equivalent to the latest Early Cretaceous, with a credibility interval ranging from 89.5 Ma (mid Late Cretaceous) to 129.5 Ma (mid Early Cretaceous). In addition, we tested the effects of changing fossil calibration priors, tree prior, different sets of calibrations and different sampling fractions but our estimate remained robust to these alternative assumptions. With 994 genera, this tree provides a comprehensive source of secondary calibrations for studies on butterflies.

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.

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.

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.

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.

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.

Preserving and vouchering butterflies and moths for large-scale museum-based molecular research.

Butterflies and moths (Lepidoptera) comprise significant portions of the world's natural history collections, but a standardized tissue preservation protocol for molecular research is largely lacking. Lepidoptera have traditionally been spread on mounting boards to display wing patterns and colors, which are often important for species identification. Many molecular phylogenetic studies have used legs from pinned specimens as the primary source for DNA in order to preserve a morphological voucher, but the amount of available tissue is often limited. Preserving an entire specimen in a cryogenic freezer is ideal for DNA preservation, but without an easily accessible voucher it can make specimen identification, verification, and morphological work difficult. Here we present a procedure that creates accessible and easily visualized "wing vouchers" of individual Lepidoptera specimens, and preserves the remainder of the insect in a cryogenic freezer for molecular research. Wings are preserved in protective holders so that both dorsal and ventral patterns and colors can be easily viewed without further damage. Our wing vouchering system has been implemented at the University of Maryland (AToL Lep Collection) and the University of Florida (Florida Museum of Natural History, McGuire Center of Lepidoptera and Biodiversity), which are among two of the largest Lepidoptera molecular collections in the world.

Nitrogen in insects: implications for trophic complexity and species diversification.

Disparities in nutrient content (nitrogen and phosphorus) between herbivores and their plant resources have lately proven to have major consequences for herbivore success, consumer-driven nutrient cycling, and the fate of primary production in ecosystems. Here we extend these findings by examining patterns of nutrient content between animals at higher trophic levels, specifically between insect herbivores and predators. Using a recently compiled database on insect nutrient content, we found that predators exhibit on average 15% greater nitrogen content than herbivores. This difference persists after accounting for variation from phylogeny and allometry. Among herbivorous insects, we also found evidence that recently derived lineages (e.g., herbivorous Diptera and Lepidoptera) have, on a relative basis, 15%-25% less body nitrogen than more ancient herbivore lineages (e.g., herbivorous Orthoptera and Hemiptera). We elaborate several testable hypotheses for the origin of differences in nitrogen content between trophic levels and among phylogenetic lineages. For example, interspecific variation in insect nitrogen content may be directly traceable to differences in dietary nitrogen (including dilution by gut contents), selected for directly in response to the differential scarcity of dietary nitrogen, or an indirect consequence of adaptation to different feeding habits. From some functional perspectives, the magnitude rather than the source of the interspecific differences in nitrogen content may be most critical. We conclude by discussing the implications of the observed patterns for both the trophic complexity of food webs and the evolutionary radiation of herbivorous insects.

EVOLUTIONARY ORIGIN OF THE CYCLORRHAPHA (DIPTERA): TEST OF ALTERNATIVE MORPHOLOGICAL HYPOTHESES.

Abstract- The higher flies, infraorder Cyclorrhapha [=Muscomorpha (McAlpine, 1989)], have undergone enormous radiation since the Cretaceous (∼100 Myr). Rapid morphological evolution in cyclorrhaphans has made their phylogenetic placement with respect to more primitive clades a long-standing problem in dipteran systematics. Of the two most plausible hypotheses, one treats the Cyclorrhapha as sister group to the orthorrhaphous superfamily Empidoidea [=Empidiformia (Hennig, 1948), Orthogenya (Brauer, 1883)], while the other places them within the empidoids. The debate over cyclorrhaphan origin has heretofore focused on homology interpretations for a few character systems, particularly the male genitalia. We provide the first attempt to assemble and quantify all of the available morphological evidence. By cladistic analysis of these data under alternative codings of genitalic features reflecting opposing homology theories, and then excluding these features altogether, we sought to judge which genitalic theory is better supported by the evidence as a whole, and how much the debate matters to resolving cyclorrhaphan origins. Using the analog of a factorial design, we also measured the effect of alternative transformation series in several other controversial characters, of outgroup choice and of successive weighting. Under all manipulations, including both genitalic codings, the Cyclorrhapha originate within the Empidoidea, near the family Atelestidae. However, trees in which the Empidoidea are constrained to be monophyletic are only 1-6 steps longer (out of ∼150), a fit not significantly worse under a permutation test for monophyly. Adult morphological data may not suffice to settle either the placement of Cyclorrhapha or the debate over genitalic homology. Moreover, the issue of genitalic homology does not appear critical to that of cyclorrhaphan origin.

Can RNA-Seq resolve the rapid radiation of advanced moths and butterflies (Hexapoda: Lepidoptera: Apoditrysia)? An exploratory study.

Recent molecular phylogenetic studies of the insect order Lepidoptera have robustly resolved family-level divergences within most superfamilies, and most divergences among the relatively species-poor early-arising superfamilies. In sharp contrast, relationships among the superfamilies of more advanced moths and butterflies that comprise the mega-diverse clade Apoditrysia (ca. 145,000 spp.) remain mostly poorly supported. This uncertainty, in turn, limits our ability to discern the origins, ages and evolutionary consequences of traits hypothesized to promote the spectacular diversification of Apoditrysia. Low support along the apoditrysian "backbone" probably reflects rapid diversification. If so, it may be feasible to strengthen resolution by radically increasing the gene sample, but case studies have been few. We explored the potential of next-generation sequencing to conclusively resolve apoditrysian relationships. We used transcriptome RNA-Seq to generate 1579 putatively orthologous gene sequences across a broad sample of 40 apoditrysians plus four outgroups, to which we added two taxa from previously published data. Phylogenetic analysis of a 46-taxon, 741-gene matrix, resulting from a strict filter that eliminated ortholog groups containing any apparent paralogs, yielded dramatic overall increase in bootstrap support for deeper nodes within Apoditrysia as compared to results from previous and concurrent 19-gene analyses. High support was restricted mainly to the huge subclade Obtectomera broadly defined, in which 11 of 12 nodes subtending multiple superfamilies had bootstrap support of 100%. The strongly supported nodes showed little conflict with groupings from previous studies, and were little affected by changes in taxon sampling, suggesting that they reflect true signal rather than artifacts of massive gene sampling. In contrast, strong support was seen at only 2 of 11 deeper nodes among the "lower", non-obtectomeran apoditrysians. These represent a much harder phylogenetic problem, for which one path to resolution might include further increase in gene sampling, together with improved orthology assignments.

A molecular phylogeny for yponomeutoidea (insecta, Lepidoptera, ditrysia) and its implications for classification, biogeography and the evolution of host plant use.

BACKGROUND: Yponomeutoidea, one of the early-diverging lineages of ditrysian Lepidoptera, comprise about 1,800 species worldwide, including notable pests and insect-plant interaction models. Yponomeutoids were one of the earliest lepidopteran clades to evolve external feeding and to extensively colonize herbaceous angiosperms. Despite the group's economic importance, and its value for tracing early lepidopteran evolution, the biodiversity and phylogeny of Yponomeutoidea have been relatively little studied. METHODOLOGY/PRINCIPAL FINDINGS: Eight nuclear genes (8 kb) were initially sequenced for 86 putative yponomeutoid species, spanning all previously recognized suprageneric groups, and 53 outgroups representing 22 families and 12 superfamilies. Eleven to 19 additional genes, yielding a total of 14.8 to 18.9 kb, were then sampled for a subset of taxa, including 28 yponomeutoids and 43 outgroups. Maximum likelihood analyses were conducted on data sets differing in numbers of genes, matrix completeness, inclusion/weighting of synonymous substitutions, and inclusion/exclusion of "rogue" taxa. Monophyly for Yponomeutoidea was supported very strongly when the 18 "rogue" taxa were excluded, and moderately otherwise. Results from different analyses are highly congruent and relationships within Yponomeutoidea are well supported overall. There is strong support overall for monophyly of families previously recognized on morphological grounds, including Yponomeutidae, Ypsolophidae, Plutellidae, Glyphipterigidae, Argyresthiidae, Attevidae, Praydidae, Heliodinidae, and Bedelliidae. We also assign family rank to Scythropiinae (Scythropiidae stat. rev.), which in our trees are strongly grouped with Bedelliidae, in contrast to all previous proposals. We present a working hypothesis of among-family relationships, and an informal higher classification. Host plant family associations of yponomeutoid subfamilies and families are non-random, but show no trends suggesting parallel phylogenesis. Our analyses suggest that previous characterizations of yponomeutoids as predominantly Holarctic were based on insufficient sampling. CONCLUSIONS/SIGNIFICANCE: We provide the first robust molecular phylogeny for Yponomeutoidea, together with a revised classification and new insights into their life history evolution and biogeography.

A large-scale, higher-level, molecular phylogenetic study of the insect order Lepidoptera (moths and butterflies).

BACKGROUND: Higher-level relationships within the Lepidoptera, and particularly within the species-rich subclade Ditrysia, are generally not well understood, although recent studies have yielded progress. We present the most comprehensive molecular analysis of lepidopteran phylogeny to date, focusing on relationships among superfamilies. METHODOLOGY PRINCIPAL FINDINGS: 483 taxa spanning 115 of 124 families were sampled for 19 protein-coding nuclear genes, from which maximum likelihood tree estimates and bootstrap percentages were obtained using GARLI. Assessment of heuristic search effectiveness showed that better trees and higher bootstrap percentages probably remain to be discovered even after 1000 or more search replicates, but further search proved impractical even with grid computing. Other analyses explored the effects of sampling nonsynonymous change only versus partitioned and unpartitioned total nucleotide change; deletion of rogue taxa; and compositional heterogeneity. Relationships among the non-ditrysian lineages previously inferred from morphology were largely confirmed, plus some new ones, with strong support. Robust support was also found for divergences among non-apoditrysian lineages of Ditrysia, but only rarely so within Apoditrysia. Paraphyly for Tineoidea is strongly supported by analysis of nonsynonymous-only signal; conflicting, strong support for tineoid monophyly when synonymous signal was added back is shown to result from compositional heterogeneity. CONCLUSIONS SIGNIFICANCE: Support for among-superfamily relationships outside the Apoditrysia is now generally strong. Comparable support is mostly lacking within Apoditrysia, but dramatically increased bootstrap percentages for some nodes after rogue taxon removal, and concordance with other evidence, strongly suggest that our picture of apoditrysian phylogeny is approximately correct. This study highlights the challenge of finding optimal topologies when analyzing hundreds of taxa. It also shows that some nodes get strong support only when analysis is restricted to nonsynonymous change, while total change is necessary for strong support of others. Thus, multiple types of analyses will be necessary to fully resolve lepidopteran phylogeny.

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

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.

Phylogenetics of eggshell morphogenesis in Antheraea (lepidoptera: saturniidae): unique origin and repeated reduction of the aeropyle crown.

Integrated phylogenetic and developmental analyses should enhance our understanding of morphological evolution and thereby improve systematists' ability to utilize morphological characters, but case studies are few. The eggshell (chorion) of Lepidoptera (Insecta) has proven especially tractable experimentally for such analyses because its morphogenesis proceeds by extracellular assembly of proteins. This study focuses on a morphological novelty, the aeropyle crown, that arises at the end of choriogenesis in the wild silkmoth genus Antheraea. Aeropyle crowns are cylindrical projections, ending in prominent prongs, that surround the openings of breathing tubes (aeropyle channels) traversing the chorion. They occur over the entire egg surface in some species, are localized to a circumferential band in many others, and in some are missing entirely, thus exhibiting variation typical of discrete characters analyzed in morphological phylogenetics. Seeking an integrated developmental-phylogenetic view, we first survey aeropyle crown variation broadly across Antheraea and related genera. We then map these observations onto a robust phylogeny, based on three nuclear genes, to test the adequacy of character codings for aeropyle crown variation and to estimate the frequency and direction of change in those characters. Thirdly, we draw on previous studies of choriogenesis, supplemented by new data on gene expression, to hypothesize developmental-genetic bases for the inferred chorion character transformations. Aeropyle crowns are inferred to arise just once, in the ancestor of Antheraea, but to undergo four or more subsequent reductions without regain, a pattern consistent with Dollo's Law. Spatial distribution shows an analogous trend, though less clear-cut, toward reduction of coverage by aeropyle crowns. These trends suggest either that there is little or no natural selection on the details of the aeropyle crown structure or that evolution toward functional optima is ongoing, although no direct evidence exists for either. Genetic, biochemical, and microscopy studies point to at least two developmental changes underlying the origin of the aeropyle crown, namely, reinitiation of deposition of chorionic lamellae after the end of normal choriogenesis (i.e., heterochrony), and sharply increased production of underlying "filler" proteins that push the nascent final lamellae upward to form the crown (i.e., heteroposy). Identification of a unique putative cis-regulatory element shared by unrelated genes involved in aeropyle crown formation suggests a possible simple mechanism for repeated evolutionary reduction and spatial restriction of aeropyle crowns.

PHYLOGENESIS OF INSECT/PLANT INTERACTIONS: HAVE PHYLLOBROTICA LEAF BEETLES (CHRYSOMELIDAE) AND THE LAMIALES DIVERSIFIED IN PARALLEL?

The relative importance of conservative versus locally adapted traits for species interactions is an increasingly common theme in evolutionary ecology. Obligate interactions such as those between parasites and hosts often exhibit such strong phylogenetic conservatism that current associations may reflect diversification in parallel. Parallel phylogenesis, documented for animal parasites, has been doubted for insect/plant interactions, but phylogenetic studies of highly specific insect/plant associations are very few. A comparison of phylogeny estimates for the strictly monophagous Phyllobrotica leaf beetles and their lamialean hostplants shows nearly complete concordance, strongly supporting the hypothesis of parallel diversification. The cladogram concordance is significant or nearly so (consensus index values equalling or exceeding the critical value) under randomization distributions based on Adams (though not Nelson) consensus trees. The one clear exception shows unusual natural history, suggesting an isolated host transfer. Insect distributions and plant fossil ages are consistent with a mid-Tertiary age for both clades, further disfavoring the alternative hypothesis of entirely subsequent evolution. The dependence of both larval and adult beetles on the hostplants, larval endophagy, and possible dependence of beetles on toxic host compounds for defense against predators are suggested to underlie the evolutionary persistence of this interaction. Current host use in these beetles appears to reflect primarily the phylogeny of the interaction, strengthening the thesis that history can play a major role in structuring insect/plant relationships.