Repeated evolution of blanched coloration in a lizard across independent white-sand habitats.

The White Sands lizards of New Mexico are a rare and classic example of convergent evolution where three species have evolved blanched coloration on the white gypsum dunes. Until now, no geological replicate of the pattern had been described. However, one of the White Sands species, the lesser earless lizard (Holbrookia maculata), has been discovered to also inhabit the Salt Basin Dunes of Texas, where it has also evolved a blanched morph. We here present a first phenotypic and genetic description of the Salt Basin Dunes population of H. maculata. Phylogenetic inference based on a housekeeping gene (ND4) and a classic candidate gene in the melanin-synthesis pathway (Melanocortin 1 Receptor; Mc1r) shows the newly discovered population as an independent lineage, with no evidence of genetic parallelism in the coding region of Mc1r. Initial morphological data suggest that while this population displays convergent evolution in blanched coloration, there are divergent patterns in limb length and habitat use behavior between the gypsum environments. Our findings present the White Sands/Salt Basin Dunes as an exceptionally promising comparative model for studies of adaptation and convergent evolution.

The fungal pathogen Batrachochytrium salamandrivorans is not detected in wild and captive amphibians from Mexico.

The recent emergence of the pathogen Batrachochytrium salamandrivorans (Bsal) is associated with rapid population declines of salamanders in Europe and its arrival to new areas could cause dramatic negative effects on other amphibian populations and species. Amphibian species, present in areas with high amphibian diversity such as Mexico, could be highly threatened due to the arrival of Bsal, particularly salamander species which are more vulnerable to chytridiomycosis caused by this pathogen. Thus, immediate surveillance is needed as a strategy to efficiently contend with this emerging infectious disease. In this study, we analyzed 490 wild and captive amphibians from 48 species across 76 sites in the North, Central, and South of Mexico to evaluate the presence of Bsal. Amphibians were sampled in sites with variable degrees of amphibian richness and suitability for Bsal according to previous studies. From the 76 sampling sites, 10 of them were located in areas with high amphibian richness and potential moderate to high Bsal habitat suitability. We did not detect Bsal in any of the samples, and no signs of the disease were observed in any individual at the time of sampling. Our results suggest that Bsal has not yet arrived at the sampled sites or could be at low prevalence within populations with low occurrence probability. This is the first study that evaluates the presence of Bsal in different regions and amphibian species in Mexico, which is the second most diverse country in salamander species in the world. We highlight the risk and the importance of continuing surveillance of Bsal in Mexico and discuss control strategies to avoid the introduction and spread of Bsal in the country.

Thermal Costs and Benefits of Replicated Color Evolution in the White Sands Desert Lizard Community.

AbstractTraits often contribute to multiple functions, complicating our understanding of the selective pressures that influence trait evolution. In the Chihuahuan Desert, predation is thought to be the primary driver of cryptic light coloration in three White Sands lizard species relative to the darker coloration of populations on adjacent dark soils. However, coloration also influences radiation absorption and thus animal body temperatures. We combined comparative physiological experiments and biophysical models to test for thermal consequences of evolving different color morphs in White Sands across the three species. While light and dark morphs have not evolved different physiological heat limits within species, differences in radiation absorption between morphs lead to body temperature differences that impact relative overheating risk and activity patterns. Moreover, for all three species, an idealized morph that matches the White Sands substrate would have considerably less activity time, by approximately 1 month, than existing light morphs. Overall, there are both benefits and costs to greater substrate matching, the balance of which may prevent the evolution of optimal crypsis. Our work highlights the importance of color in dictating thermal performance and the complexity inherent in understanding the evolution of coloration.

Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis).

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd), termed Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the site of origin for the common ancestor of modern Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swabs and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs (Rana catesbeiana) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential site of origin for Bd-GPL.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas).

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

