​Fusarium Protein Toolkit: a web-based resource for structural and variant analysis of Fusarium species.

BACKGROUND: ​​The genus Fusarium poses significant threats to food security and safety worldwide because numerous species of the fungus cause destructive diseases and/or mycotoxin contamination in crops. The adverse effects of climate change are exacerbating some existing threats and causing new problems. These challenges highlight the need for innovative solutions, including the development of advanced tools to identify targets for control strategies. DESCRIPTION: In response to these challenges, we developed the Fusarium Protein Toolkit (FPT), a web-based tool that allows users to interrogate the structural and variant landscape within the Fusarium pan-genome. The tool displays both AlphaFold and ESMFold-generated protein structure models from six Fusarium species. The structures are accessible through a user-friendly web portal and facilitate comparative analysis, functional annotation inference, and identification of related protein structures. Using a protein language model, FPT predicts the impact of over 270 million coding variants in two of the most agriculturally important species, Fusarium graminearum and F. verticillioides. To facilitate the assessment of naturally occurring genetic variation, FPT provides variant effect scores for proteins in a Fusarium pan-genome based on 22 diverse species. The scores indicate potential functional consequences of amino acid substitutions and are displayed as intuitive heatmaps using the PanEffect framework. CONCLUSION: FPT fills a knowledge gap by providing previously unavailable tools to assess structural and missense variation in proteins produced by Fusarium. FPT has the potential to deepen our understanding of pathogenic mechanisms in Fusarium, and aid the identification of genetic targets for control strategies that reduce crop diseases and mycotoxin contamination. Such targets are vital to solving the agricultural problems incited by Fusarium, particularly evolving threats resulting from climate change. Thus, FPT has the potential to contribute to improving food security and safety worldwide.

Complete genome sequence analysis of Pestalotiopsis microspora, a fungal pathogen causing kiwifruit postharvest rots.

BACKGROUND: The postharvest rot of kiwifruit is one of the most devastating diseases affecting kiwifruit quality worldwide. However, the genomic basis and pathogenicity mechanisms of kiwifruit rot pathogens are lacking. Here we report the first whole genome sequence of Pestalotiopsis microspora, one of the main pathogens causing postharvest kiwifruit rot in China. The genome of strain KFRD-2 was sequenced, de novo assembled, and analyzed. RESULTS: The genome of KFRD-2 was estimated to be approximately 50.31 Mb in size, with an overall GC content of 50.25%. Among 14,711 predicted genes, 14,423 (98.04%) exhibited significant matches to genes in the NCBI nr database. A phylogenetic analysis of 26 known pathogenic fungi, including P. microspora KFRD-2, based on conserved orthologous genes, revealed that KFRD-2's closest evolutionary relationships were to Neopestalotiopsis spp. Among KFRD-2's coding genes, 870 putative CAZy genes spanned six classes of CAZys, which play roles in degrading plant cell walls. Out of the 25 other plant pathogenic fungi, P. microspora possessed a greater number of CAZy genes than 22 and was especially enriched in GH and AA genes. A total of 845 transcription factors and 86 secondary metabolism gene clusters were predicted, representing various types. Furthermore, 28 effectors and 109 virulence-enhanced factors were identified using the PHI (pathogen host-interacting) database. CONCLUSION: This complete genome sequence analysis of the kiwifruit postharvest rot pathogen P. microspora enriches our understanding its disease pathogenesis and virulence. This study establishes a theoretical foundation for future investigations into the pathogenic mechanisms of P. microspora and the development of enhanced strategies for the efficient management of kiwifruit postharvest rots.

A putative scenario of how de novo protein-coding genes originate in the Saccharomyces cerevisiae lineage.

