Phylogenomic Analysis of a 55.1-kb 19-Gene Dataset Resolves a Monophyletic Fusarium that Includes the Fusarium solani Species Complex

Scientific communication is facilitated by a data-driven, scientifically sound taxonomy that considers the end-user's needs and established successful practice. In 2013, the Fusarium community voiced near unanimous support for a concept of Fusarium that represented a clade comprising all agriculturally and clinically important Fusarium species, including the F. solani species complex (FSSC). Subsequently, this concept was challenged in 2015 by one research group who proposed dividing the genus Fusarium into seven genera, including the FSSC described as members of the genus Neocosmospora, with subsequent justification in 2018 based on claims that the 2013 concept of Fusarium is polyphyletic. Here, we test this claim and provide a phylogeny based on exonic nucleotide sequences of 19 orthologous protein-coding genes that strongly support the monophyly of Fusarium including the FSSC. We reassert the practical and scientific argument in support of a genus Fusarium that includes the FSSC and several other basal lineages, consistent with the longstanding use of this name among plant pathologists, medical mycologists, quarantine officials, regulatory agencies, students, and researchers with a stake in its taxonomy. In recognition of this monophyly, 40 species described as genus Neocosmospora were recombined in genus Fusarium, and nine others were renamed Fusarium. Here the global Fusarium community voices strong support for the inclusion of the FSSC in Fusarium, as it remains the best scientific, nomenclatural, and practical taxonomic option available.

Genetic Differentiation of Verticillium dahliae Populations Recovered from Symptomatic and Asymptomatic Hosts

Verticillium dahliae is a soilborne fungal pathogen affecting many economically important crops that can also infect weeds and rotational crops with no apparent disease symptoms. The main research goal was to test the hypothesis that V. dahliae populations recovered from asymptomatic rotational crops and weed species are evolutionarily and genetically distinct from symptomatic hosts. We collected V. dahliae isolates from symptomatic and asymptomatic hosts growing in fields with histories of Verticillium wilt of potato in Israel and Pennsylvania (United States), and used genotyping-by-sequencing to analyze the evolutionary history and genetic differentiation between populations of different hosts. A phylogeny inferred from 26,934 single-nucleotide polymorphisms (SNPs) in 126 V. dahliae isolates displayed a highly clonal structure correlated with vegetative compatibility groups, and isolates grouped in lineages 2A, 2B⁸²⁴, 4A, and 4B, with 77% of the isolates in lineage 4B. The lineages identified in this study were differentiated by host of origin; we found 2A, 2B⁸²⁴, and 4A only in symptomatic hosts but isolates from asymptomatic hosts (weeds, oat, and sorghum) grouped exclusively within lineage 4B, and were genetically indistinguishable from 4B isolates sampled from symptomatic hosts (potato, eggplant, and avocado). Using coalescent analysis of 158 SNPs of lineage 4B, we inferred a genealogy with clades that correlated with geographic origin. In contrast, isolates from asymptomatic and symptomatic hosts shared some of the same haplotypes and were not differentiated. We conclude that asymptomatic weeds and rotational hosts may be potential reservoirs for V. dahliae populations of lineage 4B, which are pathogenic to many cultivated hosts.

Clonal Expansion and Migration of a Highly Virulent, Defoliating Lineage of Verticillium dahliae

We used a population genomics approach to test the hypothesis of clonal expansion of a highly fit genotype in populations of Verticillium dahliae. This fungal pathogen has a broad host range and can be dispersed in contaminated seed or other plant material. It has a highly clonal population structure, with several lineages having nearly worldwide distributions in agricultural crops. Isolates in lineage 1A are highly virulent and cause defoliation in cotton, okra, and olive (denoted 1A/D), whereas those in other lineages cause wilting but not defoliation (ND). We tested whether the highly virulent lineage 1A/D could have spread from the southwestern United States to the Mediterranean basin, as predicted from historical records. We found 187 single-nucleotide polymorphisms (SNPs), determined by genotyping by sequencing, among 91 isolates of lineage 1A/D and 5 isolates in the closely related lineage 1B/ND. Neighbor-joining and maximum-likelihood analyses on the 187 SNPs showed a clear divergence between 1A/D and 1B/ND haplotypes. Data for only 77 SNPs were obtained for all 96 isolates (no missing data); lineages 1A/D and 1B/ND differed by 27 of these 77 SNPs, confirming a clear divergence between the two lineages. No evidence of recombination was detected within or between these two lineages. Phylogenetic and genealogical analyses resulted in five distinct subclades of 1A/D isolates that correlated closely with geographic origins in the Mediterranean basin, consistent with the hypothesis that the D pathotype was introduced at least five times in independent founder events into this region from a relatively diverse source population. The inferred ancestral haplotype was found in two isolates sampled before 1983 from the southwestern United States, which is consistent with historical records that 1A/D originated in North America. The five subclades coalesce with the ancestral haplotype at the same time, consistent with a hypothesis of rapid population expansion in the source population during the emergence of 1A/D as a severe pathogen of cotton in the United States.

