Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

To test whether gaps resulting from sequence alignment contain phylogenetic signal concordant with those of base substitutions, we analyzed the occurrence of indel mutations upon a well-resolved, substitution-based tree for three nuclear genes in bumble bees (Bombus, Apidae: Bombini). The regions analyzed were exon and intron sequences of long-wavelength rhodopsin (LW Rh), arginine kinase (ArgK), and elongation factor-1alpha (EF-1alpha) F2 copy genes. LW Rh intron had only a few uninformative gaps, ArgK intron had relatively long gaps that were easily aligned, and EF-1alpha intron had many short gaps, resulting in multiple optimal alignments. The unambiguously aligned gaps within ArgK intron sequences showed no homoplasy upon the substitution-based tree, and phylogenetic signals within ambiguously aligned regions of EF-1alpha intron were highly congruent with those of base substitutions. We further analyzed the contribution of gap characters to phylogenetic reconstruction by incorporating them in parsimony analysis. Inclusion of gap characters consistently improved support for nodes recovered by substitutions, and inclusion of ambiguously aligned regions of EF-1alpha intron resolved several additional nodes, most of which were apical on the phylogeny. We conclude that gaps are an exceptionally reliable source of phylogenetic information that can be used to corroborate and refine phylogenies hypothesized by base substitutions, at least at lower taxonomic levels. At present, full use of gaps in phylogenetic reconstruction is best achieved in parsimony analysis, pending development of well-justified and generally applicable methods for incorporating indels in explicitly model-based methods.     

An evaluation of elongation factor 1 alpha as a phylogenetic marker for eukaryotes

Elongation factor 1 alpha (EF-1 alpha) is a highly conserved ubiquitous protein involved in translation that has been suggested to have desirable properties for phylogenetic inference. To examine the utility of EF-1 alpha as a phylogenetic marker for eukaryotes, we studied three properties of EF-1 alpha trees: congruency with other phyogenetic markers, the impact of species sampling, and the degree of substitutional saturation occurring between taxa. Our analyses indicate that the EF-1 alpha tree is congruent with some other molecular phylogenies in identifying both the deepest branches and some recent relationships in the eukaryotic line of descent. However, the topology of the intermediate portion of the EF-1 alpha tree, occupied by most of the protist lineages, differs for different phylogenetic methods, and bootstrap values for branches are low. Most problematic in this region is the failure of all phylogenetic methods to resolve the monophyly of two higher-order protistan taxa, the Ciliophora and the Alveolata. JACKMONO analyses indicated that the impact of species sampling on bootstrap support for most internal nodes of the eukaryotic EF-1 alpha tree is extreme. Furthermore, a comparison of observed versus inferred numbers of substitutions indicates that multiple overlapping substitutions have occurred, especially on the branch separating the Eukaryota from the Archaebacteria, suggesting that the rooting of the eukaryotic tree on the diplomonad lineage should be treated with caution. Overall, these results suggest that the phylogenies obtained from EF-1 alpha are congruent with other molecular phylogenies in recovering the monophyly of groups such as the Metazoa, Fungi, Magnoliophyta, and Euglenozoa. However, the interrelationships between these and other protist lineages are not well resolved. This lack of resolution may result from the combined effects of poor taxonomic sampling, relatively few informative positions, large numbers of overlapping substitutions that obscure phylogenetic signal, and lineage-specific rate increases in the EF-1 alpha data set. It is also consistent with the nearly simultaneous diversification of major eukaryotic lineages implied by the "big-bang" hypothesis of eukaryote evolution.     

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.      

Phylogenetic utility and evidence for multiple copies of elongation factor-1alpha in the spider genus Habronattus (Araneae: Salticidae)

In the continuing quest for informative genes for use in molecular systematics, the protein-coding gene Elongation factor-1alpha (EF-1alpha) has rapidly become one of the most prevalent "single-copy" nuclear genes utilized, particularly in arthropods. This paper explores the molecular evolutionary dynamics and phylogenetic utility of EF-1alpha in the salticid spider genus Habronattus. As has been reported for other arthropod lineages, our studies indicate that multiple (two) copies of EF-1alpha exist in Habronattus. These copies differ in intron structure and thus in size, making it possible to easily separate PCR amplification products. We present data for an intronless EF-1alpha copy for three Habronattus species. The presence of nonsense mutations and generally elevated rates of amino acid change suggest that this copy is evolving under relaxed functional constraints in Habronattus. A larger taxon sample (50 species plus outgroups) is presented for an EF-1alpha copy that includes both intron and exon regions. Characteristics of both regions suggest that this is a functional, orthologous copy in the species sampled. Maximum-likelihood relative-rate comparisons show that exon third codon sites are evolving more than 100 times as fast as second codon sites in these sequences and that intron sites are evolving about twice as fast as exon third sites. In combination, the EF-1alpha data provide robust, species-level phylogenetic signal that is largely congruent with morphologically well supported areas of Habronattus phylogeny. The recovery of some novel clades, and the unexpected fragmentation of others, suggests areas requiring further phylogenetic attention.     

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.     

