Nested cladistic analysis, phylogeography and speciation in the Timarcha goettingensis complex (Coleoptera, chrysomelidae)

The Timarcha goettingensis complex is a monophyletic assemblage of closely related leaf beetles (Chrysomelidae), distributed from the north half of the Iberian Peninsula to Central Europe. Oligophagy, mountainous habitat and apterism are factors which are assumed to promote speciation in these beetles. We have used cytochrome oxidase subunit II mitochondrial DNA genealogies obtained from 31 sampling localities and a nested geographical distance analysis to assess the population structure and demographic factors explaining the geographical distributions of the mtDNA haplotypes in the T. goettingensis complex. The results show that there is a significant association between genetic structuring and geography. Inferences about the historical population processes in the species complex are discussed, being in general in accordance with contiguous range expansions and past fragmentations. The use of the cohesion species concept approach suggests the existence of several systematic ranks among the different T. goettingensis populations, which is in part supported by ecological traits such as trophic selection and altitudinal distribution.     

Phylogenetic relationship of ribosomal ITS2 and mitochondrial COI among diploid and triploid Paragonimus westermani isolates

We compared patterns of intraspecific polymorphism of two markers with contrasting modes of evolution, nuclear ribosomal DNA (rDNA) and mitochondrial DNA (mtDNA), in the lung fluke, diploid and triploid Paragonimus westermani from three geographical regions of Korea. The genetic distances between three populations of Korean diploid and triploid P. westermani showed no significant difference in the nucleotide sequences of the mitochondrial cytochrome c oxidase subunit I (mtCOI) and ribosomaal second internal transcribed spacer (ITS2) genes. A highly resolved strict-consensus tree was obtained that illustrated phylogenetically useful information of the ITS2 and mtCOI sequences from diploid and triploid P. westermani.     

Testing species delimitation in sympatric species complexes: the case of an African tropical tree, Carapa spp. (Meliaceae)

Plant species delimitation within tropical ecosystems is often difficult because of the lack of diagnostic morphological characters that are clearly visible. The development of an integrated approach, which utilizes several different types of markers (both morphological and molecular), would be extremely useful in this context. Here we have addressed species delimitation of sympatric tropical tree species that belong to Carapa spp. (Meliaceae) in Central Africa. We adopted a population genetics approach, sampling numerous individuals from three locations where sympatric Carapa species are known to exist. Comparisons between morphological markers (the presence or absence of characters, leaf-shape traits) and molecular markers (chloroplast sequences, ribosomal internal transcribed spacer region (ITS) sequences, and nuclear microsatellites) demonstrated the following: (i) a strong correlation between morphological and nuclear markers; (ii) despite substantial polymorphism, the inability of chloroplast DNA to discriminate between species, suggesting that cytoplasmic markers represent ineffective DNA barcodes; (iii) lineage sorting effects when using ITS sequences; and (iv) a complex evolutionary history within the genus Carapa, which includes frequent inter-specific gene flow. Our results support the use of a population genetics approach, based on ultra-polymorphic markers, to address species delimitation within complex taxonomic groups.     

Biogeographical patterns in Artomyces pyxidatus

Artomyces pyxidatus (Auriscalpiaceae) is a lignicolous, coralloid basidiomycete found throughout temperate regions of the Northern Hemisphere. Previous studies established that populations from the eastern United States, Sweden, and China were conspecific based on mating compatibility and enzyme profiles. In this study, mating compatibility was extended to include collections from Russia, Costa Rica, Mexico, and Utah. The molecular diversity of A. pyxidatus was examined by DNA sequence and restriction site analyses of the nuclear ribosomal internally transcribed spacer region (ITS1-5.8S-ITS2). A phylogenetic analysis of twelve isolates based on ITS sequences revealed a broad geographical pattern in which Eurasian isolates comprise a sister clade to North American isolates. North American isolates appear to be further subdivided into northeastern and southwestern clades. A survey of 255 A. pyxidatus isolates using restriction enzymes revealed variable RFLP patterns that follow similar geographical patterns.     

Mitochondrial and ribosomal internal transcribed spacer (ITS2) diversity of the African malaria vector Anopheles funestus

