Ancient drainages divide cryptic species in Australia’s arid zone: Morphological and multi-gene evidence for four new species of Beaked Geckos (Rhynchoedura)

Recent speciation and adaptation to aridity in the ecologically diverse Pilbara region of Australia enabled the native tobaccos (Nicotiana; Solanaceae) to colonize all Australian deserts

Over the last 6 million years, the arid Australian Eremaean Zone (EZ) has remained as dry as it is today. A widely accepted hypothesis suggests that the flora and fauna of arid regions were more broadly distributed before aridification began. In Australia, this process started around 20 million years ago (Ma), leading to gradual speciation as the climate became increasingly arid. Here, we use genomic data to investigate the biogeography and timing of divergence of native allotetraploid tobaccos, Nicotiana section Suaveolentes (Solanaceae). The original allotetraploid migrants from South America were adapted to mesic areas of Australia and recently radiated in the EZ, including in sandy dune fields (only 1.2 Ma old), after developing drought adaptations. Coalescent and maximum likelihood analyses suggest that Nicotiana section Suaveolentes arrived on the continent around 6 Ma, with the ancestors of the Pilbara (Western Australian) lineages radiating there at the onset of extreme aridity 5 Ma by locally adapting to these various ancient, highly stable habitats. The Pilbara thus served as both a mesic refugium and cradle for adaptations to harsher conditions, due to its high topographical diversity, providing microhabitats with varying moisture levels and its proximity to the ocean, which buffers against extreme aridity. This enabled species like Nicotiana to survive in mesic refugia and subsequently adapt to more arid conditions. These results demonstrate that initially poorly adapted plant groups can develop novel adaptations in situ, permitting extensive and rapid dispersal despite the highly variable and unpredictable extreme conditions of the EZ.

Revision of zigzag geckos (Diplodactylidae: Amalosia) in eastern Australia, with description of five new species

Geckos of the genus Amalosia Wells & Wellington, 1984 occur across eastern and northern Australia. Only five species are described but additional diversity has been recognised for some time. Here we assess species diversity in eastern Australia, using morphological and genetic (ND4 mtDNA) data. We describe five new species, all morphologically distinct and highly genetically distinct (>25% divergence). Amalosia hinesi sp. nov. is found in woodlands on the western side of the Great Dividing Range in south-east Queensland and north-east New South Wales. Amalosia saxicola sp. nov. is a large species found on rocks in the MackayTownsville areas of mid-east Queensland, including on many offshore islands. Amalosia nebula sp. nov. is restricted to rocky areas in upland sclerophyll forest of the Wet Tropics region of north-east Queensland. Amalosia capensis sp. nov. is a small species found in the northern half of Cape York Peninsula. Amalosia queenslandia sp. nov. is a small species that is widespread through woodlands over much of eastern and central Queensland. These species are diagnosed from other Amalosia species in eastern Australia, including A. cf. rhombifer which occurs in north-western Queensland. Amalosia cf. rhombifer is part of the clades comprising the remainder of the A. rhombifer complex across the Northern Territory and Western Australia, which will be dealt with separately. Herein, we also we also synonymise the monotypic genus Nebulifera with Amalosia. This revision brings the number of Amalosia species to ten.

New Guinea uplift opens ecological opportunity across a continent

Sahul unites the world's largest and highest tropical island and the oldest and most arid continent on the backdrop of dynamic environmental conditions. Massive geological uplift in New Guinea is predicted to have acted as a species pump from the late Miocene onward, but the impact of this process on biogeography and diversification remains untested across Sahul as a whole. To address this, we reconstruct the assembly of a recent and diverse radiation of rodents (Murinae: Hydromyini) spanning New Guinea, Australia, and oceanic islands. Using phylogenomic data from 270 specimens, including many recently extinct and highly elusive species, we find that the orogeny and expansion of New Guinea opened ecological opportunity and triggered diversification across a continent. After a single over-water colonization from Asia ca. 8.5 Ma, ancestral Hydromyini were restricted to the tropical rainforest of proto-New Guinea for 3.5 million years. Following a shift in diversification coincident with the orogeny of New Guinea ca. 5 Ma and subsequent colonization of Australia, transitions between geographic regions (n = 24) and biomes (n = 34) become frequent. Recurrent over-water colonization between mainland and islands demonstrate how islands can play a substantial role in the assembly of continental fauna. Our results are consistent with a model of increased ecological opportunity across Sahul following major geological uplift in New Guinea ca. 5 Ma, with sustained diversification facilitated by over-water colonization from the Pleistocene to present. We show how geological processes, biome transitions, and over-water colonization collectively drove the diversification of an expansive continental radiation.

