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.

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

Compared to the number of genes available for study of both younger and older divergences, few genes have yet been identified that can strongly resolve phylogenetic splits of Mesozoic age ( approximately 65-250 mya). Thus, reconstruction of Mesozoic-age phylogenies, exemplified by basal divergences within the major orders of holometabolous insects, is likely to be especially dependent on combining multiple lines of evidence. This study tests the potential of the 18S ribosomal RNA gene for reconstructing Mesozoic-aged divergences within the insect order Lepidoptera and its ability when combined with a second, previously analyzed nuclear gene (phosphoenolpyruvate carboxykinase, PEPCK) to strongly resolve these relationships. 18S sequences were obtained for 21 taxa, representing major clades of Lepidoptera plus outgroups from the other "panorpoid orders. A well-corroborated morphology-based "test phylogeny was used to evaluate the effects of partitioning the 18S gene according to variable versus conserved domains, paired versus unpaired sites in the secondary structure, and transition versus transversion substitutions. Likelihood and unweighted parsimony analyses of the 18S data recover the "test phylogeny" almost completely, with no improvement of agreement or support provided by any form of weighting or partitioning. No conflict in signal between 18S and PEPCK was detected by the partition homogeneity test. Combined parsimony analysis yielded strong bootstrap support for nearly all relationships, much higher than for either gene alone, thereby also providing strong evidence on several hypotheses about the early evolution of lepidopteran-plant interactions. These genes in combination may be widely useful for resolving insect divergences of comparable age.

Phylogenetic utility of the nuclear gene dopa decarboxylase in noctuoid moths (Insecta: Lepidoptera: noctuoidea).

Phylogenetic utility for the nuclear gene encoding dopa decarboxylase (DDC), little used in systematics, was recently demonstrated within the noctuid moth subfamily Heliothinae. Here we extend the test of the utility of a 709-bp DDC fragment to deeper levels, analyzing 49 species representing major groups across the superfamily Noctuoidea. Parsimony, distance, and maximum-likelihood analyses recover all or nearly all of a set of "test clades" supported by clear morphological synapomorphies, spanning a wide range of taxonomic levels. DDC also upholds a recent proposal that the Noctuidae are paraphyletic. Nt3 contributes a majority of the signal and recovers the basal split between Notodontidae and all other noctuoids, despite a plateau of nt3 divergence at this level. However, nonsynonymous changes also support groups at all levels, and in contrast to nt3, amino acid divergence shows no plateau. The utility of DDC promises to extend back to the early Tertiary and Cretaceous, a time span for which few suitable genes have been identified.

Phylogenetic analysis of arthropods using two nuclear protein-encoding genes supports a crustacean + hexapod clade.

Recent phylogenetic analyses using molecular data suggest that hexapods are more closely related to crustaceans than to myriapods, a result that conflicts with long-held morphology-based hypotheses. Here we contribute additional information to this debate by conducting phylogenetic analyses on two nuclear protein-encoding genes, elongation factor-1 alpha (EF-1 alpha) and the largest subunit of RNA polymerase II (Pol II), from an extensive sample of arthropod taxa. Results were obtained from two data sets. One data set comprised 1092 nucleotides (364 amino acids) of EF-1 alpha and 372 nucleotides (124 amino acids) of Pol II from 30 arthropods and three lobopods. The other data set contained the same EF-1 alpha fragment and an expanded 1038-nucleotide (346-amino-acid) sample of Pol II from 17 arthropod taxa. Results from maximum-parsimony and maximum-likelihood analyses strongly supported the existence of a Crustacea + Hexapoda clade (Pancrustacea) over a Myriapoda + Hexapoda clade (Atelocerata). The apparent incompatibility between the molecule-based Pancrustacea hypothesis and morphology-based Atelocerata hypothesis is discussed.

More taxa or more characters revisited: combining data from nuclear protein-encoding genes for phylogenetic analyses of Noctuoidea (Insecta: Lepidoptera).

