- Research article
- Open Access
Unveiling an ancient biological invasion: molecular analysis of an old European alien, the crested porcupine (Hystrix cristata)
© Trucchi and Sbordoni; licensee BioMed Central Ltd. 2009
Received: 23 June 2008
Accepted: 18 May 2009
Published: 18 May 2009
Biological invasions can be considered one of the main threats to biodiversity, and the recognition of common ecological and evolutionary features among invaders can help developing a predictive framework to control further invasions. In particular, the analysis of successful invasive species and of their autochthonous source populations by means of genetic, phylogeographic and demographic tools can provide novel insights into the study of biological invasion patterns. Today, long-term dynamics of biological invasions are still poorly understood and need further investigations. Moreover, distribution and molecular data on native populations could contribute to the recognition of common evolutionary features of successful aliens.
We analyzed 2,195 mitochondrial base pairs, including Cytochrome b, Control Region and rRNA 12S, in 161 Italian and 27 African specimens and assessed the ancient invasive origin of Italian crested porcupine (Hystrix cristata) populations from Tunisia. Molecular coalescent-based Bayesian analyses proposed the Roman Age as a putative timeframe of introduction and suggested a retention of genetic diversity during the early phases of colonization. The characterization of the native African genetic background revealed the existence of two differentiated clades: a Mediterranean group and a Sub-Saharan one. Both standard population genetic and advanced molecular demography tools (Bayesian Skyline Plot) did not evidence a clear genetic signature of the expected increase in population size after introduction. Along with the genetic diversity retention during the bottlenecked steps of introduction, this finding could be better described by hypothesizing a multi-invasion event.
Evidences of the ancient anthropogenic invasive origin of the Italian Hystrix cristata populations were clearly shown and the native African genetic background was preliminary described. A more complex pattern than a simple demographic exponential growth from a single propagule seems to have characterized this long-term invasion.
The widespread introduction of non-native species has long been regarded as one of the major anthropogenic global changes and threats to biodiversity [1, 2]. The remarkable economic and ecological costs of biological invasions make the management of invaders one of the leading challenges in conservation biology [3, 4]. Standard methodologies implemented in population biology and new statistical tools recently developed in genetic analysis could help elucidate features and patterns relevant to invasive biology . Likewise, a comprehensive investigation of ecology and phylogeography of native source populations may suggest new insights into colonization and rapid evolutionary dynamics of invasive species . The scarcity of genetic data on the native range of successful aliens has recently been highlighted, along with the suggestion that further comparative studies could address these data in order to recognize some common evolutionary features among invaders .
Recently, a general warning on the importance of long-term perspective in invasion biology studies has been raised. Indeed, many researches have dealt with the acute phase of an invasion (i.e., the period soon after the introduction event), often lacking an adequate temporal context (see review in ). In particular, scientific efforts should be directed to better understand the long-term changes that occur in the environment and community of both invading and invaded species . In order to clarify the patterns and the evolutionary consequences of biological invasions, long-term perspective studies should focus on different aspects of the introduced populations, such as their genetic diversity and structure, population size at introduction (propagule pressure), growth rate and demographic trend . Considering the usefulness of molecular investigations in such retrospective studies , in this paper we focus on the analysis of genetic diversity, population structure and demographic trend of a putative ancient European invader, the crested porcupine Hystrix cristata.
Three species of the old world porcupines (Hystricidae, Rodentia), morphologically grouped in the Hystrix subgenus, occur all over Asia (H. indica) and Africa (H. cristata and H. africaeaustralis). The Cape porcupine (H. africaeaustralis) inhabits the South Africa region, from Cape of Good Hope to the borders of Central Africa tropical rainforests, while the crested porcupine (H. cristata) lives in North Central Africa, from the Mediterranean coast to Northern Zaire and Tanzania, and in mainland Italy and Sicily . These species are characterized by an extremely wide ecological tolerance: they can be found in arid grassland and semi-desert habitats, in Mediterranean shrub lands and deciduous forests, as well as on the Atlas and Kilimanjaro mountains, up to 3,500 m (a.s.l) of altitude. The origins of the Italian populations of H. cristata have long been debated. Based on the timescale of fossil and archaeological records  and on historical chronicles (Plinius, Naturalis Historia, VIII, 53), it was hypothesized that this species was introduced from Africa, in the late Roman Age, as a game animal. However, such hypothesis has recently been questioned, and a native European origin of the Italian porcupine has been suggested . A recent spread in North and Central-East Italy over the last century has been well documented, and the presence of archaeological and distributional records from areas where the species is now absent has also been recorded . Given a maximum density of 2–4 specimens/Km2 and the current area of occurrence, the population living in North-Central Italy can be estimated as ca. 80,000–120,000 individuals.
The main aim of this study is to highlight the geographical origin and the history of the Italian Hystrix cristata populations by means of genetic analyses, considering the hypothesis of their introduction in recent times. We also investigate the demographic parameters that characterize this putative invasion, on the basis of the genetic diversity and structure of the Italian populations. Moreover, preliminary genetic data on the African native range are shown.
The timescale of Italian invasion was inferred by means of coalescence-based Bayesian analyses, under two different tree priors which model population size changes through time. Only CR data were used in these analyses. The time to the most recent common ancestor (tMRCA) posterior densities of all introduced haplotypes and of two subgroups are shown in Figure 3b and Figure 3c. Given an intraspecific CR substitution rate of 0.4 subs/site/Myr (see Methods), the root of all samples ranges from 3.35 kyrbp (0.6 – 7.3 kyrbp HPD 95%), under the Cons_s model, to 2.9 kyrbp (0.4 – 6.7 kyrbp HPD 95%), under the Expa_g model.
Our results confirm the ancient non-native origin of the Italian H. cristata populations: i) the summary molecular statistics, derived from three mitochondrial genes, highlights the African genetic diversity as being more complex when compared with the Italian; ii) the phylogenetic analysis shows the close relationship between the Italian and the Tunisian populations; iii) the median-joining network analysis suggests that a Tunisian and the most widespread Italian haplotypes share a common ancestor. These results indicate the Tunisian area as being the most probable source of the Italian introduced populations. In our sample the native geographic range is represented by only five individuals with a global native range FST of 0.66. This supports the assignation of introduced haplotypes to a particular source region with a confidence index of 0.95 (accuracy at regional level of 0.9) as shown by Muirhead et al. . The median-joining network analysis also reveals the simple genetic structure of the introduced Italian populations i.e., a single mitogenetic clade of strictly related haplotypes. Both Central-North Italy and the Sicilian populations share the common haplotypes (Ita-A and Ita-B), while Ita-B seems to be absent from southern Italy.
Archaeological records and historical chronicles (Plinius, Naturalis Historia, VIII, 53) do not support an introduction of this species in Italy before the late Roman Empire. The first calibrated sub-fossil of H. cristata was found in South Italy (Basilicata) and dates back to 1.5 kyrbp . Since we found a close relationship between the Italian and the Tunisian haplotypes, a putative time frame for introduction could span from the early colonization of Sicily by the Phoenicians (ca. 2.5 kyrbp)  to the first Italian sub-fossil record. In the light of these considerations, H. cristata could have been introduced in Sicily or in mainland Italy during this 1 ky-time frame as a game species, as threatening beasts in circus games or as exotic animals in rich Roman country homes. Considering the coalescent-based Bayesian approach, the time-measured phylogeny of introduced haplotypes partially agrees with this historical framework. However, even if the molecular clock was calibrated using a very fast intraspecific substitution rate [19, 20], the mean value of Italian haplotypes tMRCA was dated back to ca. 3 kyrbp, farther than the above hypothesized time span. Conversely, the two recognized Italian subgroups seem to have evolved in situ, during the early phases of colonization and spread (2.1–1.1 kyrbp). Nevertheless, this scenario would need to be supported by a rare event of molecular diversity retention in the bottlenecked population during the founding event , since genetic drift reduces genetic variation in small populations . Likewise, since Central-North Italy and Sicily share the common Ita-A and Ita-B haplotypes, this rare genetic event during introduction might have occurred in both regions. The possibility of a further acceleration of the molecular substitution rates in extremely recent events as this could also be taken into account . In fact, a time-dependency of the molecular clock rates has been recently proposed and discussed [24, 25]. As for H. cristata introduction, a substitution rate of 0.5 – 0.6 subs/site/Myr in the mtDNA CR has to be accepted if we consider that the whole Italian molecular diversity was generated in situ.
A review of the most recent literature of the demographic parameters that characterize the colonization dynamics does not provide clear evidence of a common pattern during colonization. A successful invasion can take place as a multiple introduction of few individuals from genetically differentiated native populations [26, 27], or as a single founder event from extremely bottlenecked gene pools . Being an invasive species, the present population size of Italian H. cristata is the result of a demographic and spatial expansion though the genetic signature of this remarkable growth is not revealed by molecular statistics. Although there is evidence of an extremely recent demographic increase, the coalescent-based Bayesian analysis does not clearly support a population expansion model. Since this demographic analysis is based on the assumption of a single panmittic population at mutation/drift equilibrium, the results may be better explained by hypothesizing more than one event of introduction of isolated propagules.
As for the genetic background of H. cristata native range, the phylogenetic reconstruction highlighted the existence of two different haplogroups, Sub-Saharan and Mediterranean. The high level of genetic differentiation of the Sub-Saharan clade advocates for a complex evolutionary history, probably related to the broad geographical area of occurrence. Indeed, the haplotypes belonging to this clade were sampled in Central Africa as well as in Tunisia and Libya. The close relationship between the Tun-B/Lib-A and the Eritrean haplotypes suggests the existence of an ecological connection going from East to North Africa, encompassing the Sahara region. This hypothesis is supported by strong evidence of humid phases characterizing North African climate (15-8 kyrbp) . As for the second haplogroup, a low level of genetic differentiation can be detected. This could be the consequence of a less complex evolutionary history or a more recent diversification from a common ancestor when compared with the Sub-Saharan haplogroup. However, the evolution of native range H. cristata should be analyzed along with the environmental changes which occurred in Africa over the Late Pleistocene and the Holocene . As widely acknowledged in the Palaearctic Region , the distribution and the phylogeographical structure of several African taxa have largely been shaped by these past climatic shifts [32, 33]. Further investigation in the North African native range would help to clarify the relationship between the evolutionary history of the African populations and the high invasive potential of the propagules coming from such populations.
Molecular data strongly evidenced the invasive origin of the Italian Hystrix cristata populations from North Africa (Tunisia) introduced during the Roman Age. Although demographic analyses failed to support a population size expansion from a single introduced propagule, they suggest a different scenario in favour of a multi-invasion event. Moreover, this study produced preliminary data on the genetic diversity of native range populations, and provides a useful contribution for further comparative studies on successful aliens.
A total of 179 H. cristata samples (blood and quills) were collected in the field and analyzed between 2006 and 2007: 161 samples came from all over Italy and 18 samples from different African populations in Burkina Faso, Eritrea, Ethiopia Libya, Morocco, Sierra Leone, Tanzania, and Tunisia. Moreover, nine H. africaeaustralis samples from Namibia, South Africa and Zambia were analyzed as representatives of the sister taxon species. One H. indica specimen from Israel was added in the phylogenetic analyses as outgroup. Location of specimens used in the statistical analyses are given in Figure 1 [see Additional file 2].
Sequences and references of primers used in this study.
The sequences were edited and aligned using BIOEDIT. Haplotype diversity Hd, segregating sites S, nucleotide diversity averaged per site π  and, Synonymous/Non-Synonymous ratio (as regards cytochrome b gene) were calculated using DNASP. The same software was used to describe number and frequency of different haplotypes and to calculate θ per sequence from segregating sites , being this diversity value an informative synopsis of the metrics used. Net genetic distances based on the proposed GTR + Γ + I  model of sequence evolution were calculated on the whole analyzed mtDNA fragment, exporting tip-to-tip distances obtained from Maximum-Likelihood phylogeny calculated in TREEFINDER. A Maximum-Likelihood (ML) and a Bayesian (Ba) methods were used to assess the phylogenetic relationships among mitochondrial haplotypes. TREEFINDER was used to evaluate the best-fit model of sequence evolution for each analyzed data set, using the Akaike Information Criterion. As regards the whole dataset, the proposed model was a GTR + Γ + I with unequal base frequencies. TREEFINDER was used to build the ML phylogeny, while BEAST was used to reconstruct the Ba phylogeny. Robustness of the ML trees was tested with 1,000 bootstrap replicates. The relationship among the haplotypes belonging to each tree sub-clade was reconstructed by means of a Median-Joining network analysis, as implemented in NETWORK.
According to Muirhead et al. , global Fst among native populations was used to assess the confidence index to classify introduced individuals to a source geographic area. It could prevent incorrect identification of invasion source when a low number of individuals has been sampled in putative source populations. Three putative African sources were considered: Tunisia-Libya, Morocco, and East Africa comprising Eritrean, Ethiopian and Tanzanian samples. A Kimura 2-parameters distance method and 10,000 permutations was used to calculate global Fst with the software ARLEQUIN.
TMRCAs of all Italian haplotypes and of 2 sub-clades were estimate on the base of Control Region data using BEAST. Constant_size (Cons_s) and Expansion_growth (Expa_g) tree priors were compared as well as the Relaxed and the Strict molecular clock models. At least three independent runs 3 × 106 steps long were performed for each dataset. Convergence of chains, effective sample size, estimates and credible intervals for each parameter were analyzed with the program TRACER. TMRCAs were dated using a 0.4 subs/site/Myr substitution rate. This rate has been proposed by Rajabi-Maham et al.  for intraspecific polymorphism analyses in Mus musculus domesticus. Since similar mitochondrial substitution rates were recorded in Hystricidae and Muridae , this molecular rate was also considered in our study.
The mtDNA sequences were analyzed in order to find traces of demographic trend from introduction, growth rate and present population size. Several approaches were attempted. These analyses were performed only on the North Central Italy samples, which can be considered, to a large extent, as a panmittic expanding population. Only the fast-evolving mtDNA CR-I was used in these investigations. Two different neutrality tests, which have been shown to be the most powerful , were used to investigate for population expansion or bottleneck on the sampled CR-I sequences: Fu's Fs  statistics were calculated, as implemented in ARLEQUIN and R2  as implemented in DNASP. Significance of Fs was assessed by 10,000 randomization, while significance of R2 was determined by means of 10,000 coalescent simulation on the basis of observed number of segregating sites in the sample. The mismatch distribution of pairwise differences between sequences was examined to find genetic traces of the demographic trend. The sum-of-squared-differences (SSD) statistic was used to test the goodness-of-fit between the observed mismatch distribution and that expected under a sudden expansion model. The significance of SSD was assessed by 10,000 parametric bootstrap re-sampling , using ARLEQUIN. Past population demography was also investigated with a Bayesian Skyline Plot (BSP) , as implemented in BEAST. This coalescent-based approach calculates the effective breeding population size (Nef) through time directly from sampled sequence data and not from previously recognized phylogeny. This procedure can therefore account for uncertainty associated with reconstructed phylogeny. The analysis was done using a different time segmentation ranging from 5 to 10 groups (past Nef points) and substitution model (HKY + I) . Three independent runs of 50 × 106 iterations for each grouping scenarios were performed. In addition, separate runs were performed, using different coalescent tree priors which model population size through time: Cons_s, Expa_g and Exponential_growth (Expo_g). Three independent runs of 30 × 106 iterations for each scenario were done. Convergence of chains, effective sample size, estimates and credible intervals for each parameters and demographic reconstructions were analyzed with the software TRACER. The three different demographic models were then compared to each other in order to assess their relative fit to the data. Comparisons were made under TRACER workspace evaluating the marginal likelihood ratio (marginal with respect to the tree prior) of pairs of models (Bayes Factor calculation). The best approximation to the marginal likelihood comparison was found calculating the Bayes Factor on the tree Likelihood trace [52, 53].
This research was granted by Regione Lazio, Provincia di Roma and Osservatorio della Biodiversità del Lazio. Thanks to the kind collaboration of H. Azafzaf, R. Giannatelli, J. Shoshawi, G. Chiozzi, J.L. Ruiz, B. Massa, M.A. Bologna, G. Montinaro, M. Pitzalis, M. Di Cecco, M. Pellegrini, M. Visceglia, C. Bernardi, E. Aloisi Masella, A. Bonacci, U. Veken, C. Settanni, P. Parisi, A. Coluccini, E. Genovesi, F. Viviani, G. Speroni, N. Raffaelli, R. Vannucci, S. Elio, G. Cavalloni, G. Bertola, A. Coluccini, P. Politi, D'Amelia, Carotenuto, F. Scarfò, D'Amato, A. Romano, S. Marta, M. Di Rao, M. Lucarelli, M. Mattoccia, G. Gentile, Baragio, Marchionni, Schiavano, S. De Felici, A. Fabiani, A. Marroni, Todini, F. Franconi, M. Cancellieri, G. Riccarducci, S. Sanrocco, F. Terragni, A. El Sehhar, G. Tito, A. Campanaro, M. Di Domenico, G.C. Bozano, ARP Lazio and RomaNatura that have provided the biological material. We wish to thank A. Fabiani and M. Lo Ponte for the English language revision and G. Gentile and P. Gratton for their helpful comments.
- Mooney HA, Hobbs RJ: Invasive Species in a Changing World. 2000, Washington, D.C.: Island PressGoogle Scholar
- Lodge DM, Shrader-Frechette K: Nonindigenous species: ecological explanation, environmental ethics, and public policy. Conserv Biol. 2003, 17: 31-37. 10.1046/j.1523-1739.2003.02366.x.View ArticleGoogle Scholar
- Allendorf FW, Lundquist LL: Introduction: Population biology, evolution, and control of invasive species. Conserv Biol. 2003, 17: 24-30. 10.1046/j.1523-1739.2003.02365.x.View ArticleGoogle Scholar
- Shirley SM, Kark S: Amassing efforts against alien invasive species in Europe. PLoS Biol. 2006, 4 (8): e279-10.1371/journal.pbio.0040279. DOI: 10.1371/journal.pbio.0040279PubMed CentralView ArticlePubMedGoogle Scholar
- Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, McCauley DE, O'Neil P, Parker IM, Thompson JN, Weller SG: The population biology of invasive species. Annu Rev Ecol Syst. 2001, 32: 305-332. 10.1146/annurev.ecolsys.32.081501.114037.View ArticleGoogle Scholar
- Slade RW, Moritz C: Phylogeography of Bufo marinus from its natural and introduced ranges. Proc R Soc Lond B Biol Sci. 1998, 265: 769-777. 10.1098/rspb.1998.0359.View ArticleGoogle Scholar
- Suarez AV, Tsutsui ND: The evolutionary consequences of biological invasions. Mol Ecol. 2008, 17: 351-360. 10.1111/j.1365-294X.2007.03456.x.View ArticlePubMedGoogle Scholar
- Strayer DL, Eviner VT, Jeschke JM, Pace ML: Understanding the long-term effects of species invasions. Trends Ecol Evol. 2006, 21: 645-651. 10.1016/j.tree.2006.07.007.View ArticlePubMedGoogle Scholar
- Strauss SY, Lau JA, Carroll SP: Evolutionary responses of natives to introduced species: what do introductions tell us about natural communities?. Ecol Lett. 2006, 9: 357-374. 10.1111/j.1461-0248.2005.00874.x.View ArticlePubMedGoogle Scholar
- Avise JC: Phylogeography: The History and Formation of Species. 2000, Cambridge, Massachusetts: Harvard University PressGoogle Scholar
- Kingdon J: The Kingdon Field Guide to African Mammals. 1997, London: Academic pressGoogle Scholar
- Riquelme JA, Morales A: A porcupine find from roman Africa with a review of archaeozoological data from circummediterranean sites. Archaeofauna. 1997, 6: 91-95.Google Scholar
- Angelici FM, Capizzi D, Amori G, Luiselli L: Morphometric variation in the skulls of the crested porcupine Hystrix cristata from mainland Italy, Sicily and northern Africa. Mamm Biol. 2003, 68: 165-173. 10.1078/1616-5047-00078.Google Scholar
- Angelici FM, Amori G: Distribution of the Crested porcupine Hystrix cristata L., 1758 in peninsular Italy and Sicily. Mus Reg Sci Nat Boll (Torino). 1999, 16: 83-88.Google Scholar
- Niethammer J: Hystrix cristata Linnaeus, 1758 – Stachelschwein. Handbuck der Saugetiere Europas 2/I. Rodentia II. Edited by: Niethammer J, Krapp F. 1982, Wiesbaden, Academische Verlagsgesellschaft, 588-605.Google Scholar
- Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.PubMed CentralPubMedGoogle Scholar
- Muirhead JR, Gray DK, Kelly DW, Ellis SM, Heath DD, MacIsaac HJ: Identifying the source of species invasions: sampling intensity vs. genetic diversity. Mol Ecol. 2008, 17: 1020-1035. 10.1111/j.1365-294X.2008.03669.x.View ArticlePubMedGoogle Scholar
- Leighton R: Sicily before History: An Archaeological Survey from the Palaeolithic to the Iron Age. 1999, Ithaca, N.Y., Cornell University PressGoogle Scholar
- Rajabi-Maham H, Orth A, Bonhomme F: Phylogeography and postglacial expansion of Mus musculus domesticus inferred from mitochondrial DNA coalescent, from Iran to Europe. Mol Ecol. 2008, 17: 627-641.View ArticlePubMedGoogle Scholar
- Saarma U, Ho SYW, Pybus OG, Kaljuste M, Tumanov IL, Kojola I, Vorobiev AA, Markov NI, Saveljev AP, Valdmann H, Lyapunova EA, Abramov AV, Männil P, Korsten M, Vulla E, Pazetnov SV, Pazetnov VS, Putchkovskiy SV, Rõkov AM: Mitogenetic structure of brown bears (Ursus arctos L.) in northeastern Europe and a new time frame for the formation of European brown bear lineages. Mol Ecol. 2007, 16: 401-413. 10.1111/j.1365-294X.2006.03130.x.View ArticlePubMedGoogle Scholar
- Abdelkrim J, Pascal M, Samadi S: Island colonization and founder effects: the invasion of the Guadeloupe islands by ship rats (Rattus rattus). Mol Ecol. 2005, 14: 2923-2931. 10.1111/j.1365-294X.2005.02604.x.View ArticlePubMedGoogle Scholar
- Tsutsui ND, Suarez AV, Holway DA, Case TJ: Reduced genetic variation and the success of an invasive species. Proc Natl Acad Sci USA. 2000, 97: 5948-5953. 10.1073/pnas.100110397.PubMed CentralView ArticlePubMedGoogle Scholar
- Pulquerio MJF, Nichols RA: Dates from the molecular clock: how wrong can we be?. Trends Ecol Evol. 2007, 22: 180-184. 10.1016/j.tree.2006.11.013.View ArticlePubMedGoogle Scholar
- Ho SYW, Phillips MJ, Cooper A, Drummond AJ: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol. 2005, 22: 1561-1568. 10.1093/molbev/msi145.View ArticlePubMedGoogle Scholar
- Shapiro B, Drummond AJ, Rambaut A, Wilson MC, Matheus PE, Sher AV, Pybus OG, Gilbert MTP, Barnes I, Binladen J, Willerslev E, Hansen AJ, Baryshnikov GF, Burns JA, Davydov S, Driver JC, Froese DG, Harington CR, Keddie G, Kosintsev P, Kunz ML, Martin LD, Stephenson RO, Storer J, Tedford R, Zimov S, Cooper A: Rise and fall of the Beringian steppe bison. Science. 2004, 306: 1561-1565. 10.1126/science.1101074.View ArticlePubMedGoogle Scholar
- Bonhomme M, Blancher A, Cuartero S, Chikhi L, Crouau-Roy B: Origin and number of founders in an introduced insular primate: estimation from nuclear genetic data. Mol Ecol. 2008, 17: 1009-1019. 10.1111/j.1365-294X.2007.03645.x.View ArticlePubMedGoogle Scholar
- Lavergne S, Molofsky J: Increased genetic variation and evolutionary potential drive the success of an invasive grass. Proc Natl Acad Sci USA. 2007, 104: 3883-3888. 10.1073/pnas.0607324104.PubMed CentralView ArticlePubMedGoogle Scholar
- Puillandre N, Dupas S, Dangles O, Zeddam JL, Capdevielle-Dulac C, Barbin K, Torres-Leguizamon M, Silvain JF: Genetic bottleneck in invasive species: the potato tuber mothadds to the list. Biol Invasions. 2008, 10: 319-333. 10.1007/s10530-007-9132-y.View ArticleGoogle Scholar
- Gasse F: Hydrological changes in the African tropics since the Last Glacial Maximum. Quat Sci Rev. 2000, 19: 189-211. 10.1016/S0277-3791(99)00061-X.View ArticleGoogle Scholar
- deMenocal PB: African climate change and faunal evolution during the Pliocene-Pleistocene. Earth Planet Sci Lett. 2004, 220: 3-24. 10.1016/S0012-821X(04)00003-2.View ArticleGoogle Scholar
- Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Flagstad Ø, Syvertsen PO, Stenseth NC, Jakobsen KS: Environmental change and rate of evolution: the phylogeographic pattern within the hartebeest complex as related to climatic variation. Proc R Soc Lond B Biol Sci. 2001, 268: 667-677. 10.1098/rspb.2000.1416.View ArticleGoogle Scholar
- Cosson JF, Hutterer R, Libois R, Sara M, Taberlet P, Vogel P: Phylogeographical footprints of the Strait of Gibraltar and Quaternary climatic fluctuations in the Western Mediterranean: a case study with the greater white-toothed shrew, Crocidura russula (Mammalia: Soricidae). Mol Ecol. 2005, 14: 1151-1162. 10.1111/j.1365-294X.2005.02476.x.View ArticlePubMedGoogle Scholar
- Trucchi E, Gentile G, Sbodoni V: Development of primers to amplify mtDNA control region of Old World porcupines (subgenus Hystrix). Mol ecol resour. 2008, 8 (5): 1139-1141. 10.1111/j.1755-0998.2008.02201.x.View ArticlePubMedGoogle Scholar
- Nedbal MA, Allard MW, Honeycutt RL: Molecular systematic of Hystricognathi rodents: evidence from the mitochondrial 12S rRNA gene. Mol Phylogenet Evol. 1994, 3: 206-220. 10.1006/mpev.1994.1023.View ArticlePubMedGoogle Scholar
- Pompanon F, Bonin A, Bellemain E, Taberlet P: Genotyping errors: causes, consequences and solutions. Nat Rev Gen. 2005, 6: 847-846. 10.1038/nrg1707.View ArticleGoogle Scholar
- Hall TA: bioedit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, New York, Columbia University PressGoogle Scholar
- Rozas J, Sánchez-Del Barrio JC, Messeguer X, Rozas R: dnasp, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7: 256-276. 10.1016/0040-5809(75)90020-9.View ArticlePubMedGoogle Scholar
- Tavaré S: Some probabilistic and statistical problems in the analysis of DNA sequences. Lect Math Life Sc. 1986, 17: 57-86.Google Scholar
- Jobb G, von Haesler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetic. BMC Evol Biol. 2004, 4: 18-10.1186/1471-2148-4-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: beast: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Spradling TA, Hafner MS, Demastes JW: Differences in rate of Cytochrome-b evolution among species of rodents. J Mammal. 2001, 82: 65-80. 10.1644/1545-1542(2001)082<0065:DIROCB>2.0.CO;2.View ArticleGoogle Scholar
- Ramons-Onsins SE, Rozas J: Statistical Properties of New Neutrality Tests Against Population Growth. Mol Biol Evol. 2002, 19: 2092-2100.View ArticleGoogle Scholar
- Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking, and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Schneider S, Excoffier L: Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to humanmitochondrial DNA. Genetics. 1999, 152: 1079-1089.PubMed CentralPubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian Coalescent Inference of Past Population Dynamics from Molecular Sequences. Mol Biol Evol. 2005, 22: 1185-1192. 10.1093/molbev/msi103.View ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.View ArticlePubMedGoogle Scholar
- Suchard MA, Weiss RE, Sinsheimer JS: Bayesian Selection of Continuous-Time Markov Chain Evolutionary Models. Mol Biol Evol. 2001, 18: 1001-1013.View ArticlePubMedGoogle Scholar
- beast user: A general user list for discussion about the BEAST MCMC program. [http://groups.google.com/group/beast-users?hl=en]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.