Highly diverse endophytic and soil Fusarium oxysporum populations associated with field-grown tomato plants

The diversity and genetic differentiation of populations of Fusarium oxysporum associated with tomato fields, both endophytes obtained from tomato plants and isolates obtained from soil surrounding the sampled plants, were investigated. A total of 609 isolates of F. oxysporum were obtained, 295 isolates from a total of 32 asymptomatic tomato plants in two fields and 314 isolates from eight soil cores sampled from the area surrounding the plants. Included in this total were 112 isolates from the stems of all 32 plants, a niche that has not been previously included in F. oxysporum population genetics studies. Isolates were characterized using the DNA sequence of the translation elongation factor 1α gene. A diverse population of 26 sequence types was found, although two sequence types represented nearly two-thirds of the isolates studied. The sequence types were placed in different phylogenetic clades within F. oxysporum, and endophytic isolates were not monophyletic. Multiple sequence types were found in all plants, with an average of 4.2 per plant. The population compositions differed between the two fields but not between soil samples within each field. A certain degree of differentiation was observed between populations associated with different tomato cultivars, suggesting that the host genotype may affect the composition of plant-associated F. oxysporum populations. No clear patterns of genetic differentiation were observed between endophyte populations and soil populations, suggesting a lack of specialization of endophytic isolates.

Complex molecular relationship between vegetative compatibility groups (VCGs) in Verticillium dahliae: VCGs do not always align with clonal lineages

Verticillium wilts caused by the soilborne fungus Verticillium dahliae are among the most challenging diseases to control. Populations of this pathogen have been traditionally studied by means of vegetative compatibility groups (VCGs) under the assumption that VCGs comprise genetically related isolates that correlate with clonal lineages. We aimed to resolve the phylogenetic relationships among VCGs and their subgroups based on sequences of the intergenic spacer region (IGS) of the ribosomal DNA and six anonymous polymorphic sequences containing single-nucleotide polymorphisms (VdSNPs). A collection of 68 V. dahliae isolates representing the main VCGs and subgroups (VCGs 1A, 1B, 2A, 2B, 3, 4A, 4B, and 6) from different geographic origins and hosts was analyzed using the seven DNA regions. Maximum parsimony (MP) phylogenies inferred from IGS and VdSNP sequences showed five and six distinct clades, respectively. Phylogenetic analyses of individual and combined data sets indicated that certain VCG subgroups (e.g., VCGs 1A and 1B) are closely related and share a common ancestor; however, other subgroups (e.g., VCG 4B) are more closely related to members of a different VCG (e.g., VCG 2A) than to subgroups of the same VCG (VCG 4B). Furthermore, MP analyses indicated that VCG 2B is polyphyletic, with isolates placed in at least three distinct phylogenetic lineages based on IGS sequences and two lineages based on VdSNP sequences. Results from our study suggest the existence of main VCG lineages that contain VCGs 1A and 1B; VCGs 2A and 4B; and VCG 4A, for which both phylogenies agree; and the existence of other VCGs or VCG subgroups that seem to be genetically heterogeneous or show discrepancies in their phylogenetic placement: VCG 2B, VCG 3, and VCG 6. These results raise important caveats regarding the interpretation of VCG analyses: genetic homogeneity and close evolutionary relationship between members of a VCG should not be assumed.

Hidden host plant associations of soilborne fungal pathogens: an ecological perspective

Much of the current knowledge on population biology and ecology of soilborne fungal pathogens has been derived from research based on populations recovered from plants displaying disease symptoms or soil associated with symptomatic plants. Many soilborne fungal pathogens are known to cause disease on a large number of crop plants, including a variety of important agronomical, horticultural, ornamental, and forest plants species. For instance, the fungus Verticillium dahliae causes disease on >400 host plants. From a phytopathological perspective, plants on which disease symptoms have not been yet observed are considered to be nonhosts for V. dahliae. This term may be misleading because it does not provide information regarding the nature of the plant-fungus association; that is, a nonhost plant may harbor the fungus as an endophyte. Yet, there are numerous instances in the literature where V. dahliae has been isolated from asymptomatic plants; thus, these plants should be considered hosts. In this article, we synthesize scattered research that indicates that V. dahliae, aside from being a successful and significant vascular plant pathogen, may have a cryptic biology on numerous asymptomatic plants as an endophyte. Thus, we suggest here that these endophytic associations among V. dahliae and asymptomatic plants are not unusual relationships in nature. We propose to embrace the broader ecology of many fungi by differentiating between "symptomatic hosts" as those plants in which the infection and colonization by a fungus results in disease, and "asymptomatic hosts" as those plants that harbor the fungus endophytically and are different than true nonhosts that should be used for plant species that do not interact with the given fungus. In fact, if we broaden our definition of "host plant" to include asymptomatic plants that harbor the fungus as an endophyte, it is likely that the host ranges for some soilborne fungal pathogens are much larger than previously envisioned. By ignoring the potential for soilborne fungal pathogens to display endophytic relationships, we leave gaps in our knowledge about the population biology and ecology, persistence, and spread of these fungi in agroecosystems.

