Search

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host-virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts' genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Centipedes are among the most ancient groups of venomous predatory arthropods. Extant species belong to five orders, but our understanding of the composition and evolution of centipede venoms is based almost exclusively on one order, Scolopendromorpha. To gain a broader and less biased understanding we performed a comparative proteotranscriptomic analysis of centipede venoms from all five orders, including the first venom profiles for the orders Lithobiomorpha, Craterostigmomorpha, and Geophilomorpha. Our results reveal an astonishing structural diversity of venom components, with 93 phylogenetically distinct protein and peptide families. Proteomically-annotated gene trees of these putative toxin families show that centipede venom composition is highly dynamic across macroevolutionary timescales, with numerous gene duplications as well as functional recruitments and losses of toxin gene families. Strikingly, not a single family is found in the venoms of representatives of all five orders, with 67 families being unique for single orders. Ancestral state reconstructions reveal that centipede venom originated as a simple cocktail comprising just four toxin families, with very little compositional evolution happening during the approximately 50 My before the living orders had diverged. Venom complexity then increased in parallel within the orders, with scolopendromorphs evolving particularly complex venoms. Our results show that even venoms composed of toxins evolving under the strong constraint of negative selection can have striking evolutionary plasticity on the compositional level. We show that the functional recruitments and losses of toxin families that shape centipede venom arsenals are not concentrated early in their evolutionary history, but happen frequently throughout., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Activation of Rho-associated protein kinase 1 (ROCK1) and myotonic dystrophy kinase-related CDC42-binding kinase alpha (MRCKα) by caspases during apoptosis in vertebrates represents a prototypical example of co-option of kinases by proteases. How caspases acquired the ability to control these proteins during evolution of vertebrates is still unknown. Here, we report a phylogenetic and molecular study on the acquisition of caspase-cleavage sites in the family of Rho-activated kinases (RaKs). We demonstrate that the acquisition of such sites has more frequently occurred in identifiable intrinsically disordered regions (IDRs) within or flanking the coiled-coil domain. Thanks to computational identification of IDRs in protein sequences of different organisms, we predicted and validated the independent evolution of two caspase-cleavage sites in ROCK of arthropods and the loss of one of the MRCKα caspase-cleavage sites in ray-finned fishes. In conclusion, we shed light on the propensity of RaKs to evolve novel proteolytic sites, causing kinase activation and uniform subcellular distribution.

