Multilocus phylogeny and systematics of Iberian endemic Squalius (Actinopterygii,Leuciscidae)

Inferring the evolutionary history of a group of species can be challenging given the many factors involved. In recent years, the increased availability of sequences of multiple genes per species has spurred the development of new methodologies to analyse multilocus data sets. Two approaches that analyse such data are concatenated supermatrix and coalescent-based species-tree analyses. In this study, we used both of these methods to infer the phylogenetic relationships of Iberian species of the genus Squalius from one mitochondrial and six nuclear genes. We found mitonuclear discordance in the phylogenetic relationships of the group. According to the mitochondrial gene analysis, all species were recovered as monophyletic except S. pyrenaicus; besides, in the concatenated supermatrix analysis of the nuclear markers, this species resolved as polyphyletic with three divergent evolutionary lineages. The coalescent-based nuclear species-tree analysis rendered a well-resolved phylogeny compared with the supermatrix analysis, which was unable to discern between S. carolitertii, S. castellanus and one of the evolutionary lineages of S. pyrenaicus. This result is likely due to the better integration of population uncertainty in the coalescent approach. Furthermore, Bayesian multilocus species delimitation analyses based on a BPP approach strongly supported the distinct nuclear lineages as different species. Nevertheless, the supermatrix analysis was able to obtain well-supported relationships in the divergent lineages with low numbers of individuals. Our study highlights the usefulness of different analytical methodologies to obtain a more complete picture of the evolutionary history of taxa, especially when discordant patterns among genes are found.     

Phylogeny and species delimitation of North European Lumbricillus (Clitellata,Enchytraeidae)

The enchytraeid genus Lumbricillus comprises about 80 described species of clitellate worms, which are up to a few centimetres long, and they mostly inhabit the littoral zone of non‐tropical marine and brackish waters world‐wide. The phylogeny of this genus is poorly studied, but previous work has suggested that Lumbricillus is a non‐monophyletic group. In this study, species boundaries and the phylogeny of this genus is re‐assessed using more than 300 DNA‐barcoded specimens (using COI mtDNA), part of which was also sequenced for two additional mitochondrial and four nuclear molecular markers. Statistical and coalescent based applications were used for the delimitation of a total of 24 species, of which 20 were identified as belonging to 17 described morphospecies; one morphospecies was found to be a complex of four delimited species, and another four delimited species could not be matched with any described species. Furthermore, gene trees, concatenation and multispecies coalescent based species trees were estimated using Bayesian inference. The estimated phylogenies confirm a non‐monophyletic Lumbricillus as L. semifuscus is clearly excluded from the genus. Furthermore, the placement of a monophyletic clade consisting of L. arenarius, L. dubius, and an unidentified species varies between analyses; they are either found as the sister‐group to the genus Grania or as sister‐group to the remaining Lumbricillus, where the latter relationship is supported by the multispecies coalescent, which we consider as the most reliable method.     

Phylogenetic incongruence between nuclear and mitochondrial markers in the Asian colobines and the evolution of the langurs and leaf monkeys

Evidence of incongruence between mitochondrial and nuclear gene trees is now becoming documented with increasing frequency. Among the Old World monkeys, this discordance has been well demonstrated in the Cercopithecinae, but has not yet been investigated in the Colobinae. The mitochondrial relationships between the colobine genera have recently been clarified and cluster Presbytis and Trachypithecus as sister taxa to the exclusion of Semnopithecus. This is incongruent with previous morphological hypotheses that suggest the latter two are sister taxa, and perhaps even congeneric. In addition to analyzing a previously published 10,896 bp mitochondrial dataset, we sequenced and analyzed a 4297 bp fragment of the X-chromosome in order to test the competing mitochondrial and morphological phylogenetic hypotheses. The results from the mitochondrial dataset again support a Presbytis+Trachypithecus group while the X-chromosomal dataset supported a Semnopithecus+Trachypithecus group. A Shimodaira-Hasegawa test performed on both datasets indicates that the mitochondrial and X-chromosomal trees are significantly better at explaining their respective datasets than alternative topologies (p<0.05). We suggest that differential lineage sorting or ancient hybridization may be the cause of this strong discordance between the mitochondrial and X-chromosomal markers in these taxa.     

Cracking the nut: geographical adjacency of sister taxa supports vicariance in a polytomic salamander clade in the absence of node support

