Search for: All records

Most hosts contain few parasites, whereas few hosts contain many. This pattern, known as aggregation, is well-documented in macroparasites where parasite intensity distribution among hosts affects host–parasite dynamics. Infection intensity also drives fungal disease dynamics, but we lack a basic understanding of host–fungal aggregation patterns, how they compare with macroparasites and if they reflect biological processes. To begin addressing these gaps, we characterized aggregation of the fungal pathogenBatrachochytrium dendrobatidis(Bd) in amphibian hosts. Utilizing the slope of Taylor’s Power law, we found Bd intensity distributions were more aggregated than many macroparasites, conforming closely to lognormal distributions. We observed that Bd aggregation patterns are strongly correlated with known biological processes operating in amphibian populations, such as epizoological phase (i.e. invasion, post-invasion and enzootic), and intensity-dependent disease mortality. Using intensity-dependent mathematical models, we found evidence of evolution of host resistance based on aggregation shifts in systems persisting with Bd following disease-induced declines. Our results show that Bd aggregation is highly conserved across disparate systems and contains signatures of potential biological processes of amphibian–Bd systems. Our work can inform future modelling approaches and be extended to other fungal pathogens to elucidate host–fungal interactions and unite host–fungal dynamics under a common theoretical framework. 

Insights from conservation genomics have dramatically improved recovery plans for numerous endangered species. However, most taxa have yet to benefit from the full application of genomic technologies. The mountain yellow-legged frog species complex, Rana muscosa and Rana sierrae, inhabits the Sierra Nevada mountains and Transverse/Peninsular Ranges of California and Nevada. Both species have declined precipitously throughout their historical distributions. Conservation management plans outline extensive ongoing recovery efforts but are still based on the genetic structure determined primarily using a single mitochondrial sequence. Our study used two different sequencing strategies – amplicon sequencing and exome capture – to refine our understanding of the population genetics of these imperiled amphibians. We used buccal swabs, museum tissue samples, and archived skin swabs to genotype frog populations across their range. Using the amplicon sequencing and exome capture datasets separately and combined, we document five major genetic clusters. Notably, we found evidence supporting previous species boundaries within Kings Canyon National Park with some exceptions at individual sites. Though we see evidence of genetic clustering, especially in the R. muscosa clade, we also found evidence of some admixture across cluster boundaries in the R. sierrae clade, suggesting a stepping-stone model of population structure. We also find that the southern R. muscosa cluster had large runs of homozygosity and the lowest overall heterozygosity of any of the clusters, consistent with previous reports of marked declines in this area. Overall, our results clarify management unit designations across the range of an endangered species and highlight the importance of sampling the entire range of a species, even when collecting genome-scale data. 

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 Bd ASIA3, that appears to be widespread in Southeast Asia. This lineage co-occurs with the global panzootic lineage ( Bd GPL) in multiple localities. Additionally, we shed light on the global distribution of Bd GPL and highlight the expanded range of another lineage, Bd CAPE. 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.