基于孕激素受体基因的类人猿亚目的分子系统学研究

通过对类人猿亚目中部分种类的孕激素受体基因进行分析,重建类人猿亚目的 系统发育关系.扩增并测定了来源于14个属的类人猿亚目物种的孕激素受体编码区序列,并基于这一序列数据,分别采用邻接法、最大简约法和最大似然法重建了系统发育关系.除了阔鼻下目,3种方法构建的系统发生树的拓扑结构类似且各节点支持率高.重建的人猿超科和猴超科内部亲缘关系支持多数人所认可的分类系统.本研究为黑猩猩和人的姐妹群关系提供了证据,提示黑猩猩比大猩猩或其他猿猴更接近人类.阔鼻下目中蜘蛛猴科、卷尾猴科和僧面猴科三者之间的系统发育关系在本研究中未得到很好辨析.     

Phylogenetic relationships between the six superoxide dismutase proteins (FeSOD) of Trichomonas vaginalis and FeSOD6 genetic diversity

The parasitic protozoan Trichomonas vaginalis is known to contain several types of Fe-containing superoxide dismutase proteins (FeSOD). Using three different methods of phylogenetic analysis, maximum parsimony (MP), neighbor joining (NJ), and maximum likelihood (ML) methods, we examined the phylogenetic relationships among the six FeSOD (FeSOD1-FeSOD6) based on their amino acid sequences. All the analyses consistently suggested that the six proteins formed a monophyletic group implying that they probably be originated from an ancestral protein form through repeated duplication events. Although MP tree was totally unresolved, the NJ and ML trees revealed that FeSOD6 placed the most basal position and thus emerged earlier than the other five gene types during the evolution of T. vaginalis. Phylogenetic relationships among the five remaining proteins were (FeSOD2, FeSOD3), (FeSOD4, (FeSOD1, FeSOD5)) although weakly supported in terms of bootstrapping values. In addition to this, we newly designed two PCR primer specifically amplifying full-length FeSOD6 gene and examined its genetic diversity among 12 T. vaginalis isolates from five countries and three continents. They had the same nucleotide sequences except those of three Korean isolates which showed one to three different nucleotides.     

介形纲丽足目和速足目及相关类群18S rDNA分子系统发育的研究

文中分析现生介形类 (Ostracoda) 4目 2 1科 2 9属的 18SrDNA部分序列 ,采用最大似然法 (ML)、邻接法 (NJ)和最大简约法 (MP) ,尝试构建介形类的分子系统树 ;结合介形类的形态特征和化石记录 ,主要对速足目(Podocopida)、丽足目 (Myodocopida)及其超科级分类阶元的系统发生关系进行探讨。 3种分析方法均支持形态学上Podocopida ,Myodocopida和海萤超科 (Cypridinacea)的界定 ;但对Podocopida目土菱介超科 (Bairdiacea)的系统地位提出质疑 ,该类群可能不是单系发生的自然类群。上述分析显示 ,Podocopida,Myodocopida,Platycopida和Halo cypridina组成一个单系群 ;介形类在目、超科、科和属的水平上可能发生过多次辐射分化     

基于COI基因序列的九种真蝽属昆虫的分子系统学研究(半翅目:蝽科)

通过对真蝽属Pentatoma 9种昆虫线粒体COI基因约798bp的序列进行分子进化分析,并以同蝽科宽铗同蝽Acanthosoma labiduroides为外群,采用最大简约法、最大似然法和邻接法构建了分子系统树,来探讨真蝽属的系统发育关系.研究结果支持褐真蝽群P. semiannulata-group的划分,绿角真蝽Pentatoma viridicornuta应划归到褐真蝽群P. Semiannulata-group;红足真蝽群中的角肩真蝽P. angulata与红足真蝽P. rufipes遗传距离较小,它们是否为1个物种值得关注;真蝽属各群之间的系统发育关系以及是否可分为3个属值得进一步研究.     

