Molecular systematics of sciurognathi (rodentia): the mitochondrial cytochrome b and 12S rRNA genes support the Anomaluroidea (Pedetidae and Anomaluridae)

Nucleotide sequence data from the mitochondrial 12S rRNA and cytochrome b genes were used to analyze phylogenetic relationships among sciurognath rodents. Our sample taxa included representatives of 11 sciurognath and 3 hystricognath families with two marsupial species, Didelphis virginiana and Macropus robustus, as outgroups. The dataset was analyzed using both maximum-parsimony (weighted and unweighted) and likelihood methods. Three suprafamilial groupings are strongly supported: Geomyidae + Heteromyidae (Geomyoidea), Sciuridae + Aplodontidae (Sciuroidea), and Pedetidae + Anomaluridae (Anomaluroidea). Although moderately supported, two sister group relationships were identified between Gliridae and Sciuroidea and between Castor and Geomyoidea. In contrast to previous nuclear DNA evidence, the evolutionary affinities between Ctenodactylidae and Hystricognathi (Ctenohystrica) and between Muridae and Dipodidae (Myodonta) are not supported by the mitochondrial data. Molecular divergence dates based on the combined data were estimated for suprafamilial groupings and are discussed in the light of current morphological and paleontological interpretations of rodent phylogeny.     

The molecular taxonomy and evolution of the guinea pig.

On the basis of 18 protein sequences totaling 2,413 aligned amino acid sites, it is suggested that the guinea pigs and the myomorphs (rat-like rodents) are not monophyletic. Rather, the evolutionary lineage leading to the guinea pig seems to have branched off prior to the divergence among myomorphs, lagomorphs, primates, chiropterans, artiodactyls, and carnivores. It is suggested therefore that the Caviomorpha (guinea pig-like rodents) and possibly the Hystricomorpha (porcupine-like rodents) should be elevated in taxonomic rank and conferred an ordinal status distinct from the Rodentia. This suggestion calls for a reevaluation of the morphological evolution of guinea pigs and further molecular studies on the possibility of paraphyly of the order Rodentia. If the monophyly of rodents holds, it must be concluded that the pattern of molecular evolution in many guinea pig genes has been extremely unusual and that the causes for this pattern should be sought. It is also suggested that claims of large differences in the rate of molecular evolution between guinea pigs and myomorphs may have been exaggerated in many cases as a result of an erroneous phylogenetic position for the guinea pig. The average rate of amino acid replacement in the guinea pig seems to be comparable to that in the rat and the mouse. However, the data indicate that myomorph and caviomorph genes evolve, on average, about two times faster than their human counterparts. Finally, our analysis provides evidence against the hypothesis that the gundi (an African rodent) represents the most ancient rodent lineage.     

Molecular evolution of the mitochondrial 12S rRNA in Ungulata (mammalia)

