Development of a telomere vector-based approach to overcome limitations caused by lethal phenotypes in the study of essential genes in Magnaporthe oryzae.

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.

Identification of two genes essential for basidiospore formation during the postmeiotic stages in Pleurotus ostreatus.

A sporeless strain is an important breeding target in the mushroom industry. However, basidiospore production in the oyster mushroom Pleurotus ostreatus has been shown to be impaired by single-gene mutations in only two meiosis-related genes, mer3 and msh4. This study proposed a strategy for identifying the genes essential for basidiospore formation after meiotic division to determine new targets for molecular breeding. RNA-seq analysis was performed to identify P. ostreatus genes that are specifically expressed in the gill tissue of fruiting bodies, where basidiospore formation occurs. Transcriptome data during fruiting development of Coprinopsis cinerea, in which the meiotic steps progress synchronously, were then used to identify genes that are active in the postmeiotic stages. Based on these comparative analyses, five P. ostreatus genes were identified. Plasmids containing expression cassettes for hygromycin B-resistance screening, Cas9, and single-guide RNA targeting each gene were introduced into the protoplasts of dikaryotic strain, PC9×#64, to generate dikaryotic gene disruptants. Among the obtained transformants, three dikaryotic pcl1 disruptants and two cro6c disruptants did not produce basidiospores. Microscopic analyses indicated that spore formation was arrested at particular stages in these gene disruptants. These results indicate that these two genes are essential for mature spore formation in this fungus.

The Bcvic1 and Bcvic2 vegetative incompatibility genes in Botrytis cinerea encode proteins with domain architectures involved in allorecognition in other filamentous fungi.

Vegetative incompatibility is a fungal allorecognition system characterised by the inability of genetically distinct conspecific fungal strains to form a viable heterokaryon and is controlled by multiple polymorphic loci termed vic (vegetative incompatibility) or het (heterokaryon incompatibility). We have genetically identified and characterised the first vic locus in the economically important, plant-pathogenic, necrotrophic fungus Botrytis cinerea. A bulked segregant approach coupled with whole genome Illumina sequencing of near-isogenic lines of B. cinerea was used to map a vic locus to a 60-kb region of the genome. Within that locus, we identified two adjacent, highly polymorphic open reading frames, Bcvic1 and Bcvic2, which encode predicted proteins that contain domain architectures implicated in vegetative incompatibility in other filamentous fungi. Bcvic1 encodes a predicted protein containing a putative serine esterase domain, a NACHT family of NTPases domain, and several Ankyrin repeats. Bcvic2 encodes a putative syntaxin protein containing a SNARE domain; such proteins typically function in vesicular transport. Deletion of Bcvic1 and Bcvic2 individually had no effect on vegetative incompatibility. However, deletion of the region containing both Bcvic1 and Bcvic2 resulted in mutant lines that were severely restricted in growth and showed loss of vegetative incompatibility. Complementation of these mutants by ectopic expression restored the growth and vegetative incompatibility phenotype, indicating that Bcvic1 and Bcvic2 are controlling vegetative incompatibility at this vic locus.

Whole Genome Resource and Genetic Analysis of Magnaporthe oryzae from Two Field Isolates in Northeast China.

To screen candidate fungal genes that may relate to avirulence genes corresponding to the host resistance genes, we characterized two field isolates of Magnaporthe oryzae that cause rice blast disease, especially in northeast China, and performed whole-genome resequencing of these two isolates. The genome assembly and annotation data was submitted to the National Center for Biotechnology Information database. Our study unveils the predicted fungal effectors of two dominant M. oryzae isolates in northeast China, providing a resource for Avr genes to clone. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification and subcellular localization analysis of membrane protein Ycf 1 in the microsporidian Nosema bombycis.

