Phylogenetic analysis of <Emphasis Type="Italic">Metarhizium</Emphasis> spp. isolated from soil in Japan

As a result of analyzing the internal transcribed spacer (ITS) and 5′ end of the EF-1α sequence of 145 isolates of Metarhizium spp. isolated from soil in Japan using selective agar medium, eight species were identified. ITS sequence analysis divided the isolates into three clades. Two were identified as M. flavoviride var. pemphigi and M. lepidiotae, respectively. EF-1α sequence analysis identified the other clades as six species: M. anisopliae, M. brunneum, M. guizhouense, M. majus, M. pingshaense and M. robertisii. The distribution of M. flavoviride var. pemphigi was restricted to forest or wood soil, and conidial sizes of M. guizhouense and M. majus were incongruent with the phylogeny based on the sequence of the 5′ end of EF-1α.     

Molecular phylogenetics of the genus <Emphasis Type="Italic">Physalia</Emphasis> (Cnidaria: Siphonophora) in New Zealand coastal waters reveals cryptic diversity

Physalia is a genus of pelagic colonial hydrozoans often known by common names such as ‘Portuguese-man-of-war’ or ‘bluebottle’. Siphonophore systematists generally recognise only a single species in this genus, Physalia physalis, however the name Physalia utriculus is also still in common use, which has led to considerable taxonomic confusion. With some morphological variation between global regions there is the possibility that this genus holds a substantial amount of cryptic variation. We seek to examine the genetic structure of Physalia present in New Zealand coastal waters. Fifty-four specimens collected from 13 locations around New Zealand and Australia were sequenced for both mitochondrial cytochrome c oxidase I (COI) and the first internal transcribed spacer (ITS1) of the nuclear ribosomal cistron. Sequences were analysed using maximum likelihood and split decomposition neighbour networks to determine conflict between clans (the unrooted analog of clades). Three clans were identified from both the COI and ITS sequences. The results are complex and clans are not consistent between the two genes. Nevertheless, it seems that there is substantial cryptic diversity amongst Physalia present in New Zealand coastal waters.     

Morphological and phylogenetic analyses of <Emphasis Type="Italic">Uromyces appendiculatus</Emphasis> and <Emphasis Type="Italic">U. vignae</Emphasis> on legumes in Japan

Uromyces appendiculatus, inclusive of three varieties, is distinguished from U. vignae primarily by the position of urediniospore germ pores and putative host specificity. However, opinions concerning these morphological and physiological features as taxonomic characters have varied greatly, and distinction of these species has often been confused. To clarify the taxonomy of these two species, morphological features of urediniospores and teliospores of 225 rust fungus specimens on species of Phaseolus, Vigna, Apios, Lablab, and Dunbaria were examined by light microscopy and scanning electron microscopy. Forty-five specimens were subjected to molecular phylogenetic analyses. As a result, the position of germ pores in urediniospores and the teliospore-wall thickness were considered as good characters to separate three morphological groups. In molecular analyses, the specimens fell into two and three clades based on the nucleotide sequence at D1/D2 domain of LSU rDNA and ITS regions, respectively. One of the D1/D2 clades corresponded to one morphological group whereas another D1/D2 clade included two other morphological groups. In contrast, each of the three ITS clades corresponded to a separate morphological group. Neither morphological groups nor molecular clades were host limited. It is suggested that the three morphological groups that corresponded to three distinct ITS clades constitute distinct species.Contribution no. 186 from the Laboratory of Plant Parasitic Mycology, Institute of Agriculture and Forestry, University of Tsukuba, Japan     

The phylogenetic position of an <Emphasis Type="Italic">Armillaria</Emphasis> species from Amami-Oshima,a subtropical island of Japan,based on elongation factor and ITS sequences

