When do species-tree and concatenated estimates disagree? An empirical analysis with higher-level scincid lizard phylogeny

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Multi-locus phylogeny of the catfish genus Ictalurus Rafinesque, 1820 (Actinopterygii, Siluriformes) and its systematic and evolutionary implications

BACKGROUND

Ictalurus is one of the most representative groups of North American freshwater fishes. Although this group has a well-studied fossil record and has been the subject of several morphological and molecular phylogenetic studies, incomplete taxonomic sampling and insufficient taxonomic studies have produced a rather complex classification, along with intricate patterns of evolutionary history in the genus that are considered unresolved and remain under debate.

RESULTS

Based on four loci and the most comprehensive taxonomic sampling analyzed to date, including currently recognized species, previously synonymized species, undescribed taxa, and poorly studied populations, this study produced a resolved phylogenetic framework that provided plausible species delimitation and an evolutionary time framework for the genus Ictalurus.

CONCLUSIONS

Our phylogenetic hypothesis revealed that Ictalurus comprises at least 13 evolutionary units, partially corroborating the current classification and identifying populations that emerge as putative undescribed taxa. The divergence times of the species indicate that the diversification of Ictalurus dates to the early Oligocene, confirming its status as one of the oldest genera within the family Ictaluridae.

The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei

In order to adapt to diverse habitats, organisms often evolve corresponding adaptive mechanisms to cope with their survival needs. The species-rich family of Scincidae contains both limbed and limbless species, which differ fundamentally in their locomotor demands, such as relying on the movement of limbs or only body swing to move. Locomotion requires energy, and different types of locomotion have their own energy requirements. Mitochondria are the energy factories of living things, which provide a lot of energy for various physiological activities of organisms. Therefore, mitochondrial genomes could be tools to explore whether the limb loss of skinks are selected by adaptive evolution. Isopachys gyldenstolpei is a typical limbless skink. Here, we report the complete mitochondrial genomes of I. gyldenstolpei, Sphenomorphus indicus, and Tropidophorus hainanus. The latter two species were included as limbed comparator species to the limbless I. gyldenstolpei. The results showed that the full lengths of the mitochondrial genomes of I. gyldenstolpei, S. indicus, and T. hainanus were 17,210, 16,944, and 17,001 bp, respectively. Three mitochondrial genomes have typical circular double-stranded structures similar to other reptiles, including 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and the control region. Three mitochondrial genomes obtained in this study were combined with fifteen mitochondrially complete genomes of Scincidae in the NCBI database; the phylogenetic relationship between limbless I. gyldenstolpei and limbed skinks (S. indicus and T. hainanus) is discussed. Through BI and ML trees, Sphenomorphinae and Mabuyinae were monophyletic, while the paraphyly of Scincinae was also recovered. The limbless skink I. gyldenstolpei is closer to the species of Tropidophorus, which has formed a sister group with (T. hainanus + T. hangman). In the mitochondrial genome adaptations between limbless I. gyldenstolpei and limbed skinks, one positively selected site was found in the branch-site model analysis, which was located in ND2 (at position 28, BEB value = 0.907). Through analyzing the protein structure and function of the selected site, we found it was distributed in mitochondrial protein complex I. Positive selection of some mitochondrial genes in limbless skinks may be related to the requirement of energy to fit in their locomotion. Further research is still needed to confirm this conclusion though.

Evolution of diel activity patterns in skinks (Squamata: Scincidae), the world's second-largest family of terrestrial vertebrates

Many animals have strict diel activity patterns, with unique adaptations for either diurnal or nocturnal activity. Diel activity is phylogenetically conserved, yet evolutionary shifts in diel activity occur and lead to important changes in an organism's morphology, physiology, and behavior. We use phylogenetic comparative methods to examine the evolutionary history of diel activity in skinks, one of the largest families of terrestrial vertebrates. We examine how diel patterns are associated with microhabitat, ambient temperatures, and morphology. We found support for a nondiurnal ancestral skink. Strict diurnality in crown group skinks only evolved during the Paleogene. Nocturnal habits are associated with fossorial activity, limb reduction and loss, and warm temperatures. Our results shed light on the evolution of diel activity patterns in a large radiation of terrestrial ectotherms and reveal how both intrinsic biotic and extrinsic abiotic factors can shape the evolution of animal activity patterns.

Historical biogeography, systematics, and integrative taxonomy of the non-Ethiopian speckled pelage brush-furred rats (Lophuromys flavopunctatus group)

BACKGROUND

The speckled-pelage brush-furred rats (Lophuromys flavopunctatus group) have been difficult to define given conflicting genetic, morphological, and distributional records that combine to obscure meaningful accounts of its taxonomic diversity and evolution. In this study, we inferred the systematics, phylogeography, and evolutionary history of the L. flavopunctatus group using maximum likelihood and Bayesian phylogenetic inference, divergence times, historical biogeographic reconstruction, and morphometric discriminant tests. We compiled comprehensive datasets of three loci (two mitochondrial [mtDNA] and one nuclear) and two morphometric datasets (linear and geometric) from across the known range of the genus Lophuromys.

RESULTS

