Museum skins enable identification of introgression associated with cytonuclear discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale datasets and coalescent modelling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridisation. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least four operational taxonomic units (OTUs) across four species which diverged 183-278 kya. ABC modelling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Population genomic diversity and structure in the golden bandicoot: a history of isolation, extirpation, and conservation

Using genetic information to develop and implement conservation programs is vital for maintaining biodiversity and ecosystem resilience. Evaluation of the genetic variability within and among remnant populations can inform management of both natural and translocated populations to maximise species' adaptive potential, mitigate negative impacts of inbreeding, and subsequently minimise risk of extinction. Here we use reduced representation sequencing to undertake a genetic assessment of the golden bandicoot (Isoodon auratus), a threatened marsupial endemic to Australia. The currently recognised taxon consists of three subspecies distributed among multiple natural and translocated populations. After confirming the genetic distinctiveness of I. auratus from two closely related taxa, I. fusciventer and I. macrourus, we identified four genetic clusters within I. auratus. These clusters exhibited substantial genetic differentiation (pairwise FST values ranging from 0.18 to 0.65, pairwise DXY ranging from 0.1 to 0.168), reflecting long-term isolation of some populations on offshore islands and the influence of genetic drift. Mainland natural populations in the Kimberley region had the highest genetic diversity and the largest contribution to overall allelic and gene diversity compared to both natural and translocated island populations. A population translocated to Guluwuru Island in the Northern Territory had the lowest genetic diversity. Our data suggest that island populations can appear genetically unique due to genetic drift and this needs to be taken into account when considering genetic diversity in conservation efforts to maintain overall genetic diversity of the species. We effectively demonstrate how genomic information can guide practical conservation planning, especially when declining species are represented by multiple isolated populations.

Genetic mixing in conservation translocations increases diversity of a keystone threatened species, Bettongia lesueur

Translocation programmes are increasingly being informed by genetic data to monitor and enhance conservation outcomes for both natural and established populations. These data provide a window into contemporary patterns of genetic diversity, structure and relatedness that can guide managers in how to best source animals for their translocation programmes. The inclusion of historical samples, where possible, strengthens monitoring by allowing assessment of changes in genetic diversity over time and by providing a benchmark for future improvements in diversity via management practices. Here, we used reduced representation sequencing (ddRADseq) data to report on the current genetic health of three remnant and seven translocated boodie (Bettongia lesueur) populations, now extinct on the Australian mainland. In addition, we used exon capture data from seven historical mainland specimens and a subset of contemporary samples to compare pre-decline and current diversity. Both data sets showed the significant impact of population founder source (whether multiple or single) on the genetic diversity of translocated populations. Populations founded by animals from multiple sources showed significantly higher genetic diversity than the natural remnant and single-source translocation populations, and we show that by mixing the most divergent populations, exon capture heterozygosity was restored to levels close to that observed in pre-decline mainland samples. Relatedness estimates were surprisingly low across all contemporary populations and there was limited evidence of inbreeding. Our results show that a strategy of genetic mixing has led to successful conservation outcomes for the species in terms of increasing genetic diversity and provides strong rationale for mixing as a management strategy.

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

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

Linking life history to landscape for threatened species conservation in a multiuse region

Landscape-scale conservation that considers metapopulation dynamics will be essential for preventing declines of species facing multiple threats to their survival. Toward this end, we developed a novel approach that combines occurrence records, spatial-environmental data, and genetic information to model habitat, connectivity, and patterns of genetic structure and link spatial attributes to underlying ecological mechanisms. Using the threatened northern quoll (Dasyurus hallucatus) as a case study, we applied this approach to address the need for conservation decision-making tools that promote resilient metapopulations of this threatened species in the Pilbara, Western Australia, a multiuse landscape that is a hotspot for biodiversity and mining. Habitat and connectivity were predicted by different landscape characteristics. Whereas habitat suitability was overwhelmingly driven by terrain ruggedness, dispersal was facilitated by proximity to watercourses. Although there is limited evidence for major physical barriers in the Pilbara, areas with high silt and clay content (i.e., alluvial and hardpan plains) showed high resistance to dispersal. Climate subtlety shaped distributions and patterns of genetic turnover, suggesting the potential for local adaptation. By understanding these spatial-environmental associations and linking them to life-history and metapopulation dynamics, we highlight opportunities to provide targeted species management. To support this, we have created habitat, connectivity, and genetic uniqueness maps for conservation decision-making in the region. These tools have the potential to provide a more holistic approach to conservation in multiuse landscapes globally.

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.

