Phylogeography of the rare Balkan endemic Martino's vole, Dinaromys bogdanovi, reveals strong differentiation within the western Balkan Peninsula

Phytochemicals and Their Correlation with Molecular Data in Micromeria and Clinopodium (Lamiaceae) Taxa

A study of the phytochemical and molecular characteristics of ten Micromeria and six Clinopodium taxa (family Lamiaceae) distributed in the Balkan Peninsula was carried out. The phytochemicals detected in essential oils by gas chromatography, mass spectrometry, and molecular data amplified fragment length polymorphism were used to study the taxonomic relationships among the taxa and the correlations between phytochemical and molecular data. STRUCTURE analysis revealed three genetic groups, while Bayesian Analysis of Population Structure grouped the studied taxa into 11 clusters nested in the groups obtained by STRUCTURE. Principal components analysis performed with the 21 most represented compounds in the essential oils yielded results that were partly consistent with those obtained by STRUCTURE and neighbour-joining analyses. However, their geographic distributions did not support the genetic grouping of the studied taxa and populations. The Mantel test showed a significant correlation between the phytochemical and genetic data (r = 0.421, p < 0.001). Genetic distance explained 17.8% of the phytochemical distance between populations. The current taxonomic position of several of the studied taxa is yet to be satisfactorily resolved, and further studies are needed. Such future research should include nuclear and plastid DNA sequences from a larger sample of populations and individuals.

The minnow Phoxinus lumaireul (Leuciscidae) shifts the Adriatic-Black Sea basin divide in the north-western Dinaric Karst region

Karst landscapes are characterized by intermittent and sinking streams. The most common method used to study underground hydrological connections in karst is tracing tests. However, a more biologically oriented approach has been suggested: analysis of the genetic structure of aquatic organisms. Biological tracers can be sought among trogloxenes, that is, surface species that occasionally enter caves and groundwater. One such example is the fish genus Phoxinus, which exhibits high genetic diversity and complex phylogeography in the Balkan Peninsula. In the north-western Dinaric Karst, the complex hydrological network was digitalized in 2020. Contemporaneously, Phoxinus lumaireul populations in the Slovenian Dinaric Karst were intensively sampled and analysed for fragments of two mitochondrial genes and one nuclear gene. The derived phylogeographic structure and data on hydrological connections were compared to evaluate support for three alternative scenarios: The genetic structure (1) is a consequence of the ongoing geneflow through underground connections, (2) reflects a previous hydrological network or (3) is an outcome of anthropogenic translocations. The results suggest that the first two scenarios seem to have played a major role, while the third has not had profound effects on the genetic composition. Comparison between the genetic structure of Slovenian Dinaric Karst sampling sites and that of hydrologically isolated reference sampling sites indicated a greater genetic connectivity in the former. Moreover, the range of Adriatic (1a) and Black Sea (1c) haplotypes does not correspond to the Adriatic-Black Sea basin divide but is shifted northwards.

DNA barcoding for biodiversity assessment: Croatian stoneflies (Insecta: Plecoptera)

Background

The hemi-metabolous aquatic order Plecoptera (stoneflies) constitutes an indispensable part of terrestrial and aquatic food webs due to their specific life cycle and habitat requirements. Stoneflies are considered one of the most sensitive groups to environmental changes in freshwater ecosystems and anthropogenic changes have caused range contraction of many species. Given the critical threat to stoneflies, the study of their distribution, morphological variability and genetic diversity should be one of the priorities in conservation biology. However, some aspects about stoneflies, especially a fully resolved phylogeny and their patterns of distribution are not well known. A study that includes comprehensive field research and combines morphological and molecular identification of stoneflies has not been conducted in Croatia so far. Thus, the major aim of this study was to regenerate a comprehensive and taxonomically well-curated DNA barcode database for Croatian stoneflies, to highlight the morphological variability obtained for several species and to elucidate results in light of recent taxonomy.

Methods

