Multiple Routes to Color Convergence in a Radiation of Neotropical Poison Frogs.

Convergent evolution is defined as the independent evolution of similar phenotypes in different lineages. Its existence underscores the importance of external selection pressures in evolutionary history, revealing how functionally similar adaptations can evolve in response to persistent ecological challenges through a diversity of evolutionary routes. However, many examples of convergence, particularly among closely related species, involve parallel changes in the same genes or developmental pathways, raising the possibility that homology at deeper mechanistic levels is an important facilitator of phenotypic convergence. Using the genus Ranitomeya, a young, color-diverse radiation of Neotropical poison frogs, we set out to 1) provide a phylogenetic framework for this group, 2) leverage this framework to determine if color phenotypes are convergent, and 3) to characterize the underlying coloration mechanisms to test whether color convergence occurred through the same or different physical mechanisms. We generated a phylogeny for Ranitomeya using ultraconserved elements and investigated the physical mechanisms underlying bright coloration, focusing on skin pigments. Using phylogenetic comparative methods, we identified several instances of color convergence, involving several gains and losses of carotenoid and pterin pigments. We also found a compelling example of nonparallel convergence, where, in one lineage, red coloration evolved through the red pterin pigment drosopterin, and in another lineage through red ketocarotenoids. Additionally, in another lineage, "reddish" coloration evolved predominantly through structural color mechanisms. Our study demonstrates that, even within a radiation of closely related species, convergent evolution can occur through both parallel and nonparallel mechanisms, challenging the assumption that similar phenotypes among close relatives evolve through the same mechanisms.

Antifungal-Loaded Calcium Sulfate Beads as a Potential Therapeutic in Combating Candida auris.

Candida auris provides a substantial global nosocomial threat clinically. With the recent emergence that the organism can readily colonize skin niches, it will likely continue to pose a risk in health care units, particularly to patients undergoing surgery. The purpose of this study was to investigate the efficacy of antifungal-loaded calcium sulfate (CS) beads in combatting C. auris infection. We demonstrate that the CS-packed beads have the potential to interfere with planktonic and sessile C. auris.

Phylogenomic analysis of evolutionary relationships in Ranitomeya poison frogs (Family Dendrobatidae) using ultraconserved elements.

The use of genome-scale data in phylogenetics has enabled recent strides in determining the relationships between taxa that are taxonomically problematic because of extensive morphological variation. Here, we employ a phylogenomic approach to infer evolutionary relationships within Ranitomeya (Anura: Dendrobatidae), an Amazonian lineage of poison frogs consisting of 16 species with remarkable diversity in color pattern, range size, and parental care behavior. We infer phylogenies with all described species of Ranitomeya from ultraconserved nuclear genomic elements (UCEs) and also estimate divergence times. Our results differ from previous analyses regarding interspecific relationships. Notably, we find that R. toraro and R. defleri are not sister species but rather distantly related, contrary to previous analyses based on smaller genetic datasets. We recover R. uakarii as paraphyletic, designate certain populations formerly assigned to R. fantastica from Peru as R. summersi, and transfer the French Guianan and eastern Brazilian R. amazonica populations to R. variabilis. By clarifying both inter- and intraspecific relationships within Ranitomeya, our study paves the way for future tests of hypotheses on color pattern evolution and historical biogeography.

A New Method for Integrating Ecological Niche Modeling with Phylogenetics to Estimate Ancestral Distributions.