The mtDNA phylogeny supported the division of the genus Lophuromys into three primary groups with nearly equidistant pairwise differentiation: one group corresponding to the subgenus Kivumys (Kivumys group) and two groups corresponding to the subgenus Lophuromys (L. sikapusi group and L. flavopunctatus group). The L. flavopunctatus group comprised the speckled-pelage brush-furred Lophuromys endemic to Ethiopia (Ethiopian L. flavopunctatus members [ETHFLAVO]) and the non-Ethiopian ones (non-Ethiopian L. flavopunctatus members [NONETHFLAVO]) in deeply nested relationships. There were distinctly geographically structured mtDNA clades among the NONETHFLAVO, which were incongruous with the nuclear tree where several clades were unresolved. The morphometric datasets did not systematically assign samples to meaningful taxonomic units or agree with the mtDNA clades. The divergence dating and ancestral range reconstructions showed the NONETHFLAVO colonized the current ranges over two independent dispersal events out of Ethiopia in the early Pleistocene.

CONCLUSION

The phylogenetic associations and divergence times of the L. flavopunctatus group support the hypothesis that paleoclimatic impacts and ecosystem refugia during the Pleistocene impacted the evolutionary radiation of these rodents. The overlap in craniodental variation between distinct mtDNA clades among the NONETHFLAVO suggests unraveling underlying ecomorphological drivers is key to reconciling taxonomically informative morphological characters. The genus Lophuromys requires a taxonomic reassessment based on extensive genomic evidence to elucidate the patterns and impacts of genetic isolation at clade contact zones.

Partitioned Gene-Tree Analyses and Gene-Based Topology Testing Help Resolve Incongruence in a Phylogenomic Study of Host-Specialist Bees (Apidae: Eucerinae)

Incongruence among phylogenetic results has become a common occurrence in analyses of genome-scale data sets. Incongruence originates from uncertainty in underlying evolutionary processes (e.g., incomplete lineage sorting) and from difficulties in determining the best analytical approaches for each situation. To overcome these difficulties, more studies are needed that identify incongruences and demonstrate practical ways to confidently resolve them. Here, we present results of a phylogenomic study based on the analysis 197 taxa and 2,526 ultraconserved element (UCE) loci. We investigate evolutionary relationships of Eucerinae, a diverse subfamily of apid bees (relatives of honey bees and bumble bees) with >1,200 species. We sampled representatives of all tribes within the group and >80% of genera, including two mysterious South American genera, Chilimalopsis and Teratognatha. Initial analysis of the UCE data revealed two conflicting hypotheses for relationships among tribes. To resolve the incongruence, we tested concatenation and species tree approaches and used a variety of additional strategies including locus filtering, partitioned gene-trees searches, and gene-based topological tests. We show that within-locus partitioning improves gene tree and subsequent species-tree estimation, and that this approach, confidently resolves the incongruence observed in our data set. After exploring our proposed analytical strategy on eucerine bees, we validated its efficacy to resolve hard phylogenetic problems by implementing it on a published UCE data set of Adephaga (Insecta: Coleoptera). Our results provide a robust phylogenetic hypothesis for Eucerinae and demonstrate a practical strategy for resolving incongruence in other phylogenomic data sets.

A new species of Proegernia from the Namba Formation in South Australia and the early evolution and environment of Australian egerniine skinks

The diverse living Australian lizard fauna contrasts greatly with their limited Oligo-Miocene fossil record. New Oligo-Miocene fossil vertebrates from the Namba Formation (south of Lake Frome, South Australia) were uncovered from multiple expeditions from 2007 to 2018. Abundant disarticulated material of small vertebrates was concentrated in shallow lenses along the palaeolake edges, now exposed on the western of Lake Pinpa also known from Billeroo Creek 2 km northeast. The fossiliferous lens within the Namba Formation hosting the abundant aquatic (such as fish, platypus Obdurodon and waterfowl) and diverse terrestrial (such as possums, dasyuromorphs and scincids) vertebrates and is hereafter recognized as the Fish Lens. The stratigraphic provenance of these deposits in relation to prior finds in the area is also established. A new egerniine scincid taxon Proegernia mikebulli sp. nov. described herein, is based on a near-complete reconstructed mandible, maxilla, premaxilla and pterygoid. Postcranial scincid elements were also recovered with this material, but could not yet be confidently associated with P. mikebulli. This new taxon is recovered as the sister species to P. palankarinnensis, in a tip-dated total-evidence phylogenetic analysis, where both are recovered as stem Australian egerniines. These taxa also help pinpoint the timing of the arrival of scincids to Australia, with egerniines the first radiation to reach the continent.

Hydrography rather than lip morphology better explains the evolutionary relationship between Gymnogeophagus labiatus and G. lacustris in Southern Brazil (Cichlidae: Geophagini)

ABSTRACT Gymnogeophagus labiatus and G. lacustris have been long recognized as sister species exhibiting different ecological requirements. Gymnogeophagus labiatus occurs in rock bottom rivers in the hydrographic basins of Patos Lagoon (HBP) and Tramandaí River (HBT), while G. lacustris is exclusive from sand bottom coastal lagoons of the HBT. In this study, we used molecular markers, morphological measurements and data from nuptial male coloration to investigate the evolutionary relationship between these species in each hydrographic basin. We found, for all data sets, a closer relationship between G. labiatus and G. lacustris from the HBT than between G. labiatus populations from HBT and HBP. In particular, lip area had a large intraspecific plasticity, being uninformative to diagnose G. lacustris from G. labiatus. Molecular clock-based estimates suggest a recent divergence between species in the HBT (17,000 years ago), but not between G. labiatus from HBP and HBT (3.6 millions of years ago). Finally, we also found a divergent G. labiatus genetic lineage from the Camaquã River, in the HBP. These results show that the current taxonomy of G. labiatus and G. lacustris does not properly represent evolutionary lineages in these species.

