Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

A Morrison stem gekkotan reveals gecko evolution and Jurassic biogeography

Geckos are a speciose and globally distributed clade of Squamata (lizards, including snakes and amphisbaenians) that are characterized by a host of modifications for nocturnal, scansorial and insectivorous ecologies. They are among the oldest divergences in the lizard crown, so understanding the origin of geckoes (Gekkota) is essential to understanding the origin of Squamata, the most species-rich extant tetrapod clade. However, the poor fossil record of gekkotans has obscured the sequence and timing of the assembly of their distinctive morphology. Here, we describe the first North American stem gekkotan based on a three-dimensionally preserved skull from the Morrison Formation of western North America. Despite its Late Jurassic age, the new species already possesses several key characteristics of the gekkotan skull along with retained ancestral features. We show that this new stem gekkotan, and several previously named species of uncertain phylogenetic relationships, comprise a widespread clade of early crown lizards, substantiating faunal homogeneity in Laurasia during the Late Jurassic that extended across disparate ecological, body-size and physiological classes.

Embryonic development of the skull in a parthenogenetic lizard, the mourning gecko (Lepidodactylus lugubris)

Gekkotans are one of the major clades of squamate reptiles. As one of the earliest-diverging lineages, they are crucial in studying deep-level squamate phylogeny and evolution. Developmental studies can shed light on the origin of many important morphological characters, yet our knowledge of cranial development in gekkotans is very incomplete. Here, we describe the embryonic development of the skull in a parthenogenetic gekkonid, the mourning gecko (Lepidodactylus lugubris), studied using non-acidic double staining and histological sectioning. Our analysis indicates that the pterygoid is the first ossifying bone in the skull, as in almost all other studied squamates, followed closely by the surangular and prearticular. The next to appear are the dentary, frontal, parietal and squamosal. The tooth-bearing upper jaw bones, the premaxilla and maxilla, develop relatively late. In contrast to previous reports, the premaxilla starts ossifying from two distinct centres, reminiscent of the condition observed in diplodactylids and eublepharids. Only a single ossification centre of the postorbitofrontal is observed. Some of the endochondral bones of the braincase (prootic, opisthotic, supraoccipital) and the dermal parasphenoid are the last bones to appear. The skull roof is relatively poorly ossified near the time of hatching, with a large frontoparietal fontanelle still present. Many bones begin ossifying relatively later in L. lugubris than in the phyllodactylid Tarentola annularis, which suggests that its ossification sequence is heterochronic with respect to T. annularis.

Evolutionary origins of the prolonged extant squamate radiation

Squamata is the most diverse clade of terrestrial vertebrates. Although the origin of pan-squamates lies in the Triassic, the oldest undisputed members of extant clades known from nearly complete, uncrushed material come from the Cretaceous. Here, we describe three-dimensionally preserved partial skulls of two new crown lizards from the Late Jurassic of North America. Both species are placed at the base of the skink, girdled, and night lizard clade Pan-Scincoidea, which consistently occupies a position deep inside the squamate crown in both morphological and molecular phylogenies. The new lizards show that several features uniting pan-scincoids with another major lizard clade, the pan-lacertoids, in trees using morphology were convergently acquired as predicted by molecular analyses. Further, the palate of one new lizard bears a handful of ancestral saurian characteristics lost in nearly all extant squamates, revealing an underappreciated degree of complex morphological evolution in the early squamate crown. We find strong evidence for close relationships between the two new species and Cretaceous taxa from Eurasia. Together, these results suggest that early crown squamates had a wide geographic distribution and experienced complicated morphological evolution even while the Rhynchocephalia, now solely represented by the tuatara, was the dominant clade of lepidosaurs.

The Jurassic rise of squamates as supported by lepidosaur disparity and evolutionary rates

The squamates (lizards, snakes, and relatives) today comprise more than 10,000 species, and yet their sister group, the Rhynchocephalia, is represented by a single species today, the tuatara. The explosion in squamate diversity has been tracked back to the Cretaceous Terrestrial Revolution, 100 million years ago (Ma), the time when flowering plants began their takeover of terrestrial ecosystems, associated with diversification of coevolving insects and insect-eating predators such as lizards, birds, and mammals. Squamates arose much earlier, but their long pre-Cretaceous history of some 150 million years (Myr) is documented by sparse fossils. Here, we provide evidence for an initial radiation of squamate morphology in the Middle and Late Jurassic (174-145 Ma), and show that they established their key ecological roles much earlier than had been assumed, and they have not changed them much since.

