Phylogenetic resolution and systematics of the Asian tree frogs,Rhacophorus (Rhacophoridae,Amphibia)

Li, J.T., Li, Y., Murphy, R.W., Rao, D.‐Q. & Zhang, Y.‐P. (2012). Phylogenetic resolution and systematics of the Asian tree frogs, Rhacophorus (Rhacophoridae, Amphibia). —Zoologica Scripta, 41, 557–570. The treefrog genus Rhacophorus, a large genus with 80 species, has a wide range, occurring eastward from India to China, Japan, South‐east Asia, the Greater Sunda Islands and the Philippines. The phylogenetic relationships and taxonomic recognition of many species are very controversial. To stabilize the taxonomy, the phylogenetic relationships among about 52 species are investigated from 96 samples using mtDNA sequence data. Matrilineal relationships based on maximum likelihood and Bayesian inference methods resolve three well‐supported lineages (A, B and C), although the phylogenetic relationships among three lineages remain ambiguous. Analyses support recognition of two previously assigned subgenera, Leptomantis and Rhacophorus, and these correspond to lineages A and B, respectively. Given that we have three strongly supported lineages, that these lineages are morphologically distinct, and the constrained geographic distributions of these groups, we recognize each lineage as a taxon. Subgenus Leptomantis includes species mainly from Malaysia, Indonesia and the Philippines. Subgenus Rhacophorus contains a mix of species occurring in India, Indochina and southern China. Lineage C accommodates species distributed mostly in East Asia, including Japan and China. Based on genetic and morphological data from type localities, the taxonomic recognition of some species needs to be reconsidered. Rhacophorus pingbianensis and Polypedates spinus are considered as junior synonyms of Rhacophorus duboisi. Specimens of Rhacophorus rhodopus from Vietnam and Hainan, China likely represent an undescribed, cryptic species.     

Molecular phylogeny and phylogeography of the Cuban cave-fishes of the genus Lucifuga: evidence for cryptic allopatric diversity

Underground environments are increasingly recognized as reservoirs of faunal diversity. Extreme environmental conditions and limited dispersal ability of underground organisms have been acknowledged as important factors promoting divergence between species and conspecific populations. However, in many instances, there is no correlation between genetic divergence and morphological differentiation. Lucifuga Poey is a stygobiotic fish genus that lives in Cuban and Bahamian caves. In Cuba, it offers a unique opportunity to study the influence of habitat fragmentation on the genetic divergence of stygobiotic species and populations. The genus includes four species and one morphological variant that have contrasting geographical distributions. In this study, we first performed a molecular phylogenetic analysis of the Lucifuga Cuban species using mitochondrial and nuclear markers. The mitochondrial phylogeny revealed three deeply divergent clades that were supported by nuclear and morphological characters. Within two of these main clades, we identified five lineages that are candidate cryptic species and a taxonomical synonymy between Lucifuga subterranea and Lucifuga teresinarum. Secondly, phylogeographic analysis using a fragment of the cytochrome b gene was performed for Lucifuga dentata, the most widely distributed species. We found strong geographical organization of the haplotype clades at different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. At a larger temporal scale, these processes could also explain the diversification and the distribution of the different species.     

Phylogeography of the antilopine wallaroo (Macropus antilopinus) across tropical northern Australia

