Population demography of an endangered lizard, the Blue Mountains Water Skink
© Dubey et al.; licensee BioMed Central Ltd. 2013
Received: 11 October 2012
Accepted: 7 February 2013
Published: 13 February 2013
Information on the age structure within populations of an endangered species can facilitate effective management. The Blue Mountains Water Skink (Eulamprus leuraensis) is a viviparous scincid lizard that is restricted to < 40 isolated montane swamps in south-eastern Australia. We used skeletochronology of phalanges (corroborated by mark-recapture data) to estimate ages of 222 individuals from 13 populations.
These lizards grow rapidly, from neonatal size (30 mm snout-vent length) to adult size (about 70 mm SVL) within two to three years. Fecundity is low (mean 2.9 offspring per litter) and is affected by maternal body length and age. Offspring quality may decline with maternal age, based upon captive-born neonates (older females gave birth to slower offspring). In contrast to its broadly sympatric (and abundant) congener E. tympanum, E. leuraensis is short-lived (maximum 6 years, vs 15 years for E. tympanum). Litter size and offspring size are similar in the two species, but female E. leuraensis reproduce annually whereas many E. tympanum produce litters biennially. Thus, a low survival rate (rather than delayed maturation or low annual fecundity) is the key reason why E. leuraensis is endangered. Our 13 populations exhibited similar growth rates and population age structures despite substantial variation in elevation, geographic location and swamp size. However, larger populations (based on a genetic estimate of effective population size) contained older lizards, and thus a wider variance in ages.
Our study suggests that low adult survival rates, as well as specialisation on a rare and fragmented habitat type (montane swamps) contribute to the endangered status of the Blue Mountains Water Skink.
KeywordsAustralia Montane species Reptile Skeletochronology
Estimation of the age structure of endangered populations can facilitate their efficient management (e.g. [1–3]). Unfortunately, an animal’s age is difficult to calculate for many species. Mark-recapture analyses provide the most direct information, but require long-term studies, precluding any rapid conservation management plans . An alternative way to estimate an individual’s age is skeletochronology, based on histological analyses of growth marks in the skeleton [4, 5]. This technique relies upon seasonal variation in rates of skeletal growth, and thus is especially effective for ectothermic vertebrates that live in highly seasonal environments [6, 7].
The Blue Mountains Water Skink (Eulamprus leuraensis, Wells & Wellington, 1983) is a medium-sized (total length to 20 cm) viviparous scincid lizard that is restricted to less than forty small montane swamps (typically, < 2 ha) at 560 to 1,060 m elevation in the Blue Mountains and Newnes Plateau, west of the city of Sydney in south-eastern Australia. This species is classified as “endangered” under the IUCN Red List , the Threatened Species Conservation Act  and the Environmental Protection and Biodiversity Conservation Act . This lizard species is an ecological specialist, and comprises isolated small populations that are subject to considerable ongoing threats. In addition, the species’ entire distribution (montane swamps; < 2,000 km2) is listed as threatened under the Threatened Species Conservation Act  due to impacts of nearby urbanisation, (i.e., invasion by weeds, modification of the hydrological system and of bushfire regimes, pollutants and longwall mining).
These threats to population persistence are exacerbated by the low vagility of the lizards. Gene flow between populations is very limited, and animals within each of the two major parts of the species’ distribution (Blue Mountains versus Newnes Plateau) have been isolated from each other for at least a million years [11, 12]. In addition, this species has low annual fecundity (one to five neonates per annual litter; [13, 14]), further reducing the ability of populations to recover from the effects of episodes of higher-than-usual mortality [15, 16].
In the present study, we used phalangeal skeletochronology (supported by field mark-recapture data) to estimate the age structure of 13 populations of Blue Mountains Water Skinks. As well as clarifying basic issues such as growth rates and age at maturity, we took advantage of data from our earlier studies to explore the effects of maternal age on reproductive output and progeny quality (based on data from Dubey et al.,  and Dubey & Shine, ) and to compare population age structure to various site-specific parameters (including indices of genetic diversity from Dubey & Shine, [11, 12]).
Our multiple linear regression including the age and size (SVL) of gravid females as explanatory variables and the litter size as the response variable, revealed a significant effect of a female’s body size (F1,39 = 28.99, P < 0.0001) and of the interaction term between age and body size (F1,39 = 6.16, P = 0.018) on her litter size, as well as a marginally significant effect of her age (F1,39 = 3.62, P = 0.065). Older and larger females produced larger litters.
We also found a significant effect of these parameters on mean offspring size (SVL), with larger females producing larger neonates and older females smaller neonates (age: F1,34 = 8.79, P = 0.006; size: F1,34 = 4.63, P = 0.039; interaction: F1,34 = 5.59, P < 0.025). Similarly, larger females produced heavier neonates whereas older females gave birth to lighter neonates (age: F1,34 = 4.65, P = 0.039; size: F1,34 = 1.01, P = 0.32; interaction: F1,34 = 6.24, P = 0.018).
Mean age (with variance in parentheses) of Eulamprus leuraensis in our study populations
Blue Mountains Water Skinks grow fast and die young. The largest cohort of adult females comprised those in their first reproductive year post-maturation. Older females were rare, and (because they ceased growing) did not produce larger litters than their younger counterparts. Also, the offspring of older females exhibited reduced locomotor performance. In combination, these traits suggest that the reproductive capacity of a population of Blue Mountains Water Skinks rests primarily upon the output of newly-matured animals. This result provides a striking contrast to the demography of a congeneric water skink that is also found at relatively high elevations (and indeed, occurs in the eucalypt woodland of the Blue Mountains). The highland water skink E. tympanum is an abundant and widely-distributed taxon that has been the subject of detailed ecological studies, including skeletochronological work ([17–22]). Tilley  reported females living up to at least 15 years of age in a Victorian population, and Blomberg & Shine  concluded that adults lived more than 10 years in a NSW population. Mean ages of reproductive adults within E. tympanum populations were between four and seven years (in low and high elevation sites, respectively), compared to three to four years in the high-elevation E. leuraensis (present study). Offspring mass and litter size are similar in the two species. However, the proportion of females reproducing is higher in E. leuraensis than in E. tympanum (91% versus 30 to 60%, respectively; [13, 19]).
In many lizard and snake species, reproductive output increases strongly with maternal age, primarily because continued growth translates into larger females, that reproduce more often and produce more offspring when they do so (e.g., ). Blue Mountains Water Skinks do not exhibit such an increase. Instead, reproductive frequency is as high from the outset (as soon as females mature) as thereafter, because body size (the main determinant of litter size) does not increase substantially during adult life. Indeed, offspring quality (and thus, the benefits to maternal fitness) may decline in older females, based on the reduced locomotor capacities and body sizes of their offspring. Given a mean litter size of 2.9 [13, 14], the average lifetime production of neonates per female (even if she reaches six years of age) is less than 15. The short adult lifespan, coupled with the small litter size, must reduce the rate that a population can recover after its numbers have been reduced by events such as bushfires, heavy predation, or prolonged drought. Clearly, the survival rate of the newborns will have a major impact on the rate of population recovery.
The age structure of our populations suggests that many of the lizards are semelparous: that is, they reproduce in only a single year before dying. Like most cool-climate viviparous squamates, Blue Mountains Water Skinks produce only a single litter each year (although females in tropical populations of the congeneric E. quoyii can produce two litters per year: Schwarzkopf, unpublished). Although most lizards and snakes are iteroparous, semelparity is not uncommon, especially in harsh habitats such as hot and dry areas as well as at high elevations (e.g. [24–28]).
Although the 13 populations that we sampled spanned a broad range in terms of abiotic factors such as elevation, geographic location and swamp size, none of these factors had any detectable effect on the age structure of lizard populations. However, a genetic estimate of the effective population size (Theta k; based on the study of Dubey & Shine, ) was positively correlated with mean age and with intrapopulation variation in age. This pattern suggests that areas that confer higher survival of individual lizards (e.g., due to low predation rates or high food availability) thereby result in larger effective population sizes. Intuition suggests that larger swamps should contain more lizards, but our analyses show no such pattern (i.e., neither Theta k nor age structure were correlated with swamp size). Habitat quality may be more important than quantity. A low Theta k value also reflects low genetic diversity, which could directly reduce survival rates via inbreeding effects [29, 30]. Future conservation efforts could usefully explore the biotic and abiotic factors that render swamps more or less favorable to Blue Mountains Water Skinks.
Our skeletochronological data also identify a puzzling feature of this endangered lizard’s life history: its low adult survival rate. At a proximate level, this low survival rate must contribute to the precarious conservation status of the Blue Mountains Water Skink. If that low survival rate was driven by extrinsic factors such as local predation pressure or fire history, we would have expected to see at least some populations with much older animals. Thus, our data suggest that either the entire range of this species was affected by some major catastrophe about seven years prior to our study, or else low adult survival rates are a consistent feature of the biology of this species. We know of no obvious historical catastrophe that could explain the truncation of age structures across the entire range of the species. Thus, the most likely scenario is that, unlike its congeners such as E. tympanum, the Blue Mountains Water Skink is a relatively short-lived species.
The analysis above adds a demographic factor to the traits already identified as contributing to the endangerment of E. leuraensis. These lizards are habitat specialists, dependent upon a relatively scarce and highly fragmented set of highland swamps [11, 31, 32]. Their low vagility reduces gene flow among swamps, and retards recolonisation of any swamps where local populations have been extirpated [12, 13]. They live in an area that abuts Australia’s largest city, and hence are exposed to myriad threats caused by the proximity of humans [11, 31, 32]. Their dependence on moist habitats at high elevations means that they are at a high risk of habitat degradation due to climate change (southeastern Australia is predicted to become hotter and drier: ). Lastly, they exhibit the unfortunate combination of low reproductive output and low rates of adult survival. In combination, these characteristics render the Blue Mountains Water Skink highly vulnerable to any additional pressures (e.g., [13, 34, 35]), and suggest that the future conservation and management of this species will require close scrutiny and perhaps, active intervention to maintain suitable habitat conditions.
The toe-clips also provided an opportunity to mark individuals with a unique toe-clip code and hence to assess growth rates (and thus, ages relative to body size) by calculating growth increments for the 27 individuals that we recaptured in successive years. We used the von Bertalanffy growth model [36, 37] to estimate the relationship between size and age for both males and females. We fitted mark-recapture data to these models using the non-linear least-squares regression procedure in JMP 8.0.1 [37–39]. The parameters k (intrinsic growth rate) and a (asymptotic length) were seeded with initial best guesses.
Skeletochronological age estimation
Laboratory procedures for processing humeri and the second phalanx of the toes followed the standard methods of skeletochronology for amphibians and lizards [4, 40]: (1) decalcification of bones (5% formic acid, 1 h), (2) fixation in Bouin´s solution (at least 1 h), (3) Historesin™ (JUNG) embedding, (4) cross sectioning of the diaphysis at 10 μm using a JUNG RM2055 rotation microtome, (5) staining with 0.05% cresylviolet (5-10 min), (6) light microscopic examination using an Olympus BX 50 with ocular micrometer. The number of lines of arrested growth (LAG) was counted in the periosteal bone of those diaphysis sections in which the size of the medullar cavity was at its minimum and that of bone at its maximum. The number of usable diaphysis sections per individual varied between 2 and 16. All were considered for LAG counting and two of the authors independently obtained age estimates for each individual. Endosteal bone was present in all individuals, but the first LAG always remained visible (maximum resorption: ca. 30% of the complete line).
Statistical analyses were conducted using JMP 8.0 . Linear regressions were performed to compare the mean age of individuals within populations to longitude, latitude, elevation, swamp size and an estimate of the effective population size based on genetic analyses from Dubey & Shine (; Theta k, estimated from the infinite-allele equilibrium relationship  between the sample size [n], the expected number of alleles [k], and Theta [population parameter; Theta = 2Mu where M is equal to 2 N for diploid populations of size N, and u is the overall mutation rate at the haplotype level]; ).
We also performed multiple linear regressions between the age and size (SVL) of gravid females and their reproductive output (litter size and offspring size), including hatchling locomotor performance (sprint speed: see Dubey & Shine,  for methods). We omitted variables (incubation treatments, litter sizes) that did not affect locomotor speeds in this system . Mean values were calculated for each litter to avoid pseudoreplication. Data on reproductive variables were obtained from Dubey & Shine  and Dubey et al.. Based on these data and paternity analyses from Dubey et al., we also estimated the mean age of males that were known to have sired offspring.
We thank B. Nilow for technical assistance, M. Hensen for support and encouragement, J. Goudet and N. Perrin for helpful comments, and the Blue Mountains City Council, the Australian Research Council, and the Swiss National Science Foundation for funding. Tissue samples were taken under approvals from the University of Sydney Animal Care and Ethics Committee (authorization N° L04/12-2008/1/4956) and the Blue Mountains City Council.
- Chinsamy A, Valenzuela N: Skeletochronology of the endangered side-neck turtle, Podocnemis expansa. S Afr J Sci. 2008, 104: 311-314.Google Scholar
- Gillespie G: Population age structure of the spotted tree frog (Litoria spenceri): insights into population declines. Wildlife Res. 2010, 37: 19-26. 10.1071/WR08178.View ArticleGoogle Scholar
- Yang W, Liu C, Jiang J, Li C, Xie F: Age structure of females in a breeding population of Echinotriton chinhaiensis (Caudata: Salamandridae) and its conservation implication. Asian Herpetol Res. 2011, 2: 91-96. 10.3724/SP.J.1245.2011.00091.View ArticleGoogle Scholar
- Sinsch U, di Tada IE, Martino A: Longevity, demography and sex-specific growth of the Pampa de Achala toad, Bufo achalensis Cei, 1972. Stud. Neotrop. Fauna E. 2001, 36: 95-104. 10.1076/snfe.220.127.116.114.View ArticleGoogle Scholar
- Oromi N, Sanuy D, Sinsch U: Altitudinal variation of demographic life-history traits does not mimic latitudinal variation in natterjack toads (Bufo calamita). Zoology. 2012, 115: 30-7. 10.1016/j.zool.2011.08.003.View ArticlePubMedGoogle Scholar
- Castanet J, Smirina E: Introduction to the skeletochronological method in amphibians and reptiles. Ann Sci Nat, 13. ser Zool. 1990, 11: 191-196.Google Scholar
- Castanet J, Francillon-Vieillot H, Meunier FJ, de Ricqles A: Bone and individual aging. Bone growth. Edited by: Hall BK. 1993, Boca Raton: CRC Press, 245-283. 7Google Scholar
- IUCN: IUCN Red List of Threatened Species. Version. 2011, http://www.iucnredlist.org. 2011. Downloaded on 09 May 2012, .2Google Scholar
- NSW: Threatened Species Conservation Act. 1995, Sydney: Department of Environment and Climate ChangeGoogle Scholar
- Commonwealth: Environmental Protection and Biodiversity Conservation Act. 1999, Canberra: Australian GovernmentGoogle Scholar
- Dubey S, Shine R: Restricted dispersal and genetic diversity in populations of an endangered montane lizard (Eulamprus leuraensis, Scincidae). Mol Ecol. 2010, 19: 886-897. 10.1111/j.1365-294X.2010.04539.x.View ArticlePubMedGoogle Scholar
- Dubey S, Shine R: Plio-Pleistocene diversification and genetic population structure of an endangered lizard (the Blue Mountains water skink, Eulamprus leuraensis) in southeastern Australia. J Biogeogr. 2010, 37: 902-914. 10.1111/j.1365-2699.2009.02266.x.View ArticleGoogle Scholar
- Dubey S, Chevalley M, Shine R: Sexual dimorphism and sexual selection in a montane scincid lizard (Eulamprus leuraensis). Austr Ecol. 2011, 36: 68-75. 10.1111/j.1442-9993.2010.02119.x.View ArticleGoogle Scholar
- Dubey S, Shine R: Predicting the effects of climate change on an endangered montane lizard, Eulamprus leuraensis (Scincidae). Clim Chang. 2011, 17: 531-547.View ArticleGoogle Scholar
- McKinney ML: Extinction vulnerability and selectivity: combining ecological and paleontological views. Annu Rev Ecol Evol S. 1997, 28: 495-516. 10.1146/annurev.ecolsys.28.1.495.View ArticleGoogle Scholar
- Williams SE, Shoo LP, Isaac JL, Hoffmann AA, Langham G: Towards an integrated framework for assessing the vulnerability of species to climate change. PLoS Biol. 2008, 12: 2621-2626.Google Scholar
- Tilley SJ: Skeletal variation in the Australian Sphenomorphus group (Lacertilia: Scincidae). 1984, La Trobe University, Melbourne, Victoria, Australia: DissertationGoogle Scholar
- Schwarzkopf L: Sexual dimorphism in body shape without sexual dimorphism in body size in water skinks (Eulamprus quoyii). Herpetologica. 1992, 61: 116-123.View ArticleGoogle Scholar
- Schwarzkopf L: Costs of reproduction in water skinks. Ecology. 1993, 74: 1970-1981. 10.2307/1940840.View ArticleGoogle Scholar
- Doughty PD, Shine R: Detecting life-history tradeoffs: measuring energy stores in "capital" breeders reveals costs of reproduction. Oecologia. 1997, 110: 508-513. 10.1007/s004420050187.View ArticleGoogle Scholar
- Rohr DH: Demographic and life history variation in two proximate populations of a viviparous skink separated by a steep altitudinal gradient. J Anim Ecol. 1997, 66: 567-578. 10.2307/5950.View ArticleGoogle Scholar
- Blomberg S, Shine R: Modeling life history strategies with capture-recapture data: evolutionary demography of the water skink Eulamprus tympanum. Austr Ecol. 2001, 26: 349-359. 10.1046/j.1442-9993.2001.01120.x.View ArticleGoogle Scholar
- Shine R: Reproductive strategies in snakes. P. Roy. Soc. B. 2003, 270: 995-1004. 10.1098/rspb.2002.2307.View ArticleGoogle Scholar
- Shine R: "Costs" of reproduction in reptiles. Oecologia. 1980, 46: 92-100. 10.1007/BF00346972.View ArticleGoogle Scholar
- Bonnet X, Lourdais O, Shine R, Naulleau G: Reproduction in a typical capital breeder: costs, currencies, and complications in the aspic viper. Ecology. 2002, 83: 2124-2135. 10.1890/0012-9658(2002)083[2124:RIATCB]2.0.CO;2.View ArticleGoogle Scholar
- Karsten KB, Andriamandimbiarisoa LN, Fox SF, Raxworthy CJ: A unique life history among tetrapods: an annual chameleon living mostly as an egg. P Natl Acad Sci USA. 2008, 105: 8980-8984. 10.1073/pnas.0802468105.View ArticleGoogle Scholar
- Rodríguez-Romero F, Smith GR, Méndez-Sánchez F, Hernández-Gallegos O, Sánchez Nava P, de la Cruz FR M: Demography of a Semelparous, high-elevation population of Sceloporus bicanthalis (Lacertilia: Phrynosomatidae) from the Nevado de Toluca Volcano, Mexico. Southwest Nat. 2011, 56: 71-77. 10.1894/GC-193.1.View ArticleGoogle Scholar
- Bonnet X: The evolution of semelparity. Reproductive biology and phylogeny of snakes. Reproductive Biology and Phylogeny Series. Science Publishers Inc. 2011, 17: 645-672.View ArticleGoogle Scholar
- Vandewoestijne S, Schtickzelle N, Baguette M: Positive correlation between genetic diversity and fitness in a large, well-connected metapopulation. BMC Biol. 2008, 6: 46-10.1186/1741-7007-6-46.PubMed CentralView ArticlePubMedGoogle Scholar
- Valtonen TM, Roff DA, Rantala MJ: Analysis of the effects of inbreeding on lifespan and starvation resistance in Drosophila melanogaster. Genetica. 2011, 139: 525-533. 10.1007/s10709-011-9574-0.View ArticlePubMedGoogle Scholar
- Commonwealth: Environmental Protection and Biodiversity Conservation Act. 1999, Canberra: Australian GovernmentGoogle Scholar
- NSW National Parks & Wildlife Service: Blue Mountains water skink (Eulamprus leuraensis) recovery plan. 2001, Hurstville, NSW: NSW NPWSGoogle Scholar
- Dubey S, Pike D, Shine R: Predicting the impacts of climate change on genetic diversity in an endangered lizard species. Clim Chang. 2012, 10.1007/s10584-012-0540-3.Google Scholar
- Jiguet F, Gadot AS, Julliard R, Newson SE, Couvet D: Climate envelope, life history traits and the resilience of birds facing global change. Global Change Biol. 2007, 13: 1672-1684. 10.1111/j.1365-2486.2007.01386.x.View ArticleGoogle Scholar
- Forcada J, Trathan PN, Murphy EJ: Life history buffering in Antarctic mammals and birds against changing patterns of climate and environmental variation. Global Change Biol. 2008, 14: 2473-2488.Google Scholar
- Schoener TW, Schoener A: Estimating and interpreting body-size growth in some Anolis lizards. Copeia. 1978, 1978: 390-405. 10.2307/1443602.View ArticleGoogle Scholar
- Webb JK, Pike DA, Shine R: Population ecology of the velvet gecko, Oedura lesueurii in south eastern Australia: implications for the persistence of an endangered snake. Austr Ecol. 2008, 33: 839-847. 10.1111/j.1442-9993.2008.01852.x.View ArticleGoogle Scholar
- SAS Institute: User Guide. 2009, SAS Campus Drive, Cary, NC 27513: JMP, A Business Unit of SASGoogle Scholar
- Frazer NB, Gibbons JW, Greene JL: Exploring Fabens’ growth interval model with data on a long-lived vertebrate, Trachemys scripta (Reptilia: Testudinata). Copeia. 1990, 1990: 112-118. 10.2307/1445827.View ArticleGoogle Scholar
- Sinsch U, Martino A, di Tada IE: Longevity and sexual size dimorphism of the Pampa de Achala copper lizard Pristidactylus achalensis (Gallardo, 1964). Amphibia-Reptilia. 2002, 23: 177-190. 10.1163/156853802760061822.View ArticleGoogle Scholar
- Ewens WJ: The sampling theory of selectively neutral alleles. Theor Popul Biol. 1972, 3: 87-112. 10.1016/0040-5809(72)90035-4.View ArticlePubMedGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.