Testing Phylogenetic Stability with Variable Taxon Sampling

Over the past three decades, computational capabilities have grown at such a rapid rate that they have given rise to many computationally heavy science fields such as phylogenomics. As increasingly more genomes are sequenced in the three domains of life, larger and more species-complete phylogenetic tree reconstructions are leading to a better understanding of the tree of life and the evolutionary histories in deep times. However, these large datasets pose unique challenges from a modeling and computational perspective: accurately describing the evolutionary process of thousands of species is still beyond the capability of current models, while the computational burden limits our ability to test multiple hypotheses. Thus, it is common practice to reduce the size of a dataset by selecting species to represent a clade (taxon sampling). Unfortunately, this process is subjective, and comparisons of large tree of life studies show that choice and number of species used in a dataset can alter the topology obtained. Thus, taxon sampling is, in itself, a process that needs to be fully investigated to determine its effect on phylogenetic stability. Here, we present the theory and practical application of an automated pipeline that can be easily implemented to explore the effect of taxon sampling on phylogenetic reconstructions. The application of this approach was recently discussed in a study of Terrabacteria and shows its power in investigating the accuracy of deep nodes of a phylogeny.

Phenotypic characters of static homology increase phylogenetic stability under direct optimization of otherwise dynamic homology characters

Direct optimization of unaligned sequence characters provides a natural framework to explore the sensitivity of phylogenetic hypotheses to variation in analytical parameters. Phenotypic data, when combined into such analyses, are typically analyzed with static homology correspondences unlike the dynamic homology sequence data. Static homology characters may be expected to constrain the direct optimization and thus, potentially increase the similarity of phylogenetic hypotheses under different cost sets. However, whether a total-evidence approach increases the phylogenetic stability or not remains empirically largely unexplored. Here, I studied the impact of static homology data on sensitivity using six empirical data sets composed of several molecular markers and phenotypic data. The inclusion of static homology phenotypic data increased the average stability of phylogenetic hypothesis in five out of the six data sets. To investigate if any static homology characters would have similar effect, the analyses were repeated with randomized phenotypic data, and with one of the molecular markers fixed as static homology characters. These analyses had, on average, almost no effect on the phylogenetic stability, although the randomized phenotypic data sometimes resulted in even higher stability than empirical phenotypic data. The impact was related to the strength of the phylogenetic signal in the phenotypic data: higher average jackknife support of the phenotypic tree correlated with stronger stabilizing effect in the total-evidence analysis. Phenotypic data with a strong signal made the total-evidence trees topologically more similar to the phenotypic trees, thus, they constrained the dynamic homology correspondences of the sequence data. Characters that increase phylogenetic stability are particularly valuable for phylogenetic inference. These results indicate an important role and additive value of phenotypic data in increasing the stability of phylogenetic hypotheses in total-evidence analyses.

Total evidence and sensitivity phylogenetic analyses of egg-brooding frogs (Anura: Hemiphractidae)

We study the phylogenetic relationships of egg-brooding frogs, a group of 118 neotropical species, unique among anurans by having embryos with large bell-shaped gills and females carrying their eggs on the dorsum, exposed or inside a pouch. We assembled a total evidence dataset of published and newly generated data containing 51 phenotypic characters and DNA sequences of 20 loci for 143 hemiphractids and 127 outgroup terminals. We performed six analytical strategies combining different optimality criteria (parsimony and maximum likelihood), alignment methods (tree- and similarity-alignment), and three different indel coding schemes (fifth character state, unknown nucleotide, and presence/absence characters matrix). Furthermore, we analyzed a subset of the total evidence dataset to evaluate the impact of phenotypic characters on hemiphractid phylogenetic relationships. Our main results include: (i) monophyly of Hemiphractidae and its six genera for all our analyses, novel relationships among hemiphractid genera, and non-monophyly of Hemiphractinae according to our preferred phylogenetic hypothesis; (ii) non-monophyly of current supraspecific taxonomies of Gastrotheca, an updated taxonomy is provided; (iii) previous differences among studies were mainly caused by differences in analytical factors, not by differences in character/taxon sampling; (iv) optimality criteria, alignment method, and indel coding caused differences among optimal topologies, in that order of degree; (v) in most cases, parsimony analyses are more sensitive to the addition of phenotypic data than maximum likelihood analyses; (vi) adding phenotypic data resulted in an increase of shared clades for most analyses.

