Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

Combined Use of Phenotype-Based and Genome-Informed Approaches Identified a Unique Fusarium oxysporum f. sp. cubense Isolate in Hawaii

Fusarium wilt of banana, caused by Fusarium oxysporum f. sp. cubense (Foc), is a serious disease that threatens banana production worldwide. It is a long-standing problem in Hawaii, but previously, there was little knowledge of the causal pathogen. We isolated a strain of Foc, named Foc-UH, from a field experiencing the disease epidemic in Hawaii. Infection assays of a diverse panel of 26 banana clones, including varieties used for differentiating pathogen races and fruit production, revealed that Foc-UH has a race 1 pathogenic phenotype with an intermediate race 2 virulence and revealed the differential resistance of varieties to infection. Separate phylogenetic analyses using the barcoding regions of three nuclear genes, seven complete nuclear genes, and single-nucleotide polymorphisms within conserved whole-genome protein coding sequences placed Foc-UH into recently proposed taxonomic frameworks relevant to Foc and the F. oxysporum species complex. Screening of the 99.7% complete draft genome identified five secreted in xylem (SIX) gene homologs: SIX1d, SIX1f, SIX9a, SIX9b, and SIX13a. This profile is similar to that of several race 1 isolates except for the absence of SIX4 and SIX6. Foc-UH was morphologically dissimilar to the nearest related isolates. Altogether, this study identified a unique isolate that causes banana Fusarium wilt, which represents the first characterization of the causal pathogen in Hawaii. The findings and genomic resources generated in this study are expected to guide banana breeding and cultivar deployment in Hawaii and beyond and contribute to further understanding of the pathogenicity and evolutionary systematics of Foc.

The differential virulence of Fusarium strains causing corneal infections and plant diseases is associated with accessory chromosome composition

Fusarium oxysporum is a cross-kingdom pathogen. While some strains cause disseminated fusariosis and blinding corneal infections in humans, others are responsible for devastating vascular wilt diseases in plants. To better understand the distinct adaptations of F. oxysporum to animal or plant hosts, we conducted a comparative phenotypic and genetic analysis of two strains: MRL8996 (isolated from a keratitis patient) and Fol4287 (isolated from a wilted tomato [Solanum lycopersicum]). In vivo infection of mouse corneas and tomato plants revealed that, while both strains cause symptoms in both hosts, MRL8996 caused more severe corneal ulceration and perforation in mice, whereas Fol4287 induced more pronounced wilting symptoms in tomato. In vitro assays using abiotic stress treatments revealed that the human pathogen MRL8996 was better adapted to elevated temperatures, whereas the plant pathogen Fol4287 was more tolerant of osmotic and cell wall stresses. Both strains displayed broad resistance to antifungal treatment, with MRL8996 exhibiting the paradoxical effect of increased tolerance to higher concentrations of the antifungal caspofungin. We identified a set of accessory chromosomes (ACs) and protein-encoding genes with distinct transposon profiles and functions, respectively, between MRL8996 and Fol4287. Interestingly, ACs from both genomes also encode proteins with shared functions, such as chromatin remodeling and post-translational protein modifications. Our phenotypic assays and comparative genomics analyses lay the foundation for future studies correlating genotype with phenotype and for developing targeted antifungals for agricultural and clinical uses.

Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

Group I Intron as a Potential Target for Antifungal Compounds: Development of a Trans-Splicing High-Throughput Screening Strategy

The search for safe and efficient new antifungal compounds for agriculture has led to more efforts in finding new modes of action. This involves the discovery of new molecular targets, including coding and non-coding RNA. Rarely found in plants and animals but present in fungi, group I introns are of interest as their complex tertiary structure may allow selective targeting using small molecules. In this work, we demonstrate that group I introns present in phytopathogenic fungi have a self-splicing activity in vitro that can be adapted in a high-throughput screening to find new antifungal compounds. Ten candidate introns from different filamentous fungi were tested and one group ID intron found in F. oxysporum showed high self-splicing efficiency in vitro. We designed the Fusarium intron to act as a trans-acting ribozyme and used a fluorescence-based reporter system to monitor its real time splicing activity. Together, these results are opening the way to study the druggability of such introns in crop pathogen and potentially discover small molecules selectively targeting group I introns in future high-throughput screenings.

Endophytic Trichoderma species from rubber trees native to the Brazilian Amazon, including four new species

Fungi belonging to the genus Trichoderma have been widely recognized as efficient controllers of plant diseases. Although the majority of isolates currently deployed, thus far, have been isolated from soil, endophytic Trichoderma spp. is considered to be a promising option for application in biocontrol. In this study, 30 endophytic Trichoderma isolates-obtained from the leaves, stems, and roots of wild Hevea spp. in the Brazilian Amazon-were analyzed using specific DNA barcodes: sequences of internal transcribed spacers 1 and 2 of rDNA (ITS region), genes encoding translation elongation factor 1-α (TEF1-α), and the second largest subunit of RNA polymerase II (RPB2). The genealogical concordance phylogenetic species recognition (GCPSR) concept was used for species delimitation. A phylogenetic analysis showed the occurrence of Trichoderma species, such as T. erinaceum, T. ovalisporum, T. koningiopsis, T. sparsum, T. lentiforme, T. virens, and T. spirale. Molecular and morphological features resulted in the discovery of four new species, such as T. acreanum sp. nov., T. ararianum sp. nov., T. heveae sp. nov., and T. brasiliensis sp. nov. The BI and ML analyses shared a similar topology, providing high support to the final trees. The phylograms show three distinct subclades, namely, T. acreanum and T. ararianum being paraphyletic with T. koningiopsis; T. heveae with T. subviride; and T. brasiliensis with T. brevicompactum. This study adds to our knowledge of the diversity of endophytic Trichoderma species in Neotropical forests and reveals new potential biocontrol agents for the management of plant diseases.

Single-molecule Sequencing of an Animal Mitochondrial Genome Reveals Chloroplast-like Architecture and Repeat-mediated Recombination

Recent advances in long-read sequencing technology have allowed for single-molecule sequencing of entire mitochondrial genomes, opening the door for direct investigation of the mitochondrial genome architecture and recombination. We used PacBio sequencing to reassemble mitochondrial genomes from two species of New Zealand freshwater snails, Potamopyrgus antipodarum and Potamopyrgus estuarinus. These assemblies revealed a ∼1.7 kb structure within the mitochondrial genomes of both species that was previously undetected by an assembly of short reads and likely corresponding to a large noncoding region commonly present in the mitochondrial genomes. The overall architecture of these Potamopyrgus mitochondrial genomes is reminiscent of the chloroplast genomes of land plants, harboring a large single-copy (LSC) region and a small single-copy (SSC) region separated by a pair of inverted repeats (IRa and IRb). Individual sequencing reads that spanned across the Potamopyrgus IRa-SSC-IRb structure revealed the occurrence of a "flip-flop" recombination. We also detected evidence for two distinct IR haplotypes and recombination between them in wild-caught P. estuarinus, as well as extensive intermolecular recombination between single-nucleotide polymorphisms in the LSC region. The chloroplast-like architecture and repeat-mediated mitochondrial recombination we describe here raise fundamental questions regarding the origins and commonness of inverted repeats in cytoplasmic genomes and their role in mitochondrial genome evolution.

