Evaluating methods for phylogenomic analyses, and a new phylogeny for a major frog clade (Hyloidea) based on 2214 loci

Diversity and Molecular Evolution of Nonvisual Opsin Genes across Environmental, Developmental, and Morphological Adaptations in Frogs

Nonvisual opsins are transmembrane proteins expressed in the eyes and other tissues of many animals. When paired with a light-sensitive chromophore, nonvisual opsins form photopigments involved in various nonvisual, light-detection functions including circadian rhythm regulation, light-seeking behaviors, and seasonal responses. Here, we investigate the molecular evolution of nonvisual opsin genes in anuran amphibians (frogs and toads). We test several evolutionary hypotheses including the predicted loss of nonvisual opsins due to nocturnal ancestry and potential functional differences in nonvisual opsins resulting from environmental light variation across diverse anuran ecologies. Using whole-eye transcriptomes of 81 species, combined with genomes, multitissue transcriptomes, and independently annotated genes from an additional 21 species, we identify which nonvisual opsins are present in anuran genomes and those that are also expressed in the eyes, compare selective constraint among genes, and test for potential adaptive evolution by comparing selection between discrete ecological classes. At the genomic level, we recovered all 18 ancestral vertebrate nonvisual opsins, indicating that anurans demonstrate the lowest documented amount of opsin gene loss among ancestrally nocturnal tetrapods. We consistently found expression of 14 nonvisual opsins in anuran eyes and detected positive selection in a subset of these genes. We also found shifts in selective constraint acting on nonvisual opsins in frogs with differing activity periods, habitats, distributions, life histories, and pupil shapes, which may reflect functional adaptation. Although many nonvisual opsins remain poorly understood, these findings provide insight into the diversity and evolution of these genes across anurans, filling an important gap in our understanding of vertebrate opsins and setting the stage for future research on their functional evolution across taxa.

Populating a Continent: Phylogenomics Reveal the Timing of Australian Frog Diversification

The Australian continent's size and isolation make it an ideal place for studying the accumulation and evolution of biodiversity. Long separated from the ancient supercontinent Gondwana, most of Australia's plants and animals are unique and endemic, including the continent's frogs. Australian frogs comprise a remarkable ecological and morphological diversity categorized into a small number of distantly related radiations. We present a phylogenomic hypothesis based on an exon-capture dataset that spans the main clades of Australian myobatrachoid, pelodryadid hyloid, and microhylid frogs. Our time-calibrated phylogenomic-scale phylogeny identifies great disparity in the relative ages of these groups that vary from Gondwanan relics to recent immigrants from Asia and include arguably the continent's oldest living vertebrate radiation. This age stratification provides insight into the colonization of, and diversification on, the Australian continent through deep time, during periods of dramatic climatic and community changes. Contemporary Australian frog diversity highlights the adaptive capacity of anurans, particularly in response to heat and aridity, and explains why they are one of the continent's most visible faunas. [Anuran; adaptive radiation; Gondwana; phylogenetics].

Shape Evolution in Two Acts: Morphological Diversity of Larval and Adult Neoaustraranan Frogs

Phenotypic traits can evolve independently at different stages of ontogeny, optimizing adaptation to distinct ecological contexts and increasing morphological diversity in species with complex life cycles. Given the relative independence resulting from the profound changes induced by metamorphosis, niche occupation and resource utilization in tadpoles may prompt evolutionary responses that do not necessarily affect the adults. Consequently, diversity patterns observed in the larval shape may not necessarily correspond to those found in the adult shape for the same species, a premise that can be tested through the Adaptive Decoupling Hypothesis (ADH). Herein, we investigate the ADH for larval and adult shape differentiation in Neoaustrarana frogs. Neoaustrarana frogs, particularly within the Cycloramphidae family, exhibit remarkable diversity in tadpole morphology, making them an ideal model for studying adaptive decoupling. By analyzing 83 representative species across four families (Alsodidae, Batrachylidae, Cycloramphidae, and Hylodidae), we generate a morphological dataset for both larval and adult forms. We found a low correlation between larval and adult shapes, species with a highly distinct larval shape having relatively similar shape when adults. Larval morphological disparity is not a good predictor for adult morphological disparity within the group, with distinct patterns observed among families. Differences between families are notable in other aspects as well, such as the role of allometric components influencing shape and morphospace occupancy. The larval shape has higher phylogenetic structure than the adult. Evolutionary convergence emerges as a mechanism of diversification for both larval and adult shapes in the early evolution of neoaustraranans, with shape disparity of tadpoles reaching stable levels since the Oligocene. The widest occupation in morphospace involves families associated with dynamically changing environments over geological time. Our findings support the ADH driving phenotypic diversity in Neoaustrarana, underscoring the importance of considering ontogenetic stages in evolutionary studies.