A central question concerning data collection strategy for molecular phylogenies has been, is it better to increase the number of characters or the number of taxa sampled to improve the robustness of a phylogeny estimate? A recent simulation study concluded that increasing the number of taxa sampled is preferable to increasing the number of nucleotide characters, if taxa are chosen specifically to break up long branches. We explore this hypothesis by using empirical data from noctuoid moths, one of the largest superfamilies of insects. Separate studies of two nuclear genes, elongation factor-1 alpha (EF-1 alpha) and dopa decarboxylase (DDC), have yielded similar gene trees and high concordance with morphological groupings for 49 exemplar species. However, support levels were quite low for nodes deeper than the subfamily level. We tested the effects on phylogenetic signal of (1) increasing the taxon sampling by nearly 60%, to 77 species, and (2) combining data from the two genes in a single analysis. Surprisingly, the increased taxon sampling, although designed to break up long branches, generated greater disagreement between the two gene data sets and decreased support levels for deeper nodes. We appear to have inadvertently introduced new long branches, and breaking these up may require a yet larger taxon sample. Sampling additional characters (combining data) greatly increased the phylogenetic signal. To contrast the potential effect of combining data from independent genes with collection of the same total number of characters from a single gene, we simulated the latter by bootstrap augmentation of the single-gene data sets. Support levels for combined data were at least as high as those for the bootstrap-augmented data set for DDC and were much higher than those for the augmented EF-1 alpha data set. This supports the view that in obtaining additional sequence data to solve a refractory systematic problem, it is prudent to take them from an independent gene.

Evolution and phylogenetic utility of the period gene in Lepidoptera.

Evolution and phylogenetic utility of the period gene are explored through sequence analysis of a relatively conserved 909-bp fragment in 26 lepidopteran species. Taxa range from tribes to superfamilies, primarily within the putative clade Macrolepidotera plus near outgroups, and include both strongly established and problematic groupings. Their divergence dates probably range from the late Cretaceous through much of the Tertiary. Comparisons within the same set of closely related species show that amino acid substitutions in period occur 4.9 and 44 times as frequently as they do in two other nuclear genes--dopa decarboxylase and elongation factor-1 alpha, respectively. In contrast, rates of observed synonymous substitution are within 60% of each other for these three genes. Synonymous changes in period approach saturation by the family level, whereas nonsynonymous and amino acid divergences across the Macrolepidoptera are less than half the maximal values reported for this gene. Phylogenetic analyses of period strongly supported groupings at the family level and below. In contrast to previous analyses at this level with other nuclear genes, much of the information lies in nonsynonymous change. Relationships up to the superfamily level were recovered with decreasing effectiveness, and little, if any, signal was apparent regarding relationships among superfamilies. This could reflect rapid radiation of the superfamilies, however, rather than saturation in the period locus; thus, period, in combination with other genes, remains a plausible candidate for approaching the difficult problems of lepidopteran family and superfamily relationships.

Two nuclear genes yield concordant relationships within Attacini (Lepidoptera: Saturniidae).

To extend initial characterizations of their phylogenetic utility, sequences from the nuclear genes for elongation factor-1 alpha (EF-1 alpha) and dopa decarboxylase (DDC) are tested for phylogenetic concordance with each other and with previous morphological evidence within the giant silk moth tribe Attacini (Lepidoptera: Saturniidae). The sampling of DDC is expanded from the 690 basepairs of previous studies to 1051 basepairs in the current study. All nine attacine genera are sampled. EF-1 alpha and DDC agree in the placement of seven of nine genera, with placement of the other two not in strong conflict. Combination of the gene sequences results in a nearly fully resolved tree that is consistent with EF-1 alpha alone and agrees with morphology in five of eight groups. Conflict between molecules and morphology is confined to deeper-level relationships within Attacini, where node support for the molecular hypotheses, but not the morphological hypotheses, is generally very strong. A strong signal is contributed by synonymous substitutions in both genes, and by nonsynonymous change particularly in DDC. The molecular phylogeny supports a revision of attacine biogeography in that neither East Asian nor New World genera form monophyletic groups.

Molecular phylogeny of the major arthropod groups indicates polyphyly of crustaceans and a new hypothesis for the origin of hexapods.

A phylogeny of the arthropods was inferred from analyses of amino acid sequences derived from the nuclear genes encoding elongation factor-1 alpha and the largest subunit of RNA polymerase II using maximum-parsimony, neighbor-joining, and maximum-likelihood methods. Analyses of elongation factor-1 alpha from 17 arthropods and 4 outgroup taxa recovered many arthropod clades supported by previous morphological studies, including Diplopoda, Myriapoda, Insecta, Hexapoda, Branchiopoda (Crustacea), Araneae, Tetrapulmonata, Arachnida, Chelicerata, and Malacostraca (Crustacea). However, counter to previous studies, elongation factor-1 alpha placed Malacostraca as sister group to the other arthropods. Branchiopod crustaceans were found to be more closely related to hexapods and myriapods than to malacostracan crustaceans. Sequences for RNA polymerase II were obtained from 11 arthropod taxa and were analyzed separately and in combination with elongation factor-1 alpha. Results from these analyses were concordant with those derived from elongation factor-1 alpha alone and provided support for a Hexapoda/Branchiopoda clade, thus arguing against the monophyly of the traditionally defined Atelocerata (Hexapoda + Myriapoda).

