Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

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

MOTIVATION

Phylogenetics has moved into the era of genomics, incorporating enormous volumes of data to study questions at both shallow and deep scales. With this increase in information, phylogeneticists need new tools and skills to manipulate and analyze these data. To facilitate these tasks and encourage reproducibility, the community is increasingly moving toward automated workflows.

RESULTS

Here we present pipesnake, a phylogenomics pipeline written in Nextflow for the processing, assembly, and phylogenetic estimation of genomic data from short-read sequences. pipesnake is an easy to use and efficient software package designed for this next era in phylogenetics.

AVAILABILITY AND IMPLEMENTATION

pipesnake is publicly available on GitHub at https://github.com/AusARG/pipesnake and accompanied by documentation and a wiki/tutorial.

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

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.

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.

The Adrenal Gland of Squamata (Reptilia): A Comparative Overview

The adrenal gland is a complex endocrine organ composed of two components: a steroidogenic tissue, which produces steroid hormones, and a chromaffin tissue, which mainly produces norepinephrine and epinephrine. Through evolution, their relationships with each other changed. They begin as isolated chromaffin and steroidogenic cell aggregates, typical of fish, and end with the advanced compact gland, typical of mammals, which consists of an external steroidogenic cortical zone and an internal chromaffin medullary zone. The adrenal gland of reptiles is unique because, with few exceptions, it is near the gonads and genital ducts, and the chromaffin and steroidogenic tissues are closely associated. However, the degree of mixing is variable. For example, in Squamata, the mixing degree of chromaffin and steroidogenic tissues, their reciprocal position in the gland, and the relative quantities of norepinephrine and epinephrine secreted by the chromaffin cells are extremely variable. This variability could be related to the phylogenetic history of the species. After a brief discussion of the adrenal gland and its main functions in vertebrates, this overview will examine the general characteristics of the adrenal gland of squamates, the differences in morphology of the gland, and the possible relationships with the phylogeny of the different species.

Analyzing the biological traits of snakes in Avicenna's Canon of medicine and making a comparison with contemporary serpentology

Snakes are a group of reptiles in the order of squamata that have been the subject of scientific research for years. This study aimed to determine the biological characteristics of snakes that were mentioned in Avicenna's Canon of Medicine to compare with modern serpentology. Data were obtained from the Canon of Medicine using keywords concerning snakes and also appropriate literature in PubMed, Scopus, Web of Science, Scientific Information Database (SID), and IranDoc. Our results revealed that Avicenna divided snakes into three classes including highly, moderately, and slightly venomous snakes that are comparable to modern serpontology. Moreover, Avicenna clarified physiological factors such as age, gender, size, psychological state, hunger state, physical features, the type of living climate, habitat, and the time of snakebite. Considering the snake features presented in the Canon of Medicine, even though is not feasible to make a full comparison between Avicenna's snake knowledge and modern serpentology but some characteristics can be still applicable. Conclusively, the present study suggested that the criteria used for the identification and classification of snakes have changed from the medieval ages to the nowadays.

When lizards try out a more plant-based lifestyle: The macroevolution of mutualistic lizard-plant-interactions (Squamata: Sauria/Lacertilia)

Pollination and seed dispersal of plants by animals are key mutualistic processes for the conservation of plant diversity and ecosystem functioning. Although different animals frequently act as pollinators or seed dispersers, some species can provide both functions, so-called 'double mutualists', suggesting that the evolution of pollination and seed dispersal may be linked. Here, we assess the macroevolution of mutualistic behaviours in lizards (Lacertilia) by applying comparative methods to a phylogeny comprising 2,838 species. We found that both flower visitation (potential pollination) (recorded in 64 species [2.3% of total] across 9 families) and seed dispersal (recorded in 382 species [13,5% of total] across 26 families) have evolved repeatedly in Lacertilia. Furthermore, we found that seed dispersal activity pre-dated flower visitation and that the evolution of seed dispersal activity and flower visitation was correlated, illustrating a potential evolutionary mechanism behind the emergence of double mutualisms. Finally, we provide evidence that lineages with flower visitation or seed dispersal activity have higher diversification rates than lineages lacking these behaviours. Our study illustrates the repeated innovation of (double) mutualisms across Lacertilia and we argue that island settings may provide the ecological conditions under which (double) mutualisms persist over macroevolutionary timescales.

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Assessing the relative performance of fast molecular dating methods for phylogenomic data

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

Parthenogenesis doubles the rate of amino acid substitution in whiptail mitochondria