The Review of the Autotomy of Agamid Lizards with Considerations about the Types of Autotomy and Regeneration

We present a review of the data on the intervertebral autotomy and regeneration of agamid lizards based on an analysis of information obtained over a 35-year period after the publication of thorough reviews (Arnold, 1984, 1988 and Bellairs, Bryant, 1985). It is supplemented by our own studies of 869 specimens of agamid lizards (Sauria, Agamidae) stored in the herpetological collections of the Zoological Institute of the Russian Academy of Sciences (St. Petersburg, Russia) and the Zoological Museum of the Moscow State University (Moscow, Russia), represented by 31 species of 16 genera. The manifestations of the ability for autotomy and regeneration in phylogenetic lineages within the family—Leiolepidinae, Amphibolurinae, Agaminae, Draconinae—are considered. A comparative morphological analysis of the structure of the caudal vertebrae was carried out using the Computer Microtomography Methods (micro-CT) in the following ecomorphological types of agama: (1) with developed abilities to caudal autotomy and regeneration, (2) with the ability to caudal autotomy but without regeneration and (3) without the ability to autotomy. The phenomenon of intervertebral autotomy (urotomy) in snakes is considered too. Possible ways of evolution of the ability to caudal autotomy as a defense strategy against predators are discussed in the phylogenetic context.

The Making of Calibration Sausage Exemplified by Recalibrating the Transcriptomic Timetree of Jawed Vertebrates

Molecular divergence dating has the potential to overcome the incompleteness of the fossil record in inferring when cladogenetic events (splits, divergences) happened, but needs to be calibrated by the fossil record. Ideally but unrealistically, this would require practitioners to be specialists in molecular evolution, in the phylogeny and the fossil record of all sampled taxa, and in the chronostratigraphy of the sites the fossils were found in. Paleontologists have therefore tried to help by publishing compendia of recommended calibrations, and molecular biologists unfamiliar with the fossil record have made heavy use of such works (in addition to using scattered primary sources and copying from each other). Using a recent example of a large node-dated timetree inferred from molecular data, I reevaluate all 30 calibrations in detail, present the current state of knowledge on them with its various uncertainties, rerun the dating analysis, and conclude that calibration dates cannot be taken from published compendia or other secondary or tertiary sources without risking strong distortions to the results, because all such sources become outdated faster than they are published: 50 of the (primary) sources I cite to constrain calibrations were published in 2019, half of the total of 280 after mid-2016, and 90% after mid-2005. It follows that the present work cannot serve as such a compendium either; in the slightly longer term, it can only highlight known and overlooked problems. Future authors will need to solve each of these problems anew through a thorough search of the primary paleobiological and chronostratigraphic literature on each calibration date every time they infer a new timetree, and that literature is not optimized for that task, but largely has other objectives.

Embryology of the naso-palatal complex in Gekkota based on detailed 3D analysis in Lepidodactylus lugubris and Eublepharis macularius

The vomeronasal organ (VNO), nasal cavity, lacrimal duct, choanal groove, and associated parts of the superficial (soft tissue) palate are called the naso-palatal complex. Despite the morphological diversity of the squamate noses, little is known about the embryological basis of this variation. Moreover, developmental data might be especially interesting in light of the morpho-molecular discordance of squamate phylogeny, since a 'molecular scenario' implies an occurrence of unexpected scale of homoplasy also in olfactory systems. In this study, we used X-ray microtomography and light microscopy to describe morphogenesis of the naso-palatal complex in two gekkotans: Lepidodactylus lugubris (Gekkonidae) and Eublepharis macularius (Eublepharidae). Our embryological data confirmed recent findings about the nature of some developmental processes in squamates, for example, involvement of the lateral nasal prominence in the formation of the choanal groove. Moreover, our study revealed previously unknown differences between the studied gekkotans and allows us to propose redefinition of the anterior concha of Sphenodon. Interpretation of some described conditions might be problematic in the phylogenetic context, since they represent unknown: squamate, nonophidian squamate, or gekkotan features.

