Divergence in calls but not songs in the orchard oriole complex: Icterus spurius and I. fuertesi

Birdsong has important functions in attracting and competing for mates, and song characteristics are thought to diverge rapidly during the process of speciation. In contrast, other avian vocalizations that may have non‐reproductive functions, such as calls, are thought to be more evolutionarily conserved and may diverge more slowly among taxa. This study examines differences in both male song and an acoustically simpler vocalization, the ‘jeet’ call, between two closely related taxa, Icterus spurius and I. fuertesi. A previous study comparing song syllable type sharing within and between I. spurius and I. fuertesi indicated that their songs do not differ discernibly. Here we measured 18 acoustic characteristics of their songs and found strong evidence supporting this prior finding. In contrast, we measured 17 acoustic characteristics of jeet calls and found evidence of significant divergence between the two taxa in many of these characteristics. Calls in I. fuertesi have a longer duration, a larger frequency bandwidth, a lower minimum frequency, a lower beginning frequency, and greater levels of both frequency and amplitude modulation in comparison to the calls of I. spurius. In addition, I. fuertesi calls contain two distinct parts, while the calls of I. spurius have only one part. Thus, we find evidence of divergence in the calls of the two taxa but not their songs challenging the widespread assumption that complex bird song evolves more rapidly than other types of vocalizations. Understanding divergence in multiple vocalization types as well as other behavioral, morphological, and molecular traits is important to understanding the earliest stages of speciation.     

Incomplete song divergence between recently diverged taxa: syllable sharing by Orchard and Fuertes’ orioles

ABSTRACT Orchard Orioles (Icterus spurius) and Fuertes’ Orioles (I. fuertesi) recently diverged from each other, making them an ideal system for investigating trait evolution and mechanisms of reproductive isolation during the early stages of speciation. These taxa differ in adult male plumage coloration and in their breeding and wintering ranges, but quantitative comparisons of their song characteristics have revealed no discernible differences. We assessed evolutionary song divergence in this group by investigating patterns of syllable‐type sharing within and between populations. Of 529 distinct syllable types, 142 (26.8%) were shared among individuals, and sharing appeared to decrease with geographic distance. The total number of syllables shared between Orchard and Fuertes’ orioles (26; 4.9% of the total) was similar to levels of sharing between populations of Orchard Orioles. Furthermore, hierarchical cluster analyses showed individuals of the two taxa intermixed. Syllables also used as calls were shared more frequently within and between taxa, suggesting that they have evolved more slowly than those used exclusively in songs. Our results show that at least some aspects of song have not yet diverged between these incipient species, either due to cultural exchange or because songs have evolved relatively slowly compared to plumage colors.     

Coalescent analyses show isolation without migration in two closely related tropical orioles: the case of Icterus graduacauda and Icterus chrysater

The Isthmus of Tehuantepec has played an important role in shaping the avian diversity of Mexico, as well as the rest of the Western Hemisphere. It has been both a barrier and a land connector between North and South America for many groups of birds. Furthermore, climatic change over the Pleistocene has resulted in ecological fluctuations that led to periods of connection and isolation of the highlands in this area. Here we studied the divergence of two species of orioles whose distribution in the highlands is separated by the lowlands of the Isthmus of Tehuantepec: Icterus graduacauda (west of the Isthmus) and Icterus chrysater (east of the Isthmus). We sequenced multiple loci (one mitochondrial gene and six nuclear introns) and performed coalescent analyses (Isolation with Migration) to test whether their divergence resulted from prior occupancy of the ancestral area followed by a vicariant event or recent dispersal from one side or the other of this Isthmus. Results strongly indicate a vicariant event roughly 300,000 years ago in the Pleistocene followed by little or no gene flow. Both mitochondrial and nuclear genes show that the Isthmus of Tehuantepec is a strong barrier to gene flow. Thus, these two species appear to not exchange genes despite their recent divergence and the close geographic proximity of their ranges.     

