Molecular phylogeny of 21 tropical bamboo species reconstructed by integrating non-coding internal transcribed spacer (ITS1 and 2) sequences and their consensus secondary structure

The unavailability of the reproductive structure and unpredictability of vegetative characters for the identification and phylogenetic study of bamboo prompted the application of molecular techniques for greater resolution and consensus. We first employed internal transcribed spacer (ITS1, 5.8S rRNA and ITS2) sequences to construct the phylogenetic tree of 21 tropical bamboo species. While the sequence alone could grossly reconstruct the traditional phylogeny amongst the 21-tropical species studied, some anomalies were encountered that prompted a further refinement of the phylogenetic analyses. Therefore, we integrated the secondary structure of the ITS sequences to derive individual sequence-structure matrix to gain more resolution on the phylogenetic reconstruction. The results showed that ITS sequence-structure is the reliable alternative to the conventional phenotypic method for the identification of bamboo species. The best-fit topology obtained by the sequence-structure based phylogeny over the sole sequence based one underscores closer clustering of all the studied Bambusa species (Sub-tribe Bambusinae), while Melocanna baccifera, which belongs to Sub-Tribe Melocanneae, disjointedly clustered as an out-group within the consensus phylogenetic tree. In this study, we demonstrated the dependability of the combined (ITS sequence+structure-based) approach over the only sequence-based analysis for phylogenetic relationship assessment of bamboo.     

Molecular phylogenetics in 2D: ITS2 rRNA evolution and sequence-structure barcode from Veneridae to Bivalvia

In this study, we analyzed the nuclear ITS2 rRNA primary sequence and secondary structure in Veneridae and comparatively with 20 Bivalvia taxa to test the phylogenetic resolution of this marker and its suitability for molecular diagnosis at different taxonomic levels. Maximum likelihood and Bayesian trees based on primary sequences were congruent with (profile-) neighbor-joining trees based on a combined model of sequence-structure evolution. ITS2 showed higher resolution below the subfamily level, providing a phylogenetic signal comparable to (mitochondrial/nuclear) gene fragments 2-5 times longer. Structural elements of the ITS2 folding, such as specific mismatch pairing and compensatory base changes, provided further support for the monophyly of some groups and for their phylogenetic relationships. Veneridae ITS2 folding is structured in six domains (DI-VI) and shows five striking sequence-structure features. Two of them, the Basal and Apical STEMs, are common to Bivalvia, while the presence of both the Branched STEM and the Y/R stretches occurs in five superfamilies of the two Heterodonta orders Myoida and Veneroida, thus questioning their reciprocal monophyly. Our results validated the ITS2 as a suitable marker for venerids phylogenetics and taxonomy, and underlined the significance of including secondary structure information for both applications at several systematic levels within bivalves.     

The ITS2 Database

The internal transcribed spacer 2 (ITS2) has been used as a phylogenetic marker for more than two decades. As ITS2 research mainly focused on the very variable ITS2 sequence, it confined this marker to low-level phylogenetics only. However, the combination of the ITS2 sequence and its highly conserved secondary structure improves the phylogenetic resolution¹ and allows phylogenetic inference at multiple taxonomic ranks, including species delimitation^2-8.The ITS2 Database⁹ presents an exhaustive dataset of internal transcribed spacer 2 sequences from NCBI GenBank¹¹ accurately reannotated¹⁰. Following an annotation by profile Hidden Markov Models (HMMs), the secondary structure of each sequence is predicted. First, it is tested whether a minimum energy based fold¹² (direct fold) results in a correct, four helix conformation. If this is not the case, the structure is predicted by homology modeling¹³. In homology modeling, an already known secondary structure is transferred to another ITS2 sequence, whose secondary structure was not able to fold correctly in a direct fold.The ITS2 Database is not only a database for storage and retrieval of ITS2 sequence-structures. It also provides several tools to process your own ITS2 sequences, including annotation, structural prediction, motif detection and BLAST¹⁴ search on the combined sequence-structure information. Moreover, it integrates trimmed versions of 4SALE^15,16 and ProfDistS¹⁷ for multiple sequence-structure alignment calculation and Neighbor Joining¹⁸ tree reconstruction. Together they form a coherent analysis pipeline from an initial set of sequences to a phylogeny based on sequence and secondary structure.In a nutshell, this workbench simplifies first phylogenetic analyses to only a few mouse-clicks, while additionally providing tools and data for comprehensive large-scale analyses.     