An undetermined Armillaria species was collected on Amami-Oshima, a subtropical island of Japan. The phylogenetic position of the Armillaria sp. was determined using sequences of the elongation factor-1α (EF-1α) gene and the internal transcribed spacer (ITS) region (ITS1-5.8S-ITS2) of ribosomal DNA (rDNA). The phylogenetic analyses based on EF-1α and ITS sequences showed that this species differs from known Japanese taxa of Armillaria. The sequences of this species and A. novae-zelandiae from Southeast Asia were contained in a strongly supported clade, which was adjacent to a well-supported sister clade containing A. novae-zelandiae from Australia and New Zealand.     

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.     

Phylogenomic species delimitation and host‐symbiont coevolution in the fungus‐farming ant genus Sericomyrmex Mayr (Hymenoptera: Formicidae): ultraconserved elements (UCEs) resolve a recent radiation

Ants in the Neotropical genus Sericomyrmex Mayr cultivate fungi for food. Both ants and fungi are obligate, coevolved symbionts. The taxonomy of Sericomyrmex is problematic because the morphology of the worker caste is generally homogeneous across all of the species within the genus, species limits are vague, and the relationships between them are unknown. We used ultraconserved elements (UCEs) as genome‐scale markers to reconstruct evolutionary history and to infer species boundaries in Sericomyrmex. We recovered an average of ～990 UCE loci for 88 Sericomyrmex samples from across the geographical range of the genus as well as for five outgroup taxa. Using maximum likelihood and species‐tree approaches, we recovered nearly identical topologies across datasets with 50–95% matrix completeness. We identify nine species‐level lineages in Sericomyrmex, including two new species. This is less than the previously described 19 species, even accounting for two species for which we had no UCE samples, which brings the total number of Sericomyrmex species to 11. Divergence‐dating analyses recovered 4.3 Ma as the crown‐group age estimates for Sericomyrmex, indicating a recent, rapid radiation. We also sequenced mitochondrial cytochrome oxidase subunit I (COI) for 125 specimens. Resolution and support for clades in our COI phylogeny are weak, indicating that COI is not an appropriate species‐delimitation tool. However, taxa within species consistently cluster together, suggesting that COI is useful as a species identification (‘DNA barcoding’) tool. We also sequenced internal transcribed spacer (ITS) and large subunit (LSU) for 32 Sericomyrmex fungal cultivars. The fungal phylogeny confirms that Sericomyrmex fungi are generalized higher‐attine cultivars, interspersed with Trachymyrmex‐associated fungal species, indicating cultivar sharing and horizontal transfer between these two genera. Our results indicate that UCEs offer immense potential for delimiting and resolving relationships of problematic, recently diverged species.     

Adhesion forces in the self-recognition of oligosaccharide epitopes of the proteoglycan aggregation factor of the marine sponge <Emphasis Type="Italic">Microciona prolifera</Emphasis>

Cell aggregation in the marine sponge Microciona prolifera is mediated by a multimillion molecular-mass aggregation factor, termed MAF. Earlier investigations revealed that the cell aggregation activity of MAF depends on two functional domains: (i) a Ca²⁺-independent cell-binding domain and (ii) a Ca²⁺-dependent proteoglycan self-interaction domain. Structural analysis of involved carbohydrate fragments of the proteoglycan in the self-association established a sulfated disaccharide β-d-GlcpNAc3S-(1→3)-α-l-Fucp and a pyruvated trisaccharide β-d-Galp4,6(R)Pyr-(1→4)-β-d-GlcpNAc-(1→3)-α-l-Fucp. Recent UV, SPR, and TEM studies, using BSA conjugates and gold nanoparticles of the synthetic sulfated disaccharide, clearly demonstrated self-recognition on the disaccharide level in the presence of Ca²⁺-ions. To determine binding forces of the carbohydrate–carbohydrate interactions for both synthetic MAF oligosaccharides, atomic force microscopy (AFM) studies were carried out. It turned out that, in the presence of Ca²⁺-ions, the force required to separate the tip and sample coated with a self-assembling monolayer of thiol-spacer-containing β-d-GlcpNAc-(1→3)-α-l-Fucp-(1→O)(CH₂)₃S(CH₂)₆S- was found to be quantized in integer multiples of 30 ± 6 pN. No binding was observed between the two monolayers in the absence of Ca²⁺-ions. Cd²⁺-ions could partially induce the self-interaction. In contrast, similar AFM experiments with thiol-spacer-containing β-d-Galp4,6(R)Pyr-(1→4)-β-d-GlcpNAc-(1→3)-α-l-Fucp-(1→O)(CH₂)₃S(CH₂)₆S- did not show a binding in the presence of Ca²⁺-ions. Also TEM experiments of gold nanoparticles coated with the pyruvated trisaccharide could not make visible aggregation in the presence of Ca²⁺-ions. It is suggested that the self-interaction between the sulfated disaccharide fragments is stronger than that between the pyruvated trisaccharide.     

