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.

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.

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.

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.

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.

Large-scale recovery of an endangered amphibian despite ongoing exposure to multiple stressors.

Amphibians are one of the most threatened animal groups, with 32% of species at risk for extinction. Given this imperiled status, is the disappearance of a large fraction of the Earth's amphibians inevitable, or are some declining species more resilient than is generally assumed? We address this question in a species that is emblematic of many declining amphibians, the endangered Sierra Nevada yellow-legged frog (Rana sierrae). Based on >7,000 frog surveys conducted across Yosemite National Park over a 20-y period, we show that, after decades of decline and despite ongoing exposure to multiple stressors, including introduced fish, the recently emerged disease chytridiomycosis, and pesticides, R. sierrae abundance increased sevenfold during the study and at a rate of 11% per year. These increases occurred in hundreds of populations throughout Yosemite, providing a rare example of amphibian recovery at an ecologically relevant spatial scale. Results from a laboratory experiment indicate that these increases may be in part because of reduced frog susceptibility to chytridiomycosis. The disappearance of nonnative fish from numerous water bodies after cessation of stocking also contributed to the recovery. The large-scale increases in R. sierrae abundance that we document suggest that, when habitats are relatively intact and stressors are reduced in their importance by active management or species' adaptive responses, declines of some amphibians may be partially reversible, at least at a regional scale. Other studies conducted over similarly large temporal and spatial scales are critically needed to provide insight and generality about the reversibility of amphibian declines at a global scale.

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.

Diversity in growth patterns among strains of the lethal fungal pathogen Batrachochytrium dendrobatidis across extended thermal optima.

The thermal sensitivities of organisms regulate a wide range of ecological interactions, including host-parasite dynamics. The effect of temperature on disease ecology can be remarkably complex in disease systems where the hosts are ectothermic and where thermal conditions constrain pathogen reproductive rates. Amphibian chytridiomycosis, caused by the pathogen Batrachochytrium dendrobatidis (Bd), is a lethal fungal disease that is influenced by temperature. However, recent temperature studies have produced contradictory findings, suggesting that our current understanding of thermal effects on Bd may be incomplete. We investigated how temperature affects three different Bd strains to evaluate diversity in thermal responses. We quantified growth across the entire thermal range of Bd, and beyond the known thermal limits (T max and T min). Our results show that all Bd strains remained viable and grew following 24 h freeze (-12 °C) and heat shock (28 °C) treatments. Additionally, we found that two Bd strains had higher logistic growth rates (r) and carrying capacities (K) at the upper and lower extremities of the temperature range, and especially in low temperature conditions (2-3 °C). In contrast, a third strain exhibited relatively lower growth rates and carrying capacities at these same thermal extremes. Overall, our results suggest that there is considerable variation among Bd strains in thermal tolerance, and they establish a new thermal sensitivity profile for Bd. More generally, our findings point toward important questions concerning the mechanisms that dictate fungal thermal tolerances and temperature-dependent pathogenesis in other fungal disease systems.

The population genomics of rapid adaptation: disentangling signatures of selection and demography in white sands lizards.

Understanding the process of adaptation during rapid environmental change remains one of the central focal points of evolutionary biology. The recently formed White Sands system of southern New Mexico offers an outstanding example of rapid adaptation, with a variety of species having rapidly evolved blanched forms on the dunes that contrast with their close relatives in the surrounding dark soil habitat. In this study, we focus on two of the White Sands lizard species, Sceloporus cowlesi and Aspidoscelis inornata, for which previous research has linked mutations in the melanocortin-1 receptor gene (Mc1r) to blanched coloration. We sampled populations both on and off the dunes and used a custom sequence capture assay based on probed fosmid libraries to obtain >50 kb of sequence around Mc1r and hundreds of other random genomic locations. We then used model-based statistical inference methods to describe the demographic and adaptive history characterizing the colonization of White Sands. We identified a number of similarities between the two focal species, including strong evidence of selection in the blanched populations in the Mc1r region. We also found important differences between the species, suggesting different colonization times, different genetic architecture underlying the blanched phenotype and different ages of the beneficial alleles. Finally, the beneficial allele is dominant in S. cowlesi and recessive in A. inornata, allowing for a rare empirical test of theoretically expected patterns of selective sweeps under these differing models.

Complex history of the amphibian-killing chytrid fungus revealed with genome resequencing data.

