Bayesian estimation of post-Messinian divergence times in Balearic island lizards

R.P. Brown, B. Terrasa, V. Pérez-Mellado, J.A. Castro, P. Hoskisson, A. Picornell, M.M. Ramon

Research output: Contribution to journalArticle

51 Citations (Scopus)

Abstract

Phylogenetic relationships and timings of major cladogenesis events are investigated in the Balearic Island lizards Podarcis lilfordi and P. pityusensis using 2675 bp of mitochondrial and nuclear DNA sequences. Partitioned Bayesian and Maximum Parsimony analyses provided a well-resolved phylogeny with high node-support values. Bayesian MCMC estimation of node dates was investigated by comparing means of posterior distributions from different subsets of the sequence against the most robust analysis which used multiple partitions and allowed for rate heterogeneity among branches under a rate-drift model. Evolutionary rates were systematically underestimated and thus divergence times overestimated when sequences containing lower numbers of variable sites were used (based on ingroup node constraints). The following analyses allowed the best recovery of node times under the constant-rate (i.e., perfect clock) model: (i) all cytochrome b sequence (partitioned by codon position), (ii) cytochrome b (codon position 3 alone), (iii) NADH dehydrogenase (subunits 1 and 2; partitioned by codon position), (iv) cytochrome b and NADH dehydrogenase sequence together (six gene-codon partitions), (v) all unpartitioned sequence, (vi) a full multipartition analysis (nine partitions). Of these, only (iv) and (vi) performed well under the rate-drift model. These findings have significant implications for dating of recent divergence times in other taxa. The earliest P. lilfordi cladogenesis event (divergence of Menorcan populations), occurred before the end of the Pliocene, some 2.6 Ma. Subsequent events led to aWest Mallorcan lineage (2.0 Ma ago), followed 1.2 Ma ago by divergence of populations from the southern part of the Cabrera archipelago from a widely-distributed group from north Cabrera, northern and southern Mallorcan islets. Divergence within P. pityusensis is more recent with the main Ibiza and Formentera clades sharing a common ancestor at about 1.0 Ma ago. Climatic and sea level changes are likely to have initiated cladogenesis, with lineages making secondary contact during periodic landbridge formation. This oscillating cross-archipelago pattern in which ancient divergence is followed by repeated contact resembles that seen between East-West refugia populations from mainland Europe.
LanguageEnglish
Pages350-358
Number of pages8
JournalMolecular Phylogenetics and Evolution
Volume48
Issue number1
DOIs
Publication statusPublished - 2008

Fingerprint

Balearic Islands
Genetic Speciation
Lizards
Messinian
codons
lizard
Codon
Spain
Cytochromes b
lizards
divergence
cytochrome b
NADH Dehydrogenase
NADH dehydrogenase
cytochrome
Population
secondary contact
archipelago
phylogeny
Phylogeny

Keywords

  • Bayesian MCMC
  • messinian salinity crisis
  • mediterranean
  • mtDNA
  • phylogeny
  • phylogeography
  • podarcis
  • rate-drift

Cite this