Multilocus sequence evaluation for differentiating species of the trematode Family Gastrothylacidae,with a note on the utility of mitochondrial COI motifs in species identification

Amphistomiasis, a neglected trematode infectious disease of ruminants, is caused by numerous species of amphistomes belonging to six families under the Superfamily Paramphistomoidea. In the present study, four frequently used DNA markers, viz. nuclear ribosomal 28S (D1–D3 regions), 18S and ITS2 and mitochondrial COI genes, as well as sequence motifs from these genes were evaluated for their utility in species characterization of members of the amphistomes' Family Gastrothylacidae commonly prevailing in Northeast India. In sequence and phylogenetic analyses the COI gene turned out to be the most useful marker in identifying the gastrothylacid species, with the exception of Gastrothylax crumenifer, which showed a high degree of intraspecific variations among its isolates. The sequence analysis data also showed the ITS2 region to be effective for interspecies characterization, though the 28S and 18S genes were found unsuitable for the purpose. On the other hand, sequence motif analysis data revealed the motifs from the COI gene to be highly conserved and specific for their target species which allowed accurate in silico identification of the gastrothylacid species irrespective of their intraspecific differences. We propose the use of COI motifs generated in the study as a potential tool for identification of these species.     

Molecular analyses of the genus Ilex (Aquifoliaceae) in southern South America, evidence from AFLP and ITS sequence data

In order to clarify the relationships among southern South American (sSA) representatives of the genus Ilex, an amplified fragment length polymorphism (AFLP) analysis was accomplished. In addition, the phylogenetic relationships of the species were studied using ribosomal internal transcribed spacer (ITS) sequence data alone and in combination with AFLP data, taking into account the possible existence of paralogous sequences and the influence of alignment parameters. To explore stability of phylogenetic hypotheses, a sensitivity analysis was performed using 15 indel-substitution models. Within each species assayed, the AFLPs allowed the recognition of several diagnostic bands. Furthermore, the AFLP analysis revealed that individuals belonging to the same morpho-species formed coherent clades. In addition, some cases of geographical association were noted. Studies on ITS sequences revealed divergence between data obtained herein and sequence data downloaded from GenBank. The sensitivity analyses yielded different interspecific hypotheses of relationships. Notwithstanding, analyses of the ITS data alone and in combination with AFLPs, rendered clades stable to variation in the analytical parameters. Topologies obtained for the AFLPs, the ITS data alone and the combined analyses, demonstrated the existence of a group formed by I. argentina, I. brasiliensis, I. brevicuspis, I. integerrima, and I. theezans, and that I. dumosa and I. paraguariensis were distantly related to the former. Incongruence with traditional taxonomical treatments was found.     

Species delimitation of Chinese hop‐hornbeams based on molecular and morphological evidence

Species delimitation through which infers species boundaries is emerging as a major work in modern systematics. Hop‐hornbeam species in Ostrya (Betulaceae) are well known for their hard and heavy woods. Five species were described in China and their interspecific delimitations remain unclear. In this study, we firstly explored their distributions in all recorded field sites distributed in China. We then selected 110 samples from 22 natural populations of five species from this genus and one type specimen of O. yunnanensis, for molecular barcoding analyses. We sequenced four chloroplast (cp) DNA fragments (trnH–psbA, trnL–trnF, rps16, and trnG) and the nuclear internal transcribed spacer (ITS) region for all samples. Sequence variations of Ostrya from four cpDNA fragments identified three groups that showed no correspondence to any morphological delimitation because of the incomplete lineage sorting and/or possible interspecific introgression in the history. However, phylogenetic analyses of ITS sequence variations discerned four species, O. japonica, O. rehderiana, O. trichocarpa, and O. multinervis while O. yunnanensis nested within O. multinervis. Morphological clustering also discerned four species and showed the complete consistency with molecular evidence. Moreover, our phylogenetic analyses‐based ITS sequence variations suggested that O. trichocarpa comprised an isolated lineage different from the other Eurasian ones. Based on these results, hop‐hornbeams in China should be treated as four separate species. Our results further highlight the importance of ITS sequence variations in delimitating and discerning the closely related species in plants.     