The distribution of antilopine wallaroo, Macropus antilopinus, is marked by a break in the species’ range between Queensland and the Northern Territory, coinciding with the Carpentarian barrier. Previous work on M. antilopinus revealed limited genetic differentiation between the Northern Territory and Queensland M. antilopinus populations across this barrier. The study also identified a number of divergent lineages in the Northern Territory, but was unable to elucidate any geographic structure. Here, we re‐examine these results to (1) determine phylogeographic patterns across the range of M. antilopinus and (2) infer the biogeographic barriers associated with these patterns. The tropical savannahs of northern Australia: from the Cape York Peninsula in the east, to the Kimberley in the west. We examined phylogeographic patterns in M. antilopinus using a larger number of samples and three mtDNA genes: NADH dehydrogenase subunit 2, cytochrome b, and the control region. Two datasets were generated and analyzed: (1) a subset of samples with all three mtDNA regions concatenated together and (2) all samples for just control region sequences that included samples from the previous study. Analysis included generating phylogenetic trees based on Bayesian analysis and intraspecific median‐joining networks. The contemporary spatial structure of M. antilopinus mtDNA lineages revealed five shallow clades and a sixth, divergent lineage. The genetic differences that we found between Queensland and Northern Territory M. antilopinus samples confirmed the split in the geographic distribution of the species. We also found weak genetic differentiation between Northern Territory samples and those from the Kimberley region of Western Australia, possibly due to the Kimberley Plateau–Arnhem Land barrier. Within the Northern Territory, two clades appear to be parapatric in the west, while another two clades are broadly sympatric across the Northern Territory. MtDNA diversity of M. antilopinus revealed an unexpectedly complex evolutionary history involving multiple sympatric and parapatric mtDNA clades across northern Australia. These phylogeographic patterns highlight the importance of investigating genetic variation across distributions of species and integrating this information into biodiversity conservation.     

Multiple transisthmian divergences,extensive cryptic diversity,occasional long‐distance dispersal,and biogeographic patterns in a marine coastal isopod with an amphi‐American distribution

Excirolana braziliensis is a coastal intertidal isopod with a broad distribution spanning the Atlantic and Pacific tropical and temperate coasts of the American continent. Two separate regional studies (one in Panama and one in Chile) revealed the presence of highly genetically divergent lineages, implying that this taxon constitutes a cryptic species complex. The relationships among the lineages found in these two different regions and in the rest of the distribution, however, remain unknown. To better understand the phylogeographic patterns of E. braziliensis, we conducted phylogenetic analyses of specimens from much of its entire range. We obtained DNA sequences for fragments of four mitochondrial genes (16S rDNA, 12S rDNA, COI, and Cytb) and also used publicly available sequences. We conducted maximum likelihood and Bayesian phylogenetic reconstruction methods. Phylogeographic patterns revealed the following: (1) new highly divergent lineages of E. braziliensis; (2) three instances of Atlantic–Pacific divergences, some of which appear to predate the closure of the Isthmus of Panama; (3) the distributional limit of highly divergent lineages found in Brazil coincides with the boundary between two major marine coastal provinces; (4) evidence of recent long‐distance dispersal in the Caribbean; and (5) populations in the Gulf of California have closer affinities with lineages further south in the Pacific, which contrasts with the closer affinity with the Caribbean reported for other intertidal organisms. The high levels of cryptic diversity detected also bring about challenges for the conservation of this isopod and its fragile environment, the sandy shores. Our findings underscore the importance of comprehensive geographic sampling for phylogeographic and taxonomical studies of broadly distributed putative species harboring extensive cryptic diversity.     

Constrained body shape among highly genetically divergent allopatric lineages of the supralittoral isopod Ligia occidentalis (Oniscidea)

Multiple highly divergent lineages have been identified within Ligia occidentalis sensu lato, a rocky supralittoral isopod distributed along a ~3000 km latitudinal gradient that encompasses several proposed marine biogeographic provinces and ecoregions in the eastern Pacific. Highly divergent lineages have nonoverlapping geographic distributions, with distributional limits that generally correspond with sharp environmental changes. Crossbreeding experiments suggest postmating reproductive barriers exist among some of them, and surveys of mitochondrial and nuclear gene markers do not show evidence of hybridization. Populations are highly isolated, some of which appear to be very small; thus, the effects of drift are expected to reduce the efficiency of selection. Large genetic divergences among lineages, marked environmental differences in their ranges, reproductive isolation, and/or high isolation of populations may have resulted in morphological differences in L. occidentalis, not detected yet by traditional taxonomy. We used landmark‐based geometric morphometric analyses to test for differences in body shape among highly divergent lineages of L. occidentalis, and among populations within these lineages. We analyzed a total of 492 individuals from 53 coastal localities from the southern California Bight to Central Mexico, including the Gulf of California. We conducted discriminant function analyses (DFAs) on body shape morphometrics to assess morphological variation among genetically differentiated lineages and their populations. We also tested for associations between phylogeny and morphological variation, and whether genetic divergence is correlated to multivariate morphological divergence. We detected significant differences in body shape among highly divergent lineages, and among populations within these lineages. Nonetheless, neither lineages nor populations can be discriminated on the basis of body shape, because correct classification rates of cross‐validated DFAs were low. Genetic distance and phylogeny had weak to no effect on body shape variation. The supralittoral environment appears to exert strong stabilizing selection and/or strong functional constraints on body shape in L. occidentalis, thereby leading to morphological stasis in this isopod.     