Diversity and Evolution of Frog Visual Opsins: Spectral Tuning and Adaptation to Distinct Light Environments

Visual systems adapt to different light environments through several avenues including optical changes to the eye and neurological changes in how light signals are processed and interpreted. Spectral sensitivity can evolve via changes to visual pigments housed in the retinal photoreceptors through gene duplication and loss, differential and coexpression, and sequence evolution. Frogs provide an excellent, yet understudied, system for visual evolution research due to their diversity of ecologies (including biphasic aquatic-terrestrial life cycles) that we hypothesize imposed different selective pressures leading to adaptive evolution of the visual system, notably the opsins that encode the protein component of the visual pigments responsible for the first step in visual perception. Here, we analyze the diversity and evolution of visual opsin genes from 93 new eye transcriptomes plus published data for a combined dataset spanning 122 frog species and 34 families. We find that most species express the four visual opsins previously identified in frogs but show evidence for gene loss in two lineages. Further, we present evidence of positive selection in three opsins and shifts in selective pressures associated with differences in habitat and life history, but not activity pattern. We identify substantial novel variation in the visual opsins and, using microspectrophotometry, find highly variable spectral sensitivities, expanding known ranges for all frog visual pigments. Mutations at spectral-tuning sites only partially account for this variation, suggesting that frogs have used tuning pathways that are unique among vertebrates. These results support the hypothesis of adaptive evolution in photoreceptor physiology across the frog tree of life in response to varying environmental and ecological factors and further our growing understanding of vertebrate visual evolution.

Relicts in the mist: Two new frog families, genera and species highlight the role of Pantepui as a biodiversity museum throughout the Cenozoic

The iconic mountains of the Pantepui biogeographical region host many early-diverging endemic animal and plant lineages, concurring with Conan Doyle's novel about an ancient "Lost World". While this is the case of several frog lineages, others appear to have more recent origins, adding to the controversy around the diversification processes in this region. Due to its remoteness, Pantepui is challenging for biological surveys, and only a glimpse of its biodiversity has been described, which hampers comprehensive evolutionary studies in many groups. During a recent expedition to the Neblina massif on the Brazil-Venezuela border, we sampled two new frog species that could not be assigned to any known genus. Here, we perform phylogenetic analyses of mitogenomic and nuclear loci to infer the evolutionary relationships of the new taxa and support their description. We find that both species represent single lineages deeply nested within Brachycephaloidea, a major Neotropical clade of direct-developing frogs. Both species diverged >45 Ma from their closest relatives: the first is sister to all other Brachycephaloidea except for Ceuthomantis, another Pantepui endemic, and the second is sister to Brachycephalidae, endemic to the Brazilian Atlantic Forest. In addition to these considerable phylogenetic and biogeographic divergences, external morphology and osteological features support the proposition of two new family and genus-level taxa to accommodate these new branches of the amphibian tree of life. These findings add to other recently described ancient vertebrate lineages from the Neblina massif, providing a bewildering reminder that our perception of the Pantepui's biodiversity remains vastly incomplete. It also provides insights into how these mountains acted as "museums" during the diversification of Brachycephaloidea and of Neotropical biotas more broadly, in line with the influential "Plateau theory". Finally, these discoveries point at the yet unknown branches of the tree of life that may go extinct, due to global climate change and zoonotic diseases, before we even learn about their existence, amphibians living at higher elevations being particularly at risk.