Multi-site fungicides suppress banana Panama disease, caused by Fusarium oxysporum f. sp. cubense Tropical Race 4

Global banana production is currently challenged by Panama disease, caused by Fusarium oxysporum f.sp. cubense Tropical Race 4 (FocTR4). There are no effective fungicide-based strategies to control this soil-borne pathogen. This could be due to insensitivity of the pathogen to fungicides and/or soil application per se. Here, we test the effect of 12 single-site and 9 multi-site fungicides against FocTR4 and Foc Race1 (FocR1) in quantitative colony growth, and cell survival assays in purified FocTR4 macroconidia, microconidia and chlamydospores. We demonstrate that these FocTR4 morphotypes all cause Panama disease in bananas. These experiments reveal innate resistance of FocTR4 to all single-site fungicides, with neither azoles, nor succinate dehydrogenase inhibitors (SDHIs), strobilurins or benzimidazoles killing these spore forms. We show in fungicide-treated hyphae that this innate resistance occurs in a subpopulation of "persister" cells and is not genetically inherited. FocTR4 persisters respond to 3 μg ml-1 azoles or 1000 μg ml-1 strobilurins or SDHIs by strong up-regulation of genes encoding target enzymes (up to 660-fold), genes for putative efflux pumps and transporters (up to 230-fold) and xenobiotic detoxification enzymes (up to 200-fold). Comparison of gene expression in FocTR4 and Zymoseptoria tritici, grown under identical conditions, reveals that this response is only observed in FocTR4. In contrast, FocTR4 shows little innate resistance to most multi-site fungicides. However, quantitative virulence assays, in soil-grown bananas, reveals that only captan (20 μg ml-1) and all lipophilic cations (200 μg ml-1) suppress Panama disease effectively. These fungicides could help protect bananas from future yield losses by FocTR4.

Phylogenetic Diversity and Mycotoxin Potential of Emergent Phytopathogens Within the Fusarium tricinctum Species Complex

Recent studies on multiple continents indicate members of the Fusarium tricinctum species complex (FTSC) are emerging as prevalent pathogens of small-grain cereals, pulses, and other economically important crops. These understudied fusaria produce structurally diverse mycotoxins, among which enniatins (ENNs) and moniliformin (MON) are the most frequent and of greatest concern to food and feed safety. Herein a large survey of fusaria in the Fusarium Research Center and Agricultural Research Service culture collections was undertaken to assess species diversity and mycotoxin potential within the FTSC. A 151-strain collection originating from diverse hosts and substrates from different agroclimatic regions throughout the world was selected from 460 FTSC strains to represent the breadth of FTSC phylogenetic diversity. Evolutionary relationships inferred from a five-locus dataset, using maximum likelihood and parsimony, resolved the 151 strains as 24 phylogenetically distinct species, including nine that are new to science. Of the five genes analyzed, nearly full-length phosphate permease sequences contained the most phylogenetically informative characters, establishing its suitability for species-level phylogenetics within the FTSC. Fifteen of the species produced ENNs, MON, the sphingosine analog 2-amino-14,16-dimethyloctadecan-3-ol (AOD), and the toxic pigment aurofusarin (AUR) on a cracked corn kernel substrate. Interestingly, the five earliest diverging species in the FTSC phylogeny (i.e., F. iranicum, F. flocciferum, F. torulosum, and Fusarium spp. FTSC 8 and 24) failed to produce AOD and MON, but synthesized ENNs and/or AUR. Moreover, our reassessment of nine published phylogenetic studies on the FTSC identified 11 additional novel taxa, suggesting this complex comprises at least 36 species.

Identification of causal agent of wilt of common sage (Salvia officinalis L.)

Introduction: Common sage is cultivated in Europe and North America. It has strong antiviral, antibacterial and antifungal properties. This plant can be infected by different pathogenic fungi species, such as Alternaria alternata, Fusarium spp. (F. culmorum, F. equiseti, F. oxysporum), Phomopsis sclarea and Botrytis cinerea. Those species are the most frequently isolated fungi from sage stem base.

Objective: The aim of this study was to identify the causal agent of common sage wilt disease.

Methods: Studies were carried out in 2018–2020. 23 fungal isolates were identified based on their morphology and with use of PCR technique. Length and width of 100 conidia growing on SNA medium were measured after 7 days. Koch’s postulates were checked and the development of one fungus isolate (no. 13) was compared on seven media: the CMA, MEA, OA, PCA, SNA, PDA and Czapek medium. Sequences of the second largest subunit of RNA polymerase II (RPB2) were used to identify the pathogen.

Results: The fungus formed 3 kinds of spores: thin-walled, hyaline, slightly folded at the base, mostly 4-cell macroconidia, oblong, hyaline one- or two-cell microconidia and oval thick-walled chlamydospores. The Koch’s postulates were fulfilled. The fungus formed the most abundant aerial mycelium on the Czapek medium, and the least on the CMA medium. On the SNA medium, the mycelium grew into the medium and the aerial mycelium was not formed. The obtained RPB2 nucleotide sequence was 100% similar to the Fusarium oxysporum sequence deposited in GenBank (NCBI).

Conclusions: The results of this research can be used in further studies on the biological diversity of this species.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens

Exserohilum turcicum and E. rostratum, two closely related fungal species, are both economically important pathogens but have quite different target hosts (specific to plants and cross-kingdom infection, respectively). In the present study, complete circular mitochondrial genomes of the two Exserohilum species were sequenced and de novo assembled, which mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). Comparative analyses indicated that these two fungi had significant mitogenomic collinearity and consistent mitochondrial gene arrangement, yet with vastly different mitogenome sizes, 264,948 bp and 64,620 bp, respectively. By contrast with the 17 introns containing 17 intronic ORFs (one-to-one) in the E. rostratum mitogenome, E. turcicum involved far more introns (70) and intronic ORFs (126), which was considered as the main contributing factors of their mitogenome expansion/contraction. Within the generally intron-rich gene cox1, a total of 18 and 10 intron position classes (Pcls) were identified separately in the two mitogenomes. Moreover, 16.16% and 10.85% ratios of intra-mitogenomic repetitive regions were detected in E. turcicum and E. rostratum, respectively. Based on the combined mitochondrial gene dataset, we established a well-supported topology of phylogeny tree of 98 ascomycetes, implying that mitogenomes may act as an effective molecular marker for fungal phylogenetic reconstruction. Our results served as the first report on mitogenomes in the genus Exserohilum, and would have significant implications in understanding the origin, evolution and pathogenic mechanisms of this fungal lineage.