BACKGROUND: Novel protein-coding genes were considered to be born by re-organization of pre-existing genes, such as gene duplication and gene fusion. However, recent progress of genome research revealed that more protein-coding genes than expected were born de novo, that is, gene origination by accumulating mutations in non-genic DNA sequences. Nonetheless, the in-depth process (scenario) for de novo origination is not well understood. RESULTS: We have conceived bioinformatic analysis for sketching a scenario for de novo origination of protein-coding genes. For each de novo protein-coding gene, we firstly identified an edge of a given phylogenetic tree where the gene was born based on parsimony. Then, from a multiple sequence alignment of the de novo gene and its orthologous regions, we constructed ancestral DNA sequences of the gene corresponding to both end nodes of the edge. We finally revealed statistical features observed in evolution between the two ancestral sequences. In the analysis of the Saccharomyces cerevisiae lineage, we have successfully sketched a putative scenario for de novo origination of protein-coding genes. (1) In the beginning was GC-rich genome regions. (2) Neutral mutations were accumulated in the regions. (3) ORFs were extended/combined, and then (4) translation signature (Kozak consensus sequence) was recruited. Interestingly, as the scenario progresses from (2) to (4), the specificity of mutations increases. CONCLUSION: To the best of our knowledge, this is the first report outlining a scenario of de novo origination of protein-coding genes. Our bioinformatic analysis can capture events that occur during a short evolutionary time by directly observing the evolution of the ancestral sequences from non-genic to genic. This property is suitable for the analysis of fast evolving de novo genes.

Complete Genome Sequence of Candida mucifera from an Otitis Media Patient.

We describe for the first time, a high-quality genome for a rare human yeast pathogen Candida mucifera, from a patient with chronic suppurative otitis media. This pathogen exhibited reduced azole susceptibility, similar to its close relatives within the Trichomonascus ciferrii species complex.

Two near-chromosomal-level genomes of globally-distributed Macroascomycete based on single-molecule fluorescence and Hi-C methods.

Discinaceae holds significant importance within the Pezizales, representing a prominent group of macroascomycetes distributed globally. However, there is a dearth of genomic studies focusing on this family, resulting in gaps in our understanding of its evolution, development, and ecology. Here we utilized state-of-the-art genome assembly methodologies, incorporating third-generation single-molecule fluorescence and Hi-C-assisted methods, to elucidate the genomic landscapes of Gyromitra esculenta and Paragyromitra xinjiangensis. The genome sizes of two species were determined to be 47.10 Mb and 48.20 Mb, with 23 and 22 scaffolds, respectively. 10,438 and 11,469 coding proteins were identified, with functional annotations encompassing over 96.47% and 94.40%, respectively. Assessment of completeness using BUSCO revealed that 98.71% and 98.89% of the conserved proteins were identified. The application of comparative genomic technology has helped in identifying traits associated with of heterothallic life cycle traits and elucidating unique patterns of chromosomal evolution. Additionally, we identified potential saprotrophic nutritional modes and systematic phylogenetic relationships between the two species. Therefore, this study provides crucial genomic insights into the evolution, nutritional type, and ecological roles of species within the Pezizales.

Strain heterogeneity in a non-pathogenic Aspergillus fungus highlights factors associated with 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 metabolic 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 user-friendly CRISPR/Cas9 system for mutagenesis of Neurospora crassa.

As a widely used eukaryotic model organism, Neurospora crassa offers advantages in genetic studies due to its diverse biology and rapid growth. Traditional genetic manipulation methods, such as homologous recombination, require a considerable amount of time and effort. In this study, we present an easy-to-use CRIPSR/Cas9 system for N. crassa, in which the cas9 sequence is incorporated into the fungal genome and naked guide RNA is introduced via electroporation. Our approach eliminates the need for constructing multiple vectors, speeding up the mutagenesis process. Using cyclosporin-resistant-1 (csr-1) as a selectable marker gene, we achieved 100% editing efficiency under selection conditions. Furthermore, we successfully edited the non-selectable gene N-acylethanolamine amidohydrolase-2 (naa-2), demonstrating the versatility of the system. Combining gRNAs targeting csr-1 and naa-2 simultaneously increased the probability of finding mutants carrying the non-selectable mutation. The system is not only user-friendly but also effective, providing a rapid and efficient method for generating loss-of-function mutants in N. crassa compared to traditional methods.

