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

Pathogenicity is associated with population structure in a fungal pathogen of humans

Aspergillus flavus is a clinically and agriculturally important saprotrophic fungus responsible for severe human infections and extensive crop losses. We analyzed genomic data from 250 (95 clinical and 155 environmental) A. flavus isolates from 9 countries, including 70 newly sequenced clinical isolates, to examine population and pan-genome structure and their relationship to pathogenicity. We identified five A. flavus populations, including a new population, D, corresponding to distinct clades in the genome-wide phylogeny. Strikingly, > 75% of clinical isolates were from population D. Accessory genes, including genes within biosynthetic gene clusters, were significantly more common in some populations but rare in others. Population D was enriched for genes associated with zinc ion binding, lipid metabolism, and certain types of hydrolase activity. In contrast to the major human pathogen Aspergillus fumigatus, A. flavus pathogenicity in humans is strongly associated with population structure, making it a great system for investigating how population-specific genes contribute to pathogenicity.

Genome-wide patterns of noncoding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Giant transposons promote strain heterogeneity in a major fungal pathogen

Fungal infections are difficult to prevent and treat in large part due to heterogeneity in clinically relevant phenotypes. However, the genetic mechanisms driving pathogen variation remain poorly understood. Here, we determined the extent to which Starships -giant transposons capable of mobilizing numerous fungal genes-generate genetic and phenotypic variability in the human pathogen Aspergillus fumigatus . We analyzed 519 diverse strains, including 12 newly sequenced with long-read technology, to reveal 20 distinct Starships that generate genomic heterogeneity over timescales impacting experimental reproducibility. Starship -mobilized genes encode diverse functions, including 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 pathogenesis wherein Starships mediate variation in virulence-related gene content and expression. Together, our results demonstrate that Starships are a foundational mechanism generating disease-relevant genotypic and, in turn, phenotypic heterogeneity in a major human fungal pathogen.

Azole resistance mechanisms and population structure of the human pathogen Aspergillus fumigatus on retail plant products

Aspergillus fumigatus is a ubiquitous saprotroph and human-pathogenic fungus that is life-threatening to the immunocompromised. Triazole-resistant A. fumigatus was found in patients without prior treatment with azoles, leading researchers to conclude that resistance had developed in agricultural environments where azoles are used against plant pathogens. Previous studies have documented azole-resistant A. fumigatus across agricultural environments, but few have looked at retail plant products. Our objectives were to determine if azole-resistant A. fumigatus is prevalent in retail plant products produced in the United States (U.S.), as well as to identify the resistance mechanism(s) and population genetic structure of these isolates. Five hundred twenty-five isolates were collected from retail plant products and screened for azole resistance. Twenty-four isolates collected from compost, soil, flower bulbs, and raw peanuts were pan-azole resistant. These isolates had the TR34/L98H, TR46/Y121F/T289A, G448S, and H147Y cyp51A alleles, all known to underly pan-azole resistance, as well as WT alleles, suggesting that non-cyp51A mechanisms contribute to pan-azole resistance in these isolates. Minimum spanning networks showed two lineages containing isolates with TR alleles or the F46Y/M172V/E427K allele, and discriminant analysis of principle components identified three primary clusters. This is consistent with previous studies detecting three clades of A. fumigatus and identifying pan-azole-resistant isolates with TR alleles in a single clade. We found pan-azole resistance in U.S. retail plant products, particularly compost and flower bulbs, which indicates a risk of exposure to these products for susceptible populations and that highly resistant isolates are likely distributed worldwide on these products.IMPORTANCEAspergillus fumigatus has recently been designated as a critical fungal pathogen by the World Health Organization. It is most deadly to people with compromised immune systems, and with the emergence of antifungal resistance to multiple azole drugs, this disease carries a nearly 100% fatality rate without treatment or if isolates are resistant to the drugs used to treat the disease. It is important to determine the relatedness and origins of resistant A. fumigatus isolates in the environment, including plant-based retail products, so that factors promoting the development and propagation of resistant isolates can be identified.

The Narrow Footprint of Ancient Balancing Selection Revealed by Heterokaryon Incompatibility Genes in Aspergillus fumigatus

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2,142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters

BACKGROUND

Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.

RESULTS

Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.

CONCLUSION

In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Pan-azole- and multi-fungicide-resistant Aspergillus fumigatus is widespread in the United States