Persistence of a Geographically-Stable Hybrid Zone in Puerto Rican Dwarf Geckos

Determining the mechanisms that create and maintain biodiversity is a central question in ecology and evolution. Speciation is the process that creates biodiversity. Speciation is mediated by incompatibilities that lead to reproductive isolation between divergent populations and these incompatibilities can be observed in hybrid zones. Gecko lizards are a speciose clade possessing an impressive diversity of behavioral and morphological traits. In geckos, however, our understanding of the speciation process is negligible. To address this gap, we used genetic sequence data (both mitochondrial and nuclear markers) to revisit a putative hybrid zone between Sphaerodactylus nicholsi and Sphaerodactylus townsendi in Puerto Rico, initially described in 1984. First, we addressed discrepancies in the literature on the validity of both species. Second, we sampled a 10-km-wide transect across the putative hybrid zone and tested explicit predictions about its dynamics using cline models. Third, we investigated potential causes for the hybrid zone using species distribution modeling and simulations; namely, whether unique climatic variables within the hybrid zone might elicit selection for intermediate phenotypes. We find strong support for the species-level status of each species and no evidence of movement, or unique climatic variables near the hybrid zone. We suggest that this narrow hybrid zone is geographically stable and is maintained by a combination of dispersal and selection. Thus, this work has identified an extant model system within geckos that that can be used for future investigations detailing genetic mechanisms of reproductive isolation in an understudied vertebrate group.

Using high-throughput sequencing to investigate the factors structuring genomic variation of a Mediterranean grasshopper of great conservation concern

Inferring the demographic history of species is fundamental for understanding their responses to past climate/landscape alterations and improving our predictions about the future impacts of the different components of ongoing global change. Estimating the time-frame at which population fragmentation took place is also critical to determine whether such process was shaped by ancient events (e.g. past climate/geological changes) or if, conversely, it was driven by recent human activities (e.g. habitat loss). We employed genomic data (ddRAD-Seq) to determine the factors shaping contemporary patterns of genetic variation in the Iberian cross-backed grasshopper Dociostaurus crassiusculus, an endangered species with limited dispersal capacity and narrow habitat requirements. Our analyses indicate the presence of two ancient lineages and three genetic clusters resulted from historical processes of population fragmentation (~18-126 ka) that predate the Anthropocene. Landscape genetic analyses indicate that the limits of major river basins are the main geographical feature explaining large-scale patterns of genomic differentiation, with no apparent effect of human-driven habitat fragmentation. Overall, our study highlights the importance of detailed phylogeographic, demographic and spatially-explicit landscape analyses to identify evolutionary significant units and determine the relative impact of historical vs. anthropogenic factors on processes of genetic fragmentation in taxa of great conservation concern.

How the Aridification of Australia Structured the Biogeography and Influenced the Diversification of a Large Lineage of Australian Cicadas

Over the last 30 million years, Australia's landscape has undergone dramatic cooling and drying due to the establishment of the Antarctic Circumpolar Current and change in global CO$_{2}$ levels. Studies have shown that many Australian organisms went extinct during these major cooling events, while others experienced adaptive radiations and increases in diversification rates as a result of exploiting new niches in the arid zone. Despite the many studies on diversification and biogeography in Australia, few have been continent-wide and none have focused on a group of organisms adapted to feeding on plants. We studied 162 species of cicadas in the Australian Pauropsalta complex, a large generic lineage within the tribe Cicadettini. We asked whether there were changes in the diversification rate of Pauropsalta over time and if so: 1) which clades were associated with the rate change? 2) did timing of rate shifts correspond to known periods of dramatic historical climate change, 3) did increases in diversification rate along select lineages correspond to adaptive radiations with movement into the arid zone? To address these questions, we estimated a molecular phylogeny of the Pauropsalta complex using ${\sim}$5300 bp of nucleotide sequence data distributed among five loci (one mtDNA locus and four nDNA loci). We found that this large group of cicadas did not diversify at a constant rate as they spread through Australia; instead the signature of decreasing diversification rate changed roughly around the time of the expansion of the east Antarctic ice sheets ${\sim}$16 Ma and the glaciation of the northern hemisphere ${\sim}$3 Ma. Unlike other Australian taxa, the Pauropsalta complex did not explosively radiate in response to an early invasion of the arid zone. Instead multiple groups invaded the arid zone and experienced rates of diversification similar to mesic-distributed taxa. We found evidence for relictual groups, located in pre-Mesozoic habitat, that have not diversified and continue to reside on mesic hosts in isolated "habitat islands". Future work should focus on groups of similar ages with similar distribution patterns to determine whether this tempo and pattern of diversification and biogeography is consistent with evidence from other phytophagous insects.