Genomic characterization and identification of candidate genes for putative podophyllotoxin biosynthesis pathway in Penicillium herquei HGN12.1C.

Previous studies have reported the functional role, biochemical features and synthesis pathway of podophyllotoxin (PTOX) in plants. In this study, we employed combined morphological and molecular techniques to identify an endophytic fungus and extract PTOX derivatives. Based on the analysis of ITS sequences and the phylogenetic tree, the isolate was classified as Penicillium herquei HGN12.1C, with a sequence identity of 98.58%. Morphologically, the HGN12.1C strain exhibits white colonies, short-branched mycelia and densely packed hyphae. Using PacBio sequencing at an average read depth of 195×, we obtained a high-quality genome for the HGN12.1C strain, which is 34.9 Mb in size, containing eight chromosomes, one mitochondrial genome and a GC content of 46.5%. Genome analysis revealed 10 genes potentially involved in PTOX biosynthesis. These genes include VdtD, Pinoresinollariciresinol reductase (PLR), Secoisolariciresinol dehydrogenase (SDH), CYP719A23, CYP71BE54, O-methyltransferase 1 (OMT1), O-methyltransferase 3 (OMT3), 2-ODD, CYP71CU and CYP82D61. Notably, the VdtD gene in fungi shares functional similarities with the DIR gene found in plants. Additionally, we identified peltatin, a PTOX derivative, in the HGN12.1C extract. Docking analysis suggests a potential role for the 2-ODD enzyme in converting yatein to deoxypodophyllotoxin. These findings offer invaluable insights into the synthesis mechanism of PTOX in fungi, shedding light on the relationship between host plants and endophytes.

Evolution and related pathogenic genes of Pseudodiploöspora longispora on Morchella based on genomic characterization and comparative genomic analysis.

True morels (Morchella) are globally renowned medicinal and edible mushrooms. White mold disease caused by fungi is the main disease of Morchella, which has the characteristics of wide incidence and strong destructiveness. The disparities observed in the isolation rates of different pathogens indicate their varying degrees of host adaptability and competitive survival abilities. In order to elucidate its potential mechanism, this study, the pathogen of white mold disease from Dafang county, Guizhou Province was isolated and purified, identified as Pseudodiploöspora longispora by morphological, molecular biological and pathogenicity tests. Furthermore, high-quality genome of P. longisporus (40.846 Mb) was assembled N50 of 3.09 Mb, predicts 7381 protein-coding genes. Phylogenetic analysis of single-copy homologous genes showed that P. longispora and Zelopaecilomyces penicillatus have the closest evolutionary relationship, diverging into two branches approximately 50 (44.3-61.4) MYA. Additionally, compared with the other two pathogens causing Morchella disease, Z. penicillatus and Cladobotryum protrusum, it was found that they had similar proportions of carbohydrate enzyme types and encoded abundant cell wall degrading enzymes, such as chitinase and glucanase, indicating their important role in disease development. Moreover, the secondary metabolite gene clusters of P. longispora and Z. penicillatus show a high degree of similarity to leucinostatin A and leucinostatin B (peptaibols). Furthermore, a gene cluster with synthetic toxic substance Ochratoxin A was also identified in P. longispora and C. protrusum, indicating that they may pose a potential threat to food safety. This study provides valuable insights into the genome of P. longispora, contributing to pathogenicity research.

ChEC-Seq: A Comprehensive Guide for Scalable and Cost-Efficient Genome-Wide Profiling in Saccharomyces cerevisiae.