Understanding the evolutionary history of microbial pathogens is critical for mitigating the impacts of emerging infectious diseases on economically and ecologically important host species. We used a genome resequencing approach to resolve the evolutionary history of an important microbial pathogen, the chytrid Batrachochytrium dendrobatidis (Bd), which has been implicated in amphibian declines worldwide. We sequenced the genomes of 29 isolates of Bd from around the world, with an emphasis on North, Central, and South America because of the devastating effect that Bd has had on amphibian populations in the New World. We found a substantial amount of evolutionary complexity in Bd with deep phylogenetic diversity that predates observed global amphibian declines. By investigating the entire genome, we found that even the most recently evolved Bd clade (termed the global panzootic lineage) contained more genetic variation than previously reported. We also found dramatic differences among isolates and among genomic regions in chromosomal copy number and patterns of heterozygosity, suggesting complex and heterogeneous genome dynamics. Finally, we report evidence for selection acting on the Bd genome, supporting the hypothesis that protease genes are important in evolutionary transitions in this group. Bd is considered an emerging pathogen because of its recent effects on amphibians, but our data indicate that it has a complex evolutionary history that predates recent disease outbreaks. Therefore, it is important to consider the contemporary effects of Bd in a broader evolutionary context and identify specific mechanisms that may have led to shifts in virulence in this system.

Genomic transition to pathogenicity in chytrid fungi.

Understanding the molecular mechanisms of pathogen emergence is central to mitigating the impacts of novel infectious disease agents. The chytrid fungus Batrachochytrium dendrobatidis (Bd) is an emerging pathogen of amphibians that has been implicated in amphibian declines worldwide. Bd is the only member of its clade known to attack vertebrates. However, little is known about the molecular determinants of - or evolutionary transition to - pathogenicity in Bd. Here we sequence the genome of Bd's closest known relative - a non-pathogenic chytrid Homolaphlyctis polyrhiza (Hp). We first describe the genome of Hp, which is comparable to other chytrid genomes in size and number of predicted proteins. We then compare the genomes of Hp, Bd, and 19 additional fungal genomes to identify unique or recent evolutionary elements in the Bd genome. We identified 1,974 Bd-specific genes, a gene set that is enriched for protease, lipase, and microbial effector Gene Ontology terms. We describe significant lineage-specific expansions in three Bd protease families (metallo-, serine-type, and aspartyl proteases). We show that these protease gene family expansions occurred after the divergence of Bd and Hp from their common ancestor and thus are localized to the Bd branch. Finally, we demonstrate that the timing of the protease gene family expansions predates the emergence of Bd as a globally important amphibian pathogen.

Molecular and functional basis of phenotypic convergence in white lizards at White Sands.

There are many striking examples of phenotypic convergence in nature, in some cases associated with changes in the same genes. But even mutations in the same gene may have different biochemical properties and thus different evolutionary consequences. Here we dissect the molecular mechanism of convergent evolution in three lizard species with blanched coloration on the gypsum dunes of White Sands, New Mexico. These White Sands forms have rapidly evolved cryptic coloration in the last few thousand years, presumably to avoid predation. We use cell-based assays to demonstrate that independent mutations in the same gene underlie the convergent blanched phenotypes in two of the three species. Although the same gene contributes to light phenotypes in these White Sands populations, the specific molecular mechanisms leading to reduced melanin production are different. In one case, mutations affect receptor signaling and in the other, the ability of the receptor to integrate into the melanocyte membrane. These functional differences have important ramifications at the organismal level. Derived alleles in the two species show opposite dominance patterns, which in turn affect their visibility to selection and the spatial distribution of alleles across habitats. Our results demonstrate that even when the same gene is responsible for phenotypic convergence, differences in molecular mechanism can have dramatic consequences on trait expression and ultimately the adaptive trajectory.

Only skin deep: shared genetic response to the deadly chytrid fungus in susceptible frog species.

Amphibian populations around the world are threatened by an emerging infectious pathogen, the chytrid fungus Batrachochytrium dendrobatidis (Bd). How can a fungal skin infection kill such a broad range of amphibian hosts? And do different host species have a similar response to Bd infection? Here, we use a genomics approach to understand the genetic response of multiple susceptible frog species to Bd infection. We characterize the transcriptomes of two closely related endangered frog species (Rana muscosa and Rana sierrae) and analyse whole genome expression profiles from frogs in controlled Bd infection experiments. We integrate the Rana results with a comparable data set from a more distantly related susceptible species (Silurana tropicalis). We demonstrate that Bd-infected frogs show massive disruption of skin function and show no evidence of a robust immune response. The genetic response to infection is shared across the focal susceptible species, suggesting a common effect of Bd on susceptible frogs.

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.

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.

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.