Sequence Capture From Historical Museum Specimens: Maximizing Value for Population and Phylogenomic Studies

The application of high-throughput, short-read sequencing to degraded DNA has greatly increased the feasibility of generating genomic data from historical museum specimens. While many published studies report successful sequencing results from historical specimens; in reality, success and quality of sequence data can be highly variable. To examine predictors of sequencing quality, and methodological approaches to improving data accuracy, we generated and analyzed genomic sequence data from 115 historically collected museum specimens up to 180 years old. Data span both population genomic and phylogenomic scales, including historically collected specimens from 34 specimens of four species of Australian rock-wallabies (genus Petrogale) and 92 samples from 79 specimens of Australo-Papuan murine rodents (subfamily Murinae). For historical rodent specimens, where the focus was sampling for phylogenomics, we found that regardless of specimen age, DNA sequence libraries prepared from toe pad or bone subsamples performed significantly better than those taken from the skin (in terms of proportion of reads on target, number of loci captured, and data accuracy). In total, 93% of DNA libraries from toe pad or bone subsamples resulted in reliable data for phylogenetic inference, compared to 63% of skin subsamples. For skin subsamples, proportion of reads on target weakly correlated with collection year. Then using population genomic data from rock-wallaby skins as a test case, we found substantial improvement in final data quality by mapping to a high-quality “closest sister” de novo assembly from fresh tissues, compared to mapping to a sample-specific historical de novo assembly. Choice of mapping approach also affected final estimates of the number of segregating sites and Watterson's θ, both important parameters for population genomic inference. The incorporation of accurate and reliable sequence data from historical specimens has important outcomes for evolutionary studies at both population and phylogenomic scales. By assessing the outcomes of different approaches to specimen subsampling, library preparation and bioinformatic processing, our results provide a framework for increasing sequencing success for irreplaceable historical specimens.

Composition of a chemical signalling trait varies with phylogeny and precipitation across an Australian lizard radiation

The environment presents challenges to the transmission and detection of animal signalling systems, resulting in selective pressures that can drive signal divergence amongst populations in disparate environments. For chemical signals, climate is a potentially important selective force because factors such as temperature and moisture influence the persistence and detection of chemicals. We investigated an Australian lizard radiation (Heteronotia) to explore relationships between a sexually dimorphic chemical signalling trait (epidermal pore secretions) and two key climate variables: temperature and precipitation. We reconstructed the phylogeny of Heteronotia with exon capture phylogenomics, estimated phylogenetic signal in amongst-lineage chemical variation and assessed how chemical composition relates to temperature and precipitation using multivariate phylogenetic regressions. High estimates of phylogenetic signal indicate that the composition of epidermal pore secretions varies amongst lineages in a manner consistent with Brownian motion, although there are deviations to this, with stark divergences coinciding with two phylogenetic splits. Accounting for phylogenetic non-independence, we found that amongst-lineage chemical variation is associated with geographic variation in precipitation but not temperature. This contrasts somewhat with previous lizard studies, which have generally found an association between temperature and chemical composition. Our results suggest that geographic variation in precipitation can affect the evolution of chemical signalling traits, possibly influencing patterns of divergence amongst lineages and species.

Temperature predicts the rate of molecular evolution in Australian Eugongylinae skinks

Temperature differences over time and space have been hypothesized to cause variation in the rate of molecular evolution of species, but empirical evidence is mixed. To further test this hypothesis, we utilized a large exon-capture sequence data of Australian Eugongylinae skinks, exemplifying a radiation of temperature-sensitive ectotherms spanning a large latitudinal gradient. The association between temperature (and other species traits) and long-term substitution rate was assessed based on 1268 sequenced exons of 44 species pairs from the Eugongylinae subfamily using regression analyses. Temperature is the strongest, positively correlated predictor of variation in substitution rate across the Australian Eugongylinae. It explains 45% of variation in synonymous substitution rate, and 11% after controlling for all the other factors. Synonymous substitution rate is also negatively associated with body size, with a 6% variation explained by body size after controlling for the effects of temperature. Other factors are not associated with synonymous substitution rate after controlling for temperature. Overall, this study points to temperature as a strong predictor of the molecular evolution rate in the Eugongylinae subfamily, and demonstrates the power of large-scale exonic data to identify correlates of the rate of molecular evolution.

