Deep sympatric mitochondrial divergence without reproductive isolation in the common redstart Phoenicurus phoenicurus

The published complete mitochondrial genome of Blue-fronted Redstart (Phoenicurus frontalis) is a chimera and includes DNA from Pink-rumped Rosefinch Carpodacus waltoni eos (Aves: Passeriformes)

The complete mitochondrial genome of Blue-fronted Redstart (Phoenicurus frontalis), GenBank accession number MT360379 (NC_053917), was published by Li and colleages in 2020. Here we show that this mitogenome is actually a chimera containing DNA fragments of both P. frontalis (15,518 bp, 92.5%) and Pink-rumped Rosefinch (Carpodacus waltoni eos, 1258 bp, 7.5%). This mitogenome has been re-used in at least three phylogenies. Our study confirms that mitogenomes are best verified with multiple gene trees, and that any anomalies should be investigated by direct comparison of sequences.

Two new freshwater blennies from the Eastern Mediterranean basin (Teleostei: Blenniidae)

Two new species of Salariopsis are described from the Eastern Mediterranean basin. Salariopsis burcuae, new species, from the Bay of Antalya east to the Jordan, is characterised by having a short cirrus, usually not overlapping the 9th circum-orbital sensory pore, and many tiny black dots on the cheek not organised in rows or bands. The new species shows a 4.1% K2P sequence divergence on the cytochrome-c-oxidase subunit 1 (COI) barcoding region from its closest relative, S. fluviatilis. Salariopsis renatorum, new species, from the upper Ceyhan drainage and a coastal stream in Arsuz, is distinguished by having an unbranched supraocular tentacle, black lateral line pores, a short snout, and no black dots on the upper part of the flank and on the cheek. It is also distinguished from its geographically closest congener, S. burcuae, by a molecular distance of 8.8% K2P in its COI barcode region.

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Mitochondrial genomes reveal mid-Pleistocene population divergence, and post-glacial expansion, in Australasian snapper (Chrysophrys auratus)

Glacial cycles play important roles in determining the phylogeographic structure of terrestrial species, however, relatively little is known about their impacts on the distribution of marine biota. This study utilised modern (n = 350) and ancient (n = 26) mitochondrial genomes from Australasian snapper (Chrysophrys auratus) sampled in New Zealand to assess their demographic and phylogeographic history. We also tested for changes in genetic diversity using the up to 750-year-old mitochondrial genomes from pre-European archaeological sites to assess the potential impacts of human exploitation. Nucleotide diversity and haplotype diversity was high (π = 0.005, h = 0.972). There was no significant change in nucleotide diversity over the last 750 years (p = 0.343), with no detectable loss of diversity as a result of indigenous and industrial-scale fishing activity. While there was no evidence for contemporary population structure (AMOVA, p = 0.764), phylogeographic analyses identified two distinct mitochondrial clades that diverged approximately 650,000 years ago during the mid-Pleistocene, suggesting the species experienced barriers to gene flow when sea levels dropped over 120 m during previous glacial maxima. An exponential population increase was also observed around 8000 years ago consistent with a post-glacial expansion, which was likely facilitated by increased ocean temperatures and rising sea levels. This study demonstrates that glacial cycles likely played an important role in the demographic history of C. auratus and adds to our growing understanding of how dynamic climatic changes have influenced the evolution of coastal marine species.

Integrative species delimitation reveals fine-scale allopatric speciation in a good-flying insect: a case study on

Alpha taxonomy is fundamental for many biological fields. Delineation of the species boundary, however, can be challenging in a species complex, where different species share a similar morphology and diagnostic characters may not be available. In this context, integrative approaches that incorporate molecular and morphological data sets, and account for speciation history can be helpful to alpha taxonomy. Different approaches to species delimitation based on different assumptions are complementary and by integrating the results from multiple approaches we can generate a more reliable and objective taxonomic decision. In this study, we applied three molecular approaches to species delimitation and inferred the demographic history based on an isolation with migration model to test a morphologically based taxonomic hypothesis for the Cylindera pseudocylindriformis complex. We discuss the association between genetic divergence and microhabitat specialisation, and further corroborate that C. subtilis sp. nov. is a valid new species by integrating the results from model-based species delimitation and the genealogical divergence index. We argue that genetic endemism can occur at a small geographic scale, even in a winged insect like tiger beetles. Our results also indicated that there may still be undocumented species diversity of Taiwanese Cylindera remaining to be discovered. ZooBank LSID: urn:lsid:zoobank.org:pub:9DEC1432-365C-4872-8D06-73B95F30624F

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Incongruent molecular and morphological variation in the crab spider Synemaglobosum (Araneae, Thomisidae) in Europe