Phylogeny of Trichoptera (caddisflies): characterization of signal and noise within multiple datasets

Trichoptera are holometabolous insects with aquatic larvae that, together with the Lepidoptera, make up the Amphiesmenoptera. Despite extensive previous morphological work, little phylogenetic agreement has been reached about the relationship among the three suborders--Annulipalpia, Spicipalpia, and Integripalpia--or about the monophyly of Spicipalpia. In an effort to resolve this conflict, we sequenced fragments of the large and small subunit nuclear ribosomal RNAs (1078 nt; D1, D3, V4-5), the nuclear elongation factor 1 alpha gene (EF-1 alpha; 1098 nt), and a fragment of mitochondrial cytochrome oxidase I (COI; 411 nt). Seventy adult and larval morphological characters were reanalyzed and added to molecular data in a combined analysis. We evaluated signal and homoplasy in each of the molecular datasets and attempted to rank the particular datasets according to how appropriate they were for inferring relationships among suborders. This evaluation included testing for conflict among datasets, comparing tree lengths among alternative hypotheses, measuring the left-skew of tree-length distributions from maximally divergent sets of taxa, evaluating the recovery of expected clades, visualizing whether or not substitutions were accumulating with time, and estimating nucleotide compositional bias. Although all these measures cast doubt on the reliability of the deep-level signal coming from the nucleotides of the COI and EF-1 alpha genes, these data could still be included in combined analyses without overturning the results from the most conservative marker, the rRNA. The different datasets were found to be evolving under extremely different rates. A site-specific likelihood method for dealing with combined data with nonoverlapping parameters was proposed, and a similar weighting scheme under parsimony was evaluated. Among our phylogenetic conclusions, we found Annulipalpia to be the most basal of the three suborders, with Spicipalpia and Integripalpia forming a clade. Monophyly of Annulipalpia and Integripalpia was confirmed, but the relationships among spicipalpians remain equivocal.     

Exploring data interaction and nucleotide alignment in a multiple gene analysis of Ips (Coleoptera: Scolytinae)

The possibility of gene tree incongruence in a species-level phylogenetic analysis of the genus Ips (Coleoptera: Scolytidae) was investigated based on mitochondrial 16S rRNA (16S) and nuclear elongation factor-1 alpha (EF-1 alpha) sequences, and existing cytochrome oxidase I (COI) and nonmolecular data sets. Separate cladistic analyses of the data partitions resulted in partially discordant most-parsimonious trees but revealed only low conflict of the phylogenetic signal. Interactions among data partitions, which differed in the extent of sequence divergence (COI > 16S > EF-1 alpha), base composition, and homoplasy, revealed that much of the branch support emerges only in the simultaneous analysis, particularly for deeper nodes in the tree, which are almost entirely supported through "hidden support" (sensu Gatesy et al., Cladistics 15:271-313, 1999). Apparent incongruence between data partitions is in part due to suboptimal alignments and bias of character transformations, but little evidence supports invoking incongruent phylogenetic histories of genetic loci. There is also no justification for eliminating or downweighting gene partitions on the basis of their apparent homoplasy or incongruence with other partitions, because the signal emerges only in the interaction of all data. In comparison with traditional taxonomy, the pini, plastographus, and perturbatus groups are polyphyletic, whereas the grandicollis group is monophyletic except for inclusion of the (monophyletic) calligraphus group. The latidens group and some European species are distantly related and closer to other genera within Ipini. Our robust cladogram was used to revise the classification of Ips. We provide new diagnoses for Ips and four subgeneric taxa.     

Unusually high evolutionary rate of the elongation factor 1 alpha genes from the Ciliophora and its impact on the phylogeny of eukaryotes