Molecular Evolution and Phylogenetic Utility of the Internal Transcribed Spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera)

The resolution potential of internal transcribed spacer 2 (ITS2) at deeper levels remains controversial. In this study, 105 ITS2 sequences of 55 species in Calyptratae were analyzed to examine the phylogenetic utility of the spacer above the subfamily level and to further understand its evolutionary characteristics. We predicted the secondary structure of each sequence using the minimum-energy algorithm and constructed two data matrixes for phylogenetic analysis. The ITS2 regions of Calyptratae display strong A-T bias and slight variation in length. The tandem and dispersed repeats embedded in the spacers possibly resulted from replication slippage or transposition. Most foldings conformed to the four-domain model. Sequence comparison in combination with the secondary structures revealed six conserved motifs. Covariation analysis from the conserved motifs indicated that the secondary structure restrains the sequence evolution of the spacer. The deep-level phylogeny derived from the ITS2 data largely agreed with the phylogenetic hypotheses from morphologic and other molecular evidence. Our analyses suggest that the accordant resolutions generated from different analyses can be used to infer deep-level phylogenetic relations.     

ITS2 database IV: interactive taxon sampling for internal transcribed spacer 2 based phylogenies

The first step of any molecular phylogenetic analysis is the selection of the species and sequences to be included, the taxon sampling. Already here different pitfalls exist. Sequences can contain errors, annotations in databases can be inaccurate and even the taxonomic classification of a species can be wrong. Usually, these artefacts become evident only after calculation of the phylogenetic tree. Following, the taxon sampling has to be corrected iteratively. This can become tedious and time consuming, as in most cases the taxon sampling is de-coupled from the further steps of the phylogenetic analysis. Here, we present the ITS2 Workbench (http://its2.bioapps.biozentrum.uni-wuerzburg.de/), which eliminates this problem by a tight integration of taxon sampling, secondary structure prediction, multiple alignment and phylogenetic tree calculation. The ITS2 Workbench has access to more than 280,000 ITS2 sequences and their structures provided by the ITS2 database enabling sequence-structure based alignment and tree reconstruction. This allows the interactive improvement of the taxon sampling throughout the whole phylogenetic tree reconstruction process. Thus, the ITS2 Workbench enables a fast, interactive and iterative taxon sampling leading to more accurate ITS2 based phylogenies.     

Distinguishing species

Given two organisms, how can one distinguish whether they belong to the same species or not? This might be straightforward for two divergent organisms, but can be extremely difficult and laborious for closely related ones. A molecular marker giving a clear distinction would therefore be of immense benefit. The internal transcribed spacer 2 (ITS2) has been widely used for low-level phylogenetic analyses. Case studies revealed that a compensatory base change (CBC) in the helix II or helix III ITS2 secondary structure between two organisms correlated with sexual incompatibility. We analyzed more than 1300 closely related species to test whether this correlation is generally applicable. In 93%, where a CBC was found between organisms classified within the same genus, they belong to different species. Thus, a CBC in an ITS2 sequence-structure alignment is a sufficient condition to distinguish even closely related species.     

The blood alga: phylogeny of Haematococcus (Chlorophyceae) inferred from ribosomal RNA gene sequence data

The status of the green algal genera Haematococcus and Stephanosphaera has been a source of debate among algal systematists. A phylogenetic alliance between Haematococcus (sensu lato) and the colonial Stephanosphaera was affirmed by earlier molecular phylogenetic investigations. Although the data suggested that the genus Haematococcus may not be a monophyletic group, taxon sampling limited the scope of any potential taxonomic revision. Results from new molecular phylogenetic analyses of data from the 18S and 26S rRNA genes support the establishment of a separate genus, Balticola, as originally proposed by Droop in 1956. Haematococcus remains as a valid genus, with H. pluvialis as its only member. The monotypic status of H. pluvialis is supported both by molecular phylogenetic analyses of the ribosomal RNA genes and assessments of molecular evolution in the ITS2 sequences of H. pluvialis strains. The near-complete absence of compensatory base changes in a sequence-structure analysis of the highly variable ITS2 gene from more than 40 geographically diverse isolates of H. pluvialis corroborates the unity of the species inferred from molecular phylogenetic analyses of 18S and 26S rRNA gene sequence data.     

Molecular Phylogenetics and Systematics of the Bivalve Family Ostreidae Based on rRNA Sequence-Structure Models and Multilocus Species Tree

The bivalve family Ostreidae has a worldwide distribution and includes species of high economic importance. Phylogenetics and systematic of oysters based on morphology have proved difficult because of their high phenotypic plasticity. In this study we explore the phylogenetic information of the DNA sequence and secondary structure of the nuclear, fast-evolving, ITS2 rRNA and the mitochondrial 16S rRNA genes from the Ostreidae and we implemented a multi-locus framework based on four loci for oyster phylogenetics and systematics. Sequence-structure rRNA models aid sequence alignment and improved accuracy and nodal support of phylogenetic trees. In agreement with previous molecular studies, our phylogenetic results indicate that none of the currently recognized subfamilies, Crassostreinae, Ostreinae, and Lophinae, is monophyletic. Single gene trees based on Maximum likelihood (ML) and Bayesian (BA) methods and on sequence-structure ML were congruent with multilocus trees based on a concatenated (ML and BA) and coalescent based (BA) approaches and consistently supported three main clades: (i) Crassostrea, (ii) Saccostrea, and (iii) an Ostreinae-Lophinae lineage. Therefore, the subfamily Crassotreinae (including Crassostrea), Saccostreinae subfam. nov. (including Saccostrea and tentatively Striostrea) and Ostreinae (including Ostreinae and Lophinae taxa) are recognized. Based on phylogenetic and biogeographical evidence the Asian species of Crassostrea from the Pacific Ocean are assigned to Magallana gen. nov., whereas an integrative taxonomic revision is required for the genera Ostrea and Dendostrea. This study pointed out the suitability of the ITS2 marker for DNA barcoding of oyster and the relevance of using sequence-structure rRNA models and features of the ITS2 folding in molecular phylogenetics and taxonomy. The multilocus approach allowed inferring a robust phylogeny of Ostreidae providing a broad molecular perspective on their systematics.     

Thysanophora penicillioides includes multiple genetically diverged groups that coexist respectively in Abies mariesii forests in Japan

We investigated intraspecific diversity and genetic structures of a saprotrophic fungus--Thysanophora penicillioides--based on sequences of nuclear ribosomal internal transcribed spacer (ITS) in 15 discontinuous Abies mariesii forests of Japan. In such a well-defined morphological species, numerous unexpected ITS variations were revealed: 12 ITS sequence types detected in 254 isolates collected from 15 local populations were classified into five ITS sequence groups. Maximally, four ITS groups consisted of seven ITS types coexisting in one population. However, group 1 was dominant with approximately 65%; in particular, one haplotype, 1a, was most dominant with approximately 60% in respective populations. Therefore, few differences were recognized in genetic structure among local populations, implying that the gene flow of each lineage of the fungus occurs among local populations without geographic limitations. However, minor haplotypes in some ITS groups were found only in restricted areas, suggesting that they might expand steadily from their places of origin to neighboring A. mariesii forests. Aggregating sequence data of seven European strains and four North American strains from various substrates to those of Japanese strains, 18 ITS sequence types and 28 variable sites were recognized. They were clustered into nine lineages by phylogenetic analyses of the beta-tubulin and combined ITS and beta-tubulin datasets. According to phylogenetic species recognition by the concordance of genealogies, respective lineages correspond to phylogenetic species. Plural phylogenetic species coexist in a local population in an A. mariesii forest in Japan.     

Internal transcribed spacer 2 (nu ITS2 rRNA) sequence-structure phylogenetics: towards an automated reconstruction of the green algal tree of life

Background

Chloroplast-encoded genes (matK and rbcL) have been formally proposed for use in DNA barcoding efforts targeting embryophytes. Extending such a protocol to chlorophytan green algae, though, is fraught with problems including non homology (matK) and heterogeneity that prevents the creation of a universal PCR toolkit (rbcL). Some have advocated the use of the nuclear-encoded, internal transcribed spacer two (ITS2) as an alternative to the traditional chloroplast markers. However, the ITS2 is broadly perceived to be insufficiently conserved or to be confounded by introgression or biparental inheritance patterns, precluding its broad use in phylogenetic reconstruction or as a DNA barcode. A growing body of evidence has shown that simultaneous analysis of nucleotide data with secondary structure information can overcome at least some of the limitations of ITS2. The goal of this investigation was to assess the feasibility of an automated, sequence-structure approach for analysis of IT2 data from a large sampling of phylum Chlorophyta.

Methodology/Principal Findings

Sequences and secondary structures from 591 chlorophycean, 741 trebouxiophycean and 938 ulvophycean algae, all obtained from the ITS2 Database, were aligned using a sequence structure-specific scoring matrix. Phylogenetic relationships were reconstructed by Profile Neighbor-Joining coupled with a sequence structure-specific, general time reversible substitution model. Results from analyses of the ITS2 data were robust at multiple nodes and showed considerable congruence with results from published phylogenetic analyses.

Conclusions/Significance

Our observations on the power of automated, sequence-structure analyses of ITS2 to reconstruct phylum-level phylogenies of the green algae validate this approach to assessing diversity for large sets of chlorophytan taxa. Moreover, our results indicate that objections to the use of ITS2 for DNA barcoding should be weighed against the utility of an automated, data analysis approach with demonstrated power to reconstruct evolutionary patterns for highly divergent lineages.     

松萝属的5.8S-ITS2二级结构比较分析

本研究以松萝属为研究对象,构建了该属核糖体5.8S和第二转录区间（ITS2）的二级结构模型,并对种间的结构差异进行了比较。结果显示,松萝属的结构与先前发表的真核生物ITS2二级结构模型非常相似。属内种间的结构差异主要集中在第二转录区间的臂上。核酸序列分析所产生的系统树与结构特征的比较结果相一致,表明ITS2的二级结构信息可以作为一种辅助的分类手段,适用于松萝属的种的鉴定和亲缘关系研究。     

Phylogenetic utility of the ribosomal internal transcribed spacer 2 in Strumigenys spp. (Hymenoptera: Formicidae)

Sequence variations and phylogenetic relationships of Strumigenys from different localities in Taiwan were inferred from sequences of the internal transcribed spacer 2 (ITS2) of nuclear ribosomal RNA genes. The ITS2 sequences of different species of Strumigenys vary in length from 659 to 953 bp, because there are many large repeated insertion-deletion fragments, and these short tandem repeat units form microsatellites. The average GC content of the ITS2 is 60.8%. Secondary structures from ITS2 sequences are similar and present several conserved structural motifs. Eleven species and three unnamed species were analyzed using both the maximum parsimony and maximum likelihood methods. Although diversity of the ITS2 sequence is high, these data can still be a potential tool for primary analysis of molecular phylogeny and biogeography of Strumigenys at the species level on islands around Taiwan.     

The ITS1-5.8S-ITS2 sequence region in the Musaceae: structure, diversity and use in molecular phylogeny

Genes coding for 45S ribosomal RNA are organized in tandem arrays of up to several thousand copies and contain 18S, 5.8S and 26S rRNA units separated by internal transcribed spacers ITS1 and ITS2. While the rRNA units are evolutionary conserved, ITS show high level of interspecific divergence and have been used frequently in genetic diversity and phylogenetic studies. In this work we report on the structure and diversity of the ITS region in 87 representatives of the family Musaceae. We provide the first detailed information on ITS sequence diversity in the genus Musa and describe the presence of more than one type of ITS sequence within individual species. Both Sanger sequencing of amplified ITS regions and whole genome 454 sequencing lead to similar phylogenetic inferences. We show that it is necessary to identify putative pseudogenic ITS sequences, which may have negative effect on phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstruction based on ITS sequence showed that the genus Musa is divided into two distinct clades--Callimusa and Australimusa and Eumusa and Rhodochlamys. Most of the intraspecific banana hybrids analyzed contain conserved parental ITS sequences, indicating incomplete concerted evolution of rDNA loci. Independent evolution of parental rDNA in hybrids enables determination of genomic constitution of hybrids using ITS. The observation of only one type of ITS sequence in some of the presumed interspecific hybrid clones warrants further study to confirm their hybrid origin and to unravel processes leading to evolution of their genomes.     

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae)

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.     

Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