La investigación de la distribución espacial de la variación genética y fenotípica puede proporcionar información sobre los procesos evolutivos que dan forma a la diversidad en los sistemas naturales. Caracterizamos patrones de diversidad genética y fenotípica para conocer los impulsores de la diversificación de patrones de color en ranas con ojos rojos (Agalychnis callidryas) en Costa Rica. A lo largo de la costa del Pacífico, las ranas con ojos rojos tienen un patrón de color llamativo en las patas que cambia de naranja en el norte a púrpura en el sur. Medimos la variación fenotípica de las ranas en los sitios del Pacífico, con un mayor muestreo en los sitios donde ocurre la transición de naranja a púrpura. En la zona de transición, descubrimos la co-ocurrencia de múltiples morfos de patrones de color. Para explorar las posibles causas de esta variación, generamos un conjunto de datos SNP con secuenciación RAD para analizar la estructura genética de la población, medir la diversidad genética e inferir los procesos que median la dinámica genotipo-fenotipo. Investigamos cómo los patrones de parentesco genético se corresponden con medidas individuales de patrón de color a lo largo de la costa, incluidas las pruebas del papel de la hibridación en regiones geográficas donde coexisten grupos fenotípicos naranja y morado. No encontramos evidencia de que el polimorfismo del patrón de color en la zona de transición surgiera a través de una hibridación o introgresión reciente. En cambio, un fuerte patrón de aislamiento genético por distancia indica que la variación del patrón de color se retuvo a través de otros procesos, como los polimorfismos de color ancestrales, el contacto secundario antiguo o la generada por mutaciones novedosas. Descubrimos que el fenotipo de color cambia a lo largo de la costa del Pacífico más de lo que se esperaría solo por la divergencia genética y la distancia geográfica. Combinados, nuestros resultados sugieren la posibilidad de que las presiones selectivas actúen sobre el patrón de color a pequeña escala geográfica.

The recombination landscapes of spiny lizards (genus Sceloporus).

Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles-the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.

Skin bacterial metacommunities of San Francisco Bay Area salamanders are structured by host genus and habitat quality.

Host-associated microbial communities can influence physiological processes of macroorganisms, including contributing to infectious disease resistance. For instance, some bacteria that live on amphibian skin produce antifungal compounds that inhibit two lethal fungal pathogens, Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Therefore, differences in microbiome composition among host species or populations within a species can contribute to variation in susceptibility to Bd/Bsal. This study applies 16S rRNA sequencing to characterize the skin bacterial microbiomes of three widespread terrestrial salamander genera native to the western United States. Using a metacommunity structure analysis, we identified dispersal barriers for these influential bacteria between salamander families and localities. We also analysed the effects of habitat characteristics such as percent natural cover and temperature seasonality on the microbiome. We found that certain environmental variables may influence the skin microbial communities of some salamander genera more strongly than others. Each salamander family had a somewhat distinct community of putative anti-Bd skin bacteria, suggesting that salamanders may select for a functional assembly of cutaneous symbionts that could differ in its ability to protect these amphibians from disease. Our observations raise the need to consider host identity and environmental heterogeneity during the selection of probiotics to treat wildlife diseases.

Genetic variation of Batrachochytrium dendrobatidis is linked to skin bacterial diversity in the Pacific treefrog Hyliola regilla (hypochondriaca).

Symbiotic bacterial communities are crucial to combating infections and contribute to host health. The amphibian skin microbiome plays an important role in protecting their hosts against pathogens such as Batrachochytrium dendrobatidis (Bd), one of the causative agents of chytridiomycosis, which is responsible for dramatic amphibian population declines worldwide. Although symbiotic skin bacteria are known to inhibit Bd growth, an understanding of the relationship between Bd genetic variability, environmental conditions, and skin bacterial communities is limited. Therefore, we examined the associations between Bd infection load, Bd genetic diversity and skin bacterial communities in five populations of Hyliola regilla (hypochondriaca) from environmentally contrasting sites in Baja California, Mexico. We observed differences in Bd genetics and infection load among sites and environments. Genetic analysis of Bd isolates revealed patterns of spatial structure corresponding to the five sites sampled. Amphibian skin bacterial diversity and community structure differed among environments and sites. Bacterial community composition was correlated with Bd genetic differences and infection load, with specific bacterial taxa enriched on infected and un-infected frogs. Our results indicate that skin-associated bacteria and Bd strains likely interact on the host skin, with consequences for microbial community structure and Bd infection intensity.

Thermal Performance Curves of Multiple Isolates of Batrachochytrium dendrobatidis, a Lethal Pathogen of Amphibians.