Assessment of species limits in African ‘brown buntings’ (Emberiza,Passeriformes) based on mitochondrial and nuclear sequence data

We estimated a phylogeny for 10 taxa currently placed in four polytypic species that collectively encompass the African ‘brown buntings’: Cape Bunting Emberiza capensis, Cinnamon‐breasted Bunting Emberiza tahapisi, Lark‐like Bunting Emberiza impetuani and House Bunting Emberiza striolata. We made use of the mitochondrial cytochrome b gene and the nuclear introns 6–7 of ornithine decarboxylase (ODC), and intron 2 of myoglobin. There was substantial cytochrome b sequence divergence between taxa currently treated as conspecific: sahari vs. striolata (2.6–3.1% (uncorrected‐p); 3.0–3.6% (HKY + I)), and goslingi vs. tahapisi (4.4–4.7% (uncorrected‐p); 5.4–5.9% (HKY + I)). The degree of divergence of the nuclear loci among taxa was limited, and these loci lacked reciprocal monophyly, most likely as a consequence of incomplete lineage sorting. A single representative of the taxon septemstriata, generally treated as a member of the dark‐throated tahapisi group, here appears to be genetically consistent with the grey‐throated goslingi, and may be of hybrid origin. All other taxa allocated to E. striolata and E. tahapisi make up four reciprocally monophyletic groups consistent with sahari, striolata, tahapisi and goslingi, respectively. The extent of genetic evidence suggests that these taxa have been evolving as separate evolutionary lineages for a long time. This is further manifested in several morphological and vocal characteristics described previously, and we propose that these divergent taxa be treated as separate species: Cinnamon‐breasted Bunting Emberiza tahapisi, Gosling's Bunting Emberiza goslingi, Striolated Bunting Emberiza striolata and House Bunting Emberiza sahari. We do not propose any taxonomic changes regarding Emberiza impetuani or Emberiza capensis.     

Resolving species boundaries in the critically imperiled freshwater mussel species,Fusconaia mitchelli (Bivalvia: Unionidae)

Species are a fundamental unit of biology, and defining accurate species boundaries is integral to effective conservation and management of imperiled taxa. Freshwater mussels (Bivalvia: Unionidae) are among the most imperiled groups of organisms in North America, yet species boundaries remain uncertain for many taxa. The False Spike, Fusconaia mitchelli (Simpson in Dall, 1895), is a freshwater mussel considered to be endemic to central Texas (Brazos, Colorado, and Guadalupe drainages). Recent research revealed significant intraspecific genetic variation between geographically separated populations of F. mitchelli, which could be indicative of speciation; however, small sample sizes for several of the populations precluded formal taxonomic revision. Here, we increase taxon sampling and use multilocus DNA sequence data and traditional morphometrics to re‐evaluate species boundaries in F. mitchelli. We sequenced three loci: the protein‐coding mitochondrial DNA genes cytochrome c oxidase subunit 1 and NADH dehydrogenase 1, and the nuclear internal transcribed spacer 1. Phylogenetic analyses depicted deep genetic divergence between F. mitchelli in the Guadalupe and those in the Brazos and Colorado drainages, which was further supported by available biogeographic information. Morphometric analyses and coalescent‐based species delimitation models integrating both DNA sequence and morphological data provided strong support for the divergence observed between the two geographically isolated clades of F. mitchelli. Based on these results, we revise taxonomy accordingly by elevating the junior synonym Fusconaia iheringi (Wright, 1898) to represent the Brazos and Colorado populations and restrict the distribution of F. mitchelli to the Guadalupe River drainage. Our findings may impact pending management decisions to protect F. mitchelli under the U.S. Endangered Species Act.     

Extensive introgressive hybridization within the northern oriole group (Genus Icterus) revealed by three‐species isolation with migration analysis