A morphological examination of adult specimens was made using basic characteristics for distinguishing species: terminalia in males and females, head and pronotum patterns, penial morphology, and egg structures. DNA barcoding was applied to many specimens to help circumscribe known species, identify cryptic or yet undescribed species, and to construct a preliminary phylogeny for Croatian stoneflies.

Results

Sequences (658 bp in length) of 74 morphospecies from all families present in Croatia were recovered from 87% of the analysed specimens (355 of 410), with one partial sequence of 605 bp in length for Capnopsis schilleri balcanica Zwick, 1984. A total of 84% morphological species could be unambiguously identified using COI sequences. Species delineation methods confirmed the existence of five deeply divergent genetic lineages, with monophyletic origin, which also differ morphologically from their congeners and represent distinct entities. BIN (Barcode Index Number) assignment and species delineation methods clustered COI sequences into different numbers of operational taxonomic units (OTUs). ASAP delimited 76 putative species and achieved a maximum match score with morphology (97%). ABGD resulted in 62 and mPTP in 61 OTUs, indicating a more conservative approach. Most BINs were congruent with traditionally recognized species. Deep intraspecific genetic divergences in some clades highlighted the need for taxonomic revision in several species-complexes and species-groups. Research has yielded the first molecular characterization of nine species, with most having restricted distributions and confirmed the existence of several species which had been declared extinct regionally.

Morphological and molecular characterisation of the Popijač's Yellow Sally, Isoperlapopijaci sp. nov., a new stenoendemic stonefly species from Croatia (Plecoptera, Perlodidae)

A new species of the Yellow Sally genus (Isoperla Banks, 1906) is described, based on morphological (males and females adults, larval and egg) and molecular (the barcode region of the cytochrome c oxidase subunit I gene (COI)) features. Popijač’s Yellow Sally, I.popijaci Hlebec & Sivec, sp. nov. inhabits two karstic sources of the Krasulja rivulet in Croatia. Male and female of the new species are characterised by colouration patterns of the head and pronotum; the dimensions of the female subgenital plate; the medial penial armature and oval-shaped egg without collar and anchor. The larvae differ from their congeners by the uniquely coloured head and pronotum. Based on morphological characteristics I.popijacisp. nov. belongs to the I.tripartita species group. Phylogenetic and taxonomic relationships were reconstructed using three methods of phylogenetic inference and three species delimitation methods. As I.popijacisp. nov. occurs at a narrow area of the Krasulja rivulet in Krbava field, the study puts emphasis on the conservation and hotspot importance of the temporary rivers in the Dinaric karst. Furthermore, the study accentuates the necessity for further research on the genetic diversity of Plecoptera in Croatia.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in &gt; 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Genetic similarities versus morphological resemblance: Unraveling a polyploid complex in a Mediterranean biodiversity hotspot

The Balkan Peninsula is recognized as one of the hotspots of biodiversity in Europe. This area has shown since the Last Glacial Maximum appropriate conditions for species diversification and hybridization, which has led to the existence of numerous taxonomically unresolved entities. Here, we focus on the Western Balkans and explore the genetic structure and relationships among species belonging to the V. austriaca - V. orbiculata diploid-polyploid complex, including populations showing intermediate morphologies. A combination of nuclear markers (microsatellites), plastid DNA regions (trnH-psbA, ycf6-psbM) and ploidy level estimations using flow cytometry are employed to assess the genetic structure and evolutionary dynamics of this polyploid complex. To reconstruct the evolutionary history, an approximate Bayesian computation approach is combined with projections of the species distribution models onto the climatic scenarios of the Mid-Holocene (6 ka BP) and Last Glacial Maximum (22 ka BP). Four main groups were found: one well-established entity within the diploid level, V. dalmatica, a second diploid-tetraploid group which corresponds to V. orbiculata, a hexaploid cluster harboring V. austriaca subsp. jacquinii individuals, and an enigmatic tetraploid group. According to the molecular data obtained, this latter cluster represents an allopolyploid cryptic lineage −with V. orbiculata and V. dalmatica as putative parents− morphologically similar to V. orbiculata, but genetically more related to V. austriaca subsp. jacquinii. Veronica dalmatica and this “uncertain tetraploid” group are involved in the formation of the hexaploid taxon V. austriaca subsp. jacquinii, with the possibility of recent gene flow among different cytotypes. The present study supports a scenario of diversification from a diploid common ancestor leading to two different but interrelated lineages. The first one would correspond with the diploid V. orbiculata plus tetraploid individuals of this species arising through allo- and autopolyploidization, and the second one would involve all ploidy levels with allopolyploidization being prevalent.