The elongation factor 1 alpha (EF-1 alpha) has become widely employed as a phylogenetic marker for studying eukaryotic evolution. However, a disturbing problem, the artifactual polyphyly of ciliates, is always observed. It has been suggested that the addition of new sequences will help to circumvent this problem. Thus, we have determined 15 new ciliate EF-1 alpha sequences, providing for a more comprehensive taxonomic sampling of this phylum. These sequences have been analyzed together with a representation of eukaryotic sequences using distance-, parsimony-, and likelihood-based phylogenetic methods. Such analyses again failed to recover the monophyly of Ciliophora. A study of the substitution rate showed that ciliate EF-1 alpha genes exhibit a high evolutionary rate, produced in part by an increased number of variable positions. This acceleration could be related to alterations of the accessory functions acquired by this protein, likely to those involving interactions with the cytoskeleton, which is very modified in the Ciliophora. The high evolutionary rate of these sequences leads to an artificial basal emergence of some ciliates in the eukaryotic tree by effecting a long-branch attraction artifact that produces an asymmetric topology for the basal region of the tree. The use of a maximum-likelihood phylogenetic method (which is less sensitive to long-branch attraction) and the addition of sequences to break long branches allow retrieval of more symmetric topologies, which suggests that the asymmetric part of the tree is most likely artifactual. Therefore, the sole reliable part of the tree appears to correspond to the apical symmetric region. These kinds of observations suggest that the general eukaryotic evolution might have consisted of a massive radiation followed by an increase in the evolutionary rates of certain groups that emerge artificially as early branches in the asymmetric base of the tree. Ciliates in the case of the EF-1 alpha genes would offer clear evidence for this hypothesis.     

Interaction of process partitions in phylogenetic analysis: an example from the swallowtail butterfly genus Papilio

In this study, we explored how the concept of the process partition may be applied to phylogenetic analysis. Sequence data were gathered from 23 species and subspecies of the swallowtail butterfly genus Papilio, as well as from two outgroup species from the genera Eurytides and Pachliopta. Sequence data consisted of 1,010 bp of the nuclear protein-coding gene elongation factor-1 alpha (EF-1 alpha) as well as the entire sequences (a total of 2,211 bp) of the mitochondrial protein-coding genes cytochrome oxidase I and cytochrome oxidase II (COI and COII). In order to examine the interaction between the nuclear and mitochondrial partitions in a combined analysis, we used a method of visualizing branch support as a function of partition weight ratios. We demonstrated how this method may be used to diagnose error at different levels of a tree in a combined maximum-parsimony analysis. Further, we assessed patterns of evolution within and between subsets of the data by implementing a multipartition maximum-likelihood model to estimate evolutionary parameters for various putative process partitions. COI third positions have an estimated average substitution rate more than 15 times that of EF-1 alpha, while COII third positions have an estimated average substitution rate more than 22 times that of EF-1 alpha. Ultimately, we found that although the mitochondrial and nuclear data were not significantly incongruent, homoplasy in the fast-evolving mitochondrial data confounded the resolution of basal relationships in the combined unweighted parsimony analysis despite the fact that there was relatively strong support for the relationships in the nuclear data. We conclude that there may be shortcomings to the methods of "total evidence" and "conditional combination" because they may fail to detect or accommodate the type of confounding bias we found in our data.     

Elongation factor-1 alpha occurs as two copies in bees: implications for phylogenetic analysis of EF-1 alpha sequences in insects

We report the complete sequence of a paralogous copy of elongation factor-1 alpha (EF-1 alpha) in the honeybee, Apis mellifera (Hymenoptera: Apidae). This copy differs from a previously described copy in the positions of five introns and in 25% of the nucleotide sites in the coding regions. The existence of two paralogous copies of EF-1 alpha in Drosophila and Apis suggests that two copies of EF-1 alpha may be widespread in the holometabolous insect orders. To distinguish between a single, ancient gene duplication and parallel, independent fly and bee gene duplications, we performed a phylogenetic analysis of hexapod EF-1 alpha sequences. Unweighted parsimony analysis of nucleotide sequences suggests an ancient gene duplication event, whereas weighted parsimony analysis of nucleotides and unweighted parsimony analysis of amino acids suggests the contrary: that EF-1 alpha underwent parallel gene duplications in the Diptera and the Hymenoptera. The hypothesis of parallel gene duplication is supported both by congruence among nucleotide and amino acid data sets and by topology-dependent permutation tail probability (T-PTP) tests. The resulting tree topologies are also congruent with current views on the relationships among the holometabolous orders included in this study (Diptera, Hymenoptera, and Lepidoptera). More sequences, from diverse orders of holometabolous insects, will be needed to more accurately assess the historical patterns of gene duplication in EF-1 alpha.      

Assessing phylogenetic resolution among mitochondrial, nuclear, and morphological datasets in Nothonotus darters (Teleostei: Percidae)