Mass turnover and recovery dynamics of a diverse Australian continental radiation

Trends in global and local climate history have been linked to observed macroevolutionary patterns across a variety of organisms. These climatic pressures may unilaterally or asymmetrically influence the evolutionary trajectory of clades. To test and compare signatures of changing global (Eocene-Oligocene boundary cooling) and continental (Miocene aridification) environments on a continental fauna, we investigated the macroevolutionary dynamics of one of Australia's most diverse endemic radiations, pygopodoid geckos. We generated a time-calibrated phylogeny (>90% taxon coverage) to test whether (i) asymmetrical pygopodoid tree shape may be the result of mass turnover deep in the group's history, and (ii) how Miocene aridification shaped trends in biome assemblages. We find evidence of mass turnover in pygopodoids following the isolation of the Australian continental plate ∼30 million years ago, and in contrast, gradual aridification is linked to elevated speciation rates in the young arid zone. Surprisingly, our results suggest that invasion of arid habitats was not an evolutionary end point. Instead, arid Australia has acted as a source for diversity, with repeated outward dispersals having facilitated diversification of this group. This pattern contrasts trends in richness and distribution of other Australian vertebrates, illustrating the profound effects historical biome changes have on macroevolutionary patterns.

Towards a higher-level Ensifera phylogeny inferred from mitogenome sequences

Although mitogenomes are useful tools for inferring evolutionary history, only a few representative ones can be used for most Ensifera lineages. Thirty-two ensiferan mitogenomes were determined using ABI Sanger sequencing and standard primer walking of 2-3 overlapping Long-PCR fragments, or Illumina® HiSeq2000 for "shotgun" sequenced long-PCR-amplified mitochondrial or total genomic DNA. Six patterns of gene arrangements, including the novel trnR-trnS^AGN-trnA-trnN-trnG-nad3 in Lipotactes tripyrga (Lipotactinae), were identified from 59 ensiferan mitogenomes. The results suggest that trnM-trnI-trnQ and trnA-trnR-trnE-trnS^AGN-trnN-trnF rearrangements might be a shared derived character in Pseudophyllinae and Gryllidae, respectively. We found base composition biases in our dataset, which potentially complicate the inference of higher-level ensiferan phylogeny. Site-heterogeneous Bayesian inference (BI) and site-homogeneous maximum likelihood (ML) analyses recovered all ensiferan superfamilies as monophyletic. The site-homogeneous BI analysis failed to recover the monophyly of Stenopelmatoidea. As Schizodactyloidea was only represented by Comicus campestris, its monophyly could not be tested. In the Triassic/Jurassic boundary, Ensifera diverged into grylloid and non-grylloid clades. All analyses confirmed Grylloidea and Gryllotalpoidea as sister groups. Site-heterogeneous BI analysis found Schizodactyloidea as the most basal lineage and sister to the clade formed by Grylloidea and Gryllotalpoidea, but the site-homogeneous analyses placed it basally to the non-grylloid clade and recovered a sister relationship between Tettigonioidea and (Hagloidea, Rhaphidophoroidea, Stenopelmatoidea), although this clade had a low support. The site-heterogeneous BI analysis found Tettigonioidea and Hagloidea were sister groups (posterior probability (PP)=0.99), Stenopelmatoidea was sister to (Tettigonioidea, Hagloidea) (PP>0.91), and Rhaphidophoroidea was basal to the non-grylloid clade. At a lower level, all analyses divided Tettigonioidea into Phaneropteridae and Tettigoniidae.

Terrestrial mountain islands and Pleistocene climate fluctuations as motors for speciation: A case study on the genus Pseudovelia (Hemiptera: Veliidae)

This study investigated the influences of geographic isolation and climate fluctuation on the genetic diversity, speciation, and biogeography of the genus Pseudovelia (Hemiptera: Veliidae) in subtropical China and tropic Indo-China Peninsula. Species nucleotide and haplotype diversities decreased with reduction in species distribution limits. The gene tree was congruent with the taxonomy of monophyly, except for four species, P. contorta, P. extensa, P. tibialis tibialis, and P. vittiformis. The conflicts between the genes and species tree could be due to long-term isolation and incomplete lineage sorting. Diversification analysis showed that the diversification rate (0.08 sp/My shifted to 0.5 sp/My) changed at 2.1 Ma, which occurred in the early Pleistocene period. Ancestral area reconstruction suggested that subtropical species possibly evolved from the tropics region (i.e., Indo-China Peninsula). Results implied that narrow endemics harbored relatively low genetic diversity because of small effective population and genetic drift. Radiation of subtropical Pseudovelia species was rapidly promoted by Pleistocene climate fluctuations and geographic isolation. The acute rising of the Hengduan Mountain with the entire uplift of the Qinghai-Tibet Plateau induced the initial differentiation of Pseudovelia species. These results highlighted the importance of geographical isolation and climate changes in promoting speciation in mountain habitat islands.