Identification and Characterization of Fusarium nirenbergiae Associated with Saffron Corm Rot Disease

Saffron (Crocus sativus L.) is the most expensive spice plant and is distributed widely around the world. However, its production is limited by corm rot, a disastrous disease, attributed to Fusarium oxysporum in many regions of the world. In 2020, extensive surveys were carried out in Zhejiang, Shanghai, Anhui, and Guizhou provinces as saffron growing areas of China. Fourteen single-spore isolates were obtained from rotted corms and identified as F. nirenbergiae according to morphological appearance and multilocus phylogenetic analysis with translation elongation factor 1-α (tef1), DNA-directed RNA polymerase II largest subunit (rpb2), and β-tubulin (tub2). Results of the pathogenicity assay supported the conclusion that F. nirenbergiae is the pathogen responsible for corm rot. In this study, we obtained the whole genome sequence of two highly virulent F. nirenbergiae strains via the Illumina HiSeq platform. Genome sequence assemblies of approximately 52.7 and 52.2 Mb were generated for isolates WY5 and SH1, respectively. To the best of our knowledge, this is the first report of F. nirenbergiae causing C. sativus corm rot in China and indeed worldwide. Results from this research contribute to our understanding of genetic diversity, genomic information, and host determination, which will enable researchers to design appropriate management measures for this hazardous disease.

A Polyphasic Approach Reveals Novel Genotypes and Updates the Genetic Structure of the Banana Fusarium Wilt Pathogen

Fusarium oxysporum f. sp. cubense (Foc) is a soil-borne fungus that causes Fusarium wilt, a destructive plant disease that has resulted in devastating economic losses to banana production worldwide. The fungus has a complex evolutionary history and taxonomic repute and consists of three pathogenic races and at least 24 vegetative compatibility groups (VCGs). Surveys conducted in Asia, Africa, the Sultanate of Oman and Mauritius encountered isolates of F. oxysporum pathogenic to banana that were not compatible to any of the known Foc VCGs. Genetic relatedness between the undescribed and known Foc VCGs were determined using a multi-gene phylogeny and diversity array technology (DArT) sequencing. The presence of putative effector genes, the secreted in xylem (SIX) genes, were also determined. Fourteen novel Foc VCGs and 17 single-member VCGs were identified. The multi-gene tree was congruent with the DArT-seq phylogeny and divided the novel VCGs into three clades. Clustering analysis of the DArT-seq data supported the separation of Foc isolates into eight distinct clusters, with the suite of SIX genes mostly conserved within these clusters. Results from this study indicates that Foc is more diverse than hitherto assumed.

Population genomics reveals historical and ongoing recombination in the Fusarium oxysporum species complex

The Fusarium oxysporum species complex (FOSC) is a group of closely related plant pathogens long-considered strictly clonal, as sexual stages have never been recorded. Several studies have questioned whether recombination occurs in FOSC, and if it occurs its nature and frequency are unknown. We analysed 410 assembled genomes to answer whether FOSC diversified by occasional sexual reproduction interspersed with numerous cycles of asexual reproduction akin to a model of predominant clonal evolution (PCE). We tested the hypothesis that sexual reproduction occurred in the evolutionary history of FOSC by examining the distribution of idiomorphs at the mating locus, phylogenetic conflict and independent measures of recombination from genome-wide SNPs and genes. A phylogenomic dataset of 40 single copy orthologs was used to define structure a priori within FOSC based on genealogical concordance. Recombination within FOSC was tested using the pairwise homoplasy index and divergence ages were estimated by molecular dating. We called SNPs from assembled genomes using a k-mer approach and tested for significant linkage disequilibrium as an indication of PCE. We clone-corrected and tested whether SNPs were randomly associated as an indication of recombination. Our analyses provide evidence for sexual or parasexual reproduction within, but not between, clades of FOSC that diversified from a most recent common ancestor about 500 000 years ago. There was no evidence of substructure based on geography or host that might indicate how clades diversified. Competing evolutionary hypotheses for FOSC are discussed in the context of our results.

LAMP Assay for Distinguishing Fusarium oxysporum and Fusarium commune in Lotus (Nelumbo nucifera) Rhizomes

Yields of edible rhizome from cultivation of the perennial hydrophyte lotus (Nelumbo nucifera) can be severely reduced by rhizome rot disease caused by Fusarium species. There is a lack of rapid field-applicable methods for detection of these pathogens on lotus plants displaying symptoms of rhizome rot. Fusarium commune (91%) and Fusarium oxysporum (9%) were identified at different frequencies from lotus samples showing symptoms of rhizome rot. Because these two species can cause different severity of disease and their morphology is similar, molecular diagnostic-based methods to detect these two species were developed. Based on the comparison of the mitochondrial genome of the two species, three specific DNA loci targets were found. The designed primer sets for conventional PCR, quantitative PCR, and loop-mediated isothermal amplification (LAMP) precisely distinguished the above two species when isolated from lotus and other plants. The LAMP detection limits were 10 pg/μl and 1 pg/μl of total DNA for F. commune and F. oxysporum, respectively. We also carried out field-mimicked experiments on lotus seedlings and rhizomes (including inoculated samples and field-diseased samples), and the results indicated that the LAMP primer sets and the supporting portable methods are suitable for rapid diagnosis of the lotus disease in the field. The LAMP-based detection method will aid in the rapid identification of whether F. oxysporum or F. commune is infecting lotus plants with symptoms of rhizome rot and can facilitate efficient pesticide use and prevent disease spread through vegetative propagation of Fusarium-infected lotus rhizomes.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

Taxonomic Revision of the Banana Fusarium Wilt TR4 Pathogen Is Premature

Taxonomic revisions for pathogens of crops should be based on robust underpinning evidence. Recently, a substantial revision was proposed for the taxonomy of the causative agent of Fusarium wilt on banana. We reanalyzed the data on which this revision was based and discovered that the data do not robustly support the proposals. Several apparent discrepancies and errors in the published phylogenies cast further doubt on the conclusions drawn from them. Although we do not assert that the authors' conclusions are incorrect, we posit that the taxonomic changes are premature, given the data currently in the public domain.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Effectiveness of the Influence of Selected Essential Oils on the Growth of Parasitic Fusarium Isolated from Wheat Kernels from Central Europe