Systematics, biogeography and evolution of the Saharo-Arabian naked-toed geckos genus Tropiocolotes

Plate tectonics constitute one of the main mechanisms of biological diversification on Earth, often being associated with cladogenetic events at different phylogenetic levels, as well as with exchange of faunas and floras across previously isolated biogeographic regions. North Africa and Arabia share a complex geological history that dates back to the break-up of the Arabian plate from the African plate ~30-25 Mya, followed by various geological events, such as the formation of the Red Sea or the connection between the African, Arabian and Eurasian plates. Species with Saharo-Arabian distributions have shown a close association between their evolutionary history and these geological events. In this study, we investigate the systematics, biogeography and evolution of the genus Tropiocolotes, a group of small ground-dwelling geckos, comprised by 12 species distributed from the Atlantic coast of North Africa to southwestern Iran. Species delimitation analyses uncovered the existence of high levels of undescribed diversity, with forms here considered at the species level including Tropiocolotes tripolitanus (Mauritania and southern Morocco), T. nattereri (southern Israel) and T. scorteccii (Yemen and Oman). Phylogenetic and biogeographic analyses recovered two main clades, an exclusively African clade and a Saharo-Arabian clade, that split ~25 Mya following the vicariant event mediated by the separation of the Arabian and African plates. The complex geological activity around the Red Sea is associated with the diversification within the Saharo-Arabian clade, including the colonization of North Africa from a second Tropiocolotes group. Results also provide new insights into the geographic distribution of Tropiocolotes nubicus, previously considered as exclusively associated to the Nile River valley, extending its known distribution further west, up to the Central Mountains of the Sahara. Accordingly, the Nile River seems to act as a major biogeographic barrier, separating Tropiocolotes nubicus and T. steudneri in their western and eastern margins, respectively.

TREEasy: An automated workflow to infer gene trees, species trees, and phylogenetic networks from multilocus data

