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.

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.

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.

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.

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.].

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).

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.

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.

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.].

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 Wallace’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.

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 Wallace’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.

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.

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.

Barcoding utility in a mega-diverse, cross-continental genus: keeping pace with Cyrtodactylus geckos.

Over the past decade, DNA barcoding has become a staple of low-cost molecular systematic investigations. The availability of universal primers and subsidized sequencing projects (PolarBOL, SharkBOL, SpongeBOL) have driven this popularity, often without appropriate investigation into the utility of barcoding data for the taxonomic group of interest. Here, our primary aim is to determine the phylogenetic value of DNA barcoding (mitochondrial locus COI) within the gecko genus Cyrtodactylus. With >40 new species described since last systematic investigation, Cyrtodactylus represents one of the most diverse extant squamate genera, and their contemporary distribution spans the Indian subcontinent, eastward through Indochina, and into AustraloPapua. The complex biogeographic history of this group, and morphology-only designation of many species have complicated our phylogenetic understanding of Cyrtodactylus. To highlight the need for continued inclusive molecular assessment, we use Vietnamese Cyrtodactylus as a case study showing the geopolitically paraphyletic nature of their history. We compare COI to the legacy marker ND2, and discuss the value of COI as an interspecific marker, as well as its shortcomings at deeper evolutionary scales. We draw attention back to the Cold Code as a subsidized method for incorporating molecular methods into species descriptions in the effort to maintain accurate phylogenies.

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.