The pattern of sequence variation in the mitochondrial DNA cytochrome b gene (cyt-b) and ribosomal DNA internal transcribed spacer 2 (ITS2) was examined in Anopheles funestus from Senegal and Burkina Faso in West Africa and Kenya in East Africa. From both West African countries, samples included individuals hypothesized to represent reproductively isolated taxa based upon different karyotypes and behaviours. Analysis of the cyt-b data revealed high haplotypic diversity (86%) and an average pairwise difference per site of 0.42%. Sequence variation was not partitioned by geographical origin or karyotype class. The most common haplotype was sampled across Africa (approximately 6000 km). Analysis of the ITS2 data revealed one of the longest spacers yet found in anophelines (approximately 704 bp). In common with other anopheline ITS2 sequences, this one had microsatellites and frequent runs of individual nucleotides. Also in common with data from other anopheline ITS2 studies, the An. funestus sequences were almost monomorphic, with only two rare polymorphisms detected. The results from both markers are congruent and do not support the hypothesis of reproductively isolated chromosomal taxa within An. funestus. Whether the lack of support by mitochondrial DNA (mtDNA) and ribosomal DNA (rDNA) sequences is a result of the recent origin of the presumptive taxa, or of the absence of barriers to gene flow, remains to be elucidated, using more rapidly evolving markers such as microsatellites.     

Differentiation and grouping of isolates of the Ganoderma lucidum complex by random amplified polymorphic DNA-PCR compared with grouping on the basis of internal transcribed spacer sequences.

Laccate polypores of the Ganoderma lucidum species complex are widespread white rot fungi of economic importance, but isolates cannot be identified by traditional taxonomic methods. Parsimony analysis of nucleotide sequences from the internal transcribed spacers (ITS) of the ribosomal gene (rDNA) distinguished six lineages in this species complex. Each ITS lineage may represent one or more putative species. While some isolates have identical ITS sequences, all of them could be clearly differentiated by genetic fingerprinting using random amplified polymorphic DNA (RAPD). To investigate the suitability of RAPD markers for taxonomic identification and grouping of isolates of the G. lucidum complex, RAPD fragments (RAPDs) were used as phenotypic characters in numerical and parsimony analyses. Results show that data from RAPDS do not distinguish the same clades as ITS data do. Groupings based on analysis of RAPD data were very sensitive to the choice of the grouping method used, and no consistent grouping of isolates could be proposed. However, analysis with RAPDs did resolve several robust terminal clades containing putatively conspecific isolates, suggesting that RAPDs might be helpful for systematics at the lower taxonomic levels that are unresolved by ITS sequence data. The limitations of RAPDs for systematics are briefly discussed. The conclusion of this study is that ITS sequences can be used to identify isolates of the G. lucidum complex, whereas RAPDs can be used to differentiate between isolates having identical ITS sequences. The practical implications of these results are briefly illustrated.     

不同产区太子参的rDNA ITS区序列的比较

使用1对引物18SPl和26SP2对采自14个产区的太子参[Pseudostellaria heterophylla(Miq．)Pax ex Pax et Hoffm．]进行ITS基因的PCR扩增和测序。序列分析结果表明,14个产区太子参的ITSl片段长度为219—222bp,ITS2片段长度为235～236bp,5．8S片段长度为155—157bp.除江苏宜兴,江苏句容马梗,江苏南京老鹰山和江苏溧阳等4个产区的ITS序列碱基完全一致外,其他10个产区的ITS序列则有不同的变异,碱基变异数目(包括5．8S编码区)为1—17个。使用UPGMA法重建系统发生树,从分子生物学角度说明了它们的变异程度,为利用ITS区序列的差异鉴别不同产区的太子参提供了依据。     

Genetic diversity of nuclear ITS1-5.8S-ITS2 rDNA sequence in Clonorchis sinensis Cobbold, 1875 (Trematoda: Opisthorchidae) from the Russian Far East