Phylogeographical structure of the pygmy shrew: revisiting the roles of southern and northern refugia in Europe

Southern and northern glacial refugia are considered paradigms that explain the complex phylogeographical patterns and processes of European biota. Here, we provide a revisited statistical phylogeographical analysis of the pygmy shrew Sorex minutus Linnaeus, 1766 (Eulipotyphla, Soricidae), examining its genetic diversity, genetic differentiation and demographic history in the Mediterranean peninsulas and in Western and Central Europe. The results showed support for genetically distinct and diverse phylogeographical groups consistent with southern and northern glacial refugia, as expected from previous studies. We also identified geographical barriers concordant with glaciated mountain ranges during the Last Glacial Maximum (LGM), early diversification events dated between the Late Pleistocene and Early Holocene for the main phylogeographical groups, and recent (post-LGM) patterns of demographic expansions. This study is the most comprehensive investigation of this species to date, and the results have implications for the conservation of intraspecific diversity and the preservation of the evolutionary potential of S. minutus.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Small terrestrial mammals of Albania: distribution and diversity (Mammalia, Eulipotyphla, Rodentia)

In this paper new records are reported for 23 species of small terrestrial mammals (STM) of Albania collected during the field work campaigns organised in the framework of the project “Strengthening capacity in National Nature Protection – preparation for Natura 2000 network” (NaturAL) in Albania during the summer and autumn of 2016 and 2017 Data on small mammals were primarily collected through Sherman live-trapping campaigns in six high priority protected areas of Albania: Korab-Koritnik, Bredhi i Hotovës, Tomorri, Llogara-Karaburun, Divjakë-Karavasta, Liqeni i Shkodrës (Skadar lake), Lëpushë-Vermosh. Other data were obtained by analysis of owl pellets or by direct observation of individuals (dead or alive) in the field. For 21 species Erinaceusroumanicus, Neomysanomalus, Crocidurasuaveolens, Crociduraleucodon, Suncusetruscus, Talpastankovici/caeca, Myocastorcoypus, Sciurusvulgaris, Glisglis, Dryomysnitedula, Muscardinusavellanarius, Microtuslevis/arvalis, Microtussubterraneus, Microtusthomasi, Microtusfelteni, Myodesglareolus, Apodemussylvaticus, Apodemusflavicollis, Apodemusepimelas, Musmusculus, and Musmacedonicus additional records are provide and their distributions reviewed, while the presence of two new species of shrews (Sorexaraneus and Sorexminutus) for Albania is reported for the first time. A comprehensive review of the published and unpublished distribution records of STM species of the country is made, together with an updated checklist and distribution maps of the species.

Unique postglacial evolution of the hornbeam (Carpinus betulus L.) in the Carpathians and the Balkan Peninsula revealed by chloroplast DNA