Microsporidia are obligate intracellular parasites that can infect a wide range of vertebrates and invertebrates including humans and insects, such as silkworm and bees. The microsporidium Nosema bombycis can cause pebrine in Bombyx mori, which is the most destructive disease in the sericulture industry. Although membrane proteins are involved in a wide range of cellular functions and part of many important metabolic pathways, there are rare reports about the membrane proteins of microsporidia up to now. We screened a putative membrane protein Ycf 1 from the midgut transcriptome of the N. bombycis-infected silkworm. Gene cloning and bioinformatics analysis showed that the Ycf 1 gene contains a complete open reading frame (ORF) of 969 bp in length encoding a 322 amino acid polypeptide that has one signal peptide and one transmembrane domain. Indirect immunofluorescence results showed that Ycf 1 protein is distributed on the plasma membrane. Expression pattern analysis showed that the Ycf 1 gene expressed in all developmental stages of N. bombycis. Knockdown of the Ycf 1 gene by RNAi effectively inhibited the proliferation of N. bombycis. These results indicated that Ycf 1 is a membrane protein and plays an important role in the life cycle of N. bombycis.

Synonymous mutations in representative yeast genes are mostly strongly non-neutral.

Synonymous mutations in protein-coding genes do not alter protein sequences and are thus generally presumed to be neutral or nearly neutral1-5. Here, to experimentally verify this presumption, we constructed 8,341 yeast mutants each carrying a synonymous, nonsynonymous or nonsense mutation in one of 21 endogenous genes with diverse functions and expression levels and measured their fitness relative to the wild type in a rich medium. Three-quarters of synonymous mutations resulted in a significant reduction in fitness, and the distribution of fitness effects was overall similar-albeit nonidentical-between synonymous and nonsynonymous mutations. Both synonymous and nonsynonymous mutations frequently disturbed the level of mRNA expression of the mutated gene, and the extent of the disturbance partially predicted the fitness effect. Investigations in additional environments revealed greater across-environment fitness variations for nonsynonymous mutants than for synonymous mutants despite their similar fitness distributions in each environment, suggesting that a smaller proportion of nonsynonymous mutants than synonymous mutants are always non-deleterious in a changing environment to permit fixation, potentially explaining the common observation of substantially lower nonsynonymous than synonymous substitution rates. The strong non-neutrality of most synonymous mutations, if it holds true for other genes and in other organisms, would require re-examination of numerous biological conclusions about mutation, selection, effective population size, divergence time and disease mechanisms that rely on the assumption that synoymous mutations are neutral.

Transcriptome analysis of genes associated with autolysis of Coprinus comatus.

Coprinus comatus, widely known as "Jituigu", is an important commodity and food in China. The yield of C. comatus, however, is substantially reduced by the autolysis of the fruiting bodies after harvest. To gain insight into the molecular mechanism underlying this autolysis, we divided the growth of C. comatus fruiting bodies into four stages: infant stage (I), mature stage (M), discolored stage (D), and autolysis stage (A). We then subjected these stages to de novo transcriptomic analysis using high-throughput Illumina sequencing. A total of 12,946 unigenes were annotated and analyzed with the Gene Ontology (GO), Clusters of Orthologous Groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). We analyzed the differentially expressed genes (DEGs) between stages I and M, M and D, and D and A. Because the changes from M to D are thought to be related to autolysis, we focused on the DEGs between these two stages. We found that the pathways related to metabolic activity began to vary in the transition from M to D, including pathways named as autophagy-yeast, peroxisome, and starch and sucrose metabolism. This study also speculates the possible process of the autolysis of Coprinus comatus. In addition, 20 genes of interest were analyzed by quantitative real-time PCR to verify their expression profiles at the four developmental stages. This study, which is the first to describe the transcriptome of C. comatus, provides a foundation for future studies concerning the molecular basis of the autolysis of its fruiting bodies.

Improvement of laccase activity by silencing PacC in Ganoderma lucidum.

