Sequence analysis of the ribosomal DNA internal transcribed spacer 2 from populations of Anopheles nuneztovari (Diptera: Culicidae)

Sequence variation of the ribosomal DNA internal transcribed spacer 2 (ITS2) was examined for populations of the malaria vector Anopheles nuneztovari collected in Colombia, Venezuela, Bolivia, Suriname, and Brazil. Mosquitoes from Colombia and Venezuela had identical ITS2 sequences and were distinguished from sequences in other populations by three insertion/deletion events (indels) and by one transversion. The length of the ITS2 was 363-369 bp, and it had a G+C content of 55.3%- 55.7%. Variation in the length of the ITS2 between and within populations was due to indels in simple repeats. ITS2 consensus sequences were similar or identical for samples from the following three groups: (1) Colombia, Bolivia, and Venezuela; (2) Suriname and northern Brazil; and (3) eastern and central Brazil. The presence of two different consensus sequences from a single location near Manaus, Brazil, suggests that populations from eastern Brazil and those from Suriname converge in this region of the Amazon Basin. These data show that putative cryptic species of An. nuneztovari are distinguished by very minor differences in DNA sequence of the ITS2 region.      

Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae)

Extensive population structuring is known to occur in Anopheles darlingi , the primary malaria vector of the Neotropics. We analysed the phylogeographic structure of the species using the mitochondrial cytochrome oxidase I marker. Diversity is divided into six main population groups in South America: Colombia, central Amazonia, southern Brazil, south-eastern Brazil, and two groups in north-east Brazil. The ancestral distribution of the taxon is hypothesized to be central Amazonia, and there is evidence of expansion from this region during the late Pleistocene. The expansion was not a homogeneous front, however, with at least four subgroups being formed due to geographic barriers. As the species spread, populations became isolated from each other by the Amazon River and the coastal mountain ranges of south-eastern Brazil and the Andes. Analyses incorporating distances around these barriers suggest that the entire South American range of An. darlingi is at mutation–dispersal–drift equilibrium. Because the species is distributed throughout such a broad area, the limited dispersal across some landscape types promotes differentiation between otherwise proximate populations. Moreover, samples from the An. darlingi holotype location in Rio de Janeiro State are substantially derived from all other populations, implying that there may be additional genetic differences of epidemiological relevance. The results obtained contribute to our understanding of gene flow in this species and allow the formulation of human mosquito health protocols in light of the potential population differences in vector capacity or tolerance to control strategies.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 854–866.     

Phylogeny and biogeography of the alpine newt Mesotriton alpestris (Salamandridae, Caudata), inferred from mtDNA sequences

In this paper, we performed phylogenetic analyses of Mesotriton alpestris populations from the entire range of species distribution, using fragments of two mtDNA genes, cytochrome b (309bp) and 16S rRNA ( approximately 500bp). Sequence diversity patterns and phylogenetic analyses reveal the existence of a relict lineage (Clade A) of late Miocene origin, comprising populations from south-eastern Serbia. This lineage is proposed to be ancestor to a western and an eastern lineage, which diverged during the middle Pliocene. The western lineage is further divided in two clades (Clades B, C) of middle Pliocene origin that represent populations from Italy (B) and populations from central Europe and Iberia (C). Further subdivision, dated back to the middle-late Pliocene, was found within the eastern lineage, representing southern (Clade D) and central-northern (Clade E) Balkan populations, respectively. Extensive sequence divergence, implying greater isolation in multiple refugia, is found within eastern clades, while the western clades seem to have been involved in the colonization of central, western and north-eastern Europe from a hypothetical refugium in central Europe. The extent of divergence does not support the current taxonomy indicating cryptic speciation in the Balkans, while paedomorphic lineages were found to have been evolved during early-middle Pleistocene probably as a response to the ongoing dramatic climatic oscillations.     

Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species