Brown, R. P., Terrasa, B., Pérez-Mellado, V., Castro, J. A., Hoskisson, P., Picornell, A., & Ramon, M. M. (2008). Bayesian estimation of post-Messinian divergence times in Balearic island lizards. Molecular Phylogenetics and Evolution, 48(1), 350-358. https://doi.org/10.1016/j.ympev.2008.04.013
Brown, R.P. ; Terrasa, B. ; Pérez-Mellado, V. ; Castro, J.A. ; Hoskisson, P. ; Picornell, A. ; Ramon, M.M. / Bayesian estimation of post-Messinian divergence times in Balearic island lizards. In: Molecular Phylogenetics and Evolution. 2008 ; Vol. 48, No. 1. pp. 350-358.
@article{026d3f041f5b470b8cd3257be1207b22,
title = "Bayesian estimation of post-Messinian divergence times in Balearic island lizards",
abstract = "Phylogenetic relationships and timings of major cladogenesis events are investigated in the Balearic Island lizards Podarcis lilfordi and P. pityusensis using 2675 bp of mitochondrial and nuclear DNA sequences. Partitioned Bayesian and Maximum Parsimony analyses provided a well-resolved phylogeny with high node-support values. Bayesian MCMC estimation of node dates was investigated by comparing means of posterior distributions from different subsets of the sequence against the most robust analysis which used multiple partitions and allowed for rate heterogeneity among branches under a rate-drift model. Evolutionary rates were systematically underestimated and thus divergence times overestimated when sequences containing lower numbers of variable sites were used (based on ingroup node constraints). The following analyses allowed the best recovery of node times under the constant-rate (i.e., perfect clock) model: (i) all cytochrome b sequence (partitioned by codon position), (ii) cytochrome b (codon position 3 alone), (iii) NADH dehydrogenase (subunits 1 and 2; partitioned by codon position), (iv) cytochrome b and NADH dehydrogenase sequence together (six gene-codon partitions), (v) all unpartitioned sequence, (vi) a full multipartition analysis (nine partitions). Of these, only (iv) and (vi) performed well under the rate-drift model. These findings have significant implications for dating of recent divergence times in other taxa. The earliest P. lilfordi cladogenesis event (divergence of Menorcan populations), occurred before the end of the Pliocene, some 2.6 Ma. Subsequent events led to aWest Mallorcan lineage (2.0 Ma ago), followed 1.2 Ma ago by divergence of populations from the southern part of the Cabrera archipelago from a widely-distributed group from north Cabrera, northern and southern Mallorcan islets. Divergence within P. pityusensis is more recent with the main Ibiza and Formentera clades sharing a common ancestor at about 1.0 Ma ago. Climatic and sea level changes are likely to have initiated cladogenesis, with lineages making secondary contact during periodic landbridge formation. This oscillating cross-archipelago pattern in which ancient divergence is followed by repeated contact resembles that seen between East-West refugia populations from mainland Europe.",
keywords = "Bayesian MCMC, messinian salinity crisis, mediterranean, mtDNA, phylogeny, phylogeography, podarcis, rate-drift",
author = "R.P. Brown and B. Terrasa and V. P{\'e}rez-Mellado and J.A. Castro and P. Hoskisson and A. Picornell and M.M. Ramon",
year = "2008",
doi = "10.1016/j.ympev.2008.04.013",
language = "English",
volume = "48",
pages = "350--358",
journal = "Molecular Phylogenetics and Evolution",
issn = "1055-7903",
number = "1",

}

Brown, RP, Terrasa, B, Pérez-Mellado, V, Castro, JA, Hoskisson, P, Picornell, A & Ramon, MM 2008, 'Bayesian estimation of post-Messinian divergence times in Balearic island lizards' Molecular Phylogenetics and Evolution, vol. 48, no. 1, pp. 350-358. https://doi.org/10.1016/j.ympev.2008.04.013

Bayesian estimation of post-Messinian divergence times in Balearic island lizards. / Brown, R.P.; Terrasa, B.; Pérez-Mellado, V.; Castro, J.A.; Hoskisson, P.; Picornell, A.; Ramon, M.M.

In: Molecular Phylogenetics and Evolution, Vol. 48, No. 1, 2008, p. 350-358.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Bayesian estimation of post-Messinian divergence times in Balearic island lizards

AU - Brown, R.P.

AU - Terrasa, B.

AU - Pérez-Mellado, V.

AU - Castro, J.A.

AU - Hoskisson, P.

AU - Picornell, A.

AU - Ramon, M.M.

PY - 2008

Y1 - 2008

