Population structure and species boundary delimitation of cryptic Dioryctria moths: an integrative approach

Accurate delimitation of species boundaries is especially important in cryptic taxa where one or more character sources are uninformative or are in conflict. Rather than relying on a single marker to delimit species, integrative taxonomy uses multiple lines of evidence such as molecular, morphological, behavioural and geographic characters to test species limits. We examine the effectiveness of this approach by testing the delimitation of two cryptic Nearctic species of Dioryctria (Lepidoptera: Pyralidae) using three independent molecular markers [cytochrome c oxidase I (COI), second internal transcribed spacer unit (ITS2), and elongation factor 1alpha (EF1alpha)], forewing variation and larval host plant association. Although mitochondrial DNA (mtDNA) haplotypes do not form reciprocally monophyletic clades, restricted gene flow between COI haplotype groups, and concordance with ITS2 genotypes, forewing variation and host plant associations support delimitation of two Nearctic species: eastern Dioryctria reniculelloides and western Dioryctria pseudotsugella. Conversely, EF1alpha genotype variation was incongruent with the two previous markers. A case of discordance between COI and ITS2 was detected, suggesting either introgression due to hybridization or retained ancestral polymorphism due to incomplete coalescence. This study is consistent with other similar literature where molecular loci in closely related species progress from shared to fixed haplotypes/alleles, and from polyphyletic to reciprocally monophyletic relationships, although loci may vary in these characteristics despite maintenance of genomic integrity between distinct species. In particular, mtDNA in other studies generally showed a lower rate of fixation of differences than did X-linked or autosomal loci, reinforcing the need to use an integrative approach for delimiting species.     

Mitochondrial and nuclear DNA analysis revealed a cryptic species and genetic introgression in Littorina sitkana (Mollusca, Gastropoda)

We investigated mitochondrial and nuclear DNA genotypes in nominal Littorina sitkana samples from 2 localities in Eastern Hokkaido, northern Japan. Our results indicated the existence of cryptic species. In the analysis of partial mitochondrial Cytchrome b gene sequences, haplotypes of L. sitkana samples were monophyletic in a phylogenetic tree with orthologous sequences from other Littorina species, but were apparently separated in 2 clades. One included typical L. sitkana (CBa clade) samples, which formed a clade with an allopatric species, L. horikawai. The other, CBb, was independent from CBa and L. horikawai. Haplotypes of the mitochondrial 16S rRNA gene also separated into 2 clades. We additionally examined intron sequence of the heat shock cognate 70 (HSC70) nuclear gene and identified 17 haplotypes. These were also separated into 2 clades, HSCa and HSCb. Among the examined Hokkaido samples, 60% of individuals were heterozygotes. However, each heterozygote consisted of haplotypes from the same clade, HSCa or HSCb, and no admixture of HSCa and HSCb haplotypes was observed. These results indicate reproductive isolation between the 2 clades. Among the genotyped Hokkaido samples, 93% of individuals had CBa?+?HSCa or CBb?+?HSCb genotypes, and 7% had CBb?+?HSCa genotypes. The discrepancy between the mtDNA and nuclear DNA haplotypes in a few individuals may have been caused by genetic introgression due to past hybridization.     

Phylogeny and temporal diversification of darters (Percidae: Etheostomatinae)