Ganoderma lucidum is a representative white-rot fungus that has great potential to degrade lignocellulose biomass. Laccase is recognized as a class of the most important lignin-degrading enzymes in G. lucidum. However, the comprehensive regulatory mechanisms of laccase are still lacking. Based on the genome sequence of G. lucidum, 15 laccase genes were identified and their encoding proteins were analyzed in this study. All of the laccase proteins are predicted to be multicopper oxidases with conserved copper-binding domains. Most laccase proteins were secreted enzymes in addition to Lac14 in which the signal peptide could not be predicted. The activity of all laccases showed the highest level at pH 3.0 or pH 7.0, with total laccase activity of approximately 200 U/mg protein. Silencing PacC resulted in a 5.2 fold increase in laccase activity compared with WT. Five laccase genes (lac1, lac6, lac9, lac10 and lac14) showed an increased transcription levels (approximately 1.5-5.6 fold) in the PacC-silenced strains versus that in WT, while other laccase genes were downregulated or unchanged. The extracellular pH value was about 3.1, which was more acidic in the PacC-silenced strains than in the WT (pH 3.5). Moreover, maintaining the fermentation pH resulted in a downregulation of laccase activity which is induced by silencing PacC. Our findings indicate that in addition to its function in acidification of environmental pH, PacC plays an important role in regulating laccase activity in fungi.

Heterologous expression of a natural product biosynthetic gene cluster from Cordyceps militaris.

Cordyceps is a genus of ascomycete fungi widely used in old Chinese medicine, and many investigations have focus on uncovering their biological activities. Until now, only a few compounds have been identified from Cordyceps, mainly due to their poor yield. So as to make full use of Cordyceps, we used the strategy of genome mining and heterologous expression to discover natural products (NPs) from Cordyceps militaris. Analysis of the genome sequence of Cordyceps militaris CM01 showed the presence of a cryptic gene cluster encoding a highly-reducing polyketide synthetase (HR-PKS), enoyl-reductase (ER) and cytochrome P450. Heterologous expression in Aspergillus nidulans enabled the identification of two new polyketides, cordypyrone A and B. Their structures were determined by 1D and 2D NMR techniques. They showed only modest activities against pathogenic bacteria including methicillin-resistant Staphylococcus aureus (MRSA), Mycobacteria tuberculosis and Bacillus cereus.

The cold - resistance mechanism of a mutagenic Volvariella volvacea strain VH3 with outstanding traits revealed by transcriptome profiling.

BACKGROUND: The straw mushroom (Volvariella volvacea) is one of the important vegetables that is popular for its delicious taste. However, the straw mushroom is sensitive to low temperature, resulting in economic loss during transportation and storage. We obtained a novel straw mushroom strain, named VH3, via ultraviolet mutagenesis. RESULTS: Our study revealed that VH3 exhibited high cold resistance compared to an ordinary straw mushroom cultivar, V23. We found that the electrolyte leakages of VH3 were always significantly lower than that of V23 treated with 4 °C for 0 h, 2 h,4 h, 8 h, 16 h, and 24 h. Before cold treatment (0 h), there were no difference of MDA contents, SOD activities, and CAT activities between VH3 and V23. At the late stage (8 h, 26 h, and 24 h) of cold treatment, the MDA contents of VH3 were lower while both the SOD and CAT activities were higher than those of V23. To investigate the potential mechanisms of VH3 cold resistance, we performed transcriptome sequencing to detect the transcriptome profiling of VH3 and V23 after 0 h and 4 h cold treatment. Transcriptome sequencing revealed that 111 differentially expressed genes (DEG) between V23 (0 h) and VH3 (0 h) (V23-0_vs_VH3-0), consisting 50 up-regulated and 61 down-regulated DEGs. A total of 117 DEGs were obtained between V23 (4 h) and VH3(4 h) (V23-4_vs_VH3-4), containing 94 up-regulated and 23 down-regulated DEGs. Among these DEGs, VVO_00021 and VVO_00017 were up-regulated while VVO_00003, VVO_00004, VVO_00010, and VVO_00030 were down-regulated in V23-0_vs_VH3-0 and VH3-4_vs_V23-4. KEGG and GO analysis revealed that the 6 DEGs were annotated to pathways related to cold stress. Besides, the GA3 content was also decreased in VH3. CONCLUSIONS: Collectively, our study first revealed that the increased cold resistance of VH3 might be caused by the expression change of VVO_00003, VVO_00004, VVO_00017, VVO_00021, and VVO_00030, and decreased GA3.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic.

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

