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.

Metagenome assembly of high-fidelity long reads with hifiasm-meta.

De novo assembly of metagenome samples is a common approach to the study of microbial communities. Current metagenome assemblers developed for short sequence reads or noisy long reads were not optimized for accurate long reads. We thus developed hifiasm-meta, a metagenome assembler that exploits the high accuracy of recent data. Evaluated on seven empirical datasets, hifiasm-meta reconstructed tens to hundreds of complete circular bacterial genomes per dataset, consistently outperforming other metagenome assemblers.

Redefining Possible: Combining Phylogenomic and Supersparse Data in Frogs.

The data available for reconstructing molecular phylogenies have become wildly disparate. Phylogenomic studies can generate data for thousands of genetic markers for dozens of species, but for hundreds of other taxa, data may be available from only a few genes. Can these two types of data be integrated to combine the advantages of both, addressing the relationships of hundreds of species with thousands of genes? Here, we show that this is possible, using data from frogs. We generated a phylogenomic data set for 138 ingroup species and 3,784 nuclear markers (ultraconserved elements [UCEs]), including new UCE data from 70 species. We also assembled a supermatrix data set, including data from 97% of frog genera (441 total), with 1-307 genes per taxon. We then produced a combined phylogenomic-supermatrix data set (a "gigamatrix") containing 441 ingroup taxa and 4,091 markers but with 86% missing data overall. Likelihood analysis of the gigamatrix yielded a generally well-supported tree among families, largely consistent with trees from the phylogenomic data alone. All terminal taxa were placed in the expected families, even though 42.5% of these taxa each had >99.5% missing data and 70.2% had >90% missing data. Our results show that missing data need not be an impediment to successfully combining very large phylogenomic and supermatrix data sets, and they open the door to new studies that simultaneously maximize sampling of genes and taxa.

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.

Evaluation of taxonomic classification and profiling methods for long-read shotgun metagenomic sequencing datasets.

BACKGROUND: Long-read shotgun metagenomic sequencing is gaining in popularity and offers many advantages over short-read sequencing. The higher information content in long reads is useful for a variety of metagenomics analyses, including taxonomic classification and profiling. The development of long-read specific tools for taxonomic classification is accelerating, yet there is a lack of information regarding their relative performance. Here, we perform a critical benchmarking study using 11 methods, including five methods designed specifically for long reads. We applied these tools to several mock community datasets generated using Pacific Biosciences (PacBio) HiFi or Oxford Nanopore Technology sequencing, and evaluated their performance based on read utilization, detection metrics, and relative abundance estimates. RESULTS: Our results show that long-read classifiers generally performed best. Several short-read classification and profiling methods produced many false positives (particularly at lower abundances), required heavy filtering to achieve acceptable precision (at the cost of reduced recall), and produced inaccurate abundance estimates. By contrast, two long-read methods (BugSeq, MEGAN-LR & DIAMOND) and one generalized method (sourmash) displayed high precision and recall without any filtering required. Furthermore, in the PacBio HiFi datasets these methods detected all species down to the 0.1% abundance level with high precision. Some long-read methods, such as MetaMaps and MMseqs2, required moderate filtering to reduce false positives to resemble the precision and recall of the top-performing methods. We found read quality affected performance for methods relying on protein prediction or exact k-mer matching, and these methods performed better with PacBio HiFi datasets. We also found that long-read datasets with a large proportion of shorter reads (< 2 kb length) resulted in lower precision and worse abundance estimates, relative to length-filtered datasets. Finally, for classification methods, we found that the long-read datasets produced significantly better results than short-read datasets, demonstrating clear advantages for long-read metagenomic sequencing. CONCLUSIONS: Our critical assessment of available methods provides best-practice recommendations for current research using long reads and establishes a baseline for future benchmarking studies.

Phylogeography and population genetics of a widespread cold-adapted ant, Prenolepis imparis.

As species arise, evolve and diverge, they are shaped by forces that unfold across short and long timescales and at both local and vast geographical scales. It is rare, however, to be able document this history across broad sweeps of time and space in a single species. Here, we report the results of a continental-scale phylogenomic analysis across the entire range of a widespread species. We analysed sequences of 1402 orthologous ultraconserved element (UCE) loci from 75 individuals to identify population genetic structure and historical demographic patterns across the continent-wide range of a cold-adapted ant, the winter ant, Prenolepis imparis. We recovered five well-supported, genetically isolated clades representing lineages that diverged from 8.2-2.2 million years ago. These include: (i) an early diverging lineage located in Florida, (ii) a lineage that spans the southern United States, (iii) populations that extend across the midwestern and northeastern United States, (iv) populations from the western United States and (v) populations in southwestern Arizona and Mexico. Population genetic analyses revealed little or no gene flow among these lineages, but patterns consistent with more recent gene flow among populations within lineages, and localized structure with migration in the western United States. High support for five major geographical lineages and lack of evidence of contemporary gene flow indicate in situ diversification across the species' range, producing relatively ancient lineages that persisted through subsequent climate change and glaciation during the Quaternary.