Establishing species boundaries is one of the challenges taxonomists around the world have been tackling for centuries. The relation between intraspecific and interspecific variability is still under discussion and in many taxa it remains understudied. Here the hypothesis of single versus multiple species of the crab spider Synemaglobosum (Fabricius) is tested. The wide distribution range as well as its high morphological variability makes this species an interesting candidate for re-evaluation using an integrative approach. This study combines information from barcoding, phylogenetic reconstruction based on mitochondrial CO1 and ITS2 of more than 60 specimens collected over a wide range of European localities, and morphology. The findings show deep clades with up to 6% mean pairwise distance in the CO1 barcode without any biogeographical pattern. The nuclear ITS2 gene did not support the CO1 clades. Morphological assessment of somatic and genital characters in males and females and a morphometric analysis of the male palp uncovered high intraspecific variation that does not match the CO1 or ITS2 phylogenies or biogeography either. Screening for endosymbiotic Wolbachia bacteria was conducted and only a single infected specimen was found. Several scenarios might explain these inconsistent patterns. While the deep divergences in the barcoding marker might suggest cryptic or ongoing speciation or geographical isolation in the past, the lack of congruent variation in the nuclear ITS2 gene or the studied morphological character systems, especially the male palp, indicates that S.globosum might simply be highly polymorphic both in terms of its mtDNA and morphology. Therefore, more data on ecology and behaviour and full genome sequences are necessary to ultimately resolve this taxonomically intriguing case.

No hybrid snowcocks in the Altai-Hyper-variable markers can be problematic for phylogenetic inference

A recent article in Ecology and Evolution featured the discovery of hybrid snowcocks (Tetraogallus) and speculated on the hybrid origin of an extant species (T. altaicus). Comprehensive re-analyses of original data from the latter paper reliably refute the phylogenetic hypothesis taken as firm evidence of a past hybridization event in these birds. The new re-analyses showed that there is no evidence of hybridization in these snowcocks from the data available so far.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Sympatric occurrence of deeply diverged mitochondrial DNA lineages in Siberian geometrid moths (Lepidoptera: Geometridae): cryptic speciation, mitochondrial introgression, secondary admixture or effect of Wolbachia?

The divergent sympatric mitochondrial lineages within traditionally recognized species present a challenge regularly faced by taxonomists and evolutionary biologists. We encountered this problem when studying the Siberian geometrid moths, Alcis deversata and Thalera chlorosaria. Within each of these species we found two deeply diverged mitochondrial lineages that demonstrated a level of genetic differentiation exceeding the standard interspecific DNA barcode threshold. Using analyses of nuclear genes, morphology, ecological preferences and Wolbachia endosymbionts, we tested five hypotheses that might explain the mitochondrial pattern observed: cryptic speciation, ancestral polymorphism, interspecific mitochondrial introgression, secondary admixture of allopatrically evolved populations and an effect of intracellular Wolbachia endosymbionts. We demonstrate that in A. deversata and Th. chlorosaria the mitochondrial differences are not correlated with differences in nuclear genes, morphology, ecology and Wolbachia infection status, thus not supporting the hypothesis of cryptic species and an effect of Wolbachia. Mitochondrial introgression can lead to a situation in which one species has both its own mitochondrial lineage and the lineage obtained from another species. We found this situation in the species pair Alcis repandata and Alcis extinctaria. We conclude that the mitochondrial heterogeneity in A. deversata and Th. chlorosaria is most likely to be attributable to the secondary admixture of allopatrically evolved populations.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Molecular Species Delimitation of Larks (Aves: Alaudidae), and Integrative Taxonomy of the Genus Calandrella, with the Description of a Range-Restricted African Relic Taxon

