Testing phylogenetic hypotheses for reconstructing the evolutionary history of Dolichopoda cave crickets in the eastern Mediterranean

Aim To investigate the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the genus Dolichopoda (Orthoptera, Rhaphidophoridae). Location Caves in continental and insular Greece. Methods We sequenced 1967 bp of mitochondrial DNA, corresponding to three fragments of the small and large subunit of the ribosomal RNA (16S and 12S rRNA, respectively) and to the subunit I of cytochrome oxidase (COI), to reconstruct phylogenetic relationships among all 30 known Greek species of Dolichopoda. Alternative hypotheses about the colonization of the Hellenic Peninsula by Dolichopoda species were tested by comparing the degree of discordance between species trees and gene trees under four plausible biogeographical scenarios. Results The present study revealed a rather well resolved phylogeny at species level, identifying a number of clades that represent long‐separated lineages and diverse evolutionary histories within the genus Dolichopoda. Two main clades were revealed within Hellenic–Aegean species, identifying a north‐western and a south‐eastern species group. Based on Bayesian analysis, we applied a relaxed molecular clock to estimate the divergence times between the lineages. The results revealed that the origins of eastern Mediterranean lineages are much older than those of previously studied western Mediterranean Dolichopoda. Tests of alternative biogeographical hypotheses showed that a double colonization of the Hellenic Peninsula, following separate continental and trans‐Aegean routes during the Messinian stage, best accounts for the present distribution of Greek Dolichopoda species. Main conclusions Reconstruction and biogeographical hypothesis testing indicated that the colonization of Greece by Dolichopoda species comprised two episodes and two different routes. The southern lineage probably arose from a trans‐Aegean colonization during the Messinian salinity crisis (5.96–5.33 Ma). The northern lineage could be the result of dispersal from the north through the Balkan Peninsula. The opening of the Mid‐Aegean Trench could have promoted an initial diversification within the uprising Anatolian Plateau, while the Messinian marine regression offered the conditions for a rapid dispersal through the whole Aegean–Hellenic region. In addition, climatic events during the Plio‐Pleistocene may have been responsible for the speciation within each of the two different phylogeographical units, principally attributable to vicariance events.     

Inference of evolutionary patterns of the land snail Albinaria in the Aegean archipelago: is vicariance enough?

Mitochondrial DNA sequences from 16S rRNA and ATPase8 genes were used to investigate phylogeographic patterns of the land snail Albinaria (Gastropoda: Clausiliidae) in the Aegean archipelago. Forty-two populations of Albinaria were analyzed, mainly A. turrita, A. caerulea and A. brevicollis, collected from 22 Aegean islands and certain surrounding regions. Maximum parsimony, maximum likelihood and Bayesian analyses on 16S rRNA and combined datasets produced trees that share significant similarity and reveal a phylogeny with distinct branches which are in general, but not full, agreement with current taxonomy. The Aegean taxa are not monophyletic as a whole, since A. turrita does not cluster with A. caerulea and A. brevicollis. The latter form a distinct monophyletic cluster, within which two groups are evident. These groups do not readily correspond to currently accepted morphospecies; one contains the populations that inhabit the central part of the archipelago plus some eastern islands, while the other contains populations whose geographic distribution is restricted to the southeastern part of the archipelago. The divergence between these two groups is attributed to vicariance events that primarily shape contemporary distributions. Although dispersal may also be present, certain small- and large-scale vicariance events can be traced; alternative phylogeographic hypotheses are discussed in view of the historical biogeography of the region.     

Phylogeography of the Hypnea musciformis species complex (Gigartinales,Rhodophyta) with the recognition of cryptic species in the western Atlantic Ocean