Environmental resistance and habitat quality influence dispersal of the saltwater crocodile

Landscape genetics commonly focuses on the effects of environmental resistance on animal dispersal patterns, but there is an emerging focus on testing environmental effects on emigration and settlement choices. In this study, we used landscape genetics approaches to quantify dispersal patterns in the world's largest crocodilian, the saltwater crocodile (Crocodylus porosus), and demonstrated environmental influences on three processes that comprise dispersal: emigration, movement and settlement. We found that both environmental resistance and properties of the source and destination catchments (proportion of breeding habitat) were important factors influencing observed dispersal events. Our habitat quality variables related to hypotheses about resource competition and represented the ratio of breeding habitat (which limits carrying capacity), suggesting that competition for habitat influences emigration and settlement choices, together with the strong effect of environmental resistance to movement (where high-quality habitat was associated with greatest environmental permeability). Approximately 42% of crocodiles were migrants from populations other than their sampling locations and some outstandingly productive populations had a much higher proportion of emigration rather than immigration. The distance most commonly travelled between source and destination was 150-200 km although a few travelled much longer distances, up to 600-700 km. Given the extensive dispersal range, individual catchments or hydrographic regions that combine two or three adjacent catchments are an appropriate scale for population management.

Limited introgression between rock-wallabies with extensive chromosomal rearrangements

Chromosome rearrangements can result in the rapid evolution of hybrid incompatibilities. Robertsonian fusions, particularly those with monobrachial homology, can drive reproductive isolation amongst recently diverged taxa. The recent radiation of rock-wallabies (genus Petrogale) is an important model to explore the role of Robertsonian fusions in speciation. Here, we pursue that goal using an extensive sampling of populations and genomes of Petrogale from north-eastern Australia. In contrast to previous assessments using mitochondrial DNA or nuclear microsatellite loci, genomic data are able to separate the most closely related species and to resolve their divergence histories. Both phylogenetic and population genetic analyses indicate introgression between two species that differ by a single Robertsonian fusion. Based on the available data, there is also evidence for introgression between two species which share complex chromosomal rearrangements. However, the remaining results show no consistent signature of introgression amongst species pairs and where evident, indicate generally low introgression overall. X-linked loci have elevated divergence compared with autosomal loci indicating a potential role for genic evolution to produce reproductive isolation in concert with chromosome change. Our results highlight the value of genome scale data in evaluating the role of Robertsonian fusions and structural variation in divergence, speciation, and patterns of molecular evolution.

Museum Genomics

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Predictors of phylogeographic structure among codistributed taxa across the complex Australian monsoonal tropics

Differences in the geographic scale and depth of phylogeographic structure across codistributed taxa can reveal how microevolutionary processes such as population isolation and persistence drive diversification. In turn, environmental heterogeneity, species' traits, and historical biogeographic barriers may influence the potential for isolation and persistence. Using extensive SNP data and a combination of population genetic summary statistics and landscape genomic analyses, we explored predictors of the scale and depth of phylogeographic structure in codistributed lizard taxa from the topographically and climatically complex monsoonal tropics (AMT) of Australia. We first resolved intraspecific lineages and then tested whether genetic divergence across space within lineages is related to isolation by distance, resistance and/or environment and whether these factors differ across genera or between rock-related versus habitat generalist taxa. We then tested whether microevolutionary processes within lineages explain differences in the geographic scale and depth of intraspecific phylogeographic lineages. The results indicated that landscape predictors of phylogeographic structure differ between taxa. Within lineages, there was prevalent isolation by distance, but the strength of isolation by distance is independent of the taxonomic family, habitat specialization, and climate. Isolation by environment is the strongest predictor of landscape-scale genetic divergence for all taxa, with both temperature and precipitation acting as limiting factors. The strength of isolation by distance does not predict the geographic scale of the phylogeographic structure. However, more localized lineages had higher mean individual heterozygosity and less negative Tajima's D. This result implies that finer-scale phylogeographic structuring within species is associated with larger and more stable populations and, hence, persistence.

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.

Geographic restriction, genetic divergence, and morphological disparity in the Brazilian Atlantic Forests: Insights from Leposoma lizards (Gymnophthalmidae, Squamata)