Larks constitute an avian family of exceptional cryptic diversity and striking examples of convergent evolution. Therefore, traditional morphology-based taxonomy has recurrently failed to reflect evolutionary relationships. While taxonomy ideally should integrate morphology, vocalizations, behaviour, ecology, and genetics, this can be challenging for groups that span several continents including areas that are difficult to access. Here, we combine morphometrics and mitochondrial DNA to evaluate the taxonomy of Calandrella larks, with particular focus on the African C. cinerea and the Asian C. acutirostris complexes. We describe a new range-restricted West African taxon, Calandrella cinerea rufipecta ssp. nov. (type locality: Jos, Plateau State, Nigeria), with an isolated relic population 3000 km from its closest relative in the Rift Valley. We performed molecular species delimitation, employing coalescence-based multi-rate Poisson Tree Processes (mPTP) on cytochrome b sequences across 52 currently recognized lark species, including multiple taxa currently treated as subspecies. Three species-level splits were inferred within the genus Calandrella and another 13 across other genera, primarily among fragmented sub-Saharan taxa and taxa distributed from Northwest Africa to Arabia or East Africa. Previously unknown divergences date back as far as to the Miocene, indicating the presence of currently unrecognized species. However, we stress that taxonomic decisions should not be based on single datasets, such as mitochondrial DNA, although analyses of mitochondrial DNA can be a good indicator of taxa in need of further integrative taxonomic assessment.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

Hybridization between two recently diverged Neotropical passerines: The Pearly-bellied Seedeater Sporophila pileata, and the Copper Seedeater S. bouvreuil (Aves, Passeriformes, Thraupidae)

The small Neotropical finches called capuchinos are outstanding because they have experienced one of the most recent and explosive avian radiations ever documented for birds. Despite very low morphological and niche divergence among species, many of them are reproductively isolated when in sympatry due to strong sexual selection in plumage traits. However, a specific pair of mostly parapatric species, the Pearly-bellied, Sporophila pileata, and the Copper Seedeaters, S. bouvreuil, has confounded taxonomists because individuals with intermediate color patterns can be found. By analyzing diagnostic COI mtDNA sequences and adult male plumage we provide evidence for hybridization. Paternity tests using microsatellites also indicated that representatives with intermediate plumage pattern can be fertile. Our findings are consistent with the classification of S. bouvreuil and S. pileata as distinct taxa, but we demonstrate that the sexual selection mechanisms involved in the isolation of other reproductively sympatric capuchinos are not applicable to this pair of species, likely because of reduced barriers to mate recognition.

Phylogeography of Lionfishes (Pterois) Indicate Taxonomic Over Splitting and Hybrid Origin of the Invasive Pterois volitans

The lionfish is an iconic marine fish, and recently renowned for a disastrous introduction into the West Atlantic. Genetic surveys of the putative invaders (Pterois volitans and Pterois miles) in their natural Indo-Pacific range can illuminate both topics. Previous research indicated that P. volitans and P. miles are sister species that hybridize in the invasive range, but hybridization in the native range is unknown. Here, we apply mtDNA COI and 2 nuclear introns (S7 RP1 and Gpd2) from 229 lionfish including the 2 invaders and 2 closely-related taxa (44 P. miles, 91 P. volitans, 31 Pterois lunulata, and 63 Pterois russelii) from 10 locations in their native ranges. Genetic data are supplemented with key morphological characters: dorsal, anal, and pectoral fin ray counts. We observed 2 lineages (d = 4.07%, 0.89%, and 2.75% at COI, S7 RP1, and Gpd2, respectively) among the 4 putative species: an Indian Ocean lineage represented by P. miles, and a Pacific Ocean lineage represented by P. lunulata and P. russelii. All specimens of the invasive P. volitans appear to be hybrids between the Indian Ocean P. miles and a Pacific lineage encompassing P. lunulata/russelii, a conclusion supported by both genetics and morphology. The divergences between Indian and Pacific forms are within the range of species-level partitions in fishes, and we recommend retention of the names P. miles and P. russelii for Indian and Pacific forms. The hybrid origin of the Atlantic invasion invokes the possibility of heterosis as a contributing factor to invasion success.

'Ghost introgression' as a cause of deep mitochondrial divergence in a bird species complex