Sequence variation in two protein-coding genes correlates with mycelial compatibility groupings in Sclerotium rolfsii

Populations of Sclerotium rolfsii, the causal organism of Sclerotium root-rot on a wide range of hosts, can be placed into mycelial compatibility groups (MCGs). In this study, we evaluated three different molecular approaches to unequivocally identify each of 12 previously identified MCGs. These included restriction fragment length polymorphism (RFLP) patterns of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA (rDNA) and sequence analysis of two protein-coding genes: translation elongation factor 1α (EF1α) and RNA polymerase II subunit two (RPB2). A collection of 238 single-sclerotial isolates representing 12 MCGs of S. rolfsii were obtained from diseased sugar beet plants from Chile, Italy, Portugal, and Spain. ITS-RFLP analysis using four restriction enzymes (AluI, HpaII, RsaI, and MboI) displayed a low degree of variability among MCGs. Only three different restriction profiles were identified among S. rolfsii isolates, with no correlation to MCG or to geographic origin. Based on nucleotide polymorphisms, the RPB2 gene was more variable among MCGs compared with the EF1α gene. Thus, 10 of 12 MCGs could be characterized utilizing the RPB2 region only, while the EF1α region resolved 7 MCGs. However, the analysis of combined partial sequences of EF1α and RPB2 genes allowed discrimination among each of the 12 MCGs. All isolates belonging to the same MCG showed identical nucleotide sequences that differed by at least in one nucleotide from a different MCG. The consistency of our results to identify the MCG of a given S. rolfsii isolate using the combined sequences of EF1α and RPB2 genes was confirmed using blind trials. Our study demonstrates that sequence variation in the protein-coding genes EF1α and RPB2 may be exploited as a diagnostic tool for MCG typing in S. rolfsii as well as to identify previously undescribed MCGs.

Sublethal Doses of Mefenoxam Enhance Pythium Damping-off of Geranium

The effect of sublethal doses of fungicides on fungicide-resistant Pythium isolates is unknown but potentially relevant to disease management. Occasional grower reports of Pythium disease increases after fungicide applications and our observations of greater radial growth in vitro on fungicide-amended media than on nonamended media suggests that Pythium isolates may be stimulated by sublethal doses of fungicides. The objectives of this study were to determine whether Pythium isolates were stimulated by sublethal doses of mefenoxam in vitro and whether this stimulation had any influence on Pythium damping-off of geranium seedlings. A mefenoxam-resistant isolate of Pythium aphanidermatum displayed 10% mean radial growth increase in vitro with mefenoxam at 1 × 10^-10 μg/ml compared with growth on nonamended agar (nonsignificant). Geranium seedlings treated with one of eight mefenoxam concentrations were inoculated with 5-mm-diameter colonized agar plugs and evaluated for disease severity every 24 h. The area under the disease progress curve and the survival curve were estimated for each treatment and compared. Significant increases in damping-off were observed with mefenoxam at 1 × 10^-6 and 1 × 10^-10 μg/ml. Our data indicate that a Pythium isolate with resistance to mefenoxam can be stimulated by sublethal doses of the fungicide, and that this stimulation can result in significantly higher rates of Pythium damping-off of geranium seedlings.

Verticillium Wilt: A Threat to Artichoke Production

Verticillium wilt is becoming an increasing concern in artichoke production because the rapid spread of the disease to new growing areas has led to declining production. Scientists from Italy, Spain, and the United States combine to bring us up to date on diagnosis of the disease, its epidemiology and life cycle, as well as management strategies, current and forthcoming.

Temperature Response of Chickpea Cultivars to Races of Fusarium oxysporum f. sp. ciceris, Causal Agent of Fusarium Wilt