Even though lateral movements of transposons across families and even phyla within multicellular eukaryotic kingdoms have been found, little is known about transposon transfer between the kingdoms Animalia and Plantae. We discovered a novel non-LTR retrotransposon, AdLINE3, in a wild peanut species. Sequence comparisons and phylogenetic analyses indicated that AdLINE3 is a member of the RTE clade, originally identified in a nematode and rarely reported in plants. We identified RTE elements in 82 plants, spanning angiosperms to algae, including recently active elements in some flowering plants. RTE elements in flowering plants were likely derived from a single family we refer to as An-RTE. Interestingly, An-RTEs show significant DNA sequence identity with non-LTR retroelements from 42 animals belonging to four phyla. Moreover, the sequence identity of RTEs between two arthropods and two plants was higher than that of homologous genes. Phylogenetic and evolutionary analyses of RTEs from both animals and plants suggest that the An-RTE family was likely transferred horizontally into angiosperms from an ancient aphid(s) or ancestral arthropod(s). Notably, some An-RTEs were recruited as coding sequences of functional genes participating in metabolic or other biochemical processes in plants. This is the first potential example of horizontal transfer of transposons between animals and flowering plants. Our findings help to understand exchanges of genetic material between the kingdom Animalia and Plantae and suggest arthropods likely impacted on plant genome evolution., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Chemosensory-related gene (CRG) families have been studied extensively in insects, but their evolutionary history across the Arthropoda had remained relatively unexplored. Here, we address current hypotheses and prior conclusions on CRG family evolution using a more comprehensive data set. In particular, odorant receptors were hypothesized to have proliferated during terrestrial colonization by insects (hexapods), but their association with other pancrustacean clades and with independent terrestrial colonizations in other arthropod subphyla have been unclear. We also examine hypotheses on which arthropod CRG family is most ancient. Thus, we reconstructed phylogenies of CRGs, including those from new arthropod genomes and transcriptomes, and mapped CRG gains and losses across arthropod lineages. Our analysis was strengthened by including crustaceans, especially copepods, which reside outside the hexapod/branchiopod clade within the subphylum Pancrustacea. We generated the first high-resolution genome sequence of the copepod Eurytemora affinis and annotated its CRGs. We found odorant receptors and odorant binding proteins present only in hexapods (insects) and absent from all other arthropod lineages, indicating that they are not universal adaptations to land. Gustatory receptors likely represent the oldest chemosensory receptors among CRGs, dating back to the Placozoa. We also clarified and confirmed the evolutionary history of antennal ionotropic receptors across the Arthropoda. All antennal ionotropic receptors in E. affinis were expressed more highly in males than in females, suggestive of an association with male mate-recognition behavior. This study is the most comprehensive comparative analysis to date of CRG family evolution across the largest and most speciose metazoan phylum Arthropoda., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Wolbachia bacteria infect about half of all arthropods, with diverse and extreme consequences ranging from sex-ratio distortion and mating incompatibilities to protection against viruses. These phenotypic effects, combined with efficient vertical transmission from mothers to offspring, satisfactorily explain the invasion dynamics of Wolbachia within species. However, beyond the species level, the lack of congruence between the host and symbiont phylogenetic trees indicates that Wolbachia horizontal transfers and extinctions do happen and underlie its global distribution. But how often do they occur? And has the Wolbachia pandemic reached its equilibrium? Here, we address these questions by inferring recent acquisition/loss events from the distribution of Wolbachia lineages across the mitochondrial DNA tree of 3,600 arthropod specimens, spanning 1,100 species from Tahiti and surrounding islands. We show that most events occurred within the last million years, but are likely attributable to individual level variation (e.g., imperfect maternal transmission) rather than population level variation (e.g., Wolbachia extinction). At the population level, we estimate that mitochondria typically accumulate 4.7% substitutions per site during an infected episode, and 7.1% substitutions per site during the uninfected phase. Using a Bayesian time calibration of the mitochondrial tree, these numbers translate into infected and uninfected phases of approximately 7 and 9 million years. Infected species thus lose Wolbachia slightly more often than uninfected species acquire it, supporting the view that its present incidence, estimated here slightly below 0.5, represents an epidemiological equilibrium., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Transposable elements (TEs) and other repetitive DNA can accumulate in the absence of recombination, a process contributing to the degeneration of Y-chromosomes and other nonrecombining genome portions. A similar accumulation of repetitive DNA is expected for asexually reproducing species, given their entire genome is effectively nonrecombining. We tested this expectation by comparing the whole-genome TE loads of five asexual arthropod lineages and their sexual relatives, including asexual and sexual lineages of crustaceans (Daphnia water fleas), insects (Leptopilina wasps), and mites (Oribatida). Surprisingly, there was no evidence for increased TE load in genomes of asexual as compared to sexual lineages, neither for all classes of repetitive elements combined nor for specific TE families. Our study therefore suggests that nonrecombining genomes do not accumulate TEs like nonrecombining genomic regions of sexual lineages. Even if a slight but undetected increase of TEs were caused by asexual reproduction, it appears to be negligible compared to variance between species caused by processes unrelated to reproductive mode. It remains to be determined if molecular mechanisms underlying genome regulation in asexuals hamper TE activity. Alternatively, the differences in TE dynamics between nonrecombining genomes in asexual lineages versus nonrecombining genome portions in sexual species might stem from selection for benign TEs in asexual lineages because of the lack of genetic conflict between TEs and their hosts and/or because asexual lineages may only arise from sexual ancestors with particularly low TE loads., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Insects are an important model for the study of innate immune systems, but remarkably little is known about the immune system of other arthropod groups despite their importance as disease vectors, pests, and components of biological diversity. Using comparative genomics, we have characterized the immune system of all the major groups of arthropods beyond insects for the first time--studying five chelicerates, a myriapod, and a crustacean. We found clear traces of an ancient origin of innate immunity, with some arthropods having Toll-like receptors and C3-complement factors that are more closely related in sequence or structure to vertebrates than other arthropods. Across the arthropods some components of the immune system, such as the Toll signaling pathway, are highly conserved. However, there is also remarkable diversity. The chelicerates apparently lack the Imd signaling pathway and beta-1,3 glucan binding proteins--a key class of pathogen recognition receptors. Many genes have large copy number variation across species, and this may sometimes be accompanied by changes in function. For example, we find that peptidoglycan recognition proteins have frequently lost their catalytic activity and switch between secreted and intracellular forms. We also find that there has been widespread and extensive duplication of the cellular immune receptor Dscam (Down syndrome cell adhesion molecule), which may be an alternative way to generate the high diversity produced by alternative splicing in insects. In the antiviral short interfering RNAi pathway Argonaute 2 evolves rapidly and is frequently duplicated, with a highly variable copy number. Our results provide a detailed analysis of the immune systems of several important groups of animals for the first time and lay the foundations for functional work on these groups., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Orthology inference is central to phylogenomic analyses. Phylogenomic data sets commonly include transcriptomes and low-coverage genomes that are incomplete and contain errors and isoforms. These properties can severely violate the underlying assumptions of orthology inference with existing heuristics. We present a procedure that uses phylogenies for both homology and orthology assignment. The procedure first uses similarity scores to infer putative homologs that are then aligned, constructed into phylogenies, and pruned of spurious branches caused by deep paralogs, misassembly, frameshifts, or recombination. These final homologs are then used to identify orthologs. We explore four alternative tree-based orthology inference approaches, of which two are new. These accommodate gene and genome duplications as well as gene tree discordance. We demonstrate these methods in three published data sets including the grape family, Hymenoptera, and millipedes with divergence times ranging from approximately 100 to over 400 Ma. The procedure significantly increased the completeness and accuracy of the inferred homologs and orthologs. We also found that data sets that are more recently diverged and/or include more high-coverage genomes had more complete sets of orthologs. To explicitly evaluate sources of conflicting phylogenetic signals, we applied serial jackknife analyses of gene regions keeping each locus intact. The methods described here can scale to over 100 taxa. They have been implemented in python with independent scripts for each step, making it easy to modify or incorporate them into existing pipelines. All scripts are available from https://bitbucket.org/yangya/phylogenomic_dataset_construction., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Relationships between the five extant orders of centipedes have been considered solved based on morphology. Phylogenies based on samples of up to a few dozen genes have largely been congruent with the morphological tree apart from an alternative placement of one order, the relictual Craterostigmomorpha, consisting of two species in Tasmania and New Zealand. To address this incongruence, novel transcriptomic data were generated to sample all five orders of centipedes and also used as a test case for studying gene-tree incongruence. Maximum likelihood and Bayesian mixture model analyses of a data set composed of 1,934 orthologs with 45% missing data, as well as the 389 orthologs in the least saturated, stationary quartile, retrieve strong support for a sister-group relationship between Craterostigmomorpha and all other pleurostigmophoran centipedes, of which the latter group is newly named Amalpighiata. The Amalpighiata hypothesis, which shows little gene-tree incongruence and is robust to the influence of among-taxon compositional heterogeneity, implies convergent evolution in several morphological and behavioral characters traditionally used in centipede phylogenetics, such as maternal brood care, but accords with patterns of first appearances in the fossil record., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The DNA damage response (DDR) is a crucial signaling network that preserves the integrity of the genome. This network is an ensemble of distinct but often overlapping subnetworks, where different components fulfill distinct functions in precise spatial and temporal scenarios. To understand how these elements have been assembled together in humans, we performed comparative genomic analyses in 47 selected species to trace back their emergence using systematic phylogenetic analyses and estimated gene ages. The emergence of the contribution of posttranslational modifications to the complex regulation of DDR was also investigated. This is the first time a systematic analysis has focused on the evolution of DDR subnetworks as a whole. Our results indicate that a DDR core, mostly constructed around metabolic activities, appeared soon after the emergence of eukaryotes, and that additional regulatory capacities appeared later through complex evolutionary process. Potential key posttranslational modifications were also in place then, with interacting pairs preferentially appearing at the same evolutionary time, although modifications often led to the subsequent acquisition of new targets afterwards. We also found extensive gene loss in essential modules of the regulatory network in fungi, plants, and arthropods, important for their validation as model organisms for DDR studies.

The orphan nuclear receptor gene knirps and its relatives encode a small family of highly conserved proteins. We take advantage of the conservation of the family, using the recent prevalence of genomic data, to reconstruct its evolutionary history, identifying duplication events and tracing the intron-exon structure of the genes over evolution. Many arthropod species have two or three members of this family, but the orthology between members is unclear. We have analyzed the protein coding sequences of members of this family from 15 arthropod species covering all four main arthropod classes, including a total of 28 genes. All members of the family encode a highly conserved 94 amino acid core sequence, part of which is encoded by a single invariant exon. We find that many of the automated predictions of these genes contain errors, while some copies of the gene were not uncovered by automated pipelines, requiring manual corrections and curation. We use the coding sequences to present a phylogenetic analysis of the knirps family. Our analysis indicates that there was a duplication of a single ancestral gene in the lineage leading to insects, which gave rise to two paralogs, eagle and knirps-related. Descendants of this duplication can be identified by the presence or absence of a short protein-coding motif. Independent, lineage-specific duplications occurred in the two crustaceans we sampled. Within the insects, the knirps-related gene underwent further lineage-specific duplications, giving rise to--among others--the Drosophila gap gene knirps.