Ancestral range estimation and projection of niche models into the past have both become common in evolutionary studies where the ancient distributions of organisms are in question. However, these methods are hampered by complementary hurdles: discrete characterization of areas in ancestral range estimation can be overly coarse, especially at shallow timescales, and niche model projection neglects evolution. Phylogenetic niche modeling accounts for both of these issues by incorporating knowledge of evolutionary relationships into a characterization of environmental tolerances. We present a new method for phylogenetic niche modeling, implemented in R. Given past and present climate data, taxon occurrence data, and a time-calibrated phylogeny, our method constructs niche models for each extant taxon, uses ancestral character estimation to reconstruct ancestral niche models, and projects these models into paleoclimate data to provide a historical estimate of the geographic range of a lineage. Models either at nodes or along branches of the phylogeny can be estimated. We demonstrate our method on a small group of dendrobatid frogs and show that it can make inferences given species with restricted ranges and little occurrence data. We also use simulations to show that our method can reliably reconstruct the niche of a known ancestor in both geographic and environmental space. Our method brings together fields as disparate as ecological niche modeling, phylogenetics, and ancestral range estimation in a user-friendly package. [Ancestral range estimation; ancestral state reconstruction; biogeography; Dendrobatidae; ecological niche modeling; paleoclimate; phylogeography; species distribution modeling.].

Assessing the Bioactive Profile of Antifungal-Loaded Calcium Sulfate against Fungal Biofilms.

Calcium sulfate (CS) has been used clinically as a bone- or void-filling biomaterial, and its resorptive properties have provided the prospect for its use as a release mechanism for local antibiotics to control biofilms. Here, we aimed to test CS beads loaded with three antifungal drugs against planktonic and sessile fungal species to assess whether these antifungal beads could be harnessed to provide consistent release of antifungals at biofilm-inhibitory doses. A panel of different fungal species (n = 15) were selected for planktonic broth microdilution testing with fluconazole (FLZ), amphotericin B (AMB), and caspofungin (CSP). After establishing planktonic inhibition, antifungal CS beads were introduced to fungal biofilms (n = 5) to assess biofilm formation and cell viability through a combination of standard quantitative and qualitative biofilm assays. Inoculation of a hydrogel substrate, packed with antifungal CS beads, was also used to assess diffusion through a semidry material, to mimic active infection in vivo In general, antifungals released from loaded CS beads were all effective at inhibiting the pathogenic fungi over 7 days within standard MIC ranges for these fungi. We observed a significant reduction of pregrown fungal biofilms across key fungal pathogens following treatment, with visually observable changes in cell morphology and biofilm coverage provided by scanning electron microscopy. Assessment of biofilm inhibition also revealed reductions in total and viable cells across all organisms tested. These data show that antifungal-loaded CS beads produce a sustained antimicrobial effect that inhibits and kills clinically relevant fungal species in vitro as planktonic and biofilm cells.

Glycation of Host Proteins Increases Pathogenic Potential of Porphyromonas gingivalis.

The non-enzymatic addition of glucose (glycation) to circulatory and tissue proteins is a ubiquitous pathophysiological consequence of hyperglycemia in diabetes. Given the high incidence of periodontitis and diabetes and the emerging link between these conditions, it is of crucial importance to define the basic virulence mechanisms employed by periodontopathogens such as Porphyromonas gingivalis in mediating the disease process. The aim of this study was to determine whether glycated proteins are more easily utilized by P. gingivalis to stimulate growth and promote the pathogenic potential of this bacterium. We analyzed the properties of three commonly encountered proteins in the periodontal environment that are known to become glycated and that may serve as either protein substrates or easily accessible heme sources. In vitro glycated proteins were characterized using colorimetric assays, mass spectrometry, far- and near-UV circular dichroism and UV-visible spectroscopic analyses and SDS-PAGE. The interaction of glycated hemoglobin, serum albumin and type one collagen with P. gingivalis cells or HmuY protein was examined using spectroscopic methods, SDS-PAGE and co-culturing P. gingivalis with human keratinocytes. We found that glycation increases the ability of P. gingivalis to acquire heme from hemoglobin, mostly due to heme sequestration by the HmuY hemophore-like protein. We also found an increase in biofilm formation on glycated collagen-coated abiotic surfaces. We conclude that glycation might promote the virulence of P. gingivalis by making heme more available from hemoglobin and facilitating bacterial biofilm formation, thus increasing P. gingivalis pathogenic potential in vivo.