Sexual reproduction is ubiquitous in the natural world, suggesting that sex must have extensive benefits to overcome the cost of males compared to asexual reproduction. One hypothesized advantage of sex with strong theoretical support is that sex plays a role in removing deleterious mutations from the genome. Theory predicts that transitions to asexuality should lead to the suppression of recombination and segregation and, in turn, weakened natural selection, allowing for the accumulation of slightly deleterious mutations. We tested this prediction by estimating the dN/dS ratios in asexual vertebrate lineages in the genus Aspidoscelis using whole mitochondrial genomes from seven asexual and five sexual species. We found higher dN/dS ratios in asexual Aspidoscelis species, indicating that asexual whiptails accumulate nonsynonymous substitutions due to weaker purifying selection. Additionally, we estimated nucleotide diversity and found that asexuals harbor significantly less diversity. Thus, despite their recent origins, slightly deleterious mutations accumulated rapidly enough in asexual lineages to be detected. We provide empirical evidence to corroborate the connection between asexuality and increased amino acid substitutions in asexual vertebrate lineages.

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.

Phylogeography of the smooth greensnake, Opheodrys vernalis (Squamata: Colubridae): divergent lineages and variable demographics in a widely distributed yet enigmatic species

Phylogeographic studies can uncover robust details about the population structure, demographics, and diversity of species. The smooth greensnake, Opheodrys vernalis, is a small, cryptic snake occupying mesic grassland and sparsely wooded habitats. Although O. vernalis has a wide geographical range, many metapopulations are patchy and some are declining. We used mitochondrial DNA and double digest restriction-site associated DNA sequencing to construct the first phylogeographic assessment of O. vernalis. Genomic analysis of 119 individuals (mitochondrial DNA) and a subset of another 45 smooth greensnakes (nuclear DNA; N = 3031 single nucleotide polymorphisms) strongly supports two longitudinally separated lineages, with admixture in the Great Lakes region. Post-Pleistocene secondary contact best explains admixture from populations advancing northwards. Overall, populations expressed low heterozygosity, variable inbreeding rates, and moderate to high differentiation. Disjunct populations in the Rocky Mountains and central Great Plains regions might be contracting relicts, whereas northerly populations in more continuous mesic habitats (e.g., Prairie Pothole region, southern Canada) had signals of population expansion. Broadly, conservation management efforts should be focused on local populations, because habitat connectivity may facilitate gene flow and genetic diversity.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

India's biogeographic history through the eyes of blindsnakes- filling the gaps in the global typhlopoid phylogeny

The Indian subcontinent's unique geological history is reflected in the diverse assemblage of its biota. The blindsnake superfamily Typhlopoidea, with its unique mix of ancient as well as younger lineages in Asia, provides an opportunity to understand the various biotic exchange scenarios proposed for the Indian landmass. In this study, we aim to understand the biogeographic origins of the four genera of typhlopoids found in India and to decipher their times and modes of arrival in the subcontinent. Five nuclear markers were sequenced for 12 samples collected from across India, encompassing all four genera under study. Published sequences of typhlopoid genera were compiled and combined with Indian sequences to generate a global dataset. Phylogenetic relationships were reconstructed using maximum likelihood and Bayesian inference methods. Divergence times were estimated using BEAST 1.8.2. Ancestral geographical ranges were estimated using DEC + J, implemented in BioGeoBEARS. Divergence time estimates suggest that Gerrhopilus is an ancient lineage, and the lineage leading to it was present on the Indian landmass since the last 100 million years. The other three genera are more recent dispersals into India, possibly trans-oceanic. Biogeographic reconstructions suggest an East Gondwanan origin for Typhlopoidea, an African origin for Grypotyphlops and an Asian origin for Indotyphlops and Argyrophis. It appears that India harbours a combination of ancient and more recently dispersed lineages of typhlopoids. The genus Gerrhopilus is of Gondwanan origin that likely dispersed out of India into Southeast Asia. The other genera are intrusive elements that dispersed into India from Africa (Grypotyphlops) and Asia (Indotyphlops and possibly Argyrophis) post break-up of Gondwana. Thus, our study provides further evidence on the ability of blindsnakes to undergo long distance trans-oceanic dispersal. Results also suggest an Asian origin for typhlopoids from Australasia, Philippines and Wallacea.

DNA barcoding reveals distinct population of Dopasia gracilis (Squamata: Anguidae) in Mizoram, Northeast India