A new nuclear gene for insect phylogenetics: dopa decarboxylase is informative of relationships within Heliothinae (Lepidoptera: Noctuidae).

The lack of a readily accessible roster of nuclear genes informative at various taxonomic levels is a bottleneck for molecular systematics. In this report, we describe the first phylogenetic application of the sequence that encodes the enzyme dopa decarboxylase (DDC). For 14 test species within the noctuid moth subfamily Heliothinae that represent the previously best-supported groupings, a 690-bp fragment of DDC resolved relationships that are largely concordant with prior evidence from elongation factor-1 alpha (EF-1 alpha), morphology, and allozymes. Although both synonymous and nonsynonymous changes occur in DDC substantially more rapidly than they do in EF-1 alpha, DDC divergences within Heliothinae are below saturation at all codon positions. Analysis of DDC and EF-1 alpha in combination resulted in increased bootstrap support for several groupings. As a first estimate of previously unresolved relationships, DDC sequences were analyzed from 16 additional heliothines, for a total of 30 heliothine species plus outgroups. Previous relationships based on DDC were generally stable with increased taxon sampling, although a two- to eightfold downweighting of codon position 3 was required for complete concordance with the 14-species result. The weighted strict consensus trees were largely resolved and were congruent with most although not all previous hypotheses based on either morphology or EF-1 alpha. The proposed phylogeny suggests that the major agricultural pest heliothines belong to a single clade, characterized by polyphagy and associated life history traits, within this largely host-specific moth subfamily. DDC holds much promise for phylogenetic analysis of Tertiary-age animal groups.

Phylogenetic utility of elongation factor-1 alpha in noctuoidea (Insecta: Lepidoptera): the limits of synonymous substitution.

To test its phylogenetic utility, nucleotide sequence variation in a 1,240-bp fragment of the elongation factor-1 alpha (EF-1 alpha) gene was examined in 49 moth species representing the major groups of the superfamily Noctuoidea. Both parsimony and distance analyses supported the monophyly of nearly all groups for which there are clear morphological synapomorphies. Clades of subfamily rank and lower, probably mid-Tertiary and younger, were strongly supported. The third codon position contains 88% of variable sites, and approaches saturation at approximately 20% sequence divergence, possibly due to among-site rate heterogeneity and composition bias; higher divergences occur only in association with shifts in composition. Surprisingly, the few nonsynonymous changes appear no more phylogenetically reliable than synonymous changes. Signal strength for basal divergences is weak and fails to improve with character weighting; thus, dense taxon sampling is probably needed for strong inference from EF-1 alpha regarding deeper splits in Noctuoidea (probably early Tertiary). EF-1 alpha synonymous changes show promise for phylogeny reconstruction within Noctuidae and other groups of Tertiary age.

A nuclear gene for higher level phylogenetics: phosphoenolpyruvate carboxykinase tracks mesozoic-age divergences within Lepidoptera (Insecta).

The sequence of phosphoenolpyruvate carboxykinase (PEPCK) has been previously identified as a promising candidate for reconstructing Mesozoic-age divergences (Friedlander, Regier, and Mitter 1992, 1994). To test this hypothesis more rigorously, 597 nucleotides of aligned PEPCK coding sequence (approximately 30% of the coding region) were generated from 18 species representing Mesozoic-age lineages of moths (Insecta: Lepidoptera) and outgroup taxa. Relationships among basal Lepidoptera are well established by morphological analysis, providing a strong test for the utility of a gene which has not previously been used in systematics. Parsimony and other phylogenetic analyses were conducted on nucleotides by codon positions (nt1, nt2, nt3) separately and in combination, and on amino acids, for comparison to the test phylogeny. The highest concordance was achieved with nt1 + nt2, for which one of two most-parsimonious trees was identical to the test phylogeny, and with all nucleotides when nt3 was down-weighted sevenfold or higher, for which a single most-parsimonious tree identical to the test phylogeny resulted. Substitutions in nt3 approached saturation in many, but not all, pairwise comparisons and their exclusion or severe downweighting greatly increased the degree of concordance with the test phylogeny. Neighbor-joining analysis confirms this finding. The utility of PEPCK for phylogenetics is demonstrated over a time span for which few other suitable genes are currently available.

