Comparative phylogeography of the five Greek vole species infers

This article was downloaded by: [University of Patras]
On: 12 September 2012, At: 01:52
Publisher: Taylor & Francis
Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House,
37-41 Mortimer Street, London W1T 3JH, UK
Italian Journal of Zoology
Publication details, including instructions for authors and subscription information:
http://www.tandfonline.com/loi/tizo20
Comparative phylogeography of the five Greek vole
species infers the existence of multiple South Balkan
subrefugia
a
a
a
E. Thanou , G. Tryfonopoulos , B. Chondropoulos & S. Fraguedakis-Tsolis
a
a
Section of Animal Biology, Department of Biology, University of Patras, Patras, Greece
Version of record first published: 07 Feb 2012.
To cite this article: E. Thanou, G. Tryfonopoulos, B. Chondropoulos & S. Fraguedakis-Tsolis (2012): Comparative
phylogeography of the five Greek vole species infers the existence of multiple South Balkan subrefugia, Italian Journal of
Zoology, 79:3, 363-376
To link to this article: http://dx.doi.org/10.1080/11250003.2011.651163
PLEASE SCROLL DOWN FOR ARTICLE
Full terms and conditions of use: http://www.tandfonline.com/page/terms-and-conditions
This article may be used for research, teaching, and private study purposes. Any substantial or systematic
reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to
anyone is expressly forbidden.
The publisher does not give any warranty express or implied or make any representation that the contents
will be complete or accurate or up to date. The accuracy of any instructions, formulae, and drug doses should
be independently verified with primary sources. The publisher shall not be liable for any loss, actions, claims,
proceedings, demand, or costs or damages whatsoever or howsoever caused arising directly or indirectly in
connection with or arising out of the use of this material.
Italian Journal of Zoology, September 2012; 79(3): 363–376
Comparative phylogeography of the five Greek vole species infers
the existence of multiple South Balkan subrefugia
E. THANOU*, G. TRYFONOPOULOS, B. CHONDROPOULOS,
& S. FRAGUEDAKIS-TSOLIS
Section of Animal Biology, Department of Biology, University of Patras, Patras, Greece
Downloaded by [University of Patras] at 01:52 12 September 2012
(Received 15 May 2011; accepted 6 December 2011)
Abstract
Despite extensive phylogenetic studies in the genus Microtus, several of its species have not been thoroughly evaluated.
The present study focuses on a cytb molecular analysis of five vole species occurring in the southern Balkan Peninsula
(M. felteni, M. thomasi, M. subterraneus, M. levis and M. guentheri), aiming to demonstrate the importance of thorough
intraspecific sampling when the phylogeny of closely related taxa is tested. As a result, the Balkan populations of these
voles showed significant intraspecific variation that distinguishes them from other European and Asian conspecifics and in
some cases reveals distinct lineages even within the Balkan region. Their complex phylogeography suggest the existence of
multiple subrefugia, located within the southern Balkan region, which promoted diversification of these voles during the
Middle and Late Pleistocene glacial periods. The significant role of the southern Balkan Peninsula’s paleogeographical and
paleoclimatical characteristics in these small mammals’ evolution is discussed.
Keywords: Mediterranean Peninsula, Balkan subrefugia, Microtus species, glacial/interglacial cycles, paleogeography
Introduction
Populations distributed at the borders of a species’
distributional range, as those of the Mediterranean
Peninsulas, can easily be isolated or decreased in
abundance, and might show morphological and ecological adaptations because of genetic isolation and
local climatic conditions. Many of the species currently distributed in these peninsulas have experienced severe events (population fragmentation, bottleneck, etc) during the Pleistocene climatic cycles.
Their intraspecific phylogenies imply these events but
also reveal the existence of several southern refugia in the Mediterranean area from which isolated
populations recolonized northern Europe (Hewitt
2001). The Balkan Peninsula, being at the southeastern margin of Europe, is a natural border for
European populations but also a biogeographical
bridge between Asia and Europe, and bears specific
paleogeographical and paleoclimatical characteristics that could play an important role to species’
evolution. Especially during the Pleistocene, the geomorphology of this area was repeatedly altered by
eustatic movements of the sea level that were driven
by climatic change, from cold/dry to warm/humid
conditions, through the glacial/interglacial periods
(Perissoratis & Conispoliatis 2003).
The Mediterranean Peninsulas have acted as
“speciation traps”, concentrating endemics, but also
promoted intraspecific divergence. For example,
Mediterranean populations of the pygmy shrew and
the bank vole (Bilton et al. 1998; Deffontaine et al.
2005; Vega et al. 2010) show little similarity with
their conspecifics further north as well as among
each other, suggesting that the current phylogeographical profile of these species may be explained
by the assumption of multiple Mediterranean refugia. The Italian and the Balkan Peninsulas have
provided many examples of genetically differentiated
lineages within several widely distributed species,
such as the yellow-necked fieldmouse (Michaux et al.
2005), the lesser white-toothed shrew (Dubey et al.
2007a), the European ground squirrel (Kryštufek
et al. 2009a) and the common vole (Bužan et al.
2010). In the Balkan Peninsula and particularly its
*Correspondence: E. Thanou, Section of Animal Biology, Department of Biology, University of Patras, GR 26501 Patras, Greece. Tel: +30 2610 969237. Fax:
+30 2610 969218. Email: [email protected]
ISSN 1125-0003 print/ISSN 1748-5851 online © 2012 Unione Zoologica Italiana
http://dx.doi.org/10.1080/11250003.2011.651163
Downloaded by [University of Patras] at 01:52 12 September 2012
364
E. Thanou et al.
southernmost areas, such as Greece, many rodents
show higher morphological and/or genetic variation, when Greek populations are compared to
their northern European conspecifics. Such examples have been presented in immunological studies on Spermophilus citellus (Fraguedakis-Tsolis &
Ondrias 1985) and molecular ones on Mus spicilegus
(Mitsainas et al. 2009).
A recent study (Fraguedakis-Tsolis et al. 2009)
revealed distinct morphological profiles in Greek
voles: Greek M. guentheri, M. levis and M. thomasi
are distinctively smaller in size and differ in skull
measures from other European populations, while
high morphological variation exists among Greek
populations of M. thomasi, as well as M. subterraneus. Moreover, a molecular analysis on Greek
M. thomasi populations (Tryfonopoulos et al. 2008)
showed a high degree of intraspecific differentiation
and argued in favour of a new cryptic subspecies distributed in SE Central Greece, namely M. thomasi
atticus.
Recently, the phylogenetic relationships between
all currently recognized European Microtus species
were inferred, with the analysis of the cytochrome
b (cytb) mitochondrial marker (Jaarola et al. 2004).
This analysis focused mainly on the relationships
between subgenera and species, including adequate
specimens to represent all the respective taxa but with
relatively few individuals and populations per species.
This resulted in a poor representation of geographical
intraspecific variation, especially from the Balkans.
Even so, their results also revealed a wide range of
molecular variation within European species of the
genus, some of them being much more polymorphic
than others, in concordance with previous studies of American Microtus species (Conroy & Cook
2000).
It is widely accepted that the genus Microtus,
despite its recent evolutionary history and morphological homogeneity (Musser & Carleton 2005)
probably includes many hitherto unidentified cryptic
species (Graf 1982; Hellborg et al. 2005) that only
detailed within-species sampling could uncover. This
is shown in extensive phylogeographic surveys on
several vole species that revealed intraspecific diversification and complicated phylogenetic relationships
between isolated populations (Hellborg et al. 2005;
Piertney et al. 2005; Tougard et al. 2008). Detailed
within-species sampling and subsequent molecular
analyses have also unraveled a complex pattern of
intraspecific divergence that seems to be related
to the role of Mediterranean Peninsulas as glacial
refugia, but mostly to the existence of subrefugia
within them (Castiglia et al. 2008; Centeno-Cuadros
et al. 2009). Such phylogeographic patterns within
refugia are only recently addressed in Iberia and Italy,
whether respective data for the Balkans are still very
limited and refer to the northern parts of the Balkan
Peninsula (Bužan et al. 2010).
In this context, we focused on the five vole
species living in Greece: Microtus guentheri (Danford
& Alston 1880), M. levis Miller 1908 (previously
named M. rossiaemeridionalis Ognev 1924, see Wilson
& Reeder (2005) for nomenclature review of the
species), M. subterraneus (de Sélys Longchamps
1836), M. felteni (Malec & Storch 1963) and
M. thomasi Barrett-Hamilton 1903. The latter two
species are endemics to the SW Balkans, while all of
them show their southernmost European distribution
limit in Greece (Mitchell-Jones et al. 1999).
Cytb analysis was used in this study to measure the
genetic variation among Greek populations of each
of the studied species. Our data were subsequently
compared to those available in the literature, aiming
to investigate the differentiation of Greek voles compared to previously studied populations from other
areas of each species’ distributional range. In this
sense, we aim to explore the role of the Balkan
Peninsula as a refugium for ancestral vole populations and the possible existence of subrefugia from
which subsequent colonization of northern parts of
Europe has occurred. During this study, extensive
sampling was considered as a priority in order to
reveal within-species cryptic variation and provide
new data on each species’ geographical distribution
in Greece, which is not yet sufficiently known.
Materials and methods
Sampling
A total of 99 specimens, belonging to the five
Microtus species distributed in Greece, were used
(Figure 1, Table I). M. thomasi, was represented in
this study both as the nominal subspecies and as the
subpecies M. thomasi atticus. Twenty one previously
published sequences of this species (Tryfonopoulos
et al. 2008) were updated by adding one M. thomasi
thomasi from Albania, and five M. thomasi atticus from new localities which extend its previously known distribution. Moreover, five M. felteni,
16 M. subterraneus, 20 M. levis and 31 M. guentheri were sequenced. Additional sequences retrieved
from GenBank were also included in the analyses (Table I), representing populations of these five
Microtus species from other European locations, as
well as 19 other Microtus species (Table I), in order
to fully represent all the Eurasian voles that are considered closely related to the five species studied
(Jaarola et al. 2004; Kryštufek et al. 2009b). In this
Downloaded by [University of Patras] at 01:52 12 September 2012
Vole phylogeography infers multiple Balkan subrefugia
365
Figure 1. Approximate distribution of the five Microtus species living in Greece. Collection sites of individuals sequenced in this study and
their respective numbers are given in Table I.
sense, our analysis includes all recognized species
and all distinct mtDNA lineages placed within the
“arvalis” (M. arvalis “obscurus”, M. arvalis “arvalis”
and M. levis) and the “socialis” species groups
(M. socialis, M. irani, M. anatolicus, M. dogramacii,
and the “west” and “south” lineages of M. guentheri) and most of the European endemics, according to these previous studies. Chionomis nivalis was
considered to be the best fitting outgroup since
it represents the basal clade to all the Palearctic
species belonging to the genus Microtus (Jaarola et al.
2004).
DNA extraction, PCR amplification and sequencing
DNA was extracted from liver tissue, preserved
®
at -80 ◦ C, using the Invisorb Spin Tissue kit
(Invitek, Germany). A segment of the mitochondrial
cytb gene was PCR-amplified using the set of
primers and PCR conditions previously described in
Tryfonopoulos et al. (2008).
PCR products were visualized via agarose gel
®
electrophoresis, purified with the Invisorb Spin
PCRapid kit (Invitek, Germany) and the light strand
®
was sequenced on an ABI PRISM 3100 capillary
sequencer (VBC Biotech, Austria) using the primers
of the amplification procedure.
All sequences obtained in this study have been
deposited in GenBank and the respective Accession
Numbers are given in Table I.
Sequence data analyses
Sequences were aligned using CLUSTAL-W v1.4
(Thompson et al. 1994) in the BIOEDIT Sequence
Alignment Editor v5.0.9 (Hall 1999). Genetic divergence (uncorrected p-distances) was estimated in
MEGA v3.1 (Kumar et al. 2004), as the number of nucleotide substitutions per site between and
within populations, as well as net between populations’ genetic distances, after excluding the respective within - values of each of the two populations compared. Populations of each species were
geographically grouped, according to their Greek,
European or Asian origin, and net p-distances among
these groups were also calculated. The saturation
of phylogenetic information was examined using the
Xia’s test employed in DAmBE v4.5.40 (Xia 2000;
Xia & Xie 2001).
MODELTEST v3.06 (Posada & Crandall 1998)
was used to determine the best substitution model
for the Bayesian Inference (BI) analysis and
Maximum Likelihood (ML) analysis, under the
Akaike Information Criterion (AIC). Subsequently,
BI analysis was implemented using MrBayes
v3.1 (Ronquist & Huelsenbeck 2003), under the
GTR+I+G model. Four chains were run for three
million generations, sampling trees every 100th generation and “burn-in” data representing 25% of
early generations (7.5x104 trees) were discarded.
Three simultaneous and independent runs were conducted and inspected for consistency to check for
366
E. Thanou et al.
Table I. Microtus species, number of specimens, locations (site and country), codes corresponding to all sequences included in the present
study as they appear in the respective tree (Figure 2) and GenBank accession numbers. Sequences retrieved from GenBank, are marked
with an asterisk (∗ ). Localities in bold refer to new findings for the respective species.
Species
Site (as in Figure 1)
Country
Code
Ass. Nos
3
Agios Stefanos (1)
Greece
2
Afidnes (2)
Greece
1
1
1
1
1
3
Marathonas (3)
Fyli (4)
Oropos (5)
Avlona (6)
Skourta (7)
Styra, Evvoia (8)
Greece
Greece
Greece
Greece
Greece
Greece
2
Pogontas, Evvoia (9)
Greece
3
Lake Dystos, Evvoia (10)
Greece
2
Oraioi, Evvoia (11)
Greece
M. thomasi thomasi
1
1
1
1
1
1
1
M. felteni
1
3
Mt. Panachaiko (12)
Astakos (13)
Elati (14)
Balbouma (15)
Mikrokastro (16)
Mt Pieria (17)
Sarantë
Trebinje
Seli (31)
Pisoderi (32)
Greece
Greece
Greece
Greece
Greece
Greece
Albania
Bosnia
Greece
Greece
1
1
4
Nymfaio (38)
Begova Češma, Mt Pelister
Seli (31)
Greece
FYROM
Greece
1
4
Pisoderi (32)
Kato Vermio (33)
Greece
Greece
1
3
Dobro Pole, Mt Vorras (34)
Mt Pinovo (35)
Greece
Greece
1
2
Lailias, Mt Vrontous (36)
Mt Silo (37)
Greece
Greece
1
1
1
1
4
Val Piora
Pilatus
Çiğlikara
Güzyurdu
Kato Nevrokopi (27)
Switzerland
Switzerland
Turkey
Turkey
Greece
GR1.1
GR1.2
GR1.3
GR2.1
GR2.2
GR3
GR4
GR5
GR6
GR7
GR8.1
GR8.2
GR8.3
GR9.1
GR9.2
GR10.1
GR10.2
GR10.3
GR11.1
GR11.2
GR12
GR13
GR14
GR15
GR16
GR17
ALBANIA
BOSNIA
GR31.1
GR32.1
GR32.2
GR32.3
GR38
FYROM
GR31.2
GR31.3
GR31.4
GR31.5
GR32.4
GR33.1
GR33.2
GR33.3
GR33.4
GR34
GR35.1
GR35.2
GR35.3
GR36
GR37.1
GR37.2
SWITZERLAND1
SWITZERLAND2
TURKEY1
TURKEY2
GR27.1
GR27.2
GR27.3
GR27.4
EF666504∗
EF666505∗
EF666506∗
EF666510∗
EF666511∗
JN650051
JN650052
JN650053
JN650054
JN674457
EF666590∗
EF666591∗
EF666592∗
EF666588∗
EF666589∗
EF666585∗
EF666586∗
EF666587∗
EF666583∗
EF666584∗
EF666524∗
EF666536∗
EF666522∗
EF666497∗
EF666518∗
EF666949∗
JN650055
AY513844∗
FJ641181
FJ641178
FJ641179
FJ641180
JN650050
AY513798∗
FJ641162
FJ641163
FJ641164
FJ641165
FJ641173
FJ641166
FJ641167
FJ641168
FJ641169
FJ641174
FJ641179
FJ641171
FJ641172
FJ641175
FJ641176
FJ641177
AJ717745∗
AY332714∗
AY513834∗
AY513836∗
FJ641158
FJ641159
FJ641160
FJ641161
Downloaded by [University of Patras] at 01:52 12 September 2012
Microtus thomasi
atticus
M. subterraneus
M. levis
No
(Continued)
Vole phylogeography infers multiple Balkan subrefugia
367
Table I. (Continued).
Downloaded by [University of Patras] at 01:52 12 September 2012
Species
M. guentheri “west”
No
Site (as in Figure 1)
Country
Code
Ass. Nos
1
4
Olynthos (28)
Sevastiana (29)
Greece
Greece
4
Arnissa (30)
Greece
2
Seli (31)
Greece
3
Pisoderi (32)
Greece
2
Kato Vermio (33)
Greece
1
1
1
1
1
1
1
4
Kiev
Chernobyl
Kauhava
Svalbard
Istanbul
Mt Erciyes
Kangal-Sivas
Lake Yliki (18)
Ukraine
Bellarus
Finland
Norway
Turkey
Turkey
Turkey
Greece
4
Gravia (19)
Greece
2
Karyes (20)
Greece
4
Neo Monastiri (21)
Greece
4
Sykourio (22)
Greece
3
Amyntaio (23)
Greece
3
Paranesti (24)
Greece
4
Dadia (25)
Greece
3
Mikro Dereio (26)
Greece
1
2
Bitola
Kirşehir
FYROM
Turkey
Elmah
Turkey
GR28
GR29.1
GR29.2
GR29.3
GR29.4
GR30.1
GR30.2
GR30.3
GR30.4
GR31.6
GR31.7
GR32.5
GR32.6
GR32.7
GR33.5
GR33.6
UKRAINE
BELLARUS
FINLAND
NORWAY
TURKEY3
TURKEY4
TURKEY5
GR18.1
GR18.2
GR18.3
GR18.4
GR19.1
GR19.2
GR19.3
GR19.4
GR20.1
GR20.2
GR21.1
GR21.2
GR21.3
GR21.4
GR22.1
GR22.2
GR22.3
GR22.4
GR23.1
GR23.2
GR23.3
GR24.1
GR24.2
GR24.3
GR25.1
GR25.2
GR25.3
GR25.4
GR26.1
GR26.2
GR26.3
FYROM
Turkey GW3
Turkey GW4
Turkey GW5
FJ641157
FJ641142
FJ641143
FJ641144
FJ641145
FJ641146
FJ641147
FJ641148
FJ641149
FJ641150
FJ641153
FJ641154
FJ641155
FJ641156
FJ641151
FJ641152
DQ015676∗
U54495∗
AY513819∗
AY513820∗
AY513821∗
AY513822∗
AY513823∗
FJ641117
FJ641118
FJ641119
FJ641120
FJ641129
FJ641130
FJ641131
FJ641132
FJ641137
FJ641138
FJ641121
FJ641122
FJ641123
FJ641124
FJ641110
FJ641111
FJ641112
FJ641113
FJ641114
FJ641115
FJ641116
FJ641133
FJ641134
FJ641136
FJ641125
FJ641126
FJ641127
FJ641128
FJ641139
FJ641140
FJ641141
FJ767744∗
FJ767745∗
FJ767746∗
FJ767747∗
(Continued)
368
E. Thanou et al.
Table I. (Continued).
Species
Site (as in Figure 1)
Country
Code
Ass. Nos
2
Kortuteli
Turkey
Turkey GW6
Turkey GW7
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
2
Radi’ah
Aqrabat
unknown
unknown
Brandenburg
Mazarsay
Balkash
Brandenburg
Trendo
Schötz
Lauwersee
Hjerl Hede
Bialystok
Nuijamaa
Lagny
Mantet, Pyrenees
Crimea
Kavka River, Serov
Sisian
Ninotsminda
Iori River Valley
Syria
Syria
Israel
Israel
Germany
Kyrgyzstan
Kyrgyzstan
Germany
Italy
Switzerland
Netherlands
Denmark
Poland
Finland
France
Spain
Ukraine
Russia
Armenia
Georgia
Georgia
M. anatolicus
1
1
1
1
3
Reine
Ortaköy - Aksaray
Amasya
Boyali Köyü - Amasya
Cihanbeyli
Iran
Turkey
Turkey
Turkey
Turkey
M. irani
3
Balkusan
Turkey
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
Shiraz
Damar
Beniani
Bağdaşan
Trendo
Fusio
Tretie Roháčske
Smutná dolina valley
Steinberg am Rofan
unknown
Viterbo
Torino, Piedmont
Cerano, Piedmont
Fiume, Freddo
Burgos
Melgar de Fernamental
Setùbal
Algarve
Arrós, Vall d’Aran
Riba
Hecho
Trento
Pilatus
Iran
Turkey
Georgia
Turkey
Italy
Switzerland
Slovakia
Slovakia
Austria
Croatia
Italy
Italy
Italy
Italy
Spain
Spain
Portugal
Portugal
Spain
Spain
Spain
Italy
Switzerland
FJ767751∗
FJ767752∗
FJ767743∗
AY513805∗
AY513806∗
AY513807∗
DQ768147∗
AY513808∗
AY513809∗
DQ768131∗
AY220766∗
AY332711∗
AY220778∗
AY220776∗
EU439457∗
AY220770∗
AY220787∗
AY220789∗
AY220762∗
AY220764∗
AY220761∗
AY220760∗
AY513829∗
AY513830∗
AY513831∗
AY513793∗
AY513794∗
AY513795∗
FJ767740∗
FJ767741∗
FJ767742∗
FJ767748∗
FJ767749∗
FJ767750∗
FJ767739∗
DQ841703∗
AY513790∗
AY513791∗
DQ663663∗
AY332713∗
DQ841701∗
AY513838∗
DQ841693∗
EF379100∗
AY513824∗
AY513825∗
AY513826∗
AY513827∗
AY513812∗
AY513813∗
AY513796∗
AY513797∗
AY513799∗
AY513800∗
AY513801∗
DQ663668∗
AY332716∗
M. guentheri “south”
Downloaded by [University of Patras] at 01:52 12 September 2012
M. agrestis
M. juldaschi
M kirgisorum
M. arvalis “arvalis”
M. arvalis “obscurus”
M. socialis
M. dogramacii
M. majori
M. daghestanicus
M. multiplex
M. tatricus
M. bavaricus
M. liechtensteini
M. savii
M. lusitanicus
M. duodecimcostatus
M. gerbei
Chionomys nivalis
No
Vole phylogeography infers multiple Balkan subrefugia
local optima (Huelsenbeck & Bollback 2001). All
trees collected at stationarity were used to estimate
the topology of a 50% majority rule tree and its
posterior nodal probabilities. For the ML analysis
we used the PHYML online web server (Guindon
et al. 2005) under the nearest neighbour interchange (NNI) method, with 1000 bootstrap replicates and the fixed values of gamma distributed shape
parameter (1.04) and proportion of invariable sites
(0.50) that were estimated by MODELTEST.
Results
Downloaded by [University of Patras] at 01:52 12 September 2012
Sequence composition and variation
The length of the mtDNA segment analysed was
710 bp including 286 variable and 249 parsimony
informative characters (40% and 35% of total length,
respectively). Most substitutions were transitions
(72%) and the transition to transversion ratio (R) was
0.9 for the 1st codon, 2.3 for the 2nd codon, 8.8 for
the 3rd codon, while the overall ratio was 2.5. The
majority of polymorphic sites in this segment were at
the 3rd position (27.5%), followed by 1st (5.3%) and
2nd position (0.9%).
Net uncorrected p-distances among populations of
each species, geographically grouped according to
their Greek, European or Asian origin (Table II),
ranged from 0.6% (Greek M. guentheri “west” –
Turkish M. guentheri “west”) to 1.6%, the latter estimated for the pair M. thomasi thomasi M. thomasi atticus. Values twice as high (3.0–3.3%)
were observed between Greek populations of M. subterraneus and their European or Turkish conspecifics.
Phylogenetic results
The Xia’s test for the saturation of phylogenetic
information showed that the substitution saturation
index of our sequences (Iss) is significant lower
than the critical value (Iss<Iss.c, p < 0.001). Thus
our sequences can be used for phylogenetic analyses, although some saturation occurred at the third
codon position of cytb (value of 0.160 combared to
0.148 for all three codons).
The same tree topology (Figure 2), was recovered by the BI and ML analyses regarding the major
clades, with few minor differences in some terminal branches. This tree depicted the phylogenetic
relationships between most Eurasian Microtus species
in agreement with previous studies (Jaarola et al.
2004; Kryštufek et al. 2009b). The five species this
study focuses on are placed in three distinct clades.
All species of the subgenus Terricola formed one
clade further divided in two subclades: one clustered
369
M. subterraneus together with two other species
widely distributed in Asia (M. daghestanicus and
M. majori), and the other the European endemics
of the Mediterranean peninsulas, such as M. thomasi
(including M. thomasi atticus) and M. felteni. The second major clade included the arvalis group species
(M. arvalis “arvalis”, M. arvalis “obscurus” and
M. levis), and the third major clade the socialis
group species (M. socialis, M. anatolicus, M. irani,
M. dogramacii, and the west and south clades of
M. guentheri). Both M. subterraneus and M. levis were
subdivided in three subclades that separated Greek,
other European and Turkish populations, while
within the west clade of M. guentheri, Greek populations were clustered distinctly from their Turkish
conspecifics.
Discussion
The recent application of molecular analyses to several species across Europe has revealed a similar pattern of postglacial colonization from southern refugia. In many cases, hybrid zones have been described
in the Pyrenees and the Alps and DNA sequence
analysis has revealed that these were formed by
secondary contact of distinct subspecific genetic lineages emanating from refugia in Iberia, Italy and the
Balkans. The Balkan expansion came to recolonize
most of Europe, probably due to an early start in the
east after the ice age, and because the Pyrenees and
Alps hindered the spread of Iberian and Italian lineages. Nevertheless, other lineages may be also found
to occur in Greece and Turkey, and possibly further
east in the Caucasus (Hewitt 2001).
In this sense, most small mammals in Greece
are expected to share the same genome with their
northern European conspecifics, as Greek populations probably survived during glacial periods and
expanded to the north during postglacial ones.
Indeed, all Greek populations of Crocidura leucodon
belong to the European mtDNA lineage without
exceptions (Dudey et al. 2007b). Crocidura suaveolens, on the other hand, follows the same pattern,
in general, but also includes a unique mtDNA lineage from Crete Island and another discrete common lineage distributed in the NE Aegean Island
of Lesbos and W Turkey (Dudey et al. 2007a).
Another small mammal, the bank vole, shows a
distinct phylogenetic lineage distributed in the SE
Balkans and W Turkey, while central and northern
European populations descend from refugia situated
outside the Mediterranean Peninsulas (Deffontaine
et al. 2005). These examples provide evidence of
a complex history of mammalian species expansion and colonization in the southern parts of the
[1] M. guentheri “west” – Greek
[2] M. guentheri “west” – Turkish
[3] M. guentheri “south”
[4] M. dogramacii
[5] M. anatolicus
[6] M. irani
[7] M. socialis
[8] M. levis – Greek
[9] M. levis – European
[10] M. levis – Turkish
[11] M. arvalis ‘arvalis’
[12] M. arvalis ‘obscurus’
[13] M. subterraneus – Greek
[14] M. subterraneus – European
[15] M. subterraneus – Turkish
[16] M. felteni – Greek
[17] M. felteni – European
[18] M. thomasi – Greek
[19] M. thomasi – European
[20] M. thomasi atticus
0.6
6.3
3.9
9.3
6.3
6.4
11.4
11.4
11.8
11.0
11.6
13.2
13.6
14.2
13.1
12.6
10.5
10.9
11.1
[1]
5.8
3.7
8.4
6.0
6.0
11.1
11.1
11.6
10.4
11.4
13.5
13.4
13.9
13.3
12.3
10.4
10.9
10.8
[2]
5.1
7.7
7.2
6.8
11.6
11.1
11.6
9.7
10.6
15.1
15.2
14.9
14.1
13.1
11.2
11.7
12.2
[3]
7.8
5.2
5.3
10.2
10.7
10.3
9.3
10.0
13.6
12.6
12.7
11.4
10.4
9.1
9.3
9.8
[4]
5.2
6.0
12.8
13.2
12.8
10.6
11.9
15.3
15.1
15.1
13.6
12.5
11.9
12.5
12.1
[5]
3.8
10.5
10.2
10.3
8.7
8.6
12.3
11.8
12.1
11.3
9.7
9.2
9.8
9.5
[6]
11.4
11.7
11.6
10.1
11.3
13.2
12.5
12.4
11.9
10.7
9.8
10.4
10.1
[7]
1.2
1.3
5.3
4.8
13.5
13.6
13.0
14.1
13.0
11.0
11.4
11.2
[8]
1.5
5.5
4.6
14.3
14.1
13.6
14.5
13.0
11.8
12.2
11.4
[9]
5.8
4.8
13.8
13.0
12.4
13.7
12.3
10.7
11.1
11.0
[10]
2.9
12.9
12.4
11.6
13.1
11.8
10.9
11.2
10.1
[11]
11.9
12.2
11.8
13.1
12.1
11.0
11.4
11.0
[12]
3.0
3.3
13.0
12.7
11.3
12.0
11.1
[13]
2.8
13.4
12.0
10.3
11.0
10.4
[14]
14.5
13.0
11.0
11.6
9.9
[15]
1.5
8.7
9.3
9.5
[16]
7.0
7.6
7.8
[17]
0.2
1.6
[18]
2.0
[20]
Table II. The net uncorrected p-distances (%), calculated between groups of Greek, other European and Asian populations of the five vole species studied. Genetic distances of relative
species belonging to the “socialis” and the “arvalis” group are also presented.
Downloaded by [University of Patras] at 01:52 12 September 2012
370
E. Thanou et al.
371
Downloaded by [University of Patras] at 01:52 12 September 2012
Vole phylogeography infers multiple Balkan subrefugia
Figure 2. BI and ML tree topology and respective statistical support based on the cytb dataset (posterior probabilities >0.95 followed
by ML bootstrap values >75% are shown above branches). Codes for the sequences included in the analysis are listed in Table I. Maps
show the approximate worldwide range of distribution for the five vole species studied (based on Mitchell-Jones et al. 1999; Kryštufek &
Vohralik 2005; Kryštufek et al. 2009b; Wilson & Reeder 2005; IUCN 2007 European Mammal Assessment Website). Darker-shaded areas
correspond to (a) the 52-chromosome karyotype of M. subterraneus and (e) the approximate distribution of the M. guentheri “south” clade.
Balkan Peninsula. Especially Greece, may have been
colonized through W Turkey, N Balkans (connected
either with NW Europe or the Caucasus) or from
subrefugia within the S Balkan area, where isolated
populations survived and remain in recent times as
unique phylogenetic lineages.
In this study, the cytb analysis of several Greek
populations from the five species of Microtus voles, in
comparison to other European and Asiatic conspecific populations has provided a different phylogeographical pattern for each species, which is discussed
separately, under the light of possible subrefugia that
372
E. Thanou et al.
may have played a role in the colonization of the S
Balkan area.
Downloaded by [University of Patras] at 01:52 12 September 2012
M. levis
M. levis is the only representative of the “arvalis”
group in Greece and European Turkey (Kefelioğlu
1995; Mitchell-Jones et al. 1999). The present
species’ distribution in Europe and W Asia is given
in Figure 2d.
Greek M. levis voles are differentiated from
their European and Turkish conspecifics, although
phylogenetic relationships between the three resulting clades are unresolved. Genetic distances between
Greek and either European or Turkish lineages were
almost equal (Table II), and their split seems to have
occurred simultaneously. Interestingly, the genetic
variation within each of these lineages is low, despite
the fact that both European and Turkish ones include
populations from very distant geographic locations.
The other representative of the “arvalis” group,
M. arvalis, is known to have a complex evolutionary history, which results in the current existence
of many distinct phylogenetic lineages corresponding
to multiple European refugia (Tougard et al. 2008).
Our data also imply a respective phylogenetic pattern
for M. levis, but lack of paleontological and molecular data that fully represent the intraspecific variation
of European and Asian populations does not permit
safe conclusions. So far, the distinct lineage represented by Greek voles, may suggest that N Europe
was colonized through different routes, other than
the Balkans or Minor Asia. A Balkan refugium might
have acted as the gene pool for the Greek lineage and
possibly other populations in adjacent regions.
M. guentheri
As with many other species of social voles,
M. guentheri is mainly distributed in a wide area
from SE Balkans eastwards to W Turkey and southwards to parts of the Middle East (Figure 2e). In a
recent study (Kryštufek et al. 2009b), molecular
data revealed two distinct monophyletic lineages that
divided Guenther’s voles in a “west” and a “south”
clade, while the second clade shared a sister-species
relationship with M. dogramacii. Additionally, the
genetic distance between these two clades were high
enough to suggest that they may be independent
species.
In our study, we covered most of the species’
European distribution (Figure 1e), in order to provide more data on the support of the “west” and
“south” clades’ phylogeny. Our molecular results
support the placement of Greek social voles from
this study and one previously published sequence
from FYROM (Former Yugoslavian Republic of
Macedonia) in a new taxon, M. hartingi BarrettHamilton, 1903 (type loc.: Larissa, Greece), in
agreement with Kryštufek et al. (2009b). All Balkan
social voles are probably included in this taxon.
It should be noted that these voles occupy cultivated lowlands and seem to be favored rather than
disturbed by expanding agricultures. The species distribution is continuous, at least in Greece, and their
populations may easily retain gene flow, throughout
the SE Balkans.
Moreover, our data revealed further diversification
between Turkish and Balkan M. hartingi voles,
although the respective clades show moderate statistical support (0.69/65 and 0.89/72, respectively).
Divergence between these groups is low (0.6%,
Table II), presenting the lowest genetic distance
between “Greek” and “non-Greek” populations
of all studied species. Thus it is probable that
present diversification within M. hartingi resulted
from recent isolation events and the disconnection
between Greece and Anatolia, due to sea level rise
after glacial retreat (Perissoratis & Conispoliatis
2003), rather than the existence of a Balkan
refugium.
M. subterraneus
This species is distributed in a wide area, which
includes most of eastern and central Europe (Çolak
et al. 1998; Mitchell-Jones et al. 1999) (Figure 2a).
Many previous studies have reported differences in
size, allozymic variation and morphological differentiation, when M. subterraneus from the Balkans
and other European populations were compared
(Neithammer 1982b; Kryštufek et al. 1994 and
references therein; Fraguedakis-Tsolis et al. 2009).
Significant intraspecific cytb variation has also been
reported for this species (Jaarola et al. 2004),
although few specimens were tested. Nevertheless,
its phylogeographic history has only been discussed
under the light of karyological studies (Macholán
et al. 2001; Zima 2004), since central and eastern
European populations display a chromosomal number of 2n = 52, in contrast to the typical karyotype
of 2n = 54 found mainly in NW Europe and Turkey
(Figure 2a).
Although Greek and European voles studied here
belong to the 52-chromosome race (Macholan et al.
2001; Mitsainas et al. 2010), our cytb analysis placed
them in two separated and well supported mtDNA
lineages (Figure 2). According to our data, M. subterraneus populations presently distributed in Greece
are highly differentiated from those living in central
Downloaded by [University of Patras] at 01:52 12 September 2012
Vole phylogeography infers multiple Balkan subrefugia
Europe, but still more closely related to them than to
Turkish ones. The genetic distances among Greek,
European and Turkish pine voles are considerably
high (2.8–3.3%, Table II). Genetic distances of such
high values were calculated between pairs of closely
related but distinct species (M. irani – M. socialis,
M. dogramacii – M. guentheri “south”, Table II).
Sufficient representatives from northern European
and Turkish populations should be tested prior to any
suggestions regarding phylogenetic and taxonomic
relationships within M. subterraneus. Yet, our molecular analysis suggest the existence of a Balkan
refugium that sheltered a Balkan lineage currently
distributed in Greece, while Central Europe was
probably not colonized by the dispersal of this lineage.
On the other hand, Greek populations showed a
significant intraspecific variation that distinguished
three subgroups (Figure 2). The first includes Mts
Vrontous and Silo (Figure 1: locations 36, 37) which
correspond to the western and eastern parts of
the Rodopi mountain range along the Greece –
Bulgaria borders. The second include Mts Vorras
and Pinovo (Figure 1: locations 34, 35), the highest
ridges of Tzena mountain range that descends from
NW to SE between FYROM and Greece. The third
includes the remaining populations (Figure 1: locations 31, 32, 33) distributed along the eastern slopes
of Mt Pindos. This pattern implies the existence of
several subrefugia within the southern Balkans and
possibly in the numerous mountain ranges of N
Greece. All specimens of this study were found in
high elevations (over 1500 m). Although common
pine voles may be found in lower altitudes in northern Europe, in the Balkans they occur in mountainous areas, mainly in coniferous and broadleaf forests.
Their habitats are restricted and have suffered severe
fragmentation, which may promote isolation, at least
in the Southern Balkan Peninsula.
M. thomasi and M. felteni
These two species, which retain primitive dental
characters, are considered archaic species of the
European Microtus species group (Brunet-Lecompte
& Chaline 1992). A previous study that included
only one M. felteni specimen and four M. thomasi
ones (Jaarola et al. 2004), suggested a sister-species
relationship between them, which might indicate
common phylogenetic histories, as probably reflected
in their present distribution (Figures 2b, 2c). In our
study we evaluate this conclusion analyzing additional specimens representing more populations from
Greece, covering a major part of both species extant
distribution.
373
Our results demonstrate a polytomy between
clades corresponding to the European endemics,
namely the “multiplex” group, the Italian Peninsula
endemic (M. savii), the Iberian Peninsula endemics
(M. lusitanicus/M. duodecicomstatus) and the Balkan
Peninsula endemics (M. thomasi and M. felteni)
(Figure 2). Our data do not seem to support a sisterspecies relationship between M. thomasi and M. felteni (Figure 2, Table II). On the contrary, the genetic
distance between them is the greatest one for any
pair of the Mediterranean endemic species (7.9%),
while the lowest values for M. felteni and M. thomasi
are the ones exhibited between each of them and
the M. multiplex+M. bavaricus+M. liechtensteini clade
(6.5% and 4.2% respectively). Although this study
does not aim to conclude on the phylogeny between
Microtus species, our results provide evidence against
the sister-species relationship of M. thomasi and
M. felteni and indications of a possibly more complex role of the Balkan refugium on the speciation of
Balkan endemics.
Tree polytomies may result from a lack of information and can be resolved into sequential bifurcations
given additional data (“soft polytomies”) or they may
depict multiple simultaneous speciation events that
generated three or more lineages (“hard” polytomies)
(Maddison 1989; Walsh et al. 1999). Regarding the
polytomy discussed here, also observed in Jaarola
et al. (2004), the scenario of rapid radiation and
simultaneous diversification of many Microtus lineages resulting to the speciation of many European
endemics is very plausible. The two Balkan endemics
may not share a common ancestor but a common Balkan refugium. Both species survived there
through the ice ages and subsequently expanded.
This expansion was blocked to the north by the
Alps and the Dinaric Mountains, where the species
presently forming the “multiplex” group were probably isolated (Martínková et al. 2007). At the intraspecific level, M. felteni showed significant genetic differentiation (1.5%) considering its restricted geographic
distribution (Figures 1c, 2b) and in comparison to
M. thomasi (excluding M. thomasi atticus) which covers a much wider area (Figures 1b, 2c) and exhibits
much lower intraspecific differentiation (see also
Tryfonopoulos et al. 2008). High genetic differentiation implies limited gene-flow between M. felteni
populations and could be attributed to a patched
distribution within the species’ range.
Discontinued distribution of this Balkan endemic
species is also corroborated by our field observations;
thorough efforts and repeated visits during the past
decades were unsuccessful in capturing M. felteni
specimens, from most of its reported expansion
area in Greece (Ondrias 1966; Niethammer 1982a).
Downloaded by [University of Patras] at 01:52 12 September 2012
374
E. Thanou et al.
Felten’s voles were only found in a much smaller
area close to the Greek borders with Albania and
FYROM, in agreement with other recent studies
(Andera 1991; Mitsainas et al. 2010). All sampling
locations (Figure 1c) were at high altitudes, on the
east side of Mt. Pindos. Moreover, it seems that
these voles form small, discontinuous populations,
usually in sympatry with much larger populations
of M. thomasi, in Greece (authors’ personal observations) and probably in other Balkan areas (Bego
et al. 2008). Unexpectedly, M. felteni has recently
been found (Mitsainas et al. 2010; present study)
in a location outside previously known Greek distribution (Figure 1, location 31) suggesting the need
to update the biogeographical information regarding this species. Although preliminary and including
only few M. felteni populations, the present study is
the first molecular approach for this species’ phylogeny and biogeography. Our results demonstrate an
intraspecific genetic diversification probably due to
parapatric isolation in subrefugia within the mountainous S Balkan area
On the other hand, M. thomasi populations are
divided into two well supported groups, as all
populations from Evvoia and SE Sterea Ellada
were highly differentiated from Thomas’ voles living in the remaining parts of Greece and other
Balkan areas (Figures 1, 2c). These Greek populations retain a distinct phylogenetic profile, as they
form a monophyletic clade when several Microtus
species are included in the phylogenetic analyses.
Additionally, they exhibit a sister-clade relationship
with M. thomasi. Thus, our present data further support the validity of the atticus clade as a distinct taxon
and confirm its attribution to the subspecific level, as
was discussed in Tryfonopoulos et al. (2008) or even
specific level as was proposed by Rovatsos & Giagia
Athanasopoulou (2011). Moreover, new data suggest
that the distributional range of the new taxon in the
SE Sterea Ellada is wider than previously supposed,
although populations are small and discontinuous.
Despite thorough field work, we were unable to
identify a possible contact zone between the two subspecies. It is possible that the two lineages are not in
contact thus maintaining the isolation of the atticus
populations, possibly due to either the expansion and
competitive behavior of M. guentheri (Tryfonopoulos
et al. 2008) and/or to the geographical features of the
area. During the Late Pleistocene, flooded masses
of water, which resulted from melting ice in the
north, isolated SE Sterea Ellada from neighboring
areas (Palyvos 2001). Presently, major rivers with a
west-east direction, flowing from east Mt. Pindos,
still divide SE Sterea Ellada from adjoining
regions.
Conclusions
Our results have shown that all species of Greek
voles are genetically differentiated compared to their
European and/or Asiatic conspecifics, in consistency
with previous morphological studies (FraguedakisTsolis et al. 2009). The factors promoting vole
differentiation in this region probably acted during
the Middle and Late Pleistocene glacial periods, due
to the paleogeographical and paleoclimatical characteristics prevailing at that time. Yet, the pattern of
differentiation is characteristic for each of the studied
species. In the case of M. hartingi, the W Turkey and
the Balkan clade were differentiated due to parapatric
isolation when the Bosporus Trait was submerged
under sea-level disconnecting these two areas. On the
other hand, the genetic lineages found within M. levis
may suggest the existence of a Balkan refugium which
provided the genetic stock for a distinct Balkan lineage but did not expand further north to recolonize
Europe. Finally, the phylogeographic history of the
two Balkan endemics, M. thomasi and M. felteni as
well as M. subterraneus suggests that their respective intraspecific variation can be assigned to the
existence of multiple Balkan subrefugia.
Although preliminary, this study provides evidence
which imply that along with the Italian (Castiglia
et al. 2008), and the Iberian (Centeno-Cuadros
et al. 2009), the southern Balkan Peninsula not
only acted as a refugium, but it most probably
included multiple subrefugia, probably situated at
the southern mountainous parts of the Peninsula.
Moreover, it raises several interesting issues regarding the phylogeography of Eurasian vole species
(M. felteni, M. subterraneus and M. levis), which may
hold cryptic variation.
Acknowledgements
The authors wish to thank Dr. C. Stamatopoulos,
Mrs. Ch. Valassaki, Mr. A. Goulios and Mrs. Irena
Gjoni for their assistance in field work and James
Howlett for linguistic advice. This study was partly
funded by the European Social Fund (ESF),
Operational Program for Educational and Vocational
Training II (EPEAEK II, Greece), and particularly
the Program PYTHAGORAS I.
References
Andera M. 1991. First records of Microtus (Pitymys) felteni in
Albania. Mammalian Biology 56:380–387.
Bego F, Kryštufek B, Paspali G, Rogoz E. 2008. Small terrestrial
mammals of Albania: annotated list and distribution. Italian
Journal of Mammalogy 19:3–21.
Downloaded by [University of Patras] at 01:52 12 September 2012
Vole phylogeography infers multiple Balkan subrefugia
Bilton DT, Mirol PM, Mascheretti S, Fredga K, Zima J, Searle JB.
1998. Mediterranean Europe as an area of endemism for small
mammals rather than a source for northwards postglacial colonization. Proceedings of the Royal Society of London Series
B-Biological Sciences 265:1219–1226.
Brunet-Lecompte P, Chaline J. 1992. Morphological convergence
versus biochemical divergence in the holarctic ground voles:
Terricola and Pitymys (Arvicolidae, Rodentia). Neues Jahrbuch
fur Geologie und Palaontologie, Munchen 12:721–734.
Bužan EV, Förster DW, Searle JB, Kryštufek B. 2010. A new
cytochrome b phylogroup of the common vole (Microtus
arvalis) endemic to the Balkans and its implications for the
evolutionary history of the species. Biological Journal of the
Linnean Society 100:788–796.
Castiglia R, Annesia F, Aloise G, Amori G. 2008. Systematics
of the Microtus savii complex (Rodentia, Cricetidae) via
mitochondrial DNA analyses: Paraphyly and pattern of
sex chromosome evolution. Molecular Phylogenetics and
Evolution 46:1157–1164.
Centeno-Cuadros A, Delibes M, Godoy JA. 2009.
Phylogeography of Southern Water Vole (Arvicola sapidus):
evidence for refugia within the Iberian glacial refugium?
Molecular Ecology 18:3652–3667.
Çolak E, Yigit N, Sözen M, Özkurt Ş. 1998. A study on taxonomic
status of Microtus subterraneus de Selys Longchamps, 1836 and
Microtus majori Thomas, 1906 (Mammalia, Rodentia) in
Turkey. Turkish Journal of Zoology 22:119–129.
Conroy CJ, Cook JA. 2000. Molecular systematics of a holarctic rodent (Microtus Muridae). Journal of Mammalogy
81:344–359.
Deffontaine V, Libois R, Kotlík P, Sommer R, Nieberding
C, Paradis E, Searle JB, Michaux JR. 2005. Beyond the
Mediterranean peninsulas: evidence of central European
glacial refugia for a temperate forest mammal species,
the bank vole (Clethrionomys glareolus). Molecular Ecology
14:1727–1739.
Dubey S, Cosson J-F, Magnanou E, Vohralík V, Benda P,
Frynta D, Hutterer R, Vogel V, Vogel P. 2007a. Mediterranean
populations of the lesser white-toothed shrew (Crocidura suaveolens group): an unexpected puzzle of Pleistocene survivors and
prehistoric introductions. Molecular Ecology 16:3438–3452.
Dubey S, Cosson J-F, Vohralík V, Kryštufek B, Diker E,
Vogel P. 2007b. Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the
case of the bicoloured shrew (Crocidura leucodon, Soricidae).
Journal of Evolutionary Biology 20:1799–1808.
Fraguedakis-Tsolis SE, Chondropoulos BP, Stamatopoulos CV,
Giokas S. 2009. Morphological variation of the five vole species
of the genus Microtus (Mammalia, Rodentia, Arvicolinae)
occurring in Greece. Acta Zoologica 90:254–264.
Fraguedakis-Tsolis SE, Ondrias JC. 1985. Geographic variation of the ground squirrel Citellus citellus (Mammalia,
Rodentia) in Greece with the description of a new subspecies.
Säugetierkundliche Mitteilungen 32:185–198.
Graf JD. 1982. Biochemical genetics, zoogeography and taxonomy
in the family Arvicolidae (Mammalia, Rodentia). Revue Suisse
De Zoologie 89:749–787.
Guindon S, Lethiec F, Duroux P, Gascuel O. 2005. PHYML
Online: a web server for fast maximum likelihood-based
phylogenetic inference. Nucleic Acids Research (Web Server
issue):W557–W559.
Hall TA. 1999. BioEdit: a user-friendly biological sequence
alignment program for Windows 95/98/NT. Nucleic Acids
Symposium Series 41:85–98.
375
Hellborg L, Gunduz I, Jaarola M. 2005. Analysis of sexlinked sequences supports a new mammal species in Europe.
Molecular Ecology 14:2025–2031.
Hewitt GM. 2001. Speciation, hybrid zones and phylogeography
– or seeing genes in space and time. Molecular Ecology
10:537–549.
Huelsenbeck JP, Bollback JP. 2001. Empirical and hierarchical
Bayesian estimation of ancestral states. Systematic Biology
50:351–366.
IUCN 2007. Microtus felteni. Microtus guentheri. Microtus
levis. Microtus subterraneus. Microtus thomasi. In: IUCN
2007. European Mammal Assessment. hppt://ec.europa.eu/
environment/nature/conservation/species/ema/. Accessed Nov
2010 20.
Jaarola M, Martínkova N, Gündüz I, Brunhoff C, Zima J,
Nadachowski A, Amori G, Bulatova NS, Chondropoulos B,
Fraguedakis-Tsolis S, González-Esteban J, López-Fuster MJ,
Kandaurov AS, Kefelioğlu H, Mathias MD, Villate I, Searle
JB. 2004. Molecular phylogeny of the speciose vole genus
Microtus (Arvicolinae, Rodentia) inferred from mitochondrial
DNA sequences. Molecular Phylogenetics and Evolution
33:647–663.
Kefelioğlu H. 1995. The taxonomy of the genus Microtus
(Mammalia, Rodentia) and its distribution in Turkey. Turkish
Journal of Zoology 19:35–63.
Kryštufek B, Bryja J, Bužan EV. 2009a. Mitochondrial
phylogeography of the European ground squirrel, Spermophilus
citellus, yields evidence on refugia for steppic taxa in the
southern Balkans. Heredity 103:129–135.
Kryštufek B, Bužan EV, Vohralík V, Zareie R, Özkan B. 2009b.
Mitochondrial cytochrome b sequence yields new insight into
the speciation of social voles in south-west Asia. Biological
Journal of the Linnean Society 98:121–128.
Kryštufek B, Filipucci MG, Macholán M, Zima J, Vujošević
M, Simson S. 1994. Does Microtus majori occur in
Europe? Zeitschrift Fur Saugetierkunde-International Journal
of Mammalian Biology 559:349–357.
Kryštufek B, Vohralík V, editors. 2005. Mammals of Turkey and
Cyprus Volume 2. Rodentia I: Sciuridae, Dipodidae, Gliridae,
Arvicolidae. Koper: University of Primorska. pp. 149–209.
Kumar S, Tamura K, Nei M. 2004. MEGA3: Integrated Software
for Molecular Evolutionary Genetics Analysis and Sequence
Alignment. Briefings in Bioinformatic 5:150–163.
Macholàn M, Filippucci MG, Zima J. 2001. Genetic variation and
zoogeography of pine voles of the Microtus subterraneus/majori
group in Europe and Asia Minor. Journal of the Zoological
Society of London 255:31–42.
Maddison W. 1989 Reconstructing character evolution on polytomous cladograms. Cladistics 5:365–377.
Martínková N, Zima J, Jaarola M, Macholán M, Spitzenberger F.
2007. The origin and phylogenetic relationships of Microtus
bavaricus based on karyotype and mitochondrial DNA
sequences. Folia Zoologica 56:39–49.
Michaux JR, Libois M, Filippucci M-G. 2005. So close and so
different: comparative phylogeography of two small mammal
species, the Yellow-necked fieldmouse (Apodemus flavicollis)
and the Woodmouse (Apodemus sylvaticus) in the Western
Palearctic region. Heredity 94:52–63.
Mitchell-Jones AJ, Amori G, Bogdanowicz W, Krystûfek B,
Reijnders PJH, Spitzenberger F, Stubbe M, Thissen JB,
Vohralík V, Zima J. 1999. The Atlas of European mammals.
T & AD Poyser Ltd. London: Academic Press.
Mitsainas GP, Rovatsos MT, Giagia-Athanasopoulou EB.
2010. Heterochromatin study and geographical distribution
Downloaded by [University of Patras] at 01:52 12 September 2012
376
E. Thanou et al.
of Microtus species (Rodentia, Arvicolinae) from Greece.
Mammalian Biology 75:261–269.
Mitsainas GP, Tryfonopoulos GA, Thanou EG, Bisa R,
Fraguedakis-Tsolis SE, Chondropoulos BP. 2009. New data
on the distribution of Mus spicilegus Petenyi, 1882 (Rodentia,
Muridae) and a distinct mtDNA lineage in the southern
Balkans. Mammalian Biology 74:351–360.
Musser GG, Carleton MD. 2005. Superfamily Muroidea. In:
Wilson D, Reeder DM, editors. Mammal Species of the World:
A Taxonomic and Geographic Reference. Baltimore, MD:
Johns Hopkins University Press. pp. 510–756.
Niethammer J. 1982a. Microtus felteni (Malec und Storch, 1963).
In: Niethammer J, Krapp F, editors. Handbuch der Säugetiere
Europas. Wiesbaden: Akademische Verlagsgesellschaft. pp.
438–441.
Niethammer J. 1982b. Microtus subterraneus (de SélysLongchamps, 1836)-Kurzohrmaus. In: Niethammer J, Krapp
F, editors, Handbuch der Säugetiere Europas. Wiesbaden:
Akademische Verlagsgesellschaft. pp. 397–418.
Ondrias JC. 1966. The taxonomy and geographic distribution
of the rodents of Greece. Säugetierkundliche Mitteilungen
14:1–136.
Palyvos N. 2001. Geomorphological study of the broader Atalanti
area (Fthiotis, Central Greece). Ph.D. Thesis, University of
Athens, Department of Geology, 233 pp. [in Greek].
Perissoratis C, Conispoliatis N. 2003. The impacts of sea-level
changes during latest Pleistocene and Holocene times on
the morphology of the Ionian and Aegean seas (SE Alpine
Europe). Marine Geology 196:145–156.
Piertney SB, Stewart WA, Lambin X, Telfer S, Aars J, Dallas
JF. 2005. Phylogeographic structure and postglacial evolutionary history of water voles (Arvicola terrestris) in the United
Kingdom. Molecular Ecology 14:1435–1444.
Posada D, Crandall KA. 1998. MODELTEST: testing the model
of DNA substitution. Bioinformatics 14:817–818.
Ronquist F, Huelsenbeck JP. 2003. MrBayes 3: Bayesian
phylogenetic inference under mixed models. Bioinformatics
19:1572–1574
Rovatsos MTh, Giagia-Athanasopoulou EB. 2011. Taxonomical
status and phylogenetic relations between the “thomasi”
and “atticus” chromosomal races of the underground vole
Microtus thomasi (Rodentia, Arvicolinae). Mammalian Biology
(in press) doi:10.1016/j.mambio.2011.09.001.
Thompson JD, Higgins DG, Gibson TJ. 1994. CLUSTAL-W
- improving the sensitivity of progressive multiple sequence
alignment through sequence weighting, position-specific gap
penalties and weight matrix choice. Nucleic Acids Research
22:4673–4680.
Tougard C, Renvoisé E, Petitjean A, Quéré J-P. 2008. New insight
into the colonization processes of common voles: inferences
from molecular and fossil evidence. PLoS ONE 3(10):e3532.
doi:10.1371/journal.pone.0003532.
Tryfonopoulos G, Thanou E, Chondropoulos B, FraguedakisTsolis S. 2008. MtDNA analysis reveals the ongoing
speciation on Greek populations of Microtus (Terricola) thomasi
(Arvicolidae, Rodentia). Biological Journal of the Linnean
Society 95:117–130.
Vega R, Amori G, Aloise G, Cellini S, Loy A, Searle JB.
2010. Genetic and morphological variation in a Mediterranean
glacial refugium: evidence from Italian pygmy shrews, Sorex
minutus (Mammalia: Soricomorpha). Biological Journal of the
Linnean Society 100:774–787.
Walsh HE, Kidd MG, Moum T, Friesen VL. 1999. Polytomies
and the power of phylogenetic inference. Evolution
53:932–937.
Wilson DE, Reeder DM. 2005. Mammal Species of the World.
Baltimore, MD: John Hopkins University Press.
Xia X. 2000. Data analysis in molecular biology and evolution.
Boston, MA: Kluwer Academic Publishers.
Xia X, Xie Z. 2001. DAMBE: Software package for data analysis in molecular biology and evolution. Journal of Heredity
92:371–373.
Zima J. 2004. Karyotypic variation in mammals of the Balkan
Peninsula. In: Griffiths HI, Kryštufek B, Reed JM, editors,
Balkan Biodiversity: Patterns and Process in the European
Hotspot. Dordrecht: Kluwer Academik. pp. 109–133.