The urodelan genus Lyciasalamandra, which inhabits a relatively small area along the southern Turkish coast and some Aegean islands, provides an outstanding example of a diverse but phylogenetically unresolved taxon. Molecular trees contain a single basal polytomy that could be either soft or hard. We here use the information of nuclear (allozymes) and mitochondrial (fractions of the 16S rRNA and ATPase genes) datasets in combination with area relationships of lineages to resolve the phylogenetic relationships among Lyciasalamandra species in the absence of sufficient node support. We can show that neither random processes nor introgressive hybridization can be invoked to explain that the majority of pairs of sister taxa form geographically adjacent units and interpret that this pattern has been shaped by vicariant events. Topology discordance between mitochondrial and nuclear trees mainly refers to an affiliation of L. helverseni, a taxon restricted to the Karpathos archipelago, to the western-most and geographically proximate mainland taxon in the nuclear tree, while in the organelle tree it turns out to be the sister lineage to the geographically most distant eastern clade. As this discordance cannot be explained by long-branch attraction in either dataset we suppose that oversea dispersal may have accounted for a second colonization of the Karpathos archipelago. It may have initiated introgression and selection driven manifestation of alien eastern mitochondrial genomes on a western nuclear background. Our approach of testing for area relationships of sister taxa against the null hypothesis of random distribution of these taxa seems to be especially helpful in phylogenetic studies where traditional measures of phylogenetic branch support fail to reject the null hypothesis of a hard polytomy.     

Biome specificity of distinct genetic lineages within the four-striped mouse Rhabdomys pumilio (Rodentia: Muridae) from southern Africa with implications for taxonomy

Within southern Africa, a link between past climatic changes and faunal diversification has been hypothesized for a diversity of taxa. To test the hypothesis that evolutionary divergences may be correlated to vegetation changes (induced by changes in climate), we selected the widely distributed four-striped mouse, Rhabdomys, as a model. Two species are currently recognized, the mesic-adapted R. dilectus and arid-adapted R. pumilio. However, the morphology-based taxonomy and the distribution boundaries of previously described subspecies remain poorly defined. The current study, which spans seven biomes, focuses on the spatial genetic structure of the arid-adapted R. pumilio (521 specimens from 31 localities), but also includes limited sampling of the mesic-adapted R. dilectus (33 specimens from 10 localities) to act as a reference for interspecific variation within the genus. The mitochondrial COI gene and four nuclear introns (Eef1a1, MGF, SPTBN1, Bfib7) were used for the construction of gene trees. Mitochondrial DNA analyses indicate that Rhabdomys consists of four reciprocally monophyletic, geographically structured clades, with three distinct lineages present within the arid-adapted R. pumilio. These monophyletic lineages differ by at least 7.9% (±0.3) and these results are partly confirmed by a multilocus network of the combined nuclear intron dataset. Ecological niche modeling in MaxEnt supports a strong correlation between regional biomes and the distribution of distinct evolutionary lineages of Rhabdomys. A Bayesian relaxed molecular clock suggests that the geographic clades diverged between 3.09 and 4.30Ma, supporting the hypothesis that the radiation within the genus coincides with paleoclimatic changes (and the establishment of the biomes) characterizing the Miocene-Pliocene boundary. Marked genetic divergence at the mitochondrial DNA level, coupled with strong nuclear and mtDNA signals of non-monophyly of R. pumilio, support the notion that a taxonomic revision of the genus is needed.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

THE INTERFACE BETWEEN PHYLOGENETICS AND POPULATION GENETICS: INVESTIGATING GENE TREES,SPECIES TREES,AND POPULATION DYNAMICS IN THE PHYLLOPHAGA FRATERNA SPECIES GROUP

This study uses traditional and contemporary phylogenetic and population genetic analyses to assess the causes of discordance (i.e., lineage sorting and introgression) among mitochondrial and nuclear gene trees for a clade of eastern North American scarab beetles (fraterna species group, genus Phyllophaga). I estimated gene trees using individual and combined analysis of one mitochondrial and two nuclear loci in MrBayes , and inferred a species tree using a hierarchical coalescent approach based on all loci in the program Best . Because hybridization violates the assumptions of Best , I tested for introgression by comparing species monophyly between the mitochondrial and nuclear gene trees based on the prediction that cytoplasmic genomes introgress more readily than nuclear genomes. Haplotype exclusivity was identified using Bayesian tests of monophyly and the genealogical sorting index. I used the results of the phylogenetic analyses and monophyly tests to develop an explicit hypothesis of introgression that could be tested in the program IMa. Results from these analyses provided evidence for introgression across clades within the fraterna group. The tiered analytical approach used in this study demonstrated how the use of multiple methods can identify when assumptions are violated and methods are prone to yield misleading results.     