Frog phylogeny: A time-calibrated, species-level tree based on hundreds of loci and 5,242 species

Large-scale, time-calibrated phylogenies from supermatrix studies have become crucial for evolutionary and ecological studies in many groups of organisms. However, in frogs (anuran amphibians), there is a serious problem with existing supermatrix estimates. Specifically, these trees are based on a limited number of loci (15 or fewer), and the higher-level relationships estimated are discordant with recent phylogenomic estimates based on much larger numbers of loci. Here, we attempted to rectify this problem by generating an expanded supermatrix and combining this with data from phylogenomic studies. To assist in aligning ribosomal sequences for this supermatrix, we developed a new program (TaxonomyAlign) to help perform taxonomy-guided alignments. The new combined matrix contained 5,242 anuran species with data from 307 markers, but with 95% missing data overall. This dataset represented a 71% increase in species sampled relative to the previous largest supermatrix analysis of anurans (adding 2,175 species). Maximum-likelihood analyses generated a tree in which higher-level relationships (and estimated clade ages) were generally concordant with those from phylogenomic analyses but were more discordant with the previous largest supermatrix analysis. We found few obvious problems arising from the extensive missing data in most species. We also generated a set of 100 time-calibrated trees for use in comparative analyses. Overall, we provide an improved estimate of anuran phylogeny based on the largest number of combined taxa and markers to date. More broadly, we demonstrate the potential to combine phylogenomic and supermatrix analyses in other groups of organisms.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Evolutionary and Genomic Diversity of True Polyploidy in Tetrapods

True polyploid organisms have more than two chromosome sets in their somatic and germline cells. Polyploidy is a major evolutionary force and has played a significant role in the early genomic evolution of plants, different invertebrate taxa, chordates, and teleosts. However, the contribution of polyploidy to the generation of new genomic, ecological, and species diversity in tetrapods has traditionally been underestimated. Indeed, polyploidy represents an important pathway of genomic evolution, occurring in most higher-taxa tetrapods and displaying a variety of different forms, genomic configurations, and biological implications. Herein, we report and discuss the available information on the different origins and evolutionary and ecological significance of true polyploidy in tetrapods. Among the main tetrapod lineages, modern amphibians have an unparalleled diversity of polyploids and, until recently, they were considered to be the only vertebrates with closely related diploid and polyploid bisexual species or populations. In reptiles, polyploidy was thought to be restricted to squamates and associated with parthenogenesis. In birds and mammals, true polyploidy has generally been considered absent (non-tolerated). These views are being changed due to an accumulation of new data, and the impact as well as the different evolutionary and ecological implications of polyploidy in tetrapods, deserve a broader evaluation.

Glassfrogs conceal blood in their liver to maintain transparency

Transparency in animals is a complex form of camouflage involving mechanisms that reduce light scattering and absorption throughout the organism. In vertebrates, attaining transparency is difficult because their circulatory system is full of red blood cells (RBCs) that strongly attenuate light. Here, we document how glassfrogs overcome this challenge by concealing these cells from view. Using photoacoustic imaging to track RBCs in vivo, we show that resting glassfrogs increase transparency two- to threefold by removing ~89% of their RBCs from circulation and packing them within their liver. Vertebrate transparency thus requires both see-through tissues and active mechanisms that "clear" respiratory pigments from these tissues. Furthermore, glassfrogs' ability to regulate the location, density, and packing of RBCs without clotting offers insight in metabolic, hemodynamic, and blood-clot research.

A new genus and species of frog from the Kem Kem (Morocco), the second neobatrachian from Cretaceous Africa