The aim of the study was to determine the effectiveness of selected seven commercial essential oils (EsO) (grapefruit, lemongrass, tea tree (TTO), thyme, verbena, cajeput, and Litsea cubeba) on isolates of common Central European parasitic fungal species of Fusarium obtained from infected wheat kernels, and to evaluate the oils as potential natural fungicides. The study was conducted in 2 stages. At each stage, the fungicidal activity of EsO (with concentrations of 0.025; 0.05; 0.125; 0.25; 0.50; 1.0, and 2.0%) against Fusarium spp. was evaluated using the disc plate method and zones of growth inhibition were measured. At the first stage, the fungistatic activity of EsO was evaluated against four species of Fusarium from the Polish population (F. avenaceum FAPL, F. culmorum FCPL, F. graminearum FGPL and F. oxysporum FOPL). The correlation coefficient between the mycelial growth rate index (T) and the fungistatic activity (FA) was calculated. At the second stage, on the basis of the mycelium growth rate index, the effectiveness of the EsO in limiting the development of Fusarium isolates from the German population (F. culmorum FC1D, F. culmorum FC2D, F. graminearum FG1D, F. graminearum FG2D and F. poae FP0D) was assessed. The first and second stage results presented as a growth rate index were then used to indicate essential oils (as potential natural fungicides) effectively limiting the development of various common Central European parasitic species Fusarium spp. Finally, the sensitivity of four Fusarium isolates from the Polish population and five Fusarium isolates from the German population was compared. The data were compiled in STATISTICA 13.0 (StatSoft, Inc, CA, USA) at the significance level of 0.05. Fusarium isolates from the German population were generally more sensitive than those from the Polish population. The sensitivity of individual Fusarium species varied. Their vulnerability, regardless of the isolate origin, in order from the most to the least sensitive, is as follows: F. culmorum, F. graminearum, F. poae, F. avenaceum and F. oxysporum. The strongest fungicidal activity, similar to Funaben T, showed thyme oil (regardless of the concentration). Performance of citral oils (lemongrass and Litsea cubeba) was similar but at a concentration above 0.025%.

Uncovering Diagnostic Value of Mitogenome for Identification of Cryptic Species Fusarium graminearum Sensu Stricto

Fungal complexes are often composed of morphologically nearly indistinguishable species with high genetic similarity. However, despite their close relationship, they can exhibit distinct phenotypic differences in pathogenicity and production of mycotoxins. Many plant pathogenic and toxigenic fungi have been shown to consist of such cryptic species. Identification of cryptic species in economically important pathogens has added value in epidemiologic studies and provides opportunities for better control. Analysis of mitochondrial genomes or mitogenomics opens up dimensions for improved diagnostics of fungi, especially when efficient recovery of DNA is problematic. In comparison to nuclear DNA, mitochondrial DNA (mtDNA) can be amplified with improved efficacy due to its multi-copy nature. However, to date, only a few studies have demonstrated the usefulness of mtDNA for identification of cryptic species within fungal complexes. In this study, we explored the value of mtDNA for identification of one of the most important cereal pathogens Fusarium graminearum sensu stricto (F.g.). We found that homing endonucleases (HEGs), which are widely distributed in mitogenomes of fungi, display small indel polymorphism, proven to be potentially species specific. The resulting small differences in their lengths may facilitate further differentiation of F.g. from the other cryptic species belonging to F. graminearum species complex. We also explored the value of SNP analysis of the mitogenome for typing F.g. The success in identifying F.g. strains was estimated at 96%, making this tool an attractive complement to other techniques for identification of F.g.

Genetic Diversity of the Fusarium oxysporum Complex Isolated from the Grassland Biome of South Africa

The genetic diversity of pathogenic members of the Fusarium oxysporum species complex (FOSC) has been intensively studied worldwide, yet strains occurring in native soils with low anthropogenic disturbance remain poorly understood. This study focused on 355 F. oxysporum isolates from soils with low anthropogenic activity obtained from the grassland biome of South Africa. Analysis of the translation elongation factor 1-alpha (tef-1α) gene revealed high levels of sequence type diversity within the soil population in comparison with the global dataset. Phylogenetic relationships of the South African isolates revealed that four nested within FOSC clade 1. This is the first report of members of the basal clade recovered from ecosystems with low anthropogenic disturbance from Sub-Saharan Africa. The remaining strains nested within clades 2 to 5. This study contributes significantly to our understanding of the distribution of the FOSC in natural systems as we show that FOSC populations in the South African grassland biome are genetically diverse. This fills in our knowledge gap because previous studies reported only on the occurrence and diversity of the FOSC isolated from plant debris in South Africa. This is the first comprehensive survey of fusaria from grassland soils with low anthropogenic disturbance in South Africa.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

Genomic differences between the new Fusarium oxysporum f. sp. apii (Foa) race 4 on celery, the less virulent Foa races 2 and 3, and the avirulent on celery f. sp. coriandrii

BACKGROUND

Members of the F. oxysporium species complex (FOSC) in the f. sp. apii (Foa) are pathogenic on celery and those in f. sp. coriandrii (Foci) are pathogenic on coriander (=cilantro). Foci was first reported in California in 2005; a new and highly aggressive race 4 of Foa was observed in 2013 in California. Preliminary evidence indicated that Foa can also cause disease on coriander, albeit are less virulent than Foci. Comparative genomics was used to investigate the evolutionary relationships between Foa race 4, Foa race 3, and the Foci, which are all in FOSC Clade 2, and Foa race 2, which is in FOSC Clade 3.

RESULTS

A phylogenetic analysis of 2718 single-copy conserved genes and mitochondrial DNA sequence indicated that Foa races 3 and 4 and the Foci are monophyletic within FOSC Clade 2; these strains also are in a single somatic compatibility group. However, in the accessory genomes, the Foci versus Foa races 3 and 4 differ in multiple contigs. Based on significantly increased expression of Foa race 4 genes in planta vs. in vitro, we identified 23 putative effectors and 13 possible pathogenicity factors. PCR primers for diagnosis of either Foa race 2 or 4 and the Foci were identified. Finally, mixtures of conidia that were pre-stained with different fluorochromes indicated that Foa race 4 formed conidial anastomosis tubes (CATs) with Foci. Foa race 4 and Foa race 2, which are in different somatic compatibility groups, did not form CATs with each other.

CONCLUSIONS

There was no evidence that Foa race 2 was involved in the recent evolution of Foa race 4; Foa race 2 and 4 are CAT-incompatible. Although Foa races 3 and 4 and the Foci are closely related, there is no evidence that either Foci contributed to the evolution of Foa race 4, or that Foa race 4 was the recent recipient of a multi-gene chromosomal segment from another strain. However, horizontal chromosome transfer could account for the major difference in the accessory genomes of Foa race 4 and the Foci and for their differences in host range.

First Record of the Complete Mitochondrial Genome of a Saprotrophic and Opportunistic Human Pathogenic Fungus, Scopulariopsis brevicaulis

