BACKGROUND
the pilbara in northwestern australia is an ancient geological and biophysical region  that has repeatedly been identified as an important zone of biological diversity, endemism and past refugia based on congruent biogeographic patterns of both fauna and flora  <cit> .
our understanding of the biotic history of the region is still in its infancy , but the geology of the pilbara region is well understood  <cit>  and may serve as an important predictor of biological diversity.
the pilbara is underlain by sedimentary and igneous rocks ranging up to  <dig>  billion years  in age, and overlain by some of the most ancient erosion surfaces on earth  <cit> .
the pilbara craton in the northeastern portion of the region has two parts; archaean  granites and metamorphosed volcanic rocks  in the north, and, stratigraphically overlying these rocks in the south, is a group of basinally-deposited, younger  archaean to proterozoic iron-rich sedimentary rocks now exposed in the hamersley range   <cit> .
much younger terrains underlain by unmetamorphosed sedimentary rock types surround the pilbara craton.
the present pilbara land surfaces are characterized by a long history of tectonic stability and extremely slow rates of erosion  <cit> .
unlike the east coast of australia where ongoing uplift commenced  <dig> million years ago   <cit> , the pilbara has not been tectonically active since the late jurassic  when a continental sliver rifted away from the pilbara region of northern australia  <cit> .
the pilbara craton displays remarkable topographic heterogeneity compared to the coastal plains of the southern carnarvon and valley plains of the ashburton regions, or the undulating dunes of the sandy deserts to the east.
major physiographic features of the pilbara include the plateaus and ranges of the hamersley and chichester, and the heavily excised gorges of the fortescue river as well as a number of other major river systems that traverse the pilbara including the ashburton, the de grey and the shaw.
only two previous studies, both on invertebrates, have evaluated phylogenetic structure of organisms that span the pilbara and surrounding regions  <cit> .
both have demonstrated that geological structures and complexity are fundamental predictors of phylogenetic structure for subterranean fauna, but there have been no similarly detailed studies of vertebrate animals that span the pilbara.
lucasium stenodactylum is a diplodactylid gecko with a large geographic range that includes the pilbara region, and is part of a gekkotan radiation of gondwanan origin  <cit> .
we have detailed elsewhere the broader-level molecular systematics and biogeography of the species-group to which l. stenodactylum belongs, and identified an ancient split between populations distributed in the broader pilbara region and those outside the pilbara  <cit> .
we are currently assessing the taxonomic status of the pilbara populations.
here we consider l. stenodactylum populations within the pilbara based on much larger and more detailed sampling, with particular reference to testing the hypothesis that ancient surface geology, and the habitats they determine, predicts phylogenetic structure in these rock dependent animals.
methods
tissue samples were obtained for  <dig> lucasium stenodactylum, of which  <dig> were from the pilbara region .
these animals represent nearly all the l. stenodactylum specimens in tissue collections held by the western australian museum , and given the current availability of specimens, exploit the maximum possible geographical sampling from the species' range.
samples from outside the pilbara are limited in the central and northern inland parts of australia and are considered elsewhere  <cit> .
for all individuals a 1200-base pair  region of the mitochondrial genome was sequenced that included the entire mitochondrial nadh dehydrogenase subunit  <dig>  gene and the flanking transfer rna  genes trnamet , and trnatrp , trnaala  and trnaasn .
details of primers and lab protocols are outlined elsewhere  <cit> .
lucasium squarrosum and l. wombeyi were used as outgroups for all analyses  <cit> .
therefore, we used unweighted parsimony and bayesian approaches to analyze the data.
heuristic parsimony analyses were implemented with the computer program paup* <dig> b <dig>  we used tbr branch swapping and ran the parsimony analysis five times from random starting points to make sure overall tree space was well searched.
bayesian analyses were run with the computer program mrbayes v <dig> .
<dig>  we ran analyses on the full data set as one unit and also partitioned the data into four components: non-coding trna and first, second and third codons.
for both sets of analyses we allowed all parameters to be estimated from the data during the runs and used the default value of four markov chains per run.
we ran the full analyses five times to make sure overall tree-space was well sampled and to avoid getting trapped in local optima.
we ran each analysis for a total of  <dig> , <dig> generations and sampled the chain every  <dig> generations, resulting in  <dig>  sampled trees.
log-likelihood values reached a plateau after approximately  <dig>  generations for both sets of analyses , so we discarded the first  <dig>  trees as the burn-in and used the last  <dig>  trees to estimate bayesian posterior probabilities.
we used the results of  <dig>  "fast" unweighted non-parametric parsimony bootstrap replicates and bayesian posterior probabilities to assess branch support.
we also used two complementary approaches to examine fine-scale genetic patterns across and within the pilbara clades using genalex  <dig>   <cit> .
we used mantel tests of matrix correspondence following the methods of smouse et al.
<cit>  to test for isolation by distance by comparing the pairwise geographic and total genetic distance matrices across all samples, and within clades, with tests of significance by  <dig>  permutations.
we also used a powerful micro-spatial autocorrelation technique developed by smouse and peakall  <cit>  that uses pairwise geographic and genetic distance matrices to dissect out positive genetic correlation at various geographic size classes .
a significant departure from the null hypothesis of no structure  is obtained when a positive r value falls outside of the 95% confidence interval and this is interpreted as the size of the "genetic neighbourhood" based on the data at hand.
we used the combined spatial autocorrelation analysis across the five clades to provide an overall estimate of genetic neighbourhood size.
we used a distance class of  <dig> km and evaluated  <dig> size classes, giving an evaluation over  <dig> km of space.
for all  <dig> individuals from the pilbara region, we assembled a distributional data set based on latitude and longitude for mapping with arcgis desktop .
shapefiles were created for each identified pilbara clade from our phylogeny and viewed in arcmap.
distributions were then overlain on to the geological regions dataset obtained from geoscience australia  <cit> .
this dataset was derived from the  <dig> bmr 1: <dig> , <dig> geological map of australia and includes metadata such as the age, landform and stratigraphy of different geographical regions within the australian continent.
the distributions of the major pilbara clades also were examined with respect to the land types on which they occur, based on a recent inventory and condition survey of the pilbara that takes into account landforms, soils, vegetation and drainage patterns  <cit> .
these terrain types relate fundamentally to the basement geology.
RESULTS
after the exclusion of unalignable regions in the trnas, the edited alignment comprised  <dig> characters and of these  <dig>  were variable and parsimony informative.
we used modeltest  <cit>  to compare model parameter estimates with those generated for the bayesian analyses in mr. bayes.
both the hierarchical likelihood ratio tests and the akaike information criterion in modeltest supported trn+i+g model as the best-fit substitution model for the combined data set and gave a -lnl =  <dig> .
the estimated parameters were as follows: nucleotide frequencies a =  <dig> , c =  <dig> , g =  <dig> , t =  <dig> ; substitution rates a↔c  <dig> , a↔g  <dig> , a↔t  <dig> , c↔g  <dig> , c↔t  <dig> , g↔t  <dig> ; proportion of invariant sites =  <dig> ; gamma shape parameter =  <dig> .
<dig>  the bayesian analyses produced parameter estimates that were very similar to those produced by modeltest .
unweighted parsimony analyses and both unpartitioned and partitioned bayesian analyses produced very similar topologies.
figure 2a shows one of  <dig> most parsimonious trees saved during the unweighted parsimony analysis.
parsimony bootstrap support values and unpartitioned bayesian posterior probabilities are shown for major branches only, but the results of our parsimony bootstrap and unpartitioned bayesian analyses with details of each individual sample are shown in additional file 2pepperetalbootstrap and additional file 3pepperetalbayesian .
the nd2+trna data set shows strong support for the deep and old phylogenetic split between populations in the pilbara region  and non-pilbara populations , in agreement with what we found based on a reduced number of samples, but additional genes  <cit> .
this result is consistent with the stark geological contrast between the ancient and exposed surface landscape and topography in the pilbara craton relative to the surrounding dunefields of the desert regions.
for example, the archaean to proterozoic igneous  and sedimentary rocks  form rugged outcrops at the surface in the pilbara, and are elevated compared with surrounding regions.
the lower-lying area to the south and east of the pilbara craton is generally overlain by much younger sedimentary rocks of proterozoic to cainozoic age that have given rise to sandplains and dunefields  <cit> .
the nd2+trna data set also identified five major clades within the pilbara endemic l. stenodactylum clade a, with generally high bootstrap values and posterior probabilities .
uncorrected genetic divergence between clades ranged from  <dig> % between a <dig> and a <dig> to  <dig> % between a <dig> and a <dig>  and maximum uncorrected divergence within each clade were as follows, a1:  <dig> %, a2:  <dig> %, a3:  <dig> %, a4:  <dig> %, a5:  <dig> %.
our genetic data suggests that the present-day lineages within the pilbara are relatively young with inter-clade divergences dating to approximately  <dig> –  <dig>  mya and maximum intra-clade divergences dating to  <dig>  –  <dig>  mya using a standard but simplistic reptile  <dig> %/my mtdna clock  <cit> .
non-pilbara clade b exhibits deep genetic structure but this is not the focus of this paper .
our phylogeny demonstrates that clade distributions are remarkably correlated with geological  and topographic  domains at both a large geographic scale  and also at the level of fine intra-craton provinciality.
clades a <dig>  a <dig>  a <dig> and a <dig> are predominantly restricted to the rugged pilbara craton whereas clade a <dig> is mostly confined to the lower-lying non-craton, for example the carnarvon sand plains west of the craton.
however, it is notable while clade a <dig> extends eastwards of the  <dig> m contour from the sand plains, this clade does not extend onto the craton at an identical elevation ; the additional distributional controlling factor may be the contrast in exposed rock types and their associated habitats.
within the pilbara craton, clade a <dig> is restricted to the hamersley plateau and other upland regions in the archaean-proterozoic ironstone units in the south of the craton.
clade a <dig> is restricted to the sandy abydos plain developed upon the archean granite unit in the north of the craton and is bound to the east by the distinctive gorge fold ranges, a series of parallel, razor-back fold ridges, rising up to  <dig> m above the surrounding terrain  <cit>  .
clade a <dig> also is found in the far easternmost part of this granite unit, but appears to be associated with the ridges of the gorge fold ranges.
although clade a <dig> is distributed throughout both the archaean-proterozoic ironstone and the archean granite unit in the eastern portion of the craton, this clade is absent from the ranges and occupies the lowlands in the eastern part of the craton, extending its distribution up the fortescue valley.
the only clade found outside the pilbara craton is clade a <dig>  which is almost exclusively associated with much younger sandstone units that dominate the low-lying carnarvon sand plains along the northwestern australian coastline.
a small number of individuals from this clade are found within the pilbara craton on the burrup peninsula, but these individuals are geographically connected to the rest of the clade by the low-lying sandplains that stretch along the coast.
a number of major rivers traverse the pilbara, but the distribution of clades suggests these generally dry riverbeds are not current barriers to dispersal in these animals.
each clade exhibits a strong preference for a particular terrain type, each of which are related directly to underlying geological structure.
for example, clades a <dig> and a <dig> are associated with stony lower slopes and plains, while clade a <dig> is associated with sandplains, dunes and clay pans.
clade a <dig> also prefers stony lower slopes and plains while clade a <dig> is distributed between sandplains, dunes and claypans and stony lower slopes and plains.
a mantel test across all five clades showed evidence of isolation by distance in the genetic data .
strong isolation by distance also was detected in four of the five individual clades , the exception being clade a <dig>  which has too few samples to do a meaningful test.
the combined spatial autocorrelation analysis showed a significantly positive r value at all distance classes less than  <dig> km , which demonstrates that the spatial genetic structure varies on a regional scale that matches the size and distribution of available habitats.
CONCLUSIONS
together our phylogenetic, distributional, geological and habitat data provide a textbook example of ecological diversification across an ancient and heterogeneous landscape.
basement geology and superimposed topography have determined the varied micro-habitats available across the pilbara.
each of the different types of habitats occupied by these lizards will experience quite different temperature and humidity regimes, each of which should contribute strongly to distribution of clades.
while there are a variety of possible explanations, our favoured hypothesis is that ancestors of the pilbara lineages radiated into the region at the onset of aridity in australia approximately  <dig> mya and locally adapted to these ancient and highly stable terrain types and the micro-habitats derived from them.
our results are broadly consistent with patterns revealed by recent detailed phylogenetic work on subterranean diving beetles  <cit>  and amphipods  <cit>  in the pilbara region.
the pilbara is well known as important an important area for refugia with high levels of endemism in a variety of plants and animals  <cit> , but we are only just now beginning to reconstruct the biotic history of this ancient and complex landscape.
the geological history and the habitats and micro-environments derived from geological formations will form an important part of this emerging story.
authors' contributions
mp carried out all the molecular genetic work and assembled and analysed the distributional data.
pd provided access to tissue samples and participated in the design of the project.
ra provided expertise on the geology of the pilbara region.
jsk carried out the phylogenetic and spatial genetic analyses and participated in the design of the project.
mp and jsk wrote the manuscript.
all authors have read and approved the manuscript.
supplementary material
additional file 1
click here for file
additional file 2
pepperetalbootstrap
click here for file
additional file 3
pepperetalbayesian
click here for file
acknowledgements
this project was funded by australian research council grants to jsk.
we thank the calm pbs for access to tissues, christine hayes for help in the molecular lab, peter kendrick for help in the field, simon blomberg for statistical assistance and two anonymous reviewers for helpful comments that improved the manuscript.

 the pilbara in northwestern australia is an ancient geological and biophysical region that has repeatedly been identified as an important zone of biological diversity, endemism and past refugia based on congruent biogeographic patterns of both fauna and flora. our understanding of the biotic history of the region is still in its infancy, but the geology of the pilbara region is well understood and may serve as an important predictor of biological diversity. the pilbara is underlain by sedimentary and igneous rocks ranging up to billion years in age, and overlain by some of the most ancient erosion surfaces on earth. the pilbara craton in the northeastern portion of the region has two parts; archaean granites and metamorphosed volcanic rocks in the north, and, stratigraphically overlying these rocks in the south, is a group of basinally-deposited, younger archaean to proterozoic iron-rich sedimentary rocks now exposed in the hamersley range. much younger terrains underlain by unmetamorphosed sedimentary rock types surround the pilbara craton. the present pilbara land surfaces are characterized by a long history of tectonic stability and extremely slow rates of erosion. unlike the east coast of australia where ongoing uplift commenced million years ago, the pilbara has not been tectonically active since the late jurassic when a continental sliver rifted away from the pilbara region of northern australia. the pilbara craton displays remarkable topographic heterogeneity compared to the coastal plains of the southern carnarvon and valley plains of the ashburton regions, or the undulating dunes of the sandy deserts to the east. major physiographic features of the pilbara include the plateaus and ranges of the hamersley and chichester, and the heavily excised gorges of the fortescue river as well as a number of other major river systems that traverse the pilbara including the ashburton, the de grey and the shaw. only two previous studies, both on invertebrates, have evaluated phylogenetic structure of organisms that span the pilbara and surrounding regions. both have demonstrated that geological structures and complexity are fundamental predictors of phylogenetic structure for subterranean fauna, but there have been no similarly detailed studies of vertebrate animals that span the pilbara. lucasium stenodactylum is a diplodactylid gecko with a large geographic range that includes the pilbara region, and is part of a gekkotan radiation of gondwanan origin. we have detailed elsewhere the broader-level molecular systematics and biogeography of the species-group to which l. stenodactylum belongs, and identified an ancient split between populations distributed in the broader pilbara region and those outside the pilbara. we are currently assessing the taxonomic status of the pilbara populations. here we consider l. stenodactylum populations within the pilbara based on much larger and more detailed sampling, with particular reference to testing the hypothesis that ancient surface geology, and the habitats they determine, predicts phylogenetic structure in these rock dependent animals. methods tissue samples were obtained for lucasium stenodactylum, of which were from the pilbara region. these animals represent nearly all the l. stenodactylum specimens in tissue collections held by the western australian museum, and given the current availability of specimens, exploit the maximum possible geographical sampling from the species' range. samples from outside the pilbara are limited in the central and northern inland parts of australia and are considered elsewhere. for all individuals a 1200-base pair region of the mitochondrial genome was sequenced that included the entire mitochondrial nadh dehydrogenase subunit gene and the flanking transfer rna genes trnamet, and trnatrp, trnaala and trnaasn. details of primers and lab protocols are outlined elsewhere. lucasium squarrosum and l. wombeyi were used as outgroups for all analyses. therefore, we used unweighted parsimony and bayesian approaches to analyze the data. heuristic parsimony analyses were implemented with the computer program paup* b we used tbr branch swapping and ran the parsimony analysis five times from random starting points to make sure overall tree space was well searched. bayesian analyses were run with the computer program mrbayes v. we ran analyses on the full data set as one unit and also partitioned the data into four components: non-coding trna and first, second and third codons. for both sets of analyses we allowed all parameters to be estimated from the data during the runs and used the default value of four markov chains per run. we ran the full analyses five times to make sure overall tree-space was well sampled and to avoid getting trapped in local optima. we ran each analysis for a total of, generations and sampled the chain every generations, resulting in sampled trees. log-likelihood values reached a plateau after approximately generations for both sets of analyses, so we discarded the first trees as the burn-in and used the last trees to estimate bayesian posterior probabilities. we used the results of "fast" unweighted non-parametric parsimony bootstrap replicates and bayesian posterior probabilities to assess branch support. we also used two complementary approaches to examine fine-scale genetic patterns across and within the pilbara clades using genalex. we used mantel tests of matrix correspondence following the methods of smouse et al. to test for isolation by distance by comparing the pairwise geographic and total genetic distance matrices across all samples, and within clades, with tests of significance by permutations. we also used a powerful micro-spatial autocorrelation technique developed by smouse and peakall that uses pairwise geographic and genetic distance matrices to dissect out positive genetic correlation at various geographic size classes. a significant departure from the null hypothesis of no structure is obtained when a positive r value falls outside of the 95% confidence interval and this is interpreted as the size of the "genetic neighbourhood" based on the data at hand. we used the combined spatial autocorrelation analysis across the five clades to provide an overall estimate of genetic neighbourhood size. we used a distance class of km and evaluated size classes, giving an evaluation over km of space. for all individuals from the pilbara region, we assembled a distributional data set based on latitude and longitude for mapping with arcgis desktop. shapefiles were created for each identified pilbara clade from our phylogeny and viewed in arcmap. distributions were then overlain on to the geological regions dataset obtained from geoscience australia. this dataset was derived from the bmr 1:, geological map of australia and includes metadata such as the age, landform and stratigraphy of different geographical regions within the australian continent. the distributions of the major pilbara clades also were examined with respect to the land types on which they occur, based on a recent inventory and condition survey of the pilbara that takes into account landforms, soils, vegetation and drainage patterns. these terrain types relate fundamentally to the basement geology.  after the exclusion of unalignable regions in the trnas, the edited alignment comprised characters and of these were variable and parsimony informative. we used modeltest to compare model parameter estimates with those generated for the bayesian analyses in mr. bayes. both the hierarchical likelihood ratio tests and the akaike information criterion in modeltest supported trn+i+g model as the best-fit substitution model for the combined data set and gave a -lnl =. the estimated parameters were as follows: nucleotide frequencies a =, c =, g =, t =; substitution rates a↔c, a↔g, a↔t, c↔g, c↔t, g↔t; proportion of invariant sites =; gamma shape parameter =. the bayesian analyses produced parameter estimates that were very similar to those produced by modeltest. unweighted parsimony analyses and both unpartitioned and partitioned bayesian analyses produced very similar topologies. figure 2a shows one of most parsimonious trees saved during the unweighted parsimony analysis. parsimony bootstrap support values and unpartitioned bayesian posterior probabilities are shown for major branches only, but the results of our parsimony bootstrap and unpartitioned bayesian analyses with details of each individual sample are shown in additional file 2pepperetalbootstrap and additional file 3pepperetalbayesian. the nd2+trna data set shows strong support for the deep and old phylogenetic split between populations in the pilbara region and non-pilbara populations, in agreement with what we found based on a reduced number of samples, but additional genes. this result is consistent with the stark geological contrast between the ancient and exposed surface landscape and topography in the pilbara craton relative to the surrounding dunefields of the desert regions. for example, the archaean to proterozoic igneous and sedimentary rocks form rugged outcrops at the surface in the pilbara, and are elevated compared with surrounding regions. the lower-lying area to the south and east of the pilbara craton is generally overlain by much younger sedimentary rocks of proterozoic to cainozoic age that have given rise to sandplains and dunefields. the nd2+trna data set also identified five major clades within the pilbara endemic l. stenodactylum clade a, with generally high bootstrap values and posterior probabilities. uncorrected genetic divergence between clades ranged from % between a and a to % between a and a and maximum uncorrected divergence within each clade were as follows, a1: %, a2: %, a3: %, a4: %, a5: %. our genetic data suggests that the present-day lineages within the pilbara are relatively young with inter-clade divergences dating to approximately – mya and maximum intra-clade divergences dating to – mya using a standard but simplistic reptile %/my mtdna clock. non-pilbara clade b exhibits deep genetic structure but this is not the focus of this paper. our phylogeny demonstrates that clade distributions are remarkably correlated with geological and topographic domains at both a large geographic scale and also at the level of fine intra-craton provinciality. clades a a a and a are predominantly restricted to the rugged pilbara craton whereas clade a is mostly confined to the lower-lying non-craton, for example the carnarvon sand plains west of the craton. however, it is notable while clade a extends eastwards of the m contour from the sand plains, this clade does not extend onto the craton at an identical elevation; the additional distributional controlling factor may be the contrast in exposed rock types and their associated habitats. within the pilbara craton, clade a is restricted to the hamersley plateau and other upland regions in the archaean-proterozoic ironstone units in the south of the craton. clade a is restricted to the sandy abydos plain developed upon the archean granite unit in the north of the craton and is bound to the east by the distinctive gorge fold ranges, a series of parallel, razor-back fold ridges, rising up to m above the surrounding terrain. clade a also is found in the far easternmost part of this granite unit, but appears to be associated with the ridges of the gorge fold ranges. although clade a is distributed throughout both the archaean-proterozoic ironstone and the archean granite unit in the eastern portion of the craton, this clade is absent from the ranges and occupies the lowlands in the eastern part of the craton, extending its distribution up the fortescue valley. the only clade found outside the pilbara craton is clade a which is almost exclusively associated with much younger sandstone units that dominate the low-lying carnarvon sand plains along the northwestern australian coastline. a small number of individuals from this clade are found within the pilbara craton on the burrup peninsula, but these individuals are geographically connected to the rest of the clade by the low-lying sandplains that stretch along the coast. a number of major rivers traverse the pilbara, but the distribution of clades suggests these generally dry riverbeds are not current barriers to dispersal in these animals. each clade exhibits a strong preference for a particular terrain type, each of which are related directly to underlying geological structure. for example, clades a and a are associated with stony lower slopes and plains, while clade a is associated with sandplains, dunes and clay pans. clade a also prefers stony lower slopes and plains while clade a is distributed between sandplains, dunes and claypans and stony lower slopes and plains. a mantel test across all five clades showed evidence of isolation by distance in the genetic data. strong isolation by distance also was detected in four of the five individual clades, the exception being clade a which has too few samples to do a meaningful test. the combined spatial autocorrelation analysis showed a significantly positive r value at all distance classes less than km, which demonstrates that the spatial genetic structure varies on a regional scale that matches the size and distribution of available habitats.  together our phylogenetic, distributional, geological and habitat data provide a textbook example of ecological diversification across an ancient and heterogeneous landscape. basement geology and superimposed topography have determined the varied micro-habitats available across the pilbara. each of the different types of habitats occupied by these lizards will experience quite different temperature and humidity regimes, each of which should contribute strongly to distribution of clades. while there are a variety of possible explanations, our favoured hypothesis is that ancestors of the pilbara lineages radiated into the region at the onset of aridity in australia approximately mya and locally adapted to these ancient and highly stable terrain types and the micro-habitats derived from them. our results are broadly consistent with patterns revealed by recent detailed phylogenetic work on subterranean diving beetles and amphipods in the pilbara region. the pilbara is well known as important an important area for refugia with high levels of endemism in a variety of plants and animals, but we are only just now beginning to reconstruct the biotic history of this ancient and complex landscape. the geological history and the habitats and micro-environments derived from geological formations will form an important part of this emerging story. authors' contributions mp carried out all the molecular genetic work and assembled and analysed the distributional data. pd provided access to tissue samples and participated in the design of the project. ra provided expertise on the geology of the pilbara region. jsk carried out the phylogenetic and spatial genetic analyses and participated in the design of the project. mp and jsk wrote the manuscript. all authors have read and approved the manuscript. supplementary material additional file 1 click here for file additional file 2 pepperetalbootstrap click here for file additional file 3 pepperetalbayesian click here for file acknowledgements this project was funded by australian research council grants to jsk. we thank the calm pbs for access to tissues, christine hayes for help in the molecular lab, peter kendrick for help in the field, simon blomberg for statistical assistance and two anonymous reviewers for helpful comments that improved the manuscript.