Lineage differentiation, long-term persistence, and range limitation promote high levels of phylogenetic and phylogeographic endemisms and likely underlie the abundant morphologically cryptic diversity observed in the Brazilian Atlantic Forests (AF). We explore lineage differentiation and range restriction in the AF and ask if genetic divergence and morphological disparity are correlated by integrating coalescent-based species delimitation, molecular phylogenetic, and morphological analyses in the lizard genus Leposoma. We present the first species tree for Leposoma and of their tribe, the Ecpleopodini. The analyses are based on the largest dataset ever assembled for Leposoma in terms of number of species (all represented), genetic markers (12 loci), and geographic coverage (~2,500 km). The exercise allows us to robustly delimit species within the genus and phylogeographic lineages within all species. We find support for the monophyly of the genus and for the recognition of a yet undescribed species around the Baía de Todos-os-Santos, in the state of Bahia; this form is distinct from all other congeners, both genetically and morphologically. We find that L. baturitensis, from the northeastern state of Ceará, is basal to the genus - and sister to a clade of six species restricted to the AF across the eastern coast of Brazil. Relationships within this coastal clade are ((((L. annectans, Leposoma sp.), L. scincoides), L. puk) (L. nanodactylus, L. sinepollex)). Phylogenetic and phylogeographic analyses, together with precise distribution data, allowed us to update the ranges of species and phylogeographic lineages. We reveal pervasive geographic restriction of divergent lineages in Leposoma at and below species level and discuss how forest refuges and rivers might have contributed to it. We find that morphological disparity lags behind genetic divergence in the genus because although they are correlated, the first accumulates at a much slower rate than the latter. We hope to encourage new studies in the area of AF north of the Doce river; phylogeographic sampling in that region has been much less common relative to southern sites, yet it may hold the key to several important processes defining biodiversity patterns in eastern Brazil. This appears to specially apply to processes underlying geographic restriction of morphologically cryptic, yet genetic divergent lineages, as the case of Leposoma.

Comparative assessment of range-wide patterns of genetic diversity and structure with SNPs and microsatellites: A case study with Iberian amphibians

Reduced representation genome sequencing has popularized the application of single nucleotide polymorphisms (SNPs) to address evolutionary and conservation questions in nonmodel organisms. Patterns of genetic structure and diversity based on SNPs often diverge from those obtained with microsatellites to different degrees, but few studies have explicitly compared their performance under similar sampling regimes in a shared analytical framework. We compared range‐wide patterns of genetic structure and diversity in two amphibians endemic to the Iberian Peninsula: Hyla molleri and Pelobates cultripes, based on microsatellite (18 and 14 loci) and SNP (15,412 and 33,140 loci) datasets of comparable sample size and spatial extent. Model‐based clustering analyses with STRUCTURE revealed minor differences in genetic structure between marker types, but inconsistent values of the optimal number of populations (K) inferred. SNPs yielded more repeatable and less admixed ancestries with increasing K compared to microsatellites. Genetic diversity was weakly correlated between marker types, with SNPs providing a better representation of southern refugia and of gradients of genetic diversity congruent with the demographic history of both species. Our results suggest that the larger number of loci in a SNP dataset can provide more reliable inferences of patterns of genetic structure and diversity than a typical microsatellite dataset, at least at the spatial and temporal scales investigated.

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates

The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci- and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.

Presence-only and Presence-absence Data for Comparing Species Distribution Modeling Methods

Species distribution models (SDMs) are widely used to predict and study distributions of species. Many different modeling methods and associated algorithms are used and continue to emerge. It is important to understand how different approaches perform, particularly when applied to species occurrence records that were not gathered in struc­tured surveys (e.g. opportunistic records). This need motivated a large-scale, collaborative effort, published in 2006, that aimed to create objective comparisons of algorithm performance. As a benchmark, and to facilitate future comparisons of approaches, here we publish that dataset: point location records for 226 anonymized species from six regions of the world, with accompanying predictor variables in raster (grid) and point formats. A particularly interesting characteristic of this dataset is that independent presence-absence survey data are available for evaluation alongside the presence-only species occurrence data intended for modeling. The dataset is available on Open Science Framework and as an R package and can be used as a benchmark for modeling approaches and for testing new ways to evaluate the accuracy of SDMs.

A Divergent Articulavirus in an Australian Gecko Identified Using Meta-Transcriptomics and Protein Structure Comparisons

