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.     

Dynamic Insertion–Deletion of Introns in Deuterostome EF-1α Genes

To test the validity of intron–exon structure as a phylogenetic marker, the intron–exon structure of EF-1α genes was investigated for starfish, acornworms, ascidians, larvaceans, and amphioxus and compared with that of vertebrates. Of the 11 distinct intron insertion sites found within the coding regions of the deuterostome EF-1α genes, 7 are shared by several taxa, while the remainder are unique to certain taxa. Examination of the shared introns of the deuterostome EF-1α gene revealed that independent intron loss or intron insertion must have occurred in separate lineages of the deuterostome taxa. Maximum parsimony analysis of the intron–exon data matrix recovered five parsimonious trees (consistency index = 0.867). From this result, we concluded that the intron–exon structure of deuterostome EF-1α has evolved more dynamically than previously thought, rendering it unsuitable as a phylogenetic marker. We also reconstructed an evolutionary history of intron insertion–deletion events on the deuterostome phylogeny, based on several molecular phylogenetic studies. These analyses revealed that the deuterostome EF-1α gene has lost individual introns more frequently than all introns simultaneously.     

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.      

Penicillium coffeae, a new endophytic species isolated from a coffee plant and its phylogenetic relationship to P. fellutanum, P. thiersii and P. brocae based on parsimony analysis of multilocus DNA sequences

Penicillium coffeae is described as a novel endophyte isolated from a Coffea arabica L. plant in Hawaii. The species is slow growing with short, vesiculate, monoverticillate conidiophores. Phylogenetic analysis using three loci shows that P. coffeae forms a strongly supported clade with P. fellutanum, P. charlesii, P. chermesinum, P. indicum, P. phoeniceum and P. brocae. Phenotypic ally these species are quite similar but can be distinguished. The EF-1alpha gene from P. fellutanum, P. charlesii, P. chermesinum and P. indicum lack introns, P. coffeae and P. phoeniceum have a previously unknown intron at codon 20 and P. brocae and P. thiersii isolates have a single intron at codon 26. The most parsimonious interpretation of intron changes on the strongly supported phylogenetic tree requires the gain of a novel intron at position 20 and loss of intron 26 to arrive at the current distribution of introns in this gene. This is one of only a few examples of intron gain in genes.     

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.     

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.      

Introduction of a nuclear marker for phylogenetic analysis of Nepenthaceae

Nepenthaceae, the pitcher plants of the Old World tropics show a remarkable diversity in SE Asia, especially on the islands of Borneo and Sumatra. This region is considered as a secondary center of diversity. Sequence analysis of the cpDNA TRNK intron supports this hypothesis showing the species of the Malay Archipelago as neighbour group to the isolated species from Sri Lanka, the Seychelles, and Madagascar. Based on phylogenetic reconstructions an origin of recent Nepenthaceae in the Indian subcontinent is assumed. A recent investigation focused on a non-plastid, translocated copy of the TRNK intron has revealed an incongruence to tree topology based on the cpDNA TRNK intron. Although the translocated copy emerged as insufficient for phylogenetic reconstruction of Nepenthaceae some taxa showed, contrary to the cpDNA dataset, relatively high distances to the rest of the taxa. These results indicated that the phylogeny of the TRNK intron could not reflect true phylogenetic relationships. We investigated the peptide transferase 1 (PTR1), to develop a phylogenetic marker that is based on a nuclear low copy gene in Nepenthes. All sequences obtained were probably functional, indicated by the ratio of point mutations of the single codon positions in exon and intron regions. Comparative analysis showed that this locus is of similar variability as the cpDNA TRNK intron and, contrary to the translocated copy of TRNK, potential useful for phylogenetic reconstruction. While in parts congruent to the plastid TRNK intron phylogeny, a higher divergence of some sequences in PRT1 and in the previously reported, non cpDNA dataset indicates that remnants of an older species stock persisted east of Wallace's line and on the Sunda Shelf. This suggests that plastid haplotypes existing today in the main distribution center of the Nepenthaceae could be descendants of more recently dispersed lineages that had been transmitted to an old species stock.     

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.     