The present study examined the molecular organisation and sequence variation in the nuclear ribosomal DNA (rDNA) region, including the two internal transcribed spacers (ITS1 and ITS2) and the 5.8S gene of the Clonorchis sinensis from the Russian Far East. The relevant sequences from other parts of this species' area were downloaded from GenBank. The results showed 100% identity for all investigated 5.8S-ITS2 rDNA sequences. In contrast, two levels of intraspecific variations were revealed in the complete ITS1 sequences. The intra-genomic variation resulted from a C/T polymorphism in a single position. The inter-individual differences between the ITS1 sequences were both due to nucleotide and size polymorphisms resulting from a varying number of five-nucleotide repeats and followed by two ITS1 length variants. These variant frequencies correlate with the clonorchiasis level in some geographical localities. ITS1 differences, both in the mutation profile and mutation localisation, were revealed between northern and southern geographical samples. The presence of GC boxes that are identical to known regulatory motifs in eukaryotes was detected within the ITS1 sub-repeats. The predicted secondary structures for ITS1 consist of two large branches, one of which was invariable, while another depended on ITS1 length. The predicted secondary structure for ITS2 includes four helices around the core. The main differences between C. sinensis and other opisthorchids were localised on the tops of helices 2, 3, and 4. A phylogenetic MST reconstruction subdivided all ITS1 sequences into two well differentiated clusters, each with the major widespread ribotype, and showed that ribotype diversity in both Russia and Korea is much lower than in China. The results obtained demonstrate the feasibility of complete ITS1 sequences in C. sinensis population genetics and can be considered as a basis for further studies of the parasite infection because they may help to elucidate the molecular mechanisms of pathogen evolution and adaptation.     

The role of disjunction and postglacial population expansion on phylogeographical history and genetic diversity of the circumboreal plant Chamaedaphne calyculata

In the last decade a number of studies has illustrated quite different phylogeographical patterns amongst plants with a northern present‐day geographical distribution, spanning the entire circumboreal region and/or circumarctic region and southern mountains. These works, employing several marker systems, have brought to light the complex evolutionary histories of this group. Here I focus on one circumboreal plant species, Chamaedaphne calyculata (leatherleaf), to unravel its phylogeographical history and patterns of genetic diversity across its geographical range. A survey of 29 populations with combined analyses of chloroplast DNA (cpDNA), internal transcribed spacer (ITS) and AFLP markers revealed structuring into two groups: Eurasian/north‐western North American, and north‐eastern North American. The present geographical distribution of C. calyculata has resulted from colonization from two putative refugial areas: east Beringia and south‐eastern North America. The variation of chloroplast DNA (cpDNA) and ITS sequences strongly indicated that the evolutionary histories of the Eurasian/north‐western North American and the north‐eastern North American populations were independent of each other because of a geographical disjunction in the distribution area and ice‐sheet history between north‐eastern and north‐western North America. Mismatch analysis using ITS confirmed that the present‐day population structure is the result of rapid expansion, probably since the last glacial maximum. The AFLP data revealed low genetic diversity of C. calyculata (P = 19.5%, H = 0.085) over the whole geographical range, and there was no evidence of loss of genetic diversity within populations in the continuous range, either at the margins or in formerly glaciated and nonglaciated regions. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 761–775.     

Genetic diversity of Rhizoctonia solani associated with potato tubers in France

The soilborne fungus Rhizoctonia solani is a pathogen of many plants and causes severe damage in crops around the world. Strains of R. solani from the anastomosis group (AG) 3 attack potatoes, leading to great yield losses and to the downgrading of production. The study of the genetic diversity of the strains of R. solani in France allows the structure of the populations to be determined and adapted control strategies against this pathogen to be established. The diversity of 73 French strains isolated from tubers grown in the main potato seed production areas and 31 strains isolated in nine other countries was assessed by phylogenetic analyses of (i) the internal transcribed spacer sequences (ITS1 and ITS2) of ribosomal RNA (rRNA), (ii) a part of the gene tef-1α and (iii) the total DNA fingerprints of each strain established by amplified fragment length polymorphism (AFLP). The determination of the AGs of R. solani based on the sequencing of the ITS region showed three different AGs among our collection (60 AG 3 PT, 8 AG 2-1 and 5 AG 5). Grouping of the strains belonging to the same AG was confirmed by sequencing of the gene tef-1α used for the first time to study the genetic diversity of R. solani. About 42% of ITS sequences and 72% of tef-1α sequences contained polymorphic sites, suggesting that the cells of R. solani strains contain several copies of ITS and the tef-1α gene within the same nucleus or between different nuclei. Phylogenetic trees showed a greater genetic diversity within AGs in tef-1α sequences than in ITS sequences. The AFLP analyses showed an even greater diversity among the strains demonstrating that the French strains of R. solani isolated from potatoes were not a clonal population. Moreover there was no relationship between the geographical origins of the strains or the variety from which they were isolated and their genetic diversity.     