Giant Tree Frog diversification in West and Central Africa: Isolation by physical barriers, climate, and reproductive traits.

Secondary sympatry amongst sister lineages is strongly associated with genetic and ecological divergence. This pattern suggests that for closely related species to coexist in secondary sympatry, they must accumulate differences in traits that mediate ecological and/or reproductive isolation. Here, we characterized inter- and intraspecific divergence in three giant tree frog species whose distributions stretch across West and Central Africa. Using genome-wide single-nucleotide polymorphism data, we demonstrated that species-level divergence coincides temporally and geographically with a period of large-scale forest fragmentation during the late Pliocene. Our environmental niche models further supported a dynamic history of climatic suitability and stability, and indicated that all three species occupy distinct environmental niches. We found modest morphological differentiation amongst the species with significant divergence in tympanum diameter and male advertisement call. In addition, we confirmed that two species occur in secondary sympatry in Central Africa but found no evidence of hybridization. These patterns support the hypothesis that cycles of genetic exchange and isolation across West and Central Africa have contributed to globally significant biodiversity. Furthermore, divergence in both ecology and reproductive traits appear to have played important roles in maintaining distinct lineages. At the intraspecific level, we found that climatic refugia, precipitation gradients, marine incursions, and potentially riverine barriers generated phylogeographic structure throughout the Pleistocene and into the Holocene. Further studies examining phenotypic divergence and secondary contact amongst these geographically structured populations may demonstrate how smaller scale and more recent biogeographic barriers contribute to regional diversification.

La sympatrie secondaire parmi les espèces sœurs est fortement associée à la divergence génétique et écologique. Ce modèle suggère que pour que des espèces étroitement liées coexistent en sympatrie secondaire, elles doivent accumuler des différences dans les traits qui contribuent à l'isolement écologique ou reproductif. Ici, nous avons caractérisé la divergence inter- et intra-spécifique chez trois espèces de grenouilles arboricoles géantes dont les distributions s'étendent à travers l'Afrique de l'Ouest et Centrale. Avec des données génétiques, nous avons démontré que la divergence au niveau des espèces coïncide temporellement et géographiquement avec une période de fragmentation forestière à la fin du Pliocène. Nos modèles de niches environnementales ont soutenu une histoire dynamique de stabilité climatique, et ont indiqué que les trois espèces occupent des niches environnementales distinctes. Nous avons trouvé une différenciation morphologique modeste parmi les trois espèces mais une divergence significative dans le diamètre du tympan et les cris des mâles. De plus, nous avons confirmé que deux espèces sont présentes en sympatrie secondaire en Afrique Centrale mais n'avons trouvé aucune preuve d'hybridation. Ces résultats soutiennent l'hypothèse que les cycles d'échange génétique et d'isolement à travers l'Afrique de l'Ouest et Centrale ont contribué à une profonde concentration de biodiversité dans la région. De plus, la divergence des traits écologiques et reproducteurs semble avoir joué un rôle important dans le maintien de lignées distinctes. Au niveau intra-spécifique, nous avons constaté que les refuges climatiques, les gradients de précipitation, les incursions marines et potentiellement les barrières fluviales ont généré une structure phylogéographique pendant le Pléistocène et jusqu'à l'Holocène. Des études examinant la divergence phénotypique et le contact secondaire entre ces populations géographiquement structurées pourraient démontrer comment des barrières biogéographiques à échelle plus petite et plus récentes contribuent à la diversification régionale.

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

Phylogenomics of Monitor Lizards and the Role of Competition in Dictating Body Size Disparity.