Discussions aimed at resolution of the Tree of Life are most often focused on the interrelationships of major organismal lineages. In this study, we focus on the resolution of some of the most apical branches in the Tree of Life through exploration of the phylogenetic relationships of darters, a species-rich clade of North American freshwater fishes. With a near-complete taxon sampling of close to 250 species, we aim to investigate strategies for efficient multilocus data sampling and the estimation of divergence times using relaxed-clock methods when a clade lacks a fossil record. Our phylogenetic data set comprises a single mitochondrial DNA (mtDNA) gene and two nuclear genes sampled from 245 of the 248 darter species. This dense sampling allows us to determine if a modest amount of nuclear DNA sequence data can resolve relationships among closely related animal species. Darters lack a fossil record to provide age calibration priors in relaxed-clock analyses. Therefore, we use a near-complete species-sampled phylogeny of the perciform clade Centrarchidae, which has a rich fossil record, to assess two distinct strategies of external calibration in relaxed-clock divergence time estimates of darters: using ages inferred from the fossil record and molecular evolutionary rate estimates. Comparison of Bayesian phylogenies inferred from mtDNA and nuclear genes reveals that heterospecific mtDNA is present in approximately 12.5% of all darter species. We identify three patterns of mtDNA introgression in darters: proximal mtDNA transfer, which involves the transfer of mtDNA among extant and sympatric darter species, indeterminate introgression, which involves the transfer of mtDNA from a lineage that cannot be confidently identified because the introgressed haplotypes are not clearly referable to mtDNA haplotypes in any recognized species, and deep introgression, which is characterized by species diversification within a recipient clade subsequent to the transfer of heterospecific mtDNA. The results of our analyses indicate that DNA sequences sampled from single-copy nuclear genes can provide appreciable phylogenetic resolution for closely related animal species. A well-resolved near-complete species-sampled phylogeny of darters was estimated with Bayesian methods using a concatenated mtDNA and nuclear gene data set with all identified heterospecific mtDNA haplotypes treated as missing data. The relaxed-clock analyses resulted in very similar posterior age estimates across the three sampled genes and methods of calibration and therefore offer a viable strategy for estimating divergence times for clades that lack a fossil record. In addition, an informative rank-free clade-based classification of darters that preserves the rich history of nomenclature in the group and provides formal taxonomic communication of darter clades was constructed using the mtDNA and nuclear gene phylogeny. On the whole, the appeal of mtDNA for phylogeny inference among closely related animal species is diminished by the observations of extensive mtDNA introgression and by finding appreciable phylogenetic signal in a modest sampling of nuclear genes in our phylogenetic analyses of darters.     

Species boundaries in Malagasy snakes of the genus Madagascarophis (Serpentes: Colubridae sensu lato) assessed by nuclear and mitochondrial markers

We studied mitochondrial divergence in 27 individuals of colubrid snakes of the genus Madagascarophis Mertens from most of its distribution area in Madagascar. Combined analyses of 16S rRNA and cytochrome b sequences identified six major clades which only partly agreed with previously proposed classifications. Analysis of nuclear DNA sequences of the c-mos gene as well as of ISSR fingerprints revealed consistent differences only among three clades which we consider as distinct species: a widespread Madagascarophis colubrinus (Schlegel), with M. citrinus (Boettger) as a junior synonym, a southern M. meridionalis Domergue, and a presumably undescribed species from the extreme north of Madagascar. The species M. ocellatus Domergue was not available for our study. Within M. colubrinus there are two populations from the north-west, each showing two divergent haplotypes with pairwise divergences of up to 5.2% in the cytochrome b gene. Maximum divergence in this gene within M. colubrinus was 7.1%. These high values emphasise that caution needs to be applied before genetic distance values are used for species delimitation. Phylogeographically, most of the genetic variation in M. colubrinus is found in northern Madagascar, indicating that the species might have originated in this region. Later one haplotype clade colonised western and eastern Madagascar, with a putative secondary introgression into north-western populations.     

Phylogeographic analysis of mitochondrial gene flow and introgression in the salamander, Plethodon shermani

Plethodon shermani comprises a series of geographically disjunct populations occupying high-elevation mountain isolates. These populations hybridize at their borders with salamanders of the Plethodon glutinosus species complex, and past range expansions inferred from Pleistocene climatic cycles may have increased the possible genetic interactions between P. shermani and species of the P. glutinosus complex. Because mitochondrial DNA haplotypes often show introgression across species borders, we survey mtDNA variation for evidence of past and ongoing genetic interactions between P. shermani, its close relative Plethodon cheoah, and species of the P. glutinosus complex. Ongoing hybridization with the P. glutinosus-complex species Plethodon teyahalee is accompanied by extensive mitochondrial introgression in some Unicoi populations of P. shermani, but it has little genetic impact on P. shermani populations outside hybrid zones at three other isolates (Tusquitee, Wayah Bald, Standing Indian). Some Unicoi populations of P. shermani exhibit mtDNA evidence of past hybridization with diverse lineages from P. aureolus and P. glutinosus. The Tusquitee isolate of P. shermani is also characterized by mtDNA haplotypes most closely related to Plethodon aureolus and P. glutinosus, presumably introduced by past genetic contact with these species or with introgressed populations of Unicoi P. shermani. The mtDNA variation in sampled populations of the Wayah Bald and Standing Indian isolates of P. shermani appears largely unaffected by ongoing hybridization. Principal components analyses of allozymic data indicate that P. shermani isolates collectively form a genetically homogeneous unit clearly demarcated from species with which they have had current or past genetic interactions. Rapid mtDNA introgression associated with transient contacts between P. shermani and other species permits a fine-level resolution of evolutionary lineages not evident from allozymic data.     