External morphological characters are the basis of our understanding of diversity and species relationships in many darter clades. The past decade has seen the publication of many studies utilizing mtDNA sequence data to investigate darter phylogenetics, but only recently have nuclear genes been used to investigate darter relationships. Despite a long tradition of use in darter systematics few studies have examined the phylogenetic utility of external morphological characters in estimating relationships among species in darter clades. We present DNA sequence data from the mitochondrial cytochrome b (cytb) gene, the nuclear encoded S7 intron 1, and discretely coded external morphological characters for all 20 species in the darter clade Nothonotus. Bayesian phylogenetic analyses result in phylogenies that are in broad agreement with previous studies. The cytb gene tree is well resolved, while the nuclear S7 gene tree lacks phylogenetic resolution, node support, and is characterized by a lack of reciprocal monophyly for many of the Nothonotus species. The phylogenies resulting from analysis of the morphological dataset lack resolution, but nodes present are found in the cytb and S7 gene trees. The highest resolution and node support is found in the Bayesian combined data phylogeny. Based on our results we propose continued exploration of the phylogenetic utility of external morphological characters in other darter clades. Given the extensive lack of reciprocal monophyly of species observed in the S7 gene tree we predict that nuclear gene sequences may have limited utility in intraspecific phylogeographic studies of Nothonotus darters.     

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

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

LDL binds to surface-expressed human T-cadherin in transfected HEK293 cells and influences homophilic adhesive interactions

The interaction between elongation factor 1alpha (EF-1alpha) and alpha/beta-tubulins has been analyzed in vivo and in vitro. An association of both alpha- and beta-tubulins with EF-1alpha in the lysate of Tetrahymena pyriformis was detected by co-immunoprecipitation analysis. In contrast, in vitro biomolecular interaction analysis with glutathione S-transferase (GST) fusion proteins revealed that GST-beta-tubulin, but not GST-alpha-tubulin, can bind to GST-EF-1alpha. Two beta-tubulin binding sites have been identified to reside in the domains I and III of EF-1alpha. In addition, beta-tubulin itself seems to have two distinct interaction sites for each of the domains. Since domain II of EF-1alpha did not interact with beta-tubulin, we have re-evaluated the phylogenetic status of ciliates using EF-1alpha sequences devoid of domain II. The phylogenetic tree thus obtained was significantly different from that inferred from the whole sequence of EF-1alpha, suggesting the presence of functional constraints on the molecular evolution of EF-1alpha.     

The primary structure and the functional domains of an elongation factor-1 alpha from Mucor racemosus

We have determined the complete nucleotide sequence for TEF-1, one of three genes coding for elongation factor (EF)-1 alpha in Mucor racemosus. The deduced EF-1 alpha protein contains 458 amino acids encoded by two exons. The presence of an intervening sequence located near the 3' end of the gene was predicted by the nucleotide sequence data and confirmed by alkaline S1 nuclease mapping. The amino acid sequence of EF-1 alpha was compared to the published amino acid sequences of EF-1 alpha proteins from Saccharomyces cerevisiae and Artemia salina. These proteins shared nearly 85% homology. A similar comparison to the functionally analogous EF-Tu from Escherichia coli revealed several regions of amino acid homology suggesting that the functional domains are conserved in elongation factors from these diverse organisms. Secondary structure predictions indicated that alpha helix and beta sheet conformations associated with the functional domains in EF-Tu are present in the same relative location in EF-1 alpha from M. racemosus. Through this comparative structural analysis we have predicted the general location of functional domains in EF-1 alpha which interact with GTP and tRNA.     

The phylogeny of acorn weevils (genus Curculio) from mitochondrial and nuclear DNA sequences: the problem of incomplete data

We considered the contribution of two mitochondrial and two nuclear data sets for the phylogenetic reconstruction of 22 species of seed beetles in the genus Curculio (Coleoptera: Cuculionidae). A phylogenetic tree from representatives found on various hosts was inferred from a combined data set of mitochondrial DNA cytochrome oxidase subunit I, mitochondrial cytochrome b, nuclear elongation factor 1alpha, and nuclear phosphoglycerate mutase, used for the first time as a molecular marker. Separate parsimony analyses of each data set showed that individual gene trees were mainly congruent and often complementary in the support of clades but the analysis was complicated by failure of PCR amplification of nuclear genes for many taxa and hence missing data entries. When the four gene partitions were combined in a simultaneous analysis despite the missing data, this increased the resolution and taxonomic coverage compared to the individual source trees. Alternative approaches of combining the information via supertree methodology produced a comparatively less resolved tree, and hence seem inferior to combining data matrices even in cases where numerous taxa are missing. The molecular data suggest a classification of the European species into two species groups that are in accordance with morphological characteristics but the data do no support any of the previously recognised American species groups.     

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.     

