New nuclear markers and exploration of the relationships among Serraniformes (Acanthomorpha, Teleostei): the importance of working at multiple scales.

We explore the relationships within Serraniformes (Li et al., 2009) using a dense taxon sampling and seven nuclear markers. Six had already used been for teleost phylogeny (IRBP, MC1R, MLL4, Pkd1, Rhodopsin, and RNF213) at other scales, and one (MLL2) is new. The results corroborate the composition of Serraniformes described in previous publications (some Gasterosteiformes, Perciformes and Scorpaeniformes). Within the clade, Notothenioidei and Zoarcoidei are each monophyletic. Cottoidei was not monophyletic due to placement of the genus Ebinania (Psychrolutidae). Our independent data confirm the sister-group relationship of Percophidae and Notothenioidei as well as the division of Platycephaloidei in four different groups (Bembridae, Platycephalidae, Hoplichthyidae and Peristediidae with Triglidae). Within Cottoidei, Liparidae and Cyclopteridae formed a clade associated with Cottidae, the genus Cottunculus (Psychrolutidae), and Agonidae. Serranidae and Scorpaenidae are not monophyletic, with the Serranidae divided in two clades (Serraninae and Epinephelinae/Anthiinae) and Scorpaenidae including Caracanthidae and the genus Ebinania (Psychrolutidae). We discuss some morphological characters supporting clades within the Scorpaenidae.

Tracking the elusive monophyly of nototheniid fishes (Teleostei) with multiple mitochondrial and nuclear markers.

Since the first molecular study of the suborder Notothenioidei in 1994, many phylogenetic studies have been published. Among these, those with a sufficient number of taxa have all suggested that the Nototheniidae, as currently defined, is monophyletic only with the inclusion of the Channichthyidae, Artedidraconidae, Bathydraconidae and Harpagiferidae. This is corroborated by more recent studies including more taxa, but in these studies either the number of nuclear markers or the number of taxa included remained low. We obtained sequences for a large sampling covering most nototheniid genera for five markers described previously for other samplings (COI, Rhodopsin retrogene, Pkd1, HECW2, and SSRP1) and one nuclear marker never used before in phylogenetic inference (PPM1d). The topology for the combined analysis of the nuclear coding genes, as well as the topology for SSRP1 (non-coding) and the combined analysis for all markers all support the paraphyly of Nototheniidae, the genus Notothenia (including Paranotothenia) is the sister group of the clade Channichthyidae, Artedidraconidae, Bathydraconidae and Harpagiferidae, and genus Gobionotothen is a sister group to both. As in previous studies, Trematomus, Lepidonotothen and Patagonotothen form a clade that also includes Indonotothenia cyanobrancha. The position of Pleuragramma antarctica, Dissostichus species and Aethotaxis mitopteryx remains unstable and dependant on markers and analyses. We therefore propose the inclusion of the four families of the High Antarctic clade in the Nototheniidae, and their transformation into subfamilies. We transfer Paranotothenia magellanica to the genus Notothenia, as Notothenia magellanica.

Phylogenetic footprints of an Antarctic radiation: the Trematominae (Notothenioidei, Teleostei).

The teleost suborder Notothenioidei is restricted to the Southern Ocean and has been described as a species flock spanning the whole of it. Within the suborder, the subfamily Trematominae is important for coastal Antarctic ecosystems. The eleven Trematomus species occupy a large range of ecological niches. The genus is monophyletic if the genus Pagothenia (two additional species) and Cryothenia amphitreta, also nested within it, are included. Although the Trematominae have received much interest, the relationships among these fourteen species are still unclear. Several recent studies have tried to resolve these interrelationships; however no complete and clear picture has emerged, probably because of the use of a low number of insufficiently variable markers. The only common results places T. scotti as the sister-group of the rest of the subfamily and T. loennbergi close to T. lepidorhinus. We use here more variable markers. Four nuclear markers, two of which are new, and a mitochondrial marker for the biggest trematomine sampling ever gathered (14 species, 78 specimens). We found that several nuclear haplotypes are shared by several species (mostly in very closely related species). The haplotype patterns coupled with the cytogenetics of the subfamily suggest that a phenomenon of incomplete lineage sorting (ILS) is likely to be at play. Using a calibration linked to fossil evidence, we evaluate the relative ages of each clade within the Trematominae to assess the proximity of the speciation events to one another. The main trematomine diversification was recent and sudden.