The Balkan Peninsula is one of the largest and most important European glacial refugia. However, the evolutionary history and phylogeographic pattern of temperate tree species that survived in the Balkans glacial refugia and their contribution to the genetic structure of the current population in the Carpathian Mountains remains poorly understood. Using polymerase chain reaction - restriction fragment length polymorphism (PCR-RFLP), and extensive population sampling, we explored the phylogeographic pattern of Carpinus betulus in both the Balkan Peninsula and the Carpathian region. We aimed to determine the locations of potential glacial refugia, in order to delineate post-glacial colonization routes, and also to test if northern cryptic refugia had persisted during the Last Glacial Maximum (LGM). Our results provide strong support for the existence of multiple refugia similar to the 'refugia-within-refugia' scenario, which would suggest that Carpinus betulus has experienced a complex evolutionary history. In agreement with pollen data, our findings corroborate with previous hypotheses suggesting that hornbeam has a distinct postglacial evolution in the southeast of Europe (Carpathians and Balkan Peninsula) compared with that in Western Europe. Three postglacial re-colonization routes, from three distinct effective glacial refugia: (1) the Dinaric Alps; (2) the Pirin and Rhodope Mountains and (3) the Strandzha Mountains, were detected within the Balkan Peninsula. The pattern of cpDNA haplotypes distribution across the Ukrainian Carpathians revealed a "suture zone", which is a consequence of contact due to postglacial re-colonization between hornbeam populations originating from all three effective glacial refugia. The peculiar haplotype geographical structure found in C. betulus and the identified multiple glacial refugia in the Balkans may have direct consequences in the management and conservation of hornbeam forest genetic resources.

The Small and the Dead: A Review of Ancient DNA Studies Analysing Micromammal Species

The field of ancient DNA (aDNA) has recently been in a state of exponential growth, largely driven by the uptake of Next Generation Sequencing (NGS) techniques. Much of this work has focused on the mammalian megafauna and ancient humans, with comparatively less studies looking at micromammal fauna, despite the potential of these species in testing evolutionary, environmental and taxonomic theories. Several factors make micromammal fauna ideally suited for aDNA extraction and sequencing. Micromammal subfossil assemblages often include the large number of individuals appropriate for population level analyses, and, furthermore, the assemblages are frequently found in cave sites where the constant temperature and sheltered environment provide favourable conditions for DNA preservation. This review looks at studies that include the use of aDNA in molecular analysis of micromammal fauna, in order to examine the wide array of questions that can be answered in the study of small mammals using new palaeogenetic techniques. This study highlights the bias in current aDNA studies and assesses the future use of aDNA as a tool for the study of micromammal fauna.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.

Morphological, genetic and phytochemical variation of the endemic Teucrium arduini L. (Lamiaceae)

Analysis of the morphological traits of leaves, genetic variability (analyzed by AFLP) and chemical composition of essential oils (analyzed by GC-MS) was conducted on eleven populations of the endemic Illyric-Balcanic species Teucrium arduini L. in Croatia, Bosnia and Herzegovina, and Montenegro. Average blade length and width ranged from 20.00 to 31.47mm and from 11.58 to 15.66mm, respectively. Multivariate analysis (PCA, UPGMA) of morphological traits distinguished two continental Bosnian populations from the remaining populations. AFLP analysis separated the investigated populations into two groups based primarily on geographical distance. Essential oil analysis showed a total of 52 compounds, with two chemotypes distinguished based on the essential oil profile. The first was a sesquiterpene chemotype, with β-caryophyllene, germacrene D or caryophyllene oxide as the major compounds, while the second was an oxygenated monoterpene chemotype, with pulegone and piperitone oxide as the main components. The Mantel test showed a stronger correlation between the morphological traits and AFLP than between the essential oil profile and AFLP. The test also showed a stronger association between the essential oil profile and geographical position than between the morphological traits and geographical position.

Genetic structure and historical demography of Schizothorax nukiangensis (Cyprinidae) in continuous habitat

Geographic distance, different living habitats or Pleistocene climatic oscillations have frequently been found to shape population genetic structure in many species. The genetic structure of Schizothorax nukiangensis, a high altitude, valuable fish species, which is distributed throughout the Nujiang River, was investigated by mitochondrial DNA sequence analysis. The cytochrome c oxidase subunit I (COI), cytochrome b (cytb), and the mitochondrial control region (MCR) of S. nukiangensis were concatenated for examination of population structure and demographic history. The concatenated data set (2405 bp) implied a pronounced genetic population structure (overall F_ST = 0.149) and defined two population units. Strong differentiation was detected between the Sanjiangkou (SJK) population and other populations due to environmental heterogeneity, dispersal ability, and/or glacial cycles. Additional DNA sequencing of the nuclear RAG2 gene also examined significant differentiation between two units and between SJK and the upstream populations (U-unit). Recent expansion events suggest that S. nukiangensis may have undergone a rapid increase during warm interglacial periods. Surprisingly, S. nukiangensis appears to have undergone an obvious expansion during the last glaciations (LG) for cold hardiness and a sharp contraction from 1.5 ka to the present. However, two population units exhibited different reflections during the LG, which might be closely related to their living habitats and cold hardiness. A clear pattern of isolation by distance was detected in S. nukiangensis due to feeding habits, limited dispersal ability, and/or philopatry. It is vitally important that more attention be given to S. nukiangensis due to low genetic diversity, lack of gene flow, and recent population contraction.

