Diversification across biomes in a continental lizard radiation

Whip it into shape: Revision of the Demansia psammophis (Schlegel, 1837) complex (Squamata: Elapidae), with a description of a new species from central Australia

The genus Demansia Günther is the most diverse genus of Australian terrestrial elapids. A phylogenetic framework for the familiar but problematic taxa D. psammophis and D. reticulata (Gray) has been long overdue to ascertain interspecific relationships and resolve unclear taxonomic issues. Here we present an integrated molecular phylogenetic and morphological analyses to review species delineation, resulting in confirmation that both D. psammophis and D. reticulata are full species and that some populations referred to D. r. cupreiceps Storr are not distinguishable from more typical D. reticulata. We also find the widespread central Australian population (treated by most authors as part of cupreiceps) to be specifically distinct. We redescribe D. psammophis and D. reticulata to clarify morphological and geographical boundaries and describe D. cyanochasma sp. nov. based on a combination of molecular genetic markers, details of colour and pattern, adult total length and a few morphometric attributes. We also designate a lectotype for D. psammophis from the original syntype series and comment on the necessity for further taxonomic refinement of this distinctive group.

Initiation of speciation across multiple dimensions in a rock-restricted, tropical lizard

Population isolation and concomitant genetic divergence, resulting in strong phylogeographical structure, is a core aspect of speciation initiation. If and how speciation then proceeds and ultimately completes depends on multiple factors that mediate reproductive isolation, including divergence in genomes, ecology and mating traits. Here we explored these multiple dimensions in two young (Plio-Pleistocene) species complexes of gekkonid lizards (Heteronotia) from the Kimberley-Victoria River regions of tropical Australia. Using mitochondrial DNA screening and exon capture phylogenomics, we show that the rock-restricted Heteronotia planiceps exhibits exceptional fine-scale phylogeographical structure compared to the codistributed habitat generalist Heteronotia binoei. This indicates pervasive population isolation and persistence in the rock-specialist, and thus a high rate of speciation initiation across this geographically complex region, with levels of genomic divergence spanning the "grey zone" of speciation. Proximal lineages of H. planiceps were often separated by different rock substrates, suggesting a potential role for ecological isolation; however, phylogenetic incongruence and historical introgression were inferred between one such pair. Ecomorphological divergence among lineages within both H. planiceps and H. binoei was limited, except that limestone-restricted lineages of H. planiceps tended to be larger than rock-generalists. By contrast, among-lineage divergence in the chemical composition of epidermal pore secretions (putative mating trait) exceeded ecomorphology in both complexes, but with less trait overlap among lineages in H. planiceps. This system-particularly the rock-specialist H. planiceps-highlights the role of multidimensional divergence during incipient speciation, with divergence in genomes, ecomorphology and chemical signals all at play at very fine spatial scales.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Plio–Pleistocene vicariance across arid Australia in the ‘Spiny Knob-tailed Geckos’ (

Across Australia’s monsoon tropics and vast arid zone isolated regions or ‘islands’ of upland or rocky habitat are home to disjunct populations of many taxa of plants and animals. Comparative analyses of lineages that occur across these habitat islands provide opportunities to understand when and how environmental change drove isolation and diversification across arid Australia. Here we present an analysis of mitochondrial genetic diversity across disjunct populations of geckos in the Nephrurus asper group. Dating analyses suggest that disjunct and genetically divergent populations spanning the northern half of Australia diverged through the Plio–Pleistocene. Based on the timing of divergence and current habitat associations we hypothesise that species in this lineage were isolated by the expansion of unsuitable arid-zone habitats from the late Pliocene onwards. Across most areas, these barriers appear to be sandy or stony deserts. However, in eastern Australia genetically divergent populations are separated by grassland on flat vertisol-dominated soils (‘blacksoils’), suggesting that these habitats also expanded during the late Pliocene aridification. Finally, we show that western Queensland populations formerly referred to N. asper are genetically divergent and diagnosable on the basis of colour pattern and, herein, recognise these populations as a distinct species. https://zoobank.org/urn:lsid:zoobank.org:pub:9508CAAA-D014-452D-A3DA-325851615FA7

Quantitatively defining species boundaries with more efficiency and more biological realism

We introduce a widely applicable species delimitation method based on the multispecies coalescent model that is more efficient and more biologically realistic than existing methods. We extend a threshold-based method to allow the ancestral speciation rate to vary through time as a smooth piecewise function. Furthermore, we introduce the cutting-edge proposal kernels of StarBeast3 to this model, thus enabling rapid species delimitation on large molecular datasets and allowing the use of relaxed molecular clock models. We validate these methods with genomic sequence data and SNP data, and show they are more efficient than existing methods at achieving parameter convergence during Bayesian MCMC. Lastly, we apply these methods to two datasets (Hemidactylus and Galagidae) and find inconsistencies with the published literature. Our methods are powerful for rapid quantitative testing of species boundaries in large multilocus datasets and are implemented as an open source BEAST 2 package called SPEEDEMON.

Introducing SPEEDEMON, a package for BEAST 2 that better defines species boundaries based on molecular data demonstrated on gecko and loris datasets.

Phenotypic and genomic differences between biomes of the South America marsh rat, Holochilus brasiliensis

Abiotic factors can influence genetic and phenotypic divergence in several ways, and identifying the mechanisms responsible for generating this variation is challenging. However, when evaluated in combination, ecological characteristics and genetic and phenotypic information can help us to understand how habitat preferences can influence morphological and genetic patterns exhibited by taxa distributed between distinct biomes, such as the Atlantic Forest and Pampas biomes in South America. By combining distributional, environmental, phenotypic and genomic information from a habitat-specialist semi-aquatic rodent (Holochilus brasiliensis), we quantified the relationship between ecological niche differences and the phenotypic and genetic variation. The results demonstrate notable segregation among the ecological niches of H. brasiliensis within each biome, although we could not refute the hypothesis of niche similarity or equivalency. Such differences are consistent with a solid morphometric variation associated with the size of these rodents. However, the ecological and morphometric differentiation is not accompanied by the same pattern of genetic variation. Despite differences in the connectivity patterns in both biomes, the genetic differences corroborate a consistent level of migration history between biomes. Additionally, the association tests show that the environment explains a small and non- significant part of the genetic variation but a significant portion of the morphometric variation.

The living heart: Climate gradients predict desert mountain endemism

Mountain regions are centers of biodiversity endemism at a global scale but the role of arid-zone mountain ranges in shaping biodiversity patterns is poorly understood. Focusing on three guilds of taxa from a desert upland refugium in Australia, we sought to determine: (a) the relative extent to which climate, terrain or geological substrate predict endemism, and (b) whether patterns of endemism are complimentary across broad taxonomic guilds. We mapped regional endemism for plants, land snails, and vertebrates using combined Species Distribution Models (SDMs) for all endemic taxa (n = 82). We then modelled predictors of endemism using Generalised Additive Models (GAMs) and geology, terrain, and climate variables. We tested for the presence of inter- and intraguild hotspots of endemism. Many individual plant and land snail taxa were tightly linked with geology, corresponding to small distributions. Conversely, most vertebrate taxa were not constrained to specific geological substrates and occurred over larger areas. However, across all three guilds climate was the strongest predictor of regional endemism, particularly for plants wherein discrete hotspots of endemism were buffered from extreme summer temperatures. Land snail and vertebrate endemism peaked in areas with highest precipitation in the driest times of the year. Hotspots of endemism within each guild poorly predicted endemism in other guilds. We found an overarching signal that climatic gradients play a dominant role in the persistence of endemic taxa in an arid-zone mountain range system. An association with higher rainfall and cooler temperatures indicates that continuing trends toward hotter and drier climates may lead to range contractions in this, and potentially other, arid-zone mountain biotas. Contrasting patterns of endemism across guilds highlight the need to couple comprehensive regional planning for the protection of climate refugia, with targeted management of more localized and habitat specialist taxa.

Diversification of the mygalomorph spider genus Aname (Araneae: Anamidae) across the Australian arid zone: Tracing the evolution and biogeography of a continent-wide radiation

The assembly of the Australian arid zone biota has long fascinated biogeographers. Covering over two-thirds of the continent, Australia's vast arid zone biome is home to a distinctive fauna and flora, including numerous lineages which have diversified since the Eocene. Tracing the origins and speciation history of these arid zone taxa has been an ongoing endeavour since the advent of molecular phylogenetics, and an increasing number of studies on invertebrate animals are beginning to complement a rich history of research on vertebrate and plant taxa. In this study, we apply continent-wide genetic sampling and one of the largest phylogenetic data matrices yet assembled for a genus of Australian spiders, to reconstruct the phylogeny and biogeographic history of the open-holed trapdoor spider genus Aname L. Koch, 1873. This highly diverse lineage of Australian mygalomorph spiders has a distribution covering the majority of Australia west of the Great Dividing Range, but apparently excluding the high rainfall zones of eastern Australia and Tasmania. Original and legacy sequences were obtained for three mtDNA and four nuDNA markers from 174 taxa in seven genera, including 150 Aname specimen terminals belonging to 102 species-level operational taxonomic units, sampled from 32 bioregions across Australia. Reconstruction of the phylogeny and biogeographic history of Aname revealed three radiations (Tropical, Temperate-Eastern and Continental), which could be further broken into eight major inclusive clades. Ancestral area reconstruction revealed the Pilbara, Monsoon Tropics and Mid-West to be important ancestral areas for the genus Aname and its closest relatives, with the origin of Aname itself inferred in the Pilbara bioregion. From these origins in the arid north-west of Australia, our study found evidence for a series of subsequent biome transitions in separate lineages, with at least eight tertiary incursions back into the arid zone from more mesic tropical, temperate or eastern biomes, and only two major clades which experienced widespread (primary) in situ diversification within the arid zone. Based on our phylogenetic results, and results from independent legacy divergence dating studies, we further reveal the importance of climate-driven biotic change in the Miocene and Pliocene in shaping the distribution and composition of the Australian arid zone biota, and the value of continent-wide studies in revealing potentially complex patterns of arid zone diversification in dispersal-limited invertebrate taxa.

Crypsis and convergence: integrative taxonomic revision of the Gehyra australis group (Squamata: Gekkonidae) from northern Australia

For over two decades, assessments of geographic variation in mtDNA and small numbers of nuclear loci have revealed morphologically similar, but genetically divergent, intraspecific lineages in lizards from around the world. Subsequent morphological analyses often find subtle corresponding diagnostic characters to support the distinctiveness of lineages, but occasionally do not. In recent years it has become increasingly possible to survey geographic variation by sequencing thousands of loci, enabling more rigorous assessment of species boundaries across morphologically similar lineages. Here we take this approach, adding new, geographically extensive SNP data to existing mtDNA and exon capture datasets for the Gehyra australis and G. koira species complexes of gecko from northern Australia. The combination of exon-based phylogenetics with dense spatial sampling of mitochondrial DNA sequencing, SNP-based tests for introgression at lineage boundaries and newly-collected morphological evidence supports the recognition of nine species, six of which are newly described here. Detection of discrete genetic clusters using new SNP data was especially convincing where candidate taxa were continuously sampled across their distributions up to and across geographic boundaries with analyses revealing no admixture. Some species defined herein appear to be truly cryptic, showing little, if any, diagnostic morphological variation. As these SNP-based approaches are progressively applied, and with all due conservatism, we can expect to see a substantial improvement in our ability to delineate and name cryptic species, especially in taxa for which previous approaches have struggled to resolve taxonomic boundaries.

On and off the rocks: persistence and ecological diversification in a tropical Australian lizard radiation

Background

Congruent patterns in the distribution of biodiversity between regions or habitats suggest that key factors such as climatic and topographic variation may predictably shape evolutionary processes. In a number of tropical and arid biomes, genetic analyses are revealing deeper and more localised lineage diversity in rocky ranges than surrounding habitats. Two potential drivers of localised endemism in rocky areas are refugial persistence through climatic change, or ecological diversification and specialisation. Here we examine how patterns of lineage and phenotypic diversity differ across two broad habitat types (rocky ranges and open woodlands) in a small radiation of gecko lizards in the genus Gehyra (the australis group) from the Australian Monsoonal Tropics biome.

Results

Using a suite of approaches for delineating evolutionarily independent lineages, we find between 26 and 41 putative evolutionary units in the australis group (versus eight species currently recognised). Rocky ranges are home to a greater number of lineages that are also relatively more restricted in distribution, while lineages in open woodland habitats are fewer, more widely distributed, and, in one case, show evidence of range expansion. We infer at least two shifts out of rocky ranges and into surrounding woodlands. Phenotypic divergence between rocky ranges specialist and more generalist taxa is detected, but no convergent evolutionary regimes linked to ecology are inferred.

Conclusions

In climatically unstable biomes such as savannahs, rocky ranges have functioned as zones of persistence, generators of diversity and a source of colonists for surrounding areas. Phenotypic divergence can also be linked to the use of differing habitat types, however, the extent to which ecological specialisation is a primary driver or secondary outcome of localised diversification remains uncertain.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1408-1) contains supplementary material, which is available to authorized users.