Inaccurate inference of positional homologies in multiple sequence alignments and systematic errors introduced by alignment heuristics obfuscate phylogenetic inference. Alignment masking, the elimination of phylogenetically uninformative or misleading sites from an alignment before phylogenetic analysis, is a common practice in phylogenetic analysis. Although masking is often done manually, automated methods are necessary to handle the much larger data sets being prepared today. In this study, we introduce the concept of subsplits and demonstrate their use in extracting phylogenetic signal from alignments. We design a clustering approach for alignment masking where each cluster contains similar columns-similarity being defined on the basis of compatible subsplits; our approach then identifies noisy clusters and eliminates them. Trees inferred from the columns in the retained clusters are found to be topologically closer to the reference trees. We test our method on numerous standard benchmarks (both synthetic and biological data sets) and compare its performance with other methods of alignment masking. We find that our method can eliminate sites more accurately than other methods, particularly on divergent data, and can improve the topologies of the inferred trees in likelihood-based analyses. Software available upon request from the author.

Multiple visual pigments, prerequisites for color vision, are found in arthropods, but the evolutionary origin of their diversity remains obscure. In this study, we explore the opsin genes in five distantly related species of Onychophora, using deep transcriptome sequencing and screening approaches. Surprisingly, our data reveal the presence of only one opsin gene (onychopsin) in each onychophoran species, and our behavioral experiments indicate a maximum sensitivity of onychopsin to blue-green light. In our phylogenetic analyses, the onychopsins represent the sister group to the monophyletic clade of visual r-opsins of arthropods. These results concur with phylogenomic support for the sister-group status of the Onychophora and Arthropoda and provide evidence for monochromatic vision in velvet worms and in the last common ancestor of Onychophora and Arthropoda. We conclude that the diversification of visual pigments and color vision evolved in arthropods, along with the evolution of compound eyes-one of the most sophisticated visual systems known.

The rate of change in DNA is an important parameter for understanding molecular evolution and hence for inferences drawn from studies of phylogeography and phylogenetics. Most rate calibrations for mitochondrial coding regions in marine species have been made from divergence dating for fossils and vicariant events older than 1-2 My and are typically 0.5-2% per lineage per million years. Recently, calibrations made with ancient DNA (aDNA) from younger dates have yielded faster rates, suggesting that estimates of the molecular rate of change depend on the time of calibration, decaying from the instantaneous mutation rate to the phylogenetic substitution rate. aDNA methods for recent calibrations are not available for most marine taxa so instead we use radiometric dates for sea-level rise onto the Sunda Shelf following the Last Glacial Maximum (starting ∼18,000 years ago), which led to massive population expansions for marine species. Instead of divergence dating, we use a two-epoch coalescent model of logistic population growth preceded by a constant population size to infer a time in mutational units for the beginning of these expansion events. This model compares favorably to simpler coalescent models of constant population size, and exponential or logistic growth, and is far more precise than estimates from the mismatch distribution. Mean rates estimated with this method for mitochondrial coding genes in three invertebrate species are elevated in comparison to older calibration points (2.3-6.6% per lineage per million years), lending additional support to the hypothesis of calibration time dependency for molecular rates.

The wealth of available genomic data presents an unrivaled opportunity to study the molecular basis of evolution. Studies on gene family expansions and site-dependent analyses have already helped establish important insights into how proteins facilitate adaptation. However, efforts to conduct full-scale cross-genomic comparisons between species are challenged by both growing amounts of data and the inherent difficulty in accurately inferring homology between deeply rooted species. Proteins, in comparison, evolve by means of domain rearrangements, a process more amenable to study given the strength of profile-based homology inference and the lower rates with which rearrangements occur. However, adapting to a constantly changing environment can require molecular modulations beyond reach of rearrangement alone. Here, we explore rates and functional implications of novel domain emergence in contrast to domain gain and loss in 20 arthropod species of the pancrustacean clade. Emerging domains are more likely disordered in structure and spread more rapidly within their genomes than established domains. Furthermore, although domain turnover occurs at lower rates than gene family turnover, we find strong evidence that the emergence of novel domains is foremost associated with environmental adaptation such as abiotic stress response. The results presented here illustrate the simplicity with which domain-based analyses can unravel key players of nature's adaptational machinery, complementing the classical site-based analyses of adaptation.

Retroelements represent a considerable fraction of many eukaryotic genomes and are considered major drives for adaptive genetic innovations. Recent discoveries showed that despite not normally using DNA intermediates like retroviruses do, Mononegaviruses (i.e., viruses with nonsegmented, negative-sense RNA genomes) can integrate gene fragments into the genomes of their hosts. This was shown for Bornaviridae and Filoviridae, the sequences of which have been found integrated into the germ line cells of many vertebrate hosts. Here, we show that Rhabdoviridae sequences, the major Mononegavirales family, have integrated only into the genomes of arthropod species. We identified 185 integrated rhabdoviral elements (IREs) coding for nucleoproteins, glycoproteins, or RNA-dependent RNA polymerases; they were mostly found in the genomes of the mosquito Aedes aegypti and the blacklegged tick Ixodes scapularis. Phylogenetic analyses showed that most IREs in A. aegypti derived from multiple independent integration events. Since RNA viruses are submitted to much higher substitution rates as compared with their hosts, IREs thus represent fossil traces of the diversity of extinct Rhabdoviruses. Furthermore, analyses of orthologous IREs in A. aegypti field mosquitoes sampled worldwide identified an integrated polymerase IRE fragment that appeared under purifying selection within several million years, which supports a functional role in the host's biology. These results show that A. aegypti was subjected to repeated Rhabdovirus infectious episodes during its evolution history, which led to the accumulation of many integrated sequences. They also suggest that like retroviruses, integrated rhabdoviral sequences may participate actively in the evolution of their hosts.