Reanalysis of Murphy et al.'s data gives various mammalian phylogenies and suggests overcredibility of Bayesian trees

Murphy and colleagues reported that the mammalian phylogeny was resolved by Bayesian phylogenetics. However, the DNA sequences they used had many alignment gaps and undetermined nucleotide sites. We therefore reanalyzed their data by minimizing unshared nucleotide sites and retaining as many species as possible (13 species). In constructing phylogenetic trees, we used the Bayesian, maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ) methods with different substitution models. These trees were constructed by using both protein and DNA sequences. The results showed that the posterior probabilities for Bayesian trees were generally much higher than the bootstrap values for ML, MP, and NJ trees. Two different Bayesian topologies for the same set of species were sometimes supported by high posterior probabilities, implying that two different topologies can be judged to be correct by Bayesian phylogenetics. This suggests that the posterior probability in Bayesian analysis can be excessively high as an indication of statistical confidence and therefore Murphy et al.'s tree, which largely depends on Bayesian posterior probability, may not be correct.     

Molecular phylogenetic analysis of ant subfamily relationship inferred from rDNA sequences

The relationships among ant subfamilies were studied by phylogenetic analysis of rDNA sequences of 15 species from seven subfamilies. PCR primers were designed on the basis of the rDNA sequence of the Australian bulldog ant, Myrmecia croslandi, previously determined. Phylogenetic trees were constructed using sequences of a fragment of 18S rDNA (1.8 kb), a fragment of 28S rDNA (0.7 kb excluding variable regions) and a combination of the 18S and 28S rDNAs, by neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML). rDNA sequences corresponding to the same fragments from three non-ant hymenopteran species (a sawfly, a bee and a wasp) were employed as outgroups. These trees indicated that the ant subfamilies were clustered singly, and, among the seven subfamilies examined, Ponerinae and six other subfamilies are in a sister-groups relationship. The relationship among the six subfamilies, however, was not clarified. The phylogenetic trees constructed in the present study are not in contradiction to the tree from cladistic analysis of morphological data by Baroni Urbani et al. (1992) and the tree from morphological and molecular data (Ward and Brady, 2003), but are inconsistent with the traditional phylogeny. The present results thus raise a question as to the status of some traditionally employed "key" morphological characters. The present results also call for a reexamination of Amblyopone traditionally treated as a member of Ponerinae as belonging to a new subfamily.     

A preliminary molecular and morphological phylogeny of the Antarctic Epimeriidae and Iphimediidae (Crustacea, Amphipoda)

The phylogenetic relationships of 14 species of the Antarctic amphipod families Epimeriidae and Iphimediidae were investigated using 553bp of the gene for the mitochondrial cytochrome oxidase subunit I (COI) and 98 morphological characters. Both families are dominant members of the Antarctic benthic amphipod community. In contrast to previous studies, our molecular and morphological data suggest that the families Epimeriidae and Iphimediidae may not be sister taxa. Our study suggests that Iphimediidae are more closely related to Eusirus (Eusiridae) than to Epimeria (Epimeriidae). Phylogenetic analyses based on maximum parsimony (MP) and maximum likelihood (ML) indicate that the genera Iphimediella and Gnathiphimedia are not monophyletic.     

Efficiencies of different genes and different tree-building methods in recovering a known vertebrate phylogeny