Genetic subdivision and biogeography of the Danubian rheophilic barb Barbus petenyi inferred from phylogenetic analysis of mitochondrial DNA variation

The barb Barbus petenyi is a cyprinid widely distributed throughout the mountain regions in the Danube River basin. Phylogenetic analysis of the DNA sequence variation at the mitochondrial cytochrome b gene over much of this range yielded three deep-branching (5.9-9.4% average divergence), well-supported haplotype clades with mutually exclusive geographic distributions and divergence times estimated to be in the Tertiary. The clades did not form an altogether monophyletic group as the most divergent one coalesced more recently with haplotypes of phylogenetically close species than with the other B. petenyi haplotypes. This pattern was supported by bootstrap and log-likelihood Shimodaira-Hasegawa tests. The other two were sister clades, but their distinctiveness was supported by previous allozyme data. Hence, from a taxonomic point of view, the current recognition of B. petenyi is erroneous, as it does not represent a single evolutionary lineage, and we suggest that three species be recognized instead. Substantial phylogeographic differences were evident among the three putative species, the two more southerly ones displaying significant structure, which suggested that they each survived in several glacial refugia throughout the Pleistocene. The phylogeographic pattern of multiple populations of rheophilic barbs with a history of long-term persistence and separation within the Danube River basin is novel within fishes and provides a hypothesis against which phylogeographic patterns among other similarly distributed rheophilic species may be compared.     

Phylogeography of two European newt species--discordance between mtDNA and morphology

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5-1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long-term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.     

Widespread decoupling of mtDNA variation and species integrity in Grammia tiger moths (Lepidoptera: Noctuidae)

Abstract. We investigate the diversity of the North American tiger moth genus Grammia Rambur (Lepidoptera: Noctuidae) by comparing mitochondrial DNA (mtDNA) ‘barcode’ fragments of cytochrome oxidase I with non‐molecular characters such as morphology, ecology, behaviour and distribution. Mitochondrial DNA genealogy is strikingly at odds with morpho‐species taxonomy for most of the 28 sampled species, as haplotypic polyphyly not only is taxonomically widespread, but involves multiple shared haplotypes among two to four species. Morpho‐ecological traits show that those species sharing haplotypes are often not closely related. Furthermore, high mtDNA divergences occur within species. Haplotypic variation is highly discordant with species taxonomy, but variation at a continental scale reveals significant geographic structuring of haplogroups, transcending morpho‐species boundaries. A nested clade analysis and comparison of non‐molecular with mtDNA data indicate that most discordance between mtDNA and taxonomy in Grammia is explained best by taxonomically and geographically widespread ongoing hybridization events resulting in mtDNA introgression. We hypothesize that broad areas of sympatry, interspecifically compatible genitalic structure, and species overlap in pheromone components facilitate hybridization, with disparate interspecies abundances promoting mitochondrial introgression. The molecular evolution of Grammia challenges the view that interspecific gene exchange occurs rarely and is restricted to recently diverged species. These results show the value of mtDNA in detecting cryptic hybridization, while highlighting the inherent dangers of drawing taxonomic conclusions based solely on mtDNA.     

Evidence for historical introgression along a contact zone between two species of char (Pisces: Salmonidae) in northwestern North America

Phylogeographic analyses can yield valuable insights into the geographic and historical contexts of contact and hybridization between taxa. Two species of char (Salmonidae), Dolly Varden (Salvelinus malma) and bull trout (S. confluentus) have largely parapatric distributions in watersheds of northwestern North America. They are, however, sympatric in several localities and hybridization and some introgression occurs across a broad area of contact. We conducted a comparative phylogenetic analysis of Dolly Varden and bull trout to gain a historical perspective of hybridization between these species and to test for footprints of historical introgression. We resolved two major Dolly Varden mitochondrial DNA (mtDNA) clades (with 1.4-2.2% sequence divergence between haplotypes) that had different geographical distributions. Clade N is distributed across most of the range of Dolly Varden, from southern British Columbia through to the Kuril Islands in Asia. Clade S had a much more limited distribution, from Washington state, at the southern limit of the Dolly Varden range, to the middle of Vancouver Island. The distribution and inferred ages of the mtDNA clades suggested that Dolly Varden survived the Wisconsinan glaciation in a previously unsuspected refuge south of the ice sheet, and that Dolly Varden and bull trout were probably in continuous contact over most of the last 100,000 years. When bull trout were included in the phylogenetic analysis, however, the mtDNA of neither species was monophyletic: Clade S Dolly Varden clustered within the bull trout mtDNA clade. This pattern was discordant with two nuclear phylogenies produced (growth hormone 2 and rRNA internal transcribed sequence 1), in which Dolly Varden and bull trout were reciprocally monophyletic. This discordance between mtDNA- and nDNA-based phylogenies indicates that historical introgression of bull trout mtDNA into Dolly Varden occurred. Percent sequence divergence within these introgressed Dolly Varden (clade S) was 0.2-0.6%, implying that the introgression occurred prior to the most recent glaciation. Our analysis and other evidence of contact between divergent lineages in northwestern North America strongly suggests that the area may be the site of previously unsuspected suture zones of aquatic biotas.     