The use of secondary structures has been advocated to improve both the alignment and the tree reconstruction processes of ribosomal RNA (rRNA) data sets. We used simulated and empirical rRNA data to test the impact of secondary structure consideration in both steps of molecular phylogenetic analyses. A simulation approach was used to generate realistic rRNA data sets based on real 16S, 18S, and 28S sequences and structures in combination with different branch length and topologies. Alignment and tree reconstruction performance of four recent structural alignment methods was compared with exclusively sequence-based approaches. As empirical data, we used a hexapod rRNA data set to study the influence of nucleotide interdependencies in sequence alignment and tree reconstruction. Structural alignment methods delivered significantly better sequence alignments compared with pure sequence-based methods. Also, structural alignment methods delivered better trees judged by topological congruence to simulation base trees. However, the advantage of structural alignments was less pronounced and even vanished in several instances. For simulated data, application of mixed RNA/DNA models to stems and loops, respectively, led to significantly shorter branches. The application of mixed RNA/DNA models in the hexapod analyses delivered partly implausible relationships. This can be interpreted as a stronger sensitivity of mixed model setups to nonphylogenetic signal. Secondary structure consideration clearly influenced sequence alignment and tree reconstruction of ribosomal genes. Although sequence alignment quality can considerably be improved by the use of secondary structure information, the application of mixed models in tree reconstructions needs further studies to understand the observed effects.

Arthropods were the first animals to conquer land and air. They encompass more than three quarters of all described living species. This extraordinary evolutionary success is based on an astoundingly wide array of highly adaptive body organizations. A lack of robustly resolved phylogenetic relationships, however, currently impedes the reliable reconstruction of the underlying evolutionary processes. Here, we show that phylogenomic data can substantially advance our understanding of arthropod evolution and resolve several conflicts among existing hypotheses. We assembled a data set of 233 taxa and 775 genes from which an optimally informative data set of 117 taxa and 129 genes was finally selected using new heuristics and compared with the unreduced data set. We included novel expressed sequence tag (EST) data for 11 species and all published phylogenomic data augmented by recently published EST data on taxonomically important arthropod taxa. This thorough sampling reduces the chance of obtaining spurious results due to stochastic effects of undersampling taxa and genes. Orthology prediction of genes, alignment masking tools, and selection of most informative genes due to a balanced taxa-gene ratio using new heuristics were established. Our optimized data set robustly resolves major arthropod relationships. We received strong support for a sister group relationship of onychophorans and euarthropods and strong support for a close association of tardigrades and cycloneuralia. Within pancrustaceans, our analyses yielded paraphyletic crustaceans and monophyletic hexapods and robustly resolved monophyletic endopterygote insects. However, our analyses also showed for few deep splits that were recently thought to be resolved, for example, the position of myriapods, a remarkable sensitivity to methods of analyses.

In contrast to conventional splicing, which joins exons from a single primary transcript, trans-splicing links stretches of RNA from separate transcripts, derived from distinct regions of the genome. Spliced leader (SL) trans-splicing is particularly well known in trypanosomes, nematodes, and flatworms, where it provides messenger RNAs with a leader sequence and cap that allow them to be translated efficiently. One of the largest puzzles regarding SL trans-splicing is its evolutionary origin. Until now SL trans-splicing has been found in a small and disparate set of organisms (including trypanosomes, dinoflagellates, cnidarians, rotifers, nematodes, flatworms, and urochordates) but not in most other eukaryotic lineages, including well-studied groups such as fungi, plants, arthropods, and vertebrates. This patchy distribution could either suggest that trans-splicing was present in early eukaryotes/metazoans and subsequently lost in multiple lineages or that it evolved several times independently. Starting from the serendipitous discovery of SL trans-splicing in an arthropod, we undertook a comprehensive survey of this process in the animal kingdom. By surveying expressed sequence tag data from more than 70 metazoan species, we show that SL trans-splicing also occurs in at least two groups of arthropods (amphipod and copepod crustaceans), in ctenophores, and in hexactinellid sponges. However, we find no evidence for SL trans-splicing in other groups of arthropods and sponges or in 15 other phyla that we have surveyed. Although the presence of SL trans-splicing in hydrozoan cnidarians, hexactinellid sponges, and ctenophores might suggest that it was present at the base of the Metazoa, the patchy distribution that is evident at higher resolution suggests that SL trans-splicing has evolved repeatedly among metazoan lineages. In agreement with this scenario, we discuss evidence that SL precursor RNAs can readily evolve from ubiquitous small nuclear RNAs that are used for conventional splicing.

Studies on color vision in invertebrates have focused primarily on insect visual pigments, with little attention given to crustacean visual pigments. None of the blue-green-, blue-, or ultraviolet (UV)-sensitive-opsins have been identified in crustaceans. In addition, the discussion of visual pigments has been limited to long-wavelength-sensitive opsins in Pancrustacea. Here, we focused on Branchiopoda (Crustacea), which is a sister group of Hexapoda including insects. In the tadpole shrimp Triops granarius, the visual pigment chromophore was retinal. Multiple opsins were isolated from each of three branchiopod species, T. granarius, Triops longicaudatus, and the fairy shrimp Branchinella kugenumaensis (five, five, and four opsins from these species, respectively). Phylogenetic analyses and the presence of a lysine residue corresponding to position 90 in bovine rhodopsin suggested that three of the branchiopod opsins comprise UV-sensitive pigments. In addition, the phylogenetic relationships between insect and branchiopod UV-sensitive opsins revealed that the divergence of blue- and UV-sensitive pigments predates the Branchiopoda and Insecta divergence. The other branchiopod opsins show distant relationships to other known insect opsins and form novel clusters. The present results strongly suggest that the ancestral arthropod of the Chelicerata-Pancrustacea lineages possessed at least four types of opsins. The ancestors of Pancrustacea and the Insecta-Branchiopoda lineages possessed at least five and six types of opsins, respectively. Our results suggest that in the evolutionary process associated with each lineage, several opsins appeared and diversified with repeated gene duplication, of which some have been lost in some taxa.

The 7SK small nuclear RNA (snRNA) is a key player in the regulation of polymerase (pol) II transcription. The 7SK RNA was long believed to be specific to vertebrates where it is highly conserved. Homologs in basal deuterostomes and a few lophotrochozoan species were only recently reported. On longer timescales, 7SK evolves rapidly with only few conserved sequence and structure motifs. Previous attempts to identify the Drosophila homolog thus have remained unsuccessful despite considerable efforts. Here we report on the discovery of arthropod 7SK RNAs using a novel search strategy based on pol III promoters, as well as the subsequent verification of its expression. Our results demonstrate that a 7SK snRNA featuring 2 highly structured conserved domains was present already in the bilaterian ancestor.

