J. Scott Keogh

RAPID AND REPEATED ORIGIN OF INSULAR GIGANTISM AND DWARFISM IN AUSTRALIAN TIGER SNAKES

Abstract It is a well‐known phenomenon that islands can support populations of gigantic or dwarf forms of mainland conspecifics, but the variety of explanatory hypotheses for this phenomenon have been difficult to disentangle. The highly venomous Australian tiger snakes (genus Notechis) represent a well‐known and extreme example of insular body size variation. They are of special interest because there are multiple populations of dwarfs and giants and the age of the islands and thus the age of the tiger snake populations are known from detailed sea level studies. Most are 5000–7000 years old and all are less than 10,000 years old. Here we discriminate between two competing hypotheses with a molecular phylogeography dataset comprising approximately 4800 bp of mtDNA and demonstrate that populations of island dwarfs and giants have evolved five times independently. In each case the closest relatives of the giant or dwarf populations are mainland tiger snakes, and in four of the five cases, the closest relatives are also the most geographically proximate mainland tiger snakes. Moreover, these body size shifts have evolved extremely rapidly and this is reflected in the genetic divergence between island body size variants and mainland snakes. Within south eastern Australia, where populations of island giants, populations of island dwarfs, and mainland tiger snakes all occur, the maximum genetic divergence is only 0.38%. Dwarf tiger snakes are restricted to prey items that are much smaller than the prey items of mainland tiger snakes and giant tiger snakes are restricted to seasonally available prey items that are up three times larger than the prey items of mainland tiger snakes. We support the hypotheses that these body size shifts are due to strong selection imposed by the size of available prey items, rather than shared evolutionary history, and our results are consistent with the notion that adaptive plasticity also has played an important role in body size shifts. We suggest that plasticity displayed early on in the occupation of these new islands provided the flexibility necessary as the islands available prey items became more depauperate, but once the size range of available prey items was reduced, strong natural selection followed by genetic assimilation worked to optimize snake body size. The rate of body size divergence in haldanes is similar for dwarfs (hg= 0.0010) and giants (hg= 0.0020‐ 0.0025) and is in line with other studies of rapid evolution. Our data provide strong evidence for rapid and repeated morphological divergence in the wild due to similar selective pressures acting in different directions.

Molecular phylogeny and divergence dates for Australasian elapids and sea snakes (hydrophiinae): evidence from seven genes for rapid evolutionary radiations

One of the most prolific radiations of venomous snakes, the Australo‐Melanesian Hydrophiinae includes ∼100 species of Australasian terrestrial elapids plus all ∼60 species of viviparous sea snakes. Here, we estimate hydrophiine relationships based on a large data set comprising 5800 bp drawn from seven genes (mitochondrial: ND4, cytb, 12S, 16S; nuclear: rag1, cmos, myh). These data were analysed using parsimony, likelihood and Bayesian methods to better resolve hydrophiine phylogeny and provide a timescale for the terrestrial and marine radiations. Among oviparous forms, Cacophis, Furina and Demansia are basal to other Australian elapids (core oxyuranines). The Melanesian Toxicocalamus and Aspidomorphus group with Demansia, indicating multiple dispersal events between New Guinea and Australia. Oxyuranus and Pseudonaja form a robust clade. The small burrowing taxa form two separate clades, one consisting of Vermicella and Neelaps calanotus, and the other including Simoselaps, Brachyurophis and Neelaps bimaculatus. The viviparous terrestrial elapids form three separate groups: Acanthophis, the Rhinoplocephalus group and the Notechis–Hemiaspis group. True sea snakes (Hydrophiini) are robustly united with the Notechis–Hemiaspis group. Many of the retrieved groupings are consistent with previous molecular and morphological analyses, but the polyphyly of the viviparous and burrowing groups, and of Neelaps, are novel results. Bayesian relaxed clock analyses indicate very recent divergences: the ∼160 species of the core Australian radiation (including sea snakes) arose within the last 10 Myr, with most inter‐generic splits dating to between 10 and 6 Ma. The Hydrophis sea snake lineage is an exceptionally rapid radiation, with > 40 species evolving within the last 5 Myr.

Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches.