Frequent cytoplasmic exchanges between oak species that are not closely related: Quercus suber and Q. ilex in Morocco

Chloroplast (cp) and mitochondrial (mt) DNA variation were studied in 97 populations of cork oak (Quercus suber) in Morocco; in 31 of these populations, holm oak (Quercus ilex), a clearly distinct species, also occurred and was compared with Q. suber. Three cpDNA and one mtDNA primer pairs were used in the survey, each in combination with one restriction enzyme. Six haplotypes belonging to two very divergent lineages were detected; one lineage predominates in each species, and is probably ancestral, as inferred from comparisons with other oak species. In the mixed-species populations, cytoplasmic genomes were frequently shared across species, as indicated by an introgression ratio of 0.63. This index is a new measure of the propensity of species to share locally genetic markers, varying from zero (complete differentiation) to one (no differentiation). By contrast, more closely related deciduous oak species (Q. robur, Q. petraea and Q. pubescens) have introgression ratios varying from 0.82 to 0.97. The introgression events appear to have been more frequent in the direction Q. ilex (female) x Q. suber (male), a finding which seems attributable to the flowering phenology of these two species. This asymmetry may have favoured immigration of Q. suber beyond its main range, in regions already colonized by Q. ilex. There, rare hybridization and further introgression through long distance pollen flow have established populations that are morphologically indistinguishable from Q. suber but that have cytoplasmic genomes originating from the local Q. ilex populations.     

Phylogeographical patterns and genetic diversity in three species of Eurasian boreal forest beetles

Many species that occur in formerly glaciated areas of Fennoscandia have reached their current ranges from glacial refugial areas in Eurasia. Little is known of the refugia and postglacial colonization routes of insect species that are confined to boreal forests. Here, we investigate the phylogeography of three species of saproxylic beetles distributed across Eurasia: two rare boreal forest specialists, Pytho kolwensis and Pytho abieticola , and a common, less specialized species, Pytho depressus . In all species, there were two well-defined haplotype clades based on 645 bp of cytochrome oxidase subunit I gene sequence. In each species one clade was found only in China. The other clade occurred from China to north-western Europe in both P. kolwensis and P. depressus , but was apparently absent from China in P. abieticola. In spite of common phylogeographical patterns, the distribution of genetic variation differed markedly between the three species. In P. kolwensis , a highly-threatened species in old-growth forests in Fennoscandia, there was an extremely low level of genetic variation throughout Eurasia. One common haplotype, represented by 86% of the samples, dominated in all sampling localities. Levels of genetic variation were higher in both P. abieticola and P. depressus , with 31% and 58%, respectively, of the samples representing a unique haplotype. In each species, relationships between haplotypes were not well resolved, and haplotypes from one sampling locality were generally not clustered in either Neighbour-joining trees or statistical parsimony networks. These patterns in the distribution of genetic variation can be attributed to differences in the species' population sizes, ecologies, glacial refugial areas, and postglacial colonization dynamics.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 267–279.     

Phylogenetic Relationships and Species Delimitation in Pinus Section Trifoliae Inferrred from Plastid DNA

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.     

Cladistic analysis of chloroplast DNA restriction site characters in Andira (Leguminosae: Dalbergieae)