Genetic diversity and structure of Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir./Sch./Bip., Asteraceae) within the Balkan refugium

Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir. /Sch./ Bip.) is an outcrossing, perennial insecticidal plant, restricted to the eastern Adriatic coast (Mediterranean). Amplified fragment-length polymorphisms (AFLP) were used to investigate the genetic diversity and structure within and among 20 natural plant populations. The highest level of gene diversity, the number of private alleles and the frequency down-weighted marker values (DW) were found in northern Adriatic populations and gradually decreased towards the southern boundary of the species range. Genetic impoverishment of these southern populations is most likely the result of human-related activities. An analysis of molecular variance (AMOVA) indicated that most of the genetic diversity was attributed to differences among individuals within populations (85.78%), which are expected due to the outcrossing nature of the species. A Bayesian analysis of the population structure identified two dominant genetic clusters. A spatial analysis of the genetic diversity indicated that 5.6% of the genetic differentiation resulted from isolation by distance (IBD), while 12.3% of the genetic differentiation among populations followed the pattern of isolation by environmental distance (IBED). Knowledge of the genetic diversity patterns of the natural populations and the mechanism behind these patterns is required for the exploitation and possible conservation management of this endemic and economically important species.

Testing the efficiency of nested barriers to dispersal in the Mediterranean high mountain plant Edraianthus graminifolius (Campanulaceae)

Due to strong spatial heterogeneity and limited Pleistocene glaciation, the Balkan Peninsula is a major European biodiversity hot spot. Surprisingly little, however, is known about patterns and processes of intraspecific diversification of its biota in general and of high-altitude species in particular. A well-suited system to test hypotheses with respect to various isolating factors acting at different geographic scales and to explore full-range phylogeographic patterns on the Balkan Peninsula is Edraianthus graminifolius (Campanulaceae), distributed in the western Balkan mountain systems, the southwestern Carpathians and the Apennine Peninsula. To this end, we used a dense population sampling and employed amplified fragment length polymorphism (AFLP) markers and plastid DNA sequences supplemented by ecological niche modelling. The strongest splits were inferred to separate southern and northern Balkan populations from the central ones, from where range extension occurred to the Carpathians and, in more recent times, once or twice to the Apennine Peninsula. The three genetic groups in the western Balkan Peninsula were remarkably congruent among molecular markers, suggesting that the barriers to gene flow acted over long time periods facilitating allopatric differentiation. Each main group of Balkan populations contained genetically and geographically distinct subgroups, which likely are the result of local refugia during warmer periods. Evidently, the topographically highly complex and during the Last Glacial Maximum only locally glaciated Balkan Peninsula is a hot spot of species richness and endemism as well as a sanctuary of intraspecific genetic diversity, even if the underlying causes remain insufficiently understood.

Microscale vicariance and diversification of Western Balkan caddisflies linked to karstification