Scopulariopsis brevicaulis is a widely distributed soil fungus known as a common saprotroph of biodegradation. It is also an opportunistic human pathogen that can produce various secondary metabolites. Here, we report the first complete mitochondrial genome sequence of S. brevicaulis isolated from air in South Korea. Total length of the mitochondrial genome is 28,829 bp and encoded 42 genes (15 protein-coding genes, 2 rRNAs, and 25 tRNAs). Nucleotide sequence of coding region takes over 26.2%, and overall GC content is 27.6%. Phylogenetic trees present that S. brevicaulis is clustered with Lomentospora prolificans with presenting various mitochondrial genome length.

Detecting Introgression Between Members of the Fusarium fujikuroi and F. oxysporum Species Complexes by Comparative Mitogenomics

The Fusarium fujikuroi species complex (FFSC) and F. oxysporum species complex (FOSC) are two related groups of plant pathogens causing a wide diversity of diseases in agricultural crops world wide. The aims of this study are (1) to clarify the phylogeny of the FFSC, (2) to identify potential deviation from tree-like evolution, (3) to explore the value of using mitogenomes for these kinds of analyses, and (4) to better understand mitogenome evolution. In total, we have sequenced 24 species from the FFSC and a representative set of recently analyzed FOSC strains was chosen, while F. redolens was used as outgroup for the two species complexes. A species tree was constructed based on the concatenated alignment of seven nuclear genes and the mitogenome, which was contrasted to individual gene trees to identify potential conflicts. These comparisons indicated conflicts especially within the previously described African clade of the FFSC. Furthermore, the analysis of the mitogenomes revealed the presence of a variant of the large variable (LV) region in FFSC which was previously only reported for FOSC. The distribution of this variant and the results of sequence comparisons indicate horizontal genetic transfer between members of the two species complexes, most probably through introgression. In addition, a duplication of atp9 was found inside an intron of cob, which suggests that even highly conserved mitochondrial genes can have paralogs. Paralogization in turn may lead to inaccurate single gene phylogenies. In conclusion, mitochondrial genomes provide a robust basis for phylogeny. Comparative phylogenetic analysis indicated that gene flow among and between members of FFSC and FOSC has played an important role in the evolutionary history of these two groups. Since mitogenomes show greater levels of conservation and synteny than nuclear regions, they are more likely to be compatible for recombination than nuclear regions. Therefore, mitogenomes can be used as indicators to detect interspecies gene flow.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

Diversity of Mobile Genetic Elements in the Mitogenomes of Closely Related Fusarium culmorum and F. graminearum sensu stricto Strains and Its Implication for Diagnostic Purposes

Much of the mitogenome variation observed in fungal lineages seems driven by mobile genetic elements (MGEs), which have invaded their genomes throughout evolution. The variation in the distribution and nucleotide diversity of these elements appears to be the main distinction between different fungal taxa, making them promising candidates for diagnostic purposes. Fungi of the genus Fusarium display a high variation in MGE content, from MGE-poor (Fusarium oxysporum and Fusarium fujikuroi species complex) to MGE-rich mitogenomes found in the important cereal pathogens F. culmorum and F. graminearum sensu stricto. In this study, we investigated the MGE variation in these latter two species by mitogenome analysis of geographically diverse strains. In addition, a smaller set of F. cerealis and F. pseudograminearum strains was included for comparison. Forty-seven introns harboring from 0 to 3 endonucleases (HEGs) were identified in the standard set of mitochondrial protein-coding genes. Most of them belonged to the group I intron family and harbored either LAGLIDADG or GIY-YIG HEGs. Among a total of 53 HEGs, 27 were shared by all fungal strains. Most of the optional HEGs were irregularly distributed among fungal strains/species indicating ancestral mosaicism in MGEs. However, among optional MGEs, one exhibited species-specific conservation in F. culmorum. While in F. graminearum s.s. MGE patterns in cox3 and in the intergenic spacer between cox2 and nad4L may facilitate the identification of this species. Thus, our results demonstrate distinctive traits of mitogenomes for diagnostic purposes of Fusaria.

Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats

Background

Increased contamination of European and Asian wheat and barley crops with “emerging” mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises.

Results

The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)_n or (T)_n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division.

Conclusions

Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.

Population Genomic Analysis Reveals a Highly Conserved Mitochondrial Genome in Fusarium asiaticum

Fusarium asiaticum is one of the pivotal members of the Fusarium graminearum species complex (FGSC) causing Fusarium head blight (FHB) on wheat, barley and rice in large parts of Asia. Besides resulting in yield losses, FHB also causes the accumulation of mycotoxins such as nivalenol (NIV) and deoxynivalenol (DON). The aim of this study was to conduct population studies on F. asiaticum from Southern China through mitochondrial genome analyses. All strains were isolated from wheat or rice from several geographic areas in seven provinces in Southern China. Based on geographic location and host, 210 isolates were selected for next generation sequencing, and their mitogenomes were assembled by GRAbB and annotated to explore the mitochondrial genome variability of F. asiaticum. The F. asiaticum mitogenome proves extremely conserved and variation is mainly caused by absence/presence of introns harboring homing endonuclease genes. These variations could be utilized to develop molecular markers for track and trace of migrations within and between populations. This study illustrates how mitochondrial introns can be used as markers for population genetic analysis. SNP analysis demonstrate the occurrence of mitochondrial recombination in F. asiaticum as was previously found for F. oxysporum and implied for F. graminearum. Furthermore, varying degrees of genetic diversity and recombination showed a high association with different geographic regions as well as with cropping systems. The mitogenome of F. graminearum showed a much higher SNP diversity while the interspecies intron variation showed no evidence of gene flow between the two closely related and sexual compatible species.

Promising Perspectives for Detection, Identification, and Quantification of Plant Pathogenic Fungi and Oomycetes through Targeting Mitochondrial DNA

Fungi and oomycetes encompass many pathogens affecting crops worldwide. Their effective control requires screening pathogens across the local and international trade networks along with the monitoring of pathogen inocula in the field. Fundamentals to all of these concerns are their efficient detection, identification, and quantification. The use of molecular markers showed the best promise in the field of plant pathogen diagnostics. However, despite the unquestionable benefits of DNA-based methods, two significant limitations are associated with their use. The first limitation concerns the insufficient level of sensitivity due to the very low and uneven distribution of pathogens in plant material. The second limitation pertains to the inability of widely used diagnostic assays to detect cryptic species. Targeting mtDNA appears to provide a solution to these challenges. Its high copy number in microbial cells makes mtDNA an attractive target for developing highly sensitive assays. In addition, previous studies on different pathogen taxa indicated that mitogenome sequence variation could improve cryptic species delimitation accuracy. This review sheds light on the potential application of mtDNA for pathogen diagnostics. This paper covers a brief description of qPCR and DNA barcoding as two major strategies enabling the diagnostics of plant pathogenic fungi and oomycetes. Both strategies are discussed along with the potential use of mtDNA, including their strengths and weaknesses.

