Refugia within refugia – patterns in endemism and genetic divergence are linked to Late Quaternary climate stability in the Iberian Peninsula

In Europe, southern peninsulas served as major refugia during Pleistocene cold periods. However, growing evidence has revealed complex patterns of glacial survival within these southern regions, with multiple glacial refugia within each larger refugial area. We investigated the extent to which patterns of endemism and phylogeographic are concordant across animal species in the Iberian Peninsula, one of the most important unglaciated areas in Europe during the Pleistocene, can be explained in terms of climatic stability. We found that historical climatic stability (notably climate velocity measures integrating macroclimatic shifts with local spatial topoclimate gradients) was often among the most important predictors of endemic species richness for different taxonomic groups using models that also incorporated measures of modern climate. Furthermore, for some taxonomic groups, climatic stability was also correlated with patterns of spatial concordance in interpopulation genetic divergence across multiple taxa, and private haplotypes were more frequently found in relatively stable areas. Overall, our results suggest that both endemism patterns and cross‐taxa concordant phylogeographic patterns across the Iberian Peninsula to some extent are linked to spatial variation in Late Quaternary climate stability, in agreement with the proposed ‘refugia‐within‐refugia’ scenario. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 13–28.     

Phylogeny and evolution of the genus Barbus in the Iberian Peninsula as revealed by allozyme electrophoresis

Variation of 14 enzyme systems encoded by 31 presumptive loci in different barbel species was studied using starch gel electrophoresis. Eighteen samples representing 11 Barbus tetraploid taxa were taken, including 10 populations from the Iberian Peninsula, six from other parts of western and southern Europe, one from northern Africa and one diploid species as outgroup from Asia Minor. The genetic analysis reassessed of the taxonomic status of the Iberian barbels into two distinct species groups. The first group included B. bocagei, B. comiza, B. graellsii, B. gulraonis, B. microcephalus , and B. sclateri , that aligned with B. callensis from northern Africa and with B. apoensis from Asia Minor; the other group included B. haasi and B. merldlonahs that was related to the European species, B. barbus, B. plebejus and B. peloponnesius . These groups are probably not monophyletic. It is suggested that the isolation of the Iberian Peninsula from Europe since the Oligocene-Miocene may explain the genetic affinities of the Iberian barbels with those of North African rather than with the European group.     

Where are you from,stranger? The enigmatic biogeography of North African pond turtles (Emys orbicularis)

The European pond turtle (Emys orbicularis) is a Nearctic element in the African fauna and thought to have invaded North Africa from the Iberian Peninsula. All North African populations are currently identified with the subspecies E. o. occidentalis. However, a nearly range-wide sampling in North Africa used for analyses of mitochondrial and microsatellite DNA provides evidence that only Moroccan populations belong to this taxon, while eastern Algerian and Tunisian pond turtles represent an undescribed distinct subspecies. These two taxa are most closely related to E. o. galloitalica with a native distribution along the Mediterranean coast of northern Spain through southern France to western and southern Italy. This group is sister to a clade comprising several mitochondrial lineages and subspecies of E. orbicularis from Central and Eastern Europe plus Asia, and the successive sisters are E. o. hellenica and E. trinacris. Our results suggest that E. orbicularis has been present in North Africa longer than on the Iberian Peninsula and that after an initial invasion of North Africa by pond turtles from an unknown European source region, there was a phase of diversification in North Africa, followed by a later re-invasion of Europe by one of the African lineages. The differentiation of pond turtles in North Africa parallels a general phylogeographic paradigm in amphibians and reptiles, with deeply divergent lineages in the western and eastern Maghreb. Acknowledging their genetic similarity, we propose to synonymize the previously recognized Iberian subspecies E. o. fritzjuergenobsti with E. o. occidentalis sensu stricto. The seriously imperiled Moroccan populations of E. o. occidentalis represent two Management Units different in mitochondrial haplotypes and microsatellite markers. The conservation status of eastern Algerian pond turtles is unclear, while Tunisian populations are endangered. Considering that Algerian and Tunisian pond turtles represent an endemic taxon, their situation throughout the historical range should be surveyed to establish a basis for conservation measures.     

Ecological differentiation and cladogenesis of Baldellia (L.) Parl. (Alismataceae)