Chloroplast DNA (cpDNA) restriction site variation was examined in 29 accessions, representing 29 populations of 17 species of Andira (Leguminosae, Papilionoideae, tribe Dalbergieae). This sample spans the morphological and geographical diversity of this poorly known, woody genus of approximately 30 species, which occurs throughout the neotropics, with one species in Africa. Cladistic analysis of 38 restriction site characters generated 182 equally most-parsimonious trees. The strict consensus tree preserved four well-supported clades, which are incongruent with previous infrageneric classifications of Andira. None of these clades has been recognized by previous workers, and three are not marked by any major morphological innovations. There were few restriction site differences between closely related species, but two (possibly three) species showed intraspecific cpDNA polymorphism, which probably reflects processes of hybridization with subsequent introgression, or lineage sorting.     

Combining mitochondrial DNA sequences and morphological data to infer species boundaries: phylogeography of lanceheaded pitvipers in the Brazilian Atlantic forest,and the status of Bothrops pradoi (Squamata: Serpentes: Viperidae)

Phylogeographic studies using mitochondrial DNA sequence information are frequently used as the principal source of evidence to infer species boundaries. However, a critical analysis of further evidence is essential to test whether different haplotype clades identify different species. We demonstrate a hypothesis‐testing approach, using a combination of phylogeographic methods, multivariate morphometrics and matrix association tests, to investigate species boundaries in eastern Brazilian pitvipers conventionally assigned to the species Bothrops leucurus and B. pradoi. Two basal haplotype clades with partly overlapping geographical distributions are identified, which could either represent two partly sympatric species, or multiple haplotypes within one organismal lineage. We use partial Mantel matrix association tests to verify whether generalized morphology, or any of four supposedly diagnostic characters for the two species, show any association with mtDNA variation. Negative results lead to the conclusion that the haplotype clades do not denote independently evolving organismal lineages, and do not constitute separate species under any criterion.     

Hybridization and population genetics of two macaque species in Sulawesi, Indonesia

This study investigates hybridization and population genetics of two species of macaque monkey in Sulawesi, Indonesia, using molecular markers from mitochondrial, autosomal, and Y-chromosome DNA. Hybridization is the interbreeding of individuals from different parental taxa that are distinguishable by one or more heritable characteristics. Because hybridization can affect population structure of the parental taxa, it is an important consideration for conservation management. On the Indonesian island of Sulawesi an explosive diversification of macaques has occurred; seven of 19 species in the genus Macaca live on this island. The contact zone of the subjects of this study, M. maura and M. tonkeana, is located at the base of the southwestern peninsula of Sulawesi. Land conversion in Sulawesi is occurring at an alarming pace; currently two species of Sulawesi macaque, one of which is M. maura, are classified as endangered species. Results of this study indicate that hybridization among M. maura and M. tonkeana has led to different distributions of molecular variation in mitochondrial DNA and nuclear DNA in the contact zone; mitochondrial DNA shows a sharp transition from M. maura to M. tonkeana haplotypes, but nuclear DNA from the parental taxa is homogenized in a narrow hybrid zone. Similarly, within M. maura divergent mitochondrial DNA haplotypes are geographically structured but population subdivision in the nuclear genome is low or absent. In M. tonkeana, mitochondrial DNA haplotypes are geographically structured and a high level of nuclear DNA population subdivision is present in this species. These results are largely consistent with a macaque behavioral paradigm of female philopatry and obligate male dispersal, suggest that introgression between M. maura and M. tonkeana is restricted to the hybrid zone, and delineate one conservation management unit in M. maura and at least two in M. tonkeana.     

Patterns of polymorphism and gene flow of gender-associated mitochondrial DNA lineages in European mussel populations

Mussels of the genus Mytilus have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally and the F type is transmitted maternally. RFLP analysis is used to assess phylogenetic relationships and nucleotide diversity and divergence for both mtDNA genomes in European populations of M. edulis and Atlantic and Mediterranean forms of M. galloprovincialis. Ten restriction endonucleases were used to assay variation in regions of the ND2 and COIII genes for a total of 77 individuals. F and M genomes show a concordant phylogenetic split into two major divergent clades, one specific to Mediterranean M. galloprovincialis and the other containing haplotypes from the three taxa. For both genomes, the geographical distribution of mtDNA variation suggests: (i) extensive levels of mtDNA introgression; (ii) asymmetric mtDNA gene flow from Atlantic to Mediterranean populations; and (iii) recurrent historical hybridization events. Significantly higher mtDNA diversity and divergence are observed for the M than F genome in all three Mytilus taxa, although the evolutionary forces responsible for these differences cannot be resolved. The extensive mtDNA gene flow between European Mytilus taxa conflicts with the restricted mtDNA introgression observed in American mussels , implying geographical variation in the nature of nuclear/mtDNA interactions regulating biparental inheritance.     

