Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

Attenuation of phytofungal pathogenicity of Ascomycota by autophagy modulators

Autophagy in eukaryotes functions to maintain homeostasis by degradation and recycling of long-lived and unwanted cellular materials. Autophagy plays important roles in pathogenicity of various fungal pathogens, suggesting that autophagy is a novel target for development of antifungal compounds. Here, we describe bioluminescence resonance energy transfer (BRET)-based high-throughput screening (HTS) strategy to identify compounds that inhibit fungal ATG4 cysteine protease-mediated cleavage of ATG8 that is critical for autophagosome formation. We identified ebselen (EB) and its analogs ebselen oxide (EO) and 2-(4-methylphenyl)-1,2-benzisothiazol-3(2H)-one (PT) as inhibitors of fungal pathogens Botrytis cinerea and Magnaporthe oryzae ATG4-mediated ATG8 processing. The EB and its analogs inhibit spore germination, hyphal development, and appressorium formation in Ascomycota pathogens, B. cinerea, M. oryzae, Sclerotinia sclerotiorum and Monilinia fructicola. Treatment with EB and its analogs significantly reduced fungal pathogenicity. Our findings provide molecular insights to develop the next generation of antifungal compounds by targeting autophagy in important fungal pathogens.

Genetic and molecular landscapes of the generalist phytopathogen Botrytis cinerea

Botrytis cinerea Pers. Fr. (teleomorph: Botryotinia fuckeliana) is a necrotrophic fungal pathogen that attacks a wide range of plants. This updated pathogen profile explores the extensive genetic diversity of B. cinerea, highlights the progress in genome sequencing, and provides current knowledge of genetic and molecular mechanisms employed by the fungus to attack its hosts. In addition, we also discuss recent innovative strategies to combat B. cinerea.

TAXONOMY

Kingdom: Fungi, phylum: Ascomycota, subphylum: Pezizomycotina, class: Leotiomycetes, order: Helotiales, family: Sclerotiniaceae, genus: Botrytis, species: cinerea.

HOST RANGE

B. cinerea infects almost all of the plant groups (angiosperms, gymnosperms, pteridophytes, and bryophytes). To date, 1606 plant species have been identified as hosts of B. cinerea.

GENETIC DIVERSITY

This polyphagous necrotroph has extensive genetic diversity at all population levels shaped by climate, geography, and plant host variation.

PATHOGENICITY

Genetic architecture of virulence and host specificity is polygenic using multiple weapons to target hosts, including secretory proteins, complex signal transduction pathways, metabolites, and mobile small RNA.

DISEASE CONTROL STRATEGIES

Efforts to control B. cinerea, being a high-diversity generalist pathogen, are complicated. However, integrated disease management strategies that combine cultural practices, chemical and biological controls, and the use of appropriate crop varieties will lessen yield losses. Recently, studies conducted worldwide have explored the potential of small RNA as an efficient and environmentally friendly approach for combating grey mould. However, additional research is necessary, especially on risk assessment and regulatory frameworks, to fully harness the potential of this technology.

The Adaptation of Botrytis cinerea Extracellular Vesicles Proteome to Surrounding Conditions: Revealing New Tools for Its Infection Process

Extracellular vesicles (EVs) are membranous particles released by different organisms. EVs carry several sets of macromolecules implicated in cell communication. EVs have become a relevant topic in the study of pathogenic fungi due to their relationship with fungal-host interactions. One of the essential research areas in this field is the characterization protein profile of EVs since plant fungal pathogens rely heavily on secreted proteins to invade their hosts. However, EVs of Botrytis cinerea are little known, which is one of the most devastating phytopathogenic fungi. The present study has two main objectives: the characterization of B. cinerea EVs proteome changes under two pathogenic conditions and the description of their potential role during the infective process. All the experimental procedure was conducted in B. cinerea growing in a minimal salt medium supplemented with glucose as a constitutive stage and deproteinized tomato cell walls (TCW) as a virulence inductor. The isolation of EVs was performed by differential centrifugation, filtration, ultrafiltration, and sucrose cushion ultracentrifugation. EVs fractions were visualised by TEM using negative staining. Proteomic analysis of EVs cargo was addressed by LC-MS/MS. The methodology used allowed the correct isolation of B. cinerea EVs and the identification of a high number of EV proteins, including potential EV markers. The isolated EVs displayed differences in morphology under both assayed conditions. GO analysis of EV proteins showed enrichment in cell wall metabolism and proteolysis under TCW. KEGG analysis also showed the difference in EVs function under both conditions, highlighting the presence of potential virulence/pathogenic factors implicated in cell wall metabolism, among others. This work describes the first evidence of EVs protein cargo adaptation in B. cinerea, which seems to play an essential role in its infection process, sharing crucial functions with the conventional secretion pathways.

Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis-Botrytis pathosystem

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Sclerotinia sclerotiorum (Lib.) de Bary: Insights into the Pathogenomic Features of a Global Pathogen