Emerging infectious disease is a key factor in the loss of amphibian diversity. In particular, the disease chytridiomycosis has caused severe declines around the world. The lethal fungal pathogen that causes chytridiomycosis, Batrachochytrium dendrobatidis (Bd), has affected amphibians in many different environments. One primary question for researchers grappling with disease-induced losses of amphibian biodiversity is what abiotic factors drive Bd pathogenicity in different environments. To study environmental influences on Bd pathogenicity, we quantified responses of Bd phenotypic traits (e.g., viability, zoospore densities, growth rates, and carrying capacities) over a range of environmental temperatures to generate thermal performance curves. We selected multiple Bd isolates that belong to a single genetic lineage but that were collected across a latitudinal gradient. For the population viability, we found that the isolates had similar thermal optima at 21°C, but there was considerable variation among the isolates in maximum viability at that temperature. Additionally, we found the densities of infectious zoospores varied among isolates across all temperatures. Our results suggest that temperatures across geographic point of origin (latitude) may explain some of the variation in Bd viability through vertical shifts in maximal performance. However, the same pattern was not evident for other reproductive parameters (zoospore densities, growth rates, fecundity), underscoring the importance of measuring multiple traits to understand variation in pathogen responses to environmental conditions. We suggest that variation among Bd genetic variants due to environmental factors may be an important determinant of disease dynamics for amphibians across a range of diverse environments.

Divergent regional evolutionary histories of a devastating global amphibian pathogen.

Emerging infectious diseases are a pressing threat to global biological diversity. Increased incidence and severity of novel pathogens underscores the need for methodological advances to understand pathogen emergence and spread. Here, we use genetic epidemiology to test, and challenge, key hypotheses about a devastating zoonotic disease impacting amphibians globally. Using an amplicon-based sequencing method and non-invasive samples we retrospectively explore the history of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in two emblematic amphibian systems: the Sierra Nevada of California and Central Panama. The hypothesis in both regions is the hypervirulent Global Panzootic Lineage of Bd (BdGPL) was recently introduced and spread rapidly in a wave-like pattern. Our data challenge this hypothesis by demonstrating similar epizootic signatures can have radically different underlying evolutionary histories. In Central Panama, our genetic data confirm a recent and rapid pathogen spread. However, BdGPL in the Sierra Nevada has remarkable spatial structuring, high genetic diversity and a relatively older history inferred from time-dated phylogenies. Thus, this deadly pathogen lineage may have a longer history in some regions than assumed, providing insights into its origin and spread. Overall, our results highlight the importance of integrating observed wildlife die-offs with genetic data to more accurately reconstruct pathogen outbreaks.

Early presence of Batrachochytrium dendrobatidis in Mexico with a contemporary dominance of the global panzootic lineage.

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is a devastating infectious disease of amphibians. Retrospective studies using museum vouchers and genetic samples supported the hypothesis that Bd colonized Mexico from North America and then continued to spread into Central and South America, where it led to dramatic losses in tropical amphibian biodiversity (the epizootic wave hypothesis). While these studies suggest that Bd has been in Mexico since the 1970s, information regarding the historical and contemporary occurrence of different pathogen genetic lineages across the country is limited. In the current study, we investigated the historical and contemporary patterns of Bd in Mexico. We combined the swabbing of historical museum vouchers and sampling of wild amphibians with a custom Bd genotyping assay to assess the presence, prevalence, and genetic diversity of Bd over time in Mexico. We found Bd-positive museum specimens from the late 1800s, far earlier than previous records and well before recent amphibian declines. With Bd genotypes from samples collected between 1975-2019, we observed a contemporary dominance of the global panzootic lineage in Mexico and report four genetic subpopulations and potential for admixture among these populations. The observed genetic variation did not have a clear geographic signature or provide clear support for the epizootic wave hypothesis. These results provide a framework for testing new questions regarding Bd invasions and their temporal relationship to observed amphibian declines in the Americas.

Whole exome sequencing identifies the potential for genetic rescue in iconic and critically endangered Panamanian harlequin frogs.