Chromatin endogenous cleavage coupled with high-throughput sequencing (ChEC-seq) is a profiling method for protein-DNA interactions that can detect binding locations in vivo, does not require antibodies or fixation, and provides genome-wide coverage at near nucleotide resolution.The core of this method is an MNase fusion of the target protein, which allows it, when triggered by calcium exposure, to cut DNA at its binding sites and to generate small DNA fragments that can be readily separated from the rest of the genome and sequenced.Improvements since the original protocol have increased the ease, lowered the costs, and multiplied the throughput of this method to enable a scale and resolution of experiments not available with traditional methods such as ChIP-seq. This method describes each step from the initial creation and verification of the MNase-tagged yeast strains, over the ChEC MNase activation and small fragment purification procedure to the sequencing library preparation. It also briefly touches on the bioinformatic steps necessary to create meaningful genome-wide binding profiles.

Utilizing Deep Neural Networks to Fill Gaps in Small Genomes.

With the widespread adoption of next-generation sequencing technologies, the speed and convenience of genome sequencing have significantly improved, and many biological genomes have been sequenced. However, during the assembly of small genomes, we still face a series of challenges, including repetitive fragments, inverted repeats, low sequencing coverage, and the limitations of sequencing technologies. These challenges lead to unknown gaps in small genomes, hindering complete genome assembly. Although there are many existing assembly software options, they do not fully utilize the potential of artificial intelligence technologies, resulting in limited improvement in gap filling. Here, we propose a novel method, DLGapCloser, based on deep learning, aimed at assisting traditional tools in further filling gaps in small genomes. Firstly, we created four datasets based on the original genomes of Saccharomyces cerevisiae, Schizosaccharomyces pombe, Neurospora crassa, and Micromonas pusilla. To further extract effective information from the gene sequences, we also added homologous genomes to enrich the datasets. Secondly, we proposed the DGCNet model, which effectively extracts features and learns context from sequences flanking gaps. Addressing issues with early pruning and high memory usage in the Beam Search algorithm, we developed a new prediction algorithm, Wave-Beam Search. This algorithm alternates between expansion and contraction phases, enhancing efficiency and accuracy. Experimental results showed that the Wave-Beam Search algorithm improved the gap-filling performance of assembly tools by 7.35%, 28.57%, 42.85%, and 8.33% on the original results. Finally, we established new gap-filling standards and created and implemented a novel evaluation method. Validation on the genomes of Saccharomyces cerevisiae, Schizosaccharomyces pombe, Neurospora crassa, and Micromonas pusilla showed that DLGapCloser increased the number of filled gaps by 8.05%, 15.3%, 1.4%, and 7% compared to traditional assembly tools.

Characteristics of Corynespora cassiicola, the causal agent of tobacco Corynespora leaf spot, revealed by genomic and metabolic phenomic analysis.

Corynespora cassiicola is a highly diverse fungal pathogen that commonly occurs in tropical, subtropical, and greenhouse environments worldwide. In this study, the isolates were identified as C. cassiicola, and the optimum growth and sporulation were studied. The phenotypic characteristics of C. cassiicola, concerning 950 different growth conditions, were tested using Biolog PM plates 1-10. In addition, the strain of C. cassiicola DWZ from tobacco hosts was sequenced for the using Illumina PE150 and Pacbio technologies. The host resistance of tobacco Yunyan 87 with different maturity levels was investigated. In addition, the resistance evaluation of 10 common tobacco varieties was investigated. The results showed that C. cassiicola metabolized 89.47% of the tested carbon source, 100% of the nitrogen source, 100% of the phosphorus source, and 97.14% of the sulfur source. It can adapt to a variety of different osmotic pressure and pH environments, and has good decarboxylase and deaminase activities. The optimum conditions for pathogen growth and sporulation were 25-30 °C, and the growth was better on AEA and OA medium. The total length of the genome was 45.9 Mbp, the GC content was 51.23%, and a total of 13,061 protein-coding genes, 202 non-coding RNAs and 2801 and repeat sequences were predicted. Mature leaves were more susceptible than proper mature and immature leaves, and the average diameter of diseased spots reached 17.74 mm at 12 days. None of the tested ten cultivars exhibited obvious resistance to Corynespora leaf spot of tobacco, whereby all disease spot diameters reached > 10 mm and > 30 mm when at 5 and 10 days after inoculation, respectively. The phenotypic characteristics, genomic analysis of C. cassiicola and the cultivar resistance assessment of this pathogen have increased our understanding of Corynespora leaf spot of tobacco.