The relative efficiencies of different protein-coding genes of the mitochondrial genome and different tree-building methods in recovering a known vertebrate phylogeny (two whale species, cow, rat, mouse, opossum, chicken, frog, and three bony fish species) was evaluated. The tree-building methods examined were the neighbor joining (NJ), minimum evolution (ME), maximum parsimony (MP), and maximum likelihood (ML), and both nucleotide sequences and deduced amino acid sequences were analyzed. Generally speaking, amino acid sequences were better than nucleotide sequences in obtaining the true tree (topology) or trees close to the true tree. However, when only first and second codon positions data were used, nucleotide sequences produced reasonably good trees. Among the 13 genes examined, Nd5 produced the true tree in all tree-building methods or algorithms for both amino acid and nucleotide sequence data. Genes Cytb and Nd4 also produced the correct tree in most tree-building algorithms when amino acid sequence data were used. By contrast, Co2, Nd1, and Nd41 showed a poor performance. In general, large genes produced better results, and when the entire set of genes was used, all tree-building methods generated the true tree. In each tree-building method, several distance measures or algorithms were used, but all these distance measures or algorithms produced essentially the same results. The ME method, in which many different topologies are examined, was no better than the NJ method, which generates a single final tree. Similarly, an ML method, in which many topologies are examined, was no better than the ML star decomposition algorithm that generates a single final tree. In ML the best substitution model chosen by using the Akaike information criterion produced no better results than simpler substitution models. These results question the utility of the currently used optimization principles in phylogenetic construction. Relatively simple methods such as the NJ and ML star decomposition algorithms seem to produce as good results as those obtained by more sophisticated methods. The efficiencies of the NJ, ME, MP, and ML methods in obtaining the correct tree were nearly the same when amino acid sequence data were used. The most important factor in constructing reliable phylogenetic trees seems to be the number of amino acids or nucleotides used.      

REEXAMINATION OF PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA): AN ANALYSIS OF atpB AND rbcL GENE SEQUENCES

The chloroplast-encoded atp B gene was sequenced from 33 strains representing 28 species of the colonial Volvocales (the Volvocaceae and its relatives) to reexamine phylogenetic relationships as previously deduced by morphological data and rbc L gene sequence data.1128 base pairs in the coding regions of the atp B gene were analyzed by MP, NJ, and ML analyses. Although supported with relatively low bootstrap values (75% and 65% in the NJ and ML analyses, respectively), three anisogamous/oogamous volvocacean genera— Eudorina, Pleodorina, and Volvox, excluding the section Volvox (= Euvolvox, illegitimate name), constituted a large monophyletic group (Eudorina group). Outside the Eudorina group, a robust lineage composed of three species of Volvox sect. Volvox was resolved as in the rbc L gene trees, rejecting the hypothesis of the previous cladistic analysis based on morphological data that the genus Volvox is monophyletic. In addition, the NJ and ML trees suggested that Eudorina is a nonmonophyletic genus as inferred from the morphological data and rbc L gene sequences. Although phylogenetic status of the genus Gonium is ambiguous in the rbc L gene trees and the paraphyly of this genus is resolved in the cladistic analysis based on morphological data, the atp B gene sequence data suggest monophyly of Gonium with relatively low bootstrap values (56–61%) in the NJ and ML trees. On the basis of the combined sequence data (2256 base pairs) from atp B and rbc L genes, Gonium was resolved as a robust monophyletic genus in the NJ and ML trees (with 68–86% bootstrap values), and Eudorina elegans Ehrenberg represented a paraphyletic species positioned most basally within the Eudorina group. However, phylogenetic status and relationships of the families of the colonial Volvocales were still almost ambiguous even in the combined analysis.     

Molecular phylogeny of oribatid mites (Oribatida, Acari): evidence for multiple radiations of parthenogenetic lineages

