skip to main content


This content will become publicly available on June 9, 2024

Title: In host evolution of Exophiala dermatitidis in cystic fibrosis lung micro-environment
Abstract Individuals with cystic fibrosis (CF) are susceptible to chronic lung infections that lead to inflammation and irreversible lung damage. While most respiratory infections that occur in CF are caused by bacteria, some are dominated by fungi such as the slow-growing black yeast Exophiala dermatitidis. Here, we analyze isolates of E. dermatitidis cultured from two samples, collected from a single subject 2 years apart. One isolate genome was sequenced using long-read Nanopore technology as an in-population reference to use in comparative single nucleotide polymorphism and insertion–deletion variant analyses of 23 isolates. We then used population genomics and phylo-genomics to compare the isolates to each other as well as the reference genome strain E. dermatitidis NIH/UT8656. Within the CF lung population, three E. dermatitidis clades were detected, each with varying mutation rates. Overall, the isolates were highly similar suggesting that they were recently diverged. All isolates were MAT 1-1, which was consistent with their high relatedness and the absence of evidence for mating or recombination between isolates. Phylogenetic analysis grouped sets of isolates into clades that contained isolates from both early and late time points indicating there are multiple persistent lineages. Functional assessment of variants unique to each clade identified alleles in genes that encode transporters, cytochrome P450 oxidoreductases, iron acquisition, and DNA repair processes. Consistent with the genomic heterogeneity, isolates showed some stable phenotype heterogeneity in melanin production, subtle differences in antifungal minimum inhibitory concentrations, and growth on different substrates. The persistent population heterogeneity identified in lung-derived isolates is an important factor to consider in the study of chronic fungal infections, and the analysis of changes in fungal pathogens over time may provide important insights into the physiology of black yeasts and other slow-growing fungi in vivo.  more » « less
Award ID(s):
2215705
NSF-PAR ID:
10464600
Author(s) / Creator(s):
; ; ; ;
Editor(s):
Nowrousian, M
Date Published:
Journal Name:
G3: Genes, Genomes, Genetics
Volume:
13
Issue:
8
ISSN:
2160-1836
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Management of the limited number of antimicrobials currently available requires the identification of infections that contain drug-resistant isolates and the discovery of factors that promote the evolution of drug resistance. Here, we report a single fungal infection in which we have identified numerous subpopulations that differ in their alleles of a single gene that impacts drug resistance. The diversity at this locus was markedly greater than the reported heterogeneity of alleles conferring antibiotic resistance in bacterial infections. Analysis of genomes from hundreds ofClavispora(Candida)lusitaniaeisolates, through individual and pooled isolate sequencing, from a single individual with cystic fibrosis revealed at least 25 nonsynonymous mutations inMRR1, which encodes a transcription factor capable of inducing fluconazole (FLZ) resistance inCandidaspecies. Isolates with high-activity Mrr1 variants were resistant to FLZ due to elevated expression of theMDR1-encoded efflux pump. We found that high Mrr1-regulated Mdr1 activity protected against host and bacterial factors, suggesting drug resistance can be selected for indirectly and perhaps explaining the Mrr1 heterogeneity in this individual who had no prior azole exposure. Regional analysis ofC. lusitaniaepopulations from the upper and lower lobes of the right lung suggested intermingling of subpopulations throughout. Our retrospective characterization of sputum and lung populations by pooled sequencing found that alleles that confer FLZ resistance were a minority in each pool, possibly explaining why they were undetected before unsuccessful FLZ therapy. New susceptibility testing regimes may detect problematical drug-resistant subpopulations in heterogeneous single-species infections.

     
    more » « less
  3. Abstract

    Over the past few decades, large‐scale phylogenetic analyses of fungus‐gardening ants and their symbiotic fungi have depicted strong concordance among major clades of ants and their symbiotic fungi, yet within clades, fungus sharing is widespread among unrelated ant lineages. Sharing has been explained using a diffuse coevolution model within major clades. Understanding horizontal exchange within clades has been limited by conventional genetic markers that lack both interspecific and geographic variation. To examine whether reports of horizontal exchange were indeed due to symbiont sharing or the result of employing relatively uninformative molecular markers, samples ofTrachymyrmex arizonensisandTrachymyrmex pomonaeand their fungi were collected from native populations in Arizona and genotyped using conventional marker genes and genome‐wide single nucleotide polymorphisms (SNPs). Conventional markers of the fungal symbionts generally exhibited cophylogenetic patterns that were consistent with some symbiont sharing, but most fungal clades had low support. SNP analysis, in contrast, indicated that each ant species exhibited fidelity to its own fungal subclade with only one instance of a colony growing a fungus that was otherwise associated with a different ant species. This evidence supports a pattern of codivergence betweenTrachymyrmexspecies and their fungi, and thus a diffuse coevolutionary model may not accurately predict symbiont exchange. These results suggest that fungal sharing across host species in these symbioses may be less extensive than previously thought.

     
    more » « less
  4. Abstract Background

    Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogenCladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss ofAvrgenes is often reported as a means of overcoming recognition by cognate tomatoCfresistance genes. A recent near-complete reference genome assembly ofC. fulvumisolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

    Results

    Here, we obtained near-complete genome assemblies of four additionalC. fulvumisolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genesAvr4E,Avr5, andAvr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome ofC. fulvum.

    Conclusions

    Our results reveal new evolutionary aspects of theC. fulvumgenome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

     
    more » « less
  5. Abstract

    Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

     
    more » « less