Chitosan Ameliorates Candida auris Virulence in a Galleria mellonella Infection Model.

Candida auris has emerged as a multidrug-resistant nosocomial pathogen over the last decade. Outbreaks of the organism in health care facilities have resulted in life-threatening invasive candidiasis in over 40 countries worldwide. Resistance by C. auris to conventional antifungal drugs such as fluconazole and amphotericin B means that alternative therapeutics must be explored. As such, this study served to investigate the efficacy of a naturally derived polysaccharide called chitosan against aggregative (Agg) and nonaggregative (non-Agg) isolates of C. aurisin vitro and in vivo. In vitro results indicated that chitosan was effective against planktonic and sessile forms of Agg and non-Agg C. auris In a Galleria mellonella model to assess C. auris virulence, chitosan treatment was shown to ameliorate killing effects of both C. auris phenotypes (NCPF 8973 and NCPF 8978, respectively) in vivo Specifically, chitosan reduced the fungal load and increased survival rates of infected Galleria, while treatment alone was nontoxic to the larvae. Finally, chitosan treatment appeared to induce a stress-like gene expression response in NCPF 8973 in the larvae likely arising from a protective response by the organism to resist antifungal activity of the compound. Taken together, results from this study demonstrate that naturally derived compounds such as chitosan may be useful alternatives to conventional antifungals against C. auris.

Phylogenetic relationships and systematics of the Amazonian poison frog genus Ameerega using ultraconserved genomic elements.

The Amazonian poison frog genus Ameerega is one of the largest yet most understudied of the brightly colored genera in the anuran family Dendrobatidae, with 30 described species ranging throughout tropical South America. Phylogenetic analyses of Ameerega are highly discordant, lacking consistency due to variation in data types and methods, and often with limited coverage of species diversity in the genus. Here, we present a comprehensive phylogenomic reconstruction of Ameerega, utilizing state-of-the-art sequence capture techniques and phylogenetic methods. We sequenced thousands of ultraconserved elements from over 100 tissue samples, representing almost every described Ameerega species, as well as undescribed cryptic diversity. We generated topologies using maximum likelihood and coalescent methods and compared the use of maximum likelihood and Bayesian methods for estimating divergence times. Our phylogenetic inference diverged strongly from those of previous studies, and we recommend steps to bring Ameerega taxonomy in line with the new phylogeny. We place several species in a phylogeny for the first time, as well as provide evidence for six potential candidate species. We estimate that Ameerega experienced a rapid radiation approximately 7-11 million years ago and that the ancestor of all Ameerega was likely an aposematic, montane species. This study underscores the utility of phylogenomic data in improving our understanding of the phylogeny of understudied clades and making novel inferences about their evolution.

Speciation with introgression: Phylogeography and systematics of the Ameerega petersi group (Dendrobatidae).

The Tropical Andes contains exceptionally high diversity, much of it arising within the Quaternary period. The complex geology of the Andes and paleoclimate fluctuations within the Quaternary suggest complex speciation scenarios. This, in turn, has contributed to idiosyncratic speciation modes among shallowly diverged Amazonian taxa. Many relationships among these taxa remain poorly resolved. Here we use a sequence capture approach, ultraconserved elements (UCEs), to address the phylogenetic relationships among three recently diverged Peruvian Ameerega poison frog species (A. cainarachi, A. petersi, and A. smaragdina; family Dendrobatidae) and explore a possible mode of speciation in this group. We assess concordance among concatenated phylogenetic tree inference, gene-tree based species tree inference, SNP-based species tree inference, and Bayes factor lineage delimitation to resolve species boundaries. We complement these analyses with assessments of call divergence to address the presence of a prezygotic reproductive barrier. Additionally, we further explore the phylogeographic history of these species of Ameerega with demographic inference, considering evidence for admixture and population expansions. Our results support the synonymy of A. smaragdina as a junior synonym of A. petersi and we find that speciation in this group is characterized by admixture and signatures of a population bottleneck followed by expansion. We invoke the disturbance-vicariance hypothesis to explain the observed patterns and call for more, detailed investigations of in-situ speciation in the Tropical Andes.

