PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcebBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Evolutionary Biology
 
BMC Evol Biol. 2011; 11: 191.
Published online Jul 4, 2011. doi:  10.1186/1471-2148-11-191
PMCID: PMC3141439
Phylogeographic divergence in the widespread delicate skink (Lampropholis delicata) corresponds to dry habitat barriers in eastern Australia
David G Chapple,corresponding author1,2,3 Conrad J Hoskin,4,5 Stephanie NJ Chapple,2,6 and Michael B Thompson7
1School of Biological Sciences, Monash University, Clayton, Victoria 3800, Australia
2Museum Victoria, Division of Sciences, GPO Box 666, Melbourne, Victoria 3001, Australia
3Allan Wilson Centre for Molecular Ecology and Evolution, School of Biological Sciences, Victoria University of Wellington, PO Box 600, Wellington 6140, New Zealand
4Division of Evolution, Ecology & Genetics, Research School of Biology, Australian National University, Canberra, Australian Capital Territory 0200, Australia
5School of Marine and Tropical Biology, James Cook University, Townsville, Queensland 4811, Australia
6Department of Zoology, University of Melbourne, Melbourne, Victoria 3010, Australia
7School of Biological Sciences, University of Sydney, The Heydon-Laurence Building A08, New South Wales 2006, Australia
corresponding authorCorresponding author.
David G Chapple: David.Chapple/at/monash.edu; Conrad J Hoskin: Conrad.Hoskin/at/anu.edu.au; Stephanie NJ Chapple: schapple/at/museum.vic.gov.au; Michael B Thompson: Mike.Thompson/at/sydney.edu.au
Received March 14, 2011; Accepted July 4, 2011.
Background
The mesic habitats of eastern Australia harbour a highly diverse fauna. We examined the impact of climatic oscillations and recognised biogeographic barriers on the evolutionary history of the delicate skink (Lampropholis delicata), a species that occurs in moist habitats throughout eastern Australia. The delicate skink is a common and widespread species whose distribution spans 26° of latitude and nine major biogeographic barriers in eastern Australia. Sequence data were obtained from four mitochondrial genes (ND2, ND4, 12SrRNA, 16SrRNA) for 238 individuals from 120 populations across the entire native distribution of the species. The evolutionary history and diversification of the delicate skink was investigated using a range of phylogenetic (Maximum Likelihood, Bayesian) and phylogeographic analyses (genetic diversity, ΦST, AMOVA, Tajima's D, Fu's F statistic).
Results
Nine geographically structured, genetically divergent clades were identified within the delicate skink. The main clades diverged during the late Miocene-Pliocene, coinciding with the decline and fragmentation of rainforest and other wet forest habitats in eastern Australia. Most of the phylogeographic breaks within the delicate skink were concordant with dry habitat or high elevation barriers, including several recognised biogeographic barriers in eastern Australia (Burdekin Gap, St Lawrence Gap, McPherson Range, Hunter Valley, southern New South Wales). Genetically divergent populations were also located in high elevation topographic isolates inland from the main range of L. delicata (Kroombit Tops, Blackdown Tablelands, Coolah Tops). The species colonised South Australia from southern New South Wales via an inland route, possibly along the Murray River system. There is evidence for recent expansion of the species range across eastern Victoria and into Tasmania, via the Bassian Isthmus, during the late Pleistocene.
Conclusions
The delicate skink is a single widespread, but genetically variable, species. This study provides the first detailed phylogeographic investigation of a widespread species whose distribution spans virtually all of the major biogeographic barriers in eastern Australia.
The coastal regions of eastern Australia are currently dominated by wet forest and drier sclerophyllous habitats that harbour a highly diverse fauna [1,2]. While the majority (~70%) of the Australian continent is covered by arid or semi-arid vegetation, eastern Australia provides a narrow, but largely continuous expanse of habitat for mesic-adapted species [1,3,4]. These mesic habitats are generated through the presence of the Great Dividing Range (GDR), which abuts the entire length of the east coast (~2,500 km) in a north-south alignment ([5-7]; Figure Figure1).1). In the context of an expansive continent that is characterised by low topographic relief, the moderate elevation (~1000-1300 m, maximum ~2300 m) provided by the GDR generates altitudinal, climatic and environmental variation, and precipitates the required moisture to support mesic vegetation [3,5,7].
Figure 1
Figure 1
The major biogeographic barriers in eastern Australia. A description of each barrier is provided in Table 1. The location of Kroombit Tops is shown in Figure 2. Inset: The coastline of the Bass Strait region 14 kya. Tasmania has repeatedly been connected (more ...)
Although widespread glaciation never occurred in Australia [8,9], climatic oscillations have driven repeated altitudinal and distributional shifts in mesic habitats along the eastern margin of the continent (reviewed in [4]). Palaeoclimatic studies indicate that the extent and composition of the vegetation has fluctuated dramatically over the last 10 myr, although there has been a general transition from rainforest towards drier environments and sclerophyllous vegetation [1,3,8]. The rainforests that had previously dominated eastern Australia contracted between the mid- and late-Miocene, giving way to woodland and open forest vegetation that was more suited to the drier climates [3,10-12]. Lowered sea level associated with globally drier conditions facilitated the expansion of vegetation into the low lying regions of south-eastern Australia (e.g. Gippsland Basin, Murray Basin; Figure Figure1)1) that had previously been subject to marine inundation [3,11,12]. Although the extent of rainforests briefly expanded again during the early Pliocene due to a temporary return to warm and wet conditions, by the end of the Pliocene open woodlands, sclerophyllous forests and grasslands dominated the landscape of eastern Australia [3,6,8,10].
The cool-dry to warm-wet climatic fluctuations that commenced during the Pliocene intensified throughout the Pleistocene and led to the repeated expansion and contraction of mesic habitats in eastern Australia and the regular encroachment of drier habitats into the coastal fringes [4,8,13,14]. There was periodic flooding of the low lying coastal and inland basins in eastern Australia during the sea level changes associated with these climatic cycles [3,6,15], which also resulted in the connection of Tasmania (TAS) to the mainland during glacial periods by Bass Strait land bridges [16] (Figure (Figure1).1). At present, the once widespread rainforest and wet forest vegetation is restricted to small, scattered remnants within a mosaic of dry sclerophyll woodlands and open forests along the east coast [1,17].
The evolutionary history of the resident fauna of the narrow mesic strip along the east coast has been influenced by both habitat barriers and physical barriers (e.g. mountain ranges, sea straits), which led to genetic divergence and, in some cases, speciation of allopatric populations [1,18]. The most well-studied barrier in eastern Australia has been the Black Mountain Corridor (BMC) in the Wet Tropics region of north Queensland (QLD). This thin strip of rainforest currently connects the northern and southern rainforest block of the Wet Tropics but was repeatedly severed in the past by dry forest habitats during globally drier climates [18,19]. Intensive research has revealed largely concordant patterns of genetic divergence across the barrier in a wide range of rainforest taxa [e.g. [19-24]], and improved our understanding of how these barriers, in concert with climatic oscillations, have generated the high levels of biodiversity evident in eastern Australia [18,25]. However, at least nine other biogeographic barriers have been identified in eastern Australia (Tables (Tables11 and and2;2; Figure Figure1),1), several of which have yet to be investigated in detail. These include dry habitat barriers (Burdekin Gap, St Lawrence Gap, Hunter Valley), mountain ranges that act as topographic barriers (McPherson Range, southern New South Wales [NSW]), disjunct inland mountains (Kroombit Tops), sea straits (Bass Strait), and marine basins (Gippsland Basin, Murray Basin) (Tables (Tables11 and and2;2; Figure Figure11).
Table 1
Table 1
Description of the recognised biogeographic barriers in eastern Australia (see Figure 1).
Table 2
Table 2
The impact of the recognised biogeographic barriers in eastern Australia (Figure 1, Table 1) on vertebrates, invertebrates and plants.
A recent study investigated the impact of five biogeographic barriers in south-eastern Australia [26]; however, here we adopt a broader approach and examine the influence of nine biogeographic barriers (Tables (Tables11 and and2)2) throughout eastern Australia on the evolutionary history of the resident biota. In particular, we focus on the delicate skink, Lampropholis delicata (De Vis, 1888 [27]), which is unusual in that its distribution is so broad that it spans all of these barriers in eastern Australia (Figure (Figure2,2, ,3).3). The delicate skink is a small lizard (adult snout-vent length 35-51 mm) whose distribution extends across 26° of latitude from Cairns in north QLD to Hobart in TAS, with disjunct populations in far western Victoria (VIC) and south-eastern South Australia (SA) ([28]; Figure Figure2,2, ,3).3). It is a common species that occurs across a range of moist habitats, including rainforests, wet sclerophyll forests, woodland and heaths [28]. However, it also thrives in disturbed habitats and is one the most common skink species in suburban gardens along the east coast [28].
Figure 2
Figure 2
Map of eastern Australia indicating the collection localities of the Lampropholis delicata samples. The distribution of the nine major clades are indicated: clade 1 (grey solid circles), clade 2 (hollow triangles), clade 3 (black solid squares), clade (more ...)
Figure 3
Figure 3
The distribution of Lampropholis delicata clades and subclades. The distribution of each subclade identified in Figure 4 and 5 is presented (solid lines), along with the location of recognised biogeographic barriers in eastern Australia (dashed lines; (more ...)
Here we examine the phylogeography of the delicate skink using 2426 bp of mitochondrial DNA sequence data (ND2, ND4, 12SrRNA, 16SrRNA) from across the entire native range of the species (Figure (Figure2,2, ,3).3). Due to its presence in TAS and eastern VIC, Rawlinson [29] suggested that the delicate skink was a glacial relic that had occurred in southern Australia for a prolonged period of time. However, it has been implied that the delicate skink might not be native to TAS, instead reaching the state via human-assisted colonisation. This is because the delicate skink was not detected in TAS until 1963, although subsequent examination of museum collections revealed that previously mis-identified specimens had been collected during the 1920s and 1930s [30]. Unlike many reptile species whose distribution spans Bass Strait, the delicate skink does not occur on Wilsons Promontory (the most southerly projection of the Australian mainland), but it does occur on a Bass Strait Island (Flinders Island) that formed part of the Bassian Isthmus during the last glacial maxima ([28]; Figure Figure1,1, ,2).2). We conduct a range of phylogeographic analyses to examine Rawlinson's [29] hypothesis, determine the status of the Tasmanian population, and investigate the impact of historical processes on the evolutionary history of the delicate skink.
Sampling
We obtained tissue samples from 238 Lampropholis delicata, representing 120 different populations, from across the entire Australian range of the species (Figure (Figure2,2, ,3;3; Additional files 1, 2). Samples were obtained from the frozen-tissue collections of several Australian Museums (Australian Museum, South Australian Museum, CSIRO Australian National Wildlife Collection), along with our own field collections (Additional files 1, 2). We included the closely related L. guichenoti (Australian Museum NR2639) and an Australian Eugongylus-lineage skink Niveoscincus pretiosus (Australian Museum NR391) as outgroups in our study.
DNA extraction, amplification and sequencing
Total genomic DNA was extracted from liver, muscle, toe or tail-tip samples using a Qiagen DNeasy Blood and Tissue Extraction Kit (Qiagen, Hilden, Germany). For each sample we sequenced portions of four mitochondrial genes: ND2 (~600 bp), ND4 (~700 bp), 12SrRNA (~700 bp), and 16SrRNA (~500 bp). These regions were targeted because work across several taxonomic levels in squamate reptiles has indicated useful levels of variability [e.g. [26,31-34]]. The primers used to amplify and sequence these regions are provided in Additional file 3. PCR was conducted as outlined in Greaves et al. [35], except on a Corbett Research GC1-960 thermal cycler. PCR products were purified using ExoSAP-IT (USB Corporation, Cleveland, Ohio USA). The purified product was sequenced directly using a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) and then analysed on an ABI 3730XL capillary sequencer.
Sequence data were edited using CONTIGEXPRESS in VECTOR NTI ADVANCE v9.1.0 (Invitrogen), and aligned using the default parameters of CLUSTAL X v1.83[36]. We translated all coding region sequences to confirm that none contained premature stop codons. Sequence data were submitted to GenBank [GenBank:JF438009-JF438959, EF567304, EU567726, EU567768, EU567769, EU567927, EU567928, EU568019, EU568020] (Additional file 1).
Phylogenetic analyses
Maximum Likelihood (ML) and Bayesian tree building methods were used. We used MODELTEST 3.7 [37] to identify the most appropriate model of sequence evolution based on the AIC criterion. MODELTEST, conducted in PAUP* 4.0b10 [38], was also used to estimate base frequencies, substitution rates, the proportion of invariable sites (I) and the among-site substitution rate variation (G). These values were then used as settings in PhyML 3.0 [39] to generate a ML tree with 500 bootstraps.
MRBAYES 3.1.2 [40] was then used to complete Bayesian analyses. Preliminary analysis of each mtDNA region revealed congruent tree topologies. In order to evaluate partitioning strategies, we used MODELTEST to determine the most appropriate model for each partition. We then conducted a Bayesian analysis for each partitioning strategy, applying the appropriate model of evolution to each partition, and allowing among-partition rate variation. We ran each Bayesian analysis for five million generations, sampling every 100 generations (i.e. 50,000 sampled trees). We ran each analysis twice, using four heated chains per run. We discarded the first 25% of samples as burn-in and the last 37,500 trees were used to estimate the Bayesian posterior probabilities. In order to calculate the AIC and BIC scores for different partitions strategies, we calculated the number of parameters for each. Following McGuire et al. [41], for each parameter we added the number of substitution rates for the model suggested by MODELTEST for that partition (6 for GTR, or 2 for HKY), the number of free equilibrium base frequencies (3 for GTR and HKY), plus one parameter per partition where appropriate for each of I and/or G. For multi-partition strategies, we also added one parameter per partition, corresponding to the among-partition rate multiplier. To calculate the AIC and BIC scores, we used the equations: AIC = -2Li + 2ki and BIC = -2Li + (ki).(ln n) (where Li is the harmonic mean log likelihood for partitioning strategy i, ki is the total number of parameters for partitioning strategy i, and n is the total number of nucleotides). The program TRACER 1.5 [42] was used to check for chain convergence and mixing. Specifically, raw traces of sampled values versus MCMC step numbers were examined to confirm that there was no trend away from the mean and that there were no large fluctuations in the likelihood values.
Bootstrap values (500 ML bootstraps) and Bayesian posterior probabilities were used to assess branch support. We considered branches supported by bootstrap values of 70% or greater [43], and/or posterior probability values greater than or equal to 95% [44] to be supported by our data.
Molecular diversity and population divergence
Estimates of genetic diversity within L. delicata clades (number of haplotypes, h; haplotypic diversity, Hd; number of polymorphic sites, S; nucleotide diversity, π) were calculated in DNASP v4.50 [45]. Tamura-Nei (TrN)-corrected genetic distances within and among clades were calculated in MEGA 4 [46]. Genetic differentiation among clades within L. delicata was estimated in ARLEQUIN v3.5 [47]. Pairwise ΦST values (an analogue of Wright's fixation index FST) were calculated to estimate among clade differentiation. We conducted hierarchical Analysis of Molecular Variance (AMOVA; [48]) to investigate the impact of the a priori (Tables (Tables11 and and2)2) and a posteriori biogeographic barriers on the partitioning of genetic variation within L. delicata. Both tests used TrN genetic distances with gamma correction (using the value calculated from MODELTEST). Significance levels of all the estimated values were calculated by 10,000 permutations, and adjusted according to the Bonferroni correction procedure [49] for multiple pairwise comparisons as described by Holm [50].
We used Tajima's D [50], Fu's F statistic [51] (calculated in ARLEQUIN) and mismatch distributions to test for signatures of population expansion within L. delicata clades. Significant and negative Tajima's D and Fu's F statistic values are indicative of possible population expansion. Mismatch frequency histograms were plotted in DNASP to determine whether the clades exhibited evidence of spatial range expansion or a stationary population history [52]. A smooth bell shape signifies either population expansion or spatial range expansion, whereas a multimodal distribution represents a long history in situ [53-56]. To distinguish between these two types of distribution, a raggedness index (RI, sum of the squared difference between neighbouring peaks) and the sum of squared deviations (SSD) between the observed and expected mismatch were calculated using the methods of Schneider & Excoffier [57] in ARLEQUIN. The spatial expansion hypothesis (both RI and SSD) was tested using a parametric bootstrap approach (200 replicates).
As there are no suitable fossil calibration points available for Lampropholis skinks, we estimated the divergence time of L. delicata clades using an evolutionary rate of 1.3-1.63% sequence divergence per million years, based on mitochondrial DNA calibrations from other squamate reptile groups (1.3%, [58]; 1.42-1.63%, [59]; 1.55%, [60]; 1.62%, [61]; 1.63%, [62]). A strict molecular clock (0.0065-0.00815 per lineage substitution rate), implemented in BEAST v1.6.1 [63], was used to estimate the divergence times within L. delicata. The Australian Eugongylus lineage is estimated to have originated ~20 mya [64], and this information was used as the maximum age of the tree root. A GTR+I+G model of evolution was employed with a coalescent (Bayesian skyline) tree prior. The analysis was run twice, with 20 million generations per run (total 40 million generations). The output was viewed in TRACER to check that stationarity had been reached, and ensure that the effective sample size (ESS) exceeded 200 [63]. The two separate runs were then combined using LOGCOMBINER v1.6.1, with a maximum clade credibility tree generated in TREEANNOTATOR v1.6.1 and visualised in FIGTREE v1.3.1. A Bayesian skyline plot [65] was also generated in TRACER to examine the magnitude and timing of population size changes in L. delicata.
Molecular diversity and phylogeographic structure
The edited alignment comprised 2426 characters (550 bp ND2, 671 bp ND4, 708 bp 12SrRNA, 497 bp 16SrRNA; Additional file 4), of which 813 (33.5%) were variable and 587 (24.2%) were parsimony-informative. For the ingroup only, the alignment contained 638 (26.3%) variable characters, of which 543 (22.4%) were parsimony-informative. Base frequencies were unequal (A = 0.3694, T = 0.2172, C = 0.2889, G = 0.1245), but a χ2 test confirmed the homogeneity of base frequencies among sequences (df = 498, P = 1.0). The phylogenetic analyses were conducted on a dataset comprising the 165 unique haplotypes that were present within L. delicata (Table (Table33).
Table 3
Table 3
Estimates of genetic diversity for the clades present within Lampropholis delicata.
The AIC from MODELTEST supported the GTR + I + G substitution model as the most appropriate for our unpartitioned dataset. Parameters estimated under this model were: relative substitution rates (A↔C = 2.1053, A↔G = 44.1898, A↔T = 2.5350, C↔G = 1.1020, C↔T = 29.2135, relative to G↔T = 1.00), proportion of invariable sites (0.5241), and gamma distribution shape parameter (0.7015). We evaluated three partitioning strategies for our dataset (Table (Table4).4). The unpartitioned and by gene partitioning strategy were analysed using the GTR + I + G model for all nucletotides. When the data was partitioned by codon, we used a mixture of GTR and HKY models for each partition (Table (Table4).4). Both the AIC and BIC scores ranked the most highly parameterised strategy (by gene and codon) as the most appropriate. However, the topologies of the ML, unpartitioned Bayesian and partitioned Bayesian trees were congruent, therefore we present the optimal ML tree (-ln L = 14501.43296) with ML bootstrap (BS) values and unpartitioned Bayesian posterior probabilities (PP) indicating branch support (Figure (Figure4,4, ,5).5). Nine well-supported, non-overlapping clades (labelled from the top to bottom of the tree, and roughly related to their north to south geographic distribution) are present within L. delicata (Figure (Figure2,2, ,3,3, ,4,4, ,5),5), with high levels of haplotypic and nucleotide diversity within each clade (Table (Table3).3). The PP's of the main clades and subclades in the partitioned Bayesian analysis were identical to those from the unpartitioned Bayesian analysis presented in Figure Figure44 and and5,5, except that the support value for Clade 9 was lower (0.85 rather than 0.95).
Table 4
Table 4
Test of alternative partitioning strategies.
Figure 4
Figure 4
Phylogram for the delicate skink (Lampropholis delicata). The phylogeny is based on 2426 bp of mitochondrial DNA (550 bp ND2, 671 bp ND4, 708 bp 12SrRNA, 497 bp 16SrRNA). The population numbers (Figure 2, Additional file 2) are provided in parentheses. (more ...)
Figure 5
Figure 5
Continuation of Figure 4.
Clade 1 encompasses populations from coastal northern and eastern QLD, and is comprised of three subclades (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 1a includes populations north of Townsville, subclade 1b contains populations in the Mackay region, and subclade 1c stretches from the Rockhampton region to Bania State Forest (inland from Bundaberg) in south-east QLD (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 2 is restricted to Kroombit Tops (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). A complex mosaic of geographically non-overlapping clades occurs throughout south-eastern QLD and northern NSW (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 3 includes populations from the Sunshine Coast and Brisbane region of south-eastern QLD (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). There is strong support for a close affinity between Clade 2 and Clade 3 (Figure (Figure4,4, ,5).5). Subclade 3a occurs within Cooloola National Park, while subclade 3b is restricted to the Bunya Mountains (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 3c extends from the Maryborough region to the northern suburbs of Brisbane, with subclade 3d containing populations from the southern suburbs of Brisbane to Barney View on the western side of Lamington National Park (Figure (Figure2,2, ,3,3, ,4,4, ,55).
Clade 4 (96 BS) occurs along the Main Range and throughout the QLD/NSW border region, and in inland northern NSW (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 4a (93 BS) extends from Deongwar State Forest and other areas in south-eastern QLD (Lamington NP, Main Range) through the elevated regions of inland northern NSW to Riamukka State Forest, inland from Port Macquarie (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 4b is restricted to the more coastal regions of northern NSW (Mt Warning NP, Border Ranges NP, Nightcap NP, Whian Whian State Forest, Alstonville), while subclade 4c occurs further inland at Bolivia Hill, Torrington State Forest and near Armidale (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 5 extends along the northern NSW coastal region from Yamba to Cairncross State Forest near Port Macquarie, and is comprised of two subclades: subclade 5a (Yamba to Coffs Harbour) and subclade 5b (Styx River State Forest to Cairncross State Forest) (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 6 represents a disjunct inland population on the Blackdown Tableland in southern QLD (Figure (Figure2,2, ,3,3, ,4,4, ,55).
Clade 7 is geographically widespread, occurring from the Australian Capital Territory (ACT) and inland southern NSW (subclade 7a) across to western VIC (Little Desert NP) and south-eastern SA (subclade 7b) (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 8 is a disjunct population that occurs in Coolah Tops National Park in inland northern NSW (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Clade 9 is distributed from the central coast of NSW and throughout eastern VIC and TAS (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 9a encompasses the central coast of NSW and the Sydney region, subclade 9b occurs at Brayton, while subclade 9c comprises populations from coastal southern NSW (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). Subclade 9d represents a shallow clade that is distributed throughout eastern VIC and TAS (Figure (Figure2,2, ,3,3, ,4,4, ,55).
Genetic differentiation among clades and divergence time estimates
Considerable genetic differentiation was evident amongst the nine L. delicata clades, with extremely high and statistically significant pairwise ΦST values among clades (Table (Table5).5). The only comparisons that were not significant were those involving clades with low sample sizes (e.g. Clade 8). Substantial genetic distances are evident among the clades (4.3-7.4%; Table Table5),5), indicating that the divergences within L. delicata occurred during the late Miocene-Pliocene (Figure (Figure6,6, Additional file 5). The intra-clade genetic divergences in L. delicata were 0.0-2.6% (Table (Table3).3). The vast majority (97.7%) of genetic variation in L. delicata was partitioned among populations (Table (Table6).6). The nine a priori biogeographic barriers (Tables (Tables11 and and2)2) accounted for 64.1% of the genetic variation in L. delicata (Table (Table6).6). This value increased to 66.5% when the two barriers identified a posteriori (Blackdown Tableland, Coolah Tops) were included in the analysis (Table (Table66).
Table 5
Table 5
Mean Tamura-Nei corrected mtDNA genetic distances (below diagonal) and pairwise ΦST (above diagonal) among the major clades (1-9) identified in Figure 4 and 5.
Figure 6
Figure 6
BEAST maximum clade credibility tree for Lampropholis delicata. The divergence times correspond to the mean posterior estimate of their age in millions of years. The bars indicate the 95% HPD interval for the divergence time estimates. Nodes with a posterior (more ...)
Table 6
Table 6
Hierarchical Analysis of Molecular Variance (AMOVA) for Lampropholis delicata populations and biogeographic regions.
The Bayesian skyline plot indicated recent (i.e. last 0.2 myr) contraction then expansion of L. delicata populations (Figure (Figure7),7), although there was no consistent support for the model of spatial expansion in L. delicata clades or subclades. Three main clades (2, 4 and 7) and four subclades (3a, 4b, 7b and 9d) deviated significantly from the expectations of neutrality (Tajima's D, Fu's F statistic) (Table (Table3),3), suggesting recent population expansion. However, the RI and SSD values indicated that a model of population expansion could only be conclusively rejected for two main clades (1 and 4) and two subclades (1c and 4c) (Table (Table33).
Figure 7
Figure 7
Bayesian skyline plot for Lampropholis delicata. The estimated effective population size over the past 6 myr is presented. The solid line indicates the median population estimate, with the two light blue lines indicating the 95% confidence interval.
The nine main clades of Lampropholis delicata appear to have diverged during the late Miocene-Pliocene. Although the current study relied solely on mitochondrial DNA sequence data, the same tree topology is evident in a molecular phylogeny for the Lampropholis genus (including representatives from each main L. delicata clade) based on mitochondrial DNA and five nuclear genes (C. Hoskin, C. Moritz & D. Chapple, unpublished data). The divergence of L. delicata corresponds to a time when rainforest habitat in eastern Australia was in decline as a result of a drying climate, resulting in restriction of rainforest to a series of disjunct remnants that have been described as an 'archipelago of refugia' [1,3,17]. The delicate skink occurs in rainforest or rainforest fringes and therefore likely experienced similar reduction and fragmentation, resulting in genetic divergence among geographically isolated populations. Despite evidence for the expansion and contraction of some clades throughout the Pleistocene, each is geographically structured and non-overlapping (Figure (Figure2,2, ,3).3). This pattern that has been observed in a range of other taxa (Table (Table2;2; [1]), including L. guichenoti [26]. Phylogeographic breaks in the delicate skink generally correspond to dry habitat and topographic barriers (i.e. Burdekin Gap, St Lawrence Gap, Kroombit Tops, McPherson Range; Hunter Valley, southern NSW; Figure Figure1,1, ,2,2, ,3).3). However, contrary to the hypothesis of Rawlinson [29], the delicate skink appears to be a relatively recent arrival in south-eastern Australia and exhibits no evidence of restricted geneflow across the barriers in this region (e.g. East Gippsland, Bass Strait).
Phylogeographic structure in the delicate skink corresponds to dry habitat and elevational barriers
Despite its widespread distribution along the east coast of Australia, there is substantial phylogeographic structure across the native range of the delicate skink (Figure (Figure4,4, ,5).5). In many instances these breaks are concordant with dry habitat corridors, indicating that regions of drier vegetation represent effective barriers to dispersal for the mesic-adapted delicate skink. For instance, the delicate skink exhibits a moderate genetic break (4.5%, subclade 1a vs 1b; Figure Figure6)6) across the Burdekin Gap in North QLD. Equivalent Pliocene divergences between populations either side of the Burdekin Gap have been reported in open forest frogs [66,67], rainforest lizards [22], woodland lizards [68], rainforest birds (Pleistocene; [19]) and freshwater fish (Miocene; [69]) (Table (Table2).2). In contrast, the St Lawrence Gap north of Rockhampton on the central QLD coast has only been identified as a significant barrier for one lizard species (late Pleistocene-early Pliocene, [68]; Table Table2).2). Although divergence is also evident across the St Lawrence Gap in the delicate skink (4.8%, mid-late Pliocene, subclade 1b vs 1c), a more substantial break is evident a little to the south, between clades 1 and 3 (7.1%, late Miocene) in the Gladstone region (Table (Table5,5, Figure Figure2,2, ,3,3, ,4,4, ,55).
Two high elevation areas (Kroombit Tops, Blackdown Tableland) inland from the main range of L. delicata in southern QLD were found to harbour genetically divergent lineages. Both areas are remnant patches of moist forest that are surrounded by drier lowland eucalypt woodland, and are disjunct from the main distribution of the delicate skink along the east coast ([70]; Figure Figure2,2, ,3).3). Kroombit Tops (~730 m) was identified a priori as a potential habitat isolate for the delicate skink, as the region provides a cooler and wetter refuge for mesic-adapted species ([70,71]; Tables Tables11 and and2).2). The Kroombit Tops population of the delicate skink diverged from the surrounding coastal populations during the mid Pliocene (4.9%, Clade 2 vs 3; Table Table5,5, Figure Figure6),6), a pattern that has also been observed in a rainforest bird [71], two open forest frogs [66,67], and several open forest reptiles [72] (Table (Table2).2). The Blackdown Tablelands are a moderate elevation plateau (~600 m) that provides an isolated refugium for numerous mesic-adapted species. Our results for the delicate skink (4.3-6.7%, Pliocene; Figure Figure6)6) provide evidence that the fauna of this region may also be genetically divergent from the coastal populations, a pattern also seen in other open forest reptiles [72].
The Coolah Tops are a high elevation (1200 m) plateau located in inland northern NSW, just to the north of the Hunter Valley (Figure (Figure2,2, ,3).3). The refugial population of the delicate skink that occurs on the Coolah Tops was found to have diverged from the nearby populations in northern NSW during the early-mid Pliocene (4.6-4.8%; Table Table5,5, Figure Figure6).6). The dry habitat in the Hunter River Valley has been demonstrated to represent a major barrier to dispersal in both woodland and wet forest species ([73-76]; Table Table2).2). While the divergence across the Hunter Valley was estimated to have occurred in the Miocene in the congeneric L. guichenoti, which is frequently sympatric with L. delicata [26], an early-mid Pliocene break was observed across this barrier in the delicate skink (subclade 9a vs Clades 4-5; Figure Figure2,2, ,3,3, ,4,4, ,5).5). The divergence estimate for the delicate skink is consistent with those reported for most other species across the Hunter Valley (Table (Table22).
Our analyses revealed a complex mosaic of geographically structured, non-overlapping clades and subclades (Clades 3-5) in the delicate skink in south-eastern QLD and northern NSW (Figure (Figure2,2, ,3,3, ,4,4, ,5).5). The McPherson Range that occurs along the border region (Figure (Figure1)1) is concordant with the phylogeographic break between Clades 3 and 4 (Figure (Figure2,2, ,3).3). The distribution of Clade 4 extends northwards into south-east QLD to the Main Range, which runs perpendicular to the western edge of the McPherson Range (Figure (Figure1).1). Similar biogeographic patterns involving the Main and McPherson Ranges occur in wet forest (Litoria pearsoniana; [77]) and open forest frogs (Litoria fallax; [66]) (Table (Table2).2). The early-mid Pliocene split found across the McPherson/Main Ranges in the delicate skink (4.8%; Figure Figure6)6) is concordant with that observed in L. guichenoti, which also inhabits open woodlands and dry sclerophyll forest [28], but intermediate between that reported for frogs (Miocene; [66,77]) and a wet forest snake (Pleistocene; [78]) (Table (Table22).
Some relatively minor phylogeographic structure is evident among the populations in the Maryborough, Sunshine Coast and Brisbane regions of south-east QLD (subclades 3a-d; Figure Figure2,2, ,3,3, ,4,4, ,5).5). A more substantial break occurs between inland (Clade 4) and coastal (Clade 5) delicate skink populations north of the Hunter Valley in NSW (4.5%, early-mid Pliocene; Table Table5,5, Figure Figure6).6). A similar coastal vs inland divergence in northern NSW is shared with White's skink (Liopholis whitii, [74]), and reflects an equivalent pattern that is regularly observed in southern NSW (Table (Table2).2). Indeed, an analogous pattern is evident within Clade 4 in northern NSW, with subclade 4b occurring near the coastal margin, subclade 4a present in intermediate areas (with secondary contact between 4a and 4b occurring in the Border Ranges NP), and subclade 4c occurring further inland (Figure (Figure2,2, ,3).3). These biogeographic patterns, combined with the break observed in southern NSW (Figure (Figure2,2, ,3),3), indicate that high elevation areas may represent barriers to dispersal in the delicate skink.
The delicate skink is a relatively recent arrival in southern Australia
The five phylogeographic studies that have had sufficient sampling to examine the impact of the elevational and habitat barriers in southern NSW have reported a genetic break between the inland (including the ACT) and coastal regions ([26,71,73,74,79]; Table Table2).2). The impact of this barrier is pronounced in the delicate skink, as populations from the ACT and inland southern NSW are more closely related to the SA populations (3.0%, early Pleistocene-late Pliocene) than the adjacent populations along the NSW coast (5.3%, mid-Pliocene; Figure Figure2,2, ,3,3, ,4,4, ,5).5). This indicates that the delicate skink most likely reached SA from the southern NSW region via an inland route, rather than along a coastal dispersal pathway through VIC. The delicate skink may have dispersed through the mesic vegetation that is located along the Murray River, which forms the border between NSW and VIC for the majority of its length (Figure (Figure2,2, ,3).3). Indeed, the eastern water skink (Eulamprus quoyii) has a continuous distribution through the Murray-Darling River system that connects populations along the east coast (North QLD to southern NSW) to an isolated population in south-eastern SA [28]. This biogeographic pattern was not previously suspected for the delicate skink and explains the large distributional gap across western VIC between the eastern suburbs of Melbourne and Little Desert NP in north-western VIC (Figure (Figure2,2, ,3).3). Given this pattern, it was not possible to examine the impact of the Murray Basin on the delicate skink (Table (Table22).
Several frog and lizard species exhibit deep phylogeographic breaks in the East Gippsland region [26,67,74,79], a pattern that is believed to be the result of repeated marine inundation of the area since the Miocene (Tables (Tables11 and and2).2). However, the East Gippsland region does not appear to constitute a barrier to dispersal in the delicate skink, with an extremely shallow clade (intraclade divergence 0.2%) distributed across eastern VIC and across Bass Strait into TAS (Figure (Figure2,2, ,3).3). The coastal area in East Gippsland has been relatively stable since the late Pleistocene [15,80], enabling the delicate to disperse across eastern VIC from southern NSW.
The delicate skink has colonised TAS during the late Pleistocene, with the presence of shared haplotypes between populations in eastern VIC and TAS indicating a connection between these two regions until relatively recently (~12-15 kya; Table Table5,5, Figure Figure6,6, Additional file 2). The timing coincides with the inundation of the Bassian Isthmus, which connected eastern VIC (Wilsons Promotory) to north-eastern TAS between 43-14 kya ([16]; Table Table1,1, Figure Figure1).1). Although the delicate skink does not currently occur on Wilsons Promontory, it is present on Flinders Island (which formed part of the Bassian Isthmus) and is common throughout north-eastern TAS ([[28,81], DGC, personal observation]). While a second land bridge was located from western VIC, through the King Island region to western TAS from ~43-17.5 kya [16], the absence of the delicate skink from western VIC would have precluded dispersal of the species via this western route. Fossil evidence for a Nothofagus tree species on King Island 38 kya suggests that moist forest vegetation occurred along the Bass Strait land bridges [82], enabling dispersal of the delicate skink into TAS. Although some other species appear not to have used these recent land bridges (frog: Crinia signigera, [79]; reptiles: Acritoscincus duperreyi, [76]; Lerista bougainvilli, [83]; mammals: Dasyurus maculatus, [84]), others appear to have dispersed across these Bass Strait land bridges (frogs: Limnodynastes peronii and tasmaniensis, [67]; reptiles: Liopholis whitii, [73]; Notechis scutatus, [85]) (Table (Table2).2). The repeated presence of the land bridges has also restricted east-west gene flow across Bass Strait in several marine invertebrate species ([86-89]; Table Table22).
There is no evidence to support the hypothesis of Rawlinson [29] that the delicate skink is a 'glacial relic' with a relatively long presence in southern Australia. In contrast, our analyses indicate that the delicate skink only colonised VIC and TAS during the late Pleistocene from coastal southern NSW. Although the delicate skink (also known as the rainbow or plague skink in its introduced range) is a successful invasive species in the Hawaiian Islands, New Zealand, and Lord Howe Island [90], there is no strong evidence to suggest that it represents an introduced species in TAS. However, given the relatively shallow genetic divergences within subclade 9d, we are unable to completely exclude the possibility that the delicate skink reached TAS via human-associated colonisation.
We performed a detailed phylogeographic study of a species found in mesic forests down almost the entire length of eastern Australia. Lampropholis delicata is a single widespread, but genetically variable, species consisting of nine geographically structured, non-overlapping clades. This structuring is likely the result of population subdivision across dry habitat barriers (Burdekin Gap, St Lawrence Gap, Hunter Valley), topographic barriers (McPherson Range, southern NSW) and to upland habitat isolates (Kroombit Tops, Blackdown Tableland, Coolah Tops). In contrast, in the south-east of its range, the delicate skink exhibits evidence for recent dispersal into SA via an inland route, and through eastern VIC and across the Bassian Isthmus into TAS. Previous studies have demonstrated geographic variation in morphology, reproductive ecology and life history in the delicate skink [91,92]. Given the presence of multiple divergent lineages across the range, this regional variation in morphology and life history may have a genetic, as well as climatic or environmental, basis.
Authors' contributions
DGC, MBT and CJH developed the project and obtained funding for the research. The fieldwork and tissue sample collection was conducted by DGC, CJH and SNJC. DGC and SNJC completed the sequencing and analyses. All authors contributed to writing the manuscript, with each reading and approving the final manuscript.
Supplementary Material
Additional file 1
Complete collection locality table, with museum specimen and tissue voucher number and GenBank accession numbers.
Additional file 2
Clades, haplotypes, latitude and longitude for Lampropholis delicata populations sampled in the study.
Additional file 3
Oligonucleotide primers used in this study.
Additional file 4
The concatenated alignment (fasta format) for the 165 Lampropholis delicata haplotypes (ND2: 1-550, ND4: 551-1221, 12SrRNA: 1222-1929, 16SrRNA: 1930-2426).
Click here for file(411K, FASTA)
Additional file 5
Divergence time estimates for the main Lampropholis delicata clades and subclades.
Acknowledgements
We thank C. Beatson, T. Bertozzi, D. Bray, N. Clemann, S. Donnellan, K. Gray, J. Herbert, J. Melville, R. Palmer and R. Sadlier for providing tissue samples. A number of samples were also provided by collections for the south-east QLD and northern NSW Regional Forest Agreement (RFA) projects. N. Clemann, M. Driessen, T. Gordon, R. Swain, E. Wapstra and G. While provided assistance with fieldwork and/or information and advice. C. Moritz was involved in early work on L. delicata phylogeography as part of the RFAs. We thank D. Bray, R. O'Brien and J. Sumner for their assistance in lodging the specimens at Museum Victoria. This research was conducted with the approval of the Museum Victoria Animal Ethics Committee (Approval No.: 07002), and in accordance with Victorian (Research Permit: 10004254, Import Permit: 13225709) and Tasmanian (Research Permit: FA07221, Export Permit: 8594/08) scientific research permits. The research was funded by the Australian Research Council (grant to DGC, Project Number DP0771913), the National Geographic Society (grant to DGC and MBT, CRE 8085-06), and the Allan Wilson Centre for Molecular Ecology and Evolution.
  • Byrne M. Evidence for multiple refugia at different time scales during Pleistocene climatic oscillations in southern Australia inferred from phylogeography. Quaternary Science Reviews. 2008;27:2576–2585. doi: 10.1016/j.quascirev.2008.08.032. [Cross Ref]
  • Powney GD, Grenyer R, Orme CDL, Owens IPF, Meiri S. Hot, dry and different: Australian lizard richness is unlike that of mammals, amphibians and birds. Global Ecology and Biogeography. 2010;19:386–396. doi: 10.1111/j.1466-8238.2009.00521.x. [Cross Ref]
  • Martin HA. Cenozoic climate change and the development of the arid vegetation in Australia. Journal of Arid Environments. 2006;66:533–563. doi: 10.1016/j.jaridenv.2006.01.009. [Cross Ref]
  • Byrne M, Yeates DK, Joseph L, Kearney M, Bowler J, Williams MAJ, Cooper S, Donnellan SC, Keogh JS, Leijs R. et al. Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biota. Molecular Ecology. 2008;17:4398–4417. doi: 10.1111/j.1365-294X.2008.03899.x. [PubMed] [Cross Ref]
  • Keast A. In: Ecological Biogeography of Australia. Keast A, editor. The Hague: Dr W. Junk; 1981. Distributional patterns, regional biotas, and adaptations in the Australian biota: a synthesis; pp. 1891–1997.
  • Frakes LA, McGowran B, Bowler JM. In: Fauna of Australia. Dyne GR, Walton DW, editor. 1A. Canberra: Australian Government Publishing Service; 1987. Evolution of Australian environments; pp. 1–16. General Articles.
  • Taylor G. In: History of the Australian vegetation: Cretaceous to recent. Hill RS, editor. Cambridge: Cambridge University Press; 1994. Landscapes of Australia: their nature and evolution; pp. 60–79.
  • Markgraf V, McGlone M, Hope G. Neogene paleoenvironmental and paleoclimatic change in southern temperate ecosystems - a southern perspective. Trends in Ecology and Evolution. 1995;10:143–147. doi: 10.1016/S0169-5347(00)89023-0. [PubMed] [Cross Ref]
  • Barrows TT, Stone JO, Fifield LK, Cresswell RG. The timing of the last glacial maxima in Australia. Quaternary Science Reviews. 2002;21:159–173. doi: 10.1016/S0277-3791(01)00109-3. [Cross Ref]
  • Bowler JM. In: Evolution of the Flora and Fauna of Arid Australia. Barker WR, Greenslade PJM, editor. Adelaide: Peacock Publications; 1982. Aridity in the late Tertiary and Quaternary of Australia; pp. 35–45.
  • Gallagher SJ, Smith AJ, Jonasson K, Wallace MW, Holdgate GR, Daniels J, Taylor D. The Miocene palaeoenivronmental and palaeoceanographic evolution of the Gippsland Basin, southeast Australia: a record of Southern Ocean change. Palaeogeography, Palaeoclimatology, Palaeoecology. 2001;172:53–80. doi: 10.1016/S0031-0182(01)00271-1. [Cross Ref]
  • Gallagher SJ, Greenwood DR, Taylor D, Smith AJ, Wallace MW, Holdgate GR. The Pliocene climatic and environmental evolution of southeastern Australia: evidence from the marine and terrestrial realm. Palaeogeography, Palaeoclimatology, Palaeoecology. 2003;193:349–382. doi: 10.1016/S0031-0182(03)00231-1. [Cross Ref]
  • Nix H. In: Evolution of the Flora and Fauna of Arid Australia. Barker WR, Greenslade PJM, editor. Adelaide: Peacock Publications; 1982. Environmental determinants of biogeography and evolution in Terra Australis; pp. 47–66.
  • Kershaw AP, Martin HA, McEwen Mason JRC. In: History of the Australian vegetation: Cretaceous to recent. Hill RS, editor. Cambridge: Cambridge University Press; 1994. The Neogene: a period of transition; pp. 299–327.
  • Holdgate GR, Wallace MW, Gallagher SJ, Smith AJ, Keene JB, Moore D, Shafik S. Plio-Pleistocene tectonics and eustasy in the Gippsland Basin, southeast Australia: evidence from magnetic imagery and marine geological data. Australian Journal of Earth Sciences. 2003;50:403–426. doi: 10.1046/j.1440-0952.2003.01004.x. [Cross Ref]
  • Lambeck K, Chappell J. Sea level change through the last glacial cycle. Science. 2001;292:679–686. doi: 10.1126/science.1059549. [PubMed] [Cross Ref]
  • White ME. After the greening. The browning of Australia. Sydney: Kangaroo Press; 1994.
  • Moritz C, Hoskin CJ, Graham CH, Hugall A, Moussalli A. In: Phylogeny and conservation. Purvis A, Gittleman JL, Brooks T, editor. Cambridge: Cambridge University Press; 2005. Historical biogeography, diversity and conservation of Australia's tropical rainforest herpetofauna; pp. 243–264.
  • Joseph L, Moritz C. Mitochondrial DNA phylogeography of birds in eastern Australian rainforests: first fragments. Australian Journal of Zoology. 1994;42:385–403. doi: 10.1071/ZO9940385. [Cross Ref]
  • Hugall A, Moritz C, Moussalli A, Stanisic J. Reconciling paleodistribution models and comparative phylogeography in the wet tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875) Proceedings of the National Academy of Sciences of the United States of America. 2002;99:6112–6117. doi: 10.1073/pnas.092538699. [PubMed] [Cross Ref]
  • Hoskin CJ, Higgie M, McDonald KR, Moritz C. Reinforcement drives rapid allopatric speciation. Nature. 2005;437:1353–1356. doi: 10.1038/nature04004. [PubMed] [Cross Ref]
  • Dolman G, Moritz C. A multilocus perspective on refugial isolation and divergence in rainforest skinks (Carlia) Evolution. 2006;60:573–582. [PubMed]
  • Bell KL, Moritz C, Moussalli A, Yeates DK. Comparative phylogeography and speciation of dung beetles from the Australian Wet Tropics rainforest. Molecular Ecology. 2007;16:4984–4998. doi: 10.1111/j.1365-294X.2007.03533.x. [PubMed] [Cross Ref]
  • Moussalli A, Moritz C, Williams SE, Carnaval AC. Variable responses of skinks to a common history of rainforest fluctuation: concordance between phylogeography and palaeo-distribution models. Molecular Ecology. 2009;18:483–499. doi: 10.1111/j.1365-294X.2008.04035.x. [PubMed] [Cross Ref]
  • Moritz C, Patton JL, Schneider CJ, Smith TB. Diversification of rainforest faunas: an integrated molecular approach. Annual Review of Ecology and Systematics. 2000;31:533–563. doi: 10.1146/annurev.ecolsys.31.1.533. [Cross Ref]
  • Chapple DG, Chapple SNJ, Thompson MB. Biogeographic barriers in south-eastern Australia drive phylogeographic divergence in the garden skink, Lampropholis guichenoti. Journal of Biogeography. 2011.
  • de Vis CW. A contribution to the herpetology of Queensland. Proceedings of the Linnean Society of New South Wales. 1888;2:811–826.
  • Wilson S, Swan G. A complete guide to Reptiles of Australia. Third. Sydney: Reed New Holland; 2010.
  • Rawlinson PA. In: Biogeography and Ecology in Tasmania. Williams WD, editor. The Hague: Junk, W; 1974. Biogeography and ecology of the reptiles of Tasmania and the Bass Strait area; pp. 291–338.
  • Green RH. Two skink lizards newly recorded from Tasmania. Records of the Queen Victoria Museum. 1965;19:1–4.
  • Reeder TW. A phylogeny of the Australian Sphenomorphus group (Scincidae: Squamata) and the phylogenetic placement of the crocodile skinks (Tribolonotus): Bayesian approaches to assessing congruence and obtaining confidence in maximum likelihood inferred relationships. Molecular Phylogenetics and Evolution. 2003;27:384–397. doi: 10.1016/S1055-7903(02)00448-7. [PubMed] [Cross Ref]
  • Greaves SNJ, Chapple DG, Daugherty CH, Gleeson DM, Ritchie PA. Genetic divergences pre-date Pleistocene glacial cycles in the New Zealand speckled skink, Oligosoma infrapunctatum. Journal of Biogeography. 2008;35:853–864. doi: 10.1111/j.1365-2699.2007.01848.x. [Cross Ref]
  • Liggins L, Chapple DG, Daugherty CH, Ritchie PA. Origin and post-colonization evolution of the Chatham Islands skink (Oligosoma nigriplantare nigriplantare) Molecular Ecology. 2008;17:3290–3305. doi: 10.1111/j.1365-294X.2008.03832.x. [PubMed] [Cross Ref]
  • Liggins L, Chapple DG, Daugherty CH, Ritchie PA. A SINE of restricted gene flow across the Alpine Fault: phylogeography of the New Zealand common skink (Oligosoma nigriplantare polychroma) Molecular Ecology. 2008;17:3668–3683. doi: 10.1111/j.1365-294X.2008.03864.x. [PubMed] [Cross Ref]
  • Greaves SNJ, Chapple DG, Gleeson DM, Daugherty CH, Ritchie PA. Phylogeography of the spotted skink (Oligosoma lineoocellatum) and green skink (O. chloronoton) species complex (Lacertilia: Scincidae) in New Zealand reveals pre-Pleistocene divergence. Molecular Phylogenetics and Evolution. 2007;45:729–739. doi: 10.1016/j.ympev.2007.06.008. [PubMed] [Cross Ref]
  • Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The Clustal-X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research. 1997;24:4876–4882. [PMC free article] [PubMed]
  • Posada D, Crandall KA. Modeltest: testing the model of DNA substitution. Bioinformatics. 1998;14:817–818. doi: 10.1093/bioinformatics/14.9.817. [PubMed] [Cross Ref]
  • Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sunderland, Massachusetts: Sinauer Associates; 2002.
  • Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biology. 2010;59:307–321. doi: 10.1093/sysbio/syq010. [PubMed] [Cross Ref]
  • Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–1574. doi: 10.1093/bioinformatics/btg180. [PubMed] [Cross Ref]
  • McGuire JA, Witt CC, Altshuler DL, Remsen JV. Phylogenetic systematics and biogeography of hummingbirds: Bayesian and maximum likelihood analyses of partitioned data and selection of an appropriate partitioning strategy. Systematic Biology. 2007;56:837–856. doi: 10.1080/10635150701656360. [PubMed] [Cross Ref]
  • Rambaut A, Drummond AJ. Tracer v1.5. http://beast.bio.ed.ac.uk/Tracer
  • Hillis DM, Bull JJ. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Systematic Biology. 1993;42:181–192.
  • Wilcox TP, Zwickl DJ, Heath TA, Hillis DM. Phylogenetic relationships of the dwarf boas and a comparison of Bayesian and bootstrap measures of phylogenetic support. Molecular Phylogenetics and Evolution. 2002;25:361–371. doi: 10.1016/S1055-7903(02)00244-0. [PubMed] [Cross Ref]
  • Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–2497. doi: 10.1093/bioinformatics/btg359. [PubMed] [Cross Ref]
  • Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Molecular Biology and Evolution. 2007;24:1596–1599. doi: 10.1093/molbev/msm092. [PubMed] [Cross Ref]
  • Excoffier LG, Lischer HEL. ARLEQUIN suite ver 3.5: a new series of programs to perform population genetic analyses under Linux and Windows. Molecular Ecology Resources. 2010;10:564–567. doi: 10.1111/j.1755-0998.2010.02847.x. [PubMed] [Cross Ref]
  • Excoffier LG, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–491. [PubMed]
  • Rice WR. Analysing tables of statistical tests. Evolution. 1989;43:223–225. doi: 10.2307/2409177. [Cross Ref]
  • Holm S. A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics. 1979;6:65–70.
  • Tajima F. Statistical method for testing the neutral mutation hypothesis. Genetics. 1989;123:585–595. [PubMed]
  • Fu YX. New statistical tests of neutrality for DNA samples from a population. Genetics. 1996;143:557–570. [PubMed]
  • Harpending HC. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Human Biology. 1994;66:591–600. [PubMed]
  • Rogers AR, Harpending HC. Population growth makes waves in the distribution of pairwise genetic differences. Molecular Biology and Evolution. 1992;9:552–569. [PubMed]
  • Ray N, Currat M, Excoffier LG. Intra-deme molecular diversity in spatially expanding populations. Molecular Biology and Evolution. 2003;20:76–86. doi: 10.1093/molbev/msg009. [PubMed] [Cross Ref]
  • Excoffier LG. Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model. Molecular Ecology. 2004;13:853–864. doi: 10.1046/j.1365-294X.2003.02004.x. [PubMed] [Cross Ref]
  • Schneider S, Excoffier LG. Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics. 1999;152:1079–1089. [PubMed]
  • Macey JR, Schulte JA, Ananjeva NB, Larson A, Rastegar-Pouyani N, Shammakov SM, Papenfuss TJ. Phylogenetic relationships among agamid lizards of the Laudakia caucasia species group: testing hypotheses of biogeographic fragmentation and an area cladogram for the Iranian Plateau. Molecular Phylogenetics and Evolution. 1998;10:118–131. doi: 10.1006/mpev.1997.0478. [PubMed] [Cross Ref]
  • Jennings WB, Pianka ER, Donnellan SC. Systematics of the lizard family Pygopodidae with implications for the diversification of Australian temperate biotas. Systematic Biology. 2003;52:757–780. [PubMed]
  • Poulakakis N, Lymberakis P, Valakos E, Pafilis P, Zouros E, Mylonas M. Phylogeography of Balkan wall lizard (Podarcis taurica) and its relatives inferred from mitochondrial DNA sequences. Molecular Ecology. 2005;14:2433–2443. doi: 10.1111/j.1365-294X.2005.02588.x. [PubMed] [Cross Ref]
  • Shoo LP, Rose R, Doughty P, Austin JJ, Melville J. Diversification patterns of pebble-mimic dragons are consistent with historical disruption of important habitat corridors in arid Australia. Molecular Phylogenetics and Evolution. 2008;48:528–542. doi: 10.1016/j.ympev.2008.03.022. [PubMed] [Cross Ref]
  • Bryson RW, Neito-Montes de Oca A, Velasco JC. Phylogenetic position of Porthidium hespere (Viperidae: Crotalinae) and phylogeography of arid-adapted hognosed pitvipers based on mitochondrial DNA. Copeia. 2008;2008:172–178. doi: 10.1643/CH-07-043. [Cross Ref]
  • Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007;7:214. doi: 10.1186/1471-2148-7-214. [PMC free article] [PubMed] [Cross Ref]
  • Skinner A, Hugall AF, Hutchinson MN. Lygosomine phylogeny and the origins of Australian scincid lizards. Journal of Biogeography. 2011;38:1044–1058. doi: 10.1111/j.1365-2699.2010.02471.x. [Cross Ref]
  • Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Molecular Biology and Evolution. 2005;22:1185–1192. doi: 10.1093/molbev/msi103. [PubMed] [Cross Ref]
  • James CH, Moritz C. Intraspecific phylogeography in the sedge frog Litoria fallax (Hylidae) indicates pre-Pleistocene vicariance of an open forest species from eastern Australia. Molecular Ecology. 2000;9:349–358. doi: 10.1046/j.1365-294x.2000.00885.x. [PubMed] [Cross Ref]
  • Schauble CS, Moritz C. Comparative phylogeography of two open forest frogs from eastern Australia. Biological Journal of the Linnean Society. 2001;74:157–170. doi: 10.1006/bijl.2001.0562. [Cross Ref]
  • Edwards DL, Melville J. Phylogeographic analysis detects congruent biogeographic patterns between a woodland agamid and Australian wet tropics taxa despite disparate evolutionary trajectories. Journal of Biogeography. 2010;37:1543–1556.
  • Wong BBM, Keogh JS, McGlashan DJ. Current and historical patterns of drainage connectivity in eastern Australia inferred from population genetic structuring in a widespread freshwater fish Pseudomugil signifer (Pseudomugilidae) Molecular Ecology. 2004;13:391–401. doi: 10.1046/j.1365-294X.2003.02085.x. [PubMed] [Cross Ref]
  • Czechura GV. Distant exiles: frogs and reptiles recorded from Kroombit Tops, southeastern Queensland. Queensland Naturalist. 1986;27:61–67.
  • Nicholls JA, Austin JJ. Phylogeography of an east Australian wet-forest bird, the satin bowerbird (Ptilonorhynchus violaceus), derived from mtDNA, and its relationship to morphology. Molecular Ecology. 2005;14:1485–1496. doi: 10.1111/j.1365-294X.2005.02544.x. [PubMed] [Cross Ref]
  • Moritz C. In: The Biology of Biodiversity. Kato M, editor. Hong Kong: Springer-Verlag; 2000. A molecular perspective on the conservation of diversity; pp. 21–34.
  • Donnellan SC, McGuigan K, Knowles R, Mahony M, Moritz C. Genetic evidence for species boundaries in frogs of the Litoria citropa species group. Australian Journal of Zoology. 1999;47:275–293. doi: 10.1071/ZO99013. [Cross Ref]
  • Chapple DG, Keogh JS, Hutchinson MN. Substantial genetic substructuring in south-eastern and alpine Australia revealed by molecular phylogeography of the Egernia whitii (Lacertilia: Scincidae) species group. Molecular Ecology. 2005;14:1279–1292. doi: 10.1111/j.1365-294X.2005.02463.x. [PubMed] [Cross Ref]
  • Moussalli A, Hugall AF, Moritz C. A mitochondrial phylogeny of the rainforest skink genus Saproscincus, Wells and Wellington (1984) Molecular Phylogenetics and Evolution. 2005;34:190–202. doi: 10.1016/j.ympev.2004.08.022. [PubMed] [Cross Ref]
  • Dubey S, Shine R. Evolutionary diversification of the lizard genus Bassiana (Scincidae) across southern Australia. PLoS One. 2010;5:e12982. doi: 10.1371/journal.pone.0012982. [PMC free article] [PubMed] [Cross Ref]
  • McGuigan K, McDonald K, Parris K, Moritz C. Mitochondrial DNA diversity and historical biogeography of a wet forest-restricted frog (Litoria pearsoniana) from mid-east Australia. Molecular Ecology. 1998;7:175–186. doi: 10.1046/j.1365-294x.1998.00329.x. [PubMed] [Cross Ref]
  • Keogh JS, Scott IAW, Fitzgerald M, Shine R. Molecular phylogeny of the Australian venomous snake genus Hoplocephalus (Serpentes, Elapidae) and conservation genetics of the threatened H. stephensii. Conservation Genetics. 2003;4:57–65. doi: 10.1023/A:1021823423944. [Cross Ref]
  • Symula R, Keogh JS, Cannatella DC. Ancient phylogeographic divergence in southeastern Australia among populations of the widespread common froglet, Crinia signifera. Molecular Phylogenetics and Evolution. 2008;47:569–580. doi: 10.1016/j.ympev.2008.01.011. [PubMed] [Cross Ref]
  • Hodgson L. In: Gondwana to Greenhouse: Australian Environmental Science. Gostin VA, editor. Geological Society of Australia Special Publication 21; 2001. Environmental geoscience issues in the Gippsland Basin, Victoria; pp. 161–169.
  • Hutchinson MN, Swain R, Driessen M. Snakes and lizards of Tasmania. Hobart: Nature Conservation Branch DPIWE and University of Tasmania; 2001.
  • Hope GS. In: History of the Australian vegetation: Cretaceous to present. Hill RS, editor. Cambridge: Cambridge University Press; 1994. Quaternary vegetation; pp. 368–389.
  • Fairbarn J, Shine R, Moritz C, Frommer M. Phylogenetic relationships between oviparous and viviparous populations of an Australian lizard (Lerista bougainvilli, Scincidae) Molecular Phylogenetics and Evolution. 1998;10:95–103. doi: 10.1006/mpev.1997.0468. [PubMed] [Cross Ref]
  • Firestone KB, Elphinstone MS, Sherwin WB, Houlden BA. Phylogeographical population structure of tiger quolls Dasyurus maculatus (Dasyuridae: Marsupiala), an endangered carnivorous marsupial. Molecular Ecology. 1999;8:1613–1625. doi: 10.1046/j.1365-294x.1999.00745.x. [PubMed] [Cross Ref]
  • Keogh JS, Scott IAW, Hayes C. Rapid and repeated origin of insular gigantism and dwarfism in Australian tiger snakes. Evolution. 2005;59:226–233. [PubMed]
  • Dawson MN. Incipient speciation of Catostylus mosaicus (Scyphozoa, Rhizostomeae, Catostylidae), comparative phylogeography and biogeography in south-east Australia. Journal of Biogeography. 2005;32:515–533. doi: 10.1111/j.1365-2699.2004.01193.x. [Cross Ref]
  • Waters JM, King TM, O'Loughlin PM, Spencer HG. Phylogeographical disjunction in abundant high-dispersal littoral gastropods. Molecular Ecology. 2005;14:2789–2802. doi: 10.1111/j.1365-294X.2005.02635.x. [PubMed] [Cross Ref]
  • York KL, Blacket MJ, Appleton BR. The Bassian Isthmus and the major ocean currents of southeast Australia influence the phylogeography and population structure of a southern Australian intertidal barnacle Catomerus polymerus (Darwin) Molecular Ecology. 2008;17:1948–1961. doi: 10.1111/j.1365-294X.2008.03735.x. [PubMed] [Cross Ref]
  • Waters JM. Marine biogeographical disjunction in temperate Australia: historical landbridge, contemporary currents, or both? Diversity and Distributions. 2008;14:692–700. doi: 10.1111/j.1472-4642.2008.00481.x. [Cross Ref]
  • Lever C. Naturalized reptiles and amphibians of the world. Oxford: Oxford University Press; 2003.
  • Forsman A, Shine R. Parallel geographic variation in body shape and reproductive life history within the Australian scincid lizard Lampropholis delicata. Functional Ecology. 1995;9:818–828. doi: 10.2307/2389979. [Cross Ref]
  • Forsman A, Shine R. The adaptive significance of colour pattern polymorphism in the Australian scincid lizard Lampropholis delicata. Biological Journal of the Linnean Society. 1995;55:273–291. doi: 10.1111/j.1095-8312.1995.tb01066.x. [Cross Ref]
  • Ford J. Hybrid zones in Australian birds. Emu. 1987;87:158–178. doi: 10.1071/MU9870158. [Cross Ref]
  • Ford J. Minor isolates and minor geographical barriers in avian speciation in continental Australia. Emu. 1987;87:90–102. doi: 10.1071/MU9870090. [Cross Ref]
  • Cracraft J. Patterns of diversification within continental biotas: heirarchial congruence among the areas of endemism of Australian vertebrates. Australian Systematic Botany. 1991;4:211–227. doi: 10.1071/SB9910211. [Cross Ref]
  • Crisp MD, Linder HP, Weston PH. Cladistic biogeography of plants in Australia and New Guinea: congruent pattern reveals two endemic tropical tracks. Systematic Biology. 1995;44:457–473.
  • Dickinson JA, Wallace MW, Holdgate GR, Gallagher SJ, Thomas L. Origin and timing of the Miocene-Pliocene unconformity in southeast Australia. Journal of Sedimentary Research. 2002;72:288–303. doi: 10.1306/082701720288. [Cross Ref]
  • Burns EL, Eldridge MDB, Crayn DM, Houlden BA. Low phylogeographic structure in a wide spread endangered Australian frog Litoria aurea (Anura: Hylidae) Conservation Genetics. 2007;8:17–32.
  • Brown M, Cooksley H, Carthew SM, Cooper SJB. Conservation units and phylogeographic structure of an arboreal marsupial, the yellow-bellied glider (Petaurus australis) Australian Journal of Zoology. 2006;54:305–317. doi: 10.1071/ZO06034. [Cross Ref]
  • Schiffer M, Kennington WJ, Hoffmann AA, Blacket MJ. Lack of genetic structure among ecologically adapted populations of an Australian rainforest Drosophila species as indicated by microsatellite markers and mitochondrial DNA sequences. Molecular Ecology. 2007;16:1687–1700. doi: 10.1111/j.1365-294X.2006.03200.x. [PubMed] [Cross Ref]
  • Jones ME, Shepherd M, Henry RJ, Delves A. Chloroplast DNA variation and population structure in the widespread forest tree, Eucalyptus grandis. Conservation Genetics. 2006;7:691–703. doi: 10.1007/s10592-005-9104-7. [Cross Ref]
Articles from BMC Evolutionary Biology are provided here courtesy of
BioMed Central