The Iberian Peninsula is home to a large number of endemic species that often show strong genetic subdivisions which indicate possible isolation of populations in the past. Numerous phylogeographic and phylogenetic studies on Iberian flora and fauna have revealed this territory to be the origin of genetic differentiation during the Pleistocene Ice Age. To better understand the influence of Plio-Pleistocene climatic oscillations on the cladogenesis and population structuring of refugial taxa restricted to freshwaters of the Iberian Peninsula, we studied the evolutionary history of the small European genus Baldellia, both from a phylogenetic/phylogeographic perspective and regarding its ecological characteristics. Molecular analyses were based on four plastid regions and nuclear ITS2 ribosomal spacer from individuals belonging to geographically isolated Baldellia populations, while the relationships between climate and Baldellia species distribution was investigated by redundancy analysis applied to the climatic data with the occurrence of Baldellia species as external variables. The absence of a clear phylogenetic signal most likely reflects a common and recent cladogenesis for Baldellia. In contrast, the more rapidly evolving plastid sequences indicate a recent phylogeographic history of refugial repartition and ecological segregation probably occurring during the last glaciation, a scenario well supported by marked differences in species?? ecological preferences and climate analysis.     

Phylogenetic and biosystematic relationships in four highly disjunct polyploid complexes in the subgenera Ceterach and Phyllitis in Asplenium (Aspleniaceae)

Phylogenetic studies using DNA sequences of two chloroplast regions, rbcL and trnL-F, demonstrate that the proposed genus Ceterach is a small clade within the large genus Asplenium, and sister to the Phyllitis clade. The Ceterach clade is characterised by irregular anastomosing veins and often densely scaled leaf blades. Its taxonomic status as a group nested within Asplenium is confirmed, and it is accepted here as a subgenus with seven species. The Ceterach clade comprises four lineages that correspond to disjunct polyploid complexes: the A. aureum clade forming a polyploid complex (4×, 6×, 8×) in Macaronesia, the A. ceterach clade forming a polyploid complex (2×, 4×, 6×) in the Mediterranean Basin, the A. paucivenosum clade (4×, 6×) in central Asia, and the A. dalhousiae clade (2×) with a disjunct distribution in the Himalaya, Yemen and Eritrea, and southwestern North America. Asplenium paucivenosum is sister to all other members of the Ceterach clade, whereas A. dalhousiae is sister to the A. aureum clade that includes tetraploid A. aureum, hexaploid A. lolegnamense, and octoploid A. parvifolium. Asplenium ceterach and its variations – including the hexaploid A. ceterach subsp. mediterraneum subsp. nov. first described below – form a monophyletic unit, sister to a clade consisting of A. aureum and A. dalhousiae. Asplenium cordatum from Africa and A. haugthonii from the isolated atlantic island of St. Helena are not members of the Ceterach clade, which suggests that leaf blades with dense indumenta have evolved at least twice within asplenioid ferns. The allotetraploid species A. hybridum has the chloroplast DNA from A. ceterach, and therefore the latter species is the maternal ancestor of the former. The other parent of this hybrid species is A. sagittatum that is nested within the sister clade of Ceterach, the Phyllitis clade comprising A. sagittatum and A. scolopendrium. The findings suggest that the current distribution of Ceterach is either the result of long-distance dispersal or represents fragmented relicts of a previously more widely distributed species.     

Karyosystematic study of Potentilla L. subgen. Potentilla (Rosaceae) in the Iberian Peninsula

Chromosome numbers in 80 populations belonging to 18 species of Potentilla L. subgen. Potentilla from the Iberian Peninsula and two of P. maura, a North African endemic taxon, have been counted. The basic number of chromosomes is always x = 7 and these chromosomes are small (between 1 and 2 μm). For three species, the number of chromosomes is reported for the first time and, for another six, this number has been established in Iberian representatives. Moreover, new ploidy levels have been obtained for P. hispanica and P. crantzii with regard to their entire distribution area, and in P. cinerea and P. neumanniana for the Iberian Peninsula. Some taxonomic, phylogenetic and phytogeographic comments are made for several species or groups of species from the West Mediterranean region. In 13 species only one ploidy level has been found, but six species have several ploidy levels. Seven ploidy levels occur in the investigated taxa. The frequency of each ploidy level represented within Iberian Potentilla is analysed and the data are compared with those available for taxa from the rest of the distribution area of the genus.     