Neobatrachia, a clade representing the majority of extant anuran diversity, is thought to have emerged and diversified during the Cretaceous. Most of the early diversification of neobatrachians occurred in southern Gondwana, especially the regions that are today South America and Africa. Whereas five extinct neobatrachians have been described from the Cretaceous of South America in the last decade, only one is known from Africa. This difference in the known extinct diversity is linked to the lack of well-preserved specimens, understudy of fragmentary remains, and lack of known Cretaceous sites in Africa. Study of fragmentary anuran remains from Africa could allow for the identification of previously unknown neobatrachians, allowing for a better understanding of their early diversification. We reanalysed several previously described anuran specimens from the well-known Kem Kem beds, including using CT-scanning. Through our osteological study, we determined that several cranial bones and vertebrae represent a new hyperossified taxon for which we provide a formal description. Comparison to other hyperossified anurans revealed similarities and affinity of this new taxon with the neobatrachians Beelzebufo (extinct) and Ceratophrys (extant). Phylogenetic analyses supported this affinity, placing the new taxon within Neobatrachia in an unresolved clade of Ceratophryidae. This taxon is the oldest neobatrachian from Africa, and reveals that neobatrachians were already widespread throughout southern Gondwana during the earliest Late Cretaceous.

FrogCap: A modular sequence capture probe-set for phylogenomics and population genetics for all frogs, assessed across multiple phylogenetic scales

Despite the prevalence of high-throughput sequencing in phylogenetics, many relationships remain difficult to resolve because of conflicting signal among genomic regions. Selection of different types of molecular markers from different genomic regions is required to overcome these challenges. For evolutionary studies in frogs, we introduce the publicly available FrogCap suite of genomic resources, which is a large collection of ~15,000 markers that unifies previous genetic sequencing efforts. FrogCap is designed to be modular, such that subsets of markers and SNPs can be selected based on the desired phylogenetic scale. FrogCap uses a variety of marker types that include exons and introns, ultraconserved elements, and previously sequenced Sanger markers, which span up to 10,000 bp in alignment lengths; in addition, we demonstrate potential for SNP-based analyses. We tested FrogCap using 121 samples distributed across five phylogenetic scales, comparing probes designed using a consensus- or exemplar genome-based approach. Using the consensus design is more resilient to issues with sensitivity, specificity, and missing data than picking an exemplar genome sequence. We also tested the impact of different bait kit sizes (20,020 vs. 40,040) on depth of coverage and found triple the depth for the 20,020 bait  kit. We observed sequence capture success (i.e., missing data, sequenced markers/bases, marker length, and informative sites) across phylogenetic scales. The incorporation of different marker types is effective for deep phylogenetic relationships and shallow population genetics studies. Having demonstrated FrogCap's utility and modularity, we conclude that these new resources are efficacious for high-throughput sequencing projects across variable timescales.

Evolutionary analyses of visual opsin genes in frogs and toads: Diversity, duplication, and positive selection

Among major vertebrate groups, anurans (frogs and toads) are understudied with regard to their visual systems, and little is known about variation among species that differ in ecology. We sampled North American anurans representing diverse evolutionary and life histories that likely possess visual systems adapted to meet different ecological needs. Using standard molecular techniques, visual opsin genes, which encode the protein component of visual pigments, were obtained from anuran retinas. Additionally, we extracted the visual opsins from publicly available genome and transcriptome assemblies, further increasing the phylogenetic and ecological diversity of our dataset to 33 species in total. We found that anurans consistently express four visual opsin genes (RH1, LWS, SWS1, and SWS2, but not RH2) even though reported photoreceptor complements vary widely among species. The proteins encoded by these genes showed considerable sequence variation among species, including at sites known to shift the spectral sensitivity of visual pigments in other vertebrates and had conserved substitutions that may be related to dim-light adaptation. Using molecular evolutionary analyses of selection (dN/dS) we found significant evidence for positive selection at a subset of sites in the dim-light rod opsin gene RH1 and the long wavelength sensitive cone opsin LWS. The function of sites inferred to be under positive selection are largely unknown, but a few are likely to affect spectral sensitivity and other visual pigment functions based on proximity to previously identified sites in other vertebrates. We also found the first evidence of visual opsin duplication in an amphibian with the duplication of the LWS gene in the African bullfrog, which had distinct LWS copies on the sex chromosomes suggesting the possibility of sex-specific visual adaptation. Taken together, our results indicate that ecological factors, such as habitat and life history, as well as behavior, may be driving changes to anuran visual systems.

Re-evaluating the morphological evidence for the re-evolution of lost mandibular teeth in frogs