While Quaternary climatic changes are considered by some to have been a major factor promoting speciation within the neotropics, others suggest that much of the neotropical species diversity originated before the Pleistocene. Using mitochondrial and nuclear sequence data, we evaluate the relative importance of Pleistocene and pre-Pleistocene events within the evolutionary history of the Mexican iguana Ctenosaura pectinata , and related species. Results support the existence of cryptic lineages with strong mitochondrial divergence (> 4%) among them. Some of these lineages form zones of secondary contact, with one of them hybridizing with C. hemilopha . Evolutionary network analyses reveal the oldest populations of C. pectinata to be those of the northern and southern Mexican coastal regions. Inland and mid-latitudinal coastal populations are younger in age as a consequence of a history of local extinction within these regions followed by re-colonization. Estimated divergence times suggest that C. pectinata originated during the Pliocene, whereas geographically distinct mitochondrial DNA lineages first started to diverge during the Pliocene, with subsequent divergence continuing through the Pleistocene. Our results highlight the influence of both Pliocene and Pleistocene events in shaping the geographical distribution of genetic variation within neotropical lowland organisms. Areas of high genetic diversity in southern Mexico were detected, this finding plus the high levels of genetic diversity within C. pectinata , have implications for the conservation of this threatened species.     

Nuclear and mitochondrial phylogeography of the Atlantic forest endemic Xiphorhynchus fuscus (Aves: Dendrocolaptidae): biogeography and systematics implications

We studied the intraspecific evolutionary history of the South American Atlantic forest endemic Xiphorhynchusfuscus (Aves: Dendrocolaptidae) to address questions such as: Was the diversification of this bird's populations associated to areas of avian endemism? Which models of speciation (i.e., refuges, river as barriers or geotectonism) explain the diversification within X. fuscus? Does the genetic data support subspecies as independent evolutionary units (species)? We used mitochondrial (n=34) and nuclear (n=68) DNA sequences of X. fuscus to study temporal and spatial relationships within and between populations. We described four main monophyletic lineages that diverged during the Pleistocene. The subspecies taxonomy did not match all the evolutionary lineages; subspecies atlanticus was the only one that represented a monophyletic and isolated lineage. The distribution of these lineages coincided with some areas of endemism for passerines, suggesting that those areas could be regions of biotic differentiation. The ancestor of X. fuscus diverged approximately 3 million years ago from Amazonian taxa and the phylogeographic pattern suggested that X. fuscus radiated from northeastern Brazil. Neither the riverine nor the geotectonic vicariance models are supported as the primary cause for diversification of geographic lineages, but rainforest contractions and expansions (ecological vicariance) can explain most of the spatial divergence observed in this species. Finally, analyses of gene flow and divergence time estimates suggest that the endangered subspecies atlanticus (from northeastern Brazil) can be considered a full species under the general lineage species concept.     

The importance of Anopheles albitarsis E and An. darlingi in human malaria transmission in Boa Vista, state of Roraima, Brazil

In several districts of Boa Vista, state of Roraima, Brazil we found Anopheles (Nyssorhynchus) albitarsis E to be the primary vector of human malaria parasites, and during 2001-2002 it was significantly more abundant than An. darlingi (p < 0.001). Other species sampled were An. (Nys.) braziliensis, An. (Ano.) peryassui, An. (Nys.) nuneztovari, An. (Nys.) oswaldoi s.l., and An. (Nys.) triannulatus. As determined by the ELISA technique An. darlingi had a higher overall infection rate (2.1%) compared with An. albitarsis E (1.2%). However a marginally higher proportion of An. albitarsis E was infected with Plasmodium vivax compared with An. darlingi, and the An. albitarsis E biting index was also much higher These results suggest the importance of An. albitarsis E in malaria transmission in a savannah ecoregion of northern Amazonian Brazil, and reconfirm the importance of An. darlingi even if at lower abundance.     

Can Mosquito Magnet® substitute for human-landing catches to sample anopheline populations?

The efficiency of the Mosquito Magnet Liberty Plus? (MMLP) trap was evaluated in comparison to human-landing catches (HLCs) to sample anopheline populations in Jabillal, state of Bolivar, southern Venezuela. The village comprised 37 houses and a population of 101; malaria in this village is primarily due to Plasmodium vivax and the Annual Parasite Index is 316.8 per 1,000 population. A longitudinal study was conducted between June 2008-January 2009 for three nights per month every two months between 17:30 pm-21:30 pm, a time when biting mosquitoes are most active. Anopheles darlingi and Anopheles nuneztovari were the most common species collected by both methods, whereas Anopheles marajoara was more abundant according to the HLC method. The MMLP trap was more efficient for collecting An. nuneztovari [63%, confidence interval (CI): 2.53] than for collecting An. darlingi (31%, CI: 1.57). There were significant correlations (p < 0.01) between the two methods for An. darlingi [Pearson correlation (R2) = 0.65] and An. nuneztovari (R2 = 0.48). These preliminary results are encouraging for further investigations of the use of the MMLP trap for monitoring anopheline populations in remote malaria-endemic areas in the Amazon Basin.     