Use of resistant cultivars and adjustment of sowing dates are important measures for management of Fusarium wilt in chickpeas (Cicer arietinum). In this study, we examined the effect of temperature on resistance of chickpea cultivars to Fusarium wilt caused by various races of Fusarium oxysporum f. sp. ciceris. Greenhouse experiments indicated that the chickpea cultivar Ayala was moderately resistant to F. oxysporum f. sp. ciceris when inoculated plants were maintained at a day/night temperature regime of 24/21°C but was highly susceptible to the pathogen at 27/25°C. Field experiments in Israel over three consecutive years indicated that the high level of resistance of Ayala to Fusarium wilt when sown in mid- to late January differed from a moderately susceptible reaction under warmer temperatures when sowing was delayed to late February or early March. Experiments in growth chambers showed that a temperature increase of 3°C from 24 to 27°C was sufficient for the resistance reaction of cultivars Ayala and PV-1 to race 1A of the pathogen to shift from moderately or highly resistant at constant 24°C to highly susceptible at 27°C. A similar but less pronounced effect was found when Ayala plants were inoculated with F. oxysporum f. sp. ciceris race 6. Conversely, the reaction of cultivar JG-62 to races 1A and 6 was not influenced by temperature, but less disease developed on JG-62 plants inoculated with a variant of race 5 of F. oxysporum f. sp. ciceris at 27°C compared with plants inoculated at 24°C. These results indicate the importance of appropriate adjustment of temperature in tests for characterizing the resistance reactions of chickpea cultivars to the pathogen, as well as when determining the races of isolates of F. oxysporum f. sp. ciceris. Results from this study may influence choice of sowing date and use of chickpea cultivars for management of Fusarium wilt of chickpea.

Stepwise Evolution of Races in Fusarium oxysporum f. sp. ciceris Inferred from Fingerprinting with Repetitive DNA Sequences

ABSTRACT Plant pathogens often exhibit variation in virulence, the ability to cause disease on host plants with specific resistance, evident from the diversity of races observed within pathogen species. The evolution of races in asexual fungal pathogens has been hypothesized to occur in a stepwise fashion, in which mutations to virulence accumulate sequentially in clonal lineages, resulting in races capable of overcoming multiple host plant resistance genes or multiple resistant cultivars. In this study, we demonstrate a simple stepwise pattern of race evolution in Fusarium oxysporum f. sp. ciceris, the fungus that causes Fusarium wilt of chickpeas. The inferred intraspecific phylogeny of races in this fungus, based on DNA fingerprinting with repetitive sequences, shows that each of the eight races forms a monophyletic lineage. By mapping virulence to each differential cultivar (used for defining races) onto the inferred phylogeny, we show that virulence has been acquired in a simple stepwise pattern, with few parallel gains or losses. Such a clear pattern of stepwise evolution of races, to our knowledge, has not been demonstrated previously for other pathogens based on analyses of field populations. We speculate that in other systems the stepwise pattern is obscured by parallel gains or losses of virulence caused by higher mutation rates and selection by widespread deployment of resistant cultivars. Although chickpea cultivars resistant to Fusarium wilt are available, their deployment has not been extensive and the stepwise acquisition of virulence is still clearly evident.

Development of a Specific Polymerase Chain Reaction-Based Assay for the Identification of Fusarium oxysporum f. sp. ciceris and Its Pathogenic Races 0, 1A, 5, and 6

ABSTRACT Specific primers and polymerase chain reaction (PCR) assays that identify Fusarium oxysporum f. sp. ciceris and each of the F. oxysporum f. sp. ciceris pathogenic races 0, 1A, 5, and 6 were developed. F. oxysporum f. sp. ciceris- and race-specific random amplified polymorphic DNA (RAPD) markers identified in a previous study were cloned and sequenced, and sequence characterized amplified region (SCAR) primers for specific PCR were developed. Each cloned RAPD marker was characterized by Southern hybridization analysis of Eco RI-digested genomic DNA of a subset of F. oxysporum f. sp. ciceris and nonpathogenic F. oxysporum isolates. All except two cloned RAPD markers consisted of DNA sequences that were found highly repetitive in the genome of all F. oxysporum f. sp. ciceris races. F. oxysporum f. sp. ciceris isolates representing eight reported races from a wide geographic range, nonpathogenic F. oxysporum isolates, isolates of F. oxysporum f. spp. lycopersici, melonis, niveum, phaseoli, and pisi, and isolates of 47 different Fusarium spp. were tested using the SCAR markers developed. The specific primer pairs amplified a single 1,503-bp product from all F. oxysporum f. sp. ciceris isolates; and single 900- and 1,000-bp products were selectively amplified from race 0 and race 6 isolates, respectively. The specificity of these amplifications was confirmed by hybridization analysis of the PCR products. A race 5-specific identification assay was developed using a touchdown-PCR procedure. A joint use of race 0- and race 6-specific SCAR primers in a single-PCR reaction together with a PCR assay using the race 6-specific primer pair correctly identified race 1A isolates for which no RAPD marker had been found previously. All the PCR assays described herein detected up to 0.1 ng of fungal genomic DNA. The specific SCAR primers and PCR assays developed in this study clearly identify and differentiate isolates of F. oxysporum f. sp. ciceris and of each of its pathogenic races 0, 1A, 5, and 6.