The DNA barcode data of Asian Glass Lizard, Dopasia gracilis, is limited in the global database, especially from India. The present study aimed to generate a barcode sequence of morphologically identified D. gracilis from the Mizoram state in northeast India and compared with other Anguidae species. The studied species showed monophyletic clustering in the Bayesian analysis (BA) phylogeny with strong posterior probability support and also discriminated sufficient Kimura 2 parameter genetic distances. The barcode data of D. gracilis revealed high intra-species genetic variability and formed two clusters in BA phylogeny. The Templeton, Crandall, and Sing network also depicted four different haplotypes within the barcode sequences of D. gracilis. The DNA sequences generated from northeast India showed 6.5-6.6% and 7.3% genetic distances with the sequences generated from Yunnan Province and Tibetan Plateau, respectively. Considering the high genetic distances, multiple clustering, and distinct haplotypes, the present study assumed the presence of possible cryptic diversity of D. gracilis in the Indochina sub-region and a distinct population in northeast India. We recommended the generation of more DNA information from different localities to elucidate the actual diversity of D. gracilis within the known range distribution.

Lizards from the Lost World: two new species and evolutionary relationships of the Pantepui highland Riolama (Gymnophthalmidae)

The Pantepui region of northern South America harbours an endemic fauna that differs dramatically from those of the surrounding lowland rainforests and savannas. A component of this unique fauna is Riolama, a poorly known genus of microteiid lizards with four described and two undescribed species restricted to tepui mountains. We here implement an integrative approach to formally describe the two unnamed species and investigate the phylogenetic relationships and timing of diversification in Riolama using a fossil-calibrated molecular approach. Our results suggest that diversification initiated in Riolama during the Oligocene (c. 28 Mya), thereby characterizing the genus as an ancient lineage. This supports the Plateau biogeographic hypothesis to explain the diversification of the Pantepui fauna. Our divergence time estimation analysis also provides an updated temporal framework for the diversification of the highly diverse Gymnophthalmidae clade.

Optimizing the widely used nuclear protein-coding gene primers in beetle phylogenies and their application in the genus Sasajiscymnus Vandenberg (Coleoptera: Coccinellidae)

Advances in genomic biology and the increasing availability of genomic resources allow developing hundreds of nuclear protein-coding (NPC) markers, which can be used in phylogenetic research. However, for low taxonomic levels, it may be more practical to select a handful of suitable molecular loci for phylogenetic inference. Unfortunately, the presence of degenerate primers of NPC markers can be a major impediment, as the amplification success rate is low and they tend to amplify nontargeted regions. In this study, we optimized five NPC fragments widely used in beetle phylogenetics (i.e., two parts of carbamoyl-phosphate synthetase: CADXM and CADMC, Topoisomerase, Wingless and Pepck) by reducing the degenerate site of primers and the length of target genes slightly. These five NPC fragments and 6 other molecular loci were amplified to test the monophyly of the coccinellid genus Sasajiscymnus Vandenberg. The analysis of our molecular data set clearly supported the genus Sasajiscymnus may be monophyletic but confirmation with an extended sampling is required. A fossil-calibrated chronogram was generated by BEAST, indicating an origin of the genus at the end of the Cretaceous (77.87 Myr). Furthermore, a phylogenetic informativeness profile was generated to compare the phylogenetic properties of each gene more explicitly. The results showed that COI provides the strongest phylogenetic signal among all the genes, but Pepck, Topoisomerase, CADXM and CADMC are also relatively informative. Our results provide insight into the evolution of the genus Sasajiscymnus, and also enrich the molecular data resources for further study.

The Mitochondrial Genome of the Phytopathogenic Fungus Bipolaris sorokiniana and the Utility of Mitochondrial Genome to Infer Phylogeny of Dothideomycetes

A number of species in Bipolaris are important plant pathogens. Due to a limited number of synapomorphic characters, it is difficult to perform species identification and to estimate phylogeny of Bipolaris based solely on morphology. In this study, we sequenced the complete mitochondrial genome of Bipolaris sorokiniana, and presented the detailed annotation of the genome. The B. sorokiniana mitochondrial genome is 137,775 bp long, and contains two ribosomal RNA genes, 12 core protein-coding genes, 38 tRNA genes. In addition, two ribosomal protein genes (rps3 gene and rps5 gene) and the fungal mitochondrial RNase P gene (rnpB) are identified. The large genome size is mostly determined by the presence of numerous intronic and intergenic regions. A total of 28 introns are inserted in eight core protein-coding genes. Together with the published mitochondrial genome sequences, we conducted a preliminary phylogenetic inference of Dothideomycetes under various datasets and substitution models. The monophyly of Capnodiales, Botryosphaeriales and Pleosporales are consistently supported in all analyses. The Venturiaceae forms an independent lineage, with a distant phylogenetic relationship to Pleosporales. At the family level, the Mycosphaerellaceae, Botryosphaeriaceae. Phaeosphaeriaceae, and Pleosporaceae are recognized in the majority of trees.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