An examination of climate-driven flowering-time shifts at large spatial scales over 153 years in a common weedy annual.

PREMISE: Understanding species' responses to climate change is a critical challenge facing biologists today. Though many species are widespread, few studies of climate-driven shifts in flowering time have examined large continuous spatial scales for individual species. And even fewer studies have examined these shifts at time scales greater than a few decades. METHODS: We used digitized herbarium specimens and PRISM climate data to produce the spatially and temporally broadest-scale study of flowering time in a single species to date, spanning the contiguous United States and 153 years (1863-2016) for a widespread weedy annual, Triodanis perfoliata (Campanulaceae). We examined factors driving phenological shifts as well as the roles of geographic and temporal scale in understanding these trends. RESULTS: Year was a significant factor in both geospatial and climatic analyses, revealing that flowering time has advanced by ~9 days over the past ~150 years. We found that temperature as well as vapor pressure deficit, an understudied climatic parameter associated with evapotranspiration and water stress, were strongly associated with peak flowering. We also examined how sampling at different spatiotemporal scales influences the power to detect flowering-time shifts, finding that relatively large spatial and temporal scales are ideal for detecting flowering-time shifts in this widespread species. CONCLUSIONS: Our results emphasize the importance of understanding the interplay of geospatial factors at different scales to examine how species respond to climate change.

Inferring responses to climate dynamics from historical demography in neotropical forest lizards.

We apply a comparative framework to test for concerted demographic changes in response to climate shifts in the neotropical lowland forests, learning from the past to inform projections of the future. Using reduced genomic (SNP) data from three lizard species codistributed in Amazonia and the Atlantic Forest (Anolis punctatus, Anolis ortonii, and Polychrus marmoratus), we first reconstruct former population history and test for assemblage-level responses to cycles of moisture transport recently implicated in changes of forest distribution during the Late Quaternary. We find support for population shifts within the time frame of inferred precipitation fluctuations (the last 250,000 y) but detect idiosyncratic responses across species and uniformity of within-species responses across forest regions. These results are incongruent with expectations of concerted population expansion in response to increased rainfall and fail to detect out-of-phase demographic syndromes (expansions vs. contractions) across forest regions. Using reduced genomic data to infer species-specific demographical parameters, we then model the plausible spatial distribution of genetic diversity in the Atlantic Forest into future climates (2080) under a medium carbon emission trajectory. The models forecast very distinct trajectories for the lizard species, reflecting unique estimated population densities and dispersal abilities. Ecological and demographic constraints seemingly lead to distinct and asynchronous responses to climatic regimes in the tropics, even among similarly distributed taxa. Incorporating such constraints is key to improve modeling of the distribution of biodiversity in the past and future.

Predicting the genetic consequences of future climate change: The power of coupling spatial demography, the coalescent, and historical landscape changes.

PREMISE OF THE STUDY: Climate change is a widely accepted threat to biodiversity. Species distribution models (SDMs) are used to forecast whether and how species distributions may track these changes. Yet, SDMs generally fail to account for genetic and demographic processes, limiting population-level inferences. We still do not understand how predicted environmental shifts will impact the spatial distribution of genetic diversity within taxa. METHODS: We propose a novel method that predicts spatially explicit genetic and demographic landscapes of populations under future climatic conditions. We use carefully parameterized SDMs as estimates of the spatial distribution of suitable habitats and landscape dispersal permeability under present-day, past, and future conditions. We use empirical genetic data and approximate Bayesian computation to estimate unknown demographic parameters. Finally, we employ these parameters to simulate realistic and complex models of responses to future environmental shifts. We contrast parameterized models under current and future landscapes to quantify the expected magnitude of change. KEY RESULTS: We implement this framework on neutral genetic data available from Penstemon deustus. Our results predict that future climate change will result in geographically widespread declines in genetic diversity in this species. The extent of reduction will heavily depend on the continuity of population networks and deme sizes. CONCLUSIONS: To our knowledge, this is the first study to provide spatially explicit predictions of within-species genetic diversity using climatic, demographic, and genetic data. Our approach accounts for climatic, geographic, and biological complexity. This framework is promising for understanding evolutionary consequences of climate change, and guiding conservation planning.