We combined the category (CAT) mixture model (Lartillot N, Philippe H. 2004) and the nonstationary break point (BP) model (Blanquart S, Lartillot N. 2006) into a new model, CAT-BP, accounting for variations of the evolutionary process both along the sequence and across lineages. As in CAT, the model implements a mixture of distinct Markovian processes of substitution distributed among sites, thus accommodating site-specific selective constraints induced by protein structure and function. Furthermore, as in BP, these processes are nonstationary, and their equilibrium frequencies are allowed to change along lineages in a correlated way, through discrete shifts in global amino acid composition distributed along the phylogenetic tree. We implemented the CAT-BP model in a Bayesian Markov Chain Monte Carlo framework and compared its predictions with those of 3 simpler models, BP, CAT, and the site- and time-homogeneous general time-reversible (GTR) model, on a concatenation of 4 mitochondrial proteins of 20 arthropod species. In contrast to GTR, BP, and CAT, which all display a phylogenetic reconstruction artifact positioning the bees Apis mellifera and Melipona bicolor among chelicerates, the CAT-BP model is able to recover the monophyly of insects. Using posterior predictive tests, we further show that the CAT-BP combination yields better anticipations of site- and taxon-specific amino acid frequencies and that it better accounts for the homoplasies that are responsible for the artifact. Altogether, our results show that the joint modeling of heterogeneities across sites and along time results in a synergistic improvement of the phylogenetic inference, indicating that it is essential to disentangle the combined effects of both sources of heterogeneity, in order to overcome systematic errors in protein phylogenetic analyses.

There is growing interest in the use of alternative, more slowly-evolving RGCs (rare genomic changes). Recently, Rogozin and coauthors (Rogozin et al. 2007) proposed a novel phylogenetic method employing rare amino acid changes, RGC-CAMs (rare genomic changes-conserved amino acids-multiple substitutions). They applied their method to 694 sets of eukaryotic orthologs in order to distinguish the relationship between nematodes, arthropods and deuterostomes. They concluded that such rare amino acid changes were consistent with the Coelomata hypothesis, which groups arthropods and deuterostomes to the exclusion of nematodes. Here we use newly available genomic sequences from Nematostella vectensis, a basal metazoan, and from Brugia malayi, an additional nematode. We show that the apparent support for Coelomata is likely to be the result of the rapid rate of evolution leading to Caenorhabditis nematodes. Including the additional species paints a very different picture, with 13 remaining characters consistent with Ecdysozoa versus only 1 consistent with Coelomata.

As the number of sequenced genomes from diverse walks of life rapidly increases, phylogenetic analysis is entering a new era: reconstruction of the evolutionary history of organisms on the basis of full-scale comparison of their genomes. In addition to brute force, genome-wide analysis of alignments, rare genomic changes (RGCs) that are thought to comprise derived shared characters of individual clades are increasingly used in genome-wide phylogenetic studies. We propose a new type of RGCs designated RGC_CAMs (after Conserved Amino acids-Multiple substitutions), which are inferred using a genome-scale analysis of protein and underlying nucleotide sequence alignments. The RGC_CAM approach utilizes amino acid residues conserved in major eukaryotic lineages, with the exception of a few species comprising a putative clade, and selects for phylogenetic inference only those amino acid replacements that require 2 or 3 nucleotide substitutions, in order to reduce homoplasy. The RGC_CAM analysis was combined with a procedure for rigorous statistical testing of competing phylogenetic hypotheses. The RGC_CAM method is shown to be robust to branch length differences and taxon sampling. When applied to animal phylogeny, the RGC_CAM approach strongly supports the coelomate clade that unites chordates with arthropods as opposed to the ecdysozoan (molting animals) clade. This conclusion runs against the view of animal evolution that is currently prevailing in the evo-devo community. The final solution to the coelomate-ecdysozoa controversy will require a much larger set of complete genome sequences representing diverse animal taxa. It is expected that RGC_CAM and other RGC-based methods will be crucial for these future, definitive phylogenetic studies.

R2 is a non-long-terminal-repeat (LTR) retrotransposon that inserts specifically into 28S rDNA. R2 has been identified in many species of arthropods and three species of chordates. R2 may be even more widely distributed in animals, and its origin may be traceable to early animal evolution. In this study, we identified R2 elements in medaka fish, White Cloud Mountain minnow, Reeves' turtle, hagfish, sea lilies, and some arthropod species, using degenerate polymerase chain reaction methods. We also identified two R2 elements from the public genomic sequence database of the bloodfluke Schistosoma mansoni. One of the two bloodfluke R2 elements has two zinc-finger motifs at the N-terminus; this differs from other known R2 elements, which have one or three zinc-finger motifs. Phylogenetic analysis revealed that the whole phylogeny of R2 can be divided into 11 parts (subclades), in which the local R2 phylogeny and the corresponding host phylogeny are consistent. Divergence-versus-age analysis revealed that there is no reliable evidence for the horizontal transfer of R2 but supports the proposition that R2 has been vertically transferred since before the divergence of the deuterostomes and protostomes. The seeming inconsistency between the R2 phylogeny and the phylogeny of their hosts is due to the existence of paralogous lineages. The number of N-terminal zinc-finger motifs is consistent with the deep phylogeny of R2 and indicates that the common ancestor of R2 had three zinc-finger motifs at the N-terminus. This study revealed the long-term vertical inheritance and the ancient origin of sequence specificity of R2, both of which seem applicable to some other non-LTR retrotransposons.

This study combined complete 18S with partial 28S ribosomal RNA gene sequences ( approximately 2,000 nt in total) to investigate the relations of basal hexapods. Ten species of Protura, 12 of Diplura, and 10 of Collembola (representing all subgroups of these three clades) were sequenced, along with 5 true insects and 8 other arthropods, which served as out-groups. Trees were constructed with maximum parsimony, maximum likelihood, Bayesian analysis, and minimum-evolution analysis of LogDet-transformed distances. All methods yielded strong support for a clade of Protura plus Diplura, here named Nonoculata, and for monophyly of the Diplura. Parametric-bootstrapping analysis showed our data to be inconsistent with previous hypotheses (P < 0.01) that joined Protura with Collembola (Ellipura), that said Diplura are sister to true insects or are diphyletic, and that said Collembola are not hexapods. That is, our data are consistent with hexapod monophyly and Collembola grouped weakly with "Protura + Diplura" under most analytical conditions. As a caveat to the above conclusions, the sequences showed nonstationarity of nucleotide frequencies across taxa, so the CG-rich sequences of the diplurans and proturans may have grouped together artifactually; however, the fact that the LogDet method supported this group lessens this possibility. Within the basal hexapod groups, where nucleotide frequencies were stationary, traditional taxonomic subgroups generally were recovered: i.e., within Protura, the Eosentomata and Acerentomata (but Sinentomata was not monophyletic); within Collembola, the Arthropleona, Poduromorpha, and Entomobryomorpha (but Symphypleona was polyphyletic); and in Diplura, the most complete data set (> 2,100 nt) showed monophyly of Campodeoidea and of Japygoidea, and most methods united Projapygoidea with Japygoidea.

