skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The pan-genome of Aspergillus fumigatus provides a high-resolution view of its population structure revealing high levels of lineage-specific diversity driven by recombination
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
Award ID(s):
2215705
PAR ID:
10464605
Author(s) / Creator(s):
; ; ;
Editor(s):
Sil, Anita
Date Published:
Journal Name:
PLOS Biology
Volume:
20
Issue:
11
ISSN:
1545-7885
Page Range / eLocation ID:
e3001890
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Goldman, Gustavo H (Ed.)
    ABSTRACT Fungal infections are difficult to prevent and treat in large part due to strain heterogeneity, which confounds diagnostic predictability. Yet, the genetic mechanisms driving strain-to-strain variation remain poorly understood. Here, we determined the extent to whichStarships—giant transposons capable of mobilizing numerous fungal genes—generate genetic and phenotypic variability in the opportunistic human pathogenAspergillus fumigatus. We analyzed 519 diverse strains, including 11 newly sequenced with long-read technology and multiple isolates of the same reference strain, to reveal 20 distinctStarshipsthat are generating genomic heterogeneity over timescales relevant for experimental reproducibility.Starship-mobilized genes encode diverse functions, including known biofilm-related virulence factors and biosynthetic gene clusters, and many are differentially expressed during infection and antifungal exposure in a strain-specific manner. These findings support a new model of fungal evolution whereinStarshipshelp generate variation in genome structure, gene content, and expression among fungal strains. Together, our results demonstrate thatStarshipsare a previously hidden mechanism generating genotypic and, in turn, phenotypic heterogeneity in a major human fungal pathogen.IMPORTANCENo “one size fits all” option exists for treating fungal infections in large part due to genetic and phenotypic variability among strains. Accounting for strain heterogeneity is thus fundamental for developing efficacious treatments and strategies for safeguarding human health. Here, we report significant progress toward achieving this goal by uncovering a previously hidden mechanism generating heterogeneity in the human fungal pathogenAspergillus fumigatus: giant transposons, calledStarships, that span dozens of kilobases and mobilize fungal genes as cargo. By conducting a systematic investigation of these unusual transposons in a single fungal species, we demonstrate their contributions to population-level variation at the genome, pangenome, and transcriptome levels. TheStarshipcompendium we develop will not only help predict variation introduced by these elements in laboratory experiments but will serve as a foundational resource for determining howStarshipsimpact clinically relevant phenotypes, such as antifungal resistance and pathogenicity. 
    more » « less
  2. Giraud, Tatiana (Ed.)
    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
  3. Abstract Sage-grouse are two closely related iconic species of the North American West, with historically broad distributions across sagebrush-steppe habitat. Both species are dietary specialists on sagebrush during winter, with presumed adaptations to tolerate the high concentrations of toxic secondary metabolites that function as plant chemical defenses. Marked range contraction and declining population sizes since European settlement have motivated efforts to identify distinct population genetic variation, particularly that which might be associated with local genetic adaptation and dietary specialization of sage-grouse. We assembled a reference genome and performed whole-genome sequencing across sage-grouse from six populations, encompassing both species and including several populations on the periphery of the species ranges. Population genomic analyses reaffirmed genome-wide differentiation between greater and Gunnison sage-grouse, revealed pronounced intraspecific population structure, and highlighted important differentiation of a small isolated population of greater sage-grouse in the northwest of the range. Patterns of genome-wide differentiation were largely consistent with a hypothesized role of genetic drift due to limited gene flow among populations. Inferred ancient population demography suggested persistent declines in effective population sizes that have likely contributed to differentiation within and among species. Several genomic regions with single-nucleotide polymorphisms exhibiting extreme population differentiation were associated with candidate genes linked to metabolism of xenobiotic compounds. In vitro activity of enzymes isolated from sage-grouse livers supported a role for these genes in detoxification of sagebrush, suggesting that the observed interpopulation variation may underlie important local dietary adaptations, warranting close consideration for conservation strategies that link sage-grouse to the chemistry of local sagebrush. 
    more » « less
  4. null (Ed.)
    Abstract Background Antibiotic-producing Streptomyces bacteria are ubiquitous in nature, yet most studies of its diversity have focused on free-living strains inhabiting diverse soil environments and those in symbiotic relationship with invertebrates. Results We studied the draft genomes of 73 Streptomyces isolates sampled from the skin (wing and tail membranes) and fur surfaces of bats collected in Arizona and New Mexico. We uncovered large genomic variation and biosynthetic potential, even among closely related strains. The isolates, which were initially identified as three distinct species based on sequence variation in the 16S rRNA locus, could be distinguished as 41 different species based on genome-wide average nucleotide identity. Of the 32 biosynthetic gene cluster (BGC) classes detected, non-ribosomal peptide synthetases, siderophores, and terpenes were present in all genomes. On average, Streptomyces genomes carried 14 distinct classes of BGCs (range = 9–20). Results also revealed large inter- and intra-species variation in gene content (single nucleotide polymorphisms, accessory genes and singletons) and BGCs, further contributing to the overall genetic diversity present in bat-associated Streptomyces . Finally, we show that genome-wide recombination has partly contributed to the large genomic variation among strains of the same species. Conclusions Our study provides an initial genomic assessment of bat-associated Streptomyces that will be critical to prioritizing those strains with the greatest ability to produce novel antibiotics. It also highlights the need to recognize within-species variation as an important factor in genetic manipulation studies, diversity estimates and drug discovery efforts in Streptomyces . 
    more » « less
  5. ABSTRACT Yellow monkeyflowers (Mimulus guttatuscomplex, Phrymaceae) are a powerful system for studying ecological adaptation, reproductive variation, and genome evolution. To initiate pan‐genomics in this group, we present four chromosome‐scale assemblies and annotations of accessions spanning a broad evolutionary spectrum: two from a singleM. guttatuspopulation, one from the closely related selfing speciesM. nasutus, and one from a more divergent speciesM. tilingii. All assemblies are highly complete and resolve centromeric and repetitive regions. Comparative analyses reveal such extensive structural variation in repeat‐rich, gene‐poor regions that large portions of the genome are unalignable across accessions. As a result, thisMimuluspan‐genome is primarily informative in genic regions, underscoring limitations of resequencing approaches in such polymorphic taxa. We document gene presence–absence, investigate the recombination landscape using high‐resolution linkage data, and quantify nucleotide diversity. Surprisingly, pairwise differences at fourfold synonymous sites are exceptionally high—even in regions of very low recombination—reaching ~3.2% within a singleM. guttatuspopulation, ~7% within the interfertileM. guttatusspecies complex (approximately equal to SNP divergence between great apes and Old World monkeys), and ~7.4% between that complex and the reproductively isolatedM. tilingii. Genome‐wide patterns of nucleotide variation show little evidence of linked selection, and instead suggest that the concentration of genes (and likely selected sites) in high‐recombination regions may buffer diversity loss. These assemblies, annotations, and comparative analyses provide a robust genomic foundation forMimulusresearch and offer new insights into the interplay of recombination, structural variation, and molecular evolution in highly diverse plant genomes. 
    more » « less