Organismal interactions drive the accumulation of diversity by influencing species ranges, morphology, and behavior. Interactions vary from agonistic to cooperative and should result in predictable patterns in trait and range evolution. However, despite a conceptual understanding of these processes, they have been difficult to model, particularly on macroevolutionary timescales and across broad geographic spaces. Here, we investigate the influence of biotic interactions on trait evolution and community assembly in monitor lizards (Varanus). Monitors are an iconic radiation with a cosmopolitan distribution and the greatest size disparity of any living terrestrial vertebrate genus. Between the colossal Komodo dragon Varanus komodoensis and the smallest Australian dwarf goannas, Varanus length and mass vary by multiple orders of magnitude. To test the hypothesis that size variation in this genus was driven by character displacement, we extended existing phylogenetic comparative methods which consider lineage interactions to account for dynamic biogeographic history and apply these methods to Australian monitors and marsupial predators. Incorporating both exon-capture molecular and morphological data sets we use a combined evidence approach to estimate the relationships among living and extinct varaniform lizards. Our results suggest that communities of Australian Varanus show high functional diversity as a result of continent-wide interspecific competition among monitors but not with faunivorous marsupials. We demonstrate that patterns of trait evolution resulting from character displacement on continental scales are recoverable from comparative data and highlight that these macroevolutionary patterns may develop in parallel across widely distributed sympatric groups.[Character displacement; comparative methods; phylogenetics; trait evolution; Varanus.].

Cryptic diversity of a widespread global pathogen reveals expanded threats to amphibian conservation.

Biodiversity loss is one major outcome of human-mediated ecosystem disturbance. One way that humans have triggered wildlife declines is by transporting disease-causing agents to remote areas of the world. Amphibians have been hit particularly hard by disease due in part to a globally distributed pathogenic chytrid fungus (Batrachochytrium dendrobatidis [Bd]). Prior research has revealed important insights into the biology and distribution of Bd; however, there are still many outstanding questions in this system. Although we know that there are multiple divergent lineages of Bd that differ in pathogenicity, we know little about how these lineages are distributed around the world and where lineages may be coming into contact. Here, we implement a custom genotyping method for a global set of Bd samples. This method is optimized to amplify and sequence degraded DNA from noninvasive skin swab samples. We describe a divergent lineage of Bd, which we call BdASIA3, that appears to be widespread in Southeast Asia. This lineage co-occurs with the global panzootic lineage (BdGPL) in multiple localities. Additionally, we shed light on the global distribution of BdGPL and highlight the expanded range of another lineage, BdCAPE. Finally, we argue that more monitoring needs to take place where Bd lineages are coming into contact and where we know little about Bd lineage diversity. Monitoring need not use expensive or difficult field techniques but can use archived swab samples to further explore the history-and predict the future impacts-of this devastating pathogen.

Finding complexity in complexes: Assessing the causes of mitonuclear discordance in a problematic species complex of Mesoamerican toads.

Mitonuclear discordance is a frequently encountered pattern in phylogeographic studies and occurs when mitochondrial and nuclear DNA display conflicting signals. Discordance among these genetic markers can be caused by several factors including confounded taxonomies, gene flow, and incomplete lineage sorting. In this study, we present a strong case of mitonuclear discordance in a species complex of toads (Bufonidae: Incilius coccifer complex) found in the Chortís Block of Central America. To determine the cause of mitonuclear discordance in this complex, we used spatially explicit genetic data to test species limits and relationships, characterize demographic history, and quantify gene flow. We found extensive mitonuclear discordance among the three recognized species within this group, especially in populations within the Chortís Highlands of Honduras. Our data reveal nuclear introgression within the Chortís Highlands populations that was most probably driven by cyclical range expansions due to climatic fluctuations. Though we determined introgression occurred within the nuclear genome, our data suggest that it is not the key factor in driving mitonuclear discordance in the entire species complex. Rather, due to a lack of discernible geographic pattern between mitochondrial and nuclear DNA, as well as a relatively recent divergence time of this complex, we concluded that mitonuclear discordance has been caused by incomplete lineage sorting. Our study provides a framework to test sources of mitonuclear discordance and highlights the importance of using multiple marker types to test species boundaries in cryptic species.

Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians.

Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.

Macroevolutionary Patterns of Sexual Size Dimorphism Among African Tree Frogs (Family: Hyperoliidae).

Sexual size dimorphism (SSD) is shaped by multiple selective forces that drive the evolution of sex-specific body size, resulting in male or female-biased SSD. Stronger selection on one sex can result in an allometric body size scaling relationship consistent with Rensch's rule or its converse. Anurans (frogs and toads) generally display female-biased SSD, but there is variation across clades and the mechanisms driving the evolution of SSD remain poorly understood. We investigated these topics in a diverse family of African treefrogs (Hyperoliidae). Hyperoliids display traits considered rare among amphibians, including sexual dichromatism and protogynous sex change. Using phylogenetic comparative methods, we tested if adult ecology, sexual dichromatism, and sex change were predictors of body size or SSD. We also tested whether hyperoliids displayed allometric interspecific body size scaling relationships. We found a majority of hyperoliid taxa display female-biased SSD, but that adult ecology and sexual dichromatism are poor predictors of sex-specific body size and SSD. Regardless of the groupings analyzed (partitioned by clades or traits), we found support for isometric body size scaling. However, we found that sex change is a significant predictor of SSD variation. Species in the Hyperolius viridiflavus complex, which putatively display this trait, show a significant reduction in SSD and are frequently sexually monomorphic in size. Although protogynous sex change needs to be validated for several of these species, we tentatively propose this trait is a novel mechanism influencing anuran body size evolution. Beyond this association, additional factors that shape the evolution of anuran body size and SSD remain elusive.