The discovery of highly divergent RNA viruses is compromised by their limited sequence similarity to known viruses. Evolutionary information obtained from protein structural modelling offers a powerful approach to detect distantly related viruses based on the conservation of tertiary structures in key proteins such as the RNA-dependent RNA polymerase (RdRp). We utilised a template-based approach for protein structure prediction from amino acid sequences to identify distant evolutionary relationships among viruses detected in meta-transcriptomic sequencing data from Australian wildlife. The best predicted protein structural model was compared with the results of similarity searches against protein databases. Using this combination of meta-transcriptomics and protein structure prediction we identified the RdRp (PB1) gene segment of a divergent negative-sense RNA virus, denoted Lauta virus (LTAV), in a native Australian gecko (Gehyra lauta). The presence of this virus was confirmed by PCR and Sanger sequencing. Phylogenetic analysis revealed that Lauta virus likely represents a newly described genus within the family Amnoonviridae, order Articulavirales, that is most closely related to the fish virus Tilapia tilapinevirus (TiLV). These findings provide important insights into the evolution of negative-sense RNA viruses and structural conservation of the viral replicase among members of the order Articulavirales.

Colonic volvulus after laparoscopic left colectomy

Volvulus of the mobilized colon after laparoscopic left colectomy is rare. Contributing factors seem to be excessive length of the mobilized colon, absence of peritonization and absence of adhesions due to laparoscopy. Onset of colonic volvulus after laparoscopic left colectomy should lead to routine computerized tomography (CT), searching for an image suggestive of small intestinal incarceration under the neo-mesocolon, which might be an additional risk factor. Treatment consists of disincarceration of the small intestines while closing the mesocolic defect remains a subject of controversy.

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

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

Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.

Thermal physiological traits in tropical lowland amphibians: Vulnerability to climate warming and cooling

Climate change is affecting biodiversity and ecosystem function worldwide, and the lowland tropics are of special concern because organisms living in this region experience temperatures that are close to their upper thermal limits. However, it remains unclear how and whether tropical lowland species will be able to cope with the predicted pace of climate warming. Additionally, there is growing interest in examining how quickly thermal physiological traits have evolved across taxa, and whether thermal physiological traits are evolutionarily conserved or labile. We measured critical thermal maximum (CTmax) and minimum (CTmin) in 56 species of lowland Amazonian frogs to determine the extent of phylogenetic conservatism in tolerance to heat and cold, and to predict species' vulnerability to climate change. The species we studied live in sympatry and represent ~65% of the known alpha diversity at our study site. Given that critical thermal limits may have evolved differently in response to different temperature constraints, we tested whether CTmax and CTmin exhibit different rates of evolutionary change. Measuring both critical thermal traits allowed us to estimate species' thermal breadth and infer their potential to respond to abrupt changes in temperature (warming and cooling). Additionally, we assessed the contribution of life history traits and found that both critical thermal traits were correlated with species' body size and microhabitat use. Specifically, small direct-developing frogs in the Strabomantidae family appear to be at highest risk of thermal stress while tree frogs (Hylidae) and narrow mouthed frogs (Microhylidae) tolerate higher temperatures. While CTmax and CTmin had considerable variation within and among families, both critical thermal traits exhibited similar rates of evolutionary change. Our results suggest that 4% of lowland rainforest frogs assessed will experience temperatures exceeding their CTmax, 25% might be moderately affected and 70% are unlikely to experience pronounced heat stress under a hypothetical 3°C temperature increase.

Adjusting to climate: Acclimation, adaptation and developmental plasticity in physiological traits of a tropical rainforest lizard

The impact of climate change may be felt most keenly by tropical ectotherms. In these taxa, it is argued, thermal specialization means a given shift in temperature will have a larger effect on fitness. For species with limited dispersal ability, the impact of climate change depends on the capacity for their climate-relevant traits to shift. Such shifts can occur through genetic adaptation, various forms of plasticity, or a combination of these processes. Here we assess the extent and causes of shifts in 7 physiological traits in a tropical lizard, the rainforest sunskink (Lampropholis coggeri). Two populations were sampled that differ from each other in both climate and physiological traits. We compared trait values in each animal soon after field collection versus following acclimation to laboratory conditions. We also compared trait values between populations in: (i) recently field-collected animals; (ii) the same animals following laboratory acclimation; and (iii) the laboratory-reared offspring of these animals. Our results reveal high trait lability, driven primarily by acclimation and local adaptation. By contrast, developmental plasticity, resulting from incubation temperature, had little to no effect on most traits. These results suggest that, while specialized, tropical ectotherms may be capable of rapid shifts in climate-relevant traits.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Temporal genomic contrasts reveal rapid evolutionary responses in an alpine mammal during recent climate change

