Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Phylogeography and population history of the least weasel (Mustela nivalis) in the Palearctic based on multilocus analysis

The least weasel (Mustela nivalis) is one of the most widely distributed carnivorans. While previous studies have identified distinct western and eastern mitochondrial DNA (mtDNA) lineages of the species in the western Palearctic, their broader distributions across the Palearctic have remained unknown. To address the broad-scale phylogeographical structure, we expanded the sampling to populations in Eastern Europe, the Urals, the Russian Far East, and Japan, and analyzed the mtDNA control region and cytochrome b, the final intron of the zinc finger protein on Y chromosome (ZFY), and the autosomal agouti signaling protein gene (ASIP). The mtDNA data analysis exposed the previous western lineage (Clade I) but poorly supported assemblage extending across Palearctic, whereas the previous eastern lineage (Clade II) was reconfirmed and limited in the south western part of the Palearctic. The ZFY phylogeny showed a distinctive split that corresponding to the mtDNA lineage split, although less phylogeographical structure was seen in the ASIP variation. Our data concur with the previous inference of the Black Sea–Caspian Sea area having an ancestral character. The Urals region harbored high mitochondrial diversity, with an estimated coalescent time of around 100,000 years, suggesting this could have been a cryptic refugium. Based on the coalescent-based demographic reconstructions, the expansion of Clade I across the Palearctic was remarkably rapid, while Clade II was relatively stable for a longer time. It seems that Clade II has maintained a constant population size in the temperate region, and the expansive Clade I represents adaptation to the cold regions.     

Nuclear markers confirm taxonomic status and relationships among highly endangered and closely related right whale species

Right whales (genus: Eubalaena) are among the most endangered mammals, yet their taxonomy and phylogeny have been questioned. A phylogenetic hypothesis based on mitochondrial DNA (mtDNA) variation recently prompted a taxonomic revision, increasing the number of right whale species to three. We critically evaluated this hypothesis using sequence data from 13 nuclear DNA (nuDNA) loci as well as the mtDNA control region. Fixed diagnostic characters among the nuclear markers strongly support the hypothesis of three genetically distinct species, despite lack of any diagnostic morphological characters. A phylogenetics analysis of all data produced a strict consensus cladogram with strong support at nodes that define each right whale species as well as relationships among species. Results showed very little conflict among the individual partitions as well as congruence between the mtDNA and nuDNA datasets. These data clearly demonstrate the strength of using numerous independent genetic markers during a phylogenetics analysis of closely related species. In evaluating phylogenetic support contributed by individual loci, 11 of the 14 loci provided support for at least one of the nodes of interest to this study. Only a single marker (mtDNA control region) provided support at all four nodes. A study using any single nuclear marker would have failed to support the proposed phylogeny, and a strong phylogenetic hypothesis was only revealed by the simultaneous analysis of many nuclear loci. In addition, nu DNA and mtDNA data provided complementary levels of support at nodes of different evolutionary depth indicating that the combined use of mtDNA and nuDNA data is both practical and desirable.     

Molecular phylogeny and divergence time of the water penny genus Eubrianax (Coleoptera: Psephenidae) in Japan

The phylogenetic relationships among the Japanese members of the genus Eubrianax (Coleoptera: Psephenidae) were examined using the mitochondrial cytochrome oxidase subunit I (COI) gene and nuclear 28S rRNA gene sequences. Based on the molecular phylogeny as well as morphological features, the species status of Eubrianax brunneicornis Nakane, 1952 was proposed. The phylogenetic analyses recovered monophyly of the previously proposed pellucidus species group with four Japanese species, whereas a single Japanese species of the granicollis group was included in the lineage of the ramicornis group with five Japanese species. The divergence times of the species were estimated by dating the phylogenetic tree against the fossil record and a molecular clock based on the COI gene. The divergence of the Japanese species was inferred to have occurred during the Pliocene epoch.     