Until recently, studies of divergence and gene flow among closely‐related taxa were generally limited to pairs of sister taxa. However, organisms frequently exchange genes with other non‐sister taxa. The “northern oriole” group within genus Icterus exemplifies this problem. This group involves the extensively studied hybrid zone between Baltimore oriole (Icterus galbula) and Bullock's oriole (I. bullockii), an alleged hybrid zone between I. bullockii and black‐backed oriole (I. abeillei), and likely mtDNA introgression between I. galbula and I. abeillei. Here, we examine the divergence population genetics of the entire northern oriole group using a multipopulation Isolation‐with‐Migration (IM) model. In accordance with Haldane's rule, nuclear loci introgress extensively beyond the I. galbula–I. bullockii hybrid zone, while mtDNA does not. We found no evidence of introgression between I. bullockii and I. abeillei or between I. galbula and I. abeillei when all three species were analyzed together in a three‐population model. However, traditional pairwise analysis suggested some nuclear introgression from I. abeillei into I. galbula, probably reflecting genetic contributions from I. bullockii unaccounted for in a two‐population model. Thus, only by including all members of this group in the analysis was it possible to rigorously estimate the level of gene flow among these three closely related species.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Understanding the mechanisms of antitropical divergence in the seabird White‐faced Storm‐petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White‐faced Storm‐petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial Φ_ST = 0.89 (P < 0.01) and global nuclear Φ_ST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long‐term antitropical divergence in isolation. Coalescent‐based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000–1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000–600 000), both probably associated with the profound palaeo‐oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long‐distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.     

Hidden diversity and cryptic speciation refute cosmopolitan distribution in Caprella penantis (Crustacea: Amphipoda: Caprellidae)

Caprella penantis is considered a cosmopolitan species and one of the most challenging caprellids in taxonomic terms because of its remarkable intraspecific morphological variation. This study examined DNA sequences from mitochondrial (COI) and nuclear (18S) markers together with morphological data from 25 localities of C. penantis, and closely related species Caprella dilatata and Caprella andreae, all traditionally considered part of the old ‘acutifrons’ complex. The large genetic divergence and reciprocally allopatric distributions point to the existence of a species complex of at least four species, of which one is reported as a cryptic species. This study provides the first evidence of cryptic speciation in the family Caprellidae, and questions the validity of some traditional morphological characters used to delimit species in the genus Caprella. Our results are consistent with the idea that main factors were probably isolation by distance and ecological traits, promoting diversification in C. penantis. The strong genetic structure reported for this species in the Iberian Peninsula and Moroccan coasts also suggests restriction to dispersal as well as the presence of refugial areas. These results highlight the utility of the COI and 18S genes in combination with morphological characters for shedding light on systematic questions in caprellids, and patterns of genetic connectivity.     

Evidence for ecological speciation via a host shift in the holly leaf miner,Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host‐associated radiations. We investigated whether host‐associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host‐associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF‐1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.     

Mitochondrial introgressive hybridization following a demographic expansion in the tomato frogs of Madagascar,genus Dyscophus

Madagascar is a biodiversity hotspot with a unique fauna and flora largely endemic at the species level and highly threatened by habitat destruction. The processes underlying population‐level differentiation in Madagascar's biota are poorly understood and have been proposed to be related to Pleistocene climatic cycles, yet the levels of genetic divergence observed are often suggestive of ancient events. We combined molecular markers of different variability to assess the phylogeography of Madagascar's emblematic tomato frogs (Dyscophus guineti and D. antongilii) and interpret the observed pattern as resulting from ancient and recent processes. Our results suggest that the initial divergence between these taxa is probably old as reflected by protein‐coding nuclear genes and by a strong mitochondrial differentiation of the southernmost population. Dramatic changes in their demography appear to have been triggered by the end of the last glacial period and possibly by the short return of glacial conditions known as the 8K event. This dramatic change resulted in an approximately 50‐fold reduction of the effective population size in various populations of both species. We hypothesize these species' current mitochondrial DNA diversity distribution reflects a swamping of the mitochondrial genetic diversity of D. guineti by that of D. antongilii previous to the populations' bottlenecks during the Holocene, and probably as a consequence of D. antongilii demographic expansion approximately 1 million years ago. Our data support the continued recognition of D. antongilii and D. guineti as separate species and flag D. guineti as the more vulnerable species to past and probably also future environmental changes.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates

Ecological and evolutionary studies largely assume that island populations display low levels of neutral genetic variation. However, this notion has only been formally tested in a few cases involving plant taxa, and the confounding effect of selection on genetic diversity (GD) estimates based on putatively neutral markers has typically been overlooked. Here, we generated nuclear microsatellite and plastid DNA sequence data in Periploca laevigata, a plant taxon with an island–mainland distribution area, to (i) investigate whether selection affects GD estimates of populations across contrasting habitats; and (ii) test the long‐standing idea that island populations have lower GD than their mainland counterparts. Plastid data showed that colonization of the Canary Islands promoted strong lineage divergence within P. laevigata, which was accompanied by selective sweeps at several nuclear microsatellite loci. Inclusion of loci affected by strong divergent selection produced a significant downward bias in the GD estimates of the mainland lineage, but such underestimates were substantial (>14%) only when more than one loci under selection were included in the computations. When loci affected by selection were removed, we did not find evidence that insular Periploca populations have less GD than their mainland counterparts. The analysis of data obtained from a comprehensive literature survey reinforced this result, as overall comparisons of GD estimates between island and mainland populations were not significant across plant taxa (N = 66), with the only exception of island endemics with narrow distributions. This study suggests that identification and removal of markers potentially affected by selection should be routinely implemented in estimates of GD, particularly if different lineages are compared. Furthermore, it provides compelling evidence that the expectation of low GD cannot be generalized to island plant populations.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

DNA barcoding reveals taxonomic conflicts in the Herichthys bartoni species group (Pisces: Cichlidae)

Fishes of the genus Herichthys are the only representatives of the family Cichlidae to have colonized the Neartic region. In this study, we used DNA barcode of 64 individuals of the Herichthys bartoni species group to test the monophyly of the species and the efficiency of this tool to discriminate among species. The Bayesian phylogenetic tree and the neighbour joining (NJ) tree obtained from the Kimura two‐parameter model (K2P) give similar topologies. DNA barcoding resolution was very poor (25%). Additionally, the low levels of genetic divergence among taxa preclude the use of threshold values as has been suggested in earlier studies.     

Comparative phylogeographies of six species of hinged terrapins (Pelusios spp.) reveal discordant patterns and unexpected differentiation in the P. castaneus/P. chapini complex and P. rhodesianus

Using up to 2117 bp of mitochondrial DNA and up to 2012 bp of nuclear DNA, we analysed phylogeographic differentiation of six widely distributed species of African hinged terrapins (Pelusios spp.) representing different habitat types. Two taxa each live in savannahs or in forests and mesic savannahs, respectively, and the remaining two species occur in intermediate habitats. The species living in forests and mesic savannahs do not enter dry savannahs, whereas the savannah species may occur in dry and wet savannahs and even in semi‐arid steppe regions. We found no obvious correlation between habitat type and phylogeographic pattern: one savannah species (P. rhodesianus) shows phylogeographic structure, i.e. pronounced genetic differences among geographically distinct populations, and the other (P. nanus) not. One species inhabiting forests and mesic savannahs (P. carinatus) has phylogeographic structure, the other (P. gabonensis) not. The same pattern is true for the two ecologically intermediate species, with phylogeographic structure present in P. castaneus and absent in P. chapini. Nuclear evidence suggests that the latter two taxa with abutting and partially overlapping ranges are distinct, while mtDNA is only weakly differentiated. Pelusios castaneus shows pronounced phylogeographic structure, which could reflect Pleistocene range interruptions correlated with the fluctuating forest cover in West and Central Africa. Our results do not support the recognition of an extinct subspecies of P. castaneus for the Seychelles. Pelusios carinatus contains two well supported clades, which are separated by the Congo River. This species is closely related to P. rhodesianus, a taxon consisting of two deeply divergent mitochondrial clades. One of these clades is paraphyletic with respect to P. carinatus, but the two clades of P. rhodesianus are not differentiated in the studied nuclear markers and, again, paraphyletic with respect to P. carinatus. Using mtDNA sequences from the type material of P. rhodesianus, we were able to allocate this name to one of the two clades. However, owing to the confusing relationships of P. rhodesianus and P. carinatus, we refrain from taxonomic decisions.     