The organization, structure, and nucleotide variability of the ribosomal repeat unit was compared among families, genera, and species of cockroaches (Insecta:Blattodea). Sequence comparisons and molecular phylogenetic analyses were used to describe rDNA repeat unit variation at differing taxonomic levels. A reverse similar 1200 bp fragment of the 28S rDNA sequence was assessed for its potential utility in reconstructing higher-level phylogenetic relationships in cockroaches. Parsimony and maximum likelihood analyses of these data strongly support the expected pattern of relationships among cockroach groups. The examined 5' end of the 28S rDNA is shown to be an informative marker for larger studies of cockroach phylogeny. Comparative analysis of the nucleotide sequences of the rDNA internal transcribed spacers (ITS1 and ITS2) among closely related species of Blattella and Periplaneta reveals that ITS sequences can vary widely in primary sequence, length, and folding pattern. Secondary structure estimates for the ITS region of Blattella species indicate that variation in this spacer region can also influence the folding pattern of the 5.8S subunit. These results support the idea that ITS sequences play an important role in the stability and function of the rRNA cluster.     

Secondary Structure Analyses of the Nuclear rRNA Internal Transcribed Spacers and Assessment of Its Phylogenetic Utility across the Brassicaceae (Mustards)

The internal transcribed spacers of the nuclear ribosomal RNA gene cluster, termed ITS1 and ITS2, are the most frequently used nuclear markers for phylogenetic analyses across many eukaryotic groups including most plant families. The reasons for the popularity of these markers include: 1.) Ease of amplification due to high copy number of the gene clusters, 2.) Available cost-effective methods and highly conserved primers, 3.) Rapidly evolving markers (i.e. variable between closely related species), and 4.) The assumption (and/or treatment) that these sequences are non-functional, neutrally evolving phylogenetic markers. Here, our analyses of ITS1 and ITS2 for 50 species suggest that both sequences are instead under selective constraints to preserve proper secondary structure, likely to maintain complete self-splicing functions, and thus are not neutrally-evolving phylogenetic markers. Our results indicate the majority of sequence sites are co-evolving with other positions to form proper secondary structure, which has implications for phylogenetic inference. We also found that the lowest energy state and total number of possible alternate secondary structures are highly significantly different between ITS regions and random sequences with an identical overall length and Guanine-Cytosine (GC) content. Lastly, we review recent evidence highlighting some additional problematic issues with using these regions as the sole markers for phylogenetic studies, and thus strongly recommend additional markers and cost-effective approaches for future studies to estimate phylogenetic relationships.     