In the absence of nuclear-genomic differentiation between two populations, deep mitochondrial divergence (DMD) is a form of mito-nuclear discordance. Such instances of DMD are rare and might variably be explained by unusual cases of female-linked selection, by male-biased dispersal, by 'speciation reversal' or by mitochondrial capture through genetic introgression. Here we analyze DMD in an Asian Phylloscopus leaf warbler (Aves: Phylloscopidae) complex. Bioacoustic, morphological and genomic data demonstrate close similarity between the taxa affinis and occisinensis, even though DMD previously led to their classification as two distinct species. Using population genomic and comparative genomic methods on 45 whole genomes, including historical reconstructions of effective population size, genomic peaks of differentiation and genomic linkage, we infer that the form affinis is likely the product of a westward expansion in which it replaced a now-extinct congener that was the donor of its mtDNA and small portions of its nuclear genome. This study provides strong evidence of 'ghost introgression' as the cause of DMD, and we suggest that 'ghost introgression' may be a widely overlooked phenomenon in nature.

Divergent mitochondrial lineages arose within a large, panmictic population of the Savannah sparrow (Passerculus sandwichensis)

Unusual patterns of mtDNA diversity can reveal interesting aspects of a species' biology. However, making such inferences requires discerning among the many alternative scenarios that could underlie any given mtDNA pattern. Next-generation sequencing methods provide large, multilocus data sets with increased power to resolve unusual mtDNA patterns. A mtDNA-based phylogeography of the Savannah sparrow (Passerculus sandwichensis) previously identified two sympatric, but divergent (~2%) clades within the nominate subspecies group and a third clade that consisted of birds sampled from northwest Mexico. We revisited the phylogeography of this species using a population genomic data set to resolve the processes leading to the evolution of sympatric and divergent mtDNA lineages. We identified two genetic clusters in the genomic data set corresponding to (a) the nominate subspecies group and (b) northwestern Mexico birds. Following divergence, the nominate clade maintained a large, stable population, indicating that divergent mitochondrial lineages arose within a panmictic population. Simulations based on parameter estimates from this model further confirmed that this demographic history could produce observed levels of mtDNA diversity. Patterns of divergent, sympatric mtDNA lineages are frequently interpreted as admixture of historically isolated lineages. Our analyses reject this interpretation for Savannah sparrows and underscore the need for genomic data sets to resolve the evolutionary mechanisms behind anomalous, locus-specific patterns.

Extreme Mitogenomic Variation in Natural Populations of Chaetognaths

The extent of within-species genetic variation across the diversity of animal life is an underexplored problem in ecology and evolution. Although neutral genetic variation should scale positively with population size, mitochondrial diversity levels are believed to show little variation across animal species. Here, we report an unprecedented case of extreme mitochondrial diversity within natural populations of two morphospecies of chaetognaths (arrow worms). We determine that this diversity is composed of deep sympatric mitochondrial lineages, which are in some cases as divergent as human and platypus. Additionally, based on 54 complete mitogenomes, we observed mitochondrial gene order differences between several of these lineages. We examined nuclear divergence patterns (18S, 28S, and an intron) to determine the possible origin of these lineages, but did not find congruent patterns between mitochondrial and nuclear markers. We also show that extreme mitochondrial divergence in chaetognaths is not driven by positive selection. Hence, we propose that the extreme levels of mitochondrial variation could be the result of either a complex scenario of reproductive isolation, or a combination of large population size and accelerated mitochondrial mutation rate. These findings emphasize the importance of characterizing genome-wide levels of nuclear variation in these species and promote chaetognaths as a remarkable model to study mitochondrial evolution.

The absence of reproductive isolation between non-sister and deeply diverged mitochondrial lineages of the black-throated tit (Aegithalos concinnus) revealed by a multilocus genetic analysis in a contact zone

BACKGROUND

A deep divergence of mitochondrial DNA is common in species delimitated by morphological traits. Several hypotheses can explain such variations, such as cryptic species, introgression, allopatric divergence and ancestral lineage. The black-throated tit harbors several deeply divergent mitochondrial lineages. Two lineages with 5% divergence, but having a high level of gene flow, have been detected in its subspecies, A. C. concinnus and A. c. talifuenses. In this study, we conducted a genetic analysis at a contact zone of these two lineages to identify whether these lineages either reflect a high level of intraspecific variation in mitochondrial loci or represent incipient speciation. Mitochondrial ND2 and 11 microsatellite loci were used to conduct phylogenetic and population structure analyses.