Many species have experienced dramatic changes in their abundance and distribution during recent climate change, but it is often unclear whether such ecological responses are accompanied by evolutionary change. We used targeted exon sequencing of 294 museum specimens (160 historic, 134 modern) to generate independent temporal genomic contrasts spanning a century of climate change (1911-2012) for two co-distributed chipmunk species: an endemic alpine specialist (Tamias alpinus) undergoing severe range contraction and a stable mid-elevation species (T. speciosus). Using a novel analytical approach, we reconstructed the demographic histories of these populations and tested for evidence of recent positive directional selection. Only the retracting species showed substantial population genetic fragmentation through time and this was coupled with positive selection and substantial shifts in allele frequencies at a gene, Alox15, involved in regulation of inflammation and response to hypoxia. However, these rapid population and gene-level responses were not detected in an analogous temporal contrast from another area where T. alpinus has also undergone severe range contraction. Collectively, these results highlight that evolutionary responses may be variable and context dependent across populations, even when they show seemingly synchronous ecological shifts. Our results demonstrate that temporal genomic contrasts can be used to detect very recent evolutionary responses within and among contemporary populations, even in the face of complex demographic changes. Given the wealth of specimens archived in natural history museums, comparative analyses of temporal population genomic data have the potential to improve our understanding of recent and ongoing evolutionary responses to rapidly changing environments.

Current geography masks dynamic history of gene flow during speciation in northern Australian birds

Genome divergence is greatly influenced by gene flow during early stages of speciation. As populations differentiate, geographic barriers can constrain gene flow and so affect the dynamics of divergence and speciation. Current geography, specifically disjunction and continuity of ranges, is often used to predict the historical gene flow during the divergence process. We test this prediction in eight meliphagoid bird species complexes codistributed in four regions. These regions are separated by known biogeographical barriers across northern Australia and Papua New Guinea. We find that bird populations currently separated by terrestrial habitat barriers within Australia and marine barriers between Australia and Papua New Guinea have a range of divergence levels and probability of gene flow not associated with current range connectivity. Instead, geographic distance and historical range connectivity better predict divergence and probability of gene flow. In this dynamic environmental context, we also find support for a nonlinear decrease of the probability of gene flow during the divergence process. The probability of gene flow initially decreases gradually after a certain level of divergence is reached. Its decrease then accelerates until the probability is close to zero. This implies that although geographic connectivity may have more of an effect early in speciation, other factors associated with higher divergence may play a more important role in influencing gene flow midway through and later in speciation. Current geographic connectivity may then mislead inferences regarding potential for gene flow during speciation under a complex and dynamic history of geographic and reproductive isolation.

Diversification across biomes in a continental lizard radiation

Ecological opportunity is a powerful driver of evolutionary diversification, and predicts rapid lineage and phenotypic diversification following colonization of competitor-free habitats. Alternatively, topographic or environmental heterogeneity could be key to generating and sustaining diversity. We explore these hypotheses in a widespread lineage of Australian lizards: the Gehyra variegata group. This clade occurs across two biomes: the Australian monsoonal tropics (AMT), where it overlaps a separate, larger bodied clade of Gehyra and is largely restricted to rocks; and in the larger Australian arid zone (AAZ) where it has no congeners and occupies trees and rocks. New phylogenomic data and coalescent analyses of AAZ taxa resolve lineages and their relationships and reveal high diversity in the western AAZ (Pilbara region). The AMT and AAZ radiations represent separate radiations with no difference in speciation rates. Most taxa occur on rocks, with small geographic ranges relative to widespread generalist taxa across the vast central AAZ. Rock-dwelling and generalist taxa differ morphologically, but only the lineage-poor central AAZ taxa have accelerated evolution. This accords with increasing evidence that lineage and morphological diversity are poorly correlated, and suggests environmental heterogeneity and refugial dynamics have been more important than ecological release in elevating lineage diversity.

Heritability of climate-relevant traits in a rainforest skink