Primers for the amplification of nuclear introns in sea stars of the family Asteriidae

We describe polymerase chain reaction primers that consistently amplify three intron regions of approximately 2 kb in total length for two nuclear protein‐coding genes (ATP synthase beta subunit and elongation factor‐1 alpha subunit) in sea stars of the family Asteriidae. The introns are moderately polymorphic at the species level (average within‐species percentage of site differences = 0.42%, range 0–1.44%), are evolving at about 29% as fast as mitochondrial sequences in the same species and are alignable at the genus or family level, making them suitable for phylogenetic and population genetic analyses.     

The phylogenetic analysis of variable-length sequence data: elongation factor-1alpha introns in European populations of the parasitoid wasp genus Pauesia (Hymenoptera: Braconidae: Aphidiinae)

Elongation factor-1alpha (EF-1alpha) is a highly conserved nuclear coding gene that can be used to investigate recent divergences due to the presence of rapidly evolving introns. However, a universal feature of intron sequences is that even closely related species exhibit insertion and deletion events, which cause variation in the lengths of the sequences. Indels are frequently rich in evolutionary information, but most investigators ignore sites that fall within these variable regions, largely because the analytical tools and theory are not well developed. We examined this problem in the taxonomically problematic parasitoid wasp genus Pauesia (Hymenoptera: Braconidae: Aphidiinae) using congruence as a criterion for assessing a range of methods for aligning such variable-length EF-1alpha intron sequences. These methods included distance- and parsimony-based multiple-alignment programs (CLUSTAL W and MALIGN), direct optimization (POY), and two "by eye" alignment strategies. Furthermore, with one method (CLUSTAL W) we explored in detail the robustness of results to changes in the gap cost parameters. Phenetic-based alignments ("by eye" and CLUSTAL W) appeared, under our criterion, to perform as well as more readily defensible, but computationally more demanding, methods. In general, all of our alignment and tree-building strategies recovered the same basic topological structure, which means that an underlying phylogenetic signal remained regardless of the strategy chosen. However, several relationships between clades were sensitive both to alignment and to tree-building protocol. Further alignments, considering only sequences belonging to the same group, allowed us to infer a range of phylogenetic relationships that were highly robust to tree-building protocol. By comparing these topologies with those obtained by varying the CLUSTAL parameters, we generated the distribution area of congruence and taxonomic compatibility. Finally, we present the first robust estimate of the European Pauesia phylogeny by using two EF-1alpha introns and 38 taxa (plus 3 outgroups). This estimate conflicts markedly with the traditional subgeneric classification. We recommend that this classification be abandoned, and we propose a series of monophyletic species groups.     

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.      

Use of the human elongation factor 1 alpha promoter as a versatile and efficient expression system

We have characterized the promoter region of the human elongation factor 1 alpha-encoding gene (EF-1 alpha) and developed a versatile expression system which has a wide host range and a high efficiency of gene expression. To identify the promoter region of the EF-1 alpha gene necessary for efficient gene expression, we constructed four pEF-CAT plasmids that have the bacterial cat gene fused to four different sites of the human EF-1 alpha gene: (i) ligated to the end of the TATA box (pEF220-CAT); (ii) ligated in exon 1 (pEF204-CAT and pEF233-CAT), and (iii) ligated in exon 2 (pEF321-CAT). All the pEF-CAT plasmids were highly expressed in all the cell types tested, including fibroblasts and lymphoid cells. Plasmid pEF321-CAT, which contains the first exon and the first intron, gave the highest level of cat expression. Plasmids pEF204- and pEF233-CAT, which contain part of the first exon but do not contain the first intron, were less efficient in cat expression than was pEF321-CAT. Plasmid pEF220-CAT, which lacks both the first exon and the first intron, was the least efficient. Plasmid pEF321-CAT was several- to 100-fold more efficient in cat expression than plasmid pSV2-CAT depending on the recipient cell types. The promoter of pEF321 plasmid also directed the stable expression of the bacterial neo gene more efficiently than the promoter of the simian virus 40 (SV40) early gene or the long terminal repeat of Rous sarcoma virus. Using this system, the SV40 early gene and the cDNA encoding human CD4 were also expressed efficiently.     

