Plio-Pleistocene diversification and biogeographic barriers in southern Australia reflected in the phylogeography of a widespread and common lizard species

A molecular phylogeny of the Petaluridae (Odonata: Anisoptera): A 160-Million-Year-Old story of drift and extinction

Petaluridae (Odonata: Anisoptera) is a relict dragonfly family, having diverged from its sister family in the Jurassic, of eleven species that are notable among odonates (dragonflies and damselflies) for their exclusive use of fen and bog habitats, their burrowing behavior as nymphs, large body size as adults, and extended lifespans. To date, several nodes within this family remain unresolved, limiting the study of the evolution of this peculiar family. Using an anchored hybrid enrichment dataset of over 900 loci we reconstructed the species tree of Petaluridae. To estimate the temporal origin of the genera within this family, we used a set of well-vetted fossils and a relaxed molecular clock model in a divergence time estimation analysis. We estimate that Petaluridae originated in the early Cretaceous and confirm the existence of monophyletic Gondwanan and Laurasian clades within the family. Our relaxed molecular clock analysis estimated that these clades diverged from their MRCA approximately 160 mya. Extant lineages within this family were identified to have persisted from 6 (Uropetala) to 120 million years (Phenes). Our biogeographical analyses focusing on a set of key regions suggest that divergence within Petaluridae is largely correlated with continental drift, the exposure of land bridges, and the development of mountain ranges. Our results support the hypothesis that species within Petaluridae have persisted for tens of millions of years, with little fossil evidence to suggest widespread extinction in the family, despite optimal conditions for the fossilization of nymphs. Petaluridae appear to be a rare example of habitat specialists that have persisted for tens of millions of years.

A giant armoured skink from Australia expands lizard morphospace and the scope of the Pleistocene extinctions

There are more species of lizards and snakes (squamates) alive today than any other order of land vertebrates, yet their fossil record has been poorly documented compared with other groups. Here, we describe a gigantic Pleistocene skink from Australia based on extensive material that includes much of the skull and postcranial skeleton, and spans ontogenetic stages from neonate to adult. Tiliqua frangens substantially expands the known ecomorphological diversity of squamates. At approximately 2.4 kg, it was more than double the mass of any living skink, with an exceptionally broad, deep skull, squat limbs and heavy, ornamented body armour. It probably filled the armoured herbivore niche that land tortoises (testudinids), absent from Australia, occupy on other continents. Tiliqua frangens and other giant Plio-Pleistocene skinks suggest that small-bodied groups that dominate vertebrate biodiversity might have lost their largest and often most morphologically extreme representatives in the Late Pleistocene, expanding the scope of these extinctions.

A forensically validated genetic toolkit for the species and lineage identification of the highly trafficked shingleback lizard (Tiliqua rugosa)

Shingleback lizards (Tiliqua rugosa) are among the most trafficked native fauna from Australia in the illegal pet trade. There are four morphologically recognised subspecies of shinglebacks, all with differing overseas market values. Shinglebacks from different geographic locales are often trafficked and housed together, which may complicate identifying the State jurisdiction where the poaching event occurred. Additionally, shinglebacks can be housed and trafficked with other species within the same genus, which may complicate DNA analysis, especially in scenarios where indirect evidence (e.g. swabs, faeces) is taken for analysis. In this study, a forensic genetic toolkit was designed and validated to target shingleback DNA for species identification and geographic origin. To do this, field sampling across Australia was conducted to expand the phylogeographic sampling of shinglebacks across their species range and include populations suspected to be poaching hotspots. A commonly used universal reptile primer set (ND4/LEU) was then validated for use in forensic casework related to the genus Tiliqua. Two additional ND4 primer sets were designed and validated. The first primer set was designed and demonstrated to preferentially amplify an ∼510 bp region of the genus Tiliqua over other reptiles and builds on existing data to expand the available phylogeographic database. The second primer set was designed and demonstrated to solely amplify an ∼220 bp region of T. rugosa ND4 over any other reptile species. Through the validation process, all primers were demonstrated to amplify T. rugosa DNA from a variety of sample types (e.g. degraded, low quality and mixed). Two of the primer sets were able to distinguish the genetic lineage of T. rugosa from the phylogeographic database. This work provides the first forensically validated toolkit and phylogeographic genetic database for Squatmate lizards.