our phylogeny provides strong support for a deep and ancient phylogenetic split within l. stenodactylum that distinguishes populations within the pilbara region from those outside the pilbara. this is the first detailed phylogenetic examination of a vertebrate organism across the pilbara region. within the pilbara region itself, our phylogeny has identified five major clades whose distribution closely matches different surface geologies of this ancient landscape. the iconic pilbara in northwestern australia is an ancient geological and biophysical region that is an important zone of biodiversity, endemism and refugia. together our phylogenetic, distributional, geological and habitat data provide a clear example of ecological diversification across an ancient and heterogeneous landscape. our results show the geological history and the habitats derived from them will form an important part of this emerging story. our favoured hypothesis is that ancestors of the pilbara lineages radiated into the region at the onset of aridity in australia approximately mya and locally adapted to the various ancient and highly stable terrain types and the micro-habitats derived from them. in terms of specimen recovery and analysis, we are only beginning to reconstruct the biotic history of this ancient landscape. we combined phylogenetic and spatial-autocorrelation genetic analyses of mitochondrial dna data on populations of the gekkotan lizard lucasium stenodactylum within the pilbara with geological, distributional and habitat data to test the hypothesis that ancient surface geology predicts current clade-habitat associations in saxicoline animals.
 the iconic pilbara in northwestern australia is an ancient geological and biophysical region that is an important zone of biodiversity, endemism and refugia. it also is overlain by some of the oldest erosion surfaces on earth, but very little is known about the patterns of biotic diversity within the pilbara or how they relate to the landscape. we combined phylogenetic and spatial-autocorrelation genetic analyses of mitochondrial dna data on populations of the gekkotan lizard lucasium stenodactylum within the pilbara with geological, distributional and habitat data to test the hypothesis that ancient surface geology predicts current clade-habitat associations in saxicoline animals.  this is the first detailed phylogenetic examination of a vertebrate organism across the pilbara region. our phylogeny provides strong support for a deep and ancient phylogenetic split within l. stenodactylum that distinguishes populations within the pilbara region from those outside the pilbara. within the pilbara region itself, our phylogeny has identified five major clades whose distribution closely matches different surface geologies of this ancient landscape. each clade shows strong affinities with particular terrain types and topographic regions, which are directly related to different geological bedrock.  together our phylogenetic, distributional, geological and habitat data provide a clear example of ecological diversification across an ancient and heterogeneous landscape. our favoured hypothesis is that ancestors of the pilbara lineages radiated into the region at the onset of aridity in australia approximately mya and locally adapted to the various ancient and highly stable terrain types and the micro-habitats derived from them. in terms of specimen recovery and analysis, we are only beginning to reconstruct the biotic history of this ancient landscape. our results show the geological history and the habitats derived from them will form an important part of this emerging story. 