Ordering events of biochemical evolution.

Metabolic pathways exhibit structures resulting from an evolutionary process. Pathways have been inherited through time with modification, from the earliest periods of life. It is possible to compare the structure of pathways as done in comparative anatomy, i.e. for inferring ancestral pathways or parts of it (ancestral enzymatic functions), using standard phylogenetic reconstruction. Thus a phylogenetic tree of pathways provides a relative ordering of the rise of enzymatic functions. It even becomes possible to order the birth of each complete pathway in time. This particular "DNA-free" conceptual approach to evolutionary biochemistry is reviewed, gathering all the justifications given for it. Then, the method of assigning a given pathway to a time span of biochemical development is revisited. The previous method used an implicit "clock" of metabolic development that is difficult to justify. We develop a new clock-free approach, using functional biochemical arguments. Results of the two methods are not significantly different; our method is just more precise. This suggests that the clock assumed in the first method does not provoke any important artefact in describing the development of biochemical evolution. It is just unnecessary to postulate it. As a result, most of the amino acid metabolic pathways develop forwards, confirming former models of amino acid catabolism evolution, but not those for amino acid anabolism. The order of appearance of sectors of universal cellular metabolism is: (1) amino acid catabolism, (2) amino acid anabolism and closure of the urea cycle, (3) glycolysis and glycogenesis, (4) closure of the pentose-phosphate cycle, (5) closure of the Krebs cycle and fatty acids metabolism, (6) closure of the Calvin cycle.

Structural evolution of Otx genes in craniates.

Using a degenerate PCR approach, we performed an exhaustive search of Otx genes in the reedfish Erpetoichthys calabaricus, the dogfish Scyliorhinus canicula, and the hagfish Myxine glutinosa. Three novel Otx genes were identified in each of these species, and their deduced protein sequences were determined over a large C-terminal fragment located immediately downstream of the homeodomain. Like their lamprey and osteichthyan counterparts, these nine genes display a tandem duplication of a 20--25-residue C-terminal domain, which appears to be a hallmark of all craniate Otx genes identified thus far, including the highly divergent Crx gene. Phylogenetic analyses show that, together with their osteichthyan counterparts, the dogfish and reedfish genes can be classified into three gnathostome orthology classes. Two of the three genes identified in each of these species belong to the Otx1 and Otx2 orthology classes previously characterized in osteichthyans. The third one unambiguously clusters with the Otx5/Otx5b genes recently characterized in Xenopus laevis, thus defining a novel orthology class. Our results also strongly suggest that the highly divergent Crx genes identified in humans, rodents, and oxen are the mammalian representatives of this third class. The hagfish genes display no clear relationships to the three gnathostome orthology classes, but one of them appears to be closely related to the LjOtxA gene, previously identified in Lampetra japonica. Taken together, these data support the hypothesis that the Otx multigene families characterized in craniates all derive from duplications of a single ancestral gene which occurred after the splitting of cephalochordates but prior to the gnathostome radiation. Using site-by-site sequence comparisons of the gnathostome Otx proteins, we also identified structural constraints selectively acting on each of the three gnathostome orthology classes. This suggests that specialized functions for each of these orthology classes were fixed in the gnathostome lineage prior to the splitting between osteichthyans and chondrichthyans.

A combined analysis of the cystic fibrosis transmembrane conductance regulator: implications for structure and disease models.

Over the past decade, nearly 1,000 variants have been identified in the cystic fibrosis transmembrane conductance regulator (CFTR) gene in classic and atypical cystic fibrosis (CF) patients worldwide, and an enormous wealth of information concerning the structure and function of the protein has also been accumulated. These data, if evaluated together in a sequence comparison of all currently available CFTR homologs, are likely to refine the global structure-function relationship of the protein, which will, in turn, facilitate interpretation of the identified mutations in the gene. Based on such a combined analysis, we had recently defined a "functional R domain" of the CFTR protein. First, presenting two full-length cDNA sequences (termed sCFTR-I and sCFTR-II) from the Atlantic salmon (Salmo salar) and an additional partial coding sequence from the eastern gray kangaroo (Macropus giganteus), this study went further to refine the boundaries of the two nucleotide-binding domains (NBDs) and the COOH-terminal tail (C-tail), wherein NBD1 was defined as going from P439 to G646, NBD2 as going from A1225 to E1417, and the C-tail as going from E1418 to L1480. This approach also provided further insights into the differential roles of the two halves of CFTR and highlighted several well-conserved motifs that may be involved in inter- or intramolecular interactions. Moreover, a serious concern that a certain fraction of missense mutations identified in the CFTR gene may not have functional consequences was raised. Finally, phylogenetic analysis of all the full-length CFTR amino acid sequences and an extended set of exon 13--coding nucleotide sequences reinforced the idea that the rabbit may represent a better CF model than the mouse and strengthened the assertion that a long-branch attraction artifact separates the murine rodents from the rabbit and the guinea pig, the other Glires.