Sclerotinia sclerotiorum (Lib.) de Bary is a broad host-range fungus that infects an inclusive array of plant species and afflicts significant yield losses globally. Despite being a notorious pathogen, it has an uncomplicated life cycle consisting of either basal infection from myceliogenically germinated sclerotia or aerial infection from ascospores of carpogenically germinated sclerotia. This fungus is unique among necrotrophic pathogens in that it inevitably colonizes aging tissues to initiate an infection, where a saprophytic stage follows the pathogenic phase. The release of cell wall-degrading enzymes, oxalic acid, and effector proteins are considered critical virulence factors necessary for the effective pathogenesis of S. sclerotiorum. Nevertheless, the molecular basis of S. sclerotiorum pathogenesis is still imprecise and remains a topic of continuing research. Previous comprehensive sequencing of the S. sclerotiorum genome has revealed new insights into its genome organization and provided a deeper comprehension of the sophisticated processes involved in its growth, development, and virulence. This review focuses on the genetic and genomic aspects of fungal biology and molecular pathogenicity to summarize current knowledge of the processes utilized by S. sclerotiorum to parasitize its hosts. Understanding the molecular mechanisms regulating the infection process of S. sclerotiorum will contribute to devising strategies for preventing infections caused by this destructive pathogen.

Intraspecific Comparative Analysis Reveals Genomic Variation of Didymella arachidicola and Pathogenicity Factors Potentially Related to Lesion Phenotype

Didymella arachidicola is one of the most important fungal pathogens, causing foliar disease and leading to severe yield losses of peanuts (Arachis hypogaea L.) in China. Two main lesion phenotypes of peanut web blotch have been identified as reticulation type (R type) and blotch type (B type). As no satisfactory reference genome is available, the genomic variations and pathogenicity factors of D. arachidicola remain to be revealed. In the present study, we collected 41 D. arachidicola isolates from 26 geographic locations across China (33 for R type and 8 for B type). The chromosome-scale genome of the most virulent isolate (YY187) was assembled as a reference using PacBio and Hi-C technologies. In addition, we re-sequenced 40 isolates from different sampling sites. Genome-wide alignments showed high similarity among the genomic sequences from the 40 isolates, with an average mapping rate of 97.38%. An average of 3242 SNPs and 315 InDels were identified in the genomic variation analysis, which revealed an intraspecific polymorphism in D. arachidicola. The comparative analysis of the most and least virulent isolates generated an integrated gene set containing 512 differential genes. Moreover, 225 genes individually or simultaneously harbored hits in CAZy-base, PHI-base, DFVF, etc. Compared with the R type reference, the differential gene sets from all B type isolates identified 13 shared genes potentially related to lesion phenotype. Our results reveal the intraspecific genomic variation of D. arachidicola isolates and pathogenicity factors potentially related to different lesion phenotypes. This work sets a genomic foundation for understanding the mechanisms behind genomic diversity driving different pathogenic phenotypes of D. arachidicola.

Under Pressure: A Comparative Study of Botrytis cinerea Populations from Conventional and Organic Farms in Cyprus and Greece

The highly heterogeneous nature of Botrytis cinerea provides adaptive benefits to variable environmental regimes. Disentangling pathogen population structure in anthropogenic agroecosystems is crucial to designing more effective management schemes. Herein, we studied how evolutionary forces exerted in different farming systems, in terms of agrochemicals-input, shape B. cinerea populations. In total, 360 B. cinerea isolates were collected from conventional and organic, strawberry and tomato farms in Cyprus and Greece. The occurrence and frequency of sensitivities to seven botryticides were estimated. Results highlighted widespread fungicide resistance in conventional farms since only 15.5% of the isolates were sensitive. A considerable frequency of fungicide-resistant isolates was also detected in the organic farms (14.9%). High resistance frequencies were observed for boscalid (67.7%), pyraclostrobin (67.3%), cyprodinil (65.9%), and thiophanate-methyl (61.4%) in conventional farms, while high levels of multiple fungicide resistance were also evident. Furthermore, B. cinerea isolates were genotyped using a set of seven microsatellite markers (simple sequence repeat [SSR] markers). Index of association analyses (Ia and rBarD) suggest asexual reproduction of the populations, even though the mating-type idiomorphs were equally distributed, indicating frequency-dependent selection. Fungicide resistance was correlated with farming systems across countries and crops, while SSRs were able to detect population structure associated with resistance to thiophanate-methyl, pyraclostrobin, boscalid, and cyprodinil. The expected heterozygosity in organic farms was significantly higher than in conventional, suggesting the absence of selective pressure that may change the allelic abundance in organic farms. However, genetic variance among strawberry and tomato populations was high, ranking host specificity higher than other selection forces studied.

Genomic landscape of a relict fir-associated fungus reveals rapid convergent adaptation towards endophytism