The karst areas in the Dinaric region of the Western Balkan Peninsula are a hotspot of freshwater biodiversity. Many investigators have examined diversification of the subterranean freshwater fauna in these karst systems. However, diversification of surface-water fauna remains largely unexplored. We assessed local and regional diversification of surface-water species in karst systems and asked whether patterns of population differentiation could be explained by dispersal-diversification processes or allopatric diversification following karst-related microscale vicariance. We analyzed mitochondrial cytochrome c oxidase subunit I (mtCOI) sequence data of 4 caddisfly species (genus Drusus) in a phylogeographic framework to assess local and regional population genetic structure and Pliocene/Pleistocene history. We used BEAST software to assess the timing of intraspecific diversification of the target species. We compared climate envelopes of the study species and projected climatically suitable areas during the last glacial maximum (LGM) to assess differences in the species climatic niches and infer potential LGM refugia. The haplotype distribution of the 4 species (324 individuals from 32 populations) was characterized by strong genetic differentiation with few haplotypes shared among populations (16%) and deep divergence among populations of the 3 endemic species, even at local scales. Divergence among local populations of endemics often exceeded divergence among regional and continental clades of the widespread D. discolor. Major divergences among regional populations dated to 2.0 to 0.5 Mya. Species distribution model projections and genetic structure suggest that the endemic species persisted in situ and diversified locally throughout multiple Pleistocene climate cycles. The pattern for D. discolor was different and consistent with multiple invasions into the region. Patterns of population genetic structure and diversification were similar for the 3 regional endemic Drusus species and consistent with microscale vicariance after the onset of intensified karstification in the Dinaric region. Karstification may induce microscale vicariance of running surface-water habitats and probably promotes allopatric fragmentation of stream insects at small spatial scales.

Combined mitochondrial and nuclear markers revealed a deep vicariant history for Leopoldamys neilli, a cave-dwelling rodent of Thailand

BACKGROUND

Historical biogeography and evolutionary processes of cave taxa have been widely studied in temperate regions. However, Southeast Asian cave ecosystems remain largely unexplored despite their high scientific interest. Here we studied the phylogeography of Leopoldamys neilli, a cave-dwelling murine rodent living in limestone karsts of Thailand, and compared the molecular signature of mitochondrial and nuclear markers.

METHODOLOGY/PRINCIPAL FINDINGS

We used a large sampling (n = 225) from 28 localities in Thailand and a combination of mitochondrial and nuclear markers with various evolutionary rates (two intronic regions and 12 microsatellites). The evolutionary history of L. neilli and the relative role of vicariance and dispersal were investigated using ancestral range reconstruction analysis and Approximate Bayesian computation (ABC). Both mitochondrial and nuclear markers support a large-scale population structure of four main groups (west, centre, north and northeast) and a strong finer structure within each of these groups. A deep genealogical divergence among geographically close lineages is observed and denotes a high population fragmentation. Our findings suggest that the current phylogeographic pattern of this species results from the fragmentation of a widespread ancestral population and that vicariance has played a significant role in the evolutionary history of L. neilli. These deep vicariant events that occurred during Plio-Pleistocene are related to the formation of the Central Plain of Thailand. Consequently, the western, central, northern and northeastern groups of populations were historically isolated and should be considered as four distinct Evolutionarily Significant Units (ESUs).

CONCLUSIONS/SIGNIFICANCE

Our study confirms the benefit of using several independent genetic markers to obtain a comprehensive and reliable picture of L. neilli evolutionary history at different levels of resolution. The complex genetic structure of Leopoldamys neilli is supported by congruent mitochondrial and nuclear markers and has been influenced by the geological history of Thailand during Plio-Pleistocene.

Fish migrate underground: the example of Delminichthys adspersus (Cyprinidae)

Complex aquatic systems of karst harbour a rich but little-investigated biodiversity. In Croatia and Bosnia-Herzegovina karst, temporal springs are inhabited by a group of minnow-like fishes that retreat to the associated ground water during dry seasons and spend several months underground. The most abundant species in this group is Delminichthys adspersus (Heckel 1843), which also has the most fragmented distribution range. To determine the population composition and dispersal patterns, and to detect potential underground migration, a large genetic data set comprising 544 specimens of D. adspersus covering most of its distribution area was analysed. Analysis of mitochondrial cytochrome b sequences (∼1000 bp) and eight microsatellite loci showed that D. adspersus comprises at least three subpopulations with gene flow occurring among them. Coalescent-based analysis revealed a complex migration pattern, with several unidirectional dispersal paths, including between temporal springs that share no surface connection. The results of this study suggest the existence of recurrent underground migration of fish in a karst environment and demonstrate the complexity of its hydrological network. The findings are relevant to conservation strategies for endemic karst organisms and karst ecosystems as a whole.