Three reciprocally monophyletic mtDNA lineages elucidate the taxonomic status of Grant’s gazelles

The intraspecific phylogeography of Grant’s gazelles Nanger granti was assessed with mitochondrial DNA control region sequences. Samples of 177 individuals from 17 Kenyan and Tanzanian populations were analysed. Three highly divergent, reciprocally monophyletic lineages were found, with among group net nucleotide distances of 8–12%. The three lineages—notata, granti and petersii—grouped populations according to their geographic origin, encompassing populations in the north, southwest, and east, respectively. The mtDNA lineages reflected distinct evolutionary trajectories, and the data are discussed in reference to the four currently recognised subspecies. We suggest Grant’s gazelles be raised to the superspecies Nanger (granti) comprising three taxonomic units corresponding to the three mtDNA lineages. There was no evidence of gene flow between the notata and granti lineages, despite their geographic proximity, suggesting reproductive isolation. These constitute evolutionary significant units within the adaptive evolutionary framework. Due to its restricted geographic distribution and genetic and morphological distinctiveness, we suggest the petersii lineage be raised to the species Nanger (granti) petersii within the Grant’s gazelles superspecies.     

Exaggerated male forelegs are not more differentiated than wing morphology in two widespread sister species of black scavenger flies

Sexual selection represents a potent force that can drive rapid population differentiation in traits related to reproductive success. Hence, sexual traits are expected to show greater population divergence than non-sexual traits. We test this prediction by exploring patterns of morphological differentiation of the exaggerated fore femur (a male-specific sexual trait) and the wing (a non-sexual trait) among allopatric and sympatric populations of the widespread sister dung fly species Sepsis neocynipsea and Sepsis cynipsea (Diptera: Sepsidae). While both species occur in Eurasia, S. neocynipsea also abounds in North America, albeit previous studies suggest strong differentiation in morphology, behavior, and mating systems. To evaluate the degree of differentiation expected under neutrality between S. cynipsea, European S. neocynipsea, and North American S. neocynipsea, we genotyped 30 populations at nine microsatellite markers, revealing almost equal differentiation between and minor differentiation among geographic populations within the three lineages. Landmark-based analysis of 18 populations reared at constant 18 and 24°C in a laboratory common garden revealed moderate temperature-dependent phenotypic plasticity and significant heritable differentiation in size and shape of male forelegs and wings among iso-female lines of the three lineages. Following the biological species concept, there was weaker differentiation between cross-continental populations of S. neocynipsea relative to S. cynipsea, and more fore femur differentiation between the two species in sympatry versus allopatry (presumably due to character displacement). Contrary to expectation, wing morphology showed as much shape differentiation between evolutionary independent lineages as fore femora, providing no evidence for faster diversification of traits primarily engaged in mating.     

Evolution of hybridogenetic lineages in Cataglyphis ants

In most social Hymenoptera, a diploid egg develops into either a queen or a worker depending on environmental conditions. Hybridogenetic Cataglyphis ants display a bizarre genetic system, where queen‐worker caste determination is primarily determined by genetic factors. In hybridogenetic populations, all workers are F1 hybrids of two distinct lineages, whereas new queens are nearly always pure‐lineage individuals produced by clonal reproduction. The distribution and evolutionary history of these hybridogenetic populations have not yet been thoroughly analysed. Here, we studied the phylogeographic distribution of hybridogenetic populations in two closely related Spanish species: Cataglyphis humeya and Cataglyphis velox. Hybridogenesis has been previously documented in a locality of C. velox, but whether this system occurs elsewhere within the range of the two species was yet unknown. Queens and workers from 66 localities sampled across the range of the species were genotyped at 18 microsatellite markers to determine whether queens were produced by parthenogenesis and whether workers were hybrids of divergent lineages. Populations with F1 hybrid workers were identified by combining genetic, geographical and mating assortments data. In most populations of C. velox, workers were found to be hybrids of two divergent lineages. Workers were however produced via random mating in two marginal populations of C. velox, and in all populations studied of its sister species C. humeya. High‐throughput sequencing data were obtained to confirm inferences based on microsatellites and to characterize relationships between populations. Our results revealed a complicated history of reticulate evolution that may account for the origin of hybridogenetic lineages in Cataglyphis.     