Molecular phylogeographic analysis of East Asian cryptoperlan stoneflies (Insecta: Plecoptera, Peltoperlidae)

Cryptoperlan stoneflies inhabit the headwaters or upper stream areas of rivers. A molecular phylogeographic study of cryptoperlans in the Japanese archipelago and on Taiwan Island has been conducted. Altogether the mtDNA 16S rRNA region of 71 individuals from 61 populations, the mtDNA COI region of 76 individuals from 41 populations, and the nDNA Histone 3 region of 56 individuals from 52 populations were sequenced and analyzed. The respective ML, NJ, MP and Bayesian dendrograms were proposed from the sequencing data for the 16S rRNA region (362-bp), the COI region (540-bp), and the Histone 3 region (322-bp), estimated using Yoraperla uenoi as an outgroup. Based upon those data and the resulting dendrograms, it has become clear that the cryptoperlan stoneflies of the Japanese archipelago and those of Taiwan Island comprise two major clades. The first of these two major clades consists of a number of OTUs [operational taxonomic units: Cryptoperla japonica (Honshu, Shikoku and Kyushu Islands) + C. ishigakiensis (Ishigaki-jima Island) + Cryptoperla spp. (Okinawa-jima and Taiwan Islands)]. The other clade consists of the species Cryptoperla kawasawai inhabiting only Shikoku Island. Of particular note, C. kawasawai was observed to be significantly genetically differentiated from all other cryptoperlans examined. Yet, despite the fact that the specimens of C. japonica were taken from a very broad range of populations, their genetic diversity was relatively low, similar to that of C. kawasawai, which inhabits only a limited region within Shikoku Island. Furthermore, even the species C. kawasawai was revealed to be composed of two significantly genetically differentiated subclades. It is considered that this genetic structure among cryptoperlans largely reflects the geological history from the middle to upper Miocene Epoch (i.e., Tortonian stage) of the Japanese archipelago and Taiwan Island.     

Analysis of mitochondrial COI sequences of the Diachasmimorpha longicaudata (Hymenoptera: Braconidae) species complex in Thailand

Diachasmimorpha longicaudata is an Opiinae parasitoid used to control tephritid fruit flies, which cause tremendous economic losses of fruits worldwide. In Thailand, D. longicaudata is classified as three sibling species, DLA, DLB and DLBB, based on the morphological and biological species concepts but their genetic variation has not been studied. Therefore, we investigated the genetic differentiation of the mitochondrial COI gene to clarify the ambiguous taxonomy of this species complex. The 603‐bp COI region was sequenced from laboratory‐bred colonies and field‐collected specimens from seven locations representing five geographical regions in Thailand. DLA was associated with the host Bactrocera correcta while DLB and DLBB were associated with Bactrocera dorsalis. The interspecific nucleotide differences of COI sequences among the three groups ranged from 6.70% to 7.62% (Kimura 2‐parameter distance), which adequately separates species complexes within the order Hymenoptera and supports the current sibling species classification. The neighbor joining, maximum likelihood and consensus Bayesian phylogenetic trees constructed from COI sequences revealed that the three sibling species of laboratory and field‐collected D. longicaudata are monophyletic with 100% support. The high genetic variation and molecular phylogeny of the COI sequences were shown to discriminate between the D. longicaudata species examined in this study.     

Genetic admixture supports an ancient hybrid origin of the endangered Hawaiian duck