Molecular comparison of topotypic specimens confirms Anopheles (Nyssorhynchus) dunhami Causey (Diptera: Culicidae) in the Colombian Amazon

The presence of Anopheles (Nyssorhynchus) dunhami Causey in Colombia (Department of Amazonas) is confirmed for the first time through direct comparison of mtDNA cytochrome c oxidase I (COI) barcodes and nuclear rDNA second internal transcribed spacer (ITS2) sequences with topotypic specimens of An. dunhami from Tefé, Brazil. An. dunhami was identified through retrospective correlation of DNA sequences following misidentification as Anopheles nuneztovari s.l. using available morphological keys for Colombian mosquitoes. That An. dunhami occurs in Colombia and also possibly throughout the Amazon Basin, is of importance to vector control programs, as this non-vector species is morphologically similar to known malaria vectors including An. nuneztovari, Anopheles oswaldoi and Anopheles trinkae. Species identification of An. dunhami and differentiation from these closely related species are highly robust using either DNA ITS2 sequences or COI DNA barcode. DNA methods are advocated for future differentiation of these often sympatric taxa in South America.     

Biogeography and evolution of Amazonian triatomines (Heteroptera: Reduviidae): implications for Chagas disease surveillance in humid forest ecoregions

An ecological-evolutionary classification of Amazonian triatomines is proposed based on a revision of their main contemporary biogeographical patterns. Truly Amazonian triatomines include the Rhodniini, the Cavernicolini, and perhaps Eratyrus and some Bolboderini. The tribe Rhodniini comprises two major lineages (pictipes and robustus). The former gave rise to trans-Andean (pallescens) and Amazonian (pictipes) species groups, while the latter diversified within Amazonia (robustus group) and radiated to neighbouring ecoregions (Orinoco, Cerrado-Caatinga-Chaco, and Atlantic Forest). Three widely distributed Panstrongylus species probably occupied Amazonia secondarily, while a few Triatoma species include Amazonian populations that occur only in the fringes of the region. T. maculata probably represents a vicariant subset isolated from its parental lineage in the Caatinga-Cerrado system when moist forests closed a dry trans-Amazonian corridor. These diverse Amazonian triatomines display different degrees of synanthropism, defining a behavioural gradient from household invasion by adult triatomines to the stable colonisation of artificial structures. Anthropogenic ecological disturbance (driven by deforestation) is probably crucial in the onset of the process, but the fact that only a small fraction of species effectively colonises artificial environments suggests a role for evolution at the end of the gradient. Domestic infestation foci are restricted to drier subregions within Amazonia; thus, populations adapted to extremely humid rainforest microclimates may have limited chances of successfully colonising the slightly drier artificial microenvironments. These observations suggest several research avenues, from the use of climate data to map risk areas to the assessment of the synanthropic potential of individual vector species.     

High levels of population structure caused by habitat islands in the malarial vector Anopheles scanloni

The genetic structure of four populations of the malarial vector Anopheles scanloni in Thailand was studied using mitochondrial DNA sequences. Four highly divergent lineages were observed, all with signals of population expansion. Since An. scanloni is restricted to 'islands' of limestone karst habitat, this suggests there is a metapopulation-type dynamic in this species, with restricted gene flow, extinctions and drift all contributing to lineage divergence. Historical environmental change and marine transgressions may also have contributed to population extinction, expansion and divergence. Although there is some current gene flow inferred between nearby populations, it is extremely restricted between the northern and southern populations, which also differed by one fixed polymorphism at the ITS2 rDNA locus. Crossing experiments showed no post-mating barriers existing between the north and the south, but the lack of gene flow between these populations could ultimately result in speciation and has implications for malaria control strategies.     

Molecular phylogeny and plumage signal evolution in a trans Andean and circum Amazonian avian species complex

