New protein-coding genes can evolve from previously noncoding genomic regions through a process known as de novo gene emergence. Evidence suggests that this process has likely occurred throughout evolution and across the tree of life. Yet, confidently identifying de novo emerged genes remains challenging. Ancestral sequence reconstruction is a promising approach for inferring whether a gene has emerged de novo or not, as it allows us to inspect whether a given genomic locus ancestrally harbored protein-coding capacity. However, the use of ancestral sequence reconstruction in the context of de novo emergence is still in its infancy and its capabilities, limitations, and overall potential are largely unknown. Notably, it is difficult to formally evaluate the protein-coding capacity of ancestral sequences, particularly when new gene candidates are short. How well-suited is ancestral sequence reconstruction as a tool for the detection and study of de novo genes? Here, we address this question by designing an ancestral sequence reconstruction workflow incorporating different tools and sets of parameters and by introducing a formal criterion that allows to estimate, within a desired level of confidence, when protein-coding capacity originated at a particular locus. Applying this workflow on ∼2,600 short, annotated budding yeast genes (<1,000 nucleotides), we found that ancestral sequence reconstruction robustly predicts an ancient origin for the most widely conserved genes, which constitute \"easy\" cases. For less robust cases, we calculated a randomization-based empirical P-value estimating whether the observed conservation between the extant and ancestral reading frame could be attributed to chance. This formal criterion allowed us to pinpoint a branch of origin for most of the less robust cases, identifying 49 genes that can unequivocally be considered de novo originated since the split of the Saccharomyces genus, including 37 Saccharomyces cerevisiae-specific genes. We find that for the remaining equivocal cases we cannot rule out different evolutionary scenarios including rapid evolution, multiple gene losses, or a recent de novo origin. Overall, our findings suggest that ancestral sequence reconstruction is a valuable tool to study de novo gene emergence but should be applied with caution and awareness of its limitations.

Utilizing bioinformatics tools, this study expands our understanding of secondary metabolism in Botrytis cinerea, identifying novel genes within polyketide synthase (PKS), non-ribosomal peptide synthetase (NRPS), sesquiterpene cyclase (STC), diterpene cyclase (DTC), and dimethylallyltryptophan synthase (DMATS) families. These findings enrich the genetic framework associated with B. cinerea\'s pathogenicity and ecological adaptation, offering insights into uncharted metabolic pathways. Significantly, the discovery of previously unannotated genes provides new molecular targets for developing targeted antifungal strategies, promising to enhance crop protection and advance our understanding of fungal biochemistry. This research not only broadens the scope of known secondary metabolites but also opens avenues for future exploration into B. cinerea\'s biosynthetic capabilities, potentially leading to novel antifungal compounds. Our work underscores the importance of integrating bioinformatics and genomics for fungal research, paving the way for sustainable agricultural practices by pinpointing precise molecular interventions against B. cinerea. This study sets a foundation for further investigations into the fungus\'s secondary metabolism, with implications for biotechnology and crop disease management.

Reverse genetic approaches are common tools in genomics for elucidating gene functions, involving techniques such as gene deletion followed by screening for aberrant phenotypes. If the generation of gene deletion mutants fails, the question arises whether the failure stems from technical issues or because the gene of interest (GOI) is essential, meaning that the deletion causes lethality. In this report, we introduce a novel method for assessing gene essentiality using the phytopathogenic ascomycete Magnaporthe oryzae. The method is based on the observation that telomere vectors are lost in transformants during cultivation without selection pressure. We tested the hypothesis that essential genes can be identified in deletion mutants co-transformed with a telomere vector. The M. oryzae gene MoPKC, described in literature as essential, was chosen as GOI. Using CRISPR/Cas9 technology transformants with deleted GOI were generated and backed up by a telomere vector carrying a copy of the GOI and conferring fenhexamid resistance. Transformants in which the GOI deletion in the genome was not successful lost the telomere vector on media without fenhexamid. In contrast, transformants with confirmed GOI deletion retained the telomere vector even in absence of fenhexamid selection. In the latter case, the maintenance of the telomere indicates that the GOI is essential for the surveillance of the fungi, as it would have been lost otherwise. The method presented here allows to test for essentiality of genes when no mutants can be obtained from gene deletion approaches, thereby expanding the toolbox for studying gene function in ascomycetes.

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.
CONCLUSIONS: 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.

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.

Fungal secondary metabolites are not necessary for growth, but they are important for fungal metabolism and ecology because they provide selective advantages for competition, survival and interactions with the environment. These various metabolites are widely used as medicinal precursors and insecticides. Secondary metabolism genes are commonly arranged in clusters along chromosomes, which allow for the coordinate control of complete pathways. In this study, we created the Fungal Gene Cluster Database to store, retrieve, and visualize secondary metabolite gene cluster information across fungal species. The database was created by merging data from RNA sequencing, Basic Local Alignment Search Tool, genome browser, enrichment analysis and the R Shiny web framework to visualize and query putative gene clusters. This database facilitated the rapid and thorough examination of significant gene clusters across fungal species by detecting, defining and graphically displaying the architecture, organization and expression patterns of secondary metabolite gene clusters. In general, this genomic resource makes use of the tremendous chemical variety of the products of these ecologically and biotechnologically significant gene clusters to our further understanding of fungal secondary metabolism. Database URL: https://www.hebaubioinformatics.cn/FungalGeneCluster/.