RESULTS

ND2 haplotypes actually diverged into two groups within subspecies A. c. talifuenses; however, they formed a non-sister relationship when including all available GenBank ND2 sequences. Analyses of microsatellite data indicated no existing population structure and showed a pattern of isolation by distance. Individuals sampled at the contact zone were almost identified as F2 hybrids.

CONCLUSIONS

Isolation for 2.4 Ma, as suggested by a previous study, appeared to be insufficient to develop robust reproductive barriers. Reproductive barriers were weak, or even absent between the divergent lineages, highlighting that incipient speciation was unlikely to be the case. Considering the results from previous studies, the divergent lineages may be better explained by secondary contact after allopatric isolation because of Pleistocene climate changes, but other hypotheses cannot be definitively ruled out because of the lack of representative samples from the other distribution region and its relatives. Considering that divergence in morphology was hardly observed and the particular split in genetics, the two subspecies might be better merged. The current findings also highlight the likely contribution of male-biased dispersal to male-biased gene flow among mitochondrial lineages; more efforts are needed to illustrate the evolutionary history of the black-throated tit.

Integration of Genetic and Phenotypic Data in 48 Lineages of Philippine Birds Shows Heterogeneous Divergence Processes and Numerous Cryptic Species

The Philippine Islands are one of the most biologically diverse archipelagoes in the world. Current taxonomy, however, may underestimate levels of avian diversity and endemism in these islands. Although species limits can be difficult to determine among allopatric populations, quantitative methods for comparing phenotypic and genotypic data can provide useful metrics of divergence among populations and identify those that merit consideration for elevation to full species status. Using a conceptual approach that integrates genetic and phenotypic data, we compared populations among 48 species, estimating genetic divergence (p-distance) using the mtDNA marker ND2 and comparing plumage and morphometrics of museum study skins. Using conservative speciation thresholds, pairwise comparisons of genetic and phenotypic divergence suggested possible species-level divergences in more than half of the species studied (25 out of 48). In speciation process space, divergence routes were heterogeneous among taxa. Nearly all populations that surpassed high genotypic divergence thresholds were Passeriformes, and non-Passeriformes populations surpassed high phenotypic divergence thresholds more commonly than expected by chance. Overall, there was an apparent logarithmic increase in phenotypic divergence with respect to genetic divergence, suggesting the possibility that divergence among these lineages may initially be driven by divergent selection in this allopatric system. Also, genetic endemism was high among sampled islands. Higher taxonomy affected divergence in genotype and phenotype. Although broader lineage, genetic, phenotypic, and numeric sampling is needed to further explore heterogeneity among divergence processes and to accurately assess species-level diversity in these taxa, our results support the need for substantial taxonomic revisions among Philippine birds. The conservation implications are profound.

Species-Level Para- and Polyphyly in DNA Barcode Gene Trees: Strong Operational Bias in European Lepidoptera

The proliferation of DNA data is revolutionizing all fields of systematic research. DNA barcode sequences, now available for millions of specimens and several hundred thousand species, are increasingly used in algorithmic species delimitations. This is complicated by occasional incongruences between species and gene genealogies, as indicated by situations where conspecific individuals do not form a monophyletic cluster in a gene tree. In two previous reviews, non-monophyly has been reported as being common in mitochondrial DNA gene trees. We developed a novel web service “Monophylizer” to detect non-monophyly in phylogenetic trees and used it to ascertain the incidence of species non-monophyly in COI (a.k.a. cox1) barcode sequence data from 4977 species and 41,583 specimens of European Lepidoptera, the largest data set of DNA barcodes analyzed from this regard. Particular attention was paid to accurate species identification to ensure data integrity. We investigated the effects of tree-building method, sampling effort, and other methodological issues, all of which can influence estimates of non-monophyly. We found a 12% incidence of non-monophyly, a value significantly lower than that observed in previous studies. Neighbor joining (NJ) and maximum likelihood (ML) methods yielded almost equal numbers of non-monophyletic species, but 24.1% of these cases of non-monophyly were only found by one of these methods. Non-monophyletic species tend to show either low genetic distances to their nearest neighbors or exceptionally high levels of intraspecific variability. Cases of polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, as the detected cases reflected misidentifications or methodological errors. Taking into consideration variation in sampling effort, we estimate that the true incidence of non-monophyly is ∼23%, but with operational factors still being included. Within the operational factors, we separately assessed the frequency of taxonomic limitations (presence of overlooked cryptic and oversplit species) and identification uncertainties. We observed that operational factors are potentially present in more than half (58.6%) of the detected cases of non-monophyly. Furthermore, we observed that in about 20% of non-monophyletic species and entangled species, the lineages involved are either allopatric or parapatric—conditions where species delimitation is inherently subjective and particularly dependent on the species concept that has been adopted. These observations suggest that species-level non-monophyly in COI gene trees is less common than previously supposed, with many cases reflecting misidentifications, the subjectivity of species delimitation or other operational factors.