Dollo's law of irreversibility states that once a complex structure is lost, it cannot be regained in the same form. Several putative exceptions to Dollo's law have been identified using phylogenetic comparative methods, but the anatomy and development of these traits are often poorly understood. Gastrotheca guentheri is renowned as the only frog with teeth on the lower jaw. Mandibular teeth were lost in the ancestor of frogs more than 200 million years ago and subsequently regained in G. guentheri. Little is known about the teeth in this species despite being a frequent example of trait “re‐evolution,” leaving open the possibility that it may have mandibular pseudoteeth. We assessed the dental anatomy of G. guentheri using micro‐computed tomography and histology and confirmed the longstanding assumption that true mandibular teeth are present. Remarkably, the mandibular teeth of G. guentheri are nearly identical in gross morphology and development to upper jaw teeth in closely related species. The developmental genetics of tooth formation are unknown in this possibly extinct species. Our results suggest that an ancestral odontogenic pathway has been conserved but suppressed in the lower jaw since the origin of frogs, providing a possible mechanism underlying the re‐evolution of lost mandibular teeth.

Evolution of the strikingly diverse submandibular muscles in Anura

The most ventral muscles of the head (the mm. submentalis, intermandibularis, and interhyoideus) provide support to the gular region and lift the buccal floor during ventilation and feeding. These muscles show limited variation in most gnathostomes, but in Anura they exhibit a surprising diversity. The few studies that have explored this character system highlighted its potential as a source of phylogenetic information. In this paper we explored the diversity of this character system studying specimens of 567 anuran species and reviewing published data to cover a total of 1321 species, belonging to 53 of the 54 currently recognized anuran families, as well as caudates and caecilians. We defined 27 discrete characters including the number of muscle bellies, supplementary layers, hypertrophy and diversity of elastic fibres, and pigmentation, among others, and optimized them on a comprehensive phylogenetic hypothesis. We recognized 223 unambiguously optimized synapomorphies for numerous clades on different scales, including three for Anura and many for suprafamiliar clades with poor phenotypic support. Finally, we discussed the evolution of this highly diverse character system, including homology, development, and its functional role in vocalization and feeding. Interestingly, the striking levels of variation in some structures contrast with the amount of phylogenetic inertia, allowing us to recognize several general patterns. Supplementary elements of the m. intermandibularis evolved first as broad layers occuring in more than half of extant anuran species and then concentrated forming discreet bellies in several clades. The anterior portion of the gular region is not sexually dimorphic, and is likely related to ventilation and tongue protraction. Conversely, the diversity of the m. interhyoideus is strongly linked to vocal sacs, which are present only in adult males, suggesting the presence of two independent modules.

Aposematism facilitates the diversification of parental care strategies in poison frogs

Many organisms have evolved adaptations to increase the odds of survival of their offspring. Parental care has evolved several times in animals including ectotherms. In amphibians, ~ 10% of species exhibit parental care. Among these, poison frogs (Dendrobatidae) are well-known for their extensive care, which includes egg guarding, larval transport, and specialized tadpole provisioning with trophic eggs. At least one third of dendrobatids displaying aposematism by exhibiting warning coloration that informs potential predators about the presence of defensive skin toxins. Aposematism has a central role in poison frog diversification, including diet specialization, and visual and acoustic communication; and it is thought to have impacted their reproductive biology as well. We tested the latter association using multivariate phylogenetic methods at the family level. Our results show complex relationships between aposematism and certain aspects of the reproductive biology in dendrobatids. In particular, aposematic species tend to use more specialized tadpole-deposition sites, such as phytotelmata, and ferry fewer tadpoles than non-aposematic species. We propose that aposematism may have facilitated the diversification of microhabitat use in dendrobatids in the context of reproduction. Furthermore, the use of resource-limited tadpole-deposition environments may have evolved in tandem with an optimal reproductive strategy characterized by few offspring, biparental care, and female provisioning of food in the form of unfertilized eggs. We also found that in phytotelm-breeders, the rate of transition from cryptic to aposematic phenotype is 17 to 19 times higher than vice versa. Therefore, we infer that the aposematism in dendrobatids might serve as an umbrella trait for the evolution and maintenance of their complex offspring-caring activities.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Do host habitat use and cospeciation matter in the evolution of Oswaldocruzia (Nematoda, Molineidae) from neotropical amphibians?