Phylogenetic relationships among Spanish Barbus species (Pisces,Cyprinidae) shown by RAPD markers

Random amplified polymorphic DNA (RAPD) markers were used to estimate the population structure and phylogenetic relationships among the eight species of the genus Barbus that inhabit the Iberian Peninsula. Ten random oligodecamers were used to amplify DNA from 232 fish sampled from 15 populations. A total of 270 markers were detected that revealed low levels of genetic variability. The conclusions of cluster analysis indicate two main branches and three well-differentiated groups: north-eastern, Mediterranean and Atlantic. This clustering mainly reflects the evolutionary history of the genus, which is closely related to the paleogeography of the Iberian Peninsula. The contribution and application of these results to the conservation of the species, to their taxonomic status and to the process of colonization of the Iberian Peninsula by the genus Barbus are discussed.     

利用整合分类学方法进行蕨类植物复合体的物种划分: 以线裂铁角蕨复合体为例

广泛的杂交和多倍化使得铁角蕨属(Asplenium)下存在着许多分类困难的物种复合体, 针对这些类群进行整合分类学的研究, 有助于我们更加全面和深入地理解物种的界限以及形成机制。线裂铁角蕨复合体(Asplenium coenobiale complex)是铁角蕨属下一个形态多样性较高的类群, 由于缺乏全面取样和系统研究, 该复合体的物种划分长期存在争议。本研究选取线裂铁角蕨复合体中形态变异和地理分布具有代表性的个体, 通过孢粉学研究确定该类群的生殖特性, 运用流式细胞分析获取倍性信息, 同时结合叶绿体和核基因组片段系统发生分析的证据, 对该类群的系统演化关系和起源方式进行了探讨。结果表明: (1)虽然部分孢子囊败育的情况在线裂铁角蕨复合体中十分普遍, 但正常孢子囊内形成的64个孢子说明该类群植物仍能进行正常的有性生殖; (2)该复合体中存在着倍性变异, 其中多角铁角蕨(A. cornutissimum)是二倍体, 而其他成员均为四倍体; (3)依据母系遗传的叶绿体序列所构建的系统发生关系将该类群划为4个分支, 与基于核基因序列构建的系统树存在冲突, 这暗示杂交可能在该复合体的形成过程中起到了重要的推动作用。综上所述, 我们建议将线裂铁角蕨复合体划分为4个物种, 即同源四倍体新种马关铁角蕨(A. maguanense sp. nov.), 二倍体多角铁角蕨, 以及两个由同一对亲本正反交产生的异源四倍体线裂铁角蕨(A. coenobiale)和叶基宽铁角蕨(A. pulcherrimum)。     

Historical biogeography of the lineage Atopochara trivolvis PECK 1941 (Cretaceous Charophyta)

The evolutionary species Atopochara trivolvis represents a unique case of a charophyte with a well-documented, long-lasting fossil record of biostratigraphical interest and a worldwide distribution. This provides an excellent basis for analysing the historical biogeography of this important species. The origin of the lineage is recorded in the Early Berriasian of Western Europe, where an archipelago-like palaeogeography would have favoured allopatric speciation and small-sized original biogeographic ranges. The species first reached Eastern Europe during the Valanginian and expanded to China during the Hauterivian. It expanded worldwide during the Barremian–Aptian, first colonizing North Africa and South America through the Guinea Corridor and subsequently reaching North America. The Barremian–Aptian worldwide expansion was probably facilitated by the monoecious nature of A. trivolvis, which allowed for long-distance animal dispersal, as is the case with extant Chara. An initial fragmentation in the biogeographic range of A. trivolvis can be detected in the Albian–Cenomanian, when the species displayed up to 5 different morphotypes in scattered localities of Eurasia and North America. During the Turonian sea-level highstand A. trivolvis became extinct also from most of these localities, while, in the Latest Cretaceous plesiomorphic populations of A. trivolvis ulanensis survived in Mongolia and Northern China. The species became completely extinct near the K/T boundary.     

A systematic study of North American Angelica species (Apiaceae) based on nrDNA ITS and cpDNA sequences and fruit morphology

