pipesnake: generalized software for the assembly and analysis of phylogenomic datasets from conserved genomic loci.

MOTIVATION: Phylogenetics has moved into the era of genomics, incorporating enormous volumes of data to study questions at both shallow and deep scales. With this increase in information, phylogeneticists need new tools and skills to manipulate and analyze these data. To facilitate these tasks and encourage reproducibility, the community is increasingly moving toward automated workflows. RESULTS: Here we present pipesnake, a phylogenomics pipeline written in Nextflow for the processing, assembly, and phylogenetic estimation of genomic data from short-read sequences. pipesnake is an easy to use and efficient software package designed for this next era in phylogenetics. AVAILABILITY AND IMPLEMENTATION: pipesnake is publicly available on GitHub at https://github.com/AusARG/pipesnake and accompanied by documentation and a wiki/tutorial.

Populating a Continent: Phylogenomics Reveal the Timing of Australian Frog Diversification.

The Australian continent's size and isolation make it an ideal place for studying the accumulation and evolution of biodiversity. Long separated from the ancient supercontinent Gondwana, most of Australia's plants and animals are unique and endemic, including the continent's frogs. Australian frogs comprise a remarkable ecological and morphological diversity categorized into a small number of distantly related radiations. We present a phylogenomic hypothesis based on an exon-capture dataset that spans the main clades of Australian myobatrachoid, pelodryadid hyloid, and microhylid frogs. Our time-calibrated phylogenomic-scale phylogeny identifies great disparity in the relative ages of these groups which vary from Gondwanan relics to recent immigrants from Asia and include arguably the continent's oldest living vertebrate radiation. This age stratification provides insight into the colonization of⁠, and diversification on, the Australian continent through deep time, during periods of dramatic climatic and community changes. Contemporary Australian frog diversity highlights the adaptive capacity of anurans, particularly in response to heat and aridity, and explains why they are one of the continent's most visible faunas.

A return-on-investment approach for prioritization of rigorous taxonomic research needed to inform responses to the biodiversity crisis.

Global biodiversity loss is a profound consequence of human activity. Disturbingly, biodiversity loss is greater than realized because of the unknown number of undocumented species. Conservation fundamentally relies on taxonomic recognition of species, but only a fraction of biodiversity is described. Here, we provide a new quantitative approach for prioritizing rigorous taxonomic research for conservation. We implement this approach in a highly diverse vertebrate group-Australian lizards and snakes. Of 870 species assessed, we identified 282 (32.4%) with taxonomic uncertainty, of which 17.6% likely comprise undescribed species of conservation concern. We identify 24 species in need of immediate taxonomic attention to facilitate conservation. Using a broadly applicable return-on-investment framework, we demonstrate the importance of prioritizing the fundamental work of identifying species before they are lost.

Polyploidy breaks speciation barriers in Australian burrowing frogs Neobatrachus.

Polyploidy has played an important role in evolution across the tree of life but it is still unclear how polyploid lineages may persist after their initial formation. While both common and well-studied in plants, polyploidy is rare in animals and generally less understood. The Australian burrowing frog genus Neobatrachus is comprised of six diploid and three polyploid species and offers a powerful animal polyploid model system. We generated exome-capture sequence data from 87 individuals representing all nine species of Neobatrachus to investigate species-level relationships, the origin and inheritance mode of polyploid species, and the population genomic effects of polyploidy on genus-wide demography. We describe rapid speciation of diploid Neobatrachus species and show that the three independently originated polyploid species have tetrasomic or mixed inheritance. We document higher genetic diversity in tetraploids, resulting from widespread gene flow between the tetraploids, asymmetric inter-ploidy gene flow directed from sympatric diploids to tetraploids, and isolation of diploid species from each other. We also constructed models of ecologically suitable areas for each species to investigate the impact of climate on differing ploidy levels. These models suggest substantial change in suitable areas compared to past climate, which correspond to population genomic estimates of demographic histories. We propose that Neobatrachus diploids may be suffering the early genomic impacts of climate-induced habitat loss, while tetraploids appear to be avoiding this fate, possibly due to widespread gene flow. Finally, we demonstrate that Neobatrachus is an attractive model to study the effects of ploidy on the evolution of adaptation in animals.