Eukaryotic translation initiation factor 2 (eIF2) is a G protein that delivers the methionyl initiator tRNA to the small ribosomal subunit and releases it upon GTP hydrolysis after the recognition of the initiation codon. eIF2 is composed of three subunits, alpha, beta, and gamma. Subunit gamma shows the strongest conservation, and it confers both tRNA and GTP/GDP binding. Using intron positioning and protein sequence alignment, here we show that eIF2gamma is a suitable phylogenetic marker for eukaryotes. We determined or completed the sequences of 13 arthropod eIF2gamma genes. Analyzing the phylogenetic distribution of 52 different intron positions in 55 distantly related eIF2gamma genes, we identified ancient ones and shared derived introns in our data set. Obviously, intron positioning in eIF2gamma is evolutionarily conserved. However, there were episodes of complete and partial intron losses followed by intron gains. We identified 17 clusters of intron positions based on their distribution. The evolution of these clusters appears to be connected with preferred exon length and can be used to estimate the relative timing of intron gain because nearby precursor introns had to be erased from the gene before the new introns could be inserted. Moreover, we identified a putative case of intron sliding that constitutes a synapomorphic character state supporting monophyly of Coleoptera, Lepidoptera, and Diptera excluding Hymenoptera. We also performed tree reconstructions using the eIF2gamma protein sequences and intron positioning as phylogenetic information. Our results support the monophyly of Viridoplantae, Ascomycota, Homobasidiomyceta, and Apicomplexa.

Recent advances in molecular phylogenetics are continuously changing our perception of the phylogenetic relationships among the main arthropod lineages: crustaceans, hexapods, chelicerates, and myriapods. Besides the intrinsic interest in unraveling the evolution of the largest animal phylum, these studies are basic to an understanding of one of the major transitions in animal evolution-i.e., the conquest of land with all its associated structural and functional adaptations. Myriapods have been traditionally considered the closest relatives of hexapods, thus implying only one origin of terrestriality for the tracheate lineage, but this view is now challenged by molecular evidence. Sequence data available to date for centipedes and millipedes are very limited, and the taxon sampling is strongly biased. The most critical gap was the scutigeromorph centipedes, which are the sister group to all remaining Chilopoda from which they probably diverged in the Silurian if not earlier. We obtained the first complete mitochondrial sequence for a representative of this clade, the house centipede. In our phylogenetic analyses of the protein-coding genes in this mitochondrial genome, along with 16 further ones representing the other major arthropod clades plus two outgroups, the myriapods formed a clade with the chelicerates. This implies that water-to-land transition occurred at least three times (hexapods, myriapods, arachnids) during the evolution of the Arthropoda. In addition, in contrast to all previous studies, our best supported topologies favor paraphyly of the myriapods with respect to the chelicerates. This would increase to four the main events of land colonization in arthropods (once for centipedes, once for millipedes).

We report the isolation and sequencing of genomic copies of mariner transposons involved in recent horizontal transfers into the genomes of the European earwig, Forficula auricularia; the European honey bee, Apis mellifera; the Mediterranean fruit fly, Ceratitis capitata; and a blister beetle, Epicauta funebris, insects from four different orders. These elements are in the mellifera subfamily and are the second documented example of full-length mariner elements involved in this kind of phenomenon. We applied maximum likelihood methods to the coding sequences and determined that the copies in each genome were evolving neutrally, whereas reconstructed ancestral coding sequences appeared to be under selection, which strengthens our previous hypothesis that the primary selective constraint on mariner sequence evolution is the act of horizontal transfer between genomes.

Rearrangements of mitochondrial DNA gene order have been suggested as a tool for defining the pattern of evolutionary divergence in arthropod taxa. We have employed a combination of highly conserved insect-based polymerase chain reaction (PCR) primers with long-PCR to survey 14 noninsect arthropods for mitochondrial gene rearrangements. The size of the amplified fragments was used to order the primer containing genes. Five chelicerates exhibit amplicons that are consistent with the presumptive ancestral arthropod mtDNA gene order. These five species comprise two soft ticks, two prostriate hard ticks, and an opilionid. Six other chelicerates, all metastriate hard ticks, have a different arrangement that was originally discovered by this procedure and has been previously detailed in a complete mtDNA sequence. Three new arthropod mtDNA gene arrangements are described here. They were discovered in a terrestrial crustacean (Isopoda) and two myriapods (Chilopoda, centipede; Diplopoda, millipede). These rearrangements include major realignments of some of the large coding regions and two possible new positions for the tRNA(Met) (M) gene in arthropods. The long-PCR approach affords an opportunity to quickly screen divergent taxa for major rearrangements. Taxa exhibiting rearrangements can be targeted for DNA sequencing of gene boundaries to establish the details of the mtDNA organization.

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.

We determined the complete mitochondrial DNA (mtDNA) sequences of the millipedes Narceus annularus and Thyropygus sp. (Arthropoda: Diplopoda) and identified, in both genomes, all 37 genes typical for metazoan mtDNA. The arrangement of these genes is identical in the two millipedes, but differs from others found in arthropod mtDNAs in the location of at least four genes or gene blocks. This novel gene arrangement is unusual for animal mtDNA in that genes with identical transcriptional polarities are clustered in the genome, and the two clusters are separated by two noncoding regions. The only exception to this pattern is the gene for cysteine tRNA, which is located in the part of the genome that otherwise contains all genes with the opposite transcriptional polarity. We suggest that a mechanism involving complete mtDNA duplication followed by the loss of genes, predetermined by their transcriptional polarity and location in the genome, could generate this gene arrangement from the one ancestral for arthropods. The proposed mechanism has important implications for phylogenetic inferences that are drawn on the basis of gene arrangement comparisons.

Hemocyanins are copper-containing respiratory proteins of the Arthropoda that have so far been thoroughly investigated only in the Chelicerata and the Crustacea but have remained unstudied until now in the Myriapoda. Here we report the first sequence of a myriapod hemocyanin. The hemocyanin of Spirostreptus sp. (Diplopoda: Spirostreptidae) is composed of two distinct subunits that are arranged in a 6 x 6 native molecule. The cloned hemocyanin subunit cDNA codes of for a polypeptide of 653 amino acids (75.5 kDa) that includes a signal peptide of 18 amino acids. The sequence closely resembles that of the chelicerate hemocyanins. Molecular phylogenetic analyses reject with high statistical confidence the integrity of the Tracheata (i.e., Myriapoda + Insecta) but give conflicting results on the position of the myriapod hemocyanin. While distance matrix and maximum-likelihood methods support a basal position of the Spirostreptus hemocyanin with respect to the other hemocyanins, parsimony analysis suggests a sister group relationship with the chelicerate hemocyanins. The latter topology is also supported by a unique shared deletion of an alpha-helix. A common ancestry of Myriapoda and Chelicerata should be seriously considered.