Avoiding extinction in a rapidly changing environment often relies on a species' ability to quickly adapt in the face of extreme selective pressures. In Panamá, two closely related harlequin frog species (Atelopus varius and Atelopus zeteki) are threatened with extinction due to the fungal pathogen Batrachochytrium dendrobatidis (Bd). Once thought to be nearly extirpated from Panamá, A. varius have recently been rediscovered in multiple localities across their historical range; however, A. zeteki are possibly extinct in the wild. By leveraging a unique collection of 186 Atelopus tissue samples collected before and after the Bd outbreak in Panama, we describe the genetics of persistence for these species on the brink of extinction. We sequenced the transcriptome and developed an exome-capture assay to sequence the coding regions of the Atelopus genome. Using these genetic data, we evaluate the population genetic structure of historical A. varius and A. zeteki populations, describe changes in genetic diversity over time, assess the relationship between contemporary and historical individuals, and test the hypothesis that some A. varius populations have rapidly evolved to resist or tolerate Bd infection. We found a significant decrease in genetic diversity in contemporary (compared to historical) A. varius populations. We did not find strong evidence of directional allele frequency change or selection for Bd resistance genes, but we uncovered a set of candidate genes that warrant further study. Additionally, we found preliminary evidence of recent migration and gene flow in one of the largest persisting A. varius populations in Panamá, suggesting the potential for genetic rescue in this system. Finally, we propose that previous conservation units should be modified, as clear genetic breaks do not exist beyond the local population level. Our data lay the groundwork for genetically informed conservation and advance our understanding of how imperiled species might be rescued from extinction.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog.

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation Fn individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

Stepping into the past to conserve the future: Archived skin swabs from extant and extirpated populations inform genetic management of an endangered amphibian.

Moving animals on a landscape through translocations and reintroductions is an important management tool used in the recovery of endangered species, particularly for the maintenance of population genetic diversity and structure. Management of imperiled amphibian species rely heavily on translocations and reintroductions, especially for species that have been brought to the brink of extinction by habitat loss, introduced species, and disease. One striking example of amphibian declines and associated management efforts is in California's Sequoia and Kings Canyon National Parks with the mountain yellow-legged frog species complex (Rana sierrae/muscosa). Mountain yellow-legged frogs have been extirpated from more than 93% of their historic range, and limited knowledge of their population genetics has made long-term conservation planning difficult. To address this, we used 598 archived skin swabs from both extant and extirpated populations across 48 lake basins to generate a robust Illumina-based nuclear amplicon data set. We found that samples grouped into three main genetic clusters, concordant with watershed boundaries. We also found evidence for historical gene flow across watershed boundaries with a north-to-south axis of migration. Finally, our results indicate that genetic diversity is not significantly different between populations with different disease histories. Our study offers specific management recommendations for imperiled mountain yellow-legged frogs and, more broadly, provides a population genetic framework for leveraging minimally invasive samples for the conservation of threatened species.

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.

Local adaptation does not lead to genome-wide differentiation in lava flow lizards.

Adaptation can occur with or without genome-wide differentiation. If adaptive loci are linked to traits involved in reproductive isolation, genome-wide divergence is likely, and speciation is possible. However, adaptation can also lead to phenotypic differentiation without genome-wide divergence if levels of ongoing gene flow are high. Here, we use the replicated occurrence of melanism in lava flow lizards to assess the relationship between local adaptation and genome-wide differentiation. We compare patterns of phenotypic and genomic divergence among lava flow and nonlava populations for three lizard species and three lava flows in the Chihuahuan Desert. We find that local phenotypic adaptation (melanism) is not typically accompanied by genome-wide differentiation. Specifically, lava populations do not generally exhibit greater divergence from nonlava populations than expected by geography alone, regardless of whether the lava formation is 5,000 or 760,000 years old. We also infer that gene flow between lava and nonlava populations is ongoing in all lava populations surveyed. Recent work in the isolation by environment and ecological speciation literature suggests that environmentally driven genome-wide differentiation is common in nature. However, local adaptation may often simply be local adaptation rather than an early stage of ecological speciation.

Phenotypic and genetic diversity in aposematic Malagasy poison frogs (genus Mantella).

Intraspecific color variation has long fascinated evolutionary biologists. In species with bright warning coloration, phenotypic diversity is particularly compelling because many factors, including natural and sexual selection, contribute to intraspecific variation. To better understand the causes of dramatic phenotypic variation in Malagasy poison frogs, we quantified genetic structure and color and pattern variation across three closely related species, Mantella aurantiaca, Mantella crocea, and Mantella milotympanum. Although our restriction site-associated DNA (RAD) sequencing approach identified clear genetic clusters, they do not align with current species designations, which has important conservation implications for these imperiled frogs. Moreover, our results suggest that levels of intraspecific color variation within this group have been overestimated, while species diversity has been underestimated. Within major genetic clusters, we observed distinct patterns of variation including: populations that are phenotypically similar yet genetically distinct, populations where phenotypic and genetic breaks coincide, and populations that are genetically similar but have high levels of within-population phenotypic variation. We also detected admixture between two of the major genetic clusters. Our study suggests that several mechanisms-including hybridization, selection, and drift-are contributing to phenotypic diversity. Ultimately, our work underscores the need for a reevaluation of how polymorphic and polytypic populations and species are classified, especially in aposematic organisms.