There is justified concern about the impact of global warming on the persistence of tropical ectotherms. There is also growing evidence for strong selection on climate-relevant physiological traits. Understanding the evolutionary potential of populations is especially important for low dispersal organisms in isolated populations, because these populations have little choice but to adapt. Despite this, direct estimates of heritability and genetic correlations for physiological traits in ectotherms-which will determine their evolutionary responses to selection-are sparse, especially for reptiles. Here we examine the heritabilities and genetic correlations for a set of four morphological and six climate-relevant physiological traits in an isolated population of an Australian rainforest lizard, Lampropholis coggeri. These traits show considerable variation across populations in this species, suggesting local adaptation. From laboratory crosses, we estimated very low to moderate heritability of temperature-related physiological traits (h² < 0.31), but significant and higher heritability of desiccation resistance (h²~0.42). These values contrasted with uniformly higher heritabilities (h² > 0.51) for morphological traits. At the phenotypic level, there were positive associations among the morphological traits and between thermal limits. Growth rate was positively correlated with thermal limits, but there was no indication that morphology and physiology were linked in any other way. We found some support for a specialist-generalist trade-off in the thermal performance curve, but otherwise there was no evidence for evolutionary constraints, suggesting broadly labile multivariate trait structure. Our results indicate little potential to respond to selection on thermal traits in this population and provide new insights into the capacity of tropical ectotherms to adapt in situ to rapid climate change.

Insights into Platypus Population Structure and History from Whole-Genome Sequencing

The platypus is an egg-laying mammal which, alongside the echidna, occupies a unique place in the mammalian phylogenetic tree. Despite widespread interest in its unusual biology, little is known about its population structure or recent evolutionary history. To provide new insights into the dispersal and demographic history of this iconic species, we sequenced the genomes of 57 platypuses from across the whole species range in eastern mainland Australia and Tasmania. Using a highly improved reference genome, we called over 6.7 M SNPs, providing an informative genetic data set for population analyses. Our results show very strong population structure in the platypus, with our sampling locations corresponding to discrete groupings between which there is no evidence for recent gene flow. Genome-wide data allowed us to establish that 28 of the 57 sampled individuals had at least a third-degree relative among other samples from the same river, often taken at different times. Taking advantage of a sampled family quartet, we estimated the de novo mutation rate in the platypus at 7.0 × 10-9/bp/generation (95% CI 4.1 × 10-9-1.2 × 10-8/bp/generation). We estimated effective population sizes of ancestral populations and haplotype sharing between current groupings, and found evidence for bottlenecks and long-term population decline in multiple regions, and early divergence between populations in different regions. This study demonstrates the power of whole-genome sequencing for studying natural populations of an evolutionarily important species.

Phylogenomics of a rapid radiation: the Australian rainbow skinks

Background

The application of target capture with next-generation sequencing now enables phylogenomic analyses of rapidly radiating clades of species. But such analyses are complicated by extensive incomplete lineage sorting, demanding the use of methods that consider this process explicitly, such as the multispecies coalescent (MSC) model. However, the MSC makes strong assumptions about divergence history and population structure, and when using the full Bayesian implementation, current computational limits mean that relatively few loci and samples can be analysed for even modest sized radiations. We explore these issues through analyses of an extensive (> 1000 loci) dataset for the Australian rainbow skinks. This group consists of 3 genera and 41 described species, which likely diversified rapidly in Australia during the mid-late Miocene to occupy rainforest, woodland, and rocky habitats with corresponding diversity of morphology and breeding colouration. Previous phylogenetic analyses of this group have revealed short inter-nodes and high discordance among loci, limiting the resolution of inferred trees. A further complication is that many species have deep phylogeographic structure – this poses the question of how to sample individuals within species for analyses using the MSC.

Results

Phylogenies obtained using concatenation and summary coalescent species tree approaches to the full dataset are well resolved with generally consistent topology, including for previously intractable relationships near the base of the clade. As expected, branch lengths at the tips are substantially overestimated using concatenation. Comparisons of different strategies for sampling haplotypes for full Bayesian MSC analyses (for one clade and using smaller sets of loci) revealed, unexpectedly, that combining haplotypes across divergent phylogeographic lineages yielded consistent species trees.

Conclusions