Low resolution of mitochondrial COI barcodes for identifying species of the genus Larus (Charadriiformes: Laridae)

We evaluated whether cytochrome c oxidase subunit I (COI) barcodes that have been previously suggested for birds are useful for identifying species of the genus Larus, which are resident or migratory birds in Korea. We found 31 intra- or interspecific COI haplotypes from 12 of 13 Larus species in Korea. Haplotype analyses showed that the COI barcodes could not distinguish some interspecific haplotypes from 6 of 12 Larus species because there were no nucleotide substitutions among their COI haplotypes. The neighbor-joining tree formed shallow branches in the clades expected for L. saundersi, L. crassirostris, L. canus, L relictus, and L. ridbundus. In the nine Larus species, COI haplotypes were not grouped as distinct entities that were correctly assigned to their corresponding species, resulting in polytypic clades. These results indicate that the COI sequences need to be cautiously selected as a DNA barcode for identifying species of Korean Larus birds.     

Polymorphisms distinguishing different mouse species and t haplotypes.

Three anonymous chromosome 17 DNA markers, D17Tu36, D17Tu43, and D17Le66B, differentiate between house mouse species and/or between t chromosomes. The D17Tu36 probe, which maps near the Fu locus and to the In(17)4 on t chromosomes, identifies at least 15 haplotypes, each haplotype characterized by a particular combination of DNA fragments obtained after digestion with the Taq I restriction endonuclease. Ten of these haplotypes occur in Mus domesticus, while the remaining five occur in M. musculus. In each of these two species, one haplotype is borne by t chromosomes while the other haplotypes are present on non-t chromosomes. The D17Tu43 probe, which maps near the D17Leh122 locus and to the In(17)3 on t chromosomes, also identifies at least 15 haplotypes in Taq I DNA digests, of which nine occur in M. domesticus and six in M. musculus. One of the nine M. domesticus haplotypes is borne by t chromosomes, the other haplotypes are borne by non-t chromosomes; two of the six M. musculus haplotypes are borne by t chromosomes and the remaining four by non-t chromosomes. Some of the D17Tu43 haplotypes are widely distributed in a given species, while others appear to be population-specific. Exceptions to species-specificity are found only in a few mice captured near the M. domesticus-M. musculus hybrid zone or in t chromosomes that appear to be of hybrid origin. The D17Leh66B probe, which maps to the In(17)2, distinguishes three haplotypes of M. domesticus-derived t chromosomes and one haplotype of M. musculus-derived t chromosomes. Because of these characteristics, the three markers are well suited for the study of mouse population genetics in general and of t chromosome population genetics in particular. A preliminary survey of wild M. domesticus and M. musculus populations has not uncovered any evidence of widespread introgression of genes from one species to the other; possible minor introgressions were found only in the vicinity of the hybrid zone. Typing of inbred strains has revealed the contribution of only M. domesticus DNA to the chromosome 17 of the laboratory mouse.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

DNA barcoding of the endemic New Zealand leafroller moth genera, Ctenopseustis and Planotortrix

Molecular techniques such as DNA barcoding have become popular in assisting species identification especially for cryptic species complexes. We have analysed data from a 468-bp region of the mitochondrial cytochrome oxidase subunit I (COI) gene from 200 specimens of 12 species of endemic New Zealand leafroller moths (Tortricidae) from the genera Planotortrix and Ctenopseustis to assess whether the DNA barcoding region can distinguish these species. Among the 200 sequences analysed, 72 haplotypes were recovered, with each genus forming a separate major clade. Maximum likelihood phylogenetic methods were used to test whether species fell into reciprocally monophyletic clades. The optimal phylogeny showed that four species within the genus Ctenopseustis (C. obliquana, C. herana, C. filicis and C. fraterna) and three within Planotortrix (P. octo, P. excessana and P. avicenniae) are polyphyletic. Shimodaira-Hasegawa tests rejected a null hypothesis of monophyly for the species C. obliquana, C. herana, P. octo and P. excessana. Comparisons of within and between species levels of sequence divergence for the same set of seven species showed cases where maximum levels of within-species divergence were greater than some levels of between-species divergence. DNA barcoding using this region of the COI gene is able to distinguish the two genera and some species within each genus; however, many species cannot be identified using this method. Finally, we discuss the possible reasons for this polyphyly, including incomplete lineage sorting, introgression, horizontal gene transfer and incorrect taxonomy.