Datasets of fungal diversity and pseudo-chromosomal genomes of mangrove rhizosphere soil in China.

With climate change and anthropic influence on the coastal ecosystems, mangrove ecosystems are disappearing at an alarming rate. Accordingly, it becomes important to track, study, record and store the mangrove microbial community considering their ecological importance and potential for biotechnological applications. Here, we provide information on mangrove fungal community composition and diversity in mangrove ecosystems with different plant species and from various locations differing in relation to anthropic influences. We describe twelve newly assembled genomes, including four chromosomal-level genomes of fungal isolates from the mangrove ecosystems coupled with functional annotations. We envisage that these data will be of value for future studies including comparative genome analysis and large-scale temporal and/or spatial research to elucidate the potential mechanisms by which mangrove fungal communities assemble and evolve. We further anticipate that the genomes represent valuable resources for bioprospecting related to industrial or clinical uses.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome.

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Hi-C2B: Optimized Detection of Chromosomal Contacts Within Synchronized Meiotic S. cerevisiae Cells.

Hi-C, a genome-wide chromosome conformation capture assay, is a powerful tool used to study three-dimensional genome organization by converting physical pairwise interactions into counts of pairwise interactions. To study the many temporally regulated facets of meiotic recombination in S. cerevisiae, the Hi-C assay must be robust such that fine- and wide-scale comparisons between genetic datasets can be made. Here we describe an updated protocol for Hi-C (Hi-C2B) that generates reproducible libraries of interaction data with low noise and for a relatively low cost.

Genome sequencing and CAZymes repertoire analysis of Diaporthe eres P3-1W causing postharvest fruit rot of 'Hongyang' kiwifruit in China.

Postharvest rot caused by various fungal pathogens is a damaging disease affecting kiwifruit production and quality, resulting in significant annual economic losses. This study focused on isolating the strain P3-1W, identified as Diaporthe eres, as the causal agent of 'Hongyang' postharvest rot disease in China. The investigation highlighted cell wall degrading enzymes (CWDEs) as crucial pathogenic factors. Specially, the enzymatic activities of cellulase, ß-galactosidase, polygalacturonase, and pectin methylesterases peaked significantly on the second day after infection of D. eres P3-1W. To gain a comprehensive understanding of these CWDEs, the genome of this strain was sequenced using PacBio and Illumina sequencing technologies. The analysis revealed that the genome of D. eres P3-1W spans 58,489,835 bp, with an N50 of 5,939,879 bp and a GC content of 50.7%. A total of 15,407 total protein-coding genes (PCGs) were predicted and functionally annotated. Notably, 857 carbohydrate-active enzymes (CAZymes) were identified in D. eres P3-1W, with 521 CWDEs consisting of 374 glycoside hydrolases (GHs), 108 carbohydrate esterase (CEs) and 91 polysaccharide lyases (PLs). Additionally, 221 auxiliary activities (AAs), 91 glycosyltransferases (GTs), and 108 carbohydrate binding modules (CBMs) were detected. These findings offer valuable insights into the CAZymes of D. eres P3-1W.

Fragment-Guided Genome Mining of Octacyclic Cyclophane Alkaloids from Fungi.