Aspergillus fumigatus is an important global fungal pathogen of humans. Azole drugs are among the most effective treatments for A. fumigatus infection. Azoles are also widely used in agriculture as fungicides against fungal pathogens of crops. Azole-resistant A. fumigatus has been increasing in Europe and Asia for two decades where clinical resistance is thought to be driven by agricultural use of azole fungicides. The most prevalent mechanisms of azole resistance in A. fumigatus are tandem repeats (TR) in the cyp51A promoter coupled with mutations in the coding region which result in resistance to multiple azole drugs (pan-azole resistance). Azole-resistant A. fumigatus has been isolated from patients in the United States (U.S.), but little is known about its environmental distribution. To better understand the distribution of azole-resistant A. fumigatus in the U.S., we collected isolates from agricultural sites in eight states and tested 202 isolates for sensitivity to azoles. We found azole-resistant A. fumigatus in agricultural environments in seven states showing that it is widespread in the U.S. We sequenced environmental isolates representing the range of U.S. sample sites and compared them with publicly available environmental worldwide isolates in phylogenetic, principal component, and ADMIXTURE analyses. We found worldwide isolates fell into three clades, and TR-based pan-azole resistance was largely in a single clade that was strongly associated with resistance to multiple agricultural fungicides. We also found high levels of gene flow indicating recombination between clades highlighting the potential for azole-resistance to continue spreading in the U.S.IMPORTANCEAspergillus fumigatus is a fungal pathogen of humans that causes over 250,000 invasive infections each year. It is found in soils, plant debris, and compost. Azoles are the first line of defense antifungal drugs against A. fumigatus. Azoles are also used as agricultural fungicides to combat other fungi that attack plants. Azole-resistant A. fumigatus has been a problem in Europe and Asia for 20 years and has recently been reported in patients in the United States (U.S.). Until this study, we did not know much about azole-resistant A. fumigatus in agricultural settings in the U.S. In this study, we isolated azole-resistant A. fumigatus from multiple states and compared it to isolates from around the world. We show that A. fumigatus which is resistant to azoles and to other strictly agricultural fungicides is widespread in the U.S.

An endogenous DNA virus in an amphibian-killing fungus associated with pathogen genotype and virulence

The global panzootic lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis (Bd) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single-stranded (CRESS)-like DNA virus, which we named Bd DNA virus 1 (BdDV-1). We further sequenced and constructed genome assemblies with long read sequences to find that the virus is integrated into the nuclear genome in some strains. Attempts to cure virus-positive isolates were unsuccessful; however, phenotypic differences between naturally virus-positive and virus-negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo. BdDV-1 is the first-described CRESS DNA mycovirus of zoosporic true fungi, with a distribution inversely associated with the emergence of the panzootic lineage.

Genomic diversity of the pathogenic fungus Aspergillus fumigatus in Japan reveals the complex genomic basis of azole resistance

Aspergillus fumigatus is a pathogenic fungus with a global distribution. The emergence of azole-resistant A. fumigatus (ARAf) other than the TR-mutants is a problem in Japan. Additionally, the genetic diversity of A. fumigatus strains in Japan remains relatively unknown. Here we show the diversity in the A. fumigatus strains isolated in Japan as well as the complexity in the global distribution of the pathogenic strains. First, we analyzed the genome sequences of 171 strains from Japan as well as the antifungal susceptibility of these strains. Next, we conducted a population analysis of 876 strains by combining the available genomic data for strains isolated worldwide, which were grouped in six clusters. Finally, a genome-wide association study identified the genomic loci associated with ARAf strains, but not the TR-mutants. These results highlight the complexity of the genomic mechanism underlying the emergence of ARAf strains other than the TR-mutants.

Strain heterogeneity in a non-pathogenic fungus highlights factors contributing to virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolomic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Evaluation of Different Gene Prediction Tools in Coccidioides immitis

Gene prediction is required to obtain optimal biologically meaningful information from genomic sequences, but automated gene prediction software is imperfect. In this study, we compare the original annotation of the Coccidioides immitis RS genome (the reference strain of C. immitis) to annotations using the Funannotate and Augustus genome prediction pipelines. A total of 25% of the originally predicted genes (denoted CIMG) were not found in either the Funannotate or Augustus predictions. A comparison of Funannotate and Augustus predictions also found overlapping but not identical sets of genes. The predicted genes found only in the original annotation (referred to as CIMG-unique) were less likely to have a meaningful functional annotation and a lower number of orthologs and homologs in other fungi than all CIMG genes predicted by the original annotation. The CIMG-unique genes were also more likely to be lineage-specific and poorly expressed. In addition, the CIMG-unique genes were found in clusters and tended to be more frequently associated with transposable elements than all CIMG-predicted genes. The CIMG-unique genes were more likely to have experimentally determined transcription start sites that were further away from the originally predicted transcription start sites, and experimentally determined initial transcription was less likely to result in stable CIMG-unique transcripts. A sample of CIMG-unique genes that were relatively well expressed and differentially expressed in mycelia and spherules was inspected in a genome browser, and the structure of only about half of them was found to be supported by RNA-seq data. These data suggest that some of the CIMG-unique genes are not authentic gene predictions. Genes that were predicted only by the Funannotate pipeline were also less likely to have a meaningful functional annotation, be shorter, and express less well than all the genes predicted by Funannotate. C. immitis genes predicted by more than one annotation are more likely to have predicted functions, many orthologs and homologs, and be well expressed. Lineage-specific genes are relatively uncommon in this group. These data emphasize the importance and limitations of gene prediction software and suggest that improvements to the annotation of the C. immitis genome should be considered.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