Nondestructive DNA extraction from museum specimens

Natural history museums around the world hold millions of animal and plant specimens that are potentially amenable to genetic analyses. With more and more populations and species becoming extinct, the importance of these specimens for phylogenetic and phylogeographic analyses is rapidly increasing. However, as most DNA extraction methods damage the specimens, nondestructive extraction methods are useful to balance the demands of molecular biologists, morphologists, and museum curators. Here, I describe a method for nondestructive DNA extraction from bony specimens (i.e., bones and teeth). In this method, the specimens are soaked in extraction buffer, and DNA is then purified from the soaking solution using adsorption to silica. The method reliably yields mitochondrial and often also nuclear DNA. The method has been adapted to DNA extraction from other types of specimens such as arthropods.

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.

Genetic differentiation of the Western Capercaillie highlights the importance of South-eastern Europe for understanding the species phylogeography

The Western Capercaillie (Tetrao urogallus L.) is a grouse species of open boreal or high altitude forests of Eurasia. It is endangered throughout most mountain range habitat areas in Europe. Two major genetically identifiable lineages of Western Capercaillie have been described to date: the southern lineage at the species' southernmost range of distribution in Europe, and the boreal lineage. We address the question of genetic differentiation of capercaillie populations from the Rhodope and Rila Mountains in Bulgaria, across the Dinaric Mountains to the Slovenian Alps. The two lineages' contact zone and resulting conservation strategies in this so-far understudied area of distribution have not been previously determined. The results of analysis of mitochondrial DNA control region sequences of 319 samples from the studied populations show that Alpine populations were composed exclusively of boreal lineage; Dinaric populations of both, but predominantly (96%) of boreal lineage; and Rhodope-Rila populations predominantly (>90%) of southern lineage individuals. The Bulgarian mountains were identified as the core area of the southern lineage, and the Dinaric Mountains as the western contact zone between both lineages in the Balkans. Bulgarian populations appeared genetically distinct from Alpine and Dinaric populations and exhibited characteristics of a long-term stationary population, suggesting that they should be considered as a glacial relict and probably a distinct subspecies. Although all of the studied populations suffered a decline in the past, the significantly lower level of genetic diversity when compared with the neighbouring Alpine and Bulgarian populations suggests that the isolated Dinaric capercaillie is particularly vulnerable to continuing population decline. The results are discussed in the context of conservation of the species in the Balkans, its principal threats and legal protection status. Potential conservation strategies should consider the existence of the two lineages and their vulnerable Dinaric contact zone and support the specificities of the populations.

Phylogeography of the endemic Gymnocypris chilianensis (Cyprinidae): sequential westward colonization followed by allopatric evolution in response to cyclical Pleistocene glaciations on the Tibetan Plateau

The schizothoracine Gymnocypris chilianensis is restricted to the Shiyang, Ruoshui and Shule Rivers, listed from east to west, along the northeast edge of the Tibetan Plateau. This distribution provides a valuable system to test hypotheses about postglacial colonization. We used mitochondrial DNA sequence data (a control region and the cytochrome b gene; 1894 bp) to assess the phylogeographic structure of this species based on 278 specimens sampled from throughout the species' entire geographical range. We found three lineages corresponding geographically to the three rivers, suggesting three independent glacial differentiation centers within the northeast edge of the Tibetan Plateau. The phylogenetic analysis suggested that the Shiyang River population forms a lineage that separated from the other populations of G. chilianensis at the basal phylogenetic split within this species. The molecular data further demonstrated a clear pattern of decreasing genetic diversity from the eastern Shiyang River towards the central Ruoshui River and western Shule River lineages, a pattern consistent with sequential western colonization. We therefore propose a phylogeographic scenario for G. chilianensis of a gradual westerly expansion from the Shiyang River population along the northeast edge of the Tibetan Plateau, with subsequent allopatric evolution at approximately 0.37 and 0.05 million years ago (Ma), through at least two glacial maxima. Together with the genetic evidence reported in other species, our findings suggest that this common biogeographic pattern emphasizes the importance of the northeastern edge region of the Tibetan Plateau as a hotspot of genetic diversity for some taxa.