A Comprehensive Approach to Detect Hybridization Sheds Light on the Evolution of Earth's Largest Lizards.

Hybridization between species occurs more frequently in vertebrates than traditionally thought, but distinguishing ancient hybridization from other phenomena that generate similar evolutionary patterns remains challenging. Here, we used a comprehensive workflow to discover evidence of ancient hybridization between the Komodo dragon (Varanus komodoensis) from Indonesia and a common ancestor of an Australian group of monitor lizards known colloquially as sand monitors. Our data comprise $>$300 nuclear loci, mitochondrial genomes, phenotypic data, fossil and contemporary records, and past/present climatic data. We show that the four sand monitor species share more nuclear alleles with $V$. komodoensis than expected given a bifurcating phylogeny, likely as a result of hybridization between the latter species and a common ancestor of sand monitors. Sand monitors display phenotypes that are intermediate between their closest relatives and $V$. komodoensis. Biogeographic analyses suggest that $V$. komodoensis and ancestral sand monitors co-occurred in northern Australia. In agreement with the fossil record, this provides further evidence that the Komodo dragon once inhabited the Australian continent. Our study shows how different sources of evidence can be used to thoroughly characterize evolutionary histories that deviate from a treelike pattern, that hybridization can have long-lasting effects on phenotypes, and that detecting hybridization can improve our understanding of evolutionary and biogeographic patterns.[Biogeography; introgression; Komodo dragon; phylogenetic networks; phylogenomics; reticulation; Varanus.].

Phylogenomics of Monitor Lizards and the Role of Competition in Dictating Body Size Disparity.

Organismal interactions drive the accumulation of diversity by influencing species ranges, morphology, and behavior. Interactions vary from agonistic to cooperative and should result in predictable patterns in trait and range evolution. However, despite a conceptual understanding of these processes, they have been difficult to model, particularly on macroevolutionary timescales and across broad geographic spaces. Here, we investigate the influence of biotic interactions on trait evolution and community assembly in monitor lizards (Varanus). Monitors are an iconic radiation with a cosmopolitan distribution and the greatest size disparity of any living terrestrial vertebrate genus. Between the colossal Komodo dragon Varanus komodoensis and the smallest Australian dwarf goannas, Varanus length and mass vary by multiple orders of magnitude. To test the hypothesis that size variation in this genus was driven by character displacement, we extended existing phylogenetic comparative methods which consider lineage interactions to account for dynamic biogeographic history and apply these methods to Australian monitors and marsupial predators. Incorporating both exon-capture molecular and morphological data sets we use a combined evidence approach to estimate the relationships among living and extinct varaniform lizards. Our results suggest that communities of Australian Varanus show high functional diversity as a result of continent-wide interspecific competition among monitors but not with faunivorous marsupials. We demonstrate that patterns of trait evolution resulting from character displacement on continental scales are recoverable from comparative data and highlight that these macroevolutionary patterns may develop in parallel across widely distributed sympatric groups.[Character displacement; comparative methods; phylogenetics; trait evolution; Varanus.].

Evolutionary models demonstrate rapid and adaptive diversification of Australo-Papuan pythons.

Lineages may diversify when they encounter available ecological niches. Adaptive divergence by ecological opportunity often appears to follow the invasion of a new environment with open ecological space. This evolutionary process is hypothesized to explain the explosive diversification of numerous Australian vertebrate groups following the collision of the Eurasian and Australian plates 25 Mya. One of these groups is the pythons, which demonstrate their greatest phenotypic and ecological diversity in Australo-Papua (Australia and New Guinea). Here, using an updated and near complete time-calibrated phylogenomic hypothesis of the group, we show that following invasion of this region, pythons experienced a sudden burst of speciation rates coupled with multiple instances of accelerated phenotypic evolution in head and body shape and body size. These results are consistent with adaptive radiation theory with an initial rapid niche-filling phase and later slow-down approaching niche saturation. We discuss these findings in the context of other Australo-Papuan adaptive radiations and the importance of incorporating adaptive diversification systems that are not extraordinarily species-rich but ecomorphologically diverse to understand how biodiversity is generated.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace's Line.

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace's and Lyddeker's lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker's line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.].

Multilocus phylogeny of Bornean Bent-Toed geckos (Gekkonidae: Cyrtodactylus) reveals hidden diversity, taxonomic disarray, and novel biogeographic patterns.