Scrophularia arguta,a widespread annual plant in the Canary Islands: a single recent colonization event or a more complex phylogeographic pattern?

Many studies have addressed evolution and phylogeography of plant taxa in oceanic islands, but have primarily focused on endemics because of the assumption that in widespread taxa the absence of morphological differentiation between island and mainland populations is due to recent colonization. In this paper, we studied the phylogeography of Scrophularia arguta, a widespread annual species, in an attempt to determine the number and spatiotemporal origins of dispersal events to Canary Islands. Four different regions, ITS and ETS from nDNA and psbA‐trnH and psbJ‐petA from cpDNA, were used to date divergence events within S. arguta lineages and determine the phylogenetic relationships among populations. A haplotype network was obtained to elucidate the phylogenetic relationships among haplotypes. Our results support an ancient origin of S. arguta (Miocene) with expansion and genetic differentiation in the Pliocene coinciding with the aridification of northern Africa and the formation of the Mediterranean climate. Indeed, results indicate for Canary Islands three different events of colonization, including two ancient events that probably happened in the Pliocene and have originated the genetically most divergent populations into this species and, interestingly, a recent third event of colonization of Gran Canaria from mainland instead from the closest islands (Tenerife or Fuerteventura). In spite of the great genetic divergence among populations, it has not implied any morphological variation. Our work highlights the importance of nonendemic species to the genetic richness and conservation of island flora and the significance of the island populations of widespread taxa in the global biodiversity.     

BIODIVERSITY RESEARCH: Genetic diversity in two introduced biofouling amphipods (Ampithoe valida & Jassa marmorata) along the Pacific North American coast: investigation into molecular identification and cryptic diversity

Aim We investigated patterns of genetic diversity among invasive populations of Ampithoe valida and Jassa marmorata from the Pacific North American coast to assess the accuracy of morphological identification and determine whether or not cryptic diversity and multiple introductions contribute to the contemporary distribution of these species in the region. Location Native range: Atlantic North American coast; Invaded range: Pacific North American coast. Methods We assessed indices of genetic diversity based on DNA sequence data from the mitochondrial cytochrome c oxidase subunit I (COI) gene, determined the distribution of COI haplotypes among populations in both the invasive and putative native ranges of A. valida and J. marmorata and reconstructed phylogenetic relationships among COI haplotypes using both maximum parsimony and Bayesian approaches. Results Phylogenetic inference indicates that inaccurate species‐level identifications by morphological criteria are common among Jassa specimens. In addition, our data reveal the presence of three well supported but previously unrecognized clades of A. valida among specimens in the north‐eastern Pacific. Different species of Jassa and different genetic lineages of Ampithoe exhibit striking disparity in geographic distribution across the region as well as substantial differences in genetic diversity indices. Main conclusions Molecular genetic methods greatly improve the accuracy and resolution of identifications for invasive benthic marine amphipods at the species level and below. Our data suggest that multiple cryptic introductions of Ampithoe have occurred in the north‐eastern Pacific and highlight uncertainty regarding the origin and invasion histories of both Jassa and Ampithoe species. Additional morphological and genetic analyses are necessary to clarify the taxonomy and native biogeography of both amphipod genera.     

Identifying evolutionarily significant units for conservation of the endangered Malus sieversii using genome‐wide RADseq data