Phylogenetic relationship between Australian Fusarium oxysporum isolates and resolving the species complex using the multispecies coalescent model

BACKGROUND

The Fusarium oxysporum species complex (FOSC) is a ubiquitous group of fungal species readily isolated from agroecosystem and natural ecosystem soils which includes important plant and human pathogens. Genetic relatedness within the complex has been studied by sequencing either the genes or the barcoding gene regions within those genes. Phylogenetic analyses have demonstrated a great deal of diversity which is reflected in the differing number of clades identified: three, five and eight. Genetic limitation within the species in the complex has been studied through Genealogical Concordance Phylogenetic Species Recognition (GCPSR) analyses with varying number of phylogenetic 'species' identified ranging from two to 21. Such differing views have continued to confuse users of these taxonomies.

RESULTS

The phylogenetic relationships between Australian F. oxysporum isolates from both natural and agricultural ecosystems were determined using three datasets: whole genome, nuclear genes, and mitochondrial genome sequences. The phylogenies were concordant except for three isolates. There were three concordant clades from all the phylogenies suggesting similar evolutionary history for mitochondrial genome and nuclear genes for the isolates in these three clades. Applying a multispecies coalescent (MSC) model on the eight single copy nuclear protein coding genes from the nuclear gene dataset concluded that the three concordant clades correspond to three phylogenetic species within the FOSC. There was 100% posterior probability support for the formation of three species within the FOSC. This is the first report of using the MSC model to estimate species within the F. oxysporum species complex. The findings from this study were compared with previously published phylogenetics and species delimitation studies.

CONCLUSION

Phylogenetic analyses using three different gene datasets from Australian F. oxysporum isolates have all supported the formation of three major clades which delineated into three species. Species 2 (Clade 3) may be called F. oxysporum as it contains the neotype for F. oxysporum.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

Complete mitochondrial genome sequence of lettuce pathogenic fungus, Fusarium oxysporum f. sp. lactucae 09-002

Fusarium oxysporum is a famous plant pathogenic filamentous fungus. Here, we report the complete mitochondrial genome sequence of F. oxysporum f. sp. lactucae isolated from the lettuce field in Suwon area, Korea. Total length of the mitochondrial genome is 45,020 bp and it encodes 42 genes (15 protein-coding genes, two rRNAs, and 25 tRNAs). Nucleotide sequence of coding region takes over 32.7%, and overall GC content is 32.4%. Phylogenetic tree presented that F. oxysporum f. sp. lactucae 09-002 was clustered with Fusarium commune not like another F. oxysporum mitochondrial genomes, requiring further analyses.

Complete mitochondrial genome sequence of lettuce pathogenic fungus, Fusarium oxysporum f. sp. lactucae 16-086

Fusarium oxysporum f. sp. lactucae 016-086 is a plant pathogenic filamentous fungus isolated from wilted lettuce in Korea. We reported complete mitochondrial genome sequence of F. oxysporum f. sp. lactucae 016-086. Total length of this mitogenome is 41,826 bp and it encoded 42 genes (14 protein-coding genes, 2 rRNAs, and 26 tRNAs). Nucleotide sequence of coding region takes over 30.6%, and overall GC content is 32.5%. Phylogenetic tree of Fusarium mitochondrial genomes presented distinct clades along with nine formae speciales. This mitogenome will contribute distinguishing formae speciales of F. oxysoporum claearly with additional mitogenomes sequenced in the near future.

Phylogeny and species delimitations in the entomopathogenic genus Beauveria (Hypocreales, Ascomycota), including the description of B. peruviensis sp. nov

The genus Beauveria is considered a cosmopolitan anamorphic and teleomorphic genus of soilborne necrotrophic arthropod-pathogenic fungi that includes ecologically and economically important species. Species identification in Beauveria is difficult because of its structural simplicity and the lack of distinctive phenotypic variation. Therefore, the use of multi-locus sequence data is essential to establish robust species boundaries in addition to DNA-based species delimitation methods using genetic distance, coalescent, and genealogical concordance approaches (polyphasic approaches). In this regard, our study used multilocus phylogeny and five DNA-based methods to delimit species in Beauveria using three molecular makers. These polyphasic analyses allowed for the delimitation of 20-28 species in Beauveria, confirming cryptic diversity in five species (i.e. B. amorpha, B. bassiana, B. diapheromeriphila, and B. pseudobassiana) and supporting the description of B. peruviensis as a new taxon from northeastern Peru. The other five species were not evaluated as they did not have enough data (i.e. B. araneola, B. gryllotalpidicola, B. loeiensis, B. medogensis, and B. rudraprayagi). Our results demonstrate that the congruence among different methods in a polyphasic approach (e.g. genetic distance and coalescence methods) is more likely to show reliably supported species boundaries. Among the methods applied in this study, genetic distance, coalescent approaches, and multilocus phylogeny are crucial when establishing species boundaries in Beauveria.

QTL mapping of mycelial growth and aggressiveness to distinct hosts in Ceratocystis pathogens

Some species of Ceratocystis display strong host specificity, such as C. fimbriata sensu stricto that is restricted to sweet potato (Ipomoea batatas) as host. In contrast, the closely related C. manginecans, infects Acacia mangium and Mangifera indica but is not pathogenic to I. batatas. Despite the economic importance of these fungi, knowledge regarding the genetic factors that influence their pathogenicity and host specificity is limited. A recent inheritance study, based on an interspecific cross between C. fimbriata and C. manginecans and the resultant 70 F₁ progeny, confirmed that traits such as mycelial growth rate, spore production and aggressiveness on A. mangium and I. batatas are regulated by multiple genes. In the present study, a quantitative trait locus (QTL) analysis was performed to determine the genomic loci associated with these traits. All 70 progeny isolates were genotyped with SNP markers and a linkage map was constructed. The map contained 467 SNPs, distributed across nine linkage groups, with a total length of 1203 cm. Using the progeny genotypes and phenotypes, one QTL was identified on the linkage map for mycelial growth rate, one for aggressiveness to A. mangium and two for aggressiveness to I. batatas (P < 0.05). Two candidate genes, likely associated with mycelial growth rate, were identified in the QTL region. The three QTLs associated with aggressiveness to different hosts contained candidate genes involved in protein processing, detoxification and regions with effector genes and high transposable element density. The results provide a foundation for studies considering the function of genes regulating various quantitative traits in Ceratocystis.

Maternal mitochondrial inheritance in two Fusarium pathogens of prickly ash (Zanthoxylum bungeanum) in northern China