The curious case of Hermodice carunculata (Annelida: Amphinomidae): evidence for genetic homogeneity throughout the Atlantic Ocean and adjacent basins

Over the last few decades, advances in molecular techniques have led to the detection of strong geographic population structure and cryptic speciation in many benthic marine taxa, even those with long‐lived pelagic larval stages. Polychaete annelids, in particular, generally show a high degree of population divergence, especially in mitochondrial genes. Rarely have molecular studies confirmed the presence of ‘cosmopolitan’ species. The amphinomid polychaete Hermodice carunculata was long considered the sole species within its genus, with a reported distribution throughout the Atlantic and adjacent basins. However, recent studies have indicated morphological differences, primarily in the number of branchial filaments, between the East and West Atlantic populations; these differences were invoked to re‐instate Hermodice nigrolineata, formerly considered a junior synonym of H. carunculata. We utilized sequence data from two mitochondrial (cytochrome c oxidase subunit I, 16S rDNA) markers and one nuclear (internal transcribed spacer) marker to examine the genetic diversity of Hermodice throughout its distribution range in the Atlantic Ocean, including the Mediterranean Sea, the Caribbean Sea, the Gulf of Mexico and the Gulf of Guinea. Our analyses revealed generally low genetic divergences among collecting localities and between the East and West Atlantic, although phylogenetic trees based on mitochondrial data indicate the presence of a private lineage in the Mediterranean Sea. A re‐evaluation of the number of branchial filaments confirmed differences between East and West Atlantic populations; however, the differences were not diagnostic and did not reflect the observed genetic population structure. Rather, we suspect that the number of branchial filaments is a function of oxygen saturation in the environment. Our results do not support the distinction between H. carunculata in the West Atlantic and H. nigrolineata in the East Atlantic. Instead, they re‐affirm the older notion that H. carunculata is a cohesive species with a broad distribution across the Atlantic Ocean.     

A polyphasic approach to the delimitation of diatom species: a case study for the genus Pinnularia (Bacillariophyta)

Diatoms are one of the most abundant and arguably the most species‐rich group of protists. Diatom species delimitation has often been based exclusively on the recognition of morphological discontinuities without investigation of other lines of evidence. Even though DNA sequences and reproductive experiments have revealed several examples of (pseudo)cryptic diversity, our understanding of diatom species boundaries and diversity remains limited. The cosmopolitan pennate raphid diatom genus Pinnularia represents one of the most taxon‐rich diatom genera. In this study, we focused on the delimitation of species in one of the major clades of the genus, the Pinnularia subgibba group, based on 105 strains from a worldwide origin. We compared genetic distances between the sequences of seven molecular markers and selected the most variable pair, the mitochondrial cox1 and nuclear encoded LSU rDNA, to formulate a primary hypothesis on the species limits using three single‐locus automated species delimitation methods. We compared the DNA‐based primary hypotheses with morphology and with other available lines of evidence. The results indicate that our data set comprised 15 species of the P. subgibba group. The vast majority of these taxa have an uncertain taxonomic identity, suggesting that several may be unknown to science and/or members of (pseudo)cryptic species complexes within the P. subgibba group.     

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.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.