Multiple origins of green blood in New Guinea lizards

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Revision of the montane New Guinean skink genus Lobulia (Squamata: Scincidae), with the description of four new genera and nine new species

The skink genus Lobulia is endemic to New Guinea, the largest and highest tropical island in the world. Lobulia and its related genera represent an important component of the montane herpetofauna of New Guinea, but it remains understudied and poorly known. We here provide the first, large-scale, systematic revision of Lobulia, using molecular phylogenetic and morphological comparisons to assess the monophyly of the genus and the diversity of species within it. We find that Lobulia, as currently defined, is polyphyletic. The eight species currently assigned to it form three clades. Furthermore, many specimens from New Guinea of unknown specific affinity are genetically and morphologically distinct from each other. Based on these data, we re-diagnose Lobulia and two of its closely related genera, Prasinohaema and Papuascincus. We erect four new genera (Alpinoscincus gen. nov., Nubeoscincus gen. nov., Ornithuroscincus gen. nov. and Palaia gen. nov.) to address the problem of polyphyly and describe nine new species Lobulia fortis sp. nov., Lobulia huonensis sp. nov., Loublia marmorata sp. nov., Lobulia vogelkopensis sp. nov., Ornithuroscincus bengaun sp. nov., Ornithuroscincus inornatus sp. nov., Ornithuroscincus pterophilus sp. nov., Ornithuroscincus shearmani sp. nov. and Ornithuroscincus viridis sp. nov. We supplement this taxonomic revision by investigating the biogeographic history of Lobulia s.l. and find evidence for a large radiation in the accreted terranes of New Guinea, with multiple independent colonizations of montane habitats and subsequent recolonization of lowland habitats. Our study reinforces the uniqueness and richness of the montane herpetofauna of New Guinea and the importance of mountains to biodiversity in the Tropics.

New Methods to Calculate Concordance Factors for Phylogenomic Datasets

We implement two measures for quantifying genealogical concordance in phylogenomic data sets: the gene concordance factor (gCF) and the novel site concordance factor (sCF). For every branch of a reference tree, gCF is defined as the percentage of "decisive" gene trees containing that branch. This measure is already in wide usage, but here we introduce a package that calculates it while accounting for variable taxon coverage among gene trees. sCF is a new measure defined as the percentage of decisive sites supporting a branch in the reference tree. gCF and sCF complement classical measures of branch support in phylogenetics by providing a full description of underlying disagreement among loci and sites. An easy to use implementation and tutorial is freely available in the IQ-TREE software package (http://www.iqtree.org/doc/Concordance-Factor, last accessed May 13, 2020).

Elevation is a stronger predictor of morphological trait divergence than competition in a radiation of tropical lizards

Adaptations for efficient performance are expected to shape animal morphology based on selection for microhabitat use and ecological forces. The presence of competitor species is predicted to cause niches to contract and enhance trait divergence. Therefore, increased species richness is expected to lead to greater trait divergence, and to result in reduced overlap and similarity between morphologies of sympatric species. We examined patterns of morphospace occupancy and partitioning in the skink fauna of New Guinea, the world's largest tropical island. Because skink species richness is largely decoupled from elevation in New Guinea, we could examine the effects of both factors (as proxies for competition and abiotic conditions), on morphospace occupancy and partitioning. We measured 1,860 specimens from 79 species of skinks throughout Papua New Guinea, and examined their morphospace occupancy in a spatial context. We calculated, for each assemblage within equal-area cells, the volume of morphospace occupied by all skinks, the mean volume occupied per species, and the mean distance and overlap between all species pairs. We then examined whether these metrics are related to species richness and elevation. Elevation is a stronger predictor of morphospace occupancy than species richness. As elevation increases, intraspecific variation decreases and morphologies become more similar to each other such that overall morphospace occupancy decreases. Highland skinks are, on average, smaller, thinner and shorter limbed than lowland species. We hypothesise that harsh climates in the New Guinea highland habitats impose strong selection on skinks to occupy specific areas of morphospace that facilitate efficient thermoregulation in suboptimal thermal conditions. We conclude that the effect of competition on trait divergence on a community and assemblage scale is eclipsed by abiotic selection pressures in these harsh environments.

Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression

The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.

Multiple origins of green coloration in frogs mediated by a novel biliverdin-binding serpin

Many vertebrates have distinctive blue-green bones and other tissues due to unusually high biliverdin concentrations-a phenomenon called chlorosis. Despite its prevalence, the biochemical basis, biology, and evolution of chlorosis are poorly understood. In this study, we show that the occurrence of high biliverdin in anurans (frogs and toads) has evolved multiple times during their evolutionary history, and relies on the same mechanism-the presence of a class of serpin family proteins that bind biliverdin. Using a diverse combination of techniques, we purified these serpins from several species of nonmodel treefrogs and developed a pipeline that allowed us to assemble their complete amino acid and nucleotide sequences. The described proteins, hereafter named biliverdin-binding serpins (BBS), have absorption spectra that mimic those of phytochromes and bacteriophytochromes. Our models showed that physiological concentration of BBSs fine-tune the color of the animals, providing the physiological basis for crypsis in green foliage even under near-infrared light. Additionally, we found that these BBSs are most similar to human glycoprotein alpha-1-antitrypsin, but with a remarkable functional diversification. Our results present molecular and functional evidence of recurrent evolution of chlorosis, describe a biliverdin-binding protein in vertebrates, and introduce a function for a member of the serpin superfamily, the largest and most ubiquitous group of protease inhibitors.

Cryptic diversity and non-adaptive radiation of montane New Guinea skinks (Papuascincus; Scincidae)

New Guinea, the world's largest and highest tropical island, has a rich but poorly known biota. Papuascincus is a genus of skinks endemic to New Guinea's mountain regions, comprising two wide-ranging species and two species known only from their type series. The phylogeny of the genus has never been examined and the relationships among its species - as well as between it and closely related taxa - are hitherto unknown. We performed the first large-scale molecular-phylogenetic study of Papuascincus, including sampling across the genus' range in Papua New Guinea. We sequenced three mitochondrial and two nuclear markers from 65 specimens of Papuascincus and reconstructed their phylogenetic relationships. We also performed species-delimitation analyses, estimated divergence times and ancestral biogeography, and examined body-size evolution within the genus. Papuascincus was strongly supported as monophyletic. It began radiating during the mid-Miocene in the area now comprising the Central Cordillera of New Guinea, then dispersed eastward colonising the Papuan Peninsula. We found evidence of extensive cryptic diversity within the genus, with between nine and 20 supported genetic lineages. These were estimated using three methods of species delimitation and predominantly occur in allopatry. Distribution and body-size divergence patterns indicated that character displacement in size took place during the evolutionary history of Papuascincus. We conclude that the genus requires comprehensive taxonomic revision and likely represents a species-rich lineage of montane skinks.

Visible light biophotosensors using biliverdin from Antheraea yamamai

We report an endogenous photoelectric biomolecule and demonstrate that such a biomolecule can be used to detect visible light. We identify the green pigment abundantly present in natural silk cocoons of Antheraea yamamai (Japanese oak silkmoth) as biliverdin, using mass spectroscopy and optical spectroscopy. Biliverdin extracted from the green silk cocoons generates photocurrent upon light illumination with distinct colors. We further characterize the basic performance, responsiveness, and stability of the biliverdin-based biophotosensors at a photovoltaic device level using blue, green, orange, and red light illumination. Biliverdin could potentially serve as an optoelectric biomolecule toward the development of next-generation implantable photosensors and artificial photoreceptors.