Angelica is a taxonomically complex genus widespread throughout the North Temperate Zone. Previous phylogenetic studies of the genus have focused primarily on its East Asian species. The relationships among its North American members, the monophyly of these species, and the value of fruit morphology in circumscribing its taxa have yet to be examined. This study represents the most comprehensive sampling of Angelica to date (100 species) and includes all 26 species in North America. Relationships are inferred using Bayesian inference, maximum likelihood, and maximum parsimony analyses of ITS sequences and, for multiple accessions of each North American species, cpDNA ndhF-rpl32, rpl32-trnL, and psbM-psbD sequences. The fruit morphological characters examined were those considered phylogenetically important in East Asian Angelica. The results revealed that the North American species fell into three major clades: North American Angelica clade, Archangelica clade, and the Eurasian Angelica clade. Angelica dawsonii has affinities with Lomatium brandegeei. Fourteen species within the North American Angelica clade were strongly supported as monophyletic. Two paraphyletic species resulted in new combinations in A. lineariloba and A. venenosa. Conflict between the ITS-derived and cpDNA-derived phylogenies and the lack of resolution in portions of the trees may be due to chloroplast capture and rapid species radiation. Fruit morphology supported some interspecific relationships based on molecular data, and relationships revealed by ITS and cpDNA data were roughly in accordance with fruit classification type and geographic distribution region, respectively. A diagnostic key based on fruit morphology is provided for the identification of the North American Angelica taxa.     

The medical and veterinary role of Ornithodoros erraticus complex ticks (Acari: Ixodida) on the Iberian Peninsula

Argasid ticks of the Ornithodoros erraticus complex are associated with traditional pig‐farming practices on the Iberian Peninsula and are also found elsewhere in North Africa, West Africa, and western Asia. The ticks associated with pig farming on the Iberian Peninsula are the only biological vectors of African swine fever virus (ASFV) known to occur in Europe, and their ecology makes them an extremely effective reservoir of both ASFV and the Borrelia species which cause tick‐borne relapsing fever (TBRF) in humans. The recent reappearance of ASFV in the European Union, coupled with evidence that Portuguese tick populations continue to harbor Borrelia despite a lack of confirmed human infections, suggest that these populations merit closer attention. In Portugal, a series of surveys over the last twenty‐five years indicates that the number of farm sites with tick infestations has declined and suggest that populations are sensitive to changes in farm management, particularly the use of modern pig housing. Various technologies have been suggested for the control of farm‐associated Ornithodoros ticks and related species but, in our opinion, farm management changes are still the most effective strategy for population control. Furthermore, we suggest that this species could probably be eradicated from Iberian pig farms.     

Autosomal and mtDNA Markers Affirm the Distinctiveness of Lions in West and Central Africa

The evolutionary history of a species is key for understanding the taxonomy and for the design of effective management strategies for species conservation. The knowledge about the phylogenetic position of the lion (Panthera leo) in West/Central Africa is largely based on mitochondrial markers. Previous studies using mtDNA only have shown this region to hold a distinct evolutionary lineage. In addition, anthropogenic factors have led to a strong decline in West/Central African lion numbers, thus, the conservation value of these populations is particularly high. Here, we investigate whether autosomal markers are concordant with previously described phylogeographic patterns, and confirm the unique position of the West/Central African lion. Analysis of 20 microsatellites and 1,454 bp of the mitochondrial DNA in 16 lion populations representing the entire geographic range of the species found congruence in both types of markers, identifying four clusters: 1) West/Central Africa, 2) East Africa, 3) Southern Africa and 4) India. This is not in line with the current taxonomy, as defined by the IUCN, which only recognizes an African and an Asiatic subspecies. There are no indications that genetic diversity in West/Central Africa lions is lower than in either East or Southern Africa, however, given this genetic distinction and the recent declines of lion numbers in this region, we strongly recommend prioritization of conservation projects in West/Central Africa. As the current taxonomic nomenclature does not reflect the evolutionary history of the lion, we suggest that a taxonomic revision of the lion is warranted.     

Determinant variables of Iberian Peninsula Aphodiinae diversity (Coleoptera, Scarabaeoidea, Aphodiidae)