Opening the file drawer: Unexpected insights from a chytrid infection experiment.

Infection experiments are critical for understanding wildlife disease dynamics. Although infection experiments are typically designed to reduce complexity, disease outcomes still result from complex interactions between host, pathogen, and environmental factors. Cryptic variation across factors can lead to decreased repeatability of infection experiments within and between research groups and hinder research progress. Furthermore, studies with unexpected results are often relegated to the "file drawer" and potential insights gained from these experimental outcomes are lost. Here, we report unexpected results from an infection experiment studying the response of two differentially-susceptible but related frogs (American Bullfrog Rana catesbeiana and the Mountain yellow-legged frog Rana muscosa) to the amphibian-killing chytrid fungus (Batrachochytrium dendrobatidis, Bd). Despite well-documented differences in susceptibility between species, we found no evidence for antibody-mediated immune response and no Bd-related mortality in either species. Additionally, during the study, the sham-inoculated R. catesbeiana control group became unexpectedly Bd-positive. We used a custom genotyping assay to demonstrate that the aberrantly-infected R. catesbeiana carried a Bd genotype distinct from the inoculation genotype. Thus R. catesbeiana individuals were acquired with low-intensity infections that could not be detected with qPCR. In the Bd-inoculated R. catesbeiana treatment group, the inoculated genotype appeared to out-compete the cryptic infection. Thus, our results provide insight into Bd coinfection dynamics, a phenomenon that is increasingly relevant as different pathogen strains are moved around the globe. Our experiment highlights how unexpected experimental outcomes can serve as both cautionary tales and opportunities to explore unanswered research questions. We use our results as a case study to highlight common sources of anomalous results for infection experiments. We argue that understanding these factors will aid researchers in the design, execution, and interpretation of experiments to understand wildlife disease processes.

Shifts in disease dynamics in a tropical amphibian assemblage are not due to pathogen attenuation.

Infectious diseases rarely end in extinction. Yet the mechanisms that explain how epidemics subside are difficult to pinpoint. We investigated host-pathogen interactions after the emergence of a lethal fungal pathogen in a tropical amphibian assemblage. Some amphibian host species are recovering, but the pathogen is still present and is as pathogenic today as it was almost a decade ago. In addition, some species have defenses that are more effective now than they were before the epidemic. These results suggest that host recoveries are not caused by pathogen attenuation and may be due to shifts in host responses. Our findings provide insights into the mechanisms underlying disease transitions, which are increasingly important to understand in an era of emerging infectious diseases and unprecedented global pandemics.

Convergent Phenotypic Evolution despite Contrasting Demographic Histories in the Fauna of White Sands.

When are evolutionary outcomes predictable? Cases of convergent evolution can shed light on when, why, and how different species exhibit shared evolutionary trajectories. In particular, studying diverse species in a common environment can illuminate how different factors facilitate or constrain adaptive evolution. Here we integrate studies of pattern and process in the fauna at White Sands (New Mexico) to understand the determinants of convergent evolution. Numerous animal species at White Sands exhibit phenotypic convergence in response to a novel-and shared-selective environment: geologically young gypsum dunes. We synthesize 15 years of research on White Sands lizards to assess the contribution of natural selection, genetic architecture, and population demography to patterns of phenotypic evolution. We also present new data for two species of White Sands arthropods, Ammobaenetes arenicolus and Habronattus ustulatus. Overall, we find dramatic phenotypic convergence across diverse species at White Sands. Although the direction of phenotypic response is parallel, the magnitude of phenotypic response varies among species. We also find that species exhibit strikingly different demographic patterns across the ecotone. The species with the most genetic structure between White Sands and dark-soil populations generally exhibit the least phenotypic divergence, suggesting population demography as a key modulator of adaptation. Comparative studies are particularly important for understanding the determinants of convergence in natural systems.