Time-calibrated phylogeny and ecological niche models indicate Pliocene aridification drove intraspecific diversification of brushtail possums in Australia

Major aridification events in Australia during the Pliocene may have had significant impact on the distribution and structure of widespread species. To explore the potential impact of Pliocene and Pleistocene climate oscillations, we estimated the timing of population fragmentation and past connectivity of the currently isolated but morphologically similar subspecies of the widespread brushtail possum (Trichosurus vulpecula). We use ecological niche modeling (ENM) with the current fragmented distribution of brushtail possums to estimate the environmental envelope of this marsupial. We projected the ENM on models of past climatic conditions in Australia to infer the potential distribution of brushtail possums over 6 million years. D-loop haplotypes were used to describe population structure. From shotgun sequencing, we assembled whole mitochondrial DNA genomes and estimated the timing of intraspecific divergence. Our projections of ENMs suggest current possum populations were unlikely to have been in contact during the Pleistocene. Although lowered sea level during glacial periods enabled connection with habitat in Tasmania, climate fluctuation during this time would not have facilitated gene flow over much of Australia. The most recent common ancestor of sampled intraspecific diversity dates to the early Pliocene when continental aridification caused significant changes to Australian ecology and Trichosurus vulpecula distribution was likely fragmented. Phylogenetic analysis revealed that the subspecies T. v. hypoleucus (koomal; southwest), T. v. arnhemensis (langkurr; north), and T. v. vulpecula (bilda; southeast) correspond to distinct mitochondrial lineages. Despite little phenotypic differentiation, Trichosurus vulpecula populations probably experienced little gene flow with one another since the Pliocene, supporting the recognition of several subspecies and explaining their adaptations to the regional plant assemblages on which they feed.

Lineage diversity within a widespread endemic Australian skink to better inform conservation in response to regional-scale disturbance

Much attention is paid in conservation planning to the concept of a species, to ensure comparability across studies and regions when classifying taxa against criteria of endangerment and setting priorities for action. However, various jurisdictions now allow taxonomic ranks below the level of species and nontaxonomic intraspecific divisions to be factored into conservation planning-subspecies, key populations, evolutionarily significant units, or designatable units. Understanding patterns of genetic diversity and its distribution across the landscape is a key component in the identification of species boundaries and determination of substantial geographic structure within species. A total of 12,532 reliable polymorphic SNP loci were generated from 63 populations (286 individuals) covering the distribution of the Australian eastern three-lined skink, Bassiana duperreyi, to assess genetic population structure in the form of diagnosable lineages and their distribution across the landscape, with particular reference to the recent catastrophic bushfires of eastern Australia. Five well-supported diagnosable operational taxonomic units (OTUs) existed within B. duperreyi. Low levels of divergence of B. duperreyi between mainland Australia and Tasmania (no fixed allelic differences) support the notion of episodic exchange of alleles across Bass Strait (ca 60 m, 25 Kya) during periods of low sea level during the Upper Pleistocene rather than the much longer period of isolation (1.7 My) indicated by earlier studies using mitochondrial sequence variation. Our study provides foundational work for the detailed taxonomic re-evaluation of this species complex and the need for biodiversity assessment to include an examination of cryptic species and/or cryptic diversity below the level of species. Such information on lineage diversity within species and its distribution in the context of disturbance at a regional scale can be factored into conservation planning regardless of whether a decision is made to formally diagnose new species taxonomically and nomenclaturally.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Pleistocene divergence in the absence of gene flow among populations of a viviparous reptile with intraspecific variation in sex determination