Presence of a short repeat sequence in internal transcribed spacer (ITS) 1 of the rRNA gene of <Emphasis Type="Italic">Sogatella furcifera</Emphasis> (Hemiptera: Delphacidae) from geographically different populations in Asia

Nucleotide sequences of the internal transcribed spacer 1 (ITS1)–5.8S–ITS2 region of the nuclear ribosomal RNA gene were determined in the white-backed planthopper (WBPH) Sogatella furcifera (Horváth) to detect molecular variation among regional populations in Asia. We analyzed 932 sequences from 172 individuals (4–9 clones per individual) of 33 populations collected in 1987–2008 from six countries, Japan, China, Taiwan, Vietnam, Philippines, and Papua New Guinea. WBPH showed intra-individual variation in ITS1, which is mainly attributable to the frequency (0–10) of the 66-bp repeat sequence in ITS1. Among the examined clones, the sequences of 5.8S were mostly identical and those of ITS2 were similar. A single planthopper had a maximum of 6 different variants in the number of ITS1 repeats, suggesting highly varied repeat numbers in individual planthoppers. The ITS1 with four repeats was the most frequently (64%) detected. Such a repeat was not observed in two other economically important planthopper species, Nilaparvata lugens (Stål) and Laodelphax striatellus (Fallén). The ITS nucleotide sequences in the WBPH populations in Asia were genetically close and some variations in the sequences were not related to regional populations, indicating that the nucleotide sequences of the ITS region are not useful for geographical discrimination of the WBPH. This closeness seems to be caused by long distance migration and genetic exchange among populations.     

Comparative analysis of population structure across environments and geographical scales at major histocompatibility complex and microsatellite loci in Atlantic salmon (Salmo salar)

Evidence of selection acting on major histocompatibility complex (MHC) genes has been illustrated with the analysis of their nucleotide sequences and allele frequency distribution. Comparing the patterns of population differentiation at neutral markers and MHC genes in the wild may provide further insights about the relative role of selection and neutrality in shaping their diversity. In this study, we combine both methods to assess the role of selection on a MHC gene in Atlantic salmon. We compare variation at a MHC class II B locus and microsatellites among 14 samples from seven different rivers and seven subpopulations within a single river system covering a variety of habitats and different geographical scales. We show that diversifying selection is acting on the sites involved in antigen presentation and that balancing selection maintains a high level of polymorphism within populations. Despite important differences in habitat type, the comparison of the population structure at MHC and microsatellites on large geographical scales reveals a correlation between patterns of differentiation, indicating that drift and migration have been more important than selection in shaping population differentiation at the MHC locus. In contrast, strong discrepancies between patterns of population differentiation at the two types of markers provides support for the role of selection in shaping population structure within rivers. Together, these results confirm that natural selection is influencing MHC gene diversity in wild Atlantic salmon although neutral forces may also be important in their evolution.     

Primers are designed for amplification and direct sequencing of ITS region of rDNA from Myxomycetes

Four new primers were designed, based on comparison of Physarum polycephalum sequences retrieved from Genbank (primers PHYS-5 and PHYS-4) and our own sequences (primers PHYS-3 and PHYS-2), to amplify the ITS regions of rDNA, including the 5.8S gene segment from Lamproderma species. Sequencing analysis shows that Lamproderma contains ITS1-5.8S-ITS2 regions of approximately 900 bp, which is similar in size to most eukaryotes. However, the corresponding region in another common myxomycete, Fuligo septica, is more than 2000 bp due to the presence of large direct-repeat motifs in ITS1. Myxomycete rDNA ITS regions are interesting both as phylogenetic markers in taxonomic studies and as model sequences for molecular evolution.     

Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome C oxidase I sequences