Insights into skull evolution in fossorial snakes, as revealed by the cranial morphology of Atractaspis irregularis (Serpentes: Colubroidea)

Comparative osteological analyses of extant organisms provide key insight into major evolutionary transitions and phylogenetic hypotheses. This is especially true for snakes, given their unique morphology relative to other squamates and the persistent controversy regarding their evolutionary origins. However, the osteology of several major snake groups remains undescribed, thus hindering efforts to accurately reconstruct the phylogeny of snakes. One such group is the Atractaspididae, a family of fossorial colubroids. We herein present the first detailed description of the atractaspidid skull, based on fully segmented micro-computed tomography (micro-CT) scans of Atractaspis irregularis. The skull of Atractaspis presents a highly unique morphology influenced by both fossoriality and paedomorphosis. This paedomorphosis is especially evident in the jaws, palate, and suspensorium, the major elements associated with macrostomy (large-gaped feeding in snakes). Comparison to scolecophidians-a group of blind, fossorial, miniaturized snakes-in turn sheds light on current hypotheses of snake phylogeny. Features of both the naso-frontal joint and the morphofunctional system related to macrostomy refute the traditional notion that scolecophidians are fundamentally different from alethinophidians (all other extant snakes). Instead, these features support the controversial hypothesis of scolecophidians as "regressed alethinophidians," in contrast to their traditional placement as the earliest-diverging snake lineage. We propose that Atractaspis and scolecophidians fall along a morphological continuum, characterized by differing degrees of paedomorphosis. Altogether, a combination of heterochrony and miniaturization provides a mechanism for the derivation of the scolecophidian skull from an ancestral fossorial alethinophidian morphotype, exemplified by the nonminiaturized and less extreme paedomorph Atractaspis.

Development of the squamate naso-palatal complex: detailed 3D analysis of the vomeronasal organ and nasal cavity in the brown anole Anolis sagrei (Squamata: Iguania)

Background

Despite the diverse morphology of the adult squamate naso-palatal complex - consisting of the nasal cavity, vomeronasal organ (VNO), choanal groove, lacrimal duct and superficial palate - little is known about the embryology of these structures. Moreover, there are no comprehensive studies concerning development of the nasal cavity and VNO in relation to the superficial palate. In this investigation, we used X-ray microtomography and histological sections to describe embryonic development of the naso-palatal complex of iguanian lizard, the brown anole (Anolis sagrei). The purpose of the study was to describe the mechanism of formation of adult morphology in this species, which combines the peculiar anole features with typical iguanian conditions. Considering the uncertain phylogenetic position of the Iguania within Squamata, embryological data and future comparative studies may shed new light on the evolution of this large squamate clade.

Results

Development of the naso-palatal complex was divided into three phases: early, middle and late. In the early developmental phase, the vomeronasal pit originates from medial outpocketing of the nasal pit, when the facial prominences are weakly developed. In the middle developmental phase, the following events can be noted: the formation of the frontonasal mass, separation of the vestibulum, appearance of the lacrimal duct, and formation of the choanal groove, which leads to separation of the VNO from the nasal cavity. In late development, the nasal cavity and the VNO attain their adult morphology. The lacrimal duct establishes an extensive connection with the choanal groove, which eventually becomes largely separated from the oral cavity.

Conclusions

Unlike in other tetrapods, the primordium of the lacrimal duct in the brown anole develops largely beyond the nasolacrimal groove. In contrast to previous studies on squamates, the maxillary prominence is found to participate in the initial fusion with the frontonasal mass. Moreover, formation of the choanal groove occurs due to the fusion of the vomerine cushion to the subconchal fold, rather than to the choanal fold. The loss or significant reduction of the lateral nasal concha is secondary. Some features of anole adult morphology, such as the closure of the choanal groove, may constitute adaptations to vomeronasal chemoreception.

A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

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