A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths.

Molecular systematists need increased access to nuclear genes. Highly conserved, low copy number protein-encoding nuclear genes have attractive features for phylogenetic inference but have heretofore been applied mostly to very ancient divergences. By virtue of their synonymous substitutions, such genes should contain a wealth of information about lower-level taxonomic relationships as well, with the advantage that amino acid conservatism makes both alignment and primer definition straightforward. We tested this postulate for the elongation factor-1 alpha (EF-1 alpha) gene in the noctuid moth subfamily Heliothinae, which has probably diversified since the middle Tertiary. We sequenced 1,240 bp in 18 taxa representing heliothine groupings strongly supported by previous morphological and allozyme studies. The single most parsimonious gene tree and the neighbor-joining tree for all nucleotides show almost complete concordance with the morphological tree. Homoplasy and pairwise divergence levels are low, transition/transversion ratios are high, and phylogenetic information is spread evenly across gene regions. The EF-1 alpha gene and presumably other highly conserved genes hold much promise for phylogenetics of Tertiary age eukaryote groups.

Evolution of chorion gene families in lepidoptera: characterization of 15 cDNAs from the gypsy moth.

Fifteen unique chorion protein-encoding cDNAs from gypsy moth have been completely sequenced. These sequences are encoded by a family of genes, based on pairwise similarity values of 78-100% within a 225-nt region. Pairwise comparisons and maximum parsimony analysis strongly support the existence of two clusters of 11 and four sequences each, called noc1 and noc2. While noc2 consists of two subclusters, there is little character support for subclusters within noc1. The highly localized character-state distribution on the parsimony tree in gypsy moth is reminiscent of that in Bombyx mori, specifically for those chorion families that have been shown to undergo gene conversion. Gene conversion thus becomes a reasonable explanation for the homogeneity of noc1 sequences and for their distinctness from noc2. The relationship between the two major clusters of chorion sequences in gypsy moth (noc1, noc2) and Bombyx mori (Bm alpha, Bm beta) has been addressed through mixed-species tree construction. All four groups cluster separately, thus providing no direct evidence of orthologous sequences. However, the occurrence of gene conversion could have eliminated such evidence. The relationship between the chorion gene tree and the species cladogenic event is discussed, as are biases in codon usage, base composition, and nucleotide transformations.

Loss of phylogenetic information in chorion gene families of Bombyx mori gene conversion.

The silkmoth chorion has provided a stimulating model for the study of evolution and developmental regulation of gene families. Previous attempts at inferring relationships among chorion sequences have been based on pairwise comparisons of overall similarity, a potentially problematic approach. To remedy this, we identified the alignable regions of low sequence variability and then analyzed this restricted database by parsimony and neighbor-joining methods. At the deepest level, the chorion sequence tree is split into two branches, called "alpha" and "beta." Within each branch, early- and late-expressing genes each constitute monophyletic groups, while the situation with middle-expressing genes remains uncertain. The HcB gene family appears to be the most basal beta-branch group, but this conclusion is qualified because the effect of gene conversion on branching order is unknown. Previous studies by Eickbush and colleagues have strongly suggested that ErA, HcA, and HcB families undergo gene conversion within a gene family, whereas the ErB family does not. The occurrence of conversion correlates with a particular tree structure; namely, branch lengths are much greater at the base of the family than at higher internodes and terminal branches. These observations raise the possibility that chorion gene families are defined by gene conversion events (reticulate evolution) rather than by descent with modification (synapomorphy).

Cell-specific expression in the silkmoth follicle: developmental characterization of a major chorion protein, its mRNA and gene.

Choriogenesis (eggshell formation) within the silkmoth Antheraea polyphemus proceeds in parallel for the two major subpopulations of follicle cells, diverging only during the very late period when aeropyle crown surface structures form in one region but not in the other. Correlated with their appearance is the synthesis of a set of region-specific proteins. In this report, aeropyle crowns are physically isolated and their protein composition is shown to consist of those same region-specific proteins. A cDNA clone, called pcvl 16, has been selected and shown to encode a lamellar-forming, aeropyle crown-specific protein, probably of the previously described C3,4 group. These conclusions are based on hybrid-selected translation, Northern analysis, and sequence analysis. pcvl 16 was used to isolate two distinct cloned copies of the 16 gene. Both 16 genes are closely paired with another region-specific gene but the proximity of the two gene pairs to each other is uncertain. Non-region-specific chorion genes expressed at earlier times in choriogenesis surround the 16 gene pairs, suggesting that cis sequences necessary for regionalized expression may be closely linked to coding sequences. To test this hypothesis, 5'-flanking sequences from eight region-specific genes are compared and shown to share two oligonucleotide sequences. One is a known regulatory element found in virtually all moth and fly chorion genes examined. The other, located just upstream from the TATA box, is not found in non-regionally expressed chorion genes and, thus, is a candidate for specifying regional expression.