Disentangling phylogeography, polyploid evolution and taxonomy of a woodland herb (Veronica chamaedrys group, Plantaginaceae s.l.) in southeastern Europe

Southeastern Europe is a centre of European biodiversity, but very little is known about factors causing the observed richness. Here, we contribute to fill this gap by reconstructing the spatio-temporal diversification of the cytologically variable and taxonomically intricate complex of Veronica chamaedrys (Plantaginaceae s.l.), growing in open forests, forest edges and grasslands, with flow cytometry, molecular markers (AFLPs, plastid DNA sequences) and morphometry. Our results show that both diploid and tetraploid cytotypes are widespread, but diploids predominate on the southern Balkan Peninsula. Plastid sequences suggest a first split into three main lineages in the mid-Pleistocene and a continuous diversification during the last 0.4 my. Two of the identified plastid lineages coincide with geographically distinct AFLP clusters. Altogether, the genetic data suggest forest refugia on the southern-most Balkan Peninsula (Greece), in Bulgaria, Istria (Croatia and Slovenia) and maybe the southeastern Carpathians (Romania). Morphometric and genetic data show little congruence with current taxonomy.

A rapid column-based ancient DNA extraction method for increased sample throughput

Genetic analyses using museum specimens and ancient DNA from fossil samples are becoming increasingly important in phylogenetic and especially population genetic studies. Recent progress in ancient DNA sequencing technologies has substantially increased DNA sequence yields and, in combination with barcoding methods, has enabled large-scale studies using any type of DNA. Moreover, more and more studies now use nuclear DNA sequences in addition to mitochondrial ones. Unfortunately, nuclear DNA is, due to its much lower copy number in living cells compared to mitochondrial DNA, much more difficult to obtain from low-quality samples. Therefore, a DNA extraction method that optimizes DNA yields from low-quality samples and at the same time allows processing many samples within a short time frame is immediately required. In fact, the major bottleneck in the analysis process using samples containing low amounts of degraded DNA now lies in the extraction of samples, as column-based methods using commercial kits are fast but have proven to give very low yields, while more efficient methods are generally very time-consuming. Here, we present a method that combines the high DNA yield of batch-based silica extraction with the time-efficiency of column-based methods. Our results on Pleistocene cave bear samples show that DNA yields are quantitatively comparable, and in fact even slightly better than with silica batch extraction, while at the same time the number of samples that can conveniently be processed in parallel increases and both bench time and costs decrease using this method. Thus, this method is suited for harvesting the power of high-throughput sequencing using the DNA preserved in the millions of paleontological and museums specimens.

Shallow phylogeographic structuring of Vimba vimba across Europe suggests two distinct refugia during the last glaciation

Genetic variation and geographical structuring of vimba Vimba vimba were analysed across 26 sites (80 individuals) by means of mtDNA sequences (cyt b gene, mitochondrial control region) to localize hypothesized glacial refugia and to reconstruct postglacial recoloniation routes. Although genetic diversity among sequenced individuals was low, a combined analysis of the two sequenced fragments revealed a western (central and northern Europe: Danube, Elbe and lakes of Sweden) and an eastern clade (eastern Europe: Dnieper-South Bug, Don, Neman). Furthermore, a number of divergent ancestral haplotypes distributed around the Black and Caspian Seas became apparent. Mismatch analyses supported a sudden expansion model for the populations of the western clade between 50 and 10 000 bp. Overall, the study provides strong evidence for a northward and westward expansion of V. vimba from two refugial regions located in the Danubian drainage and the northern Pontic regions respectively.