Arthropod hemocyanins are members of a protein superfamily that also comprises the arthropod phenoloxidases (tyrosinases), crustacean pseudohemocyanins (cryptocyanins), and insect storage hexamerins. The evolution of these proteins was inferred by neighbor-joining, maximum-parsimony, and maximum-likelihood methods. Monte Carlo shuffling approaches provided evidence against a discernible relationship of the arthropod hemocyanin superfamily and molluscan hemocyanins or nonarthropodan tyrosinases. Within the arthropod hemocyanin superfamily, the phenoloxidase probably emerged early in the (eu-)arthropod stemline and thus form the most likely outgroup. The respiratory hemocyanins evolved from these enzymes before the radiation of the extant euarthropodan subphyla. Due to different functional constraints, replacement rates greatly vary between the clades. Divergence times were thus estimated assuming local molecular clocks using several substitution models. The results were consistent and indicated the separation of the cheliceratan and crustacean hemocyanins close to 600 MYA. The different subunit types of the multihexameric cheliceratan hemocyanin have a rather conservative structure and diversified in the arachnidan stemline between 550 and 450 MYA. By contrast, the separation of the crustacean (malacostracan) hemocyanin subunits probably occurred only about 200 MYA. The nonrespiratory pseudohemocyanins evolved within the Decapoda about 215 MYA. The insect hemocyanins and storage hexamerins emerged independently from the crustacean hemocyanins. The time of divergence of the insect proteins from the malacostracan hemocyanins was estimated to be about 430-440 MYA, providing support for the notion that the Hexapoda evolved from the same crustacean lineage as the Malacostraca.

We report a contiguous region of more than half (> 7,500 nt) of the mitochondrial genomes for Platynereis dumerii (Annelida: Polychaeta), Helobdella robusta (Annelida: Hirudinida), and Galathealinum brachiosum (Pogonophora: Perviata). The relative arrangements of all 22 genes identified for Helobdella and Galathealinum are identical to one another and to their arrangements in the mtDNA of the previously studied oligochaete annelid Lumbricus. In contrast, Platynereis differs from these taxa in the positions of several tRNA genes and in having two additional tRNA genes (trnC and trnM) and a large noncoding sequence in this region. Comparisons of relative gene arrangements and of the nucleotide and inferred amino acid sequences among these and other published taxa provide strong support for an annelid-mollusk clade that excludes arthropods, and for the inclusion of pogonophorans within Annelida, rather than giving them separate phylum status. Gene arrangement comparisons include the first use of a recently described method on previously unpublished data. Although a variety of alternative initiation codons are typically used by mitochondrial protein-encoding genes, ATG appears to be the initiator for all but one reported here. The large noncoding region (1,091 nt) identified in Platynereis has no significant sequence similarity to the noncoding region of Lumbricus, although each contains runs of TA dinucleotides and of homopolymers, which could potentially serve as signaling elements. There is strong bias for synonymous codon usage in Helobdella and especially in Galathealinum. In this latter taxon, 5 codons are completely unused, 13 are used three or fewer times, and G appears at third codon positions in only 26 of the 2,236 codons. Nucleotide composition bias appears to influence amino acid composition of the proteins.

R2 elements are non-LTR retrotransposons that insert in the 28S rRNA genes of arthropods. Partial sequence data from many species have previously suggested that these elements have been vertically inherited since the origin of this phylum. Here, we compare the complete sequences of nine R2 elements selected to represent the diversity of arthropods. All of the elements exhibited a uniform structure. Identification of their conserved sequence features, combined with our biochemical studies, allows us to make the following inferences concerning the retrotransposition mechanism of R2. While all R2 elements insert into the identical sequence of the 28S gene, it is only the location of the initial nick in the target DNA that is rigidly conserved across arthropods. Variation at the R2 5' junctions suggests that cleavage of the second strand of the target site is not conserved within or between species. The extreme 5' and 3' ends of the elements themselves are also poorly conserved, consistent with a target primed reverse transcription mechanism for attachment of the 3' end and a template switch model for the attachment of the 5' end. Comparison of the approximately 1,000-aa R2 ORF reveals that it can be divided into three domains. The central 450-aa domain can be folded by homology modeling into a tertiary structure resembling the fingers, palm, and thumb subdomains of retroviral reverse transcriptases. The carboxyl terminal end of the R2 protein appears to be the endonuclease domain, while the amino-terminal end contains zinc finger and c-myb-like DNA-binding motifs.

The entire mitochondrial genome was sequenced in a prostriate tick, Ixodes hexagonus, and a metastriate tick, Rhipicephalus sanguineus. Both genomes encode 22 tRNAs, 13 proteins, and two ribosomal RNAs. Prostriate ticks are basal members of Ixodidae and have the same gene order as Limulus polyphemus. In contrast, in R. sanguineus, a block of genes encoding NADH dehydrogenase subunit 1 (ND1), tRNA(Leu)(UUR), tRNA(Leu)(CUN), 16S rDNA, tRNA(Val), 12S rDNA, the control region, and the tRNA(Ile) and tRNA(Gln) have translocated to a position between the tRNA(Glu) and tRNA(Phe) genes. The tRNA(Cys) gene has translocated between the control region and the tRNA(Met) gene, and the tRNA(Leu)(CUN) gene has translocated between the tRNA(Ser)(UCN) gene and the control region. Furthermore, the control region is duplicated, and both copies undergo concerted evolution. Primers that flank these rearrangements confirm that this gene order is conserved in all metastriate ticks examined. Correspondence analysis of amino acid and codon use in the two ticks and in nine other arthropod mitochondrial genomes indicate a strong bias in R. sanguineus towards amino acids encoded by AT-rich codons.

To understand the early evolution of the Metazoa, it is necessary to determine the correct phylogenetic status of diploblastic animals. Despite cladistic studies of morphological characters and recent molecular phylogenetic studies, it remains uncertain whether diploblasts are monophyletic or paraphyletic, and how the phyla of diploblasts are phylogenetically related. The heat shock protein 70 (Hsp70) sequences, because of their ubiquity and high degree of conservation, could provide a useful model for phylogenetic analysis. We have sequenced almost the entire nucleic acid sequence of cytoplasmic Hsp70 from eight diploblastic species. Our data support the monophyly of diploblastic animals. However, the phylogenetic relationships of the diploblast groups were not significantly resolved. Our phylogenetic trees also support the monophyly of Metazoa with high bootstrap values, indicating that animals form an extremely robust clade.

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