Populations of the marine benthic red macroalgae Hypnea musciformis and Hypnea pseudomusciformis along the Atlantic and Pacific Oceans were tested for phylogeographic structure using the DNA barcode COI‐5P combined with rbcL for the construction of the phylogenetic tree. Strong patterns of genetic structure were detected across 210 COI‐5P DNA sequences, and 37 COI‐5P haplotypes were found, using multiple statistical approaches. Hypnea musciformis was found in the Northeast and Northwest Atlantic, the Mediterrean Sea, Namibia, and along the Pacific coast of Mexico. Two new putative species were detected, Hypnea sp. 1 in the Caribbean Sea and Hypnea sp. 2 in the Dominican Republic. Three distinct marine phylogeographic provinces were recognized in the Southern Hemisphere for H. pseudomusciformis: Uruguay, South‐Southeast Brazil, and Northeast Brazil. The degree of genetic isolation and distinctness among these provinces varied considerably. The Uruguay province was the most genetically distinct, as characterized by four unique haplotypes not shared with any of the Brazilian populations. Statistically significant results support both, isolation by distance and isolation by environment hypotheses, explaining the formation and mantainance of phylogeographic structuring along the Uruguay‐Brazil coast. Geographic, taxonomic and molecular marker concordances were found between our H. pseudomusciformis results and published studies. Furthermore, our data indicate that the Hawaiian introduced populations of H. musciformis contain Hypnea sp. 1 haplotypes, the current known distribution of which is restricted to the Caribbean.     

A multi‐locus approach resolves the phylogenetic relationships of the Simulium asakoae and Simulium ceylonicum species groups in Malaysia: evidence for distinct evolutionary lineages

A multi‐locus approach was used to examine the DNA sequences of 10 nominal species of blackfly in the Simulium subgenus Gomphostilbia (Diptera: Simuliidae) in Malaysia. Molecular data were acquired from partial DNA sequences of the mitochondria‐encoded cytochrome c oxidase subunit I (COI), 12S rRNA and 16S rRNA genes, and the nuclear‐encoded 18S rRNA and 28S rRNA genes. No single gene, nor the concatenated gene set, resolved all species or all relationships. However, all morphologically established species were supported by at least one gene. The multi‐locus sequence analysis revealed two distinct evolutionary lineages, conforming to the morphotaxonomically recognized Simulium asakoae and Simulium ceylonicum species groups.     

Global phylogeography of a cryptic copepod species complex and reproductive isolation between genetically proximate "populations"