Species with fragmented distributions are particularly useful models for investigating processes underlying biological diversification in the Neotropics. The Phaeothlypis wood-warbler complex (Aves: Parulidae) is comprised of six disjunct or parapatric populations. The geographic distribution of these six populations mirrors the classic map of Neotropical areas of endemism that were originally proposed as putative Pleistocene forest refugia, but the magnitude of mitochondrial DNA divergence between these populations suggests that they are each substantially older, with origins in the late Pliocene. Phylogenetic reconstructions based on long mtDNA coding sequences show that the Guiana Shield and Atlantic Forest populations are sister lineages, and group this combined lineage and the remaining four population-specific lineages in a five-way hard polytomy. MtDNA-based phylogenetic reconstructions provide no evidence that the three populations with conspicuous yellow rump and tail feathers currently grouped as the Buff-rumped Warbler (P. fulvicauda) form a monophyletic group. Furthermore, there is a broad discordance between mtDNA and plumage along a transect just east of the Andes, where the contact zone between highly divergent mtDNA clades is more than 1000 km north of the phenotypic hybrid zone between the bright and dark plumage forms. This discordance between mtDNA genotype and plumage phenotype is similar to patterns seen on a finer geographic scale in other avian hybrid zones and may result from asymmetric introgression of the bright plumage trait.     

Incrimination of Anopheles (Anopheles) intermedius Peryassú, An. (Nyssorhynchus) nuneztovari Gabaldón, An. (Nys.) oswaldoi Peryassú as natural vectors of Plasmodium falciparum in French Guiana

Anopheles darlingi Root is the major vector of human malaria in the Neotropics and has been considered to be the sole malaria vector in French Guiana. The presence of other potential vectors suggests that malaria may be transmitted by other species under certain conditions. From 2006-2011, all anopheline specimens collected from 11 localities were assayed to determine if the Plasmodium circumsporozoite protein was present. In addition to An. darlingi, we found Anopheles oswaldoi, Anopheles intermedius and Anopheles nuneztovari specimens that were infected with Plasmodium sp. Further investigations on the behaviour and ecology of An. oswaldoi, An. intermedius and An. nuneztovari are necessary to determine their role in malaria transmission in French Guiana.     

The hidden diversity of Coleodactylus amazonicus (Sphaerodactylinae,Gekkota) revealed by molecular data

Coleodactylus amazonicus, a small leaf-litter diurnal gecko widely distributed in Amazon Basin, has been considered a single species with no significant morphological differences between populations along its range. A recent molecular study, however, detected large genetic differences between populations of central Amazonia and those in the easternmost part of the Amazon Basin, suggesting the presence of taxonomically unrecognised diversity. In this study, DNA sequences of three mitochondrial (16S, cytb, and ND4) and two nuclear genes (RAG-1, c-mos) were used to investigate whether the species currently identified as C.amazonicus contains morphologically cryptic species lineages. The present phylogenetic analysis reveals further genetic subdivision including at least five potential species lineages, restricted to northeastern (lineage A), southeastern (lineage B), central-northern (lineage E) and central-southern (lineages C and D) parts of Amazon Basin. All clades are characterized by exclusive groups of alleles for both nuclear genes and highly divergent mitochondrial haplotype clades, with corrected pairwise net sequence divergence between sister lineages ranging from 9.1% to 20.7% for the entire mtDNA dataset. Results of this study suggest that the real diversity of “C. amazonicus” has been underestimated due to its apparent cryptic diversification.     

Amazonian Amphibian Diversity Is Primarily Derived from Late Miocene Andean Lineages

The Neotropics contains half of remaining rainforests and Earth's largest reservoir of amphibian biodiversity. However, determinants of Neotropical biodiversity (i.e., vicariance, dispersals, extinctions, and radiations) earlier than the Quaternary are largely unstudied. Using a novel method of ancestral area reconstruction and relaxed Bayesian clock analyses, we reconstructed the biogeography of the poison frog clade (Dendrobatidae). We rejected an Amazonian center-of-origin in favor of a complex connectivity model expanding over the Neotropics. We inferred 14 dispersals into and 18 out of Amazonia to adjacent regions; the Andes were the major source of dispersals into Amazonia. We found three episodes of lineage dispersal with two interleaved periods of vicariant events between South and Central America. During the late Miocene, Amazonian, and Central American-Chocoan lineages significantly increased their diversity compared to the Andean and Guianan-Venezuelan-Brazilian Shield counterparts. Significant percentage of dendrobatid diversity in Amazonia and Chocó resulted from repeated immigrations, with radiations at <10.0 million years ago (MYA), rather than in situ diversification. In contrast, the Andes, Venezuelan Highlands, and Guiana Shield have undergone extended in situ diversification at near constant rate since the Oligocene. The effects of Miocene paleogeographic events on Neotropical diversification dynamics provided the framework under which Quaternary patterns of endemism evolved.     