Molecular phylogeny of Zygomycota based on EF-1alpha and RPB1 sequences: limitations and utility of alternative markers to rDNA

Earlier molecular phylogenetic analyses based on nuclear small subunit ribosomal DNA (nSSU rDNA) suggest that the Zygomycota are polyphyletic within the Chytridiomycota. However, these analyses failed to resolve almost all interordinal relationships among basal fungi (Chytridiomycota and Zygomycota), due to lack of sufficient characters within the nSSU rDNA. To further elucidate the higher-level phylogeny of Zygomycota, we have sequenced partial RPB1 (DNA dependent RNA polymerase II largest subunit) and EF-1alpha (translation elongation factor 1 alpha) genes from 10 and 3 zygomycete fungi, respectively. Independent molecular phylogenetic analyses were performed based on each sequence by distance and maximum likelihood methods. Although deep phylogenetic relationships among basal fungi still remain poorly resolved using either gene, the RPB1-based phylogeny identified a novel monophyletic clade consisting of the Dimargaritales, Harpellales, and Kickxellales. This result suggests that regularly formed septa (cross walls that divide hyphae into segments) with a lenticular cavity are plesiomorphic for this clade, and indicates the importance of septal pore ultrastructure in zygomycete phylogeny. In addition, a peculiar mucoralean genus Mortierella, which was considered to be distantly related to the other Mucorales based on previous nSSU rDNA analyses, was resolved as the basal most divergence within the Mucorales, consistent with traditional phenotypic-based taxonomy. Although the taxa included in our analysis are restricted, the monophyly of each order suggested by nSSU rDNA phylogeny is supported by the present RPB1-based analysis. These results support the potential use of RPB1 as an alternative marker for fungal phylogenetic studies. Conversely, the overall fungal phylogeny based on EF-1alpha sequence is poorly resolved. A comparison of numbers of observed substitutions versus inferred substitutions within EF-1alpha indicates that this gene is much more saturated than RPB1. This result suggests that the EF-1alpha gene is unsuitable for resolving higher-level phylogenetic relationships within the Fungi.     

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.     

Tracing back EFL gene evolution in the cryptomonads-haptophytes assemblage: separate origins of EFL genes in haptophytes, photosynthetic cryptomonads, and goniomonads

A recently identified GTPase, elongation factor-like (EFL) protein is proposed to bear the principal functions of translation elongation factor 1alpha (EF-1alpha). Pioneering studies of EF-1alpha/EFL evolution have revealed the phylogenetically scattered distribution of EFL amongst eukaryotes, suggesting frequent eukaryote-to-eukaryote EFL gene transfer events and subsequent replacements of EF-1alpha functions by EFL. We here determined/identified seven new EFL sequences of the photosynthetic cryptomonad Cryptomonas ovata, the non-photosynthetic cryptomonad (goniomonad) Goniomonas amphinema, the foraminifer Planoglabratella opecularis, the haptophyte Chrysochromulina sp., the centroheliozoan Raphidiophrys contractilis, and two red algae Chondrus crispus and Gracilaria changii. The analyses of these EFL sequences successfully brought new insights into lateral EFL gene transfer amongst eukaryotes. Of most interest is a complex EFL evolution in a monophyletic assemblage comprised of cryptomonads and haptophytes. Since our analyses rejected any phylogenetic affinity amongst the EFL sequences from Goniomonas, photosynthetic cryptomonads, and haptophytes, the EFL genes of the three lineages most likely originated from different phylogenetic sources.     

Are Ichthyosporea animals or fungi? Bayesian phylogenetic analysis of elongation factor 1alpha of Ichthyophonus irregularis

Ichthyosporea is a recently recognized group of morphologically simple eukaryotes, many of which cause disease in aquatic organisms. Ribosomal RNA sequence analyses place Ichthyosporea near the divergence of the animal and fungal lineages, but do not allow resolution of its exact phylogenetic position. Some of the best evidence for a specific grouping of animals and fungi (Opisthokonta) has come from elongation factor 1alpha, not only phylogenetic analysis of sequences but also the presence or absence of short insertions and deletions. We sequenced the EF-1alpha gene from the ichthyosporean parasite Ichthyophonus irregularis and determined its phylogenetic position using neighbor-joining, parsimony and Bayesian methods. We also sequenced EF-1alpha genes from four chytrids to provide broader representation within fungi. Sequence analyses and the presence of a characteristic 12 amino acid insertion strongly indicate that I. irregularis is a member of Opisthokonta, but do not resolve whether I. irregularis is a specific relative of animals or of fungi. However, the EF-1alpha of I. irregularis exhibits a two amino acid deletion heretofore reported only among fungi.