Molecular systematics of the neotropical shovelnose catfish genus Pseudoplatystoma Bleeker 1862 based on nuclear and mtDNA markers

Pseudoplatystoma is a commercially important genus of Neotropical migratory catfishes widely distributed in all major river basins of South America. Historically, only three species were recognized, but a recent revision proposed eight putative morphospecies for the genus. A molecular study based on mitochondria DNA (mtDNA) provided support for recognition of only some of the species and raised questions about species boundaries in this group. We present a more encompassing analysis based on mtDNA (cytochrome b, 818bp) and nuclear DNA-based phylogenies (Rag1 intron 1, 664bp and S7 intron 1, 635bp) for a more extensive sampling (279 individuals from 42 localities) of all putative species in all major river basins. Patterns generated by individual gene genealogies and a multispecies coalescent analysis provided evidence to suggest recognition of only four distinct species in this genus: Pseudoplatystoma magdaleniatum, Pseudoplatystoma corruscans, Pseudoplatystoma tigrimun (sensu lato) and Pseudoplatystoma fasciatum (sensu lato). The species phylogeny places P. magdaleniatum as the sister group to all the other species in the genus, but the relationships among P. fasciatum s.l, P. tigrimum s.l., and P. corruscans could not be resolved with confidence.     

Diversification and biogeographical history of Neotropical plethodontid salamanders

The Neotropical bolitoglossine salamanders represent an impressive adaptive radiation, comprising roughly 40% of global salamander species diversity. Despite decades of morphological studies and molecular work, a robust multilocus phylogenetic hypothesis based on DNA sequence data is lacking for the group. We estimated species trees based on multilocus nuclear and mitochondrial data for all major lineages within the bolitoglossines, and used our new phylogenetic hypothesis to test traditional biogeographical scenarios and hypotheses of morphological evolution in the group. In contrast to previous phylogenies, our results place all Central American endemic genera in a single clade and suggest that Central America played a critical role in the early biogeographical history of the group. The large, predominantly Mexican genus Pseudoeurycea is paraphyletic, and analyses of the nuclear data place two lineages of Pseudoeurycea as the sister group of Bolitoglossa. Our phylogeny reveals extensive homoplasy in morphological characters, which may be the result of truncation or alteration of a shared developmental trajectory. We used our phylogenetic results to revise the taxonomy of the genus Pseudoeurycea. © 2015 The Linnean Society of London     

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.     

Interspecific hybridization causes long‐term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes

Classification, phylogeography and the testing of evolutionary hypotheses rely on correct estimation of species phylogeny. Early molecular phylogenies often relied on mtDNA alone, which acts as a single linkage group with one history. Over the last decade, the use of multiple nuclear sequences has often revealed conflict among gene trees. This observation can be attributed to hybridization, lineage sorting, paralogy or selection. Here, we use 54 groups of fishes from 48 studies to estimate the degree of concordance between mitochondrial and nuclear gene trees in two ecological grades of fishes: marine and freshwater. We test the hypothesis that freshwater fish phylogenies should, on average, show more discordance because of their higher propensity for hybridization in the past. In keeping with this idea, concordance between mitochondrial and nuclear gene trees (as measured by proportion of components shared) is on average 50% higher in marine fishes. We discuss why this difference almost certainly results from introgression caused by greater historical hybridization among lineages in freshwater groups, and further emphasize the need to use multiple nuclear genes, and identify conflict among them, in estimation of species phylogeny.     

Species-level phylogeny of 'Satan's perches' based on discordant gene trees (Teleostei: Cichlidae: Satanoperca Günther 1862)

Neotropical rivers are home to the largest assemblage of freshwater fishes, but little is known about the phylogeny of these fishes at the species level using multi-locus molecular markers. Here, we present a phylogeny for all known species of the genus Satanoperca, a widespread group of Neotropical cichlid fishes, based on analysis of six unlinked genetic loci. To test nominal and proposed species limits for this group, we surveyed mtDNA sequence variation among 320 individuals representing all know species. Most nominal species were supported by this approach but we determined that populations in the Xingu, Tapajós, and Araguaia+Paraná Rivers are likely undescribed species, while S. jurupari and S. mapiritensis did not show clear genetic distinction. To infer a phylogeny of these putative species, we conducted maximum likelihood and Bayesian non-clock and relaxed clock analyses of concatenated data from three genes (one mitochondrial, two nuclear). We also used a multi-species coalescent model to estimate a species tree from six unlinked loci (one mitochondrial, five nuclear). The topologies obtained were congruent with other results, but showed only minimal to moderate support for some nodes, suggesting that more loci will be needed to satisfactorily estimate the distribution of coalescent histories within Satanoperca. We determined that this variation results from topological discordance among separate gene trees, likely due to differential sorting of ancestral polymorphisms.     