Numerous complete mitochondrial DNA sequences have been determined for species within two arthropod groups, insects and crustaceans, but there are none for a third, the chelicerates. Most mitochondrial gene arrangements reported for crustaceans and insect species are identical or nearly identical to that of Drosophila yakuba. Sequences across 36 of the gene boundaries in the mitochondrial DNA (mtDNA) of a representative chelicerate. Limulus polyphemus L., also reveal an arrangement like that of Drosophila yakuba. Only the position of the tRNA(LEU)(UUR) gene differs; in Limulus it is between the genes for tRNA(LEU)(CUN) and ND1. This positioning is also found in onychophorans, mollusks, and annelids, but not in insects and crustaceans, and indicates that tRNA(LEU)(CUN)-tRNA(LEU)(UUR)-ND1 was the ancestral gene arrangement for these groups, as suggested earlier. There are no differences in the relative arrangements of protein-coding and ribosomal RNA genes between Limulus and Drosophila, and none have been observed within arthropods. The high degree of similarity of mitochondrial gene arrangements within arthropods is striking, since some taxa last shared a common ancestor before the Cambrian, and contrasts with the extensive mtDNA rearrangements occasionally observed within some other metazoan phyla (e.g., mollusks and nematodes).

To determine a possible relationship between organismal and molecular evolution, the divergence patterns of gene families were examined by taking special notice of functional difference, tissue distribution, and intracellular localization of the members. A phylogenetic analysis of 25 different gene families revealed interesting patterns of divergence of these families: Most gene duplications giving rise to different functions antedate the vertebrates-arthropods separation. On the other hand, in a group of members carrying virtually identical function to one another but differing in tissue distribution (tissue-specific isoform), most gene duplications have occurred independently in each of vertebrates and arthropods after the separation of the two animal groups. In family members encoding molecules localizing in cell compartments (compartmentalized isoforms), the gene duplications antedate the animals-fungi separation. In the cases of the Ca2+ pump and rab subfamilies, the compartmentalized isoforms were shown to have diverged during the early evolution of eukaryotes. A phylogenetic analysis of the tissue-specific isoforms from 26 different subfamilies revealed extensive gene duplications and rapid rates of amino acid substitutions in the early evolution of chordates before the separation of fishes and tetrapods. On the contrary, the genetic variations are relatively low in the later period. This pattern of evolution observed at the molecular level is correlated well with that of tissue evolution based on fossil evidence and morphological data, and thus evolution at the two levels may be related.

The complete 18S rDNA gene sequence of Macrobiotus group hufelandi (Tardigrada) was obtained and aligned with 18S rDNA and rRNA gene sequences of 24 metazoans (mainly protostomes). Discrete character (maximum-parsimony) and distance (neighbor-joining) methods were used to infer their phylogeny. The evolution of bootstrap proportions with sequence length (pattern of resolved nodes, PRN) was studied to test the resolution of the nodes in neighbor-joining trees. The results show that arthropods are monophyletic. Tardigrades represent the sister group of arthropods (in parsimony analyses) or they are related with crustaceans (distance analysis and PRN). Arthropoda are divided into two main evolutionary lines, the Hexapoda + Crustacea line (weakly supported), and the Myriapoda + Chelicerata line. The Hexapoda + Crustacea line includes Pentastomida, but the internal resolution is far from clear. The Insecta (Ectognatha) are monophyletic, but no evidence for the monophyly of Hexapoda is found. The Chelicerata are a monophyletic group and the Myriapoda cluster close to Arachnida. Overall, the results obtained represent the first molecular evidence for a Tardigrada + Arthropoda clade. In addition, the congruence between molecular phylogenies of the Arthropoda from other authors and this obtained here indicates the need to review those obtained solely on morphological characters.

Complete coding regions of the 18S rRNA gene of an enteropneust hemichordate and an echinoid and ophiuroid echinoderm were obtained and aligned with 18S rRNA gene sequences of all major chordate clades and four outgroups. Gene sequences were analyzed to test morphological character phylogenies and to assess the strength of the signal. Maximum-parsimony analysis of the sequences fails to support a monophyletic Chordata; the urochordates form the sister taxon to the hemichordates, and together this clade plus the echinoderms forms the sister taxon to the cephalochordates plus craniates. Decay, bootstrap, and tree-length distribution analyses suggest that the signal for inference of dueterostome phylogeny is weak in this molecule. Parsimony analysis of morphological plus molecular characters supports both monophyly of echinoderms plus enteropneust hemichordates and a sister group relationship of this clade to chordates. Evolutionary parsimony does not support chordate monophyly. Neighbor-joining, Fitch-Margoliash, and maximum-likelihood analyses support a chordate lineage that is the sister group to an echinoderm-plus-hemichordate lineage. The results illustrate both the limitations of the 18S rRNA molecule alone for high-level phylogeny inference and the importance of considering both molecular and morphological data in phylogeny reconstruction.

Crustacean and cheliceratan hemocyanins (oxygen-transport proteins) and insect hexamerins (storage proteins) are homologous gene products, although the latter do not bind oxygen and do not possess the copper-binding histidines present in the hemocyanins. An alignment of 19 amino acid sequences of hemocyanin subunits and insect hexamerins was made, based on the conservation of elements of secondary structure observed in X-ray structures of two hemocyanin subunits. The alignment was analyzed using parsimony and neighbor-joining methods. Results provide strong indications for grouping together the sequences of the 2 crustacean hemocyanin subunits, the 5 cheliceratan hemocyanin subunits, and the 12 insect hexamerins. Within the insect clade, four methionine-rich proteins, four arylphorins, and two juvenile hormone-suppressible proteins from Lepidoptera, as well as two dipteran proteins, form four separate groups. In the absence of an outgroup sequence, it is not possible to present information about the ancestral state from which these proteins are derived. Although this family of proteins clearly consists of homologous gene products, there remain striking differences in gene organization and site of biosynthesis of the proteins within the cell. Because studies on 18S and 12S rRNA sequences indicate a rather close relationship between insects and crustaceans, we propose that hemocyanin is the ancestral arthropod protein and that insect hexamerins lost their copper-binding capability after divergence of the insects from the crustaceans.

The polymerase chain reaction and hybridization to genomic blots were used to investigate whether the previously observed inclusion of 5S ribosomal RNA genes in the 28S-18S ribosomal DNA intergenic regions of some crustacean species (copepods) could also be detected in other arthropods. Such an arrangement was found not only in other calanoid copepod species but also in a cirriped, an euphausid, and a spider. It is interesting that species from two different calanoid copepod genera do not have this type of arrangement. We conclude that the inclusion of 5S ribosomal RNA genes within the ribosomal DNA repeats has probably occurred repeatedly during the evolution of arthropod species and that the mechanism(s) responsible for these insertions could also be responsible for their loss.