Mining Indole Alkaloid Synthesis Gene Clusters from Genomes of 53 Claviceps Strains Revealed Redundant Gene Copies and an Approximate Evolutionary Hourglass Model.

Ergot fungi (Claviceps spp.) are infamous for producing sclerotia containing a wide spectrum of ergot alkaloids (EA) toxic to humans and animals, making them nefarious villains in the agricultural and food industries, but also treasures for pharmaceuticals. In addition to three classes of EAs, several species also produce paspaline-derived indole diterpenes (IDT) that cause ataxia and staggers in livestock. Furthermore, two other types of alkaloids, i.e., loline (LOL) and peramine (PER), found in Epichloë spp., close relatives of Claviceps, have shown beneficial effects on host plants without evidence of toxicity to mammals. The gene clusters associated with the production of these alkaloids are known. We examined genomes of 53 strains of 19 Claviceps spp. to screen for these genes, aiming to understand the evolutionary patterns of these genes across the genus through phylogenetic and DNA polymorphism analyses. Our results showed (1) varied numbers of eas genes in C. sect. Claviceps and sect. Pusillae, none in sect. Citrinae, six idt/ltm genes in sect. Claviceps (except four in C. cyperi), zero to one partial (idtG) in sect. Pusillae, and four in sect. Citrinae, (2) two to three copies of dmaW, easE, easF, idt/ltmB, itd/ltmQ in sect. Claviceps, (3) frequent gene gains and losses, and (4) an evolutionary hourglass pattern in the intra-specific eas gene diversity and divergence in C. purpurea.

Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression with cell volume.

While most genes' expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiae's mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model.

Decreasing mitochondrial RNA polymerase activity reverses biased inheritance of hypersuppressive mtDNA.

Faithful inheritance of mitochondrial DNA (mtDNA) is crucial for cellular respiration/oxidative phosphorylation and mitochondrial membrane potential. However, how mtDNA is transmitted to progeny is not fully understood. We utilized hypersuppressive mtDNA, a class of respiratory deficient Saccharomyces cerevisiae mtDNA that is preferentially inherited over wild-type mtDNA (rho+), to uncover the factors governing mtDNA inheritance. We found that some regions of rho+ mtDNA persisted while others were lost after a specific hypersuppressive takeover indicating that hypersuppressive preferential inheritance may partially be due to active destruction of rho+ mtDNA. From a multicopy suppression screen, we found that overexpression of putative mitochondrial RNA exonuclease PET127 reduced biased inheritance of a subset of hypersuppressive genomes. This suppression required PET127 binding to the mitochondrial RNA polymerase RPO41 but not PET127 exonuclease activity. A temperature-sensitive allele of RPO41 improved rho+ mtDNA inheritance over a specific hypersuppressive mtDNA at semi-permissive temperatures revealing a previously unknown role for rho+ transcription in promoting hypersuppressive mtDNA inheritance.

Modulation of heterologous protein secretion in the thermotolerant methylotrophic yeast Ogataea thermomethanolica TBRC 656 by CRISPR-Cas9 system.