New clades of euthyneuran gastropods (Mollusca) from 28S rRNA sequences.

Recent morphological and molecular results on phylogeny of euthyneuran gastropods, which include opisthobranchs and pulmonates, have greatly diminished previous supposed resolution of their phylogenetic relationships. In addition to recent morphological results, sequences of the D1 and D2 domains of the 28S rRNA are here analyzed by parsimony for 31 euthyneuran species. The molecular and previous morphological data sets were not congruent according to an ILD test, and morphological and molecular data could not be analyzed simultaneously. Consequently Bremer's Combinable Component Consensus was used to obtain a new tree, with the following supported molecular results: monophyly of a new clade of opisthobranchs including actively swimming Euthyneura, i.e., pelagic Gymnosomata and Thecosomata plus benthic Anaspidea; first molecular confirmation of monophylies of Hygrophila, including Chilina, Acteonoidea, and Sacoglossa, which include both shell-bearing species and slugs; and new confirmation of the monophyly of Stylommatophora. Morphological characters which support the new clades obtained here are discussed.

Partial combination applied to phylogeny of European cyprinids using the mitochondrial control region.

Previous molecular phylogenies of European cyprinids led to some solid facts and some uncertainties. This study is based on a stretch of more than 1 kb in the mitochondrial control region newly sequenced for 35 European cyprinids and on previous cytochrome b and 16S rDNA data. The trees based on the control region are more accurate and robust than those obtained from the two other genes. Character incongruence among the three genes was tested using the incongruence length difference (ILD) test. Iterative removals of individual sequences followed by new ILD tests identified two sequences responsible for statistically significant incongruence. A partial combination was conducted, that is, a combination of the three data sets, removing the two sequences previously identified. The phylogenetic analysis of this partial combination gives a more robust and resolved picture of subfamilial interrelationships. The Rasborinae are the sister group of all other cyprinids. The monophyletic Cyprininae emerges next. Tinca tinca first and then Rhodeus are the sister groups of all the remaining nonrasborine and noncyprinine species. Gobio is the sister group of the Leuciscinae, in which the Phoxinini are the sister group of the Leuciscini. Within the Leuciscini, the genus Leuciscus and the subfamily Alburninae are both paraphyletic. The Rasborinae are the most basal cyprinid subfamily and the Tincinae are not the sister group of the Cyprininae. These two results challenge only two anatomical characters, which need to be reinterpreted or considered as homoplastic in cyprinid evolution: the modification of the first pleural rib and its parapophysis and the bony composition of the interorbital septum.

Microsatellite typing as a new tool for identification of Saccharomyces cerevisiae strains.

Since Saccharomyces cerevisiae appears to be an emerging pathogen, there is a need for a valuable molecular marker able to distinguish among strains. In this work, we investigated the potential value of microsatellite length polymorphism with a panel of 91 isolates, including 41 clinical isolates, 14 laboratory strains, and 28 strains with industrial relevance. Testing seven polymorphic regions (five trinucleotide repeats and two dinucleotide repeats) in a subgroup of 58 unrelated strains identified a total of 69 alleles (6 to 13 per locus) giving 52 different patterns with a discriminatory power of 99.03%. We found a cluster of clinical isolates sharing their genotype with a bakery strain, suggesting a digestive colonization following ingestion of this strain with diet. With the exception of this cluster of isolates and isolates collected from the same patient or from patients treated with Saccharomyces boulardii, all clinical isolates gave different and unique patterns. The genotypes are stable, and the method is reproducible. The possibility to make the method portable is of great interest for further studies using this technique. This work shows the possibility to readily identify S. boulardii (a strain increasingly isolated from invasive infections) using a unique and specific microsatellite allele.