Phylogeography of the arid-adapted Malagasy bullfrog, Laliostoma labrosum, influenced by past connectivity and habitat stability.

The rainforest biome of eastern Madagascar is renowned for its extraordinary biodiversity and restricted distribution ranges of many species, whereas the arid western region of the island is relatively species poor. We provide insight into the biogeography of western Madagascar by analyzing a multilocus phylogeographic dataset assembled for an amphibian, the widespread Malagasy bullfrog, Laliostoma labrosum. We find no cryptic species in L. labrosum (maximum 1.1% pairwise genetic distance between individuals in the 16S rRNA gene) attributable to considerable gene flow at the regional level as shown by genetic admixture in both mtDNA and three nuclear loci, especially in central Madagascar. Low breeding site fidelity, viewed as an adaptation to the unreliability of standing pools of freshwater in dry and seasonal environments, and a ubiquitous distribution within its range may underlie overall low genetic differentiation. Moreover, reductions in population size associated with periods of high aridity in western Madagascar may have purged DNA variation in this species. The mtDNA gene tree revealed seven major phylogroups within this species, five of which show mostly non-overlapping distributions. The nested positions of the northern and central mtDNA phylogroups imply a southwestern origin for all extant mtDNA lineages in L. labrosum. The current phylogeography of this species and paleo-distributions of major mtDNA lineages suggest five potential refugia in northern, western and southwestern Madagascar, likely the result of Pleistocene range fragmentation during drier and cooler climates. Lineage sorting in mtDNA and nuclear loci highlighted a main phylogeographic break between populations north and south of the Sambirano region, suggesting a role of the coastal Sambirano rainforest as a barrier to gene flow. Paleo-species distribution models and dispersal networks suggest that the persistence of some refugial populations was mainly determined by high population connectivity through space and time.

Generation of Multispecies Oral Bacteria Biofilm Models.

It is well-recognized that oral biofilms that occur in health and disease have a polymicrobial composition, though these are poorly reflected in the literature, with many studies focussing on simple mono-species biofilm model systems. The utility of polymicrobial biofilm model systems is that they more accurately reflect the oral cavity and allow researchers to ask relevant questions in basic science studies, pharmaceutical screening, and investigating inflammatory interactions. Here we describe the detailed methodology of how to sequentially construct and maintain polymicrobial biofilm models pertinent to caries, periodontal disease, and denture stomatitis.

Investigating the role of Candida albicans as a universal substrate for oral bacteria using a transcriptomic approach: implications for interkingdom biofilm control?

Candida albicans is frequently identified as a colonizer of the oral cavity in health and has recently been termed a "keystone" commensal due to its role on the bacterial communities. However, the role that C. albicans plays in such interactions is not fully understood. Therefore, this study aimed to identify the relationship between C. albicans and bacteria associated with oral symbiosis and dysbiosis. To do this, we evaluated the ability of C. albicans to support the growth of the aerobic commensal Streptococcus gordonii and the anaerobic pathogens Fusobacterium nucleatum and Porphyromonas gingivalis in the biofilm environment. RNA-Sequencing with the Illumina platform was then utilized to identify C. albicans gene expression and functional pathways involved during such interactions in dual-species and a 4-species biofilm model. Results indicated that C. albicans was capable of supporting growth of all three bacteria, with a significant increase in colony counts of each bacteria in the dual-species biofilm (p < 0.05). We identified specific functional enrichment of pathways in our 4-species community as well as transcriptional profiles unique to the F. nucleatum and S. gordonii dual-species biofilms, indicating a species-specific effect on C. albicans. Candida-related hemin acquisition and heat shock protein mediated processes were unique to the organism following co-culture with anaerobic and aerobic bacteria, respectively, suggestive that such pathways may be feasible options for therapeutic targeting to interfere with these fungal-bacterial interactions. Targeted antifungal therapy may be considered as an option for biofilm destabilization and treatment of complex communities. Moving forward, we propose that further studies must continue to investigate the role of this fungal organism in the context of the interkingdom nature of oral diseases.