Interrogating Genomic-Scale Data for Squamata (Lizards, Snakes, and Amphisbaenians) Shows no Support for Key Traditional Morphological Relationships

Genomics is narrowing uncertainty in the phylogenetic structure for many amniote groups. For one of the most diverse and species-rich groups, the squamate reptiles (lizards, snakes, and amphisbaenians), an inverse correlation between the number of taxa and loci sampled still persists across all publications using DNA sequence data and reaching a consensus on the relationships among them has been highly problematic. In this study, we use high-throughput sequence data from 289 samples covering 75 families of squamates to address phylogenetic affinities, estimate divergence times, and characterize residual topological uncertainty in the presence of genome-scale data. Importantly, we address genomic support for the traditional taxonomic groupings Scleroglossa and Macrostomata using novel machine-learning techniques. We interrogate genes using various metrics inherent to these loci, including parsimony-informative sites (PIS), phylogenetic informativeness, length, gaps, number of substitutions, and site concordance to understand why certain loci fail to find previously well-supported molecular clades and how they fail to support species-tree estimates. We show that both incomplete lineage sorting and poor gene-tree estimation (due to a few undesirable gene properties, such as an insufficient number of PIS), may account for most gene and species-tree discordance. We find overwhelming signal for Toxicofera, and also show that none of the loci included in this study supports Scleroglossa or Macrostomata. We comment on the origins and diversification of Squamata throughout the Mesozoic and underscore remaining uncertainties that persist in both deeper parts of the tree (e.g., relationships between Dibamia, Gekkota, and remaining squamates; among the three toxicoferan clades Iguania, Serpentes, and Anguiformes) and within specific clades (e.g., affinities among gekkotan, pleurodont iguanians, and colubroid families).

Phylogenetic relationships and divergence dating in the Glass Lizards (Anguinae)

The Glass Lizards are a subfamily (Anguinae) of Anguid Lizards with an elongated limbless body plan that occur throughout the Northern Hemisphere primarily in North America, Europe, and Asia, but also have a presence in North Africa and Indonesia. We used twenty-five nuclear loci (15,191 bp) and 2090 bp of the mtDNA genome to generate a phylogeny containing all known species groups to explore species relationships within the group as well as divergence dating. We also examined the group in the context of a coalescent species tree analysis and species delimitation. All major lineages were found to be monophyletic with potential cryptic diversity in some. The Anguinae first appeared in the Eocene and most lineages were present by the beginning of the Miocene. The Anguinae originated in Europe from an Anguidae ancestor that crossed the Thulean land bridge, spreading to Asia after the drying of the Turgai Sea, then across Beringia as the climate permitted. A species tree analyses found support for the major Anguinae lineages and species delimitation supported accepted species.

A Phylogenomic Framework and Divergence History of Cephalochordata Amphioxus

Amphioxus, or cephalochordates, are often used as the living invertebrate proxy of vertebrate ancestors and are widely used as evolutionary biology models of chordates. However, their phylogeny, divergence history, and speciation characteristics remain poorly understood, and phylogenomic studies to explore these problems lacking entirely from the literature. Here, we determined a new transcriptome of Branchiostoma japonicum. Combined with mass sequences of all other 18 species, a 19-way phylogeny was constructed via multiple methods (ML, BI, PhyloBayes, and ASTRAL), consistently supporting a phylogeny of [(B. belcheri + B. japonicum) + (B. lanceolatum + B. floridae) + Asymmetron lucayanum] in amphioxus. Congruent phylogenetic signals were found across mitochondrial genes, 12S RNA, and complete mitochondrial genomes according to previous reports, indicating that 12S RNA may have potential as a molecular marker for phylogenetic analysis in amphioxus. Molecular dating analysis indicated a radiation of the cephalochordates during the Cretaceous (∼104-61 million years ago), supporting an association between the diversification and speciation of cephalochordates with continental drift and associated changes in their respective habitats during this time. The identified functional enrichment analysis for species-specific domains indicated that their function mainly involves immune response, apoptosis, and lipid metabolism and utilization, signaling that pathogens and changes of energy requirements are an important driving force for amphioxus speciation. This study represents the first large-scale phylogenomic analysis of most major amphioxus genera based on phylogenomic data, providing a new perspective on both phylogeny and divergence speciation of cephalochordates.