Barcoding gap, but no support for cryptic speciation in the earthworm Aporrectodea longa (Clitellata: Lumbricidae)

DNA-barcoding, using the mitochondrial marker COI, has been found successful for the identification of specimens in many animal groups, but may not be suited for species discovery and delimitation if used alone. In this study, we investigate whether two observed COI haplogroups in the earthworm Aporrectodea longa correspond to two cryptic species or if the variation is intraspecific. This is done by complementing COI with two nuclear markers, ITS2 and Histone 3. The variation is studied using distance methods, parsimony networks and Bayesian coalescent trees, and the statistical distinctness of the groups is tested on gene trees using the genealogical sorting index, Rosenberg's P_AB and Rodrigo et al.'s P_(RD). We also applied multilocus species delimitation based on the multispecies coalescence model. The two haplogroups were found in COI, and all tests except P_(RD) found them to be significantly distinct. However, in ITS2, the same groups were not recovered in any analyses or tests. H3 was invariable in A. longa, and was, therefore, included only in the multilocus analysis, which preferred a model treating A. longa as one species over a model splitting it into two. We also compared two measurements of size, body length, and no. of segments between the groups. No difference in body length was found, and although a significant difference in no. of segments was noted the haplogroup with the lower mean showed both the highest and the lowest value. When combined, these results led us to the conclusion that there is no support for the separation of A. longa into two cryptic species. This study again highlights the importance of complementing mitochondrial barcodes with more data when establishing species boundaries.

Deeply divergent sympatric mitochondrial lineages of the earthworm Lumbricus rubellus are not reproductively isolated

Background

The accurate delimitation of species is essential to numerous areas of biological research. An unbiased assessment of the diversity, including the cryptic diversity, is of particular importance for the below ground fauna, a major component of global biodiversity. On the British Isles, the epigeic earthworm Lumbricus rubellus, which is a sentinel species in soil ecotoxicology, consists of two cryptic taxa that are differentiated in both the nuclear and the mitochondrial (mtDNA) genomes. Recently, several deeply divergent mtDNA lineages were detected in mainland Europe, but whether these earthworms also constitute cryptic species remains unclear. This information is important from an evolutionary perspective, but it is also essential for the interpretation and the design of ecotoxicological projects. In this study, we used genome-wide RADseq data to assess the reproductive isolation of the divergent mitochondrial lineages of L. rubellus that occur in sympatry in multiple localities in Central Europe.

Results

We identified five divergent (up to 16 % net p-distance) mitochondrial lineages of L. rubellus in sympatry. Because the clustering of the RADseq data was according to the population of origin and not the mtDNA lineage, reproductive isolation among the mtDNA lineages was not likely. Although each population contained multiple mtDNA lineages, subdivisions within the populations were not observed for the nuclear genome. The lack of fixed differences and sharing of the overwhelming majority of nuclear polymorphisms between localities, indicated that the populations did not constitute allopatric species. The nucleotide diversity within the populations was high, 0.7–0.8 %.

Conclusions

The deeply divergent mtDNA sympatric lineages of L. rubellus in Central Europe were not reproductively isolated groups. The earthworm L. rubellus, which is represented by several mtDNA lineages in continental Europe, apparently is a single highly polymorphic species rather than a complex of several cryptic species. This study demonstrated the critical importance of the use of multilocus nuclear data for the unbiased assessment of cryptic diversity and for the delimitation of species in soil invertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0488-9) contains supplementary material, which is available to authorized users.

Potential merger of ancient lineages in a passerine bird discovered based on evidence from host-specific ectoparasites