Mitochondrial inheritance in Fusarium zanthoxyli and F. continuum, two canker-inducing pathogens of prickly ash (Zanthoxylum bungeanum) in northern China, was investigated by genotyping ascospore progeny obtained from laboratory crosses. Polymorphic regions of the mitochondrial genomes (mitogenomes) that contained indels and single-nucleotide polymorphisms (SNPs) were identified via comparative analyses of the complete mitogenomes of the parents used in the intraspecific crosses. A reciprocal genetic cross of F. zanthoxyli NRRL 66714 × NRRL 66285, and a separate cross of F. continuum ♀ NRRL 66286 × ♂ NRRL 66218, revealed that mitochondria were only inherited from the maternal parent. In addition, the reciprocal cross demonstrated that mitochondrial inheritance is not linked to mating type. Gene order in the circular mitogenomes of the prickly ash pathogens was identical to that previously reported for other fusaria and members of the Hypocreales, except that the TRNL tRNAs were duplicated in F. zanthoxyli NRRL 66714. The genomes contained 14 polypeptide-encoding genes involved in oxidative respiration, one intron-encoded ribosomal protein (rps3) gene, two ribosomal RNA (rRNA) genes, and 26-28 tRNA genes. The F. zanthoxyli mitogenomes were 80.9 and 98.7 kb in length, whereas those of F. continuum were considerably shorter and nearly identical in length at 63.4 kb. The significant differences in mitogenome length were primarily due to variable numbers of introns and open reading frames (ORFs) encoding hypothetical proteins.

First steps towards mitochondrial pan-genomics: detailed analysis of Fusarium graminearum mitogenomes

There is a gradual shift from representing a species' genome by a single reference genome sequence to a pan-genome representation. Pan-genomes are the abstract representations of the genomes of all the strains that are present in the population or species. In this study, we employed a pan-genomic approach to analyze the intraspecific mitochondrial genome diversity of Fusarium graminearum. We present an improved reference mitochondrial genome for F. graminearum with an intron-exon annotation that was verified using RNA-seq data. Each of the 24 studied isolates had a distinct mitochondrial sequence. Length variation in the F. graminearum mitogenome was found to be largely due to variation of intron regions (99.98%). The "intronless" mitogenome length was found to be quite stable and could be informative when comparing species. The coding regions showed high conservation, while the variability of intergenic regions was highest. However, the most important variable parts are the intron regions, because they contain approximately half of the variable sites, make up more than half of the mitogenome, and show presence/absence variation. Furthermore, our analyses show that the mitogenome of F. graminearum is recombining, as was previously shown in F. oxysporum, indicating that mitogenome recombination is a common phenomenon in Fusarium. The majority of mitochondrial introns in F. graminearum belongs to group I introns, which are associated with homing endonuclease genes (HEGs). Mitochondrial introns containing HE genes may spread within populations through homing, where the endonuclease recognizes and cleaves the recognition site in the target gene. After cleavage of the "host" gene, it is replaced by the gene copy containing the intron with HEG. We propose to use introns unique to a population for tracking the spread of the given population, because introns can spread through vertical inheritance, recombination as well as via horizontal transfer. We demonstrate how pooled sequencing of strains can be used for mining mitogenome data. The usage of pooled sequencing offers a scalable solution for population analysis and for species level comparisons studies. This study may serve as a basis for future mitochondrial genome variability studies and representations.

Epitypification of Fusarium oxysporum - clearing the taxonomic chaos

Fusarium oxysporum is the most economically important and commonly encountered species of Fusarium. This soil-borne fungus is known to harbour both pathogenic (plant, animal and human) and non-pathogenic strains. However, in its current concept F. oxysporum is a species complex consisting of numerous cryptic species. Identification and naming these cryptic species is complicated by multiple subspecific classification systems and the lack of living ex-type material to serve as basic reference point for phylogenetic inference. Therefore, to advance and stabilise the taxonomic position of F. oxysporum as a species and allow naming of the multiple cryptic species recognised in this species complex, an epitype is designated for F. oxysporum. Using multi-locus phylogenetic inference and subtle morphological differences with the newly established epitype of F. oxysporum as reference point, 15 cryptic taxa are resolved in this study and described as species.

The first complete mitochondrial genome from the family Hygrophoraceae (Hygrophorus russula) by next-generation sequencing and phylogenetic implications

Hygrophorus russula (Schaeff.) Kauffman is an edible ectomycorrhizal fungus that is widely distributed in the world. In this study, the mitogenome of H. russula was sequenced and assembled. The mitogenome of H. russula is composed of circular DNA molecules, with a total size of 55,769 bp. Further analysis indicated that the frequent use of A and T in codons contributes to the high AT content (80.87%) in the H. russula mitogenome. Comparative analysis indicated that the length and base composition of the core protein-encoding genes, and the number of tRNA genes in the H. russula mitogenome varied from that of other Agaricales mitogenomes. Gene arrangement analysis revealed a novel gene order in the H. russula mitogenome. In addition, the expansion of the mitogenome in Agaricales was found to be closely related to the increase in the number of introns. Phylogenetic analysis of the combined mitochondrial gene set showed strong support for tree topologies, and H. russula was determined to be relatively distant from other Agaricales species. This study is the first report on the mitogenome of a member of genus Hygrophorus as well as family Hygrophoraceae, which improves our understanding of mitochondrial differentiation and evolution in the important ectomycorrhizal fungi Hygrophorus species.

Host-Pathogen Interaction in Fusarium oxysporum Infections: Where Do We Stand?

Fusarium oxysporum, a ubiquitous soilborne pathogen, causes devastating vascular wilt in more than 100 plant species and ranks 5th among the top 10 fungal plant pathogens. It has emerged as a human pathogen, too, causing infections in immune-compromised patients. Therefore, it is important to gain insight into the molecular processes involved in the pathogenesis of this transkingdom pathogen. A complex network comprising interconnected and overlapping signal pathways-mitogen-activated protein kinase signaling pathways, Ras proteins, G-protein signaling components and their downstream pathways, components of the velvet (LaeA/VeA/VelB) complex, and cAMP pathways-is involved in perceiving the host. This network regulates the expression of various pathogenicity genes. However, plants have evolved an elaborate protection system to combat this attack. They, too, possess intricate mechanisms at the molecular level which, once triggered by pathogen attack, transduce signals to activate defense response. This review focuses on understanding and presenting a wholistic picture of the molecular mechanisms of F. oxysporum-host interactions in plant immunity.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Updates in the Language of Histoplasma Biodiversity