Malus sieversii, a wild progenitor of the domesticated apple, is an endangered species and is assigned second conservation priority by the China Plant Red Data Book. It is urgent to carry out in situ conservation of this species, but previous studies have not identified evolutionarily significant units (ESUs) for conservation management. In this study, we investigated the genetic diversity and relationships of six M. sieversii populations from China using integrated analysis of microsatellite (nSSR) data, genome‐wide SNPs and previous results in order to propose a reasonable conservation management. The results showed that levels of genetic diversity were inconsistently reflected by our nSSR and previous studies, suggesting that indices of genetic diversity are not effective to identify priority conservation areas for M. sieversii. Based on the selection criteria of ESUs for endangered species conservation, ESUs should reflect lineage divergence, geographical separation and different adaptive variation. Our phylogenetic tree based on genome‐wide SNPs yielded a clear relationship of divergent lineages among M. sieversii populations, leading to new different from those of previous studies. Three independent lineages, including the pairs of populations Huocheng‐Yining, Gongliu‐Xinyuan and Tuoli‐Emin, were identified. The geographic distances between populations among the different phylogenetic lineages were much greater than those within the same phylogenetic lineage. A cluster analysis on environmental variables showed that the three independent lineages inhabit different environmental conditions, suggesting that they may have adapted to different environments. Based on the results, we propose that three independent ESUs should be recognized as conservation units for M. sieversii in China.     

Multiple cryptic species of sympatric generalists within the avian blood parasite Haemoproteus majoris

The avian haemosporidian parasite Haemoproteus majoris has been reported to infect a wide range of passerine birds throughout the Holarctic ecozone. Five cytochrome b (cyt b) lineages have been described as belonging to the morphological species H. majoris, and these form a tight phylogenetic cluster together with 13 undescribed lineages that differ from each other by < 1.2% in sequence divergence. Records in a database (MalAvi) that contains global findings of haemosporidian lineages generated by universal primers suggest that these lineages vary substantially in host distribution. We confirm this pattern in a data set collected at Lake Kvismaren, Sweden, where three of the generalist lineages have local transmission. However, whether these lineages represent intraspecific mitochondrial diversity or clusters of cryptic species has previously not been examined. In this study, we developed novel molecular markers to amplify the partial segments of four nuclear genes to determine the level of genetic diversity and gene phylogenies among the five morphologically described cyt b lineages of H. majoris. All five cyt b lineages were strongly associated with unique nuclear alleles at all four nuclear loci, indicating that each mitochondrial lineage represents a distinct biological species. Within lineages, there was no apparent association between nuclear alleles and host species, indicating that they form genetically unstructured populations across multiple host species.     

Locomotor responses to salt stress in native and invasive mud‐tidal gastropod populations (Batillaria)

Plasticity in salt tolerance can be crucial for successful biological invasions of novel habitats by marine gastropods. The intertidal snail Batillaria attramentaria, which is native to East Asia but invaded the western shores of North America from Japan 80 years ago, provides an opportunity to examine how environmental salinity may shape behavioral and morphological traits. In this study, we compared the movement distance of four B. attramentaria populations from native (Korea and Japan) and introduced (United States) habitats under various salinity levels (13, 23, 33, and 43 PSU) during 30 days of exposure in the lab. We sequenced a partial mitochondrial CO1 gene to infer phylogenetic relationships among populations and confirmed two divergent mitochondrial lineages constituting our sample sets. Using a statistical model‐selection approach, we investigated the effects of geographic distribution and genetic composition on locomotor performance in response to salt stress. Snails exposed to acute low salinity (13 PSU) reduced their locomotion and were unable to perform at their normal level (the moving pace of snails exposed to 33 PSU). We did not detect any meaningful differences in locomotor response to salt stress between the two genetic lineages or between the native snails (Japan vs. Korea populations), but we found significant locomotor differences between the native and introduced groups (Japan or Korea vs. the United States). We suggest that the greater magnitude of tidal salinity fluctuation at the US location may have influenced locomotor responses to salt stress in introduced snails.     

Multilocus and morphological analysis of south-eastern Iberian Wall lizards (Squamata,Podarcis)