Polymorphisms can lead to genetic isolation if there is differential mating success among conspecifics divergent for a trait. Polymorphism for sex-determining system may fall into this category, given strong selection for the production of viable males and females and the low success of heterogametic hybrids when sex chromosomes differ (Haldane's rule). Here we investigated whether populations exhibiting polymorphism for sex determination are genetically isolated, using the viviparous snow skink Carinascincus ocellatus. While a comparatively high elevation population has genotypic sex determination, in a lower elevation population there is an additional temperature component to sex determination. Based on 11,107 SNP markers, these populations appear genetically isolated. "Isolation with Migration" analysis also suggests these populations diverged in the absence of gene flow, across a period encompassing multiple Pleistocene glaciations and likely greater geographic proximity of populations. However, further experiments are required to establish whether genetic isolation may be a cause or consequence of differences in sex determination. Given the influence of temperature on sex in one lineage, we also discuss the implications for the persistence of this polymorphism under climate change.

Phylogeography of the iconic Australian pink cockatoo, Lophochroa leadbeateri

The pink cockatoo (Lophochroa leadbeateri; or Major Mitchell’s cockatoo) is one of Australia’s most iconic bird species. Two subspecies based on morphology are separated by a biogeographical divide, the Eyrean Barrier. Testing the genetic basis for this subspecies delineation, clarifying barriers to gene flow and identifying any cryptic genetic diversity will likely have important implications for conservation and management. Here, we used genome-wide single nucleotide polymorphisms (SNPs) and mitochondrial DNA data to conduct the first range-wide genetic assessment of the species. The aims were to investigate the phylogeography of the pink cockatoo, to characterize conservation units and to reassess subspecies boundaries. We found consistent but weak genetic structure between the two subspecies based on nuclear SNPs. However, phylogenetic analysis of nuclear SNPs and mitochondrial DNA sequence data did not recover reciprocally monophyletic groups, indicating incomplete evolutionary separation between the subspecies. Consequently, we have proposed that the two currently recognized subspecies be treated as separate management units rather than evolutionarily significant units. Given that poaching is suspected to be a threat to this species, we assessed the utility of our data for wildlife forensic applications. We demonstrated that a subspecies identification test could be designed using as few as 20 SNPs.

Araújo M, Kaefer IL, Lima AP. Hierarchical effects of historical and environmental factors on lizard assemblages in the upper Madeira River, Brazilian Amazonia

Investigating the role of historical and ecological factors structuring assemblages is relevant to understand mechanisms and processes affecting biodiversity across heterogeneous habitats. Considering that community assembly often involves scale-dependent processes, different spatial scales may reveal distinct factors structuring assemblages. In this study we use arboreal and leaf-litter lizard abundance data from 83 plots to investigate assemblage spatial structure at two distinct scales in southwestern Brazilian Amazonia. At a regional scale, we test the general hypothesis that the Madeira River acts as a barrier to dispersal of some lizard species, which results in distinct assemblages between river banks. At a local scale, we test the hypothesis that assemblages are not evenly distributed across heterogeneous habitats but respond to a continuum of inadequate-to-optimal portions of environmental predictors. Our results show that regional lizard assemblages are structured by the upper Madeira River acting as a regional barrier to 29.62% of the species sampled. This finding suggests species have been historically isolated at one of the river banks, or that distinct geomorphological features influence species occurrence at each river bank. At a local scale, different sets of environmental predictors affected assemblage composition between river banks or even along a river bank. These findings indicate that environmental filtering is a major cause of lizard assemblage spatial structure in the upper Madeira River, but predictor variables cannot be generalized over the extensive (nearly 500 km) study area. Based on a single study system we demonstrate that lizard assemblages along the forests near the banks of the upper Madeira River are not randomly structured but respond to multiple factors acting at different and hierarchical spatial scales.

Genomic data show little geographical structure across the naturally fragmented range of the purple-gaped honeyeater

Using single nucleotide polymorphisms and mitochondrial DNA sequences we find some evidence of genetic structure within a widespread and naturally fragmented species, the purple-gaped honeyeater (Lichenostomus cratitius), of southern Australian mallee shrublands. The very earliest stages of differentiation either side of the Nullarbor Barrier may already have been arrested by gene flow, some of which may have been anthropogenically induced.