An official website of the United States government
Abstract The amphibian chytrid fungusBatrachochytrium dendrobatidis(Bd) is a cosmopolitan pathogen with numerous distinct lineages. The global panzootic lineage (Bd-GPL) is the most widespread and virulent lineage and is responsible for many recorded amphibian declines. Mapping the extent ofBd-GPL and other more established lineages is important for predicting disease dynamics in amphibian communities. Ecuador is the most biodiverse country per unit area for amphibian taxa and, thus, a priority for studies on genotypic diversity ofBd. In this study, we tested skin swab samples collected from 464 individual amphibians across coastal, Andean montane, and Amazonian forests, for the presence of twoBdlineages known to be present in South America: the globally-distributedBd-GPL and the Brazilian-endemicBd-Asia2/Brazil lineage. By using a discriminatory qPCR SNP assay, we found widespread prevalence ofBd-GPL in Ecuador in diverse host taxa. Genotyping efficiency was 36% in this study, meaning that one in every three swabs that tested positive forBdin infection assays were successfully genotyped. Through this study, we provide further support for the presence of a singleBdlineage in this neotropical biodiversity hotspot.
more »
« less
This content will become publicly available on May 30, 2026
Comparative Genomics Reveals Intra and Inter Species Variation in the Pathogenic Fungus Batrachochytrium dendrobatidis
Abstract The Global Panzootic Lineage (GPL) of Batrachochytrium dendrobatidis (Bd) has been described as a main driver of amphibian extinctions. Pathogen studies have benefited from three Bd-GPL strain genomes, but identifying the genetic and molecular features that distinguish the B. dendrobatidis lineages requires additional high-quality genomes from diverse lineages. We sequenced and assembled genomes with Oxford Nanopore Technologies to produce assemblies of three Bd-BRAZIL isolates and one nonpathogen outgroup species Polyrhizophydium stewartii. The Bd-BRAZIL assembly sizes ranged between 22.0 and 26.1 Mb with 8,495 to 8,620 predicted protein-coding genes. We sought to categorize the pangenome of the species by identifying homologous genes across the sampled genomes as either being core and present in all strains, or accessory and shared among strains in a lineage, an analysis that has not yet been conducted on B. dendrobatidis and its lineages. We identified a core genome consisting of 6,278 gene families, and an accessory genome of 202 Bd-BRAZIL and 172 Bd-GPL specific gene families. We discovered copy number differences in pathogenicity gene families: M36 Peptidases, Crinkler Necrosis genes, Aspartyl Peptidases, Carbohydrate-Binding Module-18 genes, and S41 Proteases, between Bd-BRAZIL and Bd-GPL strains. Comparison of B. dendrobatidis and two closely related saprophytic species identified differences in protein sequence and domain counts for M36 and CBM18 families respectively. Our pangenome analysis of lineage-specific gene content led us to explore how the selection of the reference genome affects recovery of RNAseq transcripts when comparing different strains. We tested the hypothesis that genomic variation among Bd-GPL and Bd-BRAZIL lineages can impact transcript count data by comparing results with our new Bd-BRAZIL genomes as the reference genomes. Our analysis examines the genomic variation between strains in Bd-BRAZIL and Bd-GPL and offers insights into the application of these high-quality reference genomes resources for future studies.
more »
« less
- Award ID(s):
- 2215705
- PAR ID:
- 10656446
- Editor(s):
- Giraud, Tatiana
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Genome Biology and Evolution
- Volume:
- 17
- Issue:
- 6
- ISSN:
- 1759-6653
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Sil, Anita (Ed.)Aspergillus fumigatus is a deadly agent of human fungal disease where virulence heterogeneity is thought to be at least partially structured by genetic variation between strains. While population genomic analyses based on reference genome alignments offer valuable insights into how gene variants are distributed across populations, these approaches fail to capture intraspecific variation in genes absent from the reference genome. Pan-genomic analyses based on de novo assemblies offer a promising alternative to reference-based genomics with the potential to address the full genetic repertoire of a species. Here, we evaluate 260 genome sequences of A . fumigatus including 62 newly sequenced strains, using a combination of population genomics, phylogenomics, and pan-genomics. Our results offer a high-resolution assessment of population structure and recombination frequency, phylogenetically structured gene presence–absence variation, evidence for metabolic specificity, and the distribution of putative antifungal resistance genes. Although A . fumigatus disperses primarily via asexual conidia, we identified extraordinarily high levels of recombination with the lowest linkage disequilibrium decay value reported for any fungal species to date. We provide evidence for 3 primary populations of A . fumigatus , with recombination occurring only rarely between populations and often within them. These 3 populations are structured by both gene variation and distinct patterns of gene presence–absence with unique suites of accessory genes present exclusively in each clade. Accessory genes displayed functional enrichment for nitrogen and carbohydrate metabolism suggesting that populations may be stratified by environmental niche specialization. Similarly, the distribution of antifungal resistance genes and resistance alleles were often structured by phylogeny. Altogether, the pan-genome of A . fumigatus represents one of the largest fungal pan-genomes reported to date including many genes unrepresented in the Af293 reference genome. These results highlight the inadequacy of relying on a single-reference genome-based approach for evaluating intraspecific variation and the power of combined genomic approaches to elucidate population structure, genetic diversity, and putative ecological drivers of clinically relevant fungi.more » « less
-
Abstract Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species’ divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species’ pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.more » « less
-
Considerable attention has been directed to studying the infection dynamics of the fungal pathogenBatrachochytrium dendrobatidis(Bd) affecting amphibians in the high elevations of the Neotropics. However, lowland forests of the same realm remain comparatively understudied in this context. Herein, we attempt to bridge this gap by measuring the prevalence ofBdvia quantitative polymerase chain reaction (qPCR) in several anuran taxa inhabiting the Amazonian lowlands in the northeast of Ecuador. To this end, we sampled 207 anurans from 10 different families, 25 different genera, and 55 distinct host species originally collected in 2008. Taxonomy (at the family level), morphology (i.e. weight and snout-vent length), and life-long aquatic dependency of hosts (i.e. aquatic index) were also compiled to serve as potential predictors ofBdinfection status. Our findings revealed a relatively highBdprevalence of 58%, with 88% of sampled anuran families testing positive for the fungus at varying proportions. Model selection involving fitting and testing several different linear models, including mixed linear models, revealed a significant negative relationship between host weight andBdinfection status (p < 0.01). However, no significant associations were observed between taxonomy, aquatic dependency, snout-vent length, andBdinfections. In addition, we only detected the global panzootic lineage ofBd(Bd-GPL) and not theBd-Asia-2/Bd-Brazil lineage via qPCR single nucleotide polymorphism (SNP) genotyping. Our findings contribute to the understanding ofBddynamics in the Neotropical lowlands and emphasize the need for future research on the ecological factors influencingBdin the Amazon and their implications for amphibian conservation.more » « less
-
Cryptic diversity of a widespread global pathogen reveals expanded threats to amphibian conservationBiodiversity 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.more » « less