Aim The aims of this paper are to examine diversity–variability patterns for species of Aphodiinae (Coleoptera, Scarabaeoidea, Aphodiidae) on the Iberian Peninsula, and to determine the factors that influence their geographic distribution. Location Iberian Peninsula (Spain and Portugal). Methods Data from 30 studies and their bibliographies on species of Iberian Peninsula Aphodiinae were compiled. The reliability of the inventories was evaluated using parametric species richness estimators. In addition, a further 11 variables related to rarity, geographic distribution, or phylogenetic diversity were considered. Diversity variables were analysed using principal components analysis to reduce the number of dependent variables. Subsequently, the effect of differences in locality size among the 30 studies was eliminated by calculating and retaining the residuals of the curvilinear relationship of each diversity variable with the area. Generalized linear models were used to examine the relationships between diversity and 17 environmental variables. The diversity variables and their residuals were also subject to trend surface analysis in order to identify the relevance of spatially structured variables that had not been considered. The contribution of explanatory variables was determined through hierarchical variance analysis. Results Principal components analysis of biodiversity variables revealed that most of the variability could be explained using three biodiversity indexes: BI1, correlated positively with species richness, widely distributed species, frequent species, abundant species, species occurring in North Africa, Europe and the Iberian Peninsula, and phylogenetic diversity; BI2, correlated positively with numbers of infrequent and African–Iberian species; and, BI3, correlated positively with numbers of endemic, non‐abundant, European, and Iberian‐restricted species. A latitudinal disjunction emerged in BI1, with maximum scores at the north‐western and southern corners, while maximum BI2 scores were found throughout the south, and maximum BI3 scores in the north‐west. For BI1, it was climate that had the greatest influence, followed by lithology, and livestock presence. Geographic variables were the most significant for BI2, followed by climate and livestock presence. Finally, for BI3, climate variables were the most important, while geography, lithology and livestock presence had some relevance. Main conclusions The relevance of geographic variables indicates that other unaccounted‐for factors that are spatially structured could possibly explain additional variation in Aphodiinae diversity. These factors may be historic in nature, relating to the species groups, namely the Ibero‐European and the Mediterranean or Afro‐Iberian. The northern pattern could reflect the fact that the Iberian Peninsula acted as a colonization route and as a refuge during the glacial/interglacial cycles, while the southern pattern could be a consequence of the connection between the Iberian Peninsula and North Africa during the Messinian crisis, and/or a historic relationship in common, related to human activity.     

Evidence for phylogeographically informative sequence variation in the mitochondrial control region of Atlantic brown trout

Complete sequencing of the mitochondrial control region was undertaken among brown trout Salmo trutta from North Atlantic areas where previous studies, based on smaller mtDNA fragments, failed to detect any phylogeographic signal. Comparison of sequences suggests that brown trout in the Iberian Peninsula and Scandinavia belong to largely divergent evolutionary units.     

Endemic species may have complex histories: within‐refugium phylogeography of an endangered Iberian vole

Glacial refugia protected and promoted biodiversity during the Pleistocene, not only at a broader scale, but also for many endemics that contracted and expanded their ranges within refugial areas. Understanding the evolutionary history of refugial endemics is especially important in the case of endangered species to recognize the origins of their genetic structure and thus produce better informed conservation practices. The Iberian Peninsula is an important European glacial refugium, rich in endemics of conservation concern, including small mammals, such as the Cabrera vole (Microtus cabrerae). This near‐threatened rodent is characterized by an unusual suite of genetic, life history and ecological traits, being restricted to isolated geographic nuclei in fast‐disappearing Mediterranean subhumid herbaceous habitats. To reconstruct the evolutionary history of the Cabrera vole, we studied sequence variation at mitochondrial, autosomal and sex‐linked loci, using invasive and noninvasive samples. Despite low overall mitochondrial and nuclear nucleotide diversities, we observed two main well‐supported mitochondrial lineages, west and east. Phylogeographic modelling in the context of the Cabrera vole's detailed fossil record supports a demographic scenario of isolation of two populations during the Last Glacial Maximum from a single focus in the southern part of the Iberian Peninsula. In addition, our data suggest subsequent divergence within the east, and secondary contact and introgression of the expanding western population, during the late Holocene. This work emphasizes that refugial endemics may have a phylogeographic history as rich as that of more widespread species, and conservation of such endemics includes the preservation of that genetic legacy.     