The genus Oswaldocruzia Travassos, 1917 includes approximately 90 species that are parasitic on amphibians and reptiles around the world, of which 43 occur in the neotropical region. However, molecular data supporting the taxonomic status of most species of the genus are scarce, and their phylogenetic relationships remain unknown. Using specimens of four molineid taxa (Oswaldocruzia belenensis Santos, Giese, Maldonado Jr. and Lanfredi, 2008; Oswaldocruzia chabaudi Ben Slimane & Durette-Desset, 1996, Oswaldocruzia chambrieri Ben Slimane & Durette-Desset, 1996 and Kentropyxia hylae Feitosa, Furtado, Santos and Melo, 2015) from amphibian hosts collected in different regions of Pará, Brazil, we conducted morphological studies, molecular analyses and phylogenies (maximum likelihood and Bayesian inference) using the cytochrome c oxidase subunity I (Cox1) gene. The newly generated sequences were compared with those of ten publicly available Cox1 sequences of Oswaldocruzia from Mexico. Our findings demonstrated significant differences between the sequences of amazonian specimens and sequences from specimens collected in Mexico, and we suggest that host-parasite cospeciation or habitat use might be related to molineid evolution in amphibian hosts. Additionally, this work presents new hosts and new geographical records for species of Oswaldocruzia from the neotropics.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

Gene Tree Estimation Error with Ultraconserved Elements: An Empirical Study on Pseudapis Bees

Summarizing individual gene trees to species phylogenies using two-step coalescent methods is now a standard strategy in the field of phylogenomics. However, practical implementations of summary methods suffer from gene tree estimation error, which is caused by various biological and analytical factors. Greatly understudied is the choice of gene tree inference method and downstream effects on species tree estimation for empirical data sets. To better understand the impact of this method choice on gene and species tree accuracy, we compare gene trees estimated through four widely used programs under different model-selection criteria: PhyloBayes, MrBayes, IQ-Tree, and RAxML. We study their performance in the phylogenomic framework of $>$800 ultraconserved elements from the bee subfamily Nomiinae (Halictidae). Our taxon sampling focuses on the genus Pseudapis, a distinct lineage with diverse morphological features, but contentious morphology-based taxonomic classifications and no molecular phylogenetic guidance. We approximate topological accuracy of gene trees by assessing their ability to recover two uncontroversial, monophyletic groups, and compare branch lengths of individual trees using the stemminess metric (the relative length of internal branches). We further examine different strategies of removing uninformative loci and the collapsing of weakly supported nodes into polytomies. We then summarize gene trees with ASTRAL and compare resulting species phylogenies, including comparisons to concatenation-based estimates. Gene trees obtained with the reversible jump model search in MrBayes were most concordant on average and all Bayesian methods yielded gene trees with better stemminess values. The only gene tree estimation approach whose ASTRAL summary trees consistently produced the most likely correct topology, however, was IQ-Tree with automated model designation (ModelFinder program). We discuss these findings and provide practical advice on gene tree estimation for summary methods. Lastly, we establish the first phylogeny-informed classification for Pseudapis s. l. and map the distribution of distinct morphological features of the group. [ASTRAL; Bees; concordance; gene tree estimation error; IQ-Tree; MrBayes, Nomiinae; PhyloBayes; RAxML; phylogenomics; stemminess].

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Do Alignment and Trimming Methods Matter for Phylogenomic (UCE) Analyses?