Nucleotide sequences of the D3 expansion segment and its flanking regions of the 28S rDNA gene were used to evaluate phylogenetic relationships among representative sexual and asexual oribatid mites (Oribatida, Acariformes). The aim of this study was to investigate the hypothesis that oribatid mites consist of species-rich clusters of asexual species that may have radiated while being parthenogenetic. Furthermore, the systematic position of the astigmate mites (Astigmata, Acariformes) which have been hypothesised to represent a paedomorphic lineage within the oribatid mites, is investigated. This is the first phylogenetic tree for oribatid mites s.l. (incl. Astigmata) based on nucleotide sequences. Intraspecific genetic variation in the D3 region was very low, confirming the hypothesis that this region is a good species marker. Results from neighbour joining (NJ) and maximum parsimony (MP) algorithms indicate that several species-rich parthenogenetic groups like Camisiidae, Nanhermanniidae and Malaconothridae are monophyletic, consistent with the hypothesis that some oribatid mite groups diversified despite being parthenogenetic. The MP and maximum likelihood (ML) method indicated that the D3 region is a good tool for elucidating the relationship of oribatid mite species on a small scale(genera, families) but is not reliable for large-scale taxonomy, because branches from the NJ algorithm collapsed in the MP and ML tree. In all trees calculated by different algorithms the Astigmata clustered within the oribatid mites, as proposed earlier.     

基于Cyt b基因探讨羚羊亚科原羚属系统发生关系

原羚属物种在羚羊亚科中的分类地位尚存在很多争议。本文测定了原羚属的黄羊和藏原羚细胞色素b基因全序列(1140bp),并与牛科其它属31个种的同源序列进行比较,对其碱基组成变异情况及核苷酸序列差异进行了分析。基于细胞色素b基因全序列,用简约法(MP)、邻接法(NJ)和似然法(ML)构建了系统进化树。结果表明：黄羊和藏原羚的序列差异为3．78％,颠换数目近乎为0,其突变远未饱和;原羚属内黄羊和藏原羚为不同种,单系发生;原羚属与赛加羚羊属、犬羚属及跳羚属等并系发生,原羚属隶属于羚羊亚科,应为独立属;羚羊亚科组成属间多为并系起源。根据序列差异值2％／百万年的细胞色素6分子钟,推测黄羊和藏原羚分歧时间大约为1～2百万年;原羚属与羚羊亚科其它属分歧时间大约在5．7～8百万年。     

PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA) INFERRED FROM rbcL GENE SEQUENCE DATA

The chloroplast-encoded large subunit of the ribulose-1, 5-bisphosphate carboxylase / oxygenase (rbcL) gene was sequenced from 20 species of the colonial Volvocales (the Volvacaceae, Goniaceae, and Tetrabaenaceae) in order to elucidate phylogenetic relationships within the colonial Volvocales. Eleven hundred twenty-eight base pairs in the coding regions of the (rbcL) gene were analyzed by the neighbor-joining (NJ) method using three kinds of distance estimations, as well as by the maximum parsimony (MP) method. A large group comprising all the anisogamous and oogamous volvocacean species was resolved in the MP tree as well as in the NJ trees based on overall and synonymous substitutions. In all the trees constructed, Basichlamys and Tetrabaena (Tetrabaenaceae) constituted a very robust phylogenetic group. Although not supported by high bootstrap values, the MP tree and the NJ tree based on nonsynonymous substitutions indicated that the Tetrabaenaceae is the sister group to the large group comprising the Volvocaceae and the Goniaceae. In addition, the present analysis strongly suggested that Pandorina and Astrephomene are monophyletic genera whereas Eudorina is nonmonophyletic. These results are essentially consistent with the results of the recent cladistic analyses of morphological data. However, the monophyly of the Volvocaceae previously supported by four morphological synapomorphies is found only in the NJ tree based on nonsynonymous substitutions (with very low bootstrap values). The genus Volvox was clearly resolved as a polyphyletic group with V. rousseletii Pocock separated from other species of Volvox in the rbcL gene comparisons, although this genus represents a monophyletic group in the previous morphological analyses. Furthermore, none of the rbcL gene trees supported the monophyly of the Goniaceae; Astrephomene was placed in various phylogenetic positions .     

Combined nuclear and mitochondrial DNA sequences resolve generic relationships within the Cracidae (Galliformes,Aves)