Multi-locus phylogenetic analyses support the monophyly and the Miocene diversification of Iksookimia (Teleostei: Cypriniformes: Cobitidae)

The genus Iksookimia contains six species of primary freshwater fishes that are endemic to South Korea. Previous phylogenetic studies, based on DNA sequence data from three or fewer loci, have suggested non-monophyly of the genus, providing inconsistent resolutions of the relationships of Iksookimia. Our coalescent and concatenation-based phylogenetic analyses, utilizing seven unlinked nuclear-encoded genes, strongly supported Iksookimia as a monophyletic group, emphasizing the importance of multi-locus data in investigating complicated phylogenetic relationships. A relaxed molecular clock analysis using fossil calibrations, indicated that the origin of the major lineages of Iksookimia occurred between ～12 to 5 Ma, which is consistent with the Miocene uplift of the Taebaek and Sobaek Mountains and the Miocene activation of the major south-eastern faults. These palaeogeographic events may have served as vicariant events in the diversification of Iksookimia.     

Mitonuclear discordance as a confounding factor in the DNA taxonomy of monogonont rotifers

Discordance between mitochondrial and nuclear phylogenies is being increasingly recognized in animals and may confound DNA‐based taxonomy. This is especially relevant for taxa whose microscopic size often challenges any effort to distinguish between cryptic species without the assistance of molecular data. Regarding mitonuclear discordance, two strikingly contrasting scenarios have been recently demonstrated in the monogonont rotifers of the genus Brachionus. While strict mitonuclear concordance was observed in the marine B. plicatilis species complex, widespread hybridization‐driven mitonuclear discordance was revealed in the freshwater B. calyciflorus species complex. Here, we investigated the frequency of occurrence and the potential drivers of mitonuclear discordance in three additional freshwater monogonont rotifer taxa, and assessed its potential impact on the reliability of DNA taxonomy results based on commonly used single markers. We studied the cryptic species complexes of Keratella cochlearis, Polyarthra dolichoptera and Synchaeta pectinata. Phylogenetic reconstructions were based on the mitochondrial barcoding marker cytochrome c oxidase subunit I gene and the nuclear internal transcribed spacer 1 locus, which currently represent the two most typical genetic markers used in rotifer DNA taxonomy. Species were delimited according to each marker separately using a combination of tree‐based coalescent, distance‐based and allele‐sharing‐based approaches. Mitonuclear discordance was observed in all species complexes with incomplete lineage sorting and unresolved phylogenetic reconstructions recognized as the likely drivers. Evidence from additional sources, such as morphology and ecology, is thus advisable for deciding between often contrasting mitochondrial and nuclear species scenarios in these organisms.     

Considering gene flow when using coalescent methods to delimit lineages of North American pitvipers of the genus Agkistrodon

Examining species diversity and mechanisms of speciation using coalescent models provides a framework for how regional diversity is accrued, even in well‐studied areas such as the Nearctic. It is likely, that gene flow among closely‐related species with adjacent distributions may be common. However, the absence of gene flow is a primary assumption of many phylogeographical methods that produce species trees and delimit species using Bayesian or likelihood functions in a coalescent framework. In the present study, we examine delimitation when gene flow between species is present using empirical datasets from two species of North American pitvipers of the genus Agkistrodon. We also use niche modelling to determine whether these young lineages occur in distinct environmental niches. To manage the problem of gene flow between species, we first identify admixed individuals, demonstrate that gene flow has occurred, and then identify the impact of alternative population assignments of admixed individuals on delimitation posterior probabilities. In addition, we examine the influence of mitochondrial genes relative to other loci combined in coalescent analyses that delimit species. Here, we find that the copperheads (Agkistrodon contortrix) and the cottonmouths (Agkistrodon piscivorus) are each composed of two distinct species, with each occupying different niches. Importantly, we find that species can be delimited when the amount of gene flow between lineages is low, although the methods are acutely sensitive to population assignment of individuals. © 2014 The Linnean Society of London     