The merger of formerly isolated lineages is hypothesized to occur in vertebrates under certain conditions. However, despite many demonstrated instances of introgression between taxa in secondary contact, examples of lineage mergers are rare. Preliminary mtDNA sequencing of a Malagasy passerine, Xanthomixis zosterops (Passeriformes: Bernieridae), indicated a possible instance of merging lineages. We tested the hypothesis that X. zosterops lineages are merging by comparing mtDNA sequence and microsatellite data, as well as mtDNA sequence data from host-specific feather lice in the genus Myrsidea (Phthiraptera: Menoponidae). Xanthomixis zosterops comprises four deeply divergent, broadly sympatric, cryptic mtDNA clades that likely began diverging approximately 3.6 million years ago. Despite this level of divergence, the microsatellite data indicate that the X. zosterops mtDNA clades are virtually panmictic. Three major phylogroups of Myrsidea were found, supporting previous allopatry of the X. zosterops clades. In combination, the datasets from X. zosterops and its Myrsidea document a potential merger of previously allopatric lineages that likely date to the Pliocene. This represents the first report of sympatric apparent hybridization among more than two terrestrial vertebrate lineages. Further, the mtDNA phylogeographic pattern of X. zosterops, namely the syntopy of more than two deeply divergent cryptic clades, appears to be a novel scenario among vertebrates. We highlight the value of gathering multiple types of data in phylogeographic studies to contribute to the study of vertebrate speciation.

Whole plastome sequencing reveals deep plastid divergence and cytonuclear discordance between closely related balsam poplars, Populus balsamifera and P. trichocarpa (Salicaceae)

As molecular phylogenetic analyses incorporate ever-greater numbers of loci, cases of cytonuclear discordance - the phenomenon in which nuclear gene trees deviate significantly from organellar gene trees - are being reported more frequently. Plant examples of topological discordance, caused by recent hybridization between extant species, are well known. However, examples of branch-length discordance are less reported in plants relative to animals. We use a combination of de novo assembly and reference-based mapping using short-read shotgun sequences to construct a robust phylogeny of the plastome for multiple individuals of all the common Populus species in North America. We demonstrate a case of strikingly high plastome divergence, in contrast to little nuclear genome divergence, in two closely related balsam poplars, Populus balsamifera and Populus trichocarpa (Populus balsamifera ssp. trichocarpa). Previous studies with nuclear loci indicate that the two species (or subspecies) diverged since the late Pleistocene, whereas their plastomes indicate deep divergence, dating to at least the Pliocene (6-7 Myr ago). Our finding is in marked contrast to the estimated Pleistocene divergence of the nuclear genomes, previously calculated at 75 000 yr ago, suggesting plastid capture from a 'ghost lineage' of a now-extinct North American poplar.

DNA barcoding reveals 24 distinct lineages as cryptic bird species candidates in and around the Japanese Archipelago

DNA barcoding using a partial region (648 bp) of the cytochrome c oxidase I (COI) gene is a powerful tool for species identification and has revealed many cryptic species in various animal taxa. In birds, cryptic species are likely to occur in insular regions like the Japanese Archipelago due to the prevention of gene flow by sea barriers. Using COI sequences of 234 of the 251 Japanese-breeding bird species, we established a DNA barcoding library for species identification and estimated the number of cryptic species candidates. A total of 226 species (96.6%) had unique COI sequences with large genetic divergence among the closest species based on neighbour-joining clusters, genetic distance criterion and diagnostic substitutions. Eleven cryptic species candidates were detected, with distinct intraspecific deep genetic divergences, nine lineages of which were geographically separated by islands and straits within the Japanese Archipelago. To identify Japan-specific cryptic species from trans-Paleartic birds, we investigated the genetic structure of 142 shared species over an extended region covering Japan and Eurasia; 19 of these species formed two or more clades with high bootstrap values. Excluding six duplicated species from the total of 11 species within the Japanese Archipelago and 19 trans-Paleartic species, we identified 24 species that were cryptic species candidates within and surrounding the Japanese Archipelago. Repeated sea level changes during the glacial and interglacial periods may be responsible for the deep genetic divergences of Japanese birds in this insular region, which has led to inconsistencies in traditional taxonomies based on morphology.

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.