Alignment is a crucial issue in molecular phylogenetics because different alignment methods can potentially yield very different topologies for individual genes. But it is unclear if the choice of alignment methods remains important in phylogenomic analyses, which incorporate data from hundreds or thousands of genes. For example, problematic biases in alignment might be multiplied across many loci, whereas alignment errors in individual genes might become irrelevant. The issue of alignment trimming (i.e., removing poorly aligned regions or missing data from individual genes) is also poorly explored. Here, we test the impact of 12 different combinations of alignment and trimming methods on phylogenomic analyses. We compare these methods using published phylogenomic data from ultraconserved elements (UCEs) from squamate reptiles (lizards and snakes), birds, and tetrapods. We compare the properties of alignments generated by different alignment and trimming methods (e.g., length, informative sites, missing data). We also test whether these data sets can recover well-established clades when analyzed with concatenated (RAxML) and species-tree methods (ASTRAL-III), using the full data ($\sim $5000 loci) and subsampled data sets (10% and 1% of loci). We show that different alignment and trimming methods can significantly impact various aspects of phylogenomic data sets (e.g., length, informative sites). However, these different methods generally had little impact on the recovery and support values for well-established clades, even across very different numbers of loci. Nevertheless, our results suggest several "best practices" for alignment and trimming. Intriguingly, the choice of phylogenetic methods impacted the phylogenetic results most strongly, with concatenated analyses recovering significantly more well-established clades (with stronger support) than the species-tree analyses. [Alignment; concatenated analysis; phylogenomics; sequence length heterogeneity; species-tree analysis; trimming].

Higher temperatures lower rates of physiological and niche evolution

Understanding rates and patterns of change in physiological and climatic-niche variables is of urgent importance as many species are increasingly threatened by rising global temperatures. Here, we broadly test several fundamental hypotheses about physiological and niche evolution for the first time (with appropriate phylogenetic methods), using published data from 2059 vertebrate species. Our main results show that: (i) physiological tolerances to heat evolve more slowly than those to cold, (ii) the hottest climatic-niche temperatures change more slowly than the coldest climatic-niche temperatures, and (iii) physiological tolerances to heat and cold evolve more slowly than the corresponding climatic-niche variables. Physiological tolerances are significantly and positively related to the corresponding climatic-niche variables, but species often occur in climates outside the range of these tolerances. However, mismatches between climate and physiology do not necessarily mean that the climatic-niche data are misleading. Instead, some standard physiological variables used in vertebrates (i.e. critical thermal maxima and minima) may reflect when species are active (daily, seasonally) and their local-scale microhabitats (sun versus shade), rather than their large-scale climatic distributions.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

Analysis of ultraconserved elements supports African origins of narrow-mouthed frogs

Narrow-mouthed frogs (Anura: Microhylidae) are globally distributed and molecular data suggest the rapid evolution of multiple subfamilies shortly after their origin. Despite recent progress, several subfamilial relationships remain unexplored using phylogenomic data. We analysed 1,796 nuclear ultraconserved elements, a total matrix of 400,664 nucleotides, from representatives of most microhylid subfamilies. Summary method species-tree and maximum likelihood analyses unambiguously supported Hoplophryninae as the earliest diverging microhylid and confirm Chaperininae as a junior synonym of Microhylinae. Given the emerging consensus that subfamilies from mainland Africa diverged early, microhylids have likely occupied the continent for more than 66 million years.

Phylogenomic Reconstruction of the Neotropical Poison Frogs (Dendrobatidae) and Their Conservation

The evolutionary history of the Dendrobatidae, the charismatic Neotropical poison frog family, remains in flux, even after a half-century of intensive research. Understanding the evolutionary relationships between dendrobatid genera and the larger-order groups within Dendrobatidae is critical for making accurate assessments of all aspects of their biology and evolution. In this study, we provide the first phylogenomic reconstruction of Dendrobatidae with genome-wide nuclear markers known as ultraconserved elements. We performed sequence capture on 61 samples representing 33 species across 13 of the 16 dendrobatid genera, aiming for a broadly representative taxon sample. We compare topologies generated using maximum likelihood and coalescent methods and estimate divergence times using Bayesian methods. We find most of our dendrobatid tree to be consistent with previously published results based on mitochondrial and low-count nuclear data, with notable exceptions regarding the placement of Hyloxalinae and certain genera within Dendrobatinae. We also characterize how the evolutionary history and geographic distributions of the 285 poison frog species impact their conservation status. We hope that our phylogeny will serve as a backbone for future evolutionary studies and that our characterizations of conservation status inform conservation practices while highlighting taxa in need of further study.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