The Cracidae is one of the most endangered and distinctive bird families in the Neotropics, yet the higher relationships among taxa remain uncertain. The molecular phylogeny of its 11 genera was inferred using 10,678 analyzable sites (5,412 from seven different mitochondrial segments and 5,266 sites from four nuclear genes). We performed combinability tests to check conflicts in phylogenetic signals of separate genes and genomes. Phylogenetic analysis showed that the unrooted tree of ((curassows, horned guan) (guans, chachalacas)) was favored by most data partitions and that different data partitions provided support for different parts of the tree. In particular, the concatenated mitochondrial DNA (mtDNA) genes resolved shallower nodes, whereas the combined nuclear sequences resolved the basal connections among the major clades of curassows, horned guan, chachalacas, and guans. Therefore, we decided that for the Cracidae all data should be combined for phylogenetic analysis. Maximum parsimony (MP), maximum likelihood (ML), and Bayesian analyses of this large data set produced similar trees. The MP tree indicated that guans are the sister group to (horned guan, (curassows, chachalacas)), whereas the ML and Bayesian analysis recovered a tree where the horned guan is a sister clade to curassows, and these two clades had the chachalacas as a sister group. Parametric bootstrapping showed that alternative trees previously proposed for the cracid genera are significantly less likely than our estimate of their relationships. A likelihood ratio test of the hypothesis of a molecular clock for cracid mtDNA sequences using the optimal ML topology did not reject rate constancy of substitutions through time. We estimated cracids to have originated between 64 and 90 million years ago (MYA), with a mean estimate of 76 MYA. Diversification of the genera occurred approximately 41-3 MYA, corresponding with periods of global climate change and other Earth history events that likely promoted divergences of higher level taxa.     

分子系统进化关系分析的一种新方法--贝叶斯法在硬蜱属中的应用

距离矩阵邻接法、最大简约法和最大似然法是重建生物系统关系的3种主要方法。普遍认为最大似然法在原理上优于前二种方法，但其计算复杂费时。由于现行计算机的能力尚达不到其要求而实用性差，特别是在处理大数据集样本(即大于25个分类单元)时，用此方法几乎不可能。新近提出的贝叶斯法(Bayesianmethod)既保留了最大似然法的基本原理，又引进了马尔科夫链的蒙特卡洛方法，并使计算时间大大缩短。本文用贝叶斯法对硬蜱属(Ixodes)19个种的线粒体16S rDNA片段进行了系统进化分析。从总体上看，分析结果与现有的基于形态学的分类体系基本吻合。但与现存的假说相反，莱姆病的主要宿主蓖籽硬蜱复合种组并非单系。通过比较贝叶斯法与其它三种方法的结果，我们认为贝叶斯法是一种系统进化分析的好方法，它既能根据分子进化的现有理论和各种模型用概率重建系统进化关系，又克服了最大似然法计算速度慢、不适用于大数据集样本的缺陷。贝叶斯法根据后验概率直观地表示系统进化关系的分析结果，不需要用自引导法进行检验。可以预料，贝叶斯法将会被广泛地应用到系统进化分析上[动物学报49(3)：380—388，2003]。     

Molecular Phylogeny of the Genus Paramesotriton (Caudata: Salamandridae)

To elucidate the phylogeny of the genus Paramesotriton (Caudata: Salamandridae), we investigated three mitochondrial DNA gene fragments (1207 bp in total) of cytochrome b, ND2, and ND4 for its six recognized species. The phylogenetic relationships within Paramesotriton were reconstructed by maximum parsimony (MP) and maximum likelihood (ML) methods. Phylogenetic trees (MP and ML trees) that were constructed from the combined data set of the three gene fragments indicated that all six species of Paramesotriton formed a monophyletic group, with P. caudopunctatus as basal to the other five species. This result suggests that P. fuzhongensis is a valid species in Paramesotriton.     

Efficiencies of fast algorithms of phylogenetic inference under the criteria of maximum parsimony, minimum evolution, and maximum likelihood when a large number of sequences are used