Choriogenesis in the Lepidoptera: morphogenesis, protein synthesis, specific mRNA accumulation, and primary structure of a chorion cDNA from the gypsy moth.

Choriogenesis in silkmoths has been used as a model for relating developmental processes and evolutionary change. We now describe parallel studies in a nonsilkmoth lepidopteran species, the gypsy moth. Choriogenesis is described at the levels of ultrastructure, protein composition and synthesis, and specific mRNA accumulation. One complete chorion cDNA sequence is presented and features of its primary structure are discussed. This sequence is shown to be homologous to those of silkmoths.

Morphogenesis of silkmoth chorion: sequential modification of an early helicoidal framework through expansion and densification.

The silkmoth chorion is a helicoidally layered, fibrous structure which is constructed in four sequential morphogenetic modes, beginning with the assembly of a thin, low density lamellar framework. Subsequently, the framework expands in height by the insertion of additional fiber sheets into the preexisting lamellae. This expansion mode begins farthest from the follicular secretory cells and progresses in reverse. Individual fibers then grow in thickness, presumably through accretion of newly synthesized proteins, and eventually fuse. This third mode, which also begins in the most distant lamellae and proceeds in reverse, is called densification, as it results in an approximately two fold increase in overall chorion density without further lamellar expansion. Finally, lamellogenesis is recapitulated in miniature in a region of the chorion's surface, where very-late-forming lamellae are molded into prominent surface structures, the aeropyle crowns. The densification and especially the expansion modes suggest considerable fluidity in the developing chorion, consistent with its proposed cholesteric liquid crystalline structure. Such a structure is also consistent with numerous deviations from the ideal helicoidal array. These distortions and defects are described and discussed in terms of their possible origin and function.

Heterochrony and the introduction of novel modes of morphogenesis during the evolution of moth choriogenesis.

Choriogenesis in silkmoths (superfamily Bombycoidea) and in a sphingid moth (superfamily Sphingoidea) differ in major, but discrete, ways. In silkmoths, the predominant lamellar component assembles early in choriogenesis to form a thin framework. Subsequently, the lamellar framework is modified, first by expansion, and then by densification. Finally, ornate surface structures called aeropyle crowns form in some silkmoths, but they are absent in the species described here. In the sphingid, lamellar framework formation occurs throughout choriogenesis rather than largely during the early stages as in silkmoths. Lamellar densification occurs, but lamellar expansion and aeropyle crown formation do not. An evolutionary model is presented that accounts for the observed morphogenetic differences. Patterns of chorion protein synthesis in the sphingid differ from those in silkmoths in ways that are interpretable in light of the observed morphogenetic differences and the previously postulated functions of the proteins in silkmoths.

Evolutionary changes in the developmental expression of silkmoth chorion genes and their morphological consequences.

Discrete changes in silkmoth choriogenesis have occurred during evolution, as exemplified in the present report in Antheraea polyphemus and Hyalophora cecropia. At the level of morphology, the chorion of A. polyphemus has surface structures, called aeropyle crowns, that are absent from H. cecropia. Aeropyle crowns form during the very late period of choriogenesis and consist of two substructures--lamellae and filler. Filler is present in H. cecropia in greatly reduced amounts. At the level of protein synthesis, overall similarities in the two species are maintained until the very late period of choriogenesis, when synthesis of aeropyle crown components is maximal. In H. cecropia, very late period-specific proteins are reduced in number and abundance. Several of these minor proteins are candidates for E1 and E2, the components of filler. E1 and E2 RNAs are about 35 times more abundant in A. polyphemus, despite very similar gene copy numbers and times of expression in the two species. These results support the hypothesis that evolutionary changes in chorion morphology have resulted from regulatory changes in the expression of chorion genes, either at the level of transcription or mRNA decay. The hypothesis that evolutionary changes in chorion morphology are based on terminal addition onto a preexisting developmental program is discussed.