Siderophores are crucial for iron-scavenging in microorganisms. While many yeasts can uptake siderophores produced by other organisms, they are typically unable to synthesize siderophores themselves. In contrast, Wickerhamiella/Starmerella (W/S) clade yeasts gained the capacity to make the siderophore enterobactin following the remarkable horizontal acquisition of a bacterial operon enabling enterobactin synthesis. Yet, how these yeasts absorb the iron bound by enterobactin remains unresolved. Here, we demonstrate that Enb1 is the key enterobactin importer in the W/S-clade species Starmerella bombicola. Through phylogenomic analyses, we show that ENB1 is present in all W/S clade yeast species that retained the enterobactin biosynthetic genes. Conversely, it is absent in species that lost the ent genes, except for Starmerella stellata, making this species the only cheater in the W/S clade that can utilize enterobactin without producing it. Through phylogenetic analyses, we infer that ENB1 is a fungal gene that likely existed in the W/S clade prior to the acquisition of the ent genes and subsequently experienced multiple gene losses and duplications. Through phylogenetic topology tests, we show that ENB1 likely underwent horizontal gene transfer from an ancient W/S clade yeast to the order Saccharomycetales, which includes the model yeast Saccharomyces cerevisiae, followed by extensive secondary losses. Taken together, these results suggest that the fungal ENB1 and bacterial ent genes were cooperatively integrated into a functional unit within the W/S clade that enabled adaptation to iron-limited environments. This integrated fungal-bacterial circuit and its dynamic evolution determine the extant distribution of yeast enterobactin producers and cheaters.

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 of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. 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 genes Avr4E, Avr5, and Avr9, 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 of C. 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.

Groundnuts are mostly contaminated with the mold Aspergillus flavus which produces a carcinogenic mycotoxin called as aflatoxin. It is very important to understand the genetic factors underlying its pathogenicity, regulation, and biosynthesis of secondary metabolites and animal toxicities, but it still lacks useful information due to certain gaps in the era of modern technology. Therefore, the present study was considered to determine the key genes and metabolites involved in the biosynthesis of aflatoxin by using a molecular approach in a virulent strain of Aspergillus. The whole genome sequence of highly toxic and virulent Aspergillus isolates JAM-JKB-B HA-GG20 revealed 3,73,54,834 bp genome size, 2, 26, 257 number of contigs with N50 value of 49,272 bp, 12,400 genes and 48.1% of GC contained respectively. The genome sequence was compared with other known aflatoxin producing and non-producing genome of Aspergillus spp. and 61 secondary metabolite (SM) gene clusters were annotated with the toxic strain JAM-JKB-BHA-GG20 which showed similarity with other Aspergillus spp. A total number of eight genes (ver-1, AflR, pksA, uvm8, omt1, nor-1, Vha and aflP) were identified related to biosynthesis of aflatoxin and ochratoxin. Also, 69 SSR with forward and reverse primers and 137 di and tri nucleotide motifs were identified in the nucleotide sequence region related to aflatoxin gene pathway. The genes and putative metabolites identified in this study are potentially involved in host invasion and pathogenicity. As such, the genomic information obtained in this study is helpful in understanding aflatoxin gene producing pathway in comparison to other Aspergillus spp. and predicted presence of other secondary metabolites clusters viz. Nrps, T1pks etc. genes associated with a biosynthesis of OTA mycotoxin.

BACKGROUND: Penicillium chrysogenum is a filamentous fungal species with diverse habitats, yet little is known about its genetics in adapting to extreme subseafloor sedimental environments.
RESULTS: Here, we report the discovery of P. chrysogenum strain 28R-6-F01, isolated from deep coal-bearing sediments 2306 m beneath the seafloor. This strain possesses exceptional characteristics, including the ability to thrive in extreme conditions such as high temperature (45 °C), high pressure (35 Mpa), and anaerobic environments, and exhibits broad-spectrum antimicrobial activity, producing the antibiotic penicillin at a concentration of 358 μg/mL. Genome sequencing and assembly revealed a genome size of 33.19 Mb with a GC content of 48.84%, containing 6959 coding genes. Comparative analysis with eight terrestrial strains identified 88 unique genes primarily associated with penicillin and aflatoxins biosynthesis, carbohydrate degradation, viral resistance, and three secondary metabolism gene clusters. Furthermore, significant expansions in gene families related to DNA repair were observed, likely linked to the strain\'s adaptation to its environmental niche.
CONCLUSIONS: Our findings provide insights into the genomic and biological characteristics of P. chrysogenum adaptation to extreme anaerobic subseafloor sedimentary environments, such as high temperature and pressure.