Two copies of the human apolipoprotein C-I gene are linked closely to the apolipoprotein E gene

The gene for human apolipoprotein (apo) C-I was selected from human genomic cosmid and lambda libraries. Restriction endonuclease analysis showed that the gene for apoC-I is located 5.5 kilobases downstream of the gene for apoE. A copy of the apoC-I gene, apoC-I', is located 7.5 kilobases downstream of the apoC-I gene. Both genes contain four exons and three introns; the apoC-I gene is 4653 base pairs long, the apoC-I' gene 4387 base pairs. In each gene, the first intron is located 20 nucleotides upstream from the translation start signal; the second intron, within the codon of Gly-7 of the signal peptide region; and the third intron, within the codon for Arg39 of the mature plasma protein coding region. The upstream apoC-I gene encodes the known apoC-I plasma protein and differs from the downstream apoC-I' gene in about 9% of the exon nucleotide positions. The most important difference between the exons results in a change in the codon for Gln-2 of the signal peptide region, which introduces a translation stop signal in the downstream gene. Major sequence differences are found in the second and third introns of the apoC-I and apoC-I' genes, which contain 9 and 7.5 copies, respectively, of Alu family sequences. The apoC-I gene is expressed primarily in the liver, and it is activated when monocytes differentiate into macrophages. In contrast, no mRNA product of the apoC-I' gene can be detected in any tissue, suggesting that it may be a pseudogene. The similar structures and the proximity of the apoE and apoC-I genes suggest that they are derived from a common ancestor. Furthermore, they may be considered to be constituents of a family of seven apolipoprotein genes (apoE, -C-I, -C-II, -C-III, -A-I, -A-II, and -A-IV) that have a common evolutionary origin.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

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.     

Phylogenetic analysis of Myriapoda using three nuclear protein-coding genes

We assessed the ability of three nuclear protein-encoding genes-elongation factor-1alpha (EF-1alpha), RNA polymerase II (Pol II), and elongation factor-2 (EF-2)-from 59 myriapod and 12 non-myriapod species to resolve phylogenetic relationships among myriapod classes and orders. In a previous study using EF-1alpha and Pol II (2134 nt combined) from 34 myriapod taxa, Regier and Shultz recovered widely accepted classes, orders, and families but failed to resolve interclass and interordinal relationships. The result was attributed to heterogenous rates of cladogenesis (specifically, the inability of the slowly evolving sequences to capture phylogenetic signal during rapid phylogenetic diversification) but the possibility of inadequate taxon sampling or limited sequence information could not be excluded. In the present study, the myriapod taxon sample was increased by 25 taxa (73%) and sequence length per taxon was effectively doubled through addition of EF-2 (4318 nt combined). Parsimony and Bayesian analyses of the expanded data set recovered a monophyletic Myriapoda, all four myriapod classes and all multiply sampled orders, often with high node support. However, except for three diplopod clades (Colobognatha, Helminothomorpha, and a subgroup of Pentazonia), few interordinal relationships and no interclass relationships were well supported. These results are similar to those of the earlier study by Regier and Shultz, which indicates that taxon sample and sequence length alone do not readily explain the weakly supported resolution in the earlier study. We review recent paleontological evidence to further develop our proposal that heterogeneity in phylogenetic signal provided by our slowly evolving sequences is due to heterogeneity in the temporal structure of myriapod diversification.     

蕨类植物叶绿体rps12基因的分子进化研究

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用HyPhy和PAML软件对该基因进行进化速率和适应性进化研究。结果显示:位于IR区的外显子2~3,其替换率明显降低,rps12基因编码序列的替换率也随之降低,且rps12基因密码子第3位的GC含量明显升高;在蕨类植物的进化过程中,3′-rps12更倾向定位于IR区,以保持较低的替换率;rps12基因编码的123个氨基酸位点中,共检测到4个正选择位点和116个负选择位点。研究结果表明基因序列进入到IR区后,显示出降低的替换率;强烈的负选择压力表明RPS12蛋白的高度保守性以及rps12基因的功能和结构已经趋于稳定。