Combined-evidence analyses of ultraconserved elements and morphological data: an empirical example in iguanian lizards

Genomic datasets generated by next-generation sequencing are increasingly prevalent in phylogenetics, but morphological data are required to phylogenetically place fossils, corroborate molecular hypotheses and date phylogenies. Combined-evidence analyses provide an integrative assessment of tree topology. However, no attempt has been made to simultaneously analyse next-generation genomic datasets and morphological data, and the future of morphology in the context of genomic data is uncertain. I conducted combined-evidence analyses that include genomic and morphological datasets, specifically, with ultraconserved elements and two morphological matrices. In unweighted maximum-likelihood and Bayesian combined-evidence analyses, morphological signal was dwarfed by the ultraconserved elements, and some node support values were reduced relative to ultraconserved element-only analyses. Increasing the weight of morphological characters allowed those data to influence the tree, but weighting subjectivity should be considered in future analyses. More attempts should be made to simultaneously analyse genomic and morphological datasets.

Death is on Our Side: Paleontological Data Drastically Modify Phylogenetic Hypotheses

Fossils are the only remaining evidence of the majority of species that have ever existed, providing a direct window into events in evolutionary history that shaped the diversification of life on Earth. Phylogenies underpin our ability to make sense of evolution but are routinely inferred using only data available from living organisms. Although extinct taxa have been shown to add crucial information for inferring macroevolutionary patterns and processes (such as ancestral states, paleobiogeography and diversification dynamics), the role fossils play in reconstructing phylogeny is controversial. Since the early years of phylogenetic systematics, different studies have dismissed the impact of fossils due to their incompleteness, championed their ability to overturn phylogenetic hypotheses or concluded that their behavior is indistinguishable from that of extant taxa. Based on taxon addition experiments on empirical data matrices, we show that the inclusion of paleontological data has a remarkable effect in phylogenetic inference. Incorporating fossils often (yet not always) induces stronger topological changes than increasing sampling of extant taxa. Fossils also produce unique topological rearrangements, allowing the exploration of regions of treespace that are never visited by analyses of only extant taxa. Previous studies have proposed a suite of explanations for the topological behavior of fossils, such as their retention of unique morphologies or their ability to break long branches. We develop predictive models that demonstrate that the possession of distinctive character state combinations is the primary predictor of the degree of induced topological change, and that the relative impact of taxa (fossil and extant) can be predicted to some extent before any phylogenetic analysis. Our results bolster the consensus of recent empirical studies by showing the unique role of paleontological data in phylogenetic inference, and provide the first quantitative assessment of its determinants, with broad consequences for the design of taxon sampling in both morphological and total-evidence analyses. [phylogeny, morphology, fossils, parsimony, Bayesian inference.]

Exploring the impact of morphology, multiple sequence alignment and choice of optimality criteria in phylogenetic inference: a case study with the Neotropical orb-weaving spider genus Wagneriana (Araneae: Araneidae)

We present a total evidence phylogenetic analysis of the Neotropical orb-weaving spider genus Wagneriana and discuss the phylogenetic impacts of methodological choices. We analysed 167 phenotypic characters and nine loci scored for 115 Wagneriana and outgroups, including 46 newly sequenced species. We compared total evidence analyses and molecular-only analyses to evaluate the impact of phenotypic evidence, and we performed analyses using the programs POY, TNT, RAxML, GARLI, IQ-TREE and MrBayes to evaluate the effects of multiple sequence alignment and optimality criteria. In all analyses, Wagneriana carimagua and Wagneriana uropygialis were nested in the genera Parawixia and Alpaida, respectively, and the remaining species of Wagneriana fell into three main clades, none of which formed a pair of sister taxa. However, sister-group relationships among the main clades and their internal relationships were strongly influenced by methodological choices. Alignment methods had comparable topological effects to those of optimality criteria in terms of ‘subtree pruning and regrafting’ moves. The inclusion of phenotypic evidence, 2.80–3.05% of the total evidence matrices, increased support irrespective of the optimality criterion used. The monophyly of some groups was recovered only after the addition of morphological characters. A new araneid genus, Popperaneus gen. nov., is erected, and Paraverrucosa is resurrected. Four new synonymies and seven new combinations are proposed.