This study of more than 1000 loci provides a strongly-supported estimate of the phylogeny of the Australian rainbow skinks, which will inform future research on the evolution and taxonomy of this group. Our analyses suggest that species tree estimation with the MSC can be quite robust to violation of the assumption that the individuals representing a taxon are sampled from a panmictic population.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1130-4) contains supplementary material, which is available to authorized users.

Accounting for Uncertainty in Gene Tree Estimation: Summary-Coalescent Species Tree Inference in a Challenging Radiation of Australian Lizards

Accurate gene tree inference is an important aspect of species tree estimation in a summary-coalescent framework. Yet, in empirical studies, inferred gene trees differ in accuracy due to stochastic variation in phylogenetic signal between targeted loci. Empiricists should, therefore, examine the consistency of species tree inference, while accounting for the observed heterogeneity in gene tree resolution of phylogenomic data sets. Here, we assess the impact of gene tree estimation error on summary-coalescent species tree inference by screening ${\sim}2000$ exonic loci based on gene tree resolution prior to phylogenetic inference. We focus on a phylogenetically challenging radiation of Australian lizards (genus Cryptoblepharus, Scincidae) and explore effects on topology and support. We identify a well-supported topology based on all loci and find that a relatively small number of high-resolution gene trees can be sufficient to converge on the same topology. Adding gene trees with decreasing resolution produced a generally consistent topology, and increased confidence for specific bipartitions that were poorly supported when using a small number of informative loci. This corroborates coalescent-based simulation studies that have highlighted the need for a large number of loci to confidently resolve challenging relationships and refutes the notion that low-resolution gene trees introduce phylogenetic noise. Further, our study also highlights the value of quantifying changes in nodal support across locus subsets of increasing size (but decreasing gene tree resolution). Such detailed analyses can reveal anomalous fluctuations in support at some nodes, suggesting the possibility of model violation. By characterizing the heterogeneity in phylogenetic signal among loci, we can account for uncertainty in gene tree estimation and assess its effect on the consistency of the species tree estimate. We suggest that the evaluation of gene tree resolution should be incorporated in the analysis of empirical phylogenomic data sets. This will ultimately increase our confidence in species tree estimation using summary-coalescent methods and enable us to exploit genomic data for phylogenetic inference. [Coalescence; concatenation; Cryptoblepharus; exon capture; gene tree; phylogenomics; species tree.].

Convergence across a continent: adaptive diversification in a recent radiation of Australian lizards

Recent radiations are important to evolutionary biologists, because they provide an opportunity to study the mechanisms that link micro- and macroevolution. The role of ecological speciation during adaptive radiation has been intensively studied, but radiations can arise from a diversity of evolutionary processes; in particular, on large continental landmasses where allopatric speciation might frequently precede ecological differentiation. It is therefore important to establish a phylogenetic and ecological framework for recent continental-scale radiations that are species-rich and ecologically diverse. Here, we use a genomic (approx. 1 200 loci, exon capture) approach to fit branch lengths on a summary-coalescent species tree and generate a time-calibrated phylogeny for a recent and ecologically diverse radiation of Australian scincid lizards; the genus Cryptoblepharus We then combine the phylogeny with a comprehensive phenotypic dataset for over 800 individuals across the 26 species, and use comparative methods to test whether habitat specialization can explain current patterns of phenotypic variation in ecologically relevant traits. We find significant differences in morphology between species that occur in distinct environments and convergence in ecomorphology with repeated habitat shifts across the continent. These results suggest that isolated analogous habitats have provided parallel ecological opportunity and have repeatedly promoted adaptive diversification. By contrast, speciation processes within the same habitat have resulted in distinct lineages with relatively limited morphological variation. Overall, our study illustrates how alternative diversification processes might have jointly stimulated species proliferation across the continent and generated a remarkably diverse group of Australian lizards.

Biodiversity analysis in the digital era

This paper explores what the virtual biodiversity e-infrastructure will look like as it takes advantage of advances in ‘Big Data’ biodiversity informatics and e-research infrastructure, which allow integration of various taxon-level data types (genome, morphology, distribution and species interactions) within a phylogenetic and environmental framework. By overcoming the data scaling problem in ecology, this integrative framework will provide richer information and fast learning to enable a deeper understanding of biodiversity evolution and dynamics in a rapidly changing world. The Atlas of Living Australia is used as one example of the advantages of progressing towards this future. Living in this future will require the adoption of new ways of integrating scientific knowledge into societal decision making.

This article is part of the themed issue ‘From DNA barcodes to biomes’.