Phylogenetic and biogeographical relationships among some holarctic frog lung flukes (Digenea: Haematoloechidae).

A phylogenetic study of 8 North American and European species of frog lung flukes belonging to Haematoloechus was conducted using approximately 850 to 1,000 bases of the intemal transcribed spacer region (ITS 1 + 5.8S + ITS 2) and 1,250 bases of the large subunit (LSU) of the nuclear ribosomal DNA. Adequate phylogenetic resolution could not be obtained from 5.8S or ITS 2 data. Analysis of ITS 1 data produced 2 equally parsimonious trees that differed only in the position of Haematoloechus breviplexus relative to H. medioplexus and H. varioplexus. Single, identical trees were produced by analysis of both LSU sequence data and a data set comprised of all ITS and LSU data. All trees demonstrated 3 distinct evolutionary lineages within the Holarctic Haematoloechus examined. The results confirmed the taxonomic validity of H. abbreviatus and demonstrated that the presence or absence of extracecal uterine loops is not a character meaningful to the recognition of evolutionary lineages or differentiation of genera. Examination of ITS sequence data revealed almost no intraspecific variation within 5 species of Haematoloechus and demonstrated an approximately 150-base indel common to the North American H. longiplexus and the European H. asper. Two of 3 clades revealed by the phylogenetic analyses are comprised of both European and North American species, indicating that lineages of Haematoloechus arose before the breakup of Laurasia and radiated after Eurasia and North America split. Within each of 3 evolutionary lineages, members share similar patterns of arthropod host specificity distinct from patterns found in the other lineages. This suggests that second intermediate host specificity may be a trait that has been conserved through evolutionary time.