In phylogenetic inference by maximum-parsimony (MP), minimum-evolution (ME), and maximum-likelihood (ML) methods, it is customary to conduct extensive heuristic searches of MP, ME, and ML trees, examining a large number of different topologies. However, these extensive searches tend to give incorrect tree topologies. Here we show by extensive computer simulation that when the number of nucleotide sequences (m) is large and the number of nucleotides used (n) is relatively small, the simple MP or ML tree search algorithms such as the stepwise addition (SA) plus nearest neighbor interchange (NNI) search and the SA plus subtree pruning regrafting (SPR) search are as efficient as the extensive search algorithms such as the SA plus tree bisection-reconnection (TBR) search in inferring the true tree. In the case of ME methods, the simple neighbor-joining (NJ) algorithm is as efficient as or more efficient than the extensive NJ+TBR search. We show that when ME methods are used, the simple p distance generally gives better results in phylogenetic inference than more complicated distance measures such as the Hasegawa-Kishino-Yano (HKY) distance, even when nucleotide substitution follows the HKY model. When ML methods are used, the simple Jukes-Cantor (JC) model of phylogenetic inference generally shows a better performance than the HKY model even if the likelihood value for the HKY model is much higher than that for the JC model. This indicates that at least in the present case, selecting of a substitution model by using the likelihood ratio test or the AIC index is not appropriate. When n is small relative to m and the extent of sequence divergence is high, the NJ method with p distance often shows a better performance than ML methods with the JC model. However, when the level of sequence divergence is low, this is not the case.     

Evaluation of the monophyly of Fomitopsis using parsimony and MCMC methods

To evaluate the monophyly of Fomitopsis and elucidate phylogenetic relationships of its members, partial nuclear large subunit (partial 28S) ribosomal RNA genes were sequenced from 10 species of Fomitopsis and 15 related species. Phylogenetic analyses indicated that Fomitopsis was phylogenetically heterogeneous and its members were divided into three subgroups. The constrained tree excluding F. palustris (the type species of Pilatoporus) from Fomitopsis core group was rejected, thus rejecting the taxonomic concept to segregate Pilatoporus from Fomitopsis. The monophyly of taxa belonging to F. rosea complex was rejected, thus rejecting the complex definition based on morphological similarities. The exclusion of Piptoporus betulinus (the type species of Piptoporus) from Fomitopsis core group was rejected and Piptoporus proved to be heterogeneous in both best MP and MAP trees. The monophyly of F. officinalis with Fomitopsis core group also was rejected. Fomitopsis officinalis was closely related to Antrodia xantha and formed an independent lineage from Fomitopsis core group at the basal position of brown rotting fungi comprising Antrodia, Daedalea, Fomitopsis, Piptoporus and Postia. The MAP tree topologyobtained from MCMC computation of Bayesian inference was similar to the one of the best MP tree based on the parsimony analysis but showed a higher likelihood score in the Kishino-Hasegawa test and reflected better evolutionary patterns for the phylogeny of Fomitopsis.     

Phylogenetic relationships of land plants using mitochondrial small-subunit rDNA sequences

Phylogenetic relationships among embryophytes (tracheophytes, mosses, liverworts, and hornworts) were examined using 21 newly generated mitochondrial small-subunit (19S) rDNA sequences. The "core" 19S rDNA contained more phylogenetically informative sites and lower homoplasy than either nuclear 18S or plastid 16S rDNA. Results of phylogenetic analyses using parsimony (MP) and likelihood (ML) were generally congruent. Using MP, two trees were obtained that resolved either liverworts or hornworts as the basal land plant clade. The optimal ML tree showed hornworts as basal. That topology was not statistically different from the two MP trees, thus both appear to be equally viable evolutionary hypotheses. High bootstrap support was obtained for the majority of higher level embryophyte clades named in a recent morphologically based classification, e.g., Tracheophyta, Euphyllophytina, Lycophytina, and Spermatophytata. Strong support was also obtained for the following monophyletic groups: hornworts, liverworts, mosses, lycopsids, leptosporangiate and eusporangiate ferns, gymnosperms and angiosperms. This molecular analysis supported a sister relationship between Equisetum and leptosporangiate ferns and a monophyletic gymnosperms sister to angiosperms. The topologies of deeper clades were affected by taxon inclusion (particularly hornworts) as demonstrated by jackknife analyses. This study represents the first use of mitochondrial 19S rDNA for phylogenetic purposes and it appears well-suited for examining intermediate to deep evolutionary relationships among embryophytes.     