Molecular phylogeography and the evolution and conservation of Amazonian mammals

The phylogeographic structure of 15 genera of Amazonian marsupials and rodents is summarized based on comparative sequence of the mitochondrial cytochrome b gene. The data are limited in geographical coverage, with samples widely scattered throughout Amazonia from the base of the Andes in Peru to the Guianan coast and eastern Brazil. We use this approach to define species boundaries, based minimally on the principle of reciprocal monophyly, in conjunction with morphological or other genetic discontinuities. The taxa so defined are older than previously appreciated, with many lineages dating from 1 to more than 3 Myr, and thus apparently predating the early Pleistocene. We relate patterns of concordant geographical shifts with underlying tectonic history and to current positions of major rivers. Finally, we provide comments on the utility of these data and patterns to conservation, articulating a need to incorporate phylogeographic information as part of the rationale in establishing conservation priorities at the organismal and geographical area levels.     

Cryptic lineages and Pleistocene population expansion in a Brazilian Cerrado frog

Diversification of South American species endemic to open habitats has been attributed to both Tertiary events and Pleistocene climatic fluctuations. Nonetheless, phylogeographical studies of taxa in these regions are few, precluding generalizations about the timing and processes leading to differentiation and speciation. We inferred population structure of Hypsiboas albopunctatus, a frog widely distributed in the Brazilian Cerrado. Three geographically distinct lineages were recovered in our phylogeny. The Chapada dos Guimarães (CG) clade was the first to diverge from other populations and contains multiple haplotypes from a single population in western Cerrado, probably representing a cryptic species. The southeast clade (SE) includes populations along the southeastern limit of the range within the historical distribution of the Brazilian Atlantic forest. Finally, the Central Cerrado (CC) group includes haplotypes from the interior of Brazil that are paraphyletic relative to the SE clade. Analyses of historical demography indicate significant population expansion in the CC and SE populations, likely associated with colonization of newly formed open habitats. The divergence of populations in the CG clade occurred in the late Miocene, concordant with the uplift of the central Brazilian plateau. Divergence of the SE clade from the CC occurred during the mid‐Pleistocene. Thus, both Tertiary geological events and Pleistocene climatic fluctuations promoted divergences among lineages. Our study reveals a complex history of diversification in the Cerrado, a morphoclimatic domain highly threatened because of anthropogenic habitat alteration. We identified surprisingly deep divergences in a widely distributed frog, indicating that the Cerrado is not a barrier‐free habitat and that its diversity is likely underestimated.     

A phylogeographic, demographic and historical analysis of the short-tailed pit viper (Gloydius brevicaudus): evidence for early divergence and late expansion during the Pleistocene

The impact of quaternary glaciation in eastern China on local fanua and flora has been a topic of considerable interest. We used mitochondrial DNA (mtDNA) sequence data and coalescent simulations to test two general biogeographic hypothesis related to the effects of Pleistocene climatic fluctuations for a widespread ophidian species (Gloydius brevicaudus) in eastern China and Korean Peninsula. The phylogenetic analysis revealed three major lineages, the southeast Coastal, Yangtze and North Lineages. The latter two are closely related and jointly form a continental lineage. Divergence dating and coalescent simulations indicate a Late Pliocene to Early Pleistocene divergence between lineages from the southeast coast and continental interior, followed by a mid-to-late Pleistocene divergence between lineages from the north and the middle-lower Yangtze Valley across East China, suggesting that all these lineages predated the last glacial maximum. An overlapping range between the two lineages within the continental lineage and a secondary contact associated with ecological transition zones on the margins of the North China Plain were also observed. These results show that vicariance patterns dominated the history of G. brevicaudus. Though the climatic events of the Pleistocene have had a marked effect on the historical distribution and intra-specific divergence of reptiles in China, coalescent and non-coalescent demographic analyses indicate that all lineages of G. brevicaudus seem not to have been adversely affected by glacial cycles during the Late Pleistocene, presumably because of an increase in the amount of climatically mild habitat in East Asia due to a decline in elevation and the development of monsoons since the Mid-End Pleistocene.     