The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.

Miocene biome turnover drove conservative body size evolution across Australian vertebrates.

On deep time scales, changing climatic trends can have a predictable influence on macroevolution. From evidence of mass extinctions, we know that rapid climatic oscillations can indirectly open niche space and precipitate adaptive radiation, changing the course of ecological diversification. These dramatic shifts in the global climate, however, are rare events relative to extended periods of protracted climate change and biome turnover. It remains unclear whether during gradually changing periods, shifting habitats may instead promote non-adaptive speciation by facilitating allopatry and phenotypic conservatism. Using fossil-calibrated, species-level phylogenies for five Australian radiations comprising more than 800 species, we investigated temporal trends in biogeography and body size evolution. Here, we demonstrate that gradual Miocene cooling and aridification correlates with the restricted phenotypic diversification of multiple ecologically diverse vertebrate groups. This probably occurred as species ranges became fractured and isolated during continental biome restructuring, encouraging a shift towards conservatism in body size evolution. Our results provide further evidence that abiotic changes, not only biotic interactions, may act as selective forces influencing phenotypic macroevolution.

Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma.

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.

Evolutionary bursts drive morphological novelty in the world's largest skinks.

Animal phenotypes evolve and diverge as a result of differing selective pressures and drift. These processes leave unique signatures in patterns of trait evolution, impacting the tempo and mode of morphological macroevolution. While there is a broad understanding of the history of some organismal traits (e.g., body size), there is little consensus about the evolutionary mode of most others. This includes the relative contribution of prolonged (Darwinian gradualist) and episodic (Simpsonian jump) changes toward the evolution of novel morphologies. Here, we use new exon-capture and linear morphological datasets to investigate the tempo and mode of morphological evolution in Australo-Melanesian Tiliquini skinks. We generate a well-supported time-calibrated phylogenomic tree from ∼400 nuclear markers for more than 100 specimens, including undescribed diversity, and provide unprecedented resolution of the rapid Miocene diversification of these lizards. By collecting a morphological dataset that encompasses the lizard body plan (19 traits across the head, body, limb, and tail), we are able to identify that most traits evolve conservatively, but infrequent evolutionary bursts result in morphological novelty. These phenotypic discontinuities occur via rapid rate increases along individual branches, inconsistent with both gradualistic and punctuated equilibrial evolutionary modes. Instead, this "punctuated gradualism" has resulted in the rapid evolution of blue-tongued giants and armored dwarves in the ∼20 million years since colonizing Australia. These results outline the evolutionary pathway toward new morphologies and highlight the heterogeneity of evolutionary tempo and mode, even within individual traits.

Evolution across the adaptive landscape in a hyperdiverse beetle radiation.

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

Competition and geography underlie speciation and morphological evolution in Indo-Australasian monitor lizards.

How biotic and abiotic factors act together to shape biological diversity is a major question in evolutionary biology. The recent availability of large datasets and development of new methodological approaches provide new tools to evaluate the predicted effects of ecological interactions and geography on lineage diversification and phenotypic evolution. Here, we use a near complete phylogenomic-scale phylogeny and a comprehensive morphological dataset comprising more than a thousand specimens to assess the role of biotic and abiotic processes in the diversification of monitor lizards (Varanidae). This charismatic group of lizards shows striking variation in species richness among its clades and multiple instances of endemic radiation in Indo-Australasia (i.e., the Indo-Australian Archipelago and Australia), one of Earth's most biogeographically complex regions. We found heterogeneity in diversification dynamics across the family. Idiosyncratic biotic and geographic conditions appear to have driven diversification and morphological evolution in three endemic Indo-Australasian radiations. Furthermore, incumbency effects partially explain patterns in the biotic exchange between Australia and New Guinea. Our results offer insight into the dynamic history of Indo-Australasia, the evolutionary significance of competition, and the long-term consequences of incumbency effects.

Ecomorphometric Analysis of Diversity in Cranial Shape of Pygopodid Geckos.