Evolutionary implications of the frequent horizontal transfer of mismatch repair genes.

Mutation and subsequent recombination events create genetic diversity, which is subjected to natural selection. Bacterial mismatch repair (MMR) deficient mutants, exhibiting high mutation and homologous recombination rates, are frequently found in natural populations. Therefore, we have explored the possibility that MMR deficiency emerging in nature has left some "imprint" in the sequence of bacterial genomes. Comparative molecular phylogeny of MMR genes from natural Escherichia coli isolates shows that, compared to housekeeping genes, individual functional MMR genes exhibit high sequence mosaicism derived from diverse phylogenetic lineages. This apparent horizontal gene transfer correlates with hyperrecombination phenotype of MMR-deficient mutators. The sequence mosaicism of MMR genes may be a hallmark of a mechanism of adaptive evolution that involves modulation of mutation and recombination rates by recurrent losses and reacquisitions of MMR gene functions.

Phylogenetic relationships of mormyrid electric fishes (Mormyridae; Teleostei) inferred from cytochrome b sequences.

The Mormyridae are African osteoglossomorph freshwater fishes of great interest because of their electric organs. They have become an important model in studies of electrophysiology and behavior but their phylogenetic relationships are poorly known. Phylogenetic relationships among mormyrids were determined by comparing cytochrome b sequences (588 bp) of 27 species belonging to 15 genera. Results showed that the Petrocephalus species (subfamily Petrocephalinae) are the sister-group of all other mormyrids (subfamily Mormyrinae). The monophyly of the Mormyrinae was supported, as well as three original intra-Mormyrinae clades. Three genera, Marcusenius, Pollimyrus, and Brienomyrus, were found to be polyphyletic with high support. Some of these polyphylies are tentatively explained. The results confirmed that the lateral ethmoid bone was lost several times within the Mormyrinae. These findings emphasize the necessity of systematic studies and taxonomic revision of the Mormyridae. The tree obtained from the mitochondrial data showed a single rise of each electrocyte type except for electrocyte with penetrating stalk ("Pa"). Constraining the single occurrence of electrocyte type Pa did not require an excessive number of extra steps (1.86%).

Escherichia coli molecular phylogeny using the incongruence length difference test.

Molecular phylogeny of the species Escherichia coli using the E. coli reference (ECOR) collection strains has been hampered by (1) the absence of rooting in the commonly used phenogram obtained from multilocus enzyme electrophoresis (MLEE) data and (2) the existence of recombination events between strains that scramble phylogenetic trees reconstructed from the nucleotide sequences of genes. We attempted to determine the phylogeny for E. coli based on the ECOR strain data by extracting from GenBank the nucleotide sequences of 11 chromosomal structural and 2 plasmid genes for which the Salmonella enterica homologous gene sequences were available. For each of the 13 DNA data sets studied, incongruence with a nonnucleotide whole-genome data set including MLEE, random amplified polymorphic DNA, and rrn restriction fragment length polymorphism data was measured using the incongruence length difference (ILD) test of Farris et al. As previously reported, the incongruence observed between the gnd and plasmid gene data and the whole-genome data was multiple, indicating numerous horizontal transfer and/or recombination events. In five cases, the incongruence detected by the ILD test was punctual, and the donor group was identified. Congruence was not rejected for the remaining data sets. The strains responsible for incongruences with the whole-genome data set were removed, leading to a "prior-agreement" approach, i.e., the determination of a phylogeny for E. coli based on several genes, excluding (1) the genes with multiple incongruences with the whole genome data, (2) the strains responsible for punctual incongruences, and (3) the genes incongruent with each other. The obtained phylogeny shows that the most basal group of E. coli strains is the B2 group rather than the A group, as generally thought. The D group then emerges as the sister group of the rest. Finally, the A and B1 groups are sister groups. Interestingly, the most primitive taxon within E. coli in terms of branching pattern, i.e., the B2 group, includes highly virulent extraintestinal strains with derived characters (extraintestinal virulence determinants) occurring on its own branch.

[Molecular systematics of Phlebotominae: a pilot study. Paraphyly of the genus Phlebotomus]. / Systématique moléculaire des Phlebotominae: étude pilote. Paraphylie du genre Phlebotomus.