Taxonomy, systematics and geographic distribution of ground frogs (Alsodidae, Eupsophus): a comprehensive synthesis of the last six decades of research

The genus Eupsophus (ground frogs) inhabits exclusively the temperate forests of southern South America (Chile and Argentina). The current delimitation of the genus was reached in the late 1970s, when only two species were recognized, but since then the number of described species steadily increased, reaching a maximum of 11 by 2012. Subsequent studies that applied explicit species delimitation approaches decreased the number of species to six in 2017 and raised it again to 11 the following year, including an undescribed putative species. Despite these taxonomic changes, the two species groups traditionally recognized, roseus and vertebralis, have been maintained. Another recent contribution to the taxonomy of the genus was the explicit recognition of the extremely high level of external phenotypic variation exhibited by species of the roseus group, which undermines the utility of some diagnostic characters. Here we provide a critical review of the extensive taxonomic and systematic literature on the genus over the last six decades, to examine the evidence behind the recurrent taxonomic changes and advances in its systematics. We also update and complete a 2017 review of geographic information, provide additional qualitative observations of external characters commonly used in the diagnoses of species of the roseus group, and reassess the phylogenetic position of a putative new species from Tolhuaca (Chile), which was not included in the last species delimitation study. The present review shows that: 1) there is no congruence between the patterns of phenotypic and genetic/phylogenetic differentiation among species of both groups; 2) in the roseus group, the intraspecific variation in some external characters is as high as the differences described among species; 3) there is little morphological and bioacoustic differentiation within species groups, and inconsistencies in the chromosomal evidence at the genus level; 4) under the latest taxonomic proposal (2018), species of the roseus group still lack consistent and reliable diagnoses and their distribution limits are poorly defined; and 5) the population from Tolhuaca represents an additional undescribed species under the most recent taxonomic framework. Finally, we discuss the implications of these findings for the taxonomy and biogeography of the genus, pointing out some areas that require further research to understand their patterns and processes of diversification.

Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce. In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open-source UCE probe set for beetles of the suborder Adephaga.

Exploring data processing strategies in NGS target enrichment to disentangle radiations in the tribe Cardueae (Compositae)

Target enrichment is a cost-effective sequencing technique that holds promise for elucidating evolutionary relationships in fast-evolving lineages. However, potential biases and impact of bioinformatic sequence treatments in phylogenetic inference have not been thoroughly explored yet. Here, we investigate this issue with an ultimate goal to shed light into a highly diversified group of Compositae (Asteraceae) constituted by four main genera: Arctium, Cousinia, Saussurea, and Jurinea. Specifically, we compared sequence data extraction methods implemented in two easy-to-use workflows, PHYLUCE and HybPiper, and assessed the impact of two filtering practices intended to reduce phylogenetic noise. In addition, we compared two phylogenetic inference methods: (1) the concatenation approach, in which all loci were concatenated in a supermatrix; and (2) the coalescence approach, in which gene trees were produced independently and then used to construct a species tree under coalescence assumptions. Here we confirm the usefulness of the set of 1061 COS targets (a nuclear conserved orthology loci set developed for the Compositae) across a variety of taxonomic levels. Intergeneric relationships were completely resolved: there are two sister groups, Arctium-Cousinia and Saussurea-Jurinea, which are in agreement with a morphological hypothesis. Intrageneric relationships among species of Arctium, Cousinia, and Saussurea are also well defined. Conversely, conflicting species relationships remain for Jurinea. Methodological choices significantly affected phylogenies in terms of topology, branch length, and support. Across all analyses, the phylogeny obtained using HybPiper and the strictest scheme of removing fast-evolving sites was estimated as the optimal. Regarding methodological choices, we conclude that: (1) trees obtained under the coalescence approach are topologically more congruent between them than those inferred using the concatenation approach; (2) refining treatments only improved support values under the concatenation approach; and (3) branch support values are maximized when fast-evolving sites are removed in the concatenation approach, and when a higher number of loci is analyzed in the coalescence approach.