Studies of intraspecific genetic diversity of ciliates, such as population genetics and biogeography, are particularly hampered by the lack of suitable DNA markers. For example, sequences of the non-coding ribosomal internal transcribed spacer (ITS) regions are often too conserved for intraspecific analyses. We have therefore identified primers for the mitochondrial cytochrome c oxidase I (COI) gene and applied them for intraspecific investigations in Paramecium caudatum and Paramecium multimicronucleatum. Furthermore, we obtained sequences of the ITS regions from the same strains and carried out comparative sequence analyses of both data sets. The mitochondrial sequences revealed substantially higher variation in both Paramecium species, with intraspecific divergences up to 7% in P. caudatum and 9.5% in P. multimicronucleatum. Moreover, an initial survey of the population structure discovered different mitochondrial haplotypes of P. caudatum in one pond, thereby demonstrating the potential of this genetic marker for population genetic analyses. Our primers successfully amplified the COI gene of other Paramecium. This is the first report of intraspecific variation in free-living protozoans based on mitochondrial sequence data. Our results show that the high variation in mitochondrial DNA makes it a suitable marker for intraspecific and population genetic studies.     

Phylogenetic diversity of Synechococcus strains isolated from the East China Sea and the East Sea

Phylogenetic relationships among 33 Synechococcus strains isolated from the East China Sea (ECS) and the East Sea (ES) were studied based on 16S rRNA gene sequences and 16S–23S rRNA gene internal transcribed spacer (ITS) sequences. Pigment patterns of the culture strains were also examined. Based on 16S rRNA gene and ITS sequence phylogenies, the Synechococcus isolates were clustered into 10 clades, among which eight were previously identified and two were novel. Half of the culture strains belonged to clade V or VI. All strains that clustered into novel clades exhibited both phycoerythrobilin and phycourobilin. Interestingly, the pigment compositions of isolates belonging to clades V and VI differed from those reported for other oceanic regions. None of the isolates in clade V showed phycourobilin, whereas strains in clade VI exhibited both phycourobilin and phycoerythrobilin, which is in contrast to previous studies. The presence of novel lineages and the different pigment patterns in the ECS and the ES suggests the possibility that some Synechococcus lineages are distributed only in geographically restricted areas and have evolved in these regions. Therefore, further elucidation of the physiological, ecological, and genetic characteristics of the diverse Synechococcus strains is required to understand their spatial and geographical distribution.     

Identification of Paramecium bursaria syngens through molecular markers--comparative analysis of three loci in the nuclear and mitochondrial DNA