Phylogenetic relationships among Phlebotominae were inferred through a pilot study using parsimony analysis of the D2 domain of ribosomal DNA sequences: 455 pairs of bases were sequenced in nine species of Phlebotomine sandflies which belong to the genera Lutzomyia, Phlebotomus and Sergentomyia. Two taxa are used as outgroups: Psychoda sp. and Nemapalpus flavus which is the sister group of the Phlebotominae. The South American genus Lutzomyia appears to be monophyletic. The Mediterranean species Sergentomyia dentata is its sister group and is not clustered with the Old World genus Phlebotomus. The latter is a paraphyletic genus with an early individualization of the branch including the closely related subgenera Phlebotomus and Paraphlebotomus, and a late individualization of the subgenus Larroussius. These results have some consequences on the biogeography of the leishmaniasis in the Old World.

Mitochondrial phylogeny of the European cyprinids: implications for their systematics, reticulate evolution, and colonization time.

Two different mitochondrial genes, the cytochrome b and the 16S rDNA, support the same European cyprinid molecular phylogeny: the most basal subfamily is the paraphyletic Rasborinae, the Cyprininae are monophyletic, the Tincinae and Gobioninae are close to the Cyprininae or more basal lineages but not close to Leuciscinae or Alburninae, and the Leuciscinae are paraphyletic but can become monophyletic if we include the biphyletic alburninae and exclude the Phoxinini. The relationship of the Acheilognathinae remains obscure. Natural intergeneric and interspecific hybridizations are clearly demonstrated within the Leuciscinae, both from high bootstrap proportions and intermediate morphological features: Chondrostoma toxostoma and Rutilus rutilus, Scardinius erythrophthalmus and R. rutilus, and Leuciscus multicellus and Leuciscus soufia. Finally, the use of the nonsaturated and clockwise 16S mtDNA sequences have been used to infer from nonintrogressive taxa the time of the first European cyprinid cladogeneses. The estimation confirms the hypothesis of Almaça and Banarescu that European cyprinid subfamilies started to diversify 35 mya and confirms the hypothesis of Bianco on the diversification of European leuciscines in the Mediterranean area during the late Messinian (6.5 to 5.3 mya).

Weighting and congruence: a case study based on three mitochondrial genes in pitvipers.

Weighting strategies in a total evidence approach are often conducted in order to remove homoplasy, with the implicit hope to increase congruence between data partitions. Incongruence was measured using the test of Farris et al. (J. S. Farris, M. Källersjö, A. G. Kluge, and C. Bult, 1995, Cladistics 10, 315-319) within and between three mitochondrial genes from pitvipers (Crotalinae) in partitioning each codon position for the coding genes. Incongruence between codon positions of a gene can be worse than incongruence between genes at analogous positions. Third positions of codons behave very differently in terms of incongruence from one gene to another while showing similar patterns in saturation tests. Instead of removing characters in order to discard homoplasy, which is hopeless and does not increase general congruence, we advocate for the removal of those substitutions that are incongruent with the rest. The genus Calloselasma and its sister group the genus Hypnale are the most basal Crotalinae. Asiatic pitvipers are paraphyletic, while American pitvipers are monophyletic.

Low divergence in rDNA ITS sequences among five species of Fucus (Phaeophyceae) suggests a very recent radiation.

Sequences from the two ribosomal DNA internal transcribed spacers (ITS1 and ITS2) were compared among five species of Fucus. Based on the present taxon sampling, parsimony analysis showed that Fucus serratus is the sister-group of the remaining Fucus species when Ascophyllum nodosum was used as an outgroup. The topology of the tree was (Fucus serratus (F. lutarius (F. vesiculosus (F. spiralis + F. ceranoides)))). The extremely low variation observed suggests a very recent radiation of the genus which supports the view widely accepted that the Fucales are among the most evolutionarily advanced of the brown algae. We further note that sequence differences between Fucus and Ascophyllum were 28%: this does not rule out the utility of ITS sequences within the Fucaceae. The very low number of informative positions allows to demonstrate empirically that distance matrix methods group on the basis of symplesiomorphies.

The 'evolutionary signal' of homoplasy in protein-coding gene sequences and its consequences for a priori weighting in phylogeny.