The thermotolerant methylotrophic yeast Ogataea thermomethanolica TBRC 656 is a potential host strain for industrial protein production. Heterologous proteins are often retained intracellularly in yeast resulting in endoplasmic reticulum (ER) stress and poor secretion, and despite efforts to engineer protein secretory pathways, heterologous protein production is often lower than expected. We hypothesized that activation of genes involved in the secretory pathway could mitigate ER stress. In this study, we created mutants defective in protein secretory-related functions using clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) tools. Secretion of the model protein xylanase was significantly decreased in loss of function mutants for oxidative stress (sod1Δ) and vacuolar and protein sorting (vps1Δ and ypt7Δ) genes. However, xylanase secretion was unaffected in an autophagy related atg12Δ mutant. Then, we developed a system for sequence-specific activation of target gene expression (CRISPRa) in O. thermomethanolica and used it to activate SOD1, VPS1 and YPT7 genes. Production of both non-glycosylated xylanase and glycosylated phytase was enhanced in the gene activated mutants, demonstrating that CRISPR-Cas9 systems can be used as tools for understanding O. thermomethanolica genes involved in protein secretion, which could be applied for increasing heterologous protein secretion in this yeast.

Farming systems influence the compositional, structural, and functional characteristics of the sugarcane-associated microbiome.

Sugarcane (Saccharum spp.) has been produced worldwide as a relevant source of food and sustainable energy. However, the constant need to increase crop yield has led to excessive use of synthetic agrochemical inputs such as inorganic fertilizers, herbicides, and pesticides in plant cultures. It is known that these conventional practices can lead to deleterious effects on health and the environment. Organic farming emerges as a sustainable alternative to conventional systems; however, farm management influences in plant-associated microbiomes remain unclear. Here, the aim is to identify the effects of farming systems on the sugarcane microbiota. To address this issue, we sampled the microbiota from soils and plants under organic and conventional farming from two crop fields in Brazil. Then, we evaluated their compositional, structural, and functional traits through amplification and sequencing of phylogenetic markers of bacteria (16S rRNA gene, V3-V4 region) and fungi (Internal Transcribed Spacer - ITS2). The data processing and analyses by the DADA2 pipeline revealed 12,839 bacterial and 3,222 fungal sequence variants. Moreover, differences between analogous niches were detected considering the contrasting farming systems, with samples from the conventional system showing a slightly greater richness and diversity of microorganisms. The composition is also different between the farming systems, with 389 and 401 differentially abundant taxa for bacteria and fungi, respectively, including taxa capable of promoting plant growth. The microbial co-occurrence networks showed structural changes in microbial communities, where organic networks were more cohesive since they had closer taxa and less modularity by niches. Finally, the functional prediction revealed enriched metabolic pathways, including the increased presence of antimicrobial resistance in the conventional farming system. Taken together, our findings reveal functional, structural, and compositional adaptations of the microbial communities associated with sugarcane plants in the field, according to farming management. With this, we point out the need to unravel the mechanisms driving these adaptations.

Development of an efficient gene-targeting system for elucidating infection mechanisms of the fungal pathogen Trichosporon asahii.

Trichosporon asahii is a pathogenic fungus that causes severe, deep-seated fungal infections in neutropenic patients. Elucidating the infection mechanisms of T. asahii based on genetic studies requires a specific gene-targeting system. Here, we established an efficient gene-targeting system in a highly pathogenic T. asahii strain identified using the silkworm infection model. By comparing the pathogenicity of T. asahii clinical isolates in a silkworm infection model, T. asahii MPU129 was identified as a highly pathogenic strain. Using an Agrobacterium tumefaciens-mediated gene transfer system, we obtained a T. asahii MPU129 mutant lacking the ku70 gene, which encodes the Ku70 protein involved in the non-homologous end-joining repair of DNA double-strand breaks. The ku70 gene-deficient mutant showed higher gene-targeting efficiency than the wild-type strain for constructing a mutant lacking the cnb1 gene, which encodes the beta-subunit of calcineurin. The cnb1 gene-deficient mutant showed reduced pathogenicity against silkworms compared with the parental strain. These results suggest that an efficient gene-targeting system in a highly pathogenic T. asahii strain is a useful tool for elucidating the molecular mechanisms of T. asahii infection.

Genetic Analysis of the Hsm3 Protein Function in Yeast Saccharomyces cerevisiae NuB4 Complex.