The phylogenetic relationships among the wall lizards of the Podarcis hispanicus complex that inhabit the south-east (SE) of the Iberian Peninsula and other lineages of the complex remain unclear. In this study, four mitochondrial and two nuclear markers were used to study genetic relationships within this complex. The phylogenetic analyses based on mtDNA gene trees constructed with ML and BI, and a species tree using *BEAST support three divergent clades in this region: the Valencia, Galera and Albacete/Murcia lineages. These three lineages were also corroborated in species delimitation analyses based on mtDNA using bPTP, mPTP, GMYC, ABGD and BAPS. Bayesian inference species delimitation method (BPP) based on both nuclear data and a combined data set (mtDNA + nuclear) showed high posterior probabilities for these three SE lineages (≥0.94) and another Bayesian analysis (STACEY) based on combined data set recovered the same three groups in this region. Divergence time dating of the species tree provided an estimated divergence of the Galera lineage from the other SE group (Podarcis vaucheri, (Albacete/Murcia, Valencia)) at 12.48 Ma. During this period, the Betic–Rifian arc was isolated, which could have caused the isolation of the Galera form distributed to the south of the Betic Corridor. Although lizards from the Albacete/Murcia and Galera lineage are morphologically similar, they clearly represent distinct genetic lineages. The noteworthy separation of the Galera lineage enables us to conclude that this lineage must be considered as a new full species.     

Mitochondrial DNA phylogeography of the Cotesia flavipes complex of parasitic wasps (Hymenoptera: Braconidae)

The Cotesia flavipes species complex of parasitic wasps are economically important worldwide for the biological control of lepidopteran stem borers. The complex currently comprises three species: Cotesia flavipes Cameron, C. sesamiae (Cameron) and C. chilonis (Matsumura) (Hymenoptera: Braconidae), which appear morphologically similar. Despite their economic importance, little is known about the genetic diversity and phylogeography of these parasitoids. Differences in the biology of geographic populations have generally been interpreted as genetic divergence among strains, but direct genetic evidence is lacking. In Australia, several stem borer pests in neighbouring countries have been identified as significant threats to the sugar industry. However, the status of C. flavipes in Australia is unknown. To examine the genetic variation among worldwide populations of the C. flavipes complex and investigate the status of the Australian C. flavipes-like species, partial sequence data were generated for mitochondrial gene regions, 16S rRNA and COI. Parsimony, minimum evolution and Bayesian analyses based on 21 geographic populations of the complex and four outgroups supported the monophyly of the complex and the existence of genetically divergent populations of C. flavipes and C. sesamiae. The geographically isolated Australian haplotypes formed a distinct lineage within the complex and were ~3.0% divergent from the other species. The results indicated that historical biogeographic barriers and recent biological control introductions play an important role in structuring lineages within these species. This study provides a phylogeographical context for examining adaptive evolution and host range within biologically divergent strains of the C. flavipes complex.     

Patterns of morphological and ecological similarities of small-eared shrews (Soricidae,Cryptotis) in tropical montane cloud forests from Mesoamerica

It has been proposed that high morphological similarity between closely related species of small-eared shrews resulted from a recent divergence and intermittent population connectivity, presumably due to Pleistocene climatic fluctuations and associated changes in forest habitat distribution. Here we examined the morphological variation of two sister species of small-eared shrews inhabiting cloud forests from Mexico, Cryptotis obscurus and C. mexicanus. We then used ecological niche modelling to provide compelling evidence for current environmental barriers for population connectivity, and for detecting divergent ecological niches between candidate species. Our results indicated that the species boundaries in this clade should be subject to change. High morphological similarity suggested that populations of C. obscurus and C. mexicanus located west of the Isthmus of Tehuantepec, a major geographic barrier for montane species, are conspecific. Niche divergence between these two putative species was not supported indicating niche conservatism across the evolutionary history of these small-eared shrews. In addition, several barriers seem to play a main role for current lineage divergence between populations within this clade. The population located east of the Isthmus, previously referred to C. mexicanus, might prove to represent a new species based on morphological distinction and current geographic isolation. We have highlighted that estimating species’ potential distributions provides insights to evaluating the effect of geographic barriers on lineage divergence and making stronger inferences when delimiting species.     

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.     

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.