Comparative venom-gland transcriptomics and venom proteomics of four Sidewinder Rattlesnake (Crotalus cerastes) lineages reveal little differential expression despite individual variation

Changes in gene expression can rapidly influence adaptive traits in the early stages of lineage diversification. Venom is an adaptive trait comprised of numerous toxins used for prey capture and defense. Snake venoms can vary widely between conspecific populations, but the influence of lineage diversification on such compositional differences are unknown. To explore venom differentiation in the early stages of lineage diversification, we used RNA-seq and mass spectrometry to characterize Sidewinder Rattlesnake (Crotalus cerastes) venom. We generated the first venom-gland transcriptomes and complementary venom proteomes for eight individuals collected across the United States and tested for expression differences across life history traits and between subspecific, mitochondrial, and phylotranscriptomic hypotheses. Sidewinder venom was comprised primarily of hemorrhagic toxins, with few cases of differential expression attributable to life history or lineage hypotheses. However, phylotranscriptomic lineage comparisons more than doubled instances of significant expression differences compared to all other factors. Nevertheless, only 6.4% of toxins were differentially expressed overall, suggesting that shallow divergence has not led to major changes in Sidewinder venom composition. Our results demonstrate the need for consensus venom-gland transcriptomes based on multiple individuals and highlight the potential for discrepancies in differential expression between different phylogenetic hypotheses.

Systematics of Clupeiformes and testing for ecological limits on species richness in a trans-marine/freshwater clade

ABSTRACT Clupeiformes (herring, sardines, shad, anchovies and allies) are a globally distributed clade with nearly 400 marine, freshwater, and diadromous species. Although best known as filter feeding fishes that form large schools, this group occupies a diverse array of trophic guilds and habitats. Theory suggests that species richness in clades is modulated by ecological limits, which results in diversity-dependent clade growth, a pattern that most clades exhibit. As a trans-marine/freshwater clade that has undergone repeated transitions between marine and freshwaters, Clupeiformes are an excellent system for investigating the interplay between ecological diversity and macroevolutionary dynamics. In this study we review the systematics of Clupeiformes and explore discordance in phylogenetic relationships and divergence times between mitochondrial and nuclear loci. We then use comparative methods to test whether ecological limits regulate diversity in Clupeiformes. We find discordance in phylogenetic relationships at various taxonomic scales, but also considerable agreement between genomes. Our results suggest that trans-marine/freshwater clades are able to circumvent ecological limits on clade growth at regional, but not on local scales. Our study demonstrates that phylogenies are a critical link between ecology and macroevolutionary dynamics, and suggests habitat transitions can play a key role in shaping diversity patterns, particularly in the neotropics.

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

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Evolutionary Dynamics of the Gametologous CTNNB1 Gene on the Z and W Chromosomes of Snakes

Snakes exhibit genotypic sex determination with female heterogamety (ZZ males and ZW females), and the state of sex chromosome differentiation also varies among lineages. To investigate the evolutionary history of homologous genes located in the nonrecombining region of differentiated sex chromosomes in snakes, partial sequences of the gametologous CTNNB1 gene were analyzed for 12 species belonging to henophid (Cylindrophiidae, Xenopeltidae, and Pythonidae) and caenophid snakes (Viperidae, Elapidae, and Colubridae). Nonsynonymous/synonymous substitution ratios (Ka/Ks) in coding sequences were low (Ka/Ks < 1) between CTNNB1Z and CTNNB1W, suggesting that these 2 genes may have similar functional properties. However, frequencies of intron sequence substitutions and insertion–deletions were higher in CTNNB1Z than CTNNB1W, suggesting that Z-linked sequences evolved faster than W-linked sequences. Molecular phylogeny based on both intron and exon sequences showed the presence of 2 major clades: 1) Z-linked sequences of Caenophidia and 2) W-linked sequences of Caenophidia clustered with Z-linked sequences of Henophidia, which suggests that the sequence divergence between CTNNB1Z and CTNNB1W in Caenophidia may have occurred by the cessation of recombination after the split from Henophidia.

Lepidosaurian diversity in the Mesozoic-Palaeogene: the potential roles of sampling biases and environmental drivers