In the nuclear compartment of yeast, NuB4 core complex consists of three proteins, Hat1, Hat2, and Hif1, and interacts with a number of other factors. In particular, it was shown that NuB4 complex physically interacts with Hsm3p. Early we demonstrated that the gene HSM3 participates in the control of replicative and reparative spontaneous mutagenesis, and that hsm3Δ mutants increase the frequency of mutations induced by different mutagens. It was previously believed that the HSM3 gene controlled only some minor repair processes in the cell, but later it was suggested that it had a chaperone function with its participation in proteasome assembly. In this work, we analyzed the properties of three hsm3Δ, hif1Δ, and hat1Δ mutants. The results obtained showed that the Hsm3 protein may be a functional subunit of NuB4 complex. It has been shown that hsm3- and hif1-dependent UV-induced mutagenesis is completely suppressed by inactivation of the Polη polymerase. We showed a significant role of Polη for hsm3-dependent mutagenesis at non-bipyrimidine sites (NBP sites). The efficiency of expression of RNR (RiboNucleotid Reducase) genes after UV irradiation in hsm3Δ and hif1Δ mutants was several times lower than in wild-type cells. Thus, we have presented evidence that significant increase in the dNTP levels suppress hsm3- and hif1-dependent mutagenesis and Polη is responsible for hsm3- and hif1-dependent mutagenesis.

Spontaneous changes in somatic compatibility in Fusarium circinatum.

Filamentous fungi grow by the elaboration of hyphae, which may fuse to form a network as a colony develops. Fusion of hyphae can occur between genetically different individuals, provided they share a common allele at loci affecting somatic compatibility. Diversity in somatic compatibility phenotypes reduces the frequency of hyphal fusion in a population, thereby slowing the spread of deleterious genetic elements such as viruses and plasmids, which require direct cytoplasmic contact for transmission. Diverse somatic compatibility phenotypes can be generated by recombining alleles through sexual reproduction, but this mechanism may not fully account for the diversity found in nature. For example, multiple compatibility phenotypes of Fusarium circinatum were shown to be associated with the same clonal lineage, which implies they were derived by a mutation rather than recombination through sexual reproduction. Experimental tests of this hypothesis confirmed that spontaneous changes in somatic compatibility can occur at a frequency between 5 and 8 per million spores. Genomic analysis of F. circinatum strains with altered somatic compatibility revealed no consistent evidence of recombination and supported the hypothesis that a spontaneous mutation generated the observed phenotypic change. Genes known to be involved in somatic compatibility had no mutations, suggesting that mutation occurred in a gene with an as yet unexplored function in somatic compatibility.

Phylogeny and mycotoxin profile of Fusarium species isolated from sugarcane in Southern Iran.

Sugarcane is an important crop in Southern Iran for agri-food, energy, and pharmaceutical industries. Among the pathogens that colonize sugarcane, mycotoxigenic Fusarium species are reason of serious concern for both their pathogenicity on plants and ability to produce harmful mycotoxins to humans and animals. We studied 104 Fusarium strains, selected within a wider Fusarium set isolated from sugarcane in Southern Iran, for molecular identification, phylogeny and mycotoxin analyses. Most of Fusarium strains belonged to Fusarium fujikuroi Species Complex (FFSC) and identified mainly as F. proliferatum, at minor extent as F. sacchari, and rarely as F. thapsinum, and F. verticillioides. Moreover, 14 strains identified as FFSC could not be assigned to any known species, although they were phylogenetically closely related to F. andiyazi, likely representing a new phylogenetic species. A subset of FFSC strains were analyzed for in vitro production of fumonisins (FBs), beauvericin (BEA), and enniatins (ENNs). Fusarium proliferatum strains produced FBs at high amount, and, at a lesser extent, BEA, and ENNs; F.sacchari produced only BEA and B ENNs at very low level; Fusarium sp. strains produced only B ENNs. The paper provides new insights on the genetic diversity of Fusarium species and their mycotoxin profile occurring on sugarcane in Iran.