Unmasking cryptic species: morphometric and phylogenetic analyses of the Ibero‐North African Linaria incarnata complex

Linaria incarnata has been treated as a widely distributed Ibero‐North African species in the latest taxonomic reviews. Morphological and phylogenetic disparity between populations of this species has been previously reported. Here we present new morphological and phylogenetic evidence for the disintegration of L. incarnata into three distinct species: L. incarnata from the western Iberian Peninsula; L . mamorensis sp. nov. from north‐western Morocco; and L. onubensis from south‐western Spain. The relatively poor morphological differentiation between these taxa (which can be regarded as cryptic species) and their distinct phylogenetic positions indicate that characters of the L. incarnata morphotype have been acquired multiple times in the evolution of Linaria section Versicolores. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 00, 000–000.     

Global mitochondrial DNA phylogeography and biogeographic history of the antitropically and longitudinally disjunct marine bryozoan Membranipora membranacea L. (Cheilostomata): another cryptic marine sibling species complex?

The origin of disjunct distributions in high dispersal marine taxa remains an important evolutionary question as it relates to the formation of new species in an environment where barriers to gene flow are not always obvious. To reconstruct the relationships and phylogeographic history of the antitropically and longitudinally disjunct bryozoan Membranipora membranacea populations were surveyed with mtDNA cytochrome oxidase 1 (COI) sequences across its cosmopolitan range. Maximum parsimony, maximum likelihood and Bayesian genealogies revealed three deep clades in the North Pacific and one monophyletic clade each in the southeast Pacific (Chile), southwest Pacific (Australia/New Zealand), North Atlantic and southeast Atlantic (South Africa). Human-mediated dispersal has not impacted M. membranacea’s large-scale genetic structure. M. membranacea did not participate in the trans-arctic interchange. Episodic long-distance dispersal, combined with climatic vicariance can explain the disjunct distribution. Dispersal led southward across the tropics perhaps 13 mya in the East Pacific and again northwards perhaps 6 mya in the Eastern Atlantic to colonize the North Atlantic from the south, and along the West Wind Drift to colonize Australia. The clades differentiated over evolutionary time in their respective ocean region, potentially forming a sibling species complex. The taxonomic status of the clades is discussed.     

Examination of evolutionary relationships in the Cape fossorial skink species complex (Acontinae: Acontias meleagris meleagris) reveals the presence of five cryptic lineages

Variable morphological characters have obscured genealogical relationships in the Cape fossorial skink Acontias meleagris meleagris species complex. Currently the species complex contains four dubious operational taxonomic units (A. meleagris meleagris, A. m. orientalis, A. percivali tasmani and the morph lineicauda) with poorly defined species boundaries. In the present study we examine the evolutionary relationships within the species complex by sampling 24 skink populations from the known geographical distribution in the Western and Eastern Cape provinces of South Africa, representing a total of 119 specimens. We used partial sequence data derived from two mitochondrial DNA genes, 16S rRNA and COI, and one nuclear DNA gene, intron β‐fibrinogen (β‐fibint 7), to examine evolutionary relationships. Phylogenetic relationships were determined using both Maximum Parsimony (MP) and Bayesian inference (BI) from the combined mtDNA, nDNA and the total evidence data. Additionally we employed Maximum likelihood (ML) analyses on the total evidence data that comprised ～1.5 kb. Topologies derived from the combined mtDNA analyses were congruent with the total evidence analyses (mtDNA + nDNA) and retrieved five major clades with strong statistical support inferred from bootstrapping and posterior probabilities. The five clades were genealogically and geographically exclusive, diagnostic at both the mtDNA and nDNA level and characterized by pronounced sequence divergence, with no shared haplotypes between clades. Collectively these results suggest the presence of five putative cryptic operational taxonomic units within the A. meleagris meleagris species complex. Constraining the traditionally recognized taxa always retrieved a statistically worse topology suggesting that considerable taxonomic revision is required. Our results indicate that traditional morphological characters need to be reassessed to define the five novel lineages in the A. meleagris meleagris species complex. The phylogeographic pattern for the fossorial skinks we retrieved was novel compared to phylogeographic studies for codistributed above ground living taxa. These results suggest that the abiotic and biotic factors that impact subterranean taxa may differ from supraterranean taxa.