Mitogenome Variations in a Global Population of Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous, critical priority human fungal pathogen. Despite its clinical importance, there is limited knowledge regarding the variations of the genome within mitochondria, the powerhouse organelle within eukaryotic cells. In this study, we leveraged publicly available, raw, whole genome sequence data isolates from 1939 to investigate the variations in the mitochondrial genomes of A. fumigatus. These isolates were isolated from 22 countries on six continents, as well as from outer space and from within the International Space Station. In total, our analysis revealed 39 mitochondrial single nucleotide polymorphisms (mtSNPs) within this global sample, and, together, these 39 mtSNPs grouped the 1939 isolates into 79 mitochondrial multilocus genotypes (MLGs). Among the 79 MLGs, 39 were each distributed in at least two countries and 30 were each shared by at least two continents. The two most frequent MLGs were also broadly distributed: MLG11 represented 420 isolates from 11 countries and four continents and while MLG79 represented 418 isolates from 18 countries and five continents, consistent with long-distance dispersals of mitogenomes. Our population genetic analyses of the mtSNPs revealed limited differentiation among continental populations, but highly variable genetic differences among national populations, largely due to localized clonal expansions of different MLGs. Phylogenetic analysis and Discriminant Analysis of Principal Components of mtSNPs suggested the presence of at least three mitogenome clusters. Linkage disequilibrium, Index of Association, and phylogenetic incompatibility analyses collectively suggested evidence for mitogenome recombination in natural populations of A. fumigatus. In addition, sequence read depth analyses revealed an average ratio of ~20 mitogenomes per nuclear genome in this global population, but the ratios varied among strains within and between certain geographic populations. Together, our results suggest evidence for organelle dynamics, genetic differentiation, recombination, and both widespread and localized clonal expansion of the mitogenomes in the global A. fumigatus population.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

OMICS and Other Advanced Technologies in Mycological Applications

Fungi play many roles in different ecosystems. The precise identification of fungi is important in different aspects. Historically, they were identified based on morphological characteristics, but technological advancements such as polymerase chain reaction (PCR) and DNA sequencing now enable more accurate identification and taxonomy, and higher-level classifications. However, some species, referred to as "dark taxa", lack distinct physical features that makes their identification challenging. High-throughput sequencing and metagenomics of environmental samples provide a solution to identifying new lineages of fungi. This paper discusses different approaches to taxonomy, including PCR amplification and sequencing of rDNA, multi-loci phylogenetic analyses, and the importance of various omics (large-scale molecular) techniques for understanding fungal applications. The use of proteomics, transcriptomics, metatranscriptomics, metabolomics, and interactomics provides a comprehensive understanding of fungi. These advanced technologies are critical for expanding the knowledge of the Kingdom of Fungi, including its impact on food safety and security, edible mushrooms foodomics, fungal secondary metabolites, mycotoxin-producing fungi, and biomedical and therapeutic applications, including antifungal drugs and drug resistance, and fungal omics data for novel drug development. The paper also highlights the importance of exploring fungi from extreme environments and understudied areas to identify novel lineages in the fungal dark taxa.

COVID-19-Associated Pulmonary Aspergillosis Isolates Are Genomically Diverse but Similar to Each Other in Their Responses to Infection-Relevant Stresses