Lepidosauria is a speciose clade with a long evolutionary history, but there have been few attempts to explore its taxon richness through time. Here we estimate patterns of terrestrial lepidosaur genus diversity for the Triassic-Palaeogene (252-23 Ma), and compare observed and sampling-corrected richness curves generated using Shareholder Quorum Subsampling and classical rarefaction. Generalized least-squares regression (GLS) is used to investigate the relationships between richness, sampling and environmental proxies. We found low levels of richness from the Triassic until the Late Cretaceous (except in the Kimmeridgian-Tithonian of Europe). High richness is recovered for the Late Cretaceous of North America, which declined across the K-Pg boundary but remained relatively high throughout the Palaeogene. Richness decreased following the Eocene-Oligocene Grande Coupure in North America and Europe, but remained high in North America and very high in Europe compared to the Late Cretaceous; elsewhere data are lacking. GLS analyses indicate that sampling biases (particularly, the number of fossil collections per interval) are the best explanation for long-term face-value genus richness trends. The lepidosaur fossil record presents many problems when attempting to reconstruct past diversity, with geographical sampling biases being of particular concern, especially in the Southern Hemisphere.

Novel Approaches for Phylogenetic Inference from Morphological Data and Total-Evidence Dating in Squamate Reptiles (Lizards, Snakes, and Amphisbaenians)

Here, I combine previously underutilized models and priors to perform more biologically realistic phylogenetic inference from morphological data, with an example from squamate reptiles. When coding morphological characters, it is often possible to denote ordered states with explicit reference to observed or hypothetical ancestral conditions. Using this logic, we can integrate across character-state labels and estimate meaningful rates of forward and backward transitions from plesiomorphy to apomorphy. I refer to this approach as MkA, for “asymmetric.” The MkA model incorporates the biological reality of limited reversal for many phylogenetically informative characters, and significantly increases likelihoods in the empirical data sets. Despite this, the phylogeny of Squamata remains contentious. Total-evidence analyses using combined morphological and molecular data and the MkA approach tend toward recent consensus estimates supporting a nested Iguania. However, support for this topology is not unambiguous across data sets or analyses, and no mechanism has been proposed to explain the widespread incongruence between partitions, or the hidden support for various topologies in those partitions. Furthermore, different morphological data sets produced by different authors contain both different characters and different states for the same or similar characters, resulting in drastically different placements for many important fossil lineages. Effort is needed to standardize ontology for morphology, resolve incongruence, and estimate a robust phylogeny. The MkA approach provides a preliminary avenue for investigating morphological evolution while accounting for temporal evidence and asymmetry in character-state changes.

Do Macrophylogenies Yield Stable Macroevolutionary Inferences? An Example from Squamate Reptiles

Advances in the generation, retrieval, and analysis of phylogenetic data have enabled researchers to create phylogenies that contain many thousands of taxa. These "macrophylogenies"-large trees that typically derive from megaphylogeny, supermatrix, or supertree approaches-provide researchers with an unprecedented ability to conduct evolutionary analyses across broad phylogenetic scales. Many studies have now used these phylogenies to explore the dynamics of speciation, extinction, and phenotypic evolution across large swaths of the tree of life. These trees are characterized by substantial phylogenetic uncertainty on multiple levels, and the stability of macroevolutionary inferences from these data sets has not been rigorously explored. As a case study, we tested whether five recently published phylogenies for squamate reptiles-each consisting of more than 4000 species-yield congruent inferences about the processes that underlie variation in species richness across replicate evolutionary radiations of Australian snakes and lizards. We find discordance across the five focal phylogenies with respect to clade age and several diversification rate metrics, and in the effects of clade age on species richness. We also find that crown clade ages reported in the literature on these Australian groups are in conflict with all of the large phylogenies examined. Macrophylogenies offer an unprecedented opportunity to address evolutionary and ecological questions at broad phylogenetic scales, but accurately representing the uncertainty that is inherent to such analyses remains a critical challenge to our field. [Australia; macroevolution; macrophylogeny; squamates; time calibration.].

Genomic timetree and historical biogeography of Caribbean island ameiva lizards (Pholidoscelis: Teiidae)