Speciation is regarded primarily as a bifurcation from an ancestral species into two distinct taxonomic units, but gene flow can create complex signals of phylogenetic relationships, especially among different loci. We evaluated several hypotheses that could account for phylogenetic discord between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA) within Hawaiian duck (Anas wyvilliana), including stochastic lineage sorting, mtDNA capture and widespread genomic introgression. Our results best support the hypothesis that the contemporary Hawaiian duck is descended from an ancient hybridization event between the mallard (Anas platyrhynchos) and Laysan duck (Anas laysanensis). Whereas mtDNA clearly shows a sister relationship between Hawaiian duck and mallard, nuDNA is consistent with a genetic mosaic with nearly equal contributions from Laysan duck and mallard. In addition, coalescent analyses suggest that gene flow from either mallard or Laysan duck, depending on the predefined tree topology, is necessary to explain contemporary genetic diversity in Hawaiian ducks, and these estimates are more consistent with ancient, rather than contemporary, hybridization. Time since divergence estimates suggest that the genetic admixture event occurred around the Pleistocene–Holocene boundary, which is further supported by circumstantial evidence from the Hawaiian subfossil record. Although the extent of reproductive isolation from either putative parental taxon is not currently known, these species are phenotypically, genetically and ecologically different, and they meet primary criteria used in avian taxonomy for species designation. Thus, the available data are consistent with an admixed origin and support the hypothesis that the Hawaiian duck may represent a young hybrid species.     

Behavioural vs. molecular sources of conflict between nuclear and mitochondrial DNA: the role of male-biased dispersal in a Holarctic sea duck

Genetic studies of waterfowl (Anatidae) have observed the full spectrum of mitochondrial (mt) DNA population divergence, from apparent panmixia to deep, reciprocally monophyletic lineages. Yet, these studies often found weak or no nuclear (nu) DNA structure, which was often attributed to male‐biased gene flow, a common behaviour within this family. An alternative explanation for this ‘conflict’ is that the smaller effective population size and faster sorting rate of mtDNA relative to nuDNA lead to different signals of population structure. We tested these alternatives by sequencing 12 nuDNA introns for a Holarctic pair of waterfowl subspecies, the European goosander (Mergus merganser merganser) and the North American common merganser (M. m. americanus), which exhibit strong population structure in mtDNA. We inferred effective population sizes, gene flow and divergence times from published mtDNA sequences and simulated expected differentiation for nuDNA based on those histories. Between Europe and North America, nuDNA Ф_ST was 3.4‐fold lower than mtDNA Ф_ST, a result consistent with differences in sorting rates. However, despite geographically structured and monophyletic mtDNA lineages within continents, nuDNA Ф_ST values were generally zero and significantly lower than predicted. This between‐ and within‐continent contrast held when comparing mtDNA and nuDNA among published studies of ducks. Thus, male‐mediated gene flow is a better explanation than slower sorting rates for limited nuDNA differentiation within continents, which is also supported by nonmolecular data. This study illustrates the value of quantitatively testing discrepancies between mtDNA and nuDNA to reject the null hypothesis that conflict simply reflects different sorting rates.     

Nuclear and mitochondrial ribosomal RNA variability in the obscura group of Drosophila

Parts of 28S (nuclear) and 12S (mitochondrial) ribosomal RNA of Palearctic, Nearctic and African species of the obscura group have been sequenced by the direct method of sequencing. Rates of nucleotide substitutions in both molecules were compared. The nucleotide divergence is higher in the mitochondrial rRNA. Average distances of species taken in pairwise were compared to results obtained with the melanogaster subgroup: the divergence of nuclear rRNA appears lower, that of the mtDNA higher whereas genetic distances (allozymes) and sncDNA distances are similar. Noticeable variability of evolutionary rates can be observed even in low taxonomical levels. Phylogenetic trees for the obscura group are in general agreement with those obtained with other characters.     

Nuclear microsatellite and mitochondrial DNA analyses reveal the regional genetic structure and phylogeographical history of a sanguivorous land leech,Haemadipsa japonica,in Japan