Nature uses compact but functionalized biosynthetic fragments as building blocks to generate complex natural products. To leverage this strategy for the discovery of natural products with new scaffolds, we performed genome mining to identify biosynthetic gene clusters (BGCs) in fungi that embed genes that can synthesize targeted fragments. The three-enzyme pathway that biosynthesizes the strained dityrosine cyclophane in the herquline A pathway was used to identify a large number of potential BGCs that may use the cyclophane as a fragment. Characterization of a conserved BGC from fungal strains led to the isolation of octacyclin A, an octacyclic natural product with an unprecedented structure, including two hetero-[3.3.1]bicycles and a combination of fused, bridged, and macrocyclic rings. Biosynthetic steps leading to octacyclin A were fully elucidated using pathway reconstitution and enzymatic assays, unveiling intriguing chemical logic and new enzymatic reactions in building the octacyclic core. Our work demonstrates the potential utility of fragment-guided genome mining in expanding natural product chemical space.

Draft genome sequence and annotation of the polyextremotolerant polyol lipid-producing fungus aureobasidium pullulans NRRL 62042.

OBJECTIVES: The ascomycotic yeast-like fungus Aureobasidium exhibits the natural ability to synthesize several secondary metabolites, like polymalic acid, pullulan, or polyol lipids, with potential biotechnological applications. Combined with its polyextremotolerance, these properties make Aureobasidium a promising production host candidate. Hence, plenty of genomes of Aureobasidia have been sequenced recently. Here, we provide the annotated draft genome sequence of the polyol lipid-producing strain A. pullulans NRRL 62042. DATA DESCRIPTION: The genome of A. pullulans NRRL 62042 was sequenced using Illumina NovaSeq 6000. Genome assembly revealed a genome size of 24.2 Mb divided into 39 scaffolds with a GC content of 50.1%. Genome annotation using Genemark v4.68 and GenDBE yielded 9,596 genes.

MCM2-7 loading-dependent ORC release ensures genome-wide origin licensing.

Origin recognition complex (ORC)-dependent loading of the replicative helicase MCM2-7 onto replication origins in G1-phase forms the basis of replication fork establishment in S-phase. However, how ORC and MCM2-7 facilitate genome-wide DNA licensing is not fully understood. Mapping the molecular footprints of budding yeast ORC and MCM2-7 genome-wide, we discovered that MCM2-7 loading is associated with ORC release from origins and redistribution to non-origin sites. Our bioinformatic analysis revealed that origins are compact units, where a single MCM2-7 double hexamer blocks repetitive loading through steric ORC binding site occlusion. Analyses of A-elements and an improved B2-element consensus motif uncovered that DNA shape, DNA flexibility, and the correct, face-to-face spacing of the two DNA elements are hallmarks of ORC-binding and efficient helicase loading sites. Thus, our work identified fundamental principles for MCM2-7 helicase loading that explain how origin licensing is realised across the genome.

A gap-free genome assembly of Fusarium oxysporum f. sp. conglutinans, a vascular wilt pathogen.

Fusarium oxysporum is an asexual filamentous fungus that causes vascular wilt in hundreds of crop plants and poses a threat to public health through Fusariosis. F. oxysporum f. sp. conglutinans strain Fo5176, originally isolated from Brassica oleracea, is pathogenic to Arabidopsis, making it a model pathosystem for dissecting the molecular mechanisms underlying host-pathogen interactions. Assembling the F. oxysporum genome is notoriously challenging due to the presence of repeat-rich accessory chromosomes. Here, we report a gap-free genome assembly of Fo5176 using PacBio HiFi and Hi-C data. The 69.56 Mb assembly contained 18 complete chromosomes, including all centromeres and most telomeres (20/36), representing the first gap-free genome sequence of a pathogenic F. oxysporum strain. In total, 21,460 protein-coding genes were annotated, a 26.3% increase compared to the most recent assembly. This high-quality reference genome for F. oxysporum f. sp. conglutinans Fo5176 provides a valuable resource for further research into fungal pathobiology and evolution.