The phylogenetic relationships and biogeographic history of Caribbean island ameivas (Pholidoscelis) are not well-known because of incomplete sampling, conflicting datasets, and poor support for many clades. Here, we use phylogenomic and mitochondrial DNA datasets to reconstruct a well-supported phylogeny and assess historical colonization patterns in the group. We obtained sequence data from 316 nuclear loci and one mitochondrial marker for 16 of 19 extant species of the Caribbean endemic genus Pholidoscelis. Phylogenetic analyses were carried out using both concatenation and species tree approaches. To estimate divergence times, we used fossil teiids to calibrate a timetree which was used to elucidate the historical biogeography of these lizards. All phylogenetic analyses recovered four well-supported species groups (clades) recognized previously and supported novel relationships of those groups, including a (P. auberi + P. lineolatus) clade (western + central Caribbean), and a (P. exsul + P. plei) clade (eastern Caribbean). Divergence between Pholidoscelis and its sister clade was estimated to have occurred ~25 Ma, with subsequent diversification on Caribbean islands occurring over the last 11 Myr. Of the six models compared in the biogeographic analyses, the scenario which considered the distance among islands and allowed dispersal in all directions best fit the data. These reconstructions suggest that the ancestor of this group colonized either Hispaniola or Puerto Rico from Middle America. We provide a well-supported phylogeny of Pholidoscelis with novel relationships not reported in previous studies that were based on significantly smaller datasets. We propose that Pholidoscelis colonized the eastern Greater Antilles from Middle America based on our biogeographic analysis, phylogeny, and divergence time estimates. The closing of the Central American Seaway and subsequent formation of the modern Atlantic meridional overturning circulation may have promoted dispersal in this group.

Phylogenomic analyses of 539 highly informative loci dates a fully resolved time tree for the major clades of living turtles (Testudines)

Accurate time-calibrated phylogenies are the centerpiece of many macroevolutionary studies, and the relationship between the size and scale of molecular data sets and the density and accuracy of fossil calibrations is a key element of time tree studies. Here, we develop a target capture array specifically for living turtles, compare its efficiency to an ultraconserved element (UCE) dataset, and present a time-calibrated molecular phylogeny based on 539 nuclear loci sequenced from 26 species representing the breadth of living turtle diversity plus outgroups. Our gene array, based on three fully sequenced turtle genomes, is 2.4 times more variable across turtles than a recently published UCE data set for an identical subset of 13 species, confirming that taxon-specific arrays return more informative data per sequencing effort than UCEs. We used our genomic data to estimate the ages of living turtle clades including a mid-late Triassic origin for crown turtles and a mid-Carboniferous split of turtles from their sister group, Archosauria. By specifically excluding several of the earliest potential crown turtle fossils and limiting the age of fossil calibration points to the unambiguous crown lineage Caribemys oxfordiensis from the Late Jurassic (Oxfordian, 163.5-157.3Ma) we corroborate a relatively ancient age for living turtles. We also provide novel age estimates for five of the ten testudine families containing more than a single species, as well as several intrafamilial clades. Most of the diversity of crown turtles appears to date to the Paleogene, well after the Cretaceous-Paleogene mass extinction 66mya.

Arrival and diversification of mabuyine skinks (Squamata: Scincidae) in the Neotropics based on a fossil-calibrated timetree

BACKGROUND

The evolution of South American Mabuyinae skinks holds significant biogeographic interest because its sister lineage is distributed across the African continent and adjacent islands. Moreover, at least one insular species, Trachylepis atlantica, has independently reached the New World through transoceanic dispersal. To clarify the evolutionary history of both Neotropical lineages, this study aimed to infer an updated timescale using the largest species and gene sampling dataset ever assembled for this group. By extending the analysis to the Scincidae family, we could employ fossil information to estimate mabuyinae divergence times and carried out a formal statistical biogeography analysis. To unveil macroevolutionary patterns, we also inferred diversification rates for this lineage and evaluated whether the colonization of South American continent significantly altered the mode of Mabuyinae evolution.

METHODS

A time-calibrated phylogeny was inferred under the Bayesian framework employing fossil information. This timetree was used to (i) evaluate the historical biogeography of mabuiyines using the statistical approach implemented in BioGeoBEARS; (ii) estimate macroevolutionary diversification rates of the South American Mabuyinae lineages and the patterns of evolution of selected traits, namely, the mode of reproduction, body mass and snout-vent length; (iii) test the hypothesis of differential macroevolutionary patterns in South American lineages in BAMM and GeoSSE; and (iv) re-evaluate the ancestral state of the mode of reproduction of mabuyines.

RESULTS

Our results corroborated the hypothesis that the occupation of the South American continent by Mabuyinae consisted of two independent dispersion events that occurred between the Oligocene and the Miocene. We found significant differences in speciation rates between the New World and the remaining Mabuyinae clades only in GeoSSE. The influence of phenotypic traits on diversification rates was not supported by any method. Ancestral state reconstruction suggested that the ancestor of South American mabuyine was likely viviparous.