This is the first attempt to resolve the phylogenetic relationship between different syngens of Paramecium bursaria and to investigate at a molecular level the intraspecific differentiation of strains originating from very distant geographical locations. Herein we introduce a new collection of five P. bursaria syngens maintained at St Petersburg State University, as the international collection of syngens was lost in the 1960s. To analyze the degree of speciation within Paramecium bursaria, we examined 26 strains belonging to five different syngens from distant and geographically isolated localities using rDNA (ITS1-5.8S-ITS2-5'LSU) fragments, mitochondrial cytochrome c oxidase subunit I (COI), and H4 gene fragments. It was shown that P. bursaria strains of the same syngens cluster together in all three inferred molecular phylogenies. The genetic diversity among the studied P. bursaria strains based on rDNA sequences was rather low. The COI divergence of Paramecium bursaria was also definitely lower than that observed in the Paramecium aurelia complex. The nucleotide sequences of the H4 gene analyzed in the present study indicate the extent of genetic differences between the syngens of Paramecium bursaria. Our study demonstrates the diagnostic value of molecular markers, which are important tools in the identification of Paramecium bursaria syngens.     

Population genetic structure of Tomicus piniperda L. (Curculionidae: Scolytinae) on different pine species and validation of T. destruens (Woll.)

Genetic diversity and population structure of Tomicus piniperda was assessed using mitochondrial sequences on 16 populations sampled on 6 pine species in France. Amplifications of Internal transcribed space 1 (ITS1) were also performed. Our goals were to determine the taxonomic status of the Mediterranean ecotype T. piniperda destruens, and to test for host plant or geographical isolation effect on population genetic structure. We showed that T. piniperda clusters in two mtDNA haplotypic groups. Clade A corresponds to insects sampled in continental France on Pinus sylvestris, P. pinaster and P. uncinata, whereas clade B gathers the individuals sampled in Corsica on P. pinaster and P. radiata and in continental France on P. pinea and P. halepensis. Insects belonging to clade A and clade B also consistently differ in the length of ITS1. Individuals belonging to both clades were found once in sympatry on P. pinaster. Genetic distances between clades are similar to those measured between distinct species of Tomicus. We concluded that clade B actually corresponds to the destruens ecotype and forms a good species, T. destruens. Analyses of molecular variance ( amova ) were conducted separately on T. destruens and T. piniperda to test for an effect of either geographical isolation or host species. Interestingly, the effect of host plant was significant for T. piniperda only, while the effect of geographical isolation was not. Pine species therefore seems to act as a significant barrier to gene flow, even if host race formation is not observed. These results still need to be confirmed by nuclear markers.     

The Influence of Pleistocene Climatic Changes and Ocean Currents on the Phylogeography of the Southern African Barnacle,Tetraclita serrata (Thoracica; Cirripedia)

The evolutionary effects of glacial periods are poorly understood for Southern Hemisphere marine intertidal species, particularly obligatory sessile organisms. We examined this by assessing the phylogeographic patterns of the southern African volcano barnacle, Tetraclita serrata, a dominant species on rocky intertidal shores. Restricted gene flow in some geographical areas was hypothesized based on oceanic circulation patterns and known biogeographic regions. Barnacle population genetic structure was investigated using the mitochondrial cytochrome oxidase subunit 1 (COI) region for 410 individuals sampled from 20 localities spanning the South African coast. The mtDNA data were augmented by generating nuclear internal transcribed spacer 1 (ITS1) sequences from a subset of samples. Phylogenetic and population genetic analyses of mitochondrial DNA data reveal two distinct clades with mostly sympatric distributions, whereas nuclear analyses reveal only a single lineage. Shallow, but significant structure (0.0041–0.0065, P<0.01) was detected for the mtDNA data set, with the south-west African region identified as harbouring the highest levels of genetic diversity. Gene flow analyses on the mtDNA data show that individuals sampled in south-western localities experience gene flow primarily in the direction of the Benguela Current, while south and eastern localities experience bi-directional gene flow, suggesting an influence of both the inshore currents and the offshore Agulhas Current in the larval distribution of T. serrata. The mtDNA haplotype network, Bayesian Skyline Plots, mismatch distributions and time since expansion indicate that T. serrata population numbers were not severely affected by the Last Glacial Maximum (LGM), unlike other southern African marine species. The processes resulting in the two morphologically cryptic mtDNA lineages may be the result of a recent historical allopatric event followed by secondary contact or could reflect selective pressures due to differing environmental conditions.     

Internal transcribed spacer repeat-specific primers and the analysis of hybridization in the Glycine tomentella (Leguminosae) polyploid complex

Polyploid and diploid hybridization is a ubiquitous and evolutionarily important phenomenon in the plant world. Determining the parental species of a hybrid, however, is difficult. Molecular markers such as the nuclear ribosomal DNA gene complex, particularly its internal transcribed spacer (ITS) region, have proved powerful in determining hybrid parentage. In some cases, population and genomic phenomena, such as genetic drift and concerted evolution, result in the loss of all or many of the tandemly repeated copies derived from one parental species, making the recovery of hybrid history difficult or impossible. Methods such as direct sequencing and cloning are typically used to find ITS sequences contributed from parental species, but are limited in their ability to detect rare repeat types. Here we report that repeat-specific polymerase chain reaction primers can recover rare parental ITS sequences in the Glycine tomentella polyploid complex. In three allopolyploid lineages of this complex, repeat-specific primers reliably detected rare repeats that both direct sequencing and the screening of many cloned sequences failed to detect. Other strategies, such as the use of exclusion primers, may detect rare parental repeat types in hybrids when previous hypotheses regarding the second parental species are lacking.     

A molecular phylogeographic study based on DNA sequences from individual metacercariae of Paragonimus mexicanus from Guatemala and Ecuador

A molecular phylogeographic study of Paragonimus mexicanus collected from Guatemala and Ecuador was performed. Genomic DNA was extracted from individual metacercariae, and two gene regions (partial mitochondrial cytochrome c oxidase subunit 1 (CO1) and the second internal transcribed spacer of the nuclear ribosomal gene repeat (ITS2)) were amplified by the polymerase chain reaction (PCR). Sequences segregated in a phylogenetic tree according to their geographic origins. ITS2 sequences from Ecuador and Guatemala differed at only one site. Pairwise distances among CO1 sequences within a country were always lower than between countries. Nevertheless, genetic distances between countries were less than between geographical forms of P. westermani that have been suggested to be distinct species. This result suggests that populations from Guatemala and Ecuador are genetically differentiated perhaps at the level of subspecies.