N2 - Phylogenetic relationships and timings of major cladogenesis events are investigated in the Balearic Island lizards Podarcis lilfordi and P. pityusensis using 2675 bp of mitochondrial and nuclear DNA sequences. Partitioned Bayesian and Maximum Parsimony analyses provided a well-resolved phylogeny with high node-support values. Bayesian MCMC estimation of node dates was investigated by comparing means of posterior distributions from different subsets of the sequence against the most robust analysis which used multiple partitions and allowed for rate heterogeneity among branches under a rate-drift model. Evolutionary rates were systematically underestimated and thus divergence times overestimated when sequences containing lower numbers of variable sites were used (based on ingroup node constraints). The following analyses allowed the best recovery of node times under the constant-rate (i.e., perfect clock) model: (i) all cytochrome b sequence (partitioned by codon position), (ii) cytochrome b (codon position 3 alone), (iii) NADH dehydrogenase (subunits 1 and 2; partitioned by codon position), (iv) cytochrome b and NADH dehydrogenase sequence together (six gene-codon partitions), (v) all unpartitioned sequence, (vi) a full multipartition analysis (nine partitions). Of these, only (iv) and (vi) performed well under the rate-drift model. These findings have significant implications for dating of recent divergence times in other taxa. The earliest P. lilfordi cladogenesis event (divergence of Menorcan populations), occurred before the end of the Pliocene, some 2.6 Ma. Subsequent events led to aWest Mallorcan lineage (2.0 Ma ago), followed 1.2 Ma ago by divergence of populations from the southern part of the Cabrera archipelago from a widely-distributed group from north Cabrera, northern and southern Mallorcan islets. Divergence within P. pityusensis is more recent with the main Ibiza and Formentera clades sharing a common ancestor at about 1.0 Ma ago. Climatic and sea level changes are likely to have initiated cladogenesis, with lineages making secondary contact during periodic landbridge formation. This oscillating cross-archipelago pattern in which ancient divergence is followed by repeated contact resembles that seen between East-West refugia populations from mainland Europe.

AB - Phylogenetic relationships and timings of major cladogenesis events are investigated in the Balearic Island lizards Podarcis lilfordi and P. pityusensis using 2675 bp of mitochondrial and nuclear DNA sequences. Partitioned Bayesian and Maximum Parsimony analyses provided a well-resolved phylogeny with high node-support values. Bayesian MCMC estimation of node dates was investigated by comparing means of posterior distributions from different subsets of the sequence against the most robust analysis which used multiple partitions and allowed for rate heterogeneity among branches under a rate-drift model. Evolutionary rates were systematically underestimated and thus divergence times overestimated when sequences containing lower numbers of variable sites were used (based on ingroup node constraints). The following analyses allowed the best recovery of node times under the constant-rate (i.e., perfect clock) model: (i) all cytochrome b sequence (partitioned by codon position), (ii) cytochrome b (codon position 3 alone), (iii) NADH dehydrogenase (subunits 1 and 2; partitioned by codon position), (iv) cytochrome b and NADH dehydrogenase sequence together (six gene-codon partitions), (v) all unpartitioned sequence, (vi) a full multipartition analysis (nine partitions). Of these, only (iv) and (vi) performed well under the rate-drift model. These findings have significant implications for dating of recent divergence times in other taxa. The earliest P. lilfordi cladogenesis event (divergence of Menorcan populations), occurred before the end of the Pliocene, some 2.6 Ma. Subsequent events led to aWest Mallorcan lineage (2.0 Ma ago), followed 1.2 Ma ago by divergence of populations from the southern part of the Cabrera archipelago from a widely-distributed group from north Cabrera, northern and southern Mallorcan islets. Divergence within P. pityusensis is more recent with the main Ibiza and Formentera clades sharing a common ancestor at about 1.0 Ma ago. Climatic and sea level changes are likely to have initiated cladogenesis, with lineages making secondary contact during periodic landbridge formation. This oscillating cross-archipelago pattern in which ancient divergence is followed by repeated contact resembles that seen between East-West refugia populations from mainland Europe.

KW - Bayesian MCMC

KW - messinian salinity crisis

KW - mediterranean

KW - mtDNA

KW - phylogeny

KW - phylogeography

KW - podarcis

KW - rate-drift

UR - http://www.ljmu.ac.uk/BIE/BIE_docs/BROWNETAL2008.pdf

UR - http://dx.doi.org/10.1016/j.ympev.2008.04.013

U2 - 10.1016/j.ympev.2008.04.013

DO - 10.1016/j.ympev.2008.04.013

M3 - Article

VL - 48

SP - 350

EP - 358

JO - Molecular Phylogenetics and Evolution

T2 - Molecular Phylogenetics and Evolution

JF - Molecular Phylogenetics and Evolution

SN - 1055-7903

IS - 1

ER -