The copepod Eurytemora affinis has a broad geographic range within the Northern Hemisphere, inhabiting coastal regions of North America, Asia, and Europe. A phylogenetic approach was used to determine levels of genetic differentiation among populations of this species, and interpopulation crosses were performed to determine reproductive compatibility. DNA sequences from two mitochondrial genes, large subunit (16S) rRNA (450 bp) and cytochrome oxidase I (COI, 652 bp), were obtained from 38 populations spanning most of the species range and from two congeneric species, E. americana and E. herdmani. Phylogenetic analysis revealed a polytomy of highly divergent clades with maximum sequence divergences of 10% in 16S rRNA and 19% in COI. A power test (difference of a proportion) revealed that amount of sequence data collected was sufficient for resolving speciation events occurring at intervals greater than 300,000 years, but insufficient for determining whether speciation events were approximately simultaneous. Geographic and genetic distances were not correlated (Mantel's test; r = 0.023, P = 0.25), suggesting that populations had not differentiated through gradual isolation by distance. At finer spatial scales, there was almost no sharing of mtDNA haplotypes among proximate populations, indicating little genetic exchange even between nearby sites. Interpopulation crosses demonstrated reproductive incompatibility among genetically distinct populations, including those that were sympatric. Most notably, two geographically distant (4000 km) but genetically proximate (0.96% 16S, 0.15% COI) populations exhibited asymmetric reproductive isolation at the F2 generation. Large genetic divergences and reproductive isolation indicate that the morphologically conservative E. affinis constitutes a sibling species complex. Reproductive isolation between genetically proximate populations underscores the importance of using multiple measures to examine patterns of speciation.     

Contrasting geographic structure in evolutionarily divergent Lake Tanganyika catfishes

Geographic isolation is suggested to be among the most important processes in the generation of cichlid fish diversity in East Africa's Great Lakes, both through isolation by distance and fluctuating connectivity caused by changing lake levels. However, even broad scale phylogeographic patterns are currently unknown in many non‐cichlid littoral taxa from these systems. To begin to address this, we generated restriction‐site‐associated DNA sequence (RADseq) data to investigate phylogeographic structure throughout Lake Tanganyika (LT) in two broadly sympatric rocky shore catfish species from independent evolutionary radiations with differing behaviors: the mouthbrooding claroteine, Lophiobagrus cyclurus, and the brood‐parasite mochokid, Synodontis multipunctatus. Our results indicated contrasting patterns between these species, with strong lake‐wide phylogeographic signal in L. cyclurus including a deep divergence between the northern and southern lake basins. Further structuring of these populations was observed across a heterogeneous habitat over much smaller distances. Strong population growth was observed in L. cyclurus sampled from shallow shorelines, suggesting population growth associated with the colonization of new habitats following lake‐level rises. Conversely, S. multipunctatus, which occupies a broader depth range, showed little phylogeographic structure and lower rates of population growth. Our findings suggest that isolation by distance and/or habitat barriers may play a role in the divergence of non‐cichlid fishes in LT, but this effect varies by species.     

Genetic study reveals close link between Irish and Northern Spanish specimens of the protected Lusitanian slug Geomalacus maculosus

The slug Geomalacus maculosus is a prominent member of the Lusitanian fauna. As its global distribution is restricted to western Ireland and northern Iberia, it is protected under EU legislation. Nothing is known about the genetic variability and population structure of this species, so, with a special view to shedding light on the origin of the Irish G. maculosus, tissue samples from 78 specimens were collected from 13 locations within Ireland and ten locations within Iberia and partial sequences of the mitochondrial 16S rRNA and cytochrome oxidase subunit 1 (COI) and from the nuclear internal transcribed spacer 1 region (ITS‐1) were compared. The genetic diversity of the Irish G. maculosus was found to be greatly reduced compared with the Iberian populations, with only one (16S rRNA) and two (COI) mitochondrial haplotypes identified respectively. No private Irish ITS‐1 haplotype was found. Based on the COI sequences, the Irish specimens clustered closest to Spanish specimens from Northern Asturias and Cantabria, and the bGMYC analysis identified five further Iberian clades that were highly genetically differentiated suggesting long‐term allopatric divergence.     

Polymorphism of mtDNA regions in brown hare (Lepus europaeus) populations from Vojvodina (Serbia and Montenegro)

It is well known that there is heterogeneity of mitochondrial DNA (mtDNA) within and among natural populations of same species. The polymorphism level of particular regions of mtDNA gives valuable results in detection of population genetic structure. The aim of this paper was to detect polymorphism of three mtDNA regions: cytochrome oxidase I (COI), Control region, and 12S/16S rRNA, by the mtDNA RFLP-PCR method, in three Lepus europaeus populations from Vojvodina province (Serbia and Montenegro). Polymorphism was detected within the two regions, COI and Control region, while 12S/16S rRNA region was monomorphic in all 77 individuals. Eight haplotypes were detected in the brown hare population in Vojvodina, and three were unique for the Srem brown hare population.     

Molecular phylogeny of the Greek populations of the genus Ligidium (Isopoda, Oniscidea) using three mtDNA gene segments

The phylogeny of Greek populations of the terrestrial isopod genus Ligidium is reconstructed based on three mtDNA gene segments: 12S rRNA, 16S rRNA and COI. Two widely distributed European species, as well as three outgroups belonging to different isopod genera, were also included in the analyses. The samples used represent almost all Ligidium species known to occur in Greece, as well as several populations of unknown specific status plus some new records. Phylogenetic analyses of the combined data set were performed using Bayesian inference and maximum parsimony. The two main sister clades with good support indicate the sympatric differentiation of two lineages in southern continental Greece (Peloponnisos), where Ligidium populations exhibit a mosaic distribution of sibling species. The insular populations of the Aegean Islands show increased genetic divergence and form separate clades. The presence of a third lineage of Asiatic origin is strongly suggested by both the molecular phylogeny and morphology. The only presumably valid diagnostic morphological character exhibits only partial correspondence to well supported clades of the molecular phylogeny. Genetic differentiation between populations is very high, a fact that can be attributed to the strict ecological specialization of these animals that leads to increased levels of isolation even between populations that are in close proximity. As a consequence, Greek Ligidium populations, especially those present on islands, are unique genetic pools and extremely vulnerable to extinction.     

Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians

BACKGROUND: Identifying species of organisms by short sequences of DNA has been in the center of ongoing discussions under the terms DNA barcoding or DNA taxonomy. A C-terminal fragment of the mitochondrial gene for cytochrome oxidase subunit I (COI) has been proposed as universal marker for this purpose among animals. RESULTS: Herein we present experimental evidence that the mitochondrial 16S rRNA gene fulfills the requirements for a universal DNA barcoding marker in amphibians. In terms of universality of priming sites and identification of major vertebrate clades the studied 16S fragment is superior to COI. Amplification success was 100% for 16S in a subset of fresh and well-preserved samples of Madagascan frogs, while various combination of COI primers had lower success rates.COI priming sites showed high variability among amphibians both at the level of groups and closely related species, whereas 16S priming sites were highly conserved among vertebrates. Interspecific pairwise 16S divergences in a test group of Madagascan frogs were at a level suitable for assignment of larval stages to species (1-17%), with low degrees of pairwise haplotype divergence within populations (0-1%). CONCLUSION: We strongly advocate the use of 16S rRNA as standard DNA barcoding marker for vertebrates to complement COI, especially if samples a priori could belong to various phylogenetically distant taxa and false negatives would constitute a major problem.     

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

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

Comparative phylogeography of five widespread tree species: Insights into the history of western Amazonia

Various historical processes have been put forth as drivers of patterns in the spatial distribution of Amazonian trees and their population genetic variation. We tested whether five widespread tree species show congruent phylogeographic breaks and similar patterns of demographic expansion, which could be related to proposed Pleistocene refugia or the presence of geological arches in western Amazonia. We sampled Otoba parvifolia/glycycarpa (Myristicaceae), Clarisia biflora, Poulsenia armata, Ficus insipida (all Moraceae), and Jacaratia digitata (Caricaceae) across the western Amazon Basin. Plastid DNA (trnH–psbA; 674 individuals from 34 populations) and nuclear ribosomal internal transcribed spacers (ITS; 214 individuals from 30 populations) were sequenced to assess genetic diversity, genetic differentiation, population genetic structure, and demographic patterns. Overall genetic diversity for both markers varied among species, with higher values in populations of shade‐tolerant species than in pioneer species. Spatial analysis of molecular variance (SAMOVA) identified three genetically differentiated groups for the plastid marker for each species, but the areas of genetic differentiation were not concordant among species. Fewer SAMOVA groups were found for ITS, with no detectable genetic differentiation among populations in pioneers. The lack of spatially congruent phylogeographic breaks across species suggests no common biogeographic history of these Amazonian tree species. The idiosyncratic phylogeographic patterns of species could be due instead to species‐specific responses to geological and climatic changes. Population genetic patterns were similar among species with similar biological features, indicating that the ecological characteristics of species impact large‐scale phylogeography.     

Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea,inferred from mitochondria DNA markers

The Asian green mussel Perna viridis is ecologically and economically important in the coastal regions of China. In order to characterize the genetic diversity and population connectivity of P. viridis in South China Sea, a 664 bp region of mitochondrial COI gene and a 293 bp region of 16S rRNA gene were sequenced and analyzed for 78 and 92 individuals from four populations in South China Sea, respectively. A total of 15 haplotypes were defined by 14 variable nucleotide sites in COI gene, and 7 haplotypes by 6 variable nucleotide sites in 16S rRNA gene. High haplotype diversity and low nucleotide diversity were observed in COI gene, while moderate haplotype diversity and low nucleotide diversity were observed in 16S rRNA gene. Pairwise F_ST values of COI gene were all negative and those of 16S rRNA gene ranged from −0.01409 to 0.10289. The results showed that no significant genetic divergence (or shallow genetic structure) and high levels of population connectivity among the four populations of P. viridis in South China Sea.     

Genetic variation in populations of the endangered fish Ladigesocypris ghigii and its implications for conservation

1. The genetic variation of the endangered freshwater fish Ladigesocypris ghigii, endemic to the island of Rhodes (Greece), was investigated for nine populations, originating from seven different stream systems and a reservoir, both at the mtDNA and nuclear level, in order to suggest conservation actions. 2. Both restriction fragment length polymorphism analysis of five segments of mitochondrial DNA (ND‐5/6, COI and 12S‐16S rRNA) amplified by polymerase chain reaction, and random amplified polymorphic DNA analysis, revealed extremely low levels of intra‐population polymorphism. It is highly likely that the low intra‐population variability is the result of successive bottleneck events evident in shrinkage and expansion of the populations year after year, which may have led to a complete loss of several genotypes and haplotypes, and an increased degree of inbreeding. 3. Inter‐population genetic structuring was high, with fixation of haplotypes within six of the nine populations and fixation of alleles within populations originating from different waterbodies. It is probable that all haplotypes and/or alleles found were initially represented in all populations. However, because of the long time of isolation coupled with successive bottleneck and subsequent genetic drift, common mtDNA haplotypes and alleles among the populations may have become rare or extinct through stochastic lineage loss. 4. Although nucleotide divergence among haplotypes was very shallow, half of the haplotypes recorded (three of six), resulted from nucleotide changes on the 12S–16S rRNA segments, which are the most conserved part of the mitochondrial genome. This fact may indicate that the observed genetic variation did not necessarily result only from the retention of ancestral polymorphism, but may have arisen through mutation and complete lineage sorting over a relatively small number of generations, once the populations had become isolated from one another. 5. Our data suggest that two of the L. ghigii populations may be on independent evolutionary trajectories. Considering that each population appears so far well adapted within each site, all populations should be managed and conserved separately.     

Phylogeography of the Poecilimon luschani species group (Orthoptera,Tettigoniidae): a radiation strictly correlated with climatic transitions in the Pleistocene

The Quaternary biogeography of Anatolia has received considerable interest recently. Here, the genealogical history of the Anatolio?Balkan lineage of the Poecilimon luschani species group was evaluated. Using concatenated data from 16S rDNA and cytochrome c oxidase subunit I (COI) sequences, the timings of inter‐ and intraspecies radiations were estimated. The demographic history of the populations was estimated using a data set established from COI sequences. Genetic diversity was very high in almost all of the populations studied. Fixation indices suggested extreme divergence of P. luschani. A molecular chronogram estimated a radiation history for the species/subspecies over a period ranging from 1.323 to 0.440 Myr. Demographic analyses applied to 11 populations suggested departure in population size for most of the local populations. The following conclusions were reached: (1) P. luschani originated from an Anatolio‐Aegean ancestral stock and extended its range to the Balkans through Dardanelles during the Early Pleistocene; (2) the Mid‐Pleistocene Transition, the lengthening of glacial period from 41 to 100 Kyr and the initiation of intense glaciation periods are the three main events corresponding to the main nodes of the chronogram; (3) altitudinal heterogeneity played a buffer role during the glacial cycles, allowing populations to cope with severe environmental changes; (4) the effects of Pleistocene climate cycles on populations differ according to altitudinal and latitudinal location in Anatolia, and (5) habitat preferences, such as altitudinal range, may easily shift because of changes in environmental conditions. © 2014 The Linnean Society of London     

Genetic variation in the stygobitic shrimp Creaseria morleyi (Decapoda: Palaemonidae), evidence of bottlenecks and re‐invasions in the Yucatan Peninsula

The karstic nature of the Yucatan Peninsula promotes the formation of submerged caves and sinkholes that are inhabited by an endemic subterranean water fauna. By contrast with most other micro‐endemic stygobitic species, the freshwater palaemonid shrimp Creaseria morleyi is widely distributed across the northern part of the peninsula. In the present study, we investigated the phylogeographic structure of C. morleyi using two mitochondrial genes as markers, and explored hypotheses related to its evolution in the peninsula. DNA from 14 localities was extracted, and the 16S rRNA and cytochrome oxidase subunit I (COI) genes were amplified and sequenced. The different haplotypes were identified to construct a haplotype network and perform a nested clade analysis. Five haplotypes of the 16S gene were obtained, with a maximum divergence of 0.5%. One of these haplotypes is widely distributed and the most divergent is located in the north‐western section of the peninsula. Twelve haplotypes for the COI gene were found with a maximum divergence of 2%, showing the same spatial pattern. The analysis revealed two significantly different clades corresponding to populations in the centre and south‐east of the peninsula as a consequence of restricted genetic flow with isolation‐by‐distance. The divergence time of these two clades was 40–120 thousand years. The genetic variation in C. morleyi, the relationship between haplotypes and their geographic distribution, along with the geological history of the Yucatan Peninsula, may indicate that this variation is a relict of an ancient marked genetic structure reduced by changes in sea level that resulted in a series of bottlenecks. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 315–325.     

Identifying Pleistocene refugia in North American cold deserts using phylogeographic analyses and ecological niche modelling

Aim We analysed phylogeographic patterns and ecological niche models (ENMs) of the widespread velvet ant (Hymenoptera: Mutillidae) Sphaeorpthalma difficilis to understand the history of diversification in the Nearctic deserts and to identify areas that may have been cold‐desert refugia during the Pleistocene. These areas should be targeted for conservation because of their climatic stability throughout historical climate change events. Location North American arid regions. Methods The two internal transcribed spacer regions (ITS1 and ITS2) were sequenced and analysed using Bayesian techniques to uncover phylogeographic patterns of relatedness among S. difficilis populations. History of diversification was estimated using parsimony‐based and maximum likelihood character reconstructions. Molecular dating analyses were implemented in the program r8s and were calibrated with Dominican amber fossils. ENMs were developed based on current climate data and projected onto Pleistocene climate surfaces. Results The analyses suggest that S. difficilis had a complex history of Pleistocene range expansion and contraction that led to the formation of genetically distinct populations inhabiting distinct arid regions. ENMs and phylogeographic patterns indicate that several cold‐desert refugia existed in North America, particularly in the Colorado Plateau and parts of the Great Basin Desert. Main conclusions Analyses of S. difficilis are used to identify potential Pleistocene refugia in the North American cold deserts. Because these areas represent climatically stable locations, they are critical for the long‐term persistence of biodiversity. This research provides evidence that in addition to desert‐like conditions persisting through the ice age in parts of the Nearctic warm deserts, many areas maintained desert‐like characteristics in the regional cold deserts. Further work is needed to elucidate options for preserving biodiversity in these cold‐desert refugia.     

Testing for shared biogeographic history in the lower Central American freshwater fish assemblage using comparative phylogeography: concerted,independent, or multiple evolutionary responses?

A central goal of comparative phylogeography is determining whether codistributed species experienced (1) concerted evolutionary responses to past geological and climatic events, indicated by congruent spatial and temporal patterns (“concerted‐response hypothesis”); (2) independent responses, indicated by spatial incongruence (“independent‐response hypothesis”); or (3) multiple responses (“multiple‐response hypothesis”), indicated by spatial congruence but temporal incongruence (“pseudocongruence”) or spatial and temporal incongruence (“pseudoincongruence”). We tested these competing hypotheses using DNA sequence data from three livebearing fish species codistributed in the Nicaraguan depression of Central America (Alfaro cultratus, Poecilia gillii, and Xenophallus umbratilis) that we predicted might display congruent responses due to co‐occurrence in identical freshwater drainages. Spatial analyses recovered different subdivisions of genetic structure for each species, despite shared finer‐scale breaks in northwestern Costa Rica (also supported by phylogenetic results). Isolation‐with‐migration models estimated incongruent timelines of among‐region divergences, with A. cultratus and Xenophallus populations diverging over Miocene–mid‐Pleistocene while P. gillii populations diverged over mid‐late Pleistocene. Approximate Bayesian computation also lent substantial support to multiple discrete divergences over a model of simultaneous divergence across shared spatial breaks (e.g., Bayes factor [B₁₀] = 4.303 for Ψ [no. of divergences] > 1 vs. Ψ = 1). Thus, the data support phylogeographic pseudoincongruence consistent with the multiple‐response hypothesis. Model comparisons also indicated incongruence in historical demography, for example, support for intraspecific late Pleistocene population growth was unique to P. gillii, despite evidence for finer‐scale population expansions in the other taxa. Empirical tests for phylogeographic congruence indicate that multiple evolutionary responses to historical events have shaped the population structure of freshwater species codistributed within the complex landscapes in/around the Nicaraguan depression. Recent community assembly through different routes (i.e., different past distributions or colonization routes), and intrinsic ecological differences among species, has likely contributed to the unique phylogeographical patterns displayed by these Neotropical fishes.     

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA‐based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR‐amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR‐amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR‐amplified the blend using five primer sets, targeting either COI or 16S, with high‐throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.     

Phylogeographic patterns of Merodon hoverflies in the Eastern Mediterranean region: revealing connections and barriers

We investigated the phylogeographic patterns of Merodon species (Diptera, Syrphidae) in the Eastern Mediterranean. Ten species were sampled on five different islands and mainland sites as a minimum. All samples were screened for their mtDNA COI barcode haplotype diversity, and for some samples, we additionally generated genomic fingerprints. The recently established zoogeographic distribution categories classify these species as having (1) Balkan distribution; (2) Anatolian distribution; (3) continental areas and large islands distribution; and (4) with wide distribution. The ancestral haplotypes and their geographical localities were estimated with statistical parsimony (TCS). TCS networks identified as the ancestral haplotype samples that originated from localities situated within the distributional category of the species in question. Strong geographical haplotype structuring was detected for many Merodon species. We were particularly interested to test the relative importance of current (Aegean Sea) and past Mid‐Aegean Trench) barriers to dispersal for Merodon flies in the Aegean. We employed phylogenetic β‐diversity (Pβ_total) and its partition in replacement (Pβ_repl) and richness difference (Pβ_rich) to test the importance of each explanatory variable (interisland distance, MAT, and island area) in interisland differences using partial Mantel tests and hierarchical partitioning of variation. β‐Analyses confirmed the importance of both current and past barriers to dispersal on the evolution of group. Current interisland distance was particularly important to explain the replacement of haplotypes, while the MAT was driving differences in richness of haplotypes, revealing the MAT as a strong past barrier whose effects are still visible today in the phylogenetic history of the clade in the Aegean. These results support the hypothesis of a highly restricted dispersal and gene flow among Merodon populations between islands since late Pleistocene. Additionally, patterns of phylogeographic structure deduced from haplotype connections and ISSR genome fingerprinting data revealed a few putative cases of human‐mediated transfers of Merodon spp.