This content will become publicly available on June 9, 2024
- Award ID(s):
- 2215705
- NSF-PAR ID:
- 10464600
- 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
-
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
-
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 of
Clavispora (Candida )lusitaniae isolates, 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 inCandida species. 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. lusitaniae populations 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. -
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 of
Trachymyrmex arizonensis andTrachymyrmex pomonae and 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 betweenTrachymyrmex species 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. -
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 pathogen
Cladosporium fulvum , which is the first fungal plant pathogen from which avirulence (Avr ) genes were ever cloned and in which loss ofAvr genes is often reported as a means of overcoming recognition by cognate tomatoCf resistance genes. A recent near-complete reference genome assembly ofC. fulvum isolate 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 additional
C. fulvum isolates. 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 the
C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens. -
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.