Recent molecular studies have indicated that phylogeographical history of Japanese biota is likely shaped by geohistory along with biological events, such as distribution shifts, isolation, and divergence of populations. However, the genetic structure and phylogeographical history of terrestrial Annelida species, including leech species, are poorly understood. Therefore, we aimed to understand the genetic structure and phylogeographical history across the natural range of Haemadipsa japonica, a sanguivorous land leech species endemic to Japan, by using nine polymorphic nuclear microsatellites (nSSR) and cytochrome oxidase subunit one (COI) sequences of mitochondrial DNA (mtDNA). Analyses using nSSR revealed that H. japonica exhibited a stronger regional genetic differentiation among populations (G'_ST = 0.77) than other animal species, probably because of the low mobility of land leech. Analyses using mtDNA indicated that H. japonica exhibited two distinct lineages (A and B), which were estimated to have diverged in the middle Pleistocene and probably because of range fragmentation resulting from climatic change and glacial and interglacial cycles. Lineage A was widely distributed across Japan, and lineage B was found in southwestern Japan. Analyses using nSSR revealed that lineage A was roughly divided into two population groups (i.e., northeastern and southwestern Japan); these analyses also revealed a gradual decrease in genetic diversity with increasing latitude in lineage A and a strong genetic drift in populations of northeastern Japan. Combined with the largely unresolved shallow polytomies from the mtDNA phylogeny, these results implied that lineage A may have undergone a rapid northward migration, probably during the Holocene. Then, the regional genetic structure with local unique gene pools may have been formed within each lineage because of the low mobility of this leech species.     

Multilocus phylogeography of Australian teals (Anas spp.): a case study of the relationship between vagility and genetic structure

Biogeographic barriers potentially restrict gene flow but variation in dispersal or vagility can influence the effectiveness of these barriers among different species and produce characteristic patterns of population genetic structure. The objective of this study was to investigate interspecific and intraspecific genetic structure in two closely related species that differ in several life‐history characteristics. The grey teal Anas gracilis is geographically widespread throughout Australia with a distribution that crosses several recognized biogeographic barriers. This species has high vagility as its extensive movements track broad‐scale patterns in rainfall. In contrast, the closely related chestnut teal A. castanea is endemic to the mesic southeastern and southwestern regions of Australia and is more sedentary. We hypothesized that these differences in life‐history characteristics would result in more pronounced population structuring in the chestnut teal. We sequenced five nuclear loci (nuDNA) for 49 grey teal and 23 chestnut teal and compared results to published mitochondrial DNA (mtDNA) sequences. We used analysis of molecular variance to examine population structure, and applied coalescent based approaches to estimate demographic parameters. As predicted, chestnut teal were more strongly structured at both mtDNA and nuDNA (Φ_ST= 0.163 and 0.054, respectively) than were grey teal (Φ_ST < 0.0001 for both sets of loci). Surprisingly, a greater proportion of the total genetic variation was partitioned among populations within species (Φ_SC= 0.014 and 0.047 for nuDNA and mtDNA, respectively) than between the two species (Φ_CT < 0.0001 for both loci). The ‘Isolation with Migration’ coalescent model suggested a late Pleistocene divergence between the taxa, but remarkably, a deeper divergence between the southeastern and southwestern populations of chestnut teal. We conclude that dispersal potential played a prominent role in the structuring of populations within these species and that divergent selection associated with ecology and life history traits likely contributed to rapid and recent speciation in this pair.     

Morphology versus DNA barcoding: two sides of the same coin. A case study of Ceutorhynchus erysimi and C. contractus identification

Genotyping of 2 well‐known weevil species from the genus Ceutorhynchus (Coleoptera: Curculionidae) distributed in west Palearctic, C. erysimi and C. contractus, revealed phenotype versus genotype inconsistencies in a set of 56 specimens (25 C. erysimi and 31 C. contractus) collected from 25 locations in Serbia and Montenegro. An analysis of the mitochondrial cytochrome oxidase subunit I gene (COI), widely used as a barcoding region, and a nuclear gene, elongation factor‐1α (EF‐1α), revealed stable genetic divergence among these species. The average uncorrected pairwise distances for the COI and EF‐1α genes were 3.8%, and 1.3%, respectively, indicating 2 genetically well‐segregated species. However, the genetic data were not congruent with the phenotypic characteristics of the studied specimens. In the first place, C. erysimi genotypes were attached to specimens with phenotypic characteristics of C. contractus. Species‐specific PCR‐RFLP assays for the barcoding gene COI were applied for the molecular identification of 101 additional specimens of both morphospecies (33 C. erysimi and 68 C. contractus) and were found to confirm this incongruity. The discrepancy between the genetic and morphological data raises the question of the accuracy of using a barcoding approach, as it may result in misleading conclusions about the taxonomic position of the studied organism. Additionally, the typological species concept shows considerable weakness when genetic data are not supported with phenotypic characteristics as in case of asymmetric introgression, which may cause certain problems, especially in applied studies such as biological control programs in which the biological properties of the studied organisms are the main focus.     

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.     

Toward the resolution of an explosive radiation--a multilocus phylogeny of oceanic dolphins (Delphinidae)

Oceanic dolphins (Delphinidae) are the product of a rapid radiation that yielded ～36 extant species of small to medium-sized cetaceans that first emerged in the Late Miocene. Although they are a charismatic group of organisms that have become poster children for marine conservation, many phylogenetic relationships within Delphinidae remain elusive due to the slow molecular evolution of the group and the difficulty of resolving short branches from successive cladogenic events. Here I combine existing and newly generated sequences from four mitochondrial (mt) genes and 20 nuclear (nu) genes to reconstruct a well-supported phylogenetic hypothesis for Delphinidae. This study compares maximum-likelihood and Bayesian inference methods of several data sets including mtDNA, combined nuDNA, gene trees of individual nuDNA loci, and concatenated mtDNA+nuDNA. In addition, I contrast these standard phylogenetic analyses with the species tree reconstruction method of Bayesian concordance analysis (BCA). Despite finding discordance between mtDNA and individual nuDNA loci, the concatenated matrix recovers a completely resolved and robustly supported phylogeny that is also broadly congruent with BCA trees. This study strongly supports groupings such as Delphininae, Lissodelphininae, Globicephalinae, Sotalia+Delphininae, Steno+Orcaella+Globicephalinae, and Leucopleurus acutus, Lagenorhynchus albirostris, and Orcinus orca as basal delphinid taxa.     

Effects of Asymmetric Nuclear Introgression,Introgressive Mitochondrial Sweep,and Purifying Selection on Phylogenetic Reconstruction and Divergence Estimates in the Pacific Clade of Locustella Warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader''s mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species'' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.     

Geographic mode of speciation in a mountain specialist Avian family endemic to the Palearctic

Mountains host greater avian diversity than lowlands at the same latitude due to their greater diversity of habitats stratified along an elevation gradient. Here we test whether this greater ecological heterogeneity promotes sympatric speciation. We selected accentors (Prunellidae), an avian family associated with mountains of the Palearctic, as a model system. Accentors differ in their habitat/elevation preferences and south‐central Siberia and Himalayan regions each host 6 of the 13 species in the family. We used sequences of the mtDNA ND2 gene and the intron 9 of the Z chromosome specific ACO1 gene to reconstruct a complete species‐level phylogeny of Prunellidae. The tree based on joint analysis of both loci was used to reconstruct the family's biogeographic history and to date the diversification events. We also analyzed the relationship between the node age and sympatry, to determine the geographic mode of speciation in Prunellidae. Our data suggest a Miocene origin of Prunellidae in the Himalayan region. The major division between alpine species (subgenus Laiscopus) and species associated with shrubs (subgenus Prunella) and initial diversification events within the latter happened within the Himalayan region in the Miocene and Pliocene. Accentors colonized other parts of the Palearctic during the Pliocene‐Pleistocene transition. This spread across the Palearctic resulted in rapid diversification of accentors. With only a single exception dating to 0.91 Ma, lineages younger than 1.5 Ma are allopatric. In contrast, sympatry values for older nodes are >0. There was no relationship between node age and range symmetry. Allopatric speciation (not to include peripatric) is the predominant geographic mode of speciation in Prunellidae despite the favorable conditions for ecological diversification in the mountains and range overlaps among species.     

Taxonomic review and distribution of giant water bugs (Hemiptera: Belostomatidae: Lethocerinae) in the Palearctic,Oriental, and Australian regions

The aim of the present study was to analyze the morphology and COI sequences of giant water bugs (Belostomatidae: Lethocerinae) that inhabit the Palearctic, Oriental, and Australian regions (i.e., three Lethocerus and a single Kirkaldyia species). Analysis revealed that both L. patruelis and L. indicus share two light stripes on the pronotum, but L. patruelis possesses narrower stripes and L. insulanus possesses an additional longitudinal narrow dark median band, whereas K. deyrolli was characterized by not having such markings. With regards to genetic analysis, the average intraspecific genetic distances of L. patruelis and L. indicus (from Palearctic and Oriental regions are 0.3 and 0.8%, respectively), but K. deyrolli which was collected from Northeast Asia, showed much higher intraspecific genetic distance (3.7%). The genetic distance between Kirkaldyia and Lethocerus (16.4–16.8%) is similar to that found between the genera of other hemipteran taxa. This study also newly reported the extended easternmost distribution of L. patruelis up to Surat Thani Province in southern Thailand. Along with K. deyrolli, which is considered threatened in South Korea and Japan, Lethocerus species have been facing threats in Southeast Asia, where they are heavily harvested and commonly sold in local markets as food. Therefore, future conservation efforts should be directed toward Lethocerus species in tropical Southeast Asia, especially the rarely found L. patruelis.     

Phylogeny,host association and biogeographical patterns in the diverse millipede‐parasitoid genus Myriophora Brown (Diptera: Phoridae)

Myriophora is the most species‐rich group of parasitoids that attack toxic, chemically defended millipedes in the superorder Juliformia and order Polydesmida—a resource that few insect predators and parasitoids are able to exploit. Worldwide, there are an estimated 200 species of Myriophora, with the majority of the diversity centred in the Neotropical region. The phylogeny of Myriophora is unknown, biogeographical patterns are not documented, and known host associations have not been assessed in a phylogenetic context. We provide the first phylogenetic study of the genus from a data set composed of 52 taxa primarily from the Neotropical region including 10 outgroups, 40 morphological characters, and molecular data from three mitochondrial (16S, COI and ND1) and one nuclear marker (AK). We find that Myriophora dispersed from the New World to the Old World in a single event before subsequently spreading to the Afrotropical region. The ancestral hosts reconstructed for Myriophora are the benzoquinone‐producing Juliformia, and this association has been retained in the Old World clade. In the Neotropical region, Myriophora that are associated with cyanide‐producing polydesmidan millipedes are confined to a single clade that shows remarkably little genetic variation between clearly morphologically diagnosable species.     

Slow Mitochondrial COI Sequence Evolution at the Base of the Metazoan Tree and Its Implications for DNA Barcoding

The evolution rates of mtDNA in early metazoans hold important implications for DNA barcoding. Here, we present a comprehensive analysis of intra- and interspecific COI variabilities in Porifera and Cnidaria (separately as Anthozoa, Hydrozoa, and Scyphozoa) using a data set of 619 sequences from 224 species. We found variation within and between species to be much lower in Porifera and Anthozoa compared to Medusozoa (Hydrozoa and Scyphozoa), which has divergences similar to typical metazoans. Given that recent evidence has shown that fungi also exhibit limited COI divergence, slow-evolving mtDNA is likely to be plesiomorphic for the Metazoa. Higher rates of evolution could have originated independently in Medusozoa and Bilateria or been acquired in the Cnidaria + Bilateria clade and lost in the Anthozoa. Low identification success and substantial overlap between intra- and interspecific COI distances render the Anthozoa unsuitable for DNA barcoding. Caution is also advised for Porifera and Hydrozoa because of relatively low identification success rates as even threshold divergence that maximizes the “barcoding gap” does not improve identification success. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.