The complete 12S rRNA gene has been sequenced in 4 Ungulata (hoofed eutherians) and 1 marsupial and compared to 38 available mammalian sequences in order to investigate the molecular evolution of the mitochondrial small-subunit ribosomal RNA molecule. Ungulata were represented by one artiodactyl (the collared peccary, Tayassu tajacu, suborder Suiformes), two perissodactyls (the Grevy's zebra, Equus grevyi, suborder Hippomorpha; the white rhinoceros, Ceratotherium simum, suborder Ceratomorpha), and one hyracoid (the tree hyrax, Dendrohyrax dorsalis). The fifth species was a marsupial, the eastern gray kangaroo (Macropus giganteus). Several transition/transversion biases characterized the pattern of changes between mammalian 12S rRNA molecules. A bias toward transitions was found among 12S rRNA sequences of Ungulata, illustrating the general bias exhibited by ribosomal and protein-encoding genes of the mitochondrial genome. The derivation of a mammalian 12S rRNA secondary structure model from the comparison of 43 eutherian and marsupial sequences evidenced a pronounced bias against transversions in stems. Moreover, transversional compensatory changes were rare events within double-stranded regions of the ribosomal RNA. Evolutionary characteristics of the 12S rRNA were compared with those of the nuclear 18S and 28S rRNAs. From a phylogenetic point of view, transitions, transversions and indels in stems as well as transversional and indels events in loops gave congruent results for comparisons within orders. Some compensatory changes in double-stranded regions and some indels in single-stranded regions also constituted diagnostic events. The 12S rRNA molecule confirmed the monophyly of infraorder Pecora and order Cetacea and demonstrated the monophyly of suborder Suiformes. However, the monophyly of the suborder Ruminantia was not supported, and the branching pattern between Cetacea and the artiodactyl suborders Ruminantia and Suiformes was not established. The monophyly of the order Perissodactyla was evidenced, but the relationships between Artiodactyla, Cetacea, and Perissodactyla remained unresolved. Nevertheless, we found no support for a Perissodactyla + Hyracoidea clade, neither with distance approach, nor with parsimony reconstruction. The 12S rRNA was useful to solve intraordinal relationships among Ungulata, but it seemed to harbor too few informative positions to decipher the bushlike radiation of some Ungulata orders, an event which has most probably occurred in a short span of time between 55 and 70 MYA. Correspondence to: E. Douzery     

Nucleotide sequence variation in the mitochondrial 12S rRNA gene and the phylogeny of African mole-rats (Rodentia: Bathyergidae).

Mitochondrial DNA (mtDNA) sequence variation was examined in eight taxa of the African rodent family Bathyergidae, as well as in two taxa representative of the Old-World hystricognathid rodent families Petromyidae and Thryonomyidae. A total of 812 bp, constituting domains I-III of the 12S ribosomal rRNA gene, were compared for each taxon. The high levels of intrafamilial mtDNA sequence divergence observed (average 16.8, range 3.5-23.2) support an ancient origin for the five genera, 20-38 Mya. These data do not support the current subfamilial groupings of the Bathyergidae. The eastern African naked mole-rat, Heterocephalus glaber, is the most basal representative of the family, with the silvery mole-rat, Heliophobius, being the next most basal. South African forms [dune, common, and cape mole-rats (Bathyergus, Cryptomys, and Georychus, respectively)] group together. The independent origin of the common mole-rat, relative to the naked mole-rat, suggests that complex social systems evolved in parallel along different bathyergid lineages. The 12S rRNA gene is not evolving at a higher rate within the rodent lineages, relative to that seen for artiodactyls and primates. Bathyergid rodents appear to fall at an extreme end of the spectrum of mammalian variation, with respect to both transition/transversion ratios and divergence, showing much lower transition/transversion ratios than those previously reported for intrafamilial comparisons.     

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.     

Sporadic amplification of ID elements in rodents

ID sequences are members of a short interspersed element (SINE) repetitive DNA family within the rodent genome. The copy number of individual ID elements varies by up to three orders of magnitude between species. This amplification has been highly sporadic in the order Rodentia and does not follow any phylogenetic trend. Using library screening and dot-blot analysis, we estimate there are 25,000 copies of ID elements in the deer mouse, 1,500 copies in the gerbil (both cricetid rodents), and 60,000 copies of either ID or ID-like elements in a sciurid rodent (squirrel). By dot-blot analysis, we estimate there are 150,000, 4,000, 1,000, and 200 copies of ID elements in the rat, mouse, hamster, and guinea pig, respectively (which is consistent with previous reports) and 200 copies in the hystricognath rodent, nutria. Therefore, a rapid amplification took place not only after the divergence of rat and mouse but also following the deer mouse (Peromyscus) and hamster split, with no evidence of increased amplifications in hystricognath rodents. No notable variations of sequences from the BC1 genes of several myomorphic rodents were observed that would possibly explain the varied levels of ID amplification. We did observe subgenera and species-group-specific variation in the ID core sequence of the BC1 gene within the genus Peromyscus. Sequence analysis of cloned ID elements in Peromyscus show most ID elements in this genus arose prior to Peromyscus subgenus divergence. Correspondence of the consensus sequence of individual ID elements in gerbil and deer mouse further confirms BC1 as a master gene in ID amplification. Several possible mechanisms responsible for the quantitative variations are explored.The nucleotide sequences reported in this paper have been submitted to the GenBank/EMBL Data Bank with accession numbers: U33850, U33851, U33852 (BC1 sequences); and U33853, U33854, U33855, U33856, U33857, U33858, U33859, U33860, U33861, U33862, U33863, U33864, U33865, U33866, U33867 (ID sequences) Correspondence to: D.H. Kass     

Sequence and evolution of guinea pig preproinsulin DNA

Guinea pig insulin exhibits an unusually high degree of divergence from the conserved insulins of other mammals. cDNA clones encoding guinea pig preproinsulin were isolated, and their nucleic acid sequences were determined. Comparisons of the nucleic acid sequence and its predicted amino acid sequence with sequences encoding insulins of other species revealed that the gene encoding guinea pig preproinsulin evolved from the same ancestral mammalian gene as other known mammalian insulin genes.     

Human 18 S ribosomal RNA sequence inferred from DNA sequence. Variations in 18 S sequences and secondary modification patterns between vertebrates.

We have determined the DNA sequences encoding 18 S ribosomal RNA in man and in the frog, Xenopus borealis. We have also corrected the Xenopus laevis 18 S sequence: an A residue follows G-684 in the sequence. These and other available data provide a number of representative examples of variation in primary structure and secondary modification of 18 S ribosomal RNA between different groups of vertebrates. First, Xenopus laevis and Xenopus borealis 18 S ribosomal genes differ from each other by only two base substitutions, and we have found no evidence of intraspecies heterogeneity within the 18 S ribosomal DNA of Xenopus (in contrast to the Xenopus transcribed spacers). Second, the human 18 S sequence differs from that of Xenopus by approx. 6.5%. About 4% of the differences are single base changes; the remainder comprise insertions in the human sequence and other changes affecting several nucleotides. Most of these more extensive changes are clustered in a relatively short region between nucleotides 190 and 280 in the human sequence. Third, the human 18 S sequence differs from non-primate mammalian sequences by only about 1%. Fourth, nearly all of the 47 methyl groups in mammalian 18 S ribosomal RNA can be located in the sequence. The methyl group distribution corresponds closely to that in Xenopus, but there are several extra methyl groups in mammalian 18 S ribosomal RNA. Finally, minor revisions are made to the estimated numbers of pseudouridines in human and Xenopus 18 S ribosomal RNA.     

Monophyly or polyphyly of the order Rodentia: Possible conflict between morphological and molecular interpretations

Recent analyses of amino acid sequence data from selected proteins inCavia, Rattus, Homo, Bos, Sus, and a few additional rodents and other eutherians suggested that Rodentia is not a monophyletic taxon and thatCavia and other hystricognaths may have branched off earlier than the separation between Muroidea and Primates during mammalian evolution. Because this hypothesis of polyphyly is contrary to the otherwise unanimous recognition of rodent monophyly, we have reevaluated the morphological and developmental evidence from the cranium, dentition, postcranial skeleton, and fetal membranes for the taxa Hystricognathi, Muroidea, other Rodentia, Primates, Artiodactyla, and Lagomorpha, as well as for the eutherian morphotype. Our character analyses provide strong corroboration for the traditional hypothesis of rodent monophyly and lend additional support to the suggestion that Lagomorpha is the sister taxon of Rodentia. Our survey of published molecular data furnishes little or no support for the proposed hypothesis of rodent polyphyly. We conclude that this hypothesis is the result of poor sampling of sequence data from rodents and other eutherians, rather than any inherent difficulties in the use of molecular evidence for the assessment of mammalian evolution. The available molecular data suggest thatCavia differs considerably from other hystricognaths in many proteins, but the reasons for this remain to be investigated.     

Higher-level systematics of rodents (Mammalia,Rodentia): Evidence from the mitochondrial 12S rRNA gene

Phylogenetic relationships among major rodent superfamilies traditionally have been difficult to establish because of the apparent high level of convergence and parallelism seen among morphological characters and/or rapid differentiation of rodent groups in the Paleocene/Eocene. Nucleotide sequence data from the mitochondrial 12S rRNA gene were used to clarify phylogenetic relationships among the major groups of rodents as defined by Brandt (1855) and Tullberg (1899). Based on the approximately 800 bp analyzed for the 12S rRNA gene in 59 mammalian species, including 25 of the 32 extant rodent families, the major rodent groups that could be defined as monophyletic clades were the Hystricognathi, the Muroidea, and the Geomyoidea. In addition, support for superfamilial sister-group relationships was found for Aplodontoidea with Sciuroidea and Dipodoidea with Muroidea.     

Phylogeny of Tetraoninae and other galliform birds using mitochondrial 12S and ND2 genes

The avian subfamily Tetraoninae (grouse and ptarmigan) is a Holarctic group in the order Galliformes distinguished by morphological adaptations to cold environments and behavioral traits associated with elaborate courtship. Here we investigate the relationships of 17 tetraonines and 12 other galliform species using mitochondrial 12S and ND2 sequence data. We found support for the recent phylogenetic classification that separates the genus Dendragapus into two genera, Falcipennis and Dendragapus. In addition, we found support for a tetraonine clade in which the first divergence is between Bonasa umbellus and all others, followed by divergence between a Bonasa bonasia/Bonasa sewerzowi clade and the remaining tetraonines. Falcipennis canadensis is sister to a clade with four Tetrao species, and the genus Centrocercus is sister to a Dendragapus obscurus/Tympanuchus clade. Our data indicate a basal position for Cracidae and Megapodiidae among the five recognized galliform families. We also found strong support for the monophyly of Phasianidae, although the relative positions of Numididae and Odontiphoridae remains unresolved. We use a maximum likelihood approach to infer ages of 37mya for divergence of Numididae and Phasianidae and 28mya for the divergence of Tetraoninae and Meleagris gallopavo. These estimates must be viewed as tentative as they depend on tests of rates of molecular evolution and accurate fossil dates.     

The external promoter in the guinea pig 5S rRNA gene is different from the rodent promoter

The guinea pig has about 100 copies of the 5S rRNA gene per haploid genome and they are present in 2.1 kb tandem repeats. Three bona fide 5S rRNA genes and four pseudo genes were sequenced. The conserved external promoter (D box) found in rodents and primates is only partially present in the guinea pig. The "D box like" sequence in guinea pig only has eight of the 12 nucleotides in the conserved D box. The results are in accordance with investigations showing that the guinea pig is not a rodent. Conserved sequences in the non-transcribed spacer can therefore be useful in phylogenetic studies.     

Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly

The complete mitochondrial genome of the fat dormouse, Glis glis, has been sequenced (16,602 bp). A total of 23 complete mitochondrial mammalian genomes have been taken into account for phylogenetic reconstruction. Phylogenetic analyses were performed with parsimony, distance (stationary Markov model), and maximum-likelihood methods. In all cases, data strongly support the paraphyly of rodents, with dormouse and guinea pig in a different clade from rat and mouse, reaching bootstrap values of 95%. Rodent monophyly and the existence of Glires (Rodentia and Lagomorpha) are weakly supported, with maximum bootstrap values of 11% and 8.6%, respectively. This result agrees with the analyses of isochore patterns in the nuclear genome and the B2 and B2-like retroposons, which show a close relationship between dormice and guinea pigs rather than between dormice and rats and mice.      

A phylogenetic reconsideration of suctorian ciliates (Protista,Ciliophora, Phyllopharyngea) based on small subunit rRNA gene sequences

Compared with other ciliated protozoa, molecular studies of phylogenetic relationships within the subclass Suctoria are rare. In this work, phylogenetic analyses focusing on this group were performed based on all data available. In addition, the small subunit ribosomal RNA (SSU rRNA) genes of three suctorian ciliates (Acineta compressa, Acineta tuberosa and Paracineta limbata) were newly sequenced. Furthermore, the putative secondary structures of the variable region 2 of the SSU rRNA gene were predicted and compared within the Suctoria. Our results show that (i) there is support for the monophyly of the subclass Suctoria, which is a sister clade to the cyrtophorids; (ii) based on combined morphologic and molecular features, we propose the following evolutionary routine within the Suctoria: Exogenina – Evaginogenina – Endogenina; (iii) the similarities of the secondary structures of the V2 region and the SSU rRNA gene sequences within the subclass Suctoria are consistent with the branching of the phylogenetic lineages.     

Phylogenetic relationships among the Malagasy lemuriforms (Primates: Strepsirrhini) as indicated by mitochondrial sequence data from the 12S rRNA gene

Numerous phylogenetic hypotheses have been advanced for the Malagasy lemuriform radiation, drawing on data from morphology, physiology, behaviour and molecular genetics. Almost all possible relationships have been proposed and most nodes have been contested. We present a phylogenetic analysis, using several analytical methods, of a partial sequence from the 12s rRNA mitochondrial gene. This gene codes for the small ribosomal subunit, and functional constraints require that the secondary structure of the molecule is strongly conserved, which in inturn exerts constraints on the primary sequence structure. Although previous studies have suggested a very wide range of phylogenetic applicability for this molecule, our results indicate that it is most useful in strepsirrhine primates for estimating relationships among genera within families and among relatively recently diverged families (mean sequence divergence about 11%). Relationships among families separated by larger genetic distances (>12% divergence; e.g. Cheirogaleidae, Daubentoniidae, Megaladapidae) are difficult to resolve consistently. Our data show strong support for an Indridae-Lemuridae sister group and for monophyly of the Lemuridae with Varecia as the sister to all other lemurids. They also support, albeit less strongly, sister group relationships between Lemur and Hapalemur within the Lemuridae and between PmpLthecus and Avahi in the Indridae.     

The current status of the small subunit rRNA phylogeny of the coccidia (Sporozoa)

There is no current comprehensive assessment of the molecular phylogeny of the coccidia, as all recently published papers either deal with subsets of the taxa or sequence data, or provide non-robust analyses. Here, we present a comprehensive and consistent phylogenetic analysis of the available data for the small-subunit ribosomal RNA gene sequence, including a number of taxa not previously studied, based on a Bayesian tree-building analysis and the covariotide model of evolution. The assumptions of the analysis have been rigorously tested, and the benefits and limitations highlighted. Our results provide support for a number of prior conclusions, including the monophyly of the families Sarcocystidae (cyst-forming coccidia) and Eimeriidae (oocyst-forming coccidia), but with bird-host Isospora species in the Eimeriidae and mammal-host species in the Sarcocystidae. However, it is clear that a number of previously reported relationships are dependent on the evolutionary model chosen, such as the placements of Goussia janae, Lankesterella minimia and Caryospora bigenetica. Our results also confirm the monophyly of the subfamilies Toxoplasmatinae and Sarcocystinae, but only some of the previously reported groups within these subfamilies are supported by our analysis. Similarly, only some of the previously reported groups within the Eimeriidae are supported by our analysis, and the genus Eimeria is clearly paraphyletic. There are unambiguous patterns of host-parasite relationship within the coccidia, as most of the well-supported groups have a consistent and restricted range of hosts, with the exception of the Toxoplasmatinae. Furthermore, the previously reported groups for which we found no support all have a diverse range of unrelated hosts, confirming that these are unlikely to be natural groups. The most interesting unaddressed questions may relate to Isospora, which has the fewest available sequences and host-parasite relationships apparently not as straightforward as elsewhere within the suborder.     

Phylogenetic place of guinea pigs: no support of the rodent-polyphyly hypothesis from maximum-likelihood analyses of multiple protein sequences

Graur et al.'s (1991) hypothesis that the guinea pig-like rodents have an evolutionary origin within mammals that is separate from that of other rodents (the rodent-polyphyly hypothesis) was reexamined by the maximum-likelihood method for protein phylogeny, as well as by the maximum-parsimony and neighbor-joining methods. The overall evidence does not support Graur et al.'s hypothesis, which radically contradicts the traditional view of rodent monophyly. This work demonstrates that we must be careful in choosing a proper method for phylogenetic inference and that an argument based on a small data set (with respect to the length of the sequence and especially the number of species) may be unstable.      

The phylogenetic status of arthropods, as inferred from 18S rRNA sequences

Partial 18S rRNA sequences of five chelicerate arthropods plus a crustacean, myriapod, insect, chordate, echinoderm, annelid, and platyhelminth were compared. The sequence data were used to infer phylogeny by using a maximum-parsimony method, an evolutionary-distance method, and the evolutionary-parsimony method. The phylogenetic inferences generated by maximum-parsimony and distance methods support both monophyly of the Arthropoda and monophyly of the Chelicerata within the Arthropoda. These results are congruent with phylogenies based on rigorous cladistic analyses of morphological characters. Results support the inclusion of the Arthropoda within a spiralian or protostome coelomate clade that is the sister group of a deuterostome clade, refuting the hypothesis that the arthropods represent the "primitive" sister group of a protostome coelomate clade. Bootstrap analyses and consideration of all trees within 1% of the length of the most parsimonious tree suggest that relationships between the nonchelicerate arthropods and relationships within the chelicerate clade cannot be reliably inferred with the partial 18S rRNA sequence data. With the evolutionary-parsimony method, support for monophyly of the Arthropoda is found in the majority of the combinations analyzed if the coelomates are used as "outgroups." Monophyly of the Chelicerata is supported in most combinations assessed. Our analyses also indicate that the evolutionary-parsimony method, like distance and parsimony, may be biased by taxa with long branches. We suggest that a previous study's inference of the Arthropoda as paraphyletic may be the result of (a) having two few arthropod taxa available for analysis and (b) including long-branched taxa.      

Monophyly of Kelloggia Torrey ex Benth. (Rubiaceae) and evolution of its intercontinental disjunction between western North America and eastern Asia

Kelloggia Torrey ex Bentham (Rubiaceae) consists of two species disjunctly distributed in western North America (K. galioides Torrey) and the western part of eastern Asia (K. chinensis Franch.). The two species exhibit a high level of morphological divergence. To test its monophyly and to infer its biogeographic history, we estimated the phylogeny of Kelloggia and its relatives from sequences of three chloroplast DNA regions (rbcL gene, atpB-rbcL spacer, and rps16 intron). The monophyly of Kelloggia was strongly supported, and it forms a sister relationship with the tribe Rubieae. The divergence time between the two disjunct species of Kelloggia was estimated to be 5.42 ± 2.32 million years ago (mya) using the penalized likelihood method based on rbcL sequence data with fossil calibration. Our result does not support the Madrean-Tethyan hypothesis, which assumes an earlier divergence time of 20-25 mya. Ancestral area analysis, as well as dispersal-vicariance (DIVA) analysis, suggests the Asian origin of Kelloggia and the importance of Eurasia in the diversification of its close relatives in the Rubieae-Theligoneae-Paederieae group. The intercontinental disjunction in Kelloggia is suggested to have evolved via long-distance dispersal from Asia into western North America.     

Myriapod monophyly and relationships among myriapod classes based on nearly complete 28S and 18S rDNA sequences

Myriapods play a pivotal position in the arthropod phylogenetic tree. The monophyly of Myriapoda and its internal relationships have been difficult to resolve. This study combined nearly complete 28S and 18S ribosomal RNA gene sequences (3,826 nt in total) to estimate the phylogenetic position of Myriapoda and phylogenetic relationships among four myriapod classes. Our data set consists of six new myriapod sequences and homologous sequences for 18 additional species available in GenBank. Among the six new myriapod sequences, those of the one pauropod and two symphylans are very important additions because they were such difficult taxa to classify in past molecular-phylogenetic studies. Phylogenetic trees were constructed with maximum parsimony, maximum likelihood, and Bayesian analyses. All methods yielded moderate to strong support for the monophyly of Myriapoda. Symphyla grouped strongly with Pauropoda under all analytical conditions. The KH test rejected the traditional view of Dignatha and Progoneata, and the topology obtained here, though not significantly supported, was Diplopoda versus ((Symphyla + Pauropoda) + Chilopoda).