Staphylococcus aureus strains exhibit heterogenous tolerance to direct cold atmospheric plasma therapy.

The global clinical and socioeconomic impact of chronic wounds is substantial. The main difficulty that clinicians face during the treatment of chronic wounds is the risk of infection at the wound site. Infected wounds arise from an accumulation of microbial aggregates in the wound bed, leading to the formation of polymicrobial biofilms that can be largely resistant to antibiotic therapy. Therefore, it is essential for studies to identify novel therapeutics to alleviate biofilm infections. One innovative technique is the use of cold atmospheric plasma (CAP) which has been shown to possess promising antimicrobial and immunomodulatory properties. Here, different clinically relevant biofilm models will be treated with cold atmospheric plasma to assess its efficacy and killing effects. Biofilm viability was assessed using live dead qPCR, and morphological changes associated with CAP evaluated using scanning electron microscopy (SEM). Results indicated that CAP was effective against Candida albicans and Pseudomonas aeruginosa, both as mono-species biofilms and when grown in a triadic model system. CAP also significantly reduced viability in the nosocomial pathogen, Candida auris. Staphylococcus aureus Newman exhibited a level of tolerance to CAP therapy, both when grown alone or in the triadic model when grown alongside C. albicans and P. aeruginosa. However, this degree of tolerance exhibited by S. aureus was strain dependent. At a microscopic level, biofilm treatment led to subtle changes in morphology in the susceptible biofilms, with evidence of cellular deflation and shrinkage. Taken together, these results indicate a promising application of direct CAP therapy in combatting wound and skin-related biofilm infections, although biofilm composition may affect the treatment efficacy.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar.

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Spatially explicit models of dynamic histories: examination of the genetic consequences of Pleistocene glaciation and recent climate change on the American Pika.

A central goal of phylogeography is to identify and characterize the processes underlying divergence. One of the biggest impediments currently faced is how to capture the spatiotemporal dynamic under which a species evolved. Here, we described an approach that couples species distribution models (SDMs), demographic and genetic models in a spatiotemporally explicit manner. Analyses of American Pika (Ochotona princeps) from the sky islands of the central Rocky Mountains of North America are used to provide insights into key questions about integrative approaches in landscape genetics, population genetics and phylogeography. This includes (i) general issues surrounding the conversion of time-specific SDMs into simple continuous, dynamic landscapes from past to current, (ii) the utility of SDMs to inform demographic models with deme-specific carrying capacities and migration potentials as well as (iii) the contribution of the temporal dynamic of colonization history in shaping genetic patterns of contemporary populations. Our results support that the inclusion of a spatiotemporal dynamic is an important factor when studying the impact of distributional shifts on patterns of genetic data. Our results also demonstrate the utility of SDMs to generate species-specific predictions about patterns of genetic variation that account for varying degrees of habitat specialization and life history characteristics of taxa. Nevertheless, the results highlight some key issues when converting SDMs for use in demographic models. Because the transformations have direct effects on the genetic consequence of population expansion by prescribing how habitat heterogeneity and spatiotemporal variation is related to the species-specific demographic model, it is important to consider alternative transformations when studying the genetic consequences of distributional shifts.