Snails in the desert: Assessing the mitochondrial and morphological diversity and the influence of aestivation behavior on lineage differentiation in the Australian endemic Granulomelon Iredale, 1933 (Stylommatophora: Camaenidae)

Progressive aridification since the mid-Miocene has had a significant influence on the evolution of the biota in the arid zone of central Australia. Especially moisture sensitive groups, such as snails, are often restricted to topographically complex areas, which have acted as refugia in an otherwise inhospitable environment. This historical fragmentation is deemed to be a potent agent of allopatric lineage diversification. Camaenid land snails are amongst only a few terrestrial gastropods that have managed to survive in the arid zone probably due to their ability to escape desiccation through aestivation. Here, we present the first study of the mitochondrial lineage differentiation in an endemic land snail genus from the Australian 'Red Centre', Granulomelon Iredale, 1933. Exposing significant incongruence between mtDNA phylogeny and morphology-based taxonomy, we completely revise the species and genus level taxonomy of this camaenid group. We demonstrate that this genus contains three species, G. grandituberculatum, G. adcockianum and G. squamulosum, which have so far been assigned to different genera: Granulomelon Iredale, 1933 (junior synonym: Baccalena Iredale, 1937), Basedowena Iredale, 1937 and Pleuroxia Ancey, 1887. Two of these species are widespread comprising multiple divergent mitochondrial lineages. Based on a molecular clock estimate, these lineages diverged approximately during the mid-Pleistocene, a period of particularly severe aridification. The phylogeographic patterns are consistent with an isolation-by-distance model in one species but not the other. We suggest that these differences can be attributed to their distinctive aestivation behavior.

Independent Transitions between Monsoonal and Arid Biomes Revealed by Systematic Revison of a Complex of Australian Geckos (Diplodactylus; Diplodactylidae)

How the widespread expansion and intensification of aridity through the Neogene has shaped the Austral biota is a major question in Antipodean biogeography. Lineages distributed across wide aridity gradients provide opportunities to examine the timing, frequency, and direction of transitions between arid and mesic regions. Here, we use molecular genetics and morphological data to investigate the systematics and biogeography of a nominal Australian gecko species (Diplodactylus conspicillatus sensu lato) with a wide distribution spanning most of the Australian Arid Zone (AAZ) and Monsoonal Tropics (AMT). Our data support a minimum of seven genetically distinct and morphologically diagnosable taxa; we thus redefine the type species, ressurrect three names from synonymy, and describe three new species. Our inferred phylogeny suggests the history and diversification of lineages in the AAZ and AMT are intimately linked, with evidence of multiple independent interchanges since the late Miocene. However, despite this shared history, related lineages in these two regions also show evidence of broadly contrasting intra-regional responses to aridification; vicarance and speciation in older and increasingly attenuated mesic regions, versus a more dynamic history including independent colonisations and recent range expansions in the younger AAZ.

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra-a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson-Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.

Pleistocene climate change promoted rapid diversification of aquatic invertebrates in Southeast Australia

Background

The Pleistocene Ice Ages were the most recent geohistorical event of major global impact, but their consequences for most parts of the Southern hemisphere remain poorly known. We investigate a radiation of ten species of Sternopriscus, the most species-rich genus of epigean Australian diving beetles. These species are distinct based on genital morphology but cannot be distinguished readily by mtDNA and nDNA because of genotype sharing caused by incomplete lineage sorting. Their genetic similarity suggests a Pleistocene origin.

Results

We use a dataset of 3858 bp of mitochondrial and nuclear DNA to reconstruct a phylogeny of Sternopriscus using gene and species trees. Diversification analyses support the finding of a recent rapid speciation event with estimated speciation rates of up to 2.40 species per MY, which is considerably higher than the proposed average rate of 0.16 species per MY for insects. Additionally, we use ecological niche modeling and analyze data on habitat preferences to test for niche divergence between species of the recent Sternopriscus radiation. These analyses show that the species can be characterized by a set of ecological variables referring to habitat, climate and altitude.

Conclusions

Our results suggest that the repeated isolation of populations in glacial refugia might have led to divergent ecological adaptations and the fixation of morphological traits supporting reproductive isolation and therefore may have promoted speciation. The recent Sternopriscus radiation fulfills many characteristics of a species flock and would be the first described example of an aquatic insect species flock. We argue that the species of this group may represent a stage in speciation past the species flock condition because of their mostly broad and often non-overlapping ranges and preferences for different habitat types.