Phylogeny of ground beetles subgenus Nialoe (s. lat.) Tanaka (Coleoptera: Carabidae; genus Pterostichus): A molecular phylogenetic approach

We constructed a phylogeny of the ground beetle subgenus Nialoe ( s. lat. ), genus Pterostichus (Coleoptera: Carabidae) based on two mitochondrial (cytochrome oxidase I and 16S ribosomal DNA) and one nuclear (28S ribosomal DNA) gene sequences. Thirty-three representative species of the group and three outgroup species were analyzed. The resultant trees (maximum parsimonious, maximum likelihood and Bayesian trees of the combined data of the three gene sequences) indicated that there are two large and three small lineages in the group, some of which were supported by a previous morphology-based phylogeny. In all the analyses, the small lineage composed of two Korean species is sister to the rest of the subgenus, but relationships of other four lineages differed among the analyses and remained unresolved. The implications of the present results are discussed in terms of taxonomy and biogeography of the group.     

Phylogeny and cryptic diversity in geckos (Phyllopezus; Phyllodactylidae; Gekkota) from South America's open biomes

The gecko genus Phyllopezus occurs across South America's open biomes: Cerrado, Seasonally Dry Tropical Forests (SDTF, including Caatinga), and Chaco. We generated a multi-gene dataset and estimated phylogenetic relationships among described Phyllopezus taxa and related species. We included exemplars from both described Phyllopezus pollicaris subspecies, P. p. pollicaris and P. p.przewalskii. Phylogenies from the concatenated data as well as species trees constructed from individual gene trees were largely congruent. All phylogeny reconstruction methods showed Bogertia lutzae as the sister species of Phyllopezus maranjonensis, rendering Phyllopezus paraphyletic. We synonymized the monotypic genus Bogertia with Phyllopezus to maintain a taxonomy that is isomorphic with phylogenetic history. We recovered multiple, deeply divergent, cryptic lineages within P. pollicaris. These cryptic lineages possessed mtDNA distances equivalent to distances among other gekkotan sister taxa. Described P. pollicaris subspecies are not reciprocally monophyletic and current subspecific taxonomy does not accurately reflect evolutionary relationships among cryptic lineages. We highlight the conservation significance of these results in light of the ongoing habitat loss in South America's open biomes.     

Exon primed intron-crossing (EPIC) markers reveal natural hybridization and introgression in Actinidia (Actinidiaceae) with sympatric distribution

Interspecific hybridization is widespread among plants. Understanding the phylogenetic relationships among species is necessary for revealing the potential hybridization events. Actinidia, best known as kiwifruit genus found throughout a wide range in eastern Asian from Indonesia to Siberia. In this study, phylogenetic relationships of Actinidia species with sympatric distributions were investigated using three chloroplast introns (trnL-F, atpB-rbcL and rpl32-trnL) and three Exon primed intron-crossing (EPIC) markers. Chloroplast phylogeny supports non-monophyly of the five species studied excluding Actinidia fulvicoma var. lanata. The non-monophyly was also revealed by EPIC markers. Our results showed EPIC markers are more variable and informative for phylogenetic inference than that of chloroplast markers. The incongruences between loci from the plastid and nuclear DNA phylogenic trees may stem from incomplete lineage sorting or historical introgression hybridization. Incomplete lineage sorting may explain the non-monophyly between Actinidia chrysantha (section Maculatae) and other four species (section Stellatae), and introgression hybridization and high level of interspecific gene flow may explain the non-monophyly among the species of sect. Stellatae. Thus, natural hybridization and introgression may be common in Actinidia with sympatric distribution.     

Discordance of species trees with their most likely gene trees

Because of the stochastic way in which lineages sort during speciation, gene trees may differ in topology from each other and from species trees. Surprisingly, assuming that genetic lineages follow a coalescent model of within-species evolution, we find that for any species tree topology with five or more species, there exist branch lengths for which gene tree discordance is so common that the most likely gene tree topology to evolve along the branches of a species tree differs from the species phylogeny. This counterintuitive result implies that in combining data on multiple loci, the straightforward procedure of using the most frequently observed gene tree topology as an estimate of the species tree topology can be asymptotically guaranteed to produce an incorrect estimate. We conclude with suggestions that can aid in overcoming this new obstacle to accurate genomic inference of species phylogenies.