Molecular evidence for a single taxon, Anopheles nuneztovari s.l., from two endemic malaria regions in Colombia

To elucidate the Anopheles nuneztovari s.l. taxonomic status at a microgeographic level in four malaria endemic localities from Antioquia and Córdoba, Colombia, fragments of the cytochrome oxidase subunit I (COI) and the white gene were used. The COI analysis showed low genetic differentiation with fixation index (F(ST)) levels between -0.02-0.137 and Nm values between 3-∞, indicating the presence of high gene flow among An. nuneztovari s.l. populations from the four localities. The COI network showed a single most common haplotype, type 1 (n = 55), present in all localities, as the likely ancestral haplotype. Analysis of the white gene showed that An. nuneztovari s.l. populations from both departments grouped with haplotypes 19 and 20, which are part of lineage 3 reported previously. The results of the present study suggest that An. nuneztovari s.l. is a single taxon in the area of the present study.     

Diversification by host switching and dispersal shaped the diversity and distribution of avian malaria parasites in Amazonia

Understanding how pathogens and parasites diversify through time and space is fundamental to predicting emerging infectious diseases. Here, we use biogeographic, coevolutionary and phylogenetic analyses to describe the origin, diversity, and distribution of avian malaria parasites in the most diverse avifauna on Earth. We first performed phylogenetic analyses using the mitochondrial cytochrome b (cyt b) gene to determine relationships among parasite lineages. Then, we estimated divergence times and reconstructed ancestral areas to uncover how landscape evolution has shaped the diversification of Parahaemoproteus and Plasmodium in Amazonia. Finally, we assessed the coevolutionary patterns of diversification in this host–parasite system to determine how coevolution may have influenced the contemporary diversity of avian malaria parasites and their distribution among Amazonian birds. Biogeographic analysis of 324 haemosporidian parasite lineages recovered from 4178 individual birds provided strong evidence that these parasites readily disperse across major Amazonian rivers and this has occurred with increasing frequency over the last five million years. We also recovered many duplication events within areas of endemism in Amazonia. Cophylogenetic analyses of these blood parasites and their avian hosts support a diversification history dominated by host switching. The ability of avian malaria parasites to disperse geographically and shift among avian hosts has played a major role in their radiation and has shaped the current distribution and diversity of these parasites across Amazonia.     

One Species, Three Pleistocene Evolutionary Histories: Phylogeography of the Italian Crested Newt, Triturus carnifex

Phylogeographic patterns of temperate species from the Mediterranean peninsulas have been investigated intensively. Nevertheless, as more phylogeographies become available, either unique patterns or new lines of concordance continue to emerge, providing new insights on the evolution of regional biotas. Here, we investigated the phylogeography and evolutionary history of the Italian crested newt, Triturus carnifex, through phylogenetic, molecular dating and population structure analyses of two mitochondrial gene fragments (ND2 and ND4; overall 1273 bp). We found three main mtDNA lineages having parapatric distribution and estimated divergence times between Late Pliocene and Early Pleistocene. One lineage (S) was widespread south of the northern Apennine chain and was further geographically structured into five sublineages, likely of Middle Pleistocene origin. The second lineage (C) was widespread throughout the Padano-Venetian plain and did not show a clear phylogeographic structure. The third lineage (N) was observed in only two populations located on western Croatia/Slovenia. Results of analysis of molecular variance suggested that partitioning populations according to the geographic distribution of these lineages and sublineages explains 76% of the observed genetic variation. The phylogeographic structure observed within T. carnifex and divergence time estimates among its lineages, suggest that responses to Pleistocene environmental changes in this single species have been as diverse as those found previously among several codistributed temperate species combined. Consistent with the landscape heterogeneity, physiographic features, and palaeogeographical evolution of its distribution range, these responses encompass multiple refugia along the Apennine chain, lowland refugia in large peri-coastal plains, and a 'cryptic' northern refugium.