Secondary infections caused by the pulmonary fungal pathogen Aspergillus fumigatus are a significant cause of mortality in patients with severe coronavirus disease 19 (COVID-19). Even though epithelial cell damage and aberrant cytokine responses have been linked to susceptibility to COVID-19-associated pulmonary aspergillosis (CAPA), little is known about the mechanisms underpinning copathogenicity. Here, we analyzed the genomes of 11 A. fumigatus isolates from patients with CAPA in three centers from different European countries. CAPA isolates did not cluster based on geographic origin in a genome-scale phylogeny of representative A. fumigatus isolates. Phenotypically, CAPA isolates were more similar to the A. fumigatus A1160 reference strain than to the Af293 strain when grown in infection-relevant stresses, except for interactions with human immune cells wherein macrophage responses were similar to those induced by the Af293 reference strain. Collectively, our data indicate that CAPA isolates are genomically diverse but are more similar to each other in their responses to infection-relevant stresses. A larger number of isolates from CAPA patients should be studied to better understand the molecular epidemiology of CAPA and to identify genetic drivers of copathogenicity and antifungal resistance in patients with COVID-19. IMPORTANCE Coronavirus disease 2019 (COVID-19)-associated pulmonary aspergillosis (CAPA) has been globally reported as a life-threatening complication in some patients with severe COVID-19. Most of these infections are caused by the environmental mold Aspergillus fumigatus, which ranks third in the fungal pathogen priority list of the WHO. However, little is known about the molecular epidemiology of Aspergillus fumigatus CAPA strains. Here, we analyzed the genomes of 11 A. fumigatus isolates from patients with CAPA in three centers from different European countries, and carried out phenotypic analyses with a view to understanding the pathophysiology of the disease. Our data indicate that A. fumigatus CAPA isolates are genomically diverse but are more similar to each other in their responses to infection-relevant stresses.

High-Throughput Microsatellite Markers Development for Genetic Characterization of Emerging Sporothrix Species

Sporotrichosis is the main subcutaneous mycosis worldwide transmitted by animal or plant vectors and often escalates to outbreaks or epidemics. The current cat-transmitted sporotrichosis driven by Sporothrix brasiliensis has become a significant public health issue in South America. Transmission dynamics remain enigmatic due to the lack of development of polymorphic markers for molecular epidemiological analysis. This study used a high-throughput mining strategy to characterize simple sequence repeat (SSR) markers from Sporothrix genomes. A total of 118,140-143,912 SSR loci were identified (82,841-98,369 unique markers), with a 3651.55-3804.65 SSR/Mb density and a majority of dinucleotides motifs (GC/CG). We developed a panel of 15 highly polymorphic SSR markers suitable for genotyping S. brasiliensis, S. schenckii, and S. globosa. PCR amplification revealed 240 alleles in 180 Sporothrix isolates with excellent polymorphic information content (PIC = 0.9101), expected heterozygosity (H = 0.9159), and discriminating power (D = 0.7127), supporting the effectiveness of SSR markers in uncovering cryptic genetic diversity. A systematic population genetic study estimated three clusters, corresponding to S. brasiliensis (population 1, n = 97), S. schenckii (population 2, n = 49), and S. globosa (population 3, n = 34), with a weak signature of mixed ancestry between populations 1 and 2 or 3 and 2. Partitioning of genetic variation via AMOVA revealed highly structured populations (ΦPT = 0.539; Nm = 0.213; p < 0.0001), with approximately equivalent genetic variability within (46%) and between (54%) populations. Analysis of SSR diversity supports Rio de Janeiro (RJ) as the center of origin for contemporary S. brasiliensis infections. The recent emergence of cat-transmitted sporotrichosis in northeastern Brazil indicates an RJ-Northeast migration resulting in founder effects during the introduction of diseased animals into sporotrichosis-free areas. Our results demonstrated high cross-species transferability, reproducibility, and informativeness of SSR genetic markers, helping dissect deep and fine-scale genetic structures and guiding decision making to mitigate the harmful effects of the expansion of cat-transmitted sporotrichosis.

Insight into the Organization of the B10v3 Cucumber Genome by Integration of Biological and Bioinformatic Data

The availability of a well-organized and annotated reference genome is essential for genome research and the analysis of re-sequencing approaches. The B10v3 cucumber (Cucumis sativus L.) reference genome has been sequenced and assembled into 8035 contigs, a small fraction of which have been mapped to individual chromosomes. Currently, bioinformatics methods based on comparative homology have made it possible to re-order the sequenced contigs by mapping them to the reference genomes. The B10v3 genome (North-European, Borszczagowski line) was rearranged against the genomes of cucumber 9930 ('Chinese Long' line) and Gy14 (North American line). Furthermore, a better insight into the organization of the B10v3 genome was obtained by integrating the data available in the literature on the assignment of contigs to chromosomes in the B10v3 genome with the results of the bioinformatic analysis. The combination of information on the markers used in the assembly of the B10v3 genome and the results of FISH and DArT-seq experiments confirmed the reliability of the in silico assignment. Approximately 98% of the protein-coding genes within the chromosomes were assigned and a significant proportion of the repetitive fragments in the sequenced B10v3 genome were identified using the RagTag programme. In addition, BLAST analyses provided comparative information between the B10v3 genome and the 9930 and Gy14 data sets. This revealed both similarities and differences in the functional proteins found between the coding sequences region in the genomes. This study contributes to better knowledge and understanding of cucumber genome line B10v3.