In a recent article, Sepúlveda et al. (mBio 8:e01339-17, 2017, https://doi.org/10.1128/mBio.01339-17) investigated the genetic structure and evolutionary history of the human pathogen Histoplasma. Using whole-genome resequencing data, Sepúlveda et al. found that the Histoplasma genus is composed of at least four strongly differentiated lineages. Their tour de force is to use a smart combination of population genomic approaches to show that the advanced stage of intraspecific divergence observed within Histoplasma does not simply reflect population structure, but instead results from previously unidentified speciation events. The four independently evolving Histoplasma lineages are elevated to the species status and assigned names. The newly described species exhibit medically important differences in phenotype, and these findings, therefore, have important epidemiological implications. This work provides a blueprint for phylogenomic species recognition in fungi, opening the way for a new age of enlightenment in which fungal species are diagnosed using highly discriminatory tools within a hypothesis-testing framework.

Mitochondrial Recombination Reveals Mito-Mito Epistasis in Yeast

Genetic variation in mitochondrial DNA (mtDNA) provides adaptive potential although the underlying genetic architecture of fitness components within mtDNAs is not known. To dissect functional variation within mtDNAs, we first identified naturally occurring mtDNAs that conferred high or low fitness in Saccharomyces cerevisiae by comparing growth in strains containing identical nuclear genotypes but different mtDNAs. During respiratory growth under temperature and oxidative stress conditions, mitotype effects were largely independent of nuclear genotypes even in the presence of mito-nuclear interactions. Recombinant mtDNAs were generated to determine fitness components within high- and low-fitness mtDNAs. Based on phenotypic distributions of isogenic strains containing recombinant mtDNAs, we found that multiple loci contributed to mitotype fitness differences. These mitochondrial loci interacted in epistatic, nonadditive ways in certain environmental conditions. Mito-mito epistasis (i.e., nonadditive interactions between mitochondrial loci) influenced fitness in progeny from four different crosses, suggesting that mito-mito epistasis is a widespread phenomenon in yeast and other systems with recombining mtDNAs. Furthermore, we found that interruption of coadapted mito-mito interactions produced recombinant mtDNAs with lower fitness. Our results demonstrate that mito-mito epistasis results in functional variation through mitochondrial recombination in fungi, providing modes for adaptive evolution and the generation of mito-mito incompatibilities.

Gene Flow between Divergent Cereal- and Grass-Specific Lineages of the Rice Blast Fungus Magnaporthe oryzae

Delineating species and epidemic lineages in fungal plant pathogens is critical to our understanding of disease emergence and the structure of fungal biodiversity and also informs international regulatory decisions. Pyricularia oryzae (syn. Magnaporthe oryzae) is a multihost pathogen that infects multiple grasses and cereals, is responsible for the most damaging rice disease (rice blast), and is of growing concern due to the recent introduction of wheat blast to Bangladesh from South America. However, the genetic structure and evolutionary history of M. oryzae, including the possible existence of cryptic phylogenetic species, remain poorly defined. Here, we use whole-genome sequence information for 76 M. oryzae isolates sampled from 12 grass and cereal genera to infer the population structure of M. oryzae and to reassess the species status of wheat-infecting populations of the fungus. Species recognition based on genealogical concordance, using published data or extracting previously used loci from genome assemblies, failed to confirm a prior assignment of wheat blast isolates to a new species (Pyricularia graminis-tritici). Inference of population subdivisions revealed multiple divergent lineages within M. oryzae, each preferentially associated with one host genus, suggesting incipient speciation following host shift or host range expansion. Analyses of gene flow, taking into account the possibility of incomplete lineage sorting, revealed that genetic exchanges have contributed to the makeup of multiple lineages within M. oryzae These findings provide greater understanding of the ecoevolutionary factors that underlie the diversification of M. oryzae and highlight the practicality of genomic data for epidemiological surveillance in this important multihost pathogen.IMPORTANCE Infection of novel hosts is a major route for disease emergence by pathogenic microorganisms. Understanding the evolutionary history of multihost pathogens is therefore important to better predict the likely spread and emergence of new diseases. Magnaporthe oryzae is a multihost fungus that causes serious cereal diseases, including the devastating rice blast disease and wheat blast, a cause of growing concern due to its recent spread from South America to Asia. Using whole-genome analysis of 76 fungal strains from different hosts, we have documented the divergence of M. oryzae into numerous lineages, each infecting a limited number of host species. Our analyses provide evidence that interlineage gene flow has contributed to the genetic makeup of multiple M. oryzae lineages within the same species. Plant health surveillance is therefore warranted to safeguard against disease emergence in regions where multiple lineages of the fungus are in contact with one another.

Gene Flow between Divergent Cereal- and Grass-Specific Lineages of the Rice Blast Fungus Magnaporthe oryzae

Delineating species and epidemic lineages in fungal plant pathogens is critical to our understanding of disease emergence and the structure of fungal biodiversity and also informs international regulatory decisions. Pyricularia oryzae (syn. Magnaporthe oryzae) is a multihost pathogen that infects multiple grasses and cereals, is responsible for the most damaging rice disease (rice blast), and is of growing concern due to the recent introduction of wheat blast to Bangladesh from South America. However, the genetic structure and evolutionary history of M. oryzae, including the possible existence of cryptic phylogenetic species, remain poorly defined. Here, we use whole-genome sequence information for 76 M. oryzae isolates sampled from 12 grass and cereal genera to infer the population structure of M. oryzae and to reassess the species status of wheat-infecting populations of the fungus. Species recognition based on genealogical concordance, using published data or extracting previously used loci from genome assemblies, failed to confirm a prior assignment of wheat blast isolates to a new species (Pyricularia graminis-tritici). Inference of population subdivisions revealed multiple divergent lineages within M. oryzae, each preferentially associated with one host genus, suggesting incipient speciation following host shift or host range expansion. Analyses of gene flow, taking into account the possibility of incomplete lineage sorting, revealed that genetic exchanges have contributed to the makeup of multiple lineages within M. oryzae These findings provide greater understanding of the ecoevolutionary factors that underlie the diversification of M. oryzae and highlight the practicality of genomic data for epidemiological surveillance in this important multihost pathogen.IMPORTANCE Infection of novel hosts is a major route for disease emergence by pathogenic microorganisms. Understanding the evolutionary history of multihost pathogens is therefore important to better predict the likely spread and emergence of new diseases. Magnaporthe oryzae is a multihost fungus that causes serious cereal diseases, including the devastating rice blast disease and wheat blast, a cause of growing concern due to its recent spread from South America to Asia. Using whole-genome analysis of 76 fungal strains from different hosts, we have documented the divergence of M. oryzae into numerous lineages, each infecting a limited number of host species. Our analyses provide evidence that interlineage gene flow has contributed to the genetic makeup of multiple M. oryzae lineages within the same species. Plant health surveillance is therefore warranted to safeguard against disease emergence in regions where multiple lineages of the fungus are in contact with one another.