DISCUSSION

Our analyses further corroborated the existence of a transoceanic connection between Africa and South America in the Eocene/Oligocene period (Atlantogea). Following colonization of the isolated South America and subsequent dispersal through the continent by the ancestral mabuyine stock, we detected no difference in macroevolutionary regimes of New World clades. This finding argued against the ecological opportunity model as an explanation for the diversity of living mabuyines.

Evolutionary stability of sex chromosomes in snakes

Amniote vertebrates possess various mechanisms of sex determination, but their variability is not equally distributed. The large evolutionary stability of sex chromosomes in viviparous mammals and birds was believed to be connected with their endothermy. However, some ectotherm lineages seem to be comparably conserved in sex determination, but previously there was a lack of molecular evidence to confirm this. Here, we document a stability of sex chromosomes in advanced snakes based on the testing of Z-specificity of genes using quantitative PCR (qPCR) across 37 snake species (our qPCR technique is suitable for molecular sexing in potentially all advanced snakes). We discovered that at least part of sex chromosomes is homologous across all families of caenophidian snakes (Acrochordidae, Xenodermatidae, Pareatidae, Viperidae, Homalopsidae, Colubridae, Elapidae and Lamprophiidae). The emergence of differentiated sex chromosomes can be dated back to about 60 Ma and preceded the extensive diversification of advanced snakes, the group with more than 3000 species. The Z-specific genes of caenophidian snakes are (pseudo)autosomal in the members of the snake families Pythonidae, Xenopeltidae, Boidae, Erycidae and Sanziniidae, as well as in outgroups with differentiated sex chromosomes such as monitor lizards, iguanas and chameleons. Along with iguanas, advanced snakes are therefore another example of ectothermic amniotes with a long-term stability of sex chromosomes comparable with endotherms.

Why are some species older than others? A large-scale study of vertebrates

BACKGROUND

Strong variations are observed between and within taxonomic groups in the age of extant species and these differences can clarify factors that render species more vulnerable to extinction. Understanding the factors that influence the resilience of species is thus a key component of evolutionary biology, but it is also of prime importance in a context of climate change and for conservation in general. We explored the effect of extrinsic and intrinsic factors on the timing of the oldest diversification event in over 600 vertebrate species distributed worldwide. We used phylogenetic comparative methods to show that color polymorphism, latitude and reproduction (the latter through its interaction with latitude) affected the timing of the oldest diversification event within a species.

RESULTS

Species from higher latitudes tended to be younger, and colour-polymorphic species were older than monomorphic species. Mode of reproduction was important also, in that the age of oviparous species decreased with latitude, whereas no pattern was apparent for viviparous species. Organisms which have already persisted for a long time may be more likely to deal with future modifications of their environment.

CONCLUSIONS

Species that are colour polymorphic, viviparous, and/or live at low latitudes have exhibited resilience to past environmental changes, and hence may be better able to deal with current climate change.

Detecting the Anomaly Zone in Species Trees and Evidence for a Misleading Signal in Higher-Level Skink Phylogeny (Squamata: Scincidae)

The anomaly zone, defined by the presence of gene tree topologies that are more probable than the true species tree, presents a major challenge to the accurate resolution of many parts of the Tree of Life. This discrepancy can result from consecutive rapid speciation events in the species tree. Similar to the problem of long-branch attraction, including more data via loci concatenation will only reinforce the support for the incorrect species tree. Empirical phylogenetic studies often employ coalescent-based species tree methods to avoid the anomaly zone, but to this point these studies have not had a method for providing any direct evidence that the species tree is actually in the anomaly zone. In this study, we use 16 species of lizards in the family Scincidae to investigate whether nodes that are difficult to resolve place the species tree within the anomaly zone. We analyze new phylogenomic data (429 loci), using both concatenation and coalescent-based species tree estimation, to locate conflicting topological signal. We then use the unifying principle of the anomaly zone, together with estimates of ancestral population sizes and species persistence times, to determine whether the observed phylogenetic conflict is a result of the anomaly zone. We identify at least three regions of the Scincidae phylogeny that provide demographic signatures consistent with the anomaly zone, and this new information helps reconcile the phylogenetic conflict in previously published studies on these lizards. The anomaly zone presents a real problem in phylogenetics, and our new framework for identifying anomalous relationships will help empiricists leverage their resources appropriately for investigating and overcoming this challenge.