Evolutionary and biogeographic studies increasingly rely on calibrated molecular clocks to date key events. Although there has been significant recent progress in development of the techniques used for molecular dating, many issues remain. In particular, controversies abound over the appropriate use and placement of fossils for calibrating molecular clocks. Several methods have been proposed for evaluating candidate fossils; however, few studies have compared the results obtained by different approaches. Moreover, no previous study has incorporated the effects of nucleotide saturation from different data types in the evaluation of candidate fossils. In order to address these issues, we compared three approaches for evaluating fossil calibrations: the single-fossil cross-validation method of Near, Meylan, and Shaffer (2005. Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am. Nat. 165:137-146), the empirical fossil coverage method of Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am. Nat. 171:726-742), and the Bayesian multicalibration method of Sanders and Lee (2007. Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol. Lett. 3:275-279) and explicitly incorporate the effects of data type (nuclear vs. mitochondrial DNA) for identifying the most reliable or congruent fossil calibrations. We used advanced (Caenophidian) snakes as a case study; however, our results are applicable to any taxonomic group with multiple candidate fossils, provided appropriate taxon sampling and sufficient molecular sequence data are available. We found that data type strongly influenced which fossil calibrations were identified as outliers, regardless of which method was used. Despite the use of complex partitioned models of sequence evolution and multiple calibrations throughout the tree, saturation severely compressed basal branch lengths obtained from mitochondrial DNA compared with nuclear DNA. The effects of mitochondrial saturation were not ameliorated by analyzing a combined nuclear and mitochondrial data set. Although removing the third codon positions from the mitochondrial coding regions did not ameliorate saturation effects in the single-fossil cross-validations, it did in the Bayesian multicalibration analyses. Saturation significantly influenced the fossils that were selected as most reliable for all three methods evaluated. Our findings highlight the need to critically evaluate the fossils selected by data with different rates of nucleotide substitution and how data with different evolutionary rates affect the results of each method for evaluating fossils. Our empirical evaluation demonstrates that the advantages of using multiple independent fossil calibrations significantly outweigh any disadvantages.

Complex mating system and dispersal patterns in a social lizard, Egernia whitii

In contrast to the polygynous mating systems typically displayed by most reptilian taxa, long‐term genetic monogamy appears to be widespread within a lineage of group‐living Australian scincid lizards, the Egernia group. We have recently shown that Whites skink, Egernia whitii, lives in small but temporally stable social aggregations. Here, we examine the mating system, spatial organization, and dispersal patterns of E. whitii using behavioural field studies and data from four microsatellite loci. Parentage analysis of E. whitii litters revealed that its mating system is characterized by both polygyny and monogamy. Polygyny was the predominant mating system but within‐season social and genetic monogamy was common (36–45% of breeding pairs). The incidence of between‐season monogamy in E. whitii was rare compared to that reported for its congeners. Low levels of multiple paternity (12% of litters) and extra‐group paternity (16%) were detected. Social groups are generally comprised of closely related individuals, but breeding pairs were not more closely related compared to other potential mates. Spatial autocorrelation analyses revealed significant positive local genetic structure over 50 m, which was consistent for all age–sex classes. There was no clear and consistent evidence for sex‐biased dispersal, with assignment tests (mean assignment index) and relatedness analyses suggesting female‐biased dispersal, but spatial autocorrelation analyses indicating a trend for male‐biased dispersal. We discuss the implication of our results in regard to the factors promoting the evolution of monogamy within the Egernia group.

Substantial genetic substructuring in southeastern and alpine Australia revealed by molecular phylogeography of the Egernia whitii (Lacertilia: Scincidae) species group

Palaeoclimatic events and biogeographical processes since the mid‐Tertiary are believed to have strongly influenced the evolution and distribution of the terrestrial vertebrate fauna of southeastern Australia. We examined the phylogeography of the temperate‐adapted members of the Egernia whitii species group, a group of skinks that comprise both widespread low‐ to mid‐elevation (E. whitii) and montane‐restricted species (Egernia guthega, Egernia montana), in order to obtain important insights into the influence of past biogeographical processes on the herpetofauna of southeastern Australia. Sequence data were obtained from all six temperate‐adapted species within the E. whitii species group, and specifically from across the distributional ranges of E. whitii, E. guthega and E. montana. We targeted a fragment of the ND4 mitochondrial gene (696 bp) and analysed the data using maximum likelihood and Bayesian methods. Our data reveal a deep phylogeographical break in the east Gippsland region of Victoria between ‘northern’ (Queensland, New South Wales, Australian Capital Territory) and ‘southern’ (Victoria, Tasmania, South Australia) populations of E. whitii. This divergence appears to have occurred during the late Miocene–Pliocene, with the Gippsland basin possibly forming a geographical barrier to dispersal. Substantial structuring within both the ‘northern’ and the ‘southern’ clades is consistent with the effects of Plio–Pleistocene glacial‐interglacial cycles. Pleistocene glacial cycles also appear to have shaped the phylogeographical patterns observed in the alpine species, E. guthega and E. montana. We used our results to examine the biogeographical process that led to the origin and subsequent diversification of the lowland and alpine herpetofauna of southeastern Australia.

Palaeoclimate change drove diversification among isolated mountain refugia in the Australian arid zone