Phylogenetic relationships and biogeographic history of the Australian trapdoor spider genus Conothele (Araneae: Mygalomorphae: Halonoproctidae): diversification into arid habitats in an otherwise tropical radiation

In Australia, climate change and continental drift have given rise to a complex biota comprising mesic specialists, arid-adapted lineages, and taxa that have arrived on the continent from Asia. We explore the phylogenetic diversity and biogeographic history of the Australian trapdoor spider genus Conothele Thorell, 1878 that is widespread in Australia’s monsoonal tropics and arid zone. We sequenced three mtDNA and five nuDNA markers from 224 specimens. We reconstructed the phylogenetic relationships among specimens and estimated the number of operational taxonomic units (OTUs) using species delimitation methods. The timing of divergences was estimated and ancestral area reconstructions were conducted. We recovered 61 OTUs, grouped into four major clades; a single clade represented by an arboreal ecomorph, and three fossorial clades. The Australian Conothele had a crown age of ~19 million years, and ancestral area reconstructions showed a complex history with multiple transitions among the monsoonal tropics, central arid zone, south-west and Pilbara bioregion. Conothele arrived on the continent during periods of biotic exchange with Asia. Since then, Conothele has colonised much of the Australian arid and monsoonal zones, during a period of climatic instability. The Pilbara bioregion harbours high lineage diversity, emphasising the role of climate refugia.