Vanishing refuge? Testing the forest refuge hypothesis in coastal East Africa using genome-wide sequence data for seven amphibians.

High-throughput sequencing data have greatly improved our ability to understand the processes that contribute to current biodiversity patterns. The "vanishing refuge" diversification model is speculated for the coastal forests of eastern Africa, whereby some taxa have persisted and diversified between forest refugia, while others have switched to becoming generalists also present in non-forest habitats. Complex arrangements of geographical barriers (hydrology and topography) and ecological gradients between forest and non-forest habitats may have further influenced the region's biodiversity, but elucidation of general diversification processes has been limited by lack of suitable data. Here, we explicitly test alternative diversification modes in the coastal forests using genome-wide single nucleotide polymorphisms, mtDNA, spatial and environmental data for three forest (Arthroleptis xenodactyloides, Leptopelis flavomaculatus and Afrixalus sylvaticus) and four generalist (Afrixalus fornasini, A. delicatus, Leptopelis concolor and Leptopelis argenteus) amphibians. Multiple analyses provide insight about divergence times, spatial population structure, dispersal barriers, environmental stability and demographic history. We reveal highly congruent intra-specific diversity and population structure across taxa, with most divergences occurring during the late Pliocene and Pleistocene. Although stability models support the existence of some forest refugia, dispersal barriers and demographic models point towards idiosyncratic diversification modes across taxa. We identify a consistent role for riverine barriers in the diversification of generalist taxa, but mechanisms of diversification are more complex for forest taxa and potentially include topographical barriers, forest refugia and ecological gradients. Our work demonstrates the complexity of diversification processes in this region, which vary between forest and generalist taxa, but also for ecologically similar species with shared population boundaries.

Diversity and biogeography of frogs in the genus Amnirana (Anura: Ranidae) across sub-Saharan Africa.

Frogs in the genus Amnirana (family Ranidae) are widely distributed across sub-Saharan Africa and present a model system for exploring the relationship between diversification and geography across the continent. Using multiple loci from the mitochondrial (16S) and nuclear genomes (DISP2, FICD, KIAA2013, REV3L), we generated a strongly supported species-level phylogeny that provides insights into the continental biogeography of African species of Amnirana, which form a monophyletic group within the genus. Species delimitation analyses suggest that there may be as many as seven additional species of Amnirana in Africa. The biogeographic history of Amnirana is marked by several dispersal and vicariance events, including dispersal from the Lower Guinean Forest into the Congo Basin. In addition, phylogeographic patterns within two widespread species, A. albolabris and A. galamensis, reveal undescribed cryptic diversity. Populations assigned to A. albolabris in western Africa are more closely related to A. fonensis and require recognition as a distinct species. Our analyses reveal that the Lower and Upper Guinean Forest regions served as important centers of interspecific and intraspecific diversifications for Amnirana.

Evaluating mechanisms of diversification in a Guineo-Congolian tropical forest frog using demographic model selection.

The accumulation of biodiversity in tropical forests can occur through multiple allopatric and parapatric models of diversification, including forest refugia, riverine barriers and ecological gradients. Considerable debate surrounds the major diversification process, particularly in the West African Lower Guinea forests, which contain a complex geographic arrangement of topographic features and historical refugia. We used genomic data to investigate alternative mechanisms of diversification in the Gaboon forest frog, Scotobleps gabonicus, by first identifying population structure and then performing demographic model selection and spatially explicit analyses. We found that a majority of population divergences are best explained by allopatric models consistent with the forest refugia hypothesis and involve divergence in isolation with subsequent expansion and gene flow. These population divergences occurred simultaneously and conform to predictions based on climatically stable regions inferred through ecological niche modelling. Although forest refugia played a prominent role in the intraspecific diversification of S. gabonicus, we also find evidence for potential interactions between landscape features and historical refugia, including major rivers and elevational barriers such as the Cameroonian Volcanic Line. We outline the advantages of using genomewide variation in a model-testing framework to distinguish between alternative allopatric hypotheses, and the pitfalls of limited geographic and molecular sampling. Although phylogeographic patterns are often species-specific and related to life-history traits, additional comparative studies incorporating genomic data are necessary for separating shared historical processes from idiosyncratic responses to environmental, climatic and geological influences on diversification.