Refugia featured prominently in shaping evolutionary trajectories during repeated cycles of glaciation in the Quaternary, particularly in the Northern Hemisphere. The Southern Hemisphere instead experienced cycles of severe aridification but little is known about the temporal presence and role of refugia for arid‐adapted biota. Isolated mountain ranges located in the Australian arid zone likely provided refugia for many species following Mio/Pliocene (<15 Ma) aridification; however, the evolutionary consequences of the recent development of widespread sand deserts is largely unknown. To test alternative hypotheses of ancient vs. recent isolation, we generated a 10 gene data set to assess divergence history among saxicolous geckos in the genus Heteronotia that have distributions confined to major rocky ranges in the arid zone. Phylogenetic analyses show that each rocky range harbours a divergent lineage, and substantial intraspecific diversity is likely due to topographic complexity in these areas. Old divergences (∼4 Ma) among lineages pre‐date the formation of the geologically young sand deserts (<1 Ma), suggesting that Pliocene climate shifts fractured the distributions of biota long before the spread of the deserts.

Perched at the mito-nuclear crossroads: divergent mitochondrial lineages correlate with environment in the face of ongoing nuclear gene flow in an Australian bird.

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.

Molecular determination of paternity in a natural population of the multiply mating polygynous lizard Eulamprus heatwolei

We studied the mating system of the southern water skink, Eulamprus heatwolei, during spring and summer (encompassing the breeding season) in a population in southeastern Australia. We examined potential attributes that might influence the mating system and male reproductive success including home range size, physical proximity of adults and body size, and then genotyped all mothers, offspring and potential sires. Home range overlap of both sexes was extensive, with adult females sharing the greatest amount of space with each other and adult males the least amount of space with each other. However, not all adults hold home ranges. We classified approximately one quarter of adult males as home range holders and the rest as ‘floaters’. Adult females occupy home ranges more than males, with approximately three‐quarters classified as home range holders. Home range ownership is not correlated with body size for either sex, however, male body size is positively correlated with the number of adult female home ranges that his home range overlaps and adult male home ranges are larger than those of females. We used microsatellite genotyping to assign paternities to 55 offspring from 17 litters and then compared this data with our home range and behavioural observations. This species displays extreme levels of multiple paternity given the small mean clutch size of three. Multiple paternity was confirmed in 11 (64.7%) of 17 clutches but three other clutches (for a total of 82.4%) also may display multiple paternity. A total of 30 offspring from 12 litters were assigned to 10 of the 32 genotyped adult males from our study site. Of these 10 adult males, half were home range holders. Five complete clutches and a total of 25 out of the 55 offspring could not be positively assigned to any male surveyed as part of the study and were attributed to floater males or resident males adjacent to our study site that had not been genotyped. While sample sizes are small, neither male home range ownership nor body size is significantly correlated with the number of paternities a male obtained. Our study suggests a polygynous mating system for this species.

Spatial genetic analysis and long-term mark-recapture data demonstrate male-biased dispersal in a snake.

Dispersal is an important life-history trait, but it is notoriously difficult to study. The most powerful approach is to attack the problem with multiple independent sources of data. We integrated information from a 14-year demographic study with molecular data from five polymorphic microsatellite loci to test the prediction of male-biased dispersal in a common elapid species from eastern Australia, the small-eyed snake Rhinoplocephalus nigrescens. These snakes have a polygynous mating system in which males fight for access to females. Our demographic data demonstrate that males move farther than females (about twice as far on average, and about three times for maximum distances). This sex bias in adult dispersal was evident also in the genetic data, which showed a strong and significant genetic signature of male-biased dispersal. Together, the genetic and demographic data suggest that gene flow is largely mediated by males in this species.

Communal egg-laying in reptiles and amphibians: evolutionary patterns and hypotheses

Communal egg-laying is widespread among animals, occurring in insects, mollusks, fish, amphibians, reptiles, and birds, just to name a few. While some benefits of communal egg-laying may be pervasive (e.g., it saves time and energy and may ensure the survival of mothers and their offspring), the remarkable diversity in the life histories of the animals that exhibit this behavior presents a great challenge to discovering any general explanation. Reptiles and amphibians offer ideal systems for investigating communal egg-laying because they generally lack parental carea simplification that brings nest site choice behavior into sharp focus. We exhaustively reviewed the published literature for data on communal egg-laying in reptiles and amphibians. Our analysis demonstrates that the behavior is much more common than previously recognized (occurring in 481 spp.), especially among lizards (N = 255 spp.), where the behavior has evolved multiple times. Our conceptual review strongly suggests that different forces may be driving the evolution and maintenance of communal egg-laying in different taxa. Using a game theory approach, we demonstrate how a stable equilibrium may occur between solitary and communal layers, thus allowing both strategies to co-exist in some populations, and we discuss factors that may influence these proportions. We conclude by outlining future research directions for determining the proximate and ultimate causes of communal egg-laying.