Conspicuously concealed: revision of the arid clade of the Gehyra variegata (Gekkonidae) group in Western Australia using an integrative molecular and morphological approach, with the description of five cryptic species

The methods used to detect and describe morphologically cryptic species have advanced in recent years, owing to the integrative nature of molecular and morphological techniques required to elucidate them. Here we integrate recent phylogenomic work that sequenced many genes but few individuals, with new data from mtDNA and morphology from hundreds of gecko specimens of the Gehyra variegata group from the Australian arid zone. To better understand morphological and geographical boundaries among cryptic forms, we generated new sequences from 656 Gehyra individuals, largely assigned to G. variegata group members over a wide area in Western Australia, with especially dense sampling in the Pilbara region, and combined them with 566 Gehyra sequences from GenBank, resulting in a dataset of 1,222 specimens. Results indicated the existence of several cryptic species, from new species with diagnostic morphological characters, to cases when there were no useful characters to discriminate among genetically distinctive species. In addition, the cryptic species often showed counter-intuitive distributions, including broad sympatry among some forms and short range endemism in other cases. Two new species were on long branches in the phylogram and restricted to the northern Pilbara region: most records of the moderately sized G. incognita sp. nov. are near the coast with isolated inland records, whereas the small-bodied saxicoline G. unguiculata sp. nov. is only known from a small area in the extreme north of the Pilbara. Three new species were on shorter branches in the phylogram and allied to G. montium. The moderately sized G. crypta sp. nov. occurs in the western and southern Pilbara and extends south through the Murchison region; this species was distinctive genetically, but with wide overlap of characters with its sister species, G. montium. Accordingly, we provide a table of diagnostic nucleotides for this species as well as for all other species treated here. Two small-bodied species occur in isolated coastal regions: G. capensis sp. nov. is restricted to the North West Cape and G. ocellata sp. nov. occurs on Barrow Island and other neighbouring islands. The latter species showed evidence of introgression with the mtDNA of G. crypta sp. nov., possibly due to recent connectivity with the mainland owing to fluctuating sea levels. However, G. ocellata sp. nov. was more closely related to G. capensis sp. nov. in the phylogenomic data and in morphology. Our study illustrates the benefits of combining phylogenomic data with extensive screens of mtDNA to identify large numbers of individuals to the correct cryptic species. This approach was able to provide sufficient samples with which to assess morphological variation. Furthermore, determination of geographic distributions of the new cryptic species should greatly assist with identification in the field, demonstrating the utility of sampling large numbers of specimens across wide areas.