Idiosyncratic responses to climate-driven forest fragmentation and marine incursions in reed frogs from Central Africa and the Gulf of Guinea Islands.

Organismal traits interact with environmental variation to mediate how species respond to shared landscapes. Thus, differences in traits related to dispersal ability or physiological tolerance may result in phylogeographic discordance among co-distributed taxa, even when they are responding to common barriers. We quantified climatic suitability and stability, and phylogeographic divergence within three reed frog species complexes across the Guineo-Congolian forests and Gulf of Guinea archipelago of Central Africa to investigate how they responded to a shared climatic and geological history. Our species-specific estimates of climatic suitability through time are consistent with temporal and spatial heterogeneity in diversification among the species complexes, indicating that differences in ecological breadth may partly explain these idiosyncratic patterns. Likewise, we demonstrated that fluctuating sea levels periodically exposed a land bridge connecting Bioko Island with the mainland Guineo-Congolian forest and that habitats across the exposed land bridge likely enabled dispersal in some species, but not in others. We did not find evidence that rivers are biogeographic barriers across any of the species complexes. Despite marked differences in the geographic extent of stable climates and temporal estimates of divergence among the species complexes, we recovered a shared pattern of intermittent climatic suitability with recent population connectivity and demographic expansion across the Congo Basin. This pattern supports the hypothesis that genetic exchange across the Congo Basin during humid periods, followed by vicariance during arid periods, has shaped regional diversity. Finally, we identified many distinct lineages among our focal taxa, some of which may reflect incipient or unrecognized species.

Impact of species delimitation and sampling on niche models and phylogeographical inference: A case study of the East African reed frog Hyperolius substriatus Ahl, 1931.

Ecological niche models (ENMs) have been used in a wide range of ecological and evolutionary studies. In biogeographic studies these models have, among other things, helped in the discovery of new allopatric populations, and even new species. However, small sample sizes and questionable taxonomic delimitation can challenge models, often decreasing their accuracy. Herein we examine the sensitivity of ENMs to the addition of new, geographically isolated populations, and the impact of applying different taxonomic delimitations. The East African reed frog Hyperolius substriatus Ahl, 1931 was selected as a case study because it has been the subject of previous ENM predictions. Our results suggest that addition of new data and reanalysis of species lineages of H. substriatus improved our understanding of the evolutionary history of this group of frogs. ENMs provided robust predictions, even when some populations were deliberately excluded from the models. Splitting the lineages based on genetic relationships and analysing the ENMs separately provided insights about the biogeographical processes that led to the current distribution of H. substriatus.

Molecular data from contemporary and historical collections reveal a complex story of cryptic diversification in the Varanus (Polydaedalus) niloticus Species Group.

Previous studies of color pattern, tongue pigmentation, and scale counts have been used to distinguish two species of semiaquatic varanids in Africa, but these findings have yet to be tested with molecular data. The Varanus (Polydaedalus) niloticus Species Group is comprised of the Nile monitor (V. niloticus) and the Ornate monitor (V. ornatus). Due to the high rate of exploitation of both species for bushmeat, the leather industry, and the pet trade, a clear understanding of the taxonomy and genetic partitioning is necessary for effective management. Here we utilize a multilocus approach, consisting of mitochondrial and nuclear markers, totaling 4251 bp, as well as microsatellite loci to assess the taxonomic validity and intraspecific evolutionary patterns within the V. niloticus Species Group. By incorporating historical specimens from museum collections as well as contemporary samples, we obtained range-wide coverage for both species across Africa. Concordant results from various approaches all suggest that V. ornatus does not represent a distinct monophyletic group. Our analyses recovered three genetic clades within V. niloticus, representing western, northern, and southern lineages. The western clade was found to diverge first, around 7.7 mya (95% HPD: 4.6-11.0 mya) and exhibits 8.4% and 8.7% uncorrected sequence divergence between the northern and southern V. niloticus clades, respectively. This geographically separate lineage corresponds to previous descriptions of Tupinambis stellatusDaudin (1802). These findings not only call for taxonomic revision of this species group, but also shed light on the biogeographic history of Africa as well as aid in the management planning of varanids and other co-distributed African species.

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.