Phylogeny of the basal angiosperm genus Pseuduvaria (Annonaceae) inferred from five chloroplast DNA regions, with interpretation of morphological character evolution

Phylogenetic relationships within the magnoliid basal angiosperm genus Pseuduvaria (Annonaceae) are investigated using chloroplast DNA sequences from five regions: psbA-trnH spacer, trnL-F, matK, rbcL, and atpB-rbcL spacer. Over 4000 nucleotides from 51 species (of the total 53) were sequenced. The five cpDNA datasets were analyzed separately and in combination using maximum parsimony (MP), maximum likelihood (ML), and Bayesian methods. The phylogenetic trees constructed using all three phylogenetic methods, based on the combined data, strongly support the monophyly of Pseuduvaria following the inclusion of Craibella phuyensis. The trees generated using MP were less well resolved, but relationships are similar to those obtained using the other methods. ML and Bayesian analyses recovered trees with short branch lengths, showing five main clades. This study highlights the evolutionary changes in seven selected morphological characters (floral sex, stamen and carpel numbers, inner petal color, presence of inner petal glands, flowering peduncle length, and monocarp size). Although floral unisexuality is ancestral within the genus, several evolutionary lineages reveal reversal to bisexuality. Other phylogenetic transitions include the evolution of sapromyophily, and fruit-bat frugivory and seed dispersal, thus allowing a wide range of adaptations for species survival.     

New molecular data for tardigrade phylogeny, with the erection of Paramacrobiotus gen. nov.

Up to few years ago, the phylogenies of tardigrade taxa have been investigated using morphological data, but relationships within and between many taxa are still unresolved. Our aim has been to verify those relationships adding molecular analysis to morphological analysis, using nearly complete 18S ribosomal DNA gene sequences (five new) of 19 species, as well as cytochrome oxidase subunit 1 (COI) mitochondrial DNA gene sequences (15 new) from 20 species, from a total of seven families. The 18S rDNA tree was calculated by minimum evolution, maximum parsimony (MP) and maximum likelihood (ML) analyses. DNA sequences coding for COI were translated to amino acid sequences and a tree was also calculated by neighbour-joining, MP and ML analyses. For both trees (18S rDNA and COI) posterior probabilities were calculated by MrBayes. Prominent findings are as follows: the molecular data on Echiniscidae (Heterotardigrada) are in line with the phylogenetic relationships identifiable by morphological analysis. Among Eutardigrada, orders Apochela and Parachela are confirmed as sister groups. Ramazzottius (Hypsibiidae) results more related to Macrobiotidae than to the genera here considered of Hypsibiidae. Macrobiotidae and Macrobiotus result not monophyletic and confirm morphological data on the presence of at least two large groups within Macrobiotus. Using 18S rDNA and COI mtDNA genes, a new phylogenetic line has been identified within Macrobiotus , corresponding to the ' richtersi-areolatus group'. Moreover, cryptic species have been identified within the Macrobiotus ' richtersi group' and within Richtersius . Some evolutionary lines of tardigrades are confirmed, but others suggest taxonomic revision. In particular, the new genus Paramacrobiotus gen. n. has been identified, corresponding to the phylogenetic line represented by the ' richtersi-areolatus group'.