To analyse independently homoplasy for the six possible types of substitution (i.e., A-G, C-T, A-C, A-T, C-G and G-T) at each of the three codon-positions of the cytochrome b gene, two approaches were used: the first is based on the consistency index which measures the amount of homoplasy, and the second is based on the saturation analysis which describes graphically the distribution of homoplasy within the taxonomic sampling. The results obtained from a data set of 32 sequences of Artiodactyla indicate that evolution of the cytochrome b is governed by differential constraints: 1) between the six substitutions-types, 2) between the three codon-positions, and 3) between the two mtDNA strands. Moreover, we find that non-synonymous sites can be more homoplastic than synonymous sites when the possibilities of substitutions are severely restricted because of the functional requirements of hydrophobicity. Most weighting schemes applied to protein-coding genes are elaborated from unjustified assumptions. We propose to weight each substitution-type at each codon-position according to its homoplasy content evaluated either with the consistency index or with an index representing the level of mutational saturation.

Phylogenetic analysis of cystic fibrosis transmembrane conductance regulator gene in mammalian species argues for the development of a rabbit model for cystic fibrosis.

The species-specific pattern of cystic fibrosis transmembrane conductance regulator (CFTR) expression was investigated in order to identify species closely related to man which can be used as potential cystic fibrosis (CF) animal models. To this purpose, the nucleotide sequences of the CFTR promoter region of eight mammalian species representing four different orders (Primates, Artiodactyla, Lagomorpha and Rodentia) were analyzed. Distance matrices and unrooted trees of the CFTR promoter region sequences yielded two deeply separated groups, one including man (Homo sapiens), nonhuman primates (Hylobates lar, Macaca fascicularis, Saimiri sciureus), cow (Bos taurus), and rabbit (Oryctolagus cuniculus) and the other including the rodents (Rattus norvegicus, Mus musculus). Divergences between rodent and nonrodent groups have been observed in putative cis transcriptional regulatory elements and can be involved in the differences of pattern of expression between these two groups. Comparison of the available CFTR cDNA sequences enabled us to root the tree with a noneutherian outgroup and to perform a phylogenetic analysis. This analysis did not detect any base composition bias and supported polyphyletic Glires. Although a long-branch attraction artifact cannot be completely excluded, these findings converge toward the recent statement (Graur, Duret, and Gouy 1996) that Lagomorpha is more closely related to Primates than to Rodentia. In addition, the phenylalanine residue in exon 10 involved in the most common CF mutation in man is conserved in rabbit. These phylogenetic analyses as well as anatomical and developmental data suggest that, once rabbit embryonic stem cells become available, the rabbit will provide a suitable tool for both gene transfer and pharmacological investigations and could lead to a better CF model than the current murine models.

Emergence and scattering of multiple neurofibromatosis (NF1)-related sequences during hominoid evolution suggest a process of pericentromeric interchromosomal transposition.

Type 1 neurofibromatosis (NF1) gene encodes for a member of the GTPase activating protein family and is considered to be a tumor suppressor gene. Its very high rate of de novo mutation in humans led us to study a specific feature of this gene: the presence of numerous NF1-related sequences. According to our results, the human genome contains at least 11 NF1-related sequences, nine of which are scattered near centromeric sequences of seven different chromosomes. These NF1-related sequences, whose extent is quite varied according to loci, are unprocessed copies of the NF1 gene, and bear numerous mutations. A phylogenetic analysis of the six largest sequences indicates that they are all derived from a common ancestor, which would have appeared 22-33 million years ago, and was subsequently duplicated several times during hominoid evolution. The most recent duplication and interchromosomal transposition occurred in the last million years suggesting that the process could still be ongoing. Intriguing similarities between the evolution of alpha-satellite DNA and NF1-related sequences suggest the involvement of a common genetic mechanism for the generation and pericentric spreading of these NF1 partial copies.

Phylogenetic analysis of the Monogenea and their relationships with Digenea and Eucestoda inferred from 28S rDNA sequences.

Platyhelminth phylogeny is controversial. Phylogenetic analyses of the partial domain C1 and the full domains D1 and C2 (358 nucleotides) from the 28S ribosomal RNA gene for 21 species from the Monogenea, Digenea, Cestoda, and, as the outgroup, Tricladida reveal major departures from prevailing theory. The Digenea and not the Monogenea (Monopisthocotylea and Polyopisthocotylea) form the sister group of the cestodes; the Monopisthocotylea and Polyopisthocotylea are each monophyletic, but the Monogenea do not form a monophylum; the sister group of the Digenea + Cestoda is the Polyopisthocotylea; and Monopisthocotylea are the sister group of all other parasitic flatworms.