Multilocus genomic data sets can be used to infer a rich set of information about the evolutionary history of a lineage, including gene trees, species trees, and phylogenetic networks. However, user-friendly tools to run such integrated analyses are lacking, and workflows often require tedious reformatting and handling time to shepherd data through a series of individual programs. Here, we present a tool written in Python-TREEasy-that performs automated sequence alignment (with MAFFT), gene tree inference (with IQ-Tree), species inference from concatenated data (with IQ-Tree and RaxML-NG), species tree inference from gene trees (with ASTRAL, MP-EST, and STELLS2), and phylogenetic network inference (with SNaQ and PhyloNet). The tool only requires FASTA files and nine parameters as inputs. The tool can be run as command line or through a Graphical User Interface (GUI). As examples, we reproduced a recent analysis of staghorn coral evolution, and performed a new analysis on the evolution of the "WGD clade" of yeast. The latter revealed novel patterns that were not identified by previous analyses. TREEasy represents a reliable and simple tool to accelerate research in systematic biology (https://github.com/MaoYafei/TREEasy).

Another one bites the dust: taxonomic sampling of a key genus in phylogenomic datasets reveals more non-monophyletic groups in traditional scorpion classification

Historically, morphological characters have been used to support the monophyly, composition, and phylogenetic relationships of scorpion families. Although recent phylogenomic analyses have recovered most of these traditional higher-level relationships as non-monophyletic, certain key taxa have yet to be sampled using a phylogenomic approach. Salient among these is the monotypic genus Caraboctonus Pocock,1893, the type species of the family Caraboctonidae Kraepelin, 1905. Here, we examined the putative monophyly and phylogenetic placement of this family, sampling the library of C. keyserlingi Pocock, 1893 using high throughput transcriptomic sequencing. Our phylogenomic analyses recovered Caraboctonidae as polyphyletic due to the distant placement of the genera Caraboctonus and Hadrurus Thorell, 1876. Caraboctonus was stably recovered as the sister-group of the monotypic family Superstitioniidae Stahnke, 1940, whereas Hadrurus formed an unstable relationship with Uroctonus Thorell, 1876and Belisarius Simon, 1879. Four-cluster likelihood mapping revealed that the instability inherent to the placement of Hadrurus, Uroctonus and Belisarius was attributable to significant gene tree conflict in the internodes corresponding to their divergences. To redress the polyphyly of Caraboctonidae, the following systematic actions have been taken: (1) the family Caraboctonidae has been delimited to consist of 23 species in the genera Caraboctonus and Hadruroides Pocock, 1893; (2) Caraboctonidae, previously included in the superfamily Iuroidea Thorell, 1876 or as incertae sedis, is transferred to the superfamily Caraboctonoidea (new rank); (3) the superfamily Hadruroidea (new rank) is established and the status of Hadrurinae Stahnke, 1973 is elevated to family (Hadruridae new status) including 9 species in the genera Hadrurus and Hoffmannihadrurus Fet & Soleglad, 2004 and (4) we treat Uroctonus and Belisarius as insertae sedis with respect to superfamilial placement. Our systematic actions engender the monophyly of both Iuroidea and Caraboctonidae. Future phylogenomic investigations should target similar taxon-poor and understudied lineages of potential phylogenetic significance, which are anticipated to reveal additional non-monophyletic groups.

Hemipenial morphology of nine South American species of Mabuya (Scincidae: Lygosominae) with comments on the morphology of the family

Species of the genus Mabuya are barely recognizable because of their highly conservative external morphology, resulting in unstable systematic and taxonomic status. Considering that hemipenial morphology has been extensively used as a source of systematic information for several groups of Squamata, we described the hemipenial morphology of nine species of Mabuya, searching for diagnostic features useful to taxonomic and systematic studies. To have comparative parameters within the family, we also presented the description of the hemipenes of two other Scincidae species (i.e., Trachylepis atlantica and Tiliqua rugosa). This information, together with the available literature, was analyzed in a comparative framework, allowing us to define some variation for the genus and family. The hemipenial morphology of the species of Mabuya studied is very conservative; features that show subtle variation were the shape of the protuberances present on the asulcate face of the hemipenial body and the presence and type of a lateral fold, on the body. Also, within the family, some variation was found in the shape and length of lobes, length of branches of the sulcus spermaticus, and in the kind of folds and protuberances in both faces of the lobes and hemipenial body. Although these features were not tested in a systematic context, this information suggested variations in hemipenial morphology that could contribute to taxonomic diagnostics and phylogenetic characters. The importance of obtaining proper preparations for gathering accurate information is emphasized, and the use of fresh-killed specimens that shows fewer difficulties for the procedures of preparation is highly recommended.

Dark offshoot: Phylogenomic data sheds light on the evolutionary history of a new species of cave brittle star

Caves are a useful system for testing evolutionary and biogeographic hypotheses, as they are isolated, and their environmental conditions have resulted in adaptive selection across different taxa. Although in recent years many more cave species have been discovered, cave-dwelling members of the class Ophiuroidea (brittle stars) remain scarce. Out of the more than two thousand species of brittle stars described to date, only three are regarded as true cave-dwellers. These occurrences represent rare colonising events, compared to other groups that are known to have successfully diversified in these systems. A third species from an anchihaline cave system in the Yucatan Peninsula, Mexico, has been previously identified from cytochrome oxidase I (COI) barcodes. In this study, we reassess the species boundaries of this putative cave species using a phylogenomic dataset (20 specimens in 13 species, 100 exons, 18.7 kbp). We perform species delimitation analyses using robust full-coalescent methods for discovery and validation of hypotheses on species boundaries, as well as infer its phylogenetic relationships with species distributed in adjacent marine regions, in order to investigate the origin of this cave-adapted species. We assess which hypotheses on the origin of subterranean taxa can be applied to this species by taking into account its placement within the genus Ophionereis and its demographic history. We provide a detailed description of Ophionereis commutabilis n. sp., and evaluate its morphological characters in the light of its successful adaptation to life in caves.

Multilocus phylogeny and a new classification for African, Asian and Indian supple and writhing skinks (Scincidae: Lygosominae)

The genera Lepidothyris, Lygosoma and Mochlus comprise the writhing or supple skinks, a group of semi-fossorial, elongate-bodied skinks distributed across the Old World Tropics. Due to their generalized morphology and lack of diagnostic characters, species- and clade-level relationships have long been debated. Recent molecular phylogenetic studies of the group have provided some clarification of species-level relationships, but a number of issues regarding higher level relationships among genera still remain. Here we present a phylogenetic estimate of relationships among species in Lygosoma, Mochlus and Lepidothyris generated by concatenated and species tree analyses of multilocus data using the most extensive taxonomic sampling of the group to date. We also use multivariate statistics to examine species and clade distributions in morpho space. Our results reject the monophyly of Lygosoma s.l., Lygosoma s.s. and Mochlus, which highlights the instability of the current taxonomic classification of the group. We, therefore, revise the taxonomy of the writhing skinks to better reflect the evolutionary history of Lygosoma s.l. by restricting Lygosoma for Southeast Asia, resurrecting the genus Riopa for a clade of Indian and Southeast Asian species, expanding the genus Mochlus to include all African species of writhing skinks and describing a new genus in Southeast Asia.

Signal, Uncertainty, and Conflict in Phylogenomic Data for a Diverse Lineage of Microbial Eukaryotes (Diatoms, Bacillariophyta)

Diatoms (Bacillariophyta) are a species-rich group of eukaryotic microbes diverse in morphology, ecology, and metabolism. Previous reconstructions of the diatom phylogeny based on one or a few genes have resulted in inconsistent resolution or low support for critical nodes. We applied phylogenetic paralog pruning techniques to a data set of 94 diatom genomes and transcriptomes to infer perennially difficult species relationships, using concatenation and summary-coalescent methods to reconstruct species trees from data sets spanning a wide range of thresholds for taxon and column occupancy in gene alignments. Conflicts between gene and species trees decreased with both increasing taxon occupancy and bootstrap cutoffs applied to gene trees. Concordance between gene and species trees was lowest for short internodes and increased logarithmically with increasing edge length, suggesting that incomplete lineage sorting disproportionately affects species tree inference at short internodes, which are a common feature of the diatom phylogeny. Although species tree topologies were largely consistent across many data treatments, concatenation methods appeared to outperform summary-coalescent methods for sparse alignments. Our results underscore that approaches to species-tree inference based on few loci are likely to be misled by unrepresentative sampling of gene histories, particularly in lineages that may have diversified rapidly. In addition, phylogenomic studies of diatoms, and potentially other hyperdiverse groups, should maximize the number of gene trees with high taxon occupancy, though there is clearly a limit to how many of these genes will be available.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof et al. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

StarBEAST2 Brings Faster Species Tree Inference and Accurate Estimates of Substitution Rates

Fully Bayesian multispecies coalescent (MSC) methods like *BEAST estimate species trees from multiple sequence alignments. Today thousands of genes can be sequenced for a given study, but using that many genes with *BEAST is intractably slow. An alternative is to use heuristic methods which compromise accuracy or completeness in return for speed. A common heuristic is concatenation, which assumes that the evolutionary history of each gene tree is identical to the species tree. This is an inconsistent estimator of species tree topology, a worse estimator of divergence times, and induces spurious substitution rate variation when incomplete lineage sorting is present. Another class of heuristics directly motivated by the MSC avoids many of the pitfalls of concatenation but cannot be used to estimate divergence times. To enable fuller use of available data and more accurate inference of species tree topologies, divergence times, and substitution rates, we have developed a new version of *BEAST called StarBEAST2. To improve convergence rates we add analytical integration of population sizes, novel MCMC operators and other optimizations. Computational performance improved by 13.5× and 13.8× respectively when analyzing two empirical data sets, and an average of 33.1× across 30 simulated data sets. To enable accurate estimates of per-species substitution rates, we introduce species tree relaxed clocks, and show that StarBEAST2 is a more powerful and robust estimator of rate variation than concatenation. StarBEAST2 is available through the BEAUTi package manager in BEAST 2.4 and above.

RADseq dataset with 90% missing data fully resolves recent radiation of Petalidium (Acanthaceae) in the ultra-arid deserts of Namibia

Deserts, even those at tropical latitudes, often have strikingly low levels of plant diversity, particularly within genera. One remarkable exception to this pattern is the genus Petalidium (Acanthaceae), in which 37 of 40 named species occupy one of the driest environments on Earth, the Namib Desert of Namibia and neighboring Angola. To contribute to understanding this enigmatic diversity, we generated RADseq data for 47 accessions of Petalidium representing 22 species. We explored the impacts of 18 different combinations of assembly parameters in de novo assembly of the data across nine levels of missing data plus a best practice assembly using a reference Acanthaceae genome for a total of 171 sequence datasets assembled. RADseq data assembled at several thresholds of missing data, including 90% missing data, yielded phylogenetic hypotheses of Petalidium that were confidently and nearly fully resolved, which is notable given that divergence time analyses suggest a crown age for African species of 3.6-1.4 Ma. De novo assembly of our data yielded the most strongly supported and well-resolved topologies; in contrast, reference-based assembly performed poorly, perhaps due in part to moderate phylogenetic divergence between the reference genome, Ruellia speciosa, and the ingroup. Overall, we found that Petalidium, despite the harshness of the environment in which species occur, shows a net diversification rate (0.8-2.1 species per my) on par with those of diverse genera in tropical, Mediterranean, and alpine environments.

Comparative skull anatomy of terrestrial and crevice-dwelling Trachylepis skinks (Squamata: Scincidae) with a survey of resources in scincid cranial osteology

Skinks account for more than 25% of all lizard species; however, representatives of fewer than a quarter of all species have been characterized osteologically. All but a few of the available cranial descriptions concentrate solely on characters that can be seen externally on the intact skull. Mabuyid skinks of the genus Trachylepis are the dominant, fully limbed skinks in Sub-Saharan Africa, and nearly all species have the same generalized body plan. Although a few rock crevice-dwelling species possess slight body depression, extreme dorsoventral depression is observed only in Trachylepis laevis. We investigated the detailed skull anatomy of three Trachylepis skinks (T. laevis, T. sulcata, and T. gonwouoi, a recently described species allied to T. affinis) using high-resolution X-ray micro-computed tomography. Our goals were to review the scincid cranial osteology literature in a phylogenetic context, provide a detailed anatomical atlas for the mabuyid lineage, and investigate the morphological adaptations of the highly modified T. laevis. Our results demonstrate that there is significant morphological variation between these three taxa, including the loss and fusion of structures, as well as changes in the shape, scale, and relationship between individual elements. Trachylepis laevis possesses several osteological modifications that have produced a reducton in head depth that are likely functional consequences of extreme rupicolous habits, including a flat skull roof, many strongly recumbent elements, and a depressed neurocranium.We hypothesize these modifications may correspond to descreased bite force and increased capabilities of cranial kinesis. Our study is the first element-by-element description of a skink using computed tomography technology.

The phylogenetic relationships and molecular systematics of scincid lizards of the genus Heremites (Sauria, Scincidae) in the Middle East based on mtDNA sequences

The taxonomic status of species included in the genus Heremites in Iran and Iraq is uncertain. Three of these species have been assigned to the genus based on morphology: Heremites auratus transcaucasica, H. vittatus, and H. septemtaeniatus. We examined the phylogenetic relationships and taxonomic status of the Iranian and Iraqi species of Heremites by performing phylogenetic analyses using mitochondrial DNA sequences (cytochrome b and 16S rRNA). Phylogenetic relationships and estimated genetic distances indicated that the Heremites populations of the area (Iran and Iraq) form five distinct clades. Three of these clades are found only in Iran, specifically in: (1) Fars and Hormozgan provinces; (2) Northeastern Khuzestan; and (3) Khorasan and Isfahan provinces. The fourth clade (H. septemtaeniatus) is found in west and Mahshahr in Iran as well as in eastern and northern parts of Iraq. The fifth clade, Heremites vittatus, is found in Iran and Iraq. We also confirm the absence of H. auratus in Iran and Iraq. It also indicated that H. vittatus is sister taxon to the other groups that our analyses estimate the divergence of this clade in the Middle Miocene (15.9 Mya). The clade containing the Fars-Hormozgan and Khuzestan populations diverged at the end of the Miocene (8.5 Mya). The Isfahan and Khorasan populations separated at the Pliocene (4.2 Mya) from the western Iranian group, the group in Mahshahr, Iran and the groups in northern and eastern Iraq.

Macrosystematics of eutherian mammals combining HTS data to expand taxon coverage

In the last few years high-throughput sequencing technologies have permitted significant advances in mammalian phylogenetic studies from a genomic perspective. However, these studies have been restricted to a sparse number of species with available reference genomes. Thus, several issues inside the eutherian mammals phylogeny remain unresolved. This may be due in part to limited taxon sampling, as taxonomic density is known to affect phylogenetic resolution. In this context, we present a protocol to increase taxon coverage using high-throughput sequencing data (RNA or DNA) generated for other biological studies and available in public databases. Following this procedure we addressed pending or controversial issues concerning the phylogenetic position of Dermoptera, Pholidota and Chiroptera, considering multiple and independent loci. Also for Chiroptera and Arctoidea we evaluated the relationships of the lineages that compose it. Although the maximum number of genes used is moderate (95), in some cases taxon coverage doubles that of previous related studies. Globally, all coalescent-based (STAR, MP-EST and ASTRAL) and concatenated (IQ-TREE and BEAST2) methods used for species tree reconstruction were consistent to each other and most of interrogated nodes received high statistical support.

A Species-Level Phylogeny of Extant Snakes with Description of a New Colubrid Subfamily and Genus

BACKGROUND

With over 3,500 species encompassing a diverse range of morphologies and ecologies, snakes make up 36% of squamate diversity. Despite several attempts at estimating higher-level snake relationships and numerous assessments of generic- or species-level phylogenies, a large-scale species-level phylogeny solely focusing on snakes has not been completed. Here, we provide the largest-yet estimate of the snake tree of life using maximum likelihood on a supermatrix of 1745 taxa (1652 snake species + 7 outgroup taxa) and 9,523 base pairs from 10 loci (5 nuclear, 5 mitochondrial), including previously unsequenced genera (2) and species (61).

RESULTS

Increased taxon sampling resulted in a phylogeny with a new higher-level topology and corroborate many lower-level relationships, strengthened by high nodal support values (> 85%) down to the species level (73.69% of nodes). Although the majority of families and subfamilies were strongly supported as monophyletic with > 88% support values, some families and numerous genera were paraphyletic, primarily due to limited taxon and loci sampling leading to a sparse supermatrix and minimal sequence overlap between some closely-related taxa. With all rogue taxa and incertae sedis species eliminated, higher-level relationships and support values remained relatively unchanged, except in five problematic clades.

CONCLUSION

Our analyses resulted in new topologies at higher- and lower-levels; resolved several previous topological issues; established novel paraphyletic affiliations; designated a new subfamily, Ahaetuliinae, for the genera Ahaetulla, Chrysopelea, Dendrelaphis, and Dryophiops; and appointed Hemerophis (Coluber) zebrinus to a new genus, Mopanveldophis. Although we provide insight into some distinguished problematic nodes, at the deeper phylogenetic scale, resolution of these nodes may require sampling of more slowly-evolving nuclear genes.

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.

Molecular phylogeny of Panaspis and Afroablepharus skinks (Squamata: Scincidae) in the savannas of sub-Saharan Africa

African snake-eyed skinks are relatively small lizards of the genera Panaspis and Afroablepharus. Species allocation of these genera frequently changed during the 20th century based on morphology, ecology, and biogeography. Members of these genera occur primarily in savanna habitats throughout sub-Saharan Africa and include species whose highly conserved morphology poses challenges for taxonomic studies. We sequenced two mitochondrial (16S and cyt b) and two nuclear genes (PDC and RAG1) from 76 Panaspis and Afroablepharus samples from across eastern, central, and southern Africa. Concatenated gene-tree and divergence-dating analyses were conducted to infer phylogenies and biogeographic patterns. Molecular data sets revealed several cryptic lineages, with most radiations occurring during the mid-Miocene to Pliocene. We infer that rifting processes (including the formation of the East African Rift System) and climatic oscillations contributed to the expansion and contraction of savannas, and caused cladogenesis in snake-eyed skinks. Species in Panaspis and Afroablepharus used in this study, including type species for both genera, formed a monophyletic group. As a result, the latter genus should be synonymized with the former, which has priority. Conservatively, we continue to include the West African species P. breviceps and P. togoensis within an expanded Panaspis, but note that they occur in relatively divergent clades, and their taxonomic status may change with improved taxon sampling. Divergence estimates and cryptic speciation patterns of snake-eyed skinks were consistent with previous studies of other savanna vertebrate lineages from the same areas examined in this study.

Phylogenomic analyses reveal novel relationships among snake families

Snakes are a diverse and important group of vertebrates. However, relationships among the major groups of snakes have remained highly uncertain, with recent studies hypothesizing very different (and typically weakly supported) relationships. Here, we address family-level snake relationships with new phylogenomic data from 3776 nuclear loci from ultraconserved elements (1.40million aligned base pairs, 52% missing data overall) sampled from 29 snake species that together represent almost all families, a dataset ∼100 times larger than used in previous studies. We found relatively strong support from species-tree analyses (NJst) for most relationships, including three largely novel clades: (1) a clade uniting the boas, pythons and their relatives, (2) a clade placing cylindrophiids and uropeltids with this clade, and (3) a clade uniting bolyeriids (Round Island boas) with pythonids and their relatives (xenopeltids and loxocemids). Relationships among families of advanced snakes (caenophidians) were also strongly supported. The results show the potential for phylogenomic analyses to resolve difficult groups, but also show a surprising sensitivity of the analyses to the inclusion or exclusion of outgroups.

Resolving evolutionary relationships in lichen-forming fungi using diverse phylogenomic datasets and analytical approaches

Evolutionary histories are now being inferred from unprecedented, genome-scale datasets for a broad range of organismal groups. While phylogenomic data has helped in resolving a number of difficult, long-standing questions, constructing appropriate datasets from genomes is not straightforward, particularly in non-model groups. Here we explore the utility of phylogenomic data to infer robust phylogenies for a lineage of closely related lichen-forming fungal species. We assembled multiple, distinct nuclear phylogenomic datasets, ranging from ca. 25 Kb to 16.8 Mb and inferred topologies using both concatenated gene tree approaches and species tree methods based on the multispecies coalescent model. In spite of evidence for rampant incongruence among individual loci, these genome-scale datasets provide a consistent, well-supported phylogenetic hypothesis using both concatenation and multispecies coalescent approaches (ASTRAL-II and SVDquartets). However, the popular full hierarchical coalescent approach implemented in *BEAST provided inconsistent inferences, both in terms of nodal support and topology, with smaller subsets of the phylogenomic data. While comparable, well-supported topologies can be accurately inferred with only a small fraction of the overall genome, consistent results across a variety of datasets and methodological approaches provide reassurance that phylogenomic data can effectively be used to provide robust phylogenies for closely related lichen-forming fungal lineages.

Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies

Phylogenetic resolution of ancient rapid radiations has remained problematic despite major advances in statistical approaches and DNA sequencing technologies. Here we report on a combined phylogenetic approach utilizing transcriptome data in conjunction with Sanger sequence data to investigate a tandem of ancient divergences in the harvestmen superfamily Ischyropsalidoidea (Arachnida, Opiliones, Dyspnoi). We rely on Sanger sequences to resolve nodes within and between closely related genera, and use RNA-seq data from a subset of taxa to resolve a short and ancient internal branch. We use several analytical approaches to explore this succession of ancient diversification events, including concatenated and coalescent-based analyses and maximum likelihood gene trees for each locus. We evaluate the robustness of phylogenetic inferences using a randomized locus sub-sampling approach, and find congruence across these methods despite considerable incongruence across gene trees. Incongruent gene trees are not recovered in frequencies expected from a simple multispecies coalescent model, and we reject incomplete lineage sorting as the sole contributor to gene tree conflict. Using these approaches we attain robust support for higher-level phylogenetic relationships within Ischyropsalidoidea.

Comparing species tree estimation with large anchored phylogenomic and small Sanger-sequenced molecular datasets: an empirical study on Malagasy pseudoxyrhophiine snakes

BACKGROUND

Using molecular data generated by high throughput next generation sequencing (NGS) platforms to infer phylogeny is becoming common as costs go down and the ability to capture loci from across the genome goes up. While there is a general consensus that greater numbers of independent loci should result in more robust phylogenetic estimates, few studies have compared phylogenies resulting from smaller datasets for commonly used genetic markers with the large datasets captured using NGS. Here, we determine how a 5-locus Sanger dataset compares with a 377-locus anchored genomics dataset for understanding the evolutionary history of the pseudoxyrhophiine snake radiation centered in Madagascar. The Pseudoxyrhophiinae comprise ~86 % of Madagascar's serpent diversity, yet they are poorly known with respect to ecology, behavior, and systematics. Using the 377-locus NGS dataset and the summary statistics species-tree methods STAR and MP-EST, we estimated a well-supported species tree that provides new insights concerning intergeneric relationships for the pseudoxyrhophiines. We also compared how these and other methods performed with respect to estimating tree topology using datasets with varying numbers of loci.

METHODS

Using Sanger sequencing and an anchored phylogenomics approach, we sequenced datasets comprised of 5 and 377 loci, respectively, for 23 pseudoxyrhophiine taxa. For each dataset, we estimated phylogenies using both gene-tree (concatenation) and species-tree (STAR, MP-EST) approaches. We determined the similarity of resulting tree topologies from the different datasets using Robinson-Foulds distances. In addition, we examined how subsets of these data performed compared to the complete Sanger and anchored datasets for phylogenetic accuracy using the same tree inference methodologies, as well as the program *BEAST to determine if a full coalescent model for species tree estimation could generate robust results with fewer loci compared to the summary statistics species tree approaches. We also examined the individual gene trees in comparison to the 377-locus species tree using the program MetaTree.

RESULTS

Using the full anchored dataset under a variety of methods gave us the same, well-supported phylogeny for pseudoxyrhophiines. The African pseudoxyrhophiine Duberria is the sister taxon to the Malagasy pseudoxyrhophiines genera, providing evidence for a monophyletic radiation in Madagascar. In addition, within Madagascar, the two major clades inferred correspond largely to the aglyphous and opisthoglyphous genera, suggesting that feeding specializations associated with tooth venom delivery may have played a major role in the early diversification of this radiation. The comparison of tree topologies from the concatenated and species-tree methods using different datasets indicated the 5-locus dataset cannot beused to infer a correct phylogeny for the pseudoxyrhophiines under any method tested here and that summary statistics methods require 50 or more loci to consistently recover the species-tree inferred using the complete anchored dataset. However, as few as 15 loci may infer the correct topology when using the full coalescent species tree method *BEAST. MetaTree analyses of each gene tree from the Sanger and anchored datasets found that none of the individual gene trees matched the 377-locus species tree, and that no gene trees were identical with respect to topology.

CONCLUSIONS

Our results suggest that ≥50 loci may be necessary to confidently infer phylogenies when using summaryspecies-tree methods, but that the coalescent-based method *BEAST consistently recovers the same topology using only 15 loci. These results reinforce that datasets with small numbers of markers may result in misleading topologies, and further, that the method of inference used to generate a phylogeny also has a major influence on the number of loci necessary to infer robust species trees.

The systematics of carnivorous sponges

Carnivorous sponges are characterized by their unique method of capturing mesoplanktonic prey coupled with the complete or partial reduction of the aquiferous system characteristic of the phylum Porifera. Current systematics place the vast majority of carnivorous sponges within Cladorhizidae, with certain species assigned to Guitarridae and Esperiopsidae. Morphological characters have not been able to show whether this classification is evolutionary accurate, and whether carnivory has evolved once or in several lineages. In the present paper we present the first comprehensive molecular phylogeny of the carnivorous sponges, interpret these results in conjunction with morphological characters, and propose a revised classification of the group. Molecular phylogenies were inferred using 18S rDNA and a combined dataset of partial 28S rDNA, COI and ALG11 sequences. The results recovered carnivorous sponges as a clade closely related to the families Mycalidae and Guitarridae, showing family Cladorhizidae to be monophyletic and also including carnivorous species currently placed in other families. The genus Lycopodina is resurrected for species currently placed in the paraphyletic subgenus Asbestopluma (Asbestopluma) featuring forceps spicules and lacking sigmas or sigmancistras. The genera Chondrocladia and Cladorhiza are found to be monophyletic. However, results indicate that the subgenus Chondrocladia is polyphyletic with respect to the subgenera Meliiderma and Symmetrocladia. Euchelipluma, formerly Guitarridae, is retained, but transferred to Cladorhizidae. The four known carnivorous species currently in Esperiopsis are transferred to Abyssocladia. Neocladia is a junior homonym and is here renamed Koltunicladia. Our results provide strong evidence in support of the hypothesis that carnivory in sponges has evolved only once. While spicule characters mostly reflect monophyletic groups at the generic level, differences between genera represent evolution within family Cladorhizidae rather than evolution of carnivory in separate lineages. Conflicting spicule characters can be reinterpreted to support the inclusion of all carnivorous sponges within Cladorhizidae, and a carnivorous habit should thus be considered the main diagnostic character in systematic classification.

How Should Genes and Taxa be Sampled for Phylogenomic Analyses with Missing Data? An Empirical Study in Iguanian Lizards

Targeted sequence capture is becoming a widespread tool for generating large phylogenomic data sets to address difficult phylogenetic problems. However, this methodology often generates data sets in which increasing the number of taxa and loci increases amounts of missing data. Thus, a fundamental (but still unresolved) question is whether sampling should be designed to maximize sampling of taxa or genes, or to minimize the inclusion of missing data cells. Here, we explore this question for an ancient, rapid radiation of lizards, the pleurodont iguanians. Pleurodonts include many well-known clades (e.g., anoles, basilisks, iguanas, and spiny lizards) but relationships among families have proven difficult to resolve strongly and consistently using traditional sequencing approaches. We generated up to 4921 ultraconserved elements with sampling strategies including 16, 29, and 44 taxa, from 1179 to approximately 2.4 million characters per matrix and approximately 30% to 60% total missing data. We then compared mean branch support for interfamilial relationships under these 15 different sampling strategies for both concatenated (maximum likelihood) and species tree (NJst) approaches (after showing that mean branch support appears to be related to accuracy). We found that both approaches had the highest support when including loci with up to 50% missing taxa (matrices with ~40-55% missing data overall). Thus, our results show that simply excluding all missing data may be highly problematic as the primary guiding principle for the inclusion or exclusion of taxa and genes. The optimal strategy was somewhat different for each approach, a pattern that has not been shown previously. For concatenated analyses, branch support was maximized when including many taxa (44) but fewer characters (1.1 million). For species-tree analyses, branch support was maximized with minimal taxon sampling (16) but many loci (4789 of 4921). We also show that the choice of these sampling strategies can be critically important for phylogenomic analyses, since some strategies lead to demonstrably incorrect inferences (using the same method) that have strong statistical support. Our preferred estimate provides strong support for most interfamilial relationships in this important but phylogenetically challenging group.

Multi-locus fossil-calibrated phylogeny of Atheriniformes (Teleostei, Ovalentaria)

Phylogenetic relationships among families within the order Atheriniformes have been difficult to resolve on the basis of morphological evidence. Molecular studies so far have been fragmentary and based on a small number taxa and loci. In this study, we provide a new phylogenetic hypothesis based on sequence data collected for eight molecular markers for a representative sample of 103 atheriniform species, covering 2/3 of the genera in this order. The phylogeny is calibrated with six carefully chosen fossil taxa to provide an explicit timeframe for the diversification of this group. Our results support the subdivision of Atheriniformes into two suborders (Atherinopsoidei and Atherinoidei), the nesting of Notocheirinae within Atherinopsidae, and the monophyly of tribe Menidiini, among others. We propose taxonomic changes for Atherinopsoidei, but a few weakly supported nodes in our phylogeny suggests that further study is necessary to support a revised taxonomy of Atherinoidei. The time-calibrated phylogeny was used to infer ancestral habitat reconstructions to explain the current distribution of marine and freshwater taxa. Based on these results, the current distribution of Atheriniformes is likely due to widespread marine dispersal along the margins of continents, infrequent trans-oceanic dispersal, and repeated invasion of freshwater habitats. This conclusion is supported by post-Gondwanan divergence times among families within the order, and a high probability of a marine ancestral habitat.