Pygopodids are elongate, functionally limbless geckos found throughout Australia. The clade presents low taxonomic diversity (∼45 spp.), but a variety of cranial morphologies, habitat use, and locomotor abilities that vary between and within genera. In order to assess potential relationships between cranial morphology and ecology, computed tomography scans of 29 species were used for 3D geometric morphometric analysis. A combination of 24 static landmarks and 20 sliding semi-landmarks were subjected to Generalized Procrustes Alignment. Disparity in cranial shape was visualized through Principal Component Analysis, and a multivariate analysis of variance (MANOVA) was used to test for an association between shape, habitat, and diet. A subset of 27 species with well-resolved phylogenetic relationships was used to generate a phylomorphospace and conduct phylogeny-corrected MANOVA. Similar analyses were done solely on Aprasia taxa to explore species-level variation. Most of the variation across pygopodids was described by principal component (PC) 1(54%: cranial roof width, parabasisphenoid, and occipital length), PC2 (12%: snout elongation and braincase width), and PC3 (6%: elongation and shape of the palate and rostrum). Without phylogenetic correction, both habitat and diet were significant influencers of variation in cranial morphology. However, in the phylogeny-corrected MANOVA, habitat remained weakly significant, but not diet, which can be explained by generic-level differences in ecology rather than among species. Our results demonstrate that at higher levels, phylogeny has a strong effect on morphology, but that influence may be due to small sample size when comparing genera. However, because some closely related taxa occupy distant regions of morphospace, diverging diets, and use of fossorial habitats may contribute to variation seen in these geckos.

What's under the hood? Phylogeny and taxonomy of the snake genera Parasuta Worrell and Suta Worrell (Squamata: Elapidae), with a description of a new species from the Pilbara, Western Australia.

Despite decades of phylogenetic studies, the generic and species-level relationships of some Australian elapid snakes remain problematic. The morphologically conservative genus Parasuta comprises small nocturnal snakes with a particularly obfuscated taxonomic history. Here we provide a molecular phylogenetic analysis of all currently recognised species including members of the sister genus Suta and provide new morphological data that lead to a taxonomic revision of generic and species boundaries. We failed to find support for monophyly of Parasuta or Suta, instead supporting previous evidence that these two genera should be combined. Our species-level investigations revise the boundaries between P. gouldii (Gray) and P. spectabilis (Krefft) resulting in recognition that both P. spectabilis bushi (Storr) and P. spectabilis nullarbor (Storr) are conspecific with P. gouldii. We also find the Pilbara population of P. monachus (Storr) to be specifically distinct. As a consequence of this information, we synonymise Parasuta with its senior synonym Suta, redescribe S. gouldii, S. monachus and S. spectabilis to clarify morphological and geographical boundaries and describe S. gaikhorstorum sp. nov., which differs from all other described Suta species, including the geographically proximate and similar-looking S. monachus, by a combination of molecular genetic markers, morphometric attributes, details of colouration and scalation. The recognition of S. gaikhorstorum sp. nov. adds to the growing list of the many endemic reptiles from this exceptionally diverse biotic region. We also designate a lectotype for S. spectabilis from the original syntype series, highlight a distinctive population from the Great Victoria Desert in Western Australia and comment on further unresolved issues regarding the relationships between S. dwyeri (Worrell) and S. nigriceps (Gȕnther).

How mountains shape biodiversity: The role of the Andes in biogeography, diversification, and reproductive biology in South America's most species-rich lizard radiation (Squamata: Liolaemidae).

Testing hypotheses on drivers of clade evolution and trait diversification provides insight into many aspects of evolutionary biology. Often, studies investigate only intrinsic biological properties of organisms as the causes of diversity, however, extrinsic properties of a clade's environment, particularly geological history, may also offer compelling explanations. The Andes are a young mountain chain known to have shaped many aspects of climate and diversity of South America. The Liolaemidae are a radiation of South American reptiles with over 300 species found across most biomes and with similar numbers of egg-laying and live-bearing species. Using the most complete dated phylogeny of the family, we tested the role of Andean uplift in biogeography, diversification patterns, and parity mode of the Liolaemidae. We find that the Andes promoted lineage diversification and acted as a species pump into surrounding biomes. We also find strong support for the role of Andean uplift in boosting the species diversity of these lizards via allopatric fragmentation. Finally, we find repeated shifts in parity mode associated with changing thermal niches, with live-bearing favored in cold climates and egg-laying favored in warm climates. Importantly, we find evidence for possible reversals to oviparity, an evolutionary transition believed to be extremely rare.

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.