Characterization of the <Emphasis Type="Italic">Pragia fontium</Emphasis> Lipopolysaccharides

Lipopolysaccharides (LPSs) of two strains Pragia fontium 97U116 and 27480 were isolated and characterized; they were close to those of other representatives of the family Enterobacteriaceae in fatty acid composition and contained, respectively, 3-hydroxytetradecanoic acid as the predominant component (45.8 and 45.1%), tetradecanoic (23.5 and 28.9%), hexadecanoic (12.6 and 7.9%), hexadecenoic (12.6 and 7.9%), and dodecanoic (4.9 and 4.2%) fatty acids. The O-specific polysaccharides consisted of linear penta- and tetrasaccharide repeating units: →2)-α-D-Galf-(1→3)-α-L-Rhap2Ac-(1→4)-α-D-GlcpNAc-(1→2)-α-L-Rhap-(1→3)-β-D-GlcpNAc-(1→ →4)-β-D-ManpNAc3NAcA-(1→2)-α-L-Rhap-(1→3)-β-L-Rhap-(1→4)-α-D-GlcpNAc-(1→ The LPSs of P. fontium 97U116 and 27480 were serologically active and belonged to different serogroups; they were less toxic than those of strain E. coli O55:B5, but more pyrogenic than the Pyrogenal preparation.     

Genetic relationships between Oeneis urda and O. mongolica (Nymphalidae: Lepidoptera)

The species status of Oeneis urda (Eversmann) and O. mongolica (Oberthür) has been argued based on morphological characters. Reexamination of their major morphological characters has shown a slight differentiation in the two species. Sequences of three mitochondrial genes (COI, ND6, and ND1) and one nuclear region (internal transcribed spacer 2, ITS2) from two O. urda populations (Yangyang and Mt. Hanla) and one O. mongolica population (Uljin) were performed for phylogenetic and population genetic inferences. Sharing of identical sequences in the ND6 gene and ITS2, minimal sequence divergence in the COI and ND1 genes, and phylogenetically undividable sequence types in all mitochondrial genes and ITS2 suggest genetic continuity between the two species. Nevertheless, significant F_ST estimates (P < 0.05) were found for the COI gene in comparisons between Yangyang (O. urda) and Uljin (O. mongolica), between Yangyang (O. urda) and Mt. Hanla (O. urda), and between Uljin (O. mongolica) and Mt. Hanla (O. urda) populations. These F_ST estimates, along with other gene‐based analyses collectively suggest isolation of the two species at some point in the past, but incomplete separation between the two species on the mainland (Yangyang and Uljin) and biogeographically forced isolation of the O. urda population on Mt. Hanla collectively appear to complicate species status of these two species that were once further clearly separated.     

Asymmetric synthesis of (<Emphasis Type="Italic">R</Emphasis>)-3-fluoroalanine from 3-fluoropyruvate using omega-transaminase

In this study, (R)-3-fluoroalanine was asymmetrically synthesized from 3-fluoropyruvate (F-pyruvate) and (S)-α-methylbenzylamine (MBA) using recombinant ω-transaminase (TA) from Vibrio fluvialis JS17. The reaction was severely inhibited by acetophenone (deaminated product of α-MBA). In the presence of 5 mM acetophenone, the reactivity of the enzyme towards F-pyruvate decreased by 78%. To overcome the product inhibition by acetophenone, a biphasic reaction was successfully used. The conversion of F-pyruvate into (R)-3-fluoroalanine (enatiomeric exess (e.e.) > 99%) was about 95% in the biphasic system (75 mM F-pyruvate, 100 mM (S)-α-MBA, and 3.0 U/mL), whereas 31% was obtained without product extraction. The use of racemic α-MBA as an amino donor instead of (S)-α-MBA can reduce the reaction cost and also produce chiral amines through kinetic resolution. When the kinetic resolution of racemic α-MBA (40 mM) was carried out with F-pyruvate (30 mM) and ω-TA (3.0 U/mL) in 100 mM phosphate buffer (pH 7.0), the e.e. of (R)-α-MBA reached 98.4% with 52.2% conversion for 10 h and 21 mM (R)-3-fluoroalanine was produced with 70% conversion and an e.e. > 99%.     

The distribution and molecular diversity of the Eastern Atlantic and Mediterranean chthamalids (Crustacea, Cirripedia)

The three chthamalids Chthamalus stellatus , C. montagui and Euraphia depressa are common inhabitants of the intertidal zone in the Eastern Atlantic, Mediterranean Sea and Black Sea. In this study, we investigated the occurrence of these barnacles in a wide range of their distribution. Population divergences of these two species have been inferred using three molecular markers — internal transcribed spacer (ITS), elongation factor 1α (EF-1α) and cytochrome oxidase subunit I (COI). ITS sequences of C. stellatus were identical throughout the species range, whereas ITS sequences of C. montagui indicated that the Black Sea and Mediterranean populations are isolated from the Atlantic population. The COI and EF-1α sequences were the most variable and informative. They indicated a high genetic divergence between Atlantic, Mediterranean and Black Sea populations for C. montagui . In addition significant genetic structure was found among the populations of C. stellatus based on EF-1α but not COI. Interestingly, our molecular dating analysis correlated the pattern of diversification in C. montagui to major geological changes that occurred in the Mediterranean during the end of the Messinian and Pleiocene periods. We suggest that palaeohistory shaped the divergences between Chthamalus populations that have probably been maintained by current hydrographic conditions. Finally, COI phylogenetic analysis placed the genus Euraphia within the Chthamalus clade, suggesting the need for a taxonomic revision of Euraphia . This study represents the most detailed phylogeographical analysis of intertidal Mediterranean species to date, and shows that geological events have strongly shaped the current diversity pattern of this fauna.     

Relationships within aphids Cinara (Cupressobium) (Hemiptera) based on mitochondrial and nuclear DNA sequences

The relationships between Cinara (Cupressobium) aphids inhabiting woody parts and leaves of conifers belonging to Cupressaceae have been studied using a mitochondrial gene (COI) and a nuclear gene (EF1-α). Based on the COI sequences, genetic distances between species ranged from 5.6 % between Cinara (C.) tujafilina (del Guercio) and Cinara (C.) juniperi (De Geer) to 10.5 % between C. (C.) tujafilina and Cinara (C.) mordvilkoi (Pa?ek). Genetic distances among EF1-α sequences were lower and showed from 0.1 % between C. cupressi and C. juniperi to 2.3 % between C. tujafilina and C. mordvilkoi. Molecular phylogenetic trees were constructed using the Bayesian inference (BI) phylogenetic analysis and maximum parsimony (MP) criterion. Phylogenetic trees obtained based on COI and EF1-α marker genes created two sister clades. Our results indicate that Cinara (Cupressobium) are a monophyletic group of aphids. Phylogenetic relationships amongst Cupressobium aphids do not result from the association with the host plant, but from the feeding site on the host plant or an ability to change the microhabitat on the plant. As closely related species inhabit similar microhabitats on different host plants, it suggests that the host switching is the main mode of speciation in this subgenus.     

Conservation genetics of highly isolated populations of the xerothermic beetle Crioceris quatuordecimpunctata (Chrysomelidae)

Xerothermic species are rare and threatened in central and eastern Europe. In light of the continuing loss of steppe‐like habitats due to anthropogenic fragmentation and degradation, the evaluation of genetic variation in populations inhabiting them is of immediate importance if appropriate conservation measures are to be undertaken. Here we report on the genetic diversity of the rare leaf beetle Crioceris quatuordecimpunctata, whose populations in central and eastern Europe inhabit highly geographically isolated areas. All of the studied populations (in Poland, Ukraine, and Slovakia) were differentiated at the mitochondrial marker COI. However, with respect to the nuclear marker ITS1, Polish populations were monomorphic, but distinct from all other populations. The distinctiveness of the studied populations was confirmed by Wolbachia screening, which showed that all populations carried different strains (one or two), which were probably transferred independently from other insects. On the other hand, no diversity was found in any marker within particular populations, which could be caused (at least for mtDNA) by a Wolbachia selective sweep. Crioceris quatuordecimpunctata probably consists of isolated populations, which went through narrow bottlenecks leading to a drastic reduction in their genetic diversity. As these populations are reciprocally monophyletic for mtDNA haplotypes and show a significant divergence of allele frequencies at nuclear loci, they could be classified as evolutionarily significant units (ESUs). In addition, DNA barcodes were used to identify Asparagus officinalis as the host plant for members of all studied populations. These data should be valuable in efforts to conserve populations of C. quatuordecimpunctata (e.g., for guiding reintroductions).     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

The phylogeny of <Emphasis Type="Italic">Schistidium</Emphasis> (Bryophyta,Grimmiaceae) based on the primary and secondary structure of nuclear rDNA internal transcribed spacers

Phylogeny of Schistidium (Bryophyta, Grimmiaceae) was studied by comparing the nucleotide sequences of internal transcribed spacers ITS1-2 of nuclear rDNA and the trnT-trnD region of chloroplast DNA. Phylogenetic trees constructed based on nuclear and chloroplast sequences were consistent, comprising a basal grade and two large clades. Morphological characteristics specific for these clades were described. Secondary structures of ITS1 and ITS2 Schistidium species were modeled using thermodynamic criteria. Four different structures of the longest ITS1 hairpin were identified. These results were used to analyze possible paths of Schistidium evolution. Characteristics of the ITS2 secondary structure support the two major clades recognized in the phylogenetic trees.     

Diversity in the genus Rhabdias (Nematoda,Rhabdiasidae): Evidence for cryptic speciation

Lungworms from the genus Rhabdias are common parasites of amphibians and reptiles distributed worldwide. To assess the diversity of Rhabdias spp., we performed molecular analyses of 35 specimens sampled in different regions of Brazil. Molecular analyses were based on the internal transcribed spacer (ITS), large subunit (28S) ribosomal and the cytochrome oxidase I (COI) mitochondrial genes. DNA sequence divergence was compared among ribosomal and mitochondrial genes, analyses using the general mixed Yule‐coalescent (GMYC) method based on the COI gene were used to identify possible cryptic diversity, and phylogenetic analyses using concatenated ITS and 28S ribosomal genes were used to test the monophyly of Rhabdiasidae. We revealed five morphospecies: R. cf. stenocephala, R. breviensis, R. pseudosphaerocephala and two new species, Rhabdias sp.4 and Rhabdias sp.5. DNA sequence levels of divergence among genes ITS, 28S and COI were compared, and the efficiency of the molecular markers to identify species (ITS and COI) and lineages (COI) was tested. GMYC was assigned to 17 well‐supported clades (i.e., 17 species), and cryptic diversity was detected in the Neotropical region as evidenced by the multiple lineages in R. breviensis and R. pseudosphaerocephala. In addition, our results suggest evidence for host–parasite cophylogeny in the R. pseudosphaerocephala complex and dispersal events among their populations. Phylogenetic analyses supported the monophyly of Rhabdiasidae and improved the resolution of main clades. Rhabdias breviensis is closely related to Rhabdias cf. africanus, Rhabdias cf. stenocephala, R. pseudosphaerocephala, Rhabdias sp.4 and Rhabdias sp.5 grouping together in a main clade with Neotropical‐related species. The large geographical distribution appeared to be a phylogenetic pattern among the species of Rhabdias from the neotropics.     

Structural peculiarities of the O-specific polysaccharides of <Emphasis Type="Italic">Azospirillum</Emphasis> bacteria of serogroup III

Lipopolysaccharides and O-specific polysaccharides were isolated from the outer membrane of bacterial cells of three strains belonging to two Azospirillum species, and their structures were established by monosaccharide analysis including determination of the absolute configurations, methylation analysis, and one- and two-dimensional NMR spectroscopy. It was shown that while having the identical composition, the O-polysaccharides have different branched tetrasaccharide repeating units. Two neutral polysaccharides were found in the lipopolysaccharide of A. brasilense 54, and the structure for the predominant O-polysaccharide was determined. The structural data, together with results of serological studies, enabled assignment of strains examined to a novel serogroup, III. The chemical basis for the serological relatedness among the azospirilla of this serogroup is presumably the presence of a common →3)-α-L-Rhap-(1→2)-α-L-Rhap-(1→3)-α-L-Rhap-(1→oligosaccharide motif in their O-polysaccharides.     

Molecular phylogeny of the subterranean genus Niphargus (Crustacea: Amphipoda) in the Middle East: a comparison with European Niphargids

The subterranean genus Niphargus is one of the most species‐rich genera among freshwater amphipods in the world, distributed in the Western Palearctic. Thus far, taxonomic and phylogenetic research has focused mainly on the European half of the genus range. In this study, 25 populations of Niphargus from Iran, Lebanon and the Crimean Peninsula were investigated. Bayesian inference based on 28S, H3 and COI gene sequences suggests that populations from the area belong to four different clades. Three species from Crimea and one species from Iran are nested at basal nodes, indicating their rather ancient origin. The rest of the species are younger and belong to two separate clades. One Crimean species is a sister‐species to east Romanian species. The second clade includes one species from Lebanon and all but one population from Iran. The origin of this clade corresponds to marine transgression between the Black Sea and Mediterranean approximately 12 Mya. This clade was further investigated taxonomically. Revision of qualitative morphological traits and unilocus species delimitation methods using COI suggest that this clade comprises 12–16 species, of which only three have been described so far. Multilocus coalescence delimitation methods (using fragments of COI, 28S, H3 and ITS) strongly supported 11 of these species. The remaining populations comprise at least two species complexes that require further and more detailed taxonomic research. © 2015 The Linnean Society of London     

Phylogeographic history of the New Zealand stick insect Niveaphasma annulata (Phasmatodea) estimated from mitochondrial and nuclear loci

We have assessed the utility of a single-copy nuclear locus and mitochondrial DNA (mtDNA) in a phylogeographic study of the New Zealand stick insect Niveaphasma annulata (Hutton). We amplified sequences from the mitochondrial cytochrome oxidase subunit I (COI) gene and the single-copy nuclear gene elongation factor-1α (EF1α) from 97 individuals. Allelic phase at the EF1α locus was determined using Denaturing Gradient Gel Electrophoresis. Phylogenetic analyses showed broad congruence between the geographic distribution of three major COI clades and EF1α alleles, which suggested that the phylogenetic patterns reflect population history rather than lineage sorting. However, the geographic boundaries of these clades were not always in exact agreement between the two loci. Our data indicate that Niveaphasma annulata was most likely separated into a number of refugia during Pleistocene glacial advances. Subsequent to glacial retreat these refugial populations have expanded and now form a number of zones of secondary contact. We contrast these patterns with those observed from other New Zealand taxa. Our study offers compelling evidence for the use of nuclear genes alongside mtDNA for future phylogeographic studies.