Comparative and pan-genomic analyses of the endophytic fungus Pezicula neosporulosa (Helotiales, Ascomycota) from needles of the relict fir, Abies beshanzuensis, showed expansions of carbohydrate metabolism and secondary metabolite biosynthetic genes characteristic for unrelated plant-beneficial helotialean, such as dark septate endophytes and ericoid mycorrhizal fungi. The current species within the relatively young Pliocene genus Pezicula are predominantly saprotrophic, while P. neosporulosa lacks such features. To understand the genomic background of this putatively convergent evolution, we performed population analyses of 77 P. neosporulosa isolates. This revealed a mosaic structure of a dozen non-recombining and highly genetically polymorphic subpopulations with a unique mating system structure. We found that one idiomorph of a probably duplicated mat1-2 gene was found in putatively heterothallic isolates, while the other co-occurred with mat1-1 locus suggesting homothallic reproduction for these strains. Moreover, 24 and 81 genes implicated in plant cell-wall degradation and secondary metabolite biosynthesis, respectively, showed signatures of the balancing selection. These findings highlight the evolutionary pattern of the two gene families for allowing the fungus a rapid adaptation towards endophytism and facilitating diverse symbiotic interactions.

Botrytis cinerea Loss and Restoration of Virulence during In Vitro Culture Follows Flux in Global DNA Methylation

Pathogenic fungi can lose virulence after protracted periods of culture, but little is known of the underlying mechanisms. Here, we present the first analysis of DNA methylation flux at a single-base resolution for the plant pathogen B. cinerea and identify differentially methylated genes/genomic regions associated with virulence erosion during in vitro culture. Cultures were maintained for eight months, with subcultures and virulence testing every month. Methylation-sensitive amplified polymorphisms were performed at monthly intervals to characterise global changes to the pathogen’s genome during culture and also on DNA from mycelium inoculated onto Arabidopsis thaliana after eight months in culture. Characterisation of culture-induced epialleles was assessed by whole-genome re-sequencing and whole-genome bisulfite sequencing. Virulence declined with time in culture and recovered after inoculation on A. thaliana. Variation detected by methylation-sensitive amplified polymorphisms followed virulence changes during culture. Whole-genome (bisulfite) sequencing showed marked changes in global and local methylation during culture but no significant genetic changes. We imply that virulence is a non-essential plastic character that is at least partly modified by the changing levels of DNA methylation during culture. We hypothesise that changing DNA methylation during culture may be responsible for the high virulence/low virulence transition in B. cinerea and speculate that this may offer fresh opportunities to control pathogen virulence.

Population Genomics Reveals Molecular Determinants of Specialization to Tomato in the Polyphagous Fungal Pathogen Botrytis cinerea in France

Many fungal plant pathogens encompass multiple populations specialized on different plant species. Understanding the factors underlying pathogen adaptation to their hosts is a major challenge of evolutionary microbiology, and it should help to prevent the emergence of new specialized pathogens on novel hosts. Previous studies have shown that French populations of the gray mold pathogen Botrytis cinerea parasitizing tomato and grapevine are differentiated from each other, and have higher aggressiveness on their host of origin than on other hosts, indicating some degree of host specialization in this polyphagous pathogen. Here, we aimed at identifying the genomic features underlying the specialization of B. cinerea populations to tomato and grapevine. Based on whole genome sequences of 32 isolates, we confirmed the subdivision of B. cinerea pathogens into two genetic clusters on grapevine and another, single cluster on tomato. Levels of genetic variation in the different clusters were similar, suggesting that the tomato-specific cluster has not recently emerged following a bottleneck. Using genome scans for selective sweeps and divergent selection, tests of positive selection based on polymorphism and divergence at synonymous and nonsynonymous sites, and analyses of presence and absence variation, we identified several candidate genes that represent possible determinants of host specialization in the tomato-associated population. This work deepens our understanding of the genomic changes underlying the specialization of fungal pathogen populations.

Deciphering the Dynamics of Signaling Cascades and Virulence Factors of B. cinerea during Tomato Cell Wall Degradation

The ascomycete Botrytis cinerea is one of the most relevant plant pathogenic fungi, affecting fruits, flowers, and greenhouse-grown crops. The infection strategy used by the fungus comprises a magnificent set of tools to penetrate and overcome plant defenses. In this context, the plant-pathogen communication through membrane receptors and signal transduction cascades is essential to trigger specific routes and the final success of the infection. In previous reports, proteomics approaches to B. cinerea signal transduction cascades changes in response to different carbon source and plant-based elicitors have been performed. Analyzing the secretome, membranome, phosphoproteome, and the phosphomembranome. Moreover, phenotypic changes in fungal biology was analyzed, specifically toxin production. To obtain the whole picture of the process and reveal the network from a system biology approach, this proteomic information has been merged with the phenotypic characterization, to be analyzed using several bioinformatics algorithms (GO, STRING, MCODE) in order to unravel key points in the signal transduction regulation crucial to overcome plant defenses, as well as new virulence/pathogenicity factors that could be used as therapeutic targets in the control of the gray mold rot disease. A total of 1721 and 663 exclusive or overexpressed proteins were identified under glucose (GLU) and deproteinized tomato cell walls (TCW), summarizing all of the protein identifications under phenotypic characterized stages. Under GO analysis, there are more biological process and molecular functions described in GLU, highlighting the increase in signaling related categories. These results agree with the high number of total identified proteins in GLU, probably indicating a more varied and active metabolism of the fungus. When analyzing only GO annotations related with signal transduction, it was revealed that there were proteins related to TOR signaling, the phosphorelay signal transduction system, and inositol lipid-mediated signaling, only under GLU conditions. On the contrary, calcium-mediated signaling GO annotation is only present between the proteins identified under TCW conditions. To establish a potential relationship between expressed proteins, cluster analyses showed 41 and 14 clusters under GLU and TCW conditions, confirming an increase in biological activity in GLU, where we identified a larger number of clusters related to transcription, translation, and cell division, between others. From these analyses, clusters related to signal transduction and clusters related to mycotoxin production were found, which correlated with the phenotypic characterization. The identification of the proteins encompassed in each condition and signal transduction cascade would provide the research community with new information about the B. cinerea infection process and potential candidates of pathogenicity/virulence factors, overcoming plant defenses, and new therapeutic targets.

Quantitative interactions: the disease outcome of Botrytis cinerea across the plant kingdom

Botrytis cinerea is a fungal pathogen that causes necrotic disease on more than a thousand known hosts widely spread across the plant kingdom. How B. cinerea interacts with such extensive host diversity remains largely unknown. To address this question, we generated an infectivity matrix of 98 strains of B. cinerea on 90 genotypes representing eight host plants. This experimental infectivity matrix revealed that the disease outcome is largely explained by variations in either the host resistance or pathogen virulence. However, the specific interactions between host and pathogen account for 16% of the disease outcome. Furthermore, the disease outcomes cluster among genotypes of a species but are independent of the relatedness between hosts. When analyzing the host specificity and virulence of B. cinerea, generalist strains are predominant. In this fungal necrotroph, specialization may happen by a loss in virulence on most hosts rather than an increase of virulence on a specific host. To uncover the genetic architecture of Botrytis host specificity and virulence, a genome-wide association study (GWAS) was performed and revealed up to 1492 genes of interest. The genetic architecture of these traits is widespread across the B. cinerea genome. The complexity of the disease outcome might be explained by hundreds of functionally diverse genes putatively involved in adjusting the infection to diverse hosts.

Eco-friendly management of postharvest fungal decays in kiwifruit

Kiwifruit is purchased by consumers worldwide and is increasing in demand. Unfortunately, kiwifruit is susceptible to postharvest decay caused by a variety of fungal pathogens, including Botrytis cinerea, Penicillium expansum, Alternaria alternata, Botryosphaeria dothidea, and Diaporthe spp. Among these pathogens, B. cinerea is the most prevalent and devastating. Infections by these fungal pathogens result in a deterioration in fruit quality and a reduction in marketable yield. Eco-friendly methods to control kiwifruit postharvest decay have been explored as alternatives to the use of synthetic fungicides. In this review, we provide an overview and discuss the virulence and pathogenesis of fungi that are causal agents of kiwifruit decay, especially B. cinerea, including recent molecular and genomic studies. Advances in pre- and postharvest measures for postharvest decay management, including biological control, physical applications, the use of natural compounds and plant hormones, and the use of combined methods, are also reviewed. Eco-friendly control measures are a critical component of an integrated management approach for sustainable production of kiwifruit. The need for further research on the use of microbial consortia for the management of postharvest diseases of kiwifruit is also discussed.

A Major Effect Gene Controlling Development and Pathogenicity in Botrytis cinerea Identified Through Genetic Analysis of Natural Mycelial Non-pathogenic Isolates

Botrytis cinerea is a necrotrophic plant pathogenic fungus with a wide host range. Its natural populations are phenotypically and genetically very diverse. A survey of B. cinerea isolates causing gray mold in the vineyards of Castilla y León, Spain, was carried out and as a result eight non-pathogenic natural variants were identified. Phenotypically these isolates belong to two groups. The first group consists of seven isolates displaying a characteristic mycelial morphotype, which do not sporulate and is unable to produce sclerotia. The second group includes one isolate, which sporulates profusely and does not produce sclerotia. All of them are unresponsive to light. Crosses between a representative mycelial non-pathogenic isolate and a highly aggressive field isolate revealed that the phenotypic differences regarding pathogenicity, sporulation and production of sclerotia cosegregated in the progeny and are determined by a single genetic locus. By applying a bulked segregant analysis strategy based on the comparison of the two parental genomes the locus was mapped to a 110 kb region in chromosome 4. Subcloning and transformation experiments revealed that the polymorphism is an SNP affecting gene Bcin04g03490 in the reference genome of B. cinerea. Genetic complementation analysis and sequencing of the Bcin04g03490 alleles demonstrated that the mutations in the mycelial isolates are allelic and informed about the nature of the alterations causing the phenotypes observed. Integration of the allele of the pathogenic isolate into the non-pathogenic isolate fully restored the ability to infect, to sporulate and to produce sclerotia. Therefore, it is concluded that a major effect gene controlling differentiation and developmental processes as well as pathogenicity has been identified in B. cinerea. It encodes a protein with a GAL4-like Zn(II)2Cys6 binuclear cluster DNA binding domain and an acetyltransferase domain, suggesting a role in regulation of gene expression through a mechanism involving acetylation of specific substrates.

Two sides of the same story in grapevine-pathogen interactions

Proteases are an integral part of plant defence systems, and their role in plant-pathogen interactions is unequivocal. Emerging evidence suggests that different protease families contribute to the establishment not only of hypersensitive response, priming, and signalling, but also of recognition events through complex proteolytic cascades. Moreover, they play a crucial role in pathogen/microbe-associated molecular pattern (PAMP/MAMP)-triggered immunity as well as in effector-triggered immunity. However, despite important advances in our understanding of the role of proteases in plant defence, the contribution of proteases to pathogen defence in grapevine remains poorly understood. In this review, we summarize current knowledge of the main grapevine pathosystems and explore the role of serine, cysteine, and aspartic proteases from both the host and pathogen point of views.

The Destructive Fungal Pathogen Botrytis cinerea-Insights from Genes Studied with Mutant Analysis

Botrytis cinerea is one of the most destructive fungal pathogens affecting numerous plant hosts, including many important crop species. As a molecularly under-studied organism, its genome was only sequenced at the beginning of this century and it was recently updated with improved gene annotation and completeness. In this review, we summarize key molecular studies on B. cinerea developmental and pathogenesis processes, specifically on genes studied comprehensively with mutant analysis. Analyses of these studies have unveiled key genes in the biological processes of this pathogen, including hyphal growth, sclerotial formation, conidiation, pathogenicity and melanization. In addition, our synthesis has uncovered gaps in the present knowledge regarding development and virulence mechanisms. We hope this review will serve to enhance the knowledge of the biological mechanisms behind this notorious fungal pathogen.

Influence of genetic diversity of seventeen Beauveria bassiana isolates from different hosts on virulence by comparative genomics

Background

Beauveria bassiana (B. bassiana) is a famous entomopathogenic fungus that could parasitize on hundreds of insect species, which are being used as an environmentally friendly mycoinsecticide. Nevertheless, the possible effect of genetic diversity of these B. bassiana isolates from different hosts on virulence has not been explored before. In order to explore that issue, we compared the genome sequences among seventeen B. bassiana isolates from 17 different insects using whole genome re-sequencing, with B. bassiana strain ARSEF 2860 as the reference genome.

Results

There were a total of 10,098 missense mutated genes, 720 positively selected genes were identified in 17 strains of B. bassiana. Among these, two genes with high frequency mutations encode the toxin-producing non-ribosomal peptide synthase (NRPS) protein. Seven genes undergoing positive selection were enriched in the two-component signaling pathway that is known to regulate the fungal toxicity. In addition, the domain changes of three positively selected genes are also directly related to the virulence plasticity. Besides, the functional categorization of mutated genes showed that most of them involved in the biological functions of toxic proteins involved in.

Conclusions

Based on our data, our results indicate that several mutated genes and positively selected genes may underpin virulence of B. bassiana towards hosts during infection process, which provide an insight into the potential effects of natural variation on the virulence of B. bassiana, which will be useful in screening out potential virulence factors in B. bassiana.

Population Genetics and Fungicide Resistance of Botrytis cinerea on Vitis and Prunus spp. in California

Botrytis cinerea, the causal agent of gray mold, has high genetic diversity and a broad host range. In Vitis sp. and Prunus spp., B. cinerea causes pre- and postharvest diseases, and fungicides are routinely applied to prevent yield loss. In total, 535 isolates of B. cinerea collected from Vitis sp. and Prunus spp. in 2012, 2016, and 2017 were genotyped using 18 microsatellite markers and the transposable elements (TEs) Boty and Flipper. Only nine of the polymorphic markers and the two TEs were considered informative and retained for the final analyses. Of the 532 isolates, 297 were tested for resistance to seven fungicides representing six Fungicide Resistance Action Committee classes. After clone correction, 295 multilocus genotype groups were retained across the 3 years in 326 individuals, and four genetic subpopulations were detected. High levels of clonality were observed across the dataset. Significant pairwise differentiation was detected among years, locations, and TE composition. However, most of the diversity observed was within a subpopulation and not among subpopulations. No genetic differentiation was detected among resistant and sensitive isolates for individual fungicide classes. When resistance to the total number of fungicides was compared, regardless of the fungicide class, significant differentiation was detected among isolates that are resistant to two fungicide classes and those resistant to three or four fungicide groups. Fungicide resistance frequencies were stable for most chemistries evaluated with the exception of fluopyram, which increased from 2012 to 2016/2017.

Prediction of pathogenicity genes involved in adaptation to a lupin host in the fungal pathogens Botrytis cinerea and Sclerotinia sclerotiorum via comparative genomics

Background

Narrow-leafed lupin is an emerging crop of significance in agriculture, livestock feed and human health food. However, its susceptibility to various diseases is a major obstacle towards increased adoption. Sclerotinia sclerotiorum and Botrytis cinerea – both necrotrophs with broad host-ranges - are reported among the top 10 lupin pathogens. Whole-genome sequencing and comparative genomics are useful tools to discover genes responsible for interactions between pathogens and their hosts.

Results

Genomes were assembled for one isolate of B. cinerea and two isolates of S. sclerotiorum, which were isolated from either narrow-leafed or pearl lupin species. Comparative genomics analysis between lupin-derived isolates and others isolated from alternate hosts was used to predict between 94 to 98 effector gene candidates from among their respective non-conserved gene contents.

Conclusions

Detection of minor differences between relatively recently-diverged isolates, originating from distinct regions and with hosts, may highlight novel or recent gene mutations and losses resulting from host adaptation in broad host-range fungal pathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5774-2) contains supplementary material, which is available to authorized users.

Plant-necrotroph co-transcriptome networks illuminate a metabolic battlefield

A central goal of studying host-pathogen interaction is to understand how host and pathogen manipulate each other to promote their own fitness in a pathosystem. Co-transcriptomic approaches can simultaneously analyze dual transcriptomes during infection and provide a systematic map of the cross-kingdom communication between two species. Here we used the Arabidopsis-B. cinerea pathosystem to test how plant host and fungal pathogen interact at the transcriptomic level. We assessed the impact of genetic diversity in pathogen and host by utilization of a collection of 96 isolates infection on Arabidopsis wild-type and two mutants with jasmonate or salicylic acid compromised immunities. We identified ten B. cinereagene co-expression networks (GCNs) that encode known or novel virulence mechanisms. Construction of a dual interaction network by combining four host- and ten pathogen-GCNs revealed potential connections between the fungal and plant GCNs. These co-transcriptome data shed lights on the potential mechanisms underlying host-pathogen interaction.

Phenotypic and Genetic Characterization of Botrytis cinerea Population from Kiwifruit in Sichuan Province, China

Botrytis cinerea (anamorph of Botryotinia fuckeliana) causes gray mold on numerous plants, including kiwifruit. The primary aim of this study was to investigate the phenotypic and genetic characteristics of the Botrytis cinerea population from kiwifruit in Sichuan Province, China. In all, 176 isolates were collected from kiwifruit orchards from eight geographic regions in Sichuan. All isolates were identified as B. cinerea sensu stricto based on the combined datasets, including morphological criteria, determination of the Bc-hch allele, and phylogenetic analysis of the genes RPB2, G3PDH, and HSP60. Three colony types (i.e., sclerotial, mycelial, and conidial) were observed on potato dextrose agar after 2 weeks, with sclerotial isolates, the predominant category, accounting for 40.91%. No obvious differences in microscopic characteristics were observed among the three types. Three genotypes of transposable elements were identified in the B. cinerea population: boty, flipper, and transposa types. The most prevalent genotype from different geographic populations of B. cinerea was transposa; in contrast, the flipper genotype accounted for only 3.98% of the total population, whereas the vacuma genotype was absent. According to MAT locus amplification, 87 and 89 isolates are MAT1-1 and MAT1-2 type, respectively, and the two mating types were found to be balanced overall in the population. Forty-eight representative isolates were all able to cause gray mold to some extent, and disease severities were significantly different between the cultivars Hongyang and Hort16A (P < 0.01). Disease severity was significantly greater on young leaves than on mature leaves (P < 0.01). No significant relationship was found between pathogenicity and geographical region, colony type, or transposon distribution. The results obtained in the present study suggest a relatively uniform species diversity of Botrytis but rich phenotypic and genetic differentiation within the B. cinerea population on kiwifruit in China. Utilizing resistant cultivars and rain-shelter cultivation instead of fungicides may be an effective approach to delaying pathogen variability.

Interactions of Tomato and Botrytis cinerea Genetic Diversity: Parsing the Contributions of Host Differentiation, Domestication, and Pathogen Variation

Although the impacts of crop domestication on specialist pathogens are well known, less is known about the interaction of crop variation and generalist pathogens. To study how genetic variation within a crop affects plant resistance to generalist pathogens, we infected a collection of wild and domesticated tomato accessions with a genetically diverse population of the generalist pathogen Botrytis cinerea We quantified variation in lesion size of 97 B. cinerea genotypes (isolates) on six domesticated tomato genotypes (Solanum lycopersicum) and six wild tomato genotypes (Solanum pimpinellifolium). Lesion size was significantly affected by large effects of the host and pathogen's genotype, with a much smaller contribution of domestication. This pathogen collection also enables genome-wide association mapping of B. cinerea Genome-wide association mapping of the pathogen showed that virulence is highly polygenic and involves a diversity of mechanisms. Breeding against this pathogen would likely require the use of diverse isolates to capture all possible mechanisms. Critically, we identified a subset of B. cinerea genes where allelic variation was linked to altered virulence against wild versus domesticated tomato, as well as loci that could handle both groups. This generalist pathogen already has a large collection of allelic variation that must be considered when designing a breeding program.

Digital Imaging Combined with Genome-Wide Association Mapping Links Loci to Plant-Pathogen Interaction Traits

Plant resistance to generalist pathogens with broad host ranges, such as Botrytis cinerea (Botrytis), is typically quantitative and highly polygenic. Recent studies have begun to elucidate the molecular genetic basis of plant-pathogen interactions using commonly measured traits, including lesion size and/or pathogen biomass. However, with the advent of digital imaging and high-throughput phenomics, there are a large number of additional traits available to study quantitative resistance. In this study, we used high-throughput digital imaging analysis to investigate previously poorly characterized visual traits of plant-pathogen interactions related to disease resistance using the Arabidopsis (Arabidopsis thaliana)/Botrytis pathosystem. From a large collection of visual lesion trait measurements, we focused on color, shape, and size to test how these aspects of the Arabidopsis/Botrytis interaction are genetically related. Through genome-wide association mapping in Arabidopsis, we show that lesion color and shape are genetically separable traits associated with plant disease resistance. Moreover, by employing defined mutants in 23 candidate genes identified from the genome-wide association mapping, we demonstrate links between loci and each of the different plant-pathogen interaction traits. These results expand our understanding of the functional mechanisms driving plant disease resistance.

Many Shades of Grey in Botrytis-Host Plant Interactions

The grey mould Botrytis cinerea causes disease in more than 1000 plant species, including important crops. The interaction between Botrytis and its (potential) hosts is determined by quantitative susceptibility and virulence traits in both interacting partners, resulting in a greyscale of disease outcomes. Fungal infection was long thought to rely mainly on its capacity to kill the host plant and degrade plant tissue. Recent research has revealed that Botrytis exploits two crucial biological processes in host plants for its own success. We highlight recent findings that illustrate that the interactions between Botrytis and its host plants are subtle and we discuss the molecular and cellular mechanisms controlling the many shades of grey during these interactions.

Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem

To respond to pathogen attack, selection and associated evolution has led to the creation of plant immune system that are a highly effective and inducible defense system. Central to this system are the plant defense hormones jasmonic acid (JA) and salicylic acid (SA) and crosstalk between the two, which may play an important role in defense responses to specific pathogens or even genotypes. Here, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem to test how the host's defense system functions against genetic variation in a pathogen. We measured defense-related phenotypes and transcriptomic responses in Arabidopsis wild-type Col-0 and JA- and SA-signaling mutants, coi1-1 and npr1-1, individually challenged with 96 diverse B. cinerea isolates. Those data showed genetic variation in the pathogen influences on all components within the plant defense system at the transcriptional level. We identified four gene coexpression networks and two vectors of defense variation triggered by genetic variation in B. cinerea This showed that the JA and SA signaling pathways functioned to constrain/canalize the range of virulence in the pathogen population, but the underlying transcriptomic response was highly plastic. These data showed that plants utilize major defense hormone pathways to buffer disease resistance, but not the metabolic or transcriptional responses to genetic variation within a pathogen.

Quantitative Resistance: More Than Just Perception of a Pathogen

Molecular plant pathology has focused on studying large-effect qualitative resistance loci that predominantly function in detecting pathogens and/or transmitting signals resulting from pathogen detection. By contrast, less is known about quantitative resistance loci, particularly the molecular mechanisms controlling variation in quantitative resistance. Recent studies have provided insight into these mechanisms, showing that genetic variation at hundreds of causal genes may underpin quantitative resistance. Loci controlling quantitative resistance contain some of the same causal genes that mediate qualitative resistance, but the predominant mechanisms of quantitative resistance extend beyond pathogen recognition. Indeed, most causal genes for quantitative resistance encode specific defense-related outputs such as strengthening of the cell wall or defense compound biosynthesis. Extending previous work on qualitative resistance to focus on the mechanisms of quantitative resistance, such as the link between perception of microbe-associated molecular patterns and growth, has shown that the mechanisms underlying these defense outputs are also highly polygenic. Studies that include genetic variation in the pathogen have begun to highlight a potential need to rethink how the field considers broad-spectrum resistance and how it is affected by genetic variation within pathogen species and between pathogen species. These studies are broadening our understanding of quantitative resistance and highlighting the potentially vast scale of the genetic basis of quantitative resistance.

Landscape of genomic diversity and host adaptation in Fusarium graminearum

BACKGROUND

Fusarium graminearum is one of the main causal agents of the Fusarium Head Blight, a worldwide disease affecting cereal cultures, whose presence can lead to contaminated grains with chemically stable and harmful mycotoxins. Resistant cultivars and fungicides are frequently used to control this pathogen, and several observations suggest an adaptation of F. graminearum that raises concerns regarding the future of current plant disease management strategies. To understand the genetic basis as well as the extent of its adaptive potential, we investigated the landscape of genomic diversity among six French isolates of F. graminearum, at single-nucleotide resolution using whole-genome re-sequencing.

RESULTS

A total of 242,756 high-confidence genetic variants were detected when compared to the reference genome, among which 96% are single nucleotides polymorphisms. One third of these variants were observed in all isolates. Seventy-seven percent of the total polymorphism is located in 32% of the total length of the genome, comprising telomeric/subtelomeric regions as well as discrete interstitial sections, delineating clear variant enriched genomic regions- 7.5 times in average. About 80% of all the F. graminearum protein-coding genes were found polymorphic. Biological functions are not equally affected: genes potentially involved in host adaptation are preferentially located within polymorphic islands and show greater diversification rate than genes fulfilling basal functions. We further identified 29 putative effector genes enriched with non-synonymous effect mutation.

CONCLUSIONS

Our results highlight a remarkable level of polymorphism in the genome of F. graminearum distributed in a specific pattern. Indeed, the landscape of genomic diversity follows a bi-partite organization of the genome according to polymorphism and biological functions. We measured, for the first time, the level of sequence diversity for the entire gene repertoire of F. graminearum and revealed that the majority are polymorphic. Those assumed to play a role in host-pathogen interaction are discussed, in the light of the subsequent consequences for host adaptation. The annotated genetic variants discovered for this major pathogen are valuable resources for further genetic and genomic studies.

Expansive Phenotypic Landscape of Botrytis cinerea Shows Differential Contribution of Genetic Diversity and Plasticity

The modern evolutionary synthesis suggests that both environmental variation and genetic diversity are critical determinants of pathogen success. However, the relative contribution of these two sources of variation is not routinely measured. To estimate the relative contribution of plasticity and genetic diversity for virulence-associated phenotypes in a generalist plant pathogen, we grew a population of 15 isolates of Botrytis cinerea from throughout the world, under a variety of in vitro and in planta conditions. Under in planta conditions, phenotypic differences between the isolates were determined by the combination of genotypic variation within the pathogen and environmental variation. In contrast, phenotypic differences between the isolates under in vitro conditions were predominantly determined by genetic variation in the pathogen. Using a correlation network approach, we link the phenotypic variation under in vitro experimental conditions to phenotypic variation during plant infection. This study indicates that there is a high level of phenotypic variation within B. cinerea that is controlled by a mixture of genetic variation, environment, and genotype × environment. This argues that future experiments into the pathogenicity of B. cinerea must account for the genetic and environmental variation within the pathogen to better sample the potential phenotypic space of the pathogen.

The Quantitative Basis of the Arabidopsis Innate Immune System to Endemic Pathogens Depends on Pathogen Genetics

The most established model of the eukaryotic innate immune system is derived from examples of large effect monogenic quantitative resistance to pathogens. However, many host-pathogen interactions involve many genes of small to medium effect and exhibit quantitative resistance. We used the Arabidopsis-Botrytis pathosystem to explore the quantitative genetic architecture underlying host innate immune system in a population of Arabidopsis thaliana. By infecting a diverse panel of Arabidopsis accessions with four phenotypically and genotypically distinct isolates of the fungal necrotroph B. cinerea, we identified a total of 2,982 genes associated with quantitative resistance using lesion area and 3,354 genes associated with camalexin production as measures of the interaction. Most genes were associated with resistance to a specific Botrytis isolate, which demonstrates the influence of pathogen genetic variation in analyzing host quantitative resistance. While known resistance genes, such as receptor-like kinases (RLKs) and nucleotide-binding site leucine-rich repeat proteins (NLRs), were found to be enriched among associated genes, they only account for a small fraction of the total genes associated with quantitative resistance. Using publically available co-expression data, we condensed the quantitative resistance associated genes into co-expressed gene networks. GO analysis of these networks implicated several biological processes commonly connected to disease resistance, including defense hormone signaling and ROS production, as well as novel processes, such as leaf development. Validation of single gene T-DNA knockouts in a Col-0 background demonstrate a high success rate (60%) when accounting for differences in environmental and Botrytis genetic variation. This study shows that the genetic architecture underlying host innate immune system is extremely complex and is likely able to sense and respond to differential virulence among pathogen genotypes.

'Omics' and Plant Responses to Botrytis cinerea

Botrytis cinerea is a dangerous plant pathogenic fungus with wide host ranges. This aggressive pathogen uses multiple weapons to invade and cause serious damages on its host plants. The continuing efforts of how to solve the "puzzle" of the multigenic nature of B. cinerea's pathogenesis and plant defense mechanisms against the disease caused by this mold, the integration of omic approaches, including genomics, transcriptomics, proteomics and metabolomics, along with functional analysis could be a potential solution. Omic studies will provide a foundation for development of genetic manipulation and breeding programs that will eventually lead to crop improvement and protection. In this mini-review, we will highlight the current progresses in research in plant stress responses to B. cinerea using high-throughput omic technologies. We also discuss the opportunities that omic technologies can provide to research on B. cinerea-plant interactions as an example showing the impacts of omics on agricultural research.

Isolate Dependency of Brassica rapa Resistance QTLs to Botrytis cinerea

Generalist necrotrophic pathogens including Botrytis cinerea cause significant yield and financial losses on Brassica crops. However, there is little knowledge about the mechanisms underlying the complex interactions encoded by both host and pathogen genomes in this interaction. This potentially includes multiple layers of plant defense and pathogen virulence mechanisms that could complicate in breeding broad spectrum resistance within Brassica species. Glucosinolates (GSLs) are a diverse group of defense metabolites that play a key role in interaction between Brassica and biotic attackers. In this study, we utilized a collection of diverse B. cinerea isolates to investigate resistance within the Brassica rapa R500 × IMB211 recombinant inbred line population. We tested variation on lesion development and glucosinolate accumulation in parental lines and all population lines. We then mapped quantitative trait loci (QTL) for both resistances to B. cinerea and defense metabolites in this population. Phenotypic analysis and QTL mapping demonstrate that the genetic basis of resistance to B. cinerea in B. rapa is isolate specific and polygenic with transgressive segregation that both parents contribute resistance alleles. QTLs controlling defensive GSLs are highly dependent on pathogen infection. An overlap of two QTLs identified between resistance to B. cinerea and defense metabolites also showed isolate specific effects. This work suggests that directly searching for resistance loci may not be the best approach at improving resistance in B. rapa to necrotrophic pathogen.