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

Challenges and Opportunities Arising from Host-Botrytis cinerea Interactions to Outline Novel and Sustainable Control Strategies: The Key Role of RNA Interference

The necrotrophic plant pathogenic fungus Botrytis cinerea (Pers., 1794), the causative agent of gray mold disease, causes significant losses in agricultural production. Control of this fungal pathogen is quite difficult due to its wide host range and environmental persistence. Currently, the management of the disease is still mainly based on chemicals, which can have harmful effects not only on the environment and on human health but also because they favor the development of strains resistant to fungicides. The flexibility and plasticity of B. cinerea in challenging plant defense mechanisms and its ability to evolve strategies to escape chemicals require the development of new control strategies for successful disease management. In this review, some aspects of the host-pathogen interactions from which novel and sustainable control strategies could be developed (e.g., signaling pathways, molecules involved in plant immune mechanisms, hormones, post-transcriptional gene silencing) were analyzed. New biotechnological tools based on the use of RNA interference (RNAi) are emerging in the crop protection scenario as versatile, sustainable, effective, and environmentally friendly alternatives to the use of chemicals. RNAi-based fungicides are expected to be approved soon, although they will face several challenges before reaching the market.

Evaluation of biocontrol efficacy of rhizosphere dwelling bacteria for management of Fusarium wilt and Botrytis gray mold of chickpea

BACKGROUND

Chickpea (Cicer arietinum L.) production is affected by many biotic factors, among them Fusarium wilt caused by Fusarium oxysporum f. sp. ciceri and Botrytis gray mold caused by Botrytis cinerea led to severe losses. As fungicide application is not advisable, biological management is the best alternative for plant protection. The rhizosphere-dwelling antagonistic bacteria are one of the important successful alternative strategy to manage these diseases of chickpea. Rhizosphere dwelling bacteria serve as biocontrol agents by different mechanisms like producing antibiotics, different enzymes, siderophores against pathogens and thereby reducing the growth of pathogens.

RESULTS

The present study aimed to isolate rhizospheric bacteria from the soils of different chickpea fields to evaluate biocontrol efficacy of the isolated bacteria to manage Fusarium wilt and Botrytis gray mold in chickpea. A total of 67 bacteria were isolated from chickpea rhizosphere from Bundelkhand region of India. Study revealed the isolated bacteria could reduce the Fusarium oxysporum f. sp. ciceris and Botrytis cinerea infection in chickpea between 17.29 and 75.29%. After screening of all the bacteria for their biocontrol efficacy, 13 most promising bacterial isolates were considered for further study out of which, three bacterial isolates (15d, 9c and 14a) have shown the maximum in vitro antagonistic effects against Fusarium oxysporum f. sp. ciceri and Botrytis cinerea comparable to in vivo effects. However, Isolate (15d) showed highest 87.5% and 82.69% reduction in disease against Fusarium wilt and Botrytis gray mold respectively, under pot condition. Three most potential isolates were characterized at molecular level using 16S rRNA gene and found to be Priestia megaterium (9c and 14a) and Serratia marcescens (15d).

CONCLUSION

This study identified two native biocontrol agents Priestia megaterium and Serratia marcescens from the rhizospheric soils of Bundelkhand region of India for control of Fusarium wilt, Botrytis gray mold. In future, efforts should be made to further validate the biocontrol agents in conjugation with nanomaterials for enhancing the synergistic effects in managing the fungal diseases in chickpea. This study will definitely enhance our understanding of these bioagents, and to increase their performance by developing effective formulations, application methods, and integrated strategies.

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.

Plant immune inducer ZNC promotes rutin accumulation and enhances resistance to Botrytis cinerea in tomato

Gray mold is a destructive disease caused by Botrytis cinerea, a pervasive plant pathogen, which poses a threat to both tomato growth and postharvest storage. The utilization of induced resistance presents a potential strategy for combating plant pathogenic attacks. ZNC (zhinengcong), an extract derived from the endophytic fungus Paecilomyces variotii, has been discovered to play a vital role in preventing diverse forms of bacterial infections. Nevertheless, the precise mechanism behind its ability to enhance tomato resistance to fungi remains unclear. In this study, we found that the exogenous spraying of ZNC could significantly improve the resistance of tomato plants to B. cinerea. The results of both the metabolomic analysis and high-performance liquid chromatography (HPLC) demonstrated that tomato plants responded to ZNC treatment by accumulating high levels of rutin. Additional transcriptome analysis uncovered that rutin enhances tomato resistance possible by initiating the generation of reactive oxygen species (ROS) and phosphorylation of mitogen-activated protein kinases (MPKs) related genes expression during the initial phase of invasion by B. cinerea. In addition, we also found that rutin might activate plant immunity by eliciting ethylene (ET) and jasmonic acid (JA)-mediated pathways. Therefore, plant immune inducer ZNC and rutin has bright application prospects and high utilization value to control gray mold.

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.

Laminarin-triggered defence responses are geographically dependent in natural populations of Solanum chilense

Natural plant populations are polymorphic and show intraspecific variation in resistance properties against pathogens. The activation of the underlying defence responses can depend on variation in perception of pathogen-associated molecular patterns or elicitors. To dissect such variation, we evaluated the responses induced by laminarin (a glucan, representing an elicitor from oomycetes) in the wild tomato species Solanum chilense and correlated this to observed infection frequencies of Phytophthora infestans. We measured reactive oxygen species burst and levels of diverse phytohormones upon elicitation in 83 plants originating from nine populations. We found high diversity in basal and elicitor-induced levels of each component. Further we generated linear models to explain the observed infection frequency of P. infestans. The effect of individual components differed dependent on the geographical origin of the plants. We found that the resistance in the southern coastal region, but not in the other regions, was directly correlated to ethylene responses and confirmed this positive correlation using ethylene inhibition assays. Our findings reveal high diversity in the strength of defence responses within a species and the involvement of different components with a quantitatively different contribution of individual components to resistance in geographically separated populations of a wild plant species.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

Morphological, ultrastructural and molecular variations in susceptible and resistant genotypes of chickpea infected with Botrytis grey mould

Biotic stress due to fungal infection is detrimental to the growth and development of chickpea. In our study, two chickpea genotypes viz Cicer pinnatifidum (resistant) and PBG5 (susceptible) were inoculated with (1 × 104 spore mL−1) of nectrotrophic fungus Botrytis cinerea at seedling stage. These seedlings were evaluated for morphological, ultrastructural, and molecular differences after 3, 5 and 7 days post inoculation (dpi). Visual symptoms were recorded in terms of water-soaked lesions, rotten pods and twigs with fungal colonies. Light and scanning electron microscopy (SEM) revealed the differences in number of stomata, hyphal network and extent of topographical damage in resistant (C. pinnatifidum) and susceptible (PBG5) genotypes, which were validated by stomatal index studies done by using fluorescence microscopy in the infection process of B. cinerea in leaves of both chickpea genotypes. In case of control (water inoculated) samples, there were differences in PCR analysis done using five primers for screening the genetic variations between two genotypes. The presence of a Botrytis responsive gene (LrWRKY) of size ~300 bp was observed in uninoculated resistant genotype which might have a role in resistance against Botrytis grey mould. The present investigation provides information about the variation in the infection process of B. cinerea in two genotypes which can be further exploited to develop robust and effective strategies to manage grey mould disease.

Quantitative Genetic Analysis of Interactions in the Pepper-Phytophthora capsici Pathosystem

The development of pepper cultivars with durable resistance to the oomycete Phytophthora capsici has been challenging due to differential interactions between the species that allow certain pathogen isolates to cause disease on otherwise resistant host genotypes. Currently, little is known about the pathogen genes involved in these interactions. To investigate the genetic basis of P. capsici virulence on individual pepper genotypes, we inoculated sixteen pepper accessions, representing commercial varieties, sources of resistance, and host differentials, with 117 isolates of P. capsici, for a total of 1,864 host-pathogen combinations. Analysis of disease outcomes revealed a significant effect of inter-species genotype-by-genotype interactions, although these interactions were quantitative rather than qualitative in scale. Isolates were classified into five pathogen subpopulations, as determined by their genotypes at over 60,000 single-nucleotide polymorphisms (SNPs). While absolute virulence levels on certain pepper accessions significantly differed between subpopulations, a multivariate phenotype reflecting relative virulence levels on certain pepper genotypes compared with others showed the strongest association with pathogen subpopulation. A genome-wide association study (GWAS) identified four pathogen loci significantly associated with virulence, two of which colocalized with putative RXLR effector genes and another with a polygalacturonase gene cluster. All four loci appeared to represent broad-spectrum virulence genes, as significant SNPs demonstrated consistent effects regardless of the host genotype tested. Host genotype-specific virulence variants in P. capsici may be difficult to map via GWAS with all but excessively large sample sizes, perhaps controlled by genes of small effect or by multiple allelic variants that have arisen independently. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Elucidating the effect of tomato leaf surface microstructure on Botrytis cinerea using synthetic systems

For some pathogenic fungi, sensing surface topography is part of their infection strategy. Their directional growth and transformation to a new developmental stage is influenced by contact with topographic features, which is referred to as thigmo-response, the exact functionality of which is not fully understood. Research on thigmo-responses is often performed on biomimetically patterned surfaces (BPS). Polydimethylsiloxane (PDMS) is especially suitable for fabrication of BPS. Here, we used synthetic BPS surfaces, mimicking tomato leaf surface, made from PDMS with the pathogenic fungus Botrytis cinerea to study the influence of structural features of the leaf surface on the fungus behavior. As a control, a PDMS surface without microstructure was fabricated to maintain the same chemical properties. Pre-penetration processes of B. cinerea, including the distribution of conidia on the surface, germination, and germ tube growth were observed on both leaf-patterned and flat PDMS. Microstructure affected the location of immediate attachment of conidia. Additionally, the microstructure of the plant host stimulated the development of germ tube in B. cinerea, at a higher rate than that observed on flat surface, suggesting that microstructure plays a role in fungus attachment and development.

Contrasting Roles of Ethylene Response Factors in Pathogen Response and Ripening in Fleshy Fruit

Fleshy fruits are generally hard and unpalatable when unripe; however, as they mature, their quality is transformed by the complex and dynamic genetic and biochemical process of ripening, which affects all cell compartments. Ripening fruits are enriched with nutrients such as acids, sugars, vitamins, attractive volatiles and pigments and develop a pleasant taste and texture and become attractive to eat. Ripening also increases sensitivity to pathogens, and this presents a crucial problem for fruit postharvest transport and storage: how to enhance pathogen resistance while maintaining ripening quality. Fruit development and ripening involve many changes in gene expression regulated by transcription factors (TFs), some of which respond to hormones such as auxin, abscisic acid (ABA) and ethylene. Ethylene response factor (ERF) TFs regulate both fruit ripening and resistance to pathogen stresses. Different ERFs regulate fruit ripening and/or pathogen responses in both fleshy climacteric and non-climacteric fruits and function cooperatively or independently of other TFs. In this review, we summarize the current status of studies on ERFs that regulate fruit ripening and responses to infection by several fungal pathogens, including a systematic ERF transcriptome analysis of fungal grey mould infection of tomato caused by Botrytis cinerea. This deepening understanding of the function of ERFs in fruit ripening and pathogen responses may identify novel approaches for engineering transcriptional regulation to improve fruit quality and pathogen resistance.

Genetic composition and diversity of Arabica coffee in the crop's center of origin and its impact on four major fungal diseases

Conventional wisdom states that genetic variation reduces disease levels in plant populations. Nevertheless, crop species have been subject to a gradual loss of genetic variation through selection for specific traits during breeding, thereby increasing their vulnerability to biotic stresses such as pathogens. We explored how genetic variation in Arabica coffee sites in southwestern Ethiopia was related to the incidence of four major fungal diseases. Sixty sites were selected along a gradient of management intensity, ranging from nearly wild to intensively managed coffee stands. We used genotyping-by-sequencing of pooled leaf samples (pool-GBS) derived from 16 individual coffee shrubs in each of the sixty sites to assess the variation in genetic composition (multivariate: reference allele frequency) and genetic diversity (univariate: mean expected heterozygosity) between sites. We found that genetic composition had a clear spatial pattern and that genetic diversity was higher in less managed sites. The incidence of the four fungal diseases was related to the genetic composition of the coffee stands, but in a specific way for each disease. In contrast, genetic diversity was only related to the within-site variation of coffee berry disease, but not to the mean incidence of any of the four diseases across sites. Given that fungal diseases are major challenges of Arabica coffee in its native range, our findings that genetic composition of coffee sites impacted the major fungal diseases may serve as baseline information to study the molecular basis of disease resistance in coffee. Overall, our study illustrates the need to consider both host genetic composition and genetic diversity when investigating the genetic basis for variation in disease levels.

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.

Multiomics analyses reveal the roles of the ASR1 transcription factor in tomato fruits

The transcription factor ASR1 (ABA, STRESS, RIPENING 1) plays multiple roles in plant responses to abiotic stresses as well as being involved in the regulation of central metabolism in several plant species. However, despite the high expression of ASR1 in tomato fruits, large scale analyses to uncover its function in fruits are still lacking. In order to study its function in the context of fruit ripening, we performed a multiomics analysis of ASR1-antisense transgenic tomato fruits at the transcriptome and metabolome levels. Our results indicate that ASR1 is involved in several pathways implicated in the fruit ripening process, including cell wall, amino acid, and carotenoid metabolism, as well as abiotic stress pathways. Moreover, we found that ASR1-antisense fruits are more susceptible to the infection by the necrotrophic fungus Botrytis cinerea. Given that ASR1 could be regulated by fruit ripening regulators such as FRUITFULL1/FRUITFULL2 (FUL1/FUL2), NON-RIPENING (NOR), and COLORLESS NON-RIPENING (CNR), we positioned it in the regulatory cascade of red ripe tomato fruits. These data extend the known range of functions of ASR1 as an important auxiliary regulator of tomato fruit ripening.

Tracking cell wall changes in wine and table grapes undergoing Botrytis cinerea infection using glycan microarrays

BACKGROUND AND AIMS

The necrotrophic fungus Botrytis cinerea infects a broad range of fruit crops including domesticated grapevine Vitis vinifera cultivars. Damage caused by this pathogen is severely detrimental to the table and wine grape industries and results in substantial crop losses worldwide. The apoplast and cell wall interface is an important setting where many plant-pathogen interactions take place and where some defence-related messenger molecules are generated. Limited studies have investigated changes in grape cell wall composition upon infection with B. cinerea, with much being inferred from studies on other fruit crops.

METHODS

In this study, comprehensive microarray polymer profiling in combination with monosaccharide compositional analysis was applied for the first time to investigate cell wall compositional changes in the berries of wine (Sauvignon Blanc and Cabernet Sauvignon) and table (Dauphine and Barlinka) grape cultivars during Botrytis infection and tissue maceration. This was used in conjunction with scanning electron microscopy (SEM) and X-ray computed tomography (CT) to characterize infection progression.

KEY RESULTS

Grapes infected at veraison did not develop visible infection symptoms, whereas grapes inoculated at the post-veraison and ripe stages showed evidence of significant tissue degradation. The latter was characterized by a reduction in signals for pectin epitopes in the berry cell walls, implying the degradation of pectin polymers. The table grape cultivars showed more severe infection symptoms, and corresponding pectin depolymerization, compared with wine grape cultivars. In both grape types, hemicellulose layers were largely unaffected, as was the arabinogalactan protein content, whereas in moderate to severely infected table grape cultivars, evidence of extensin epitope deposition was present.

CONCLUSIONS

Specific changes in the grape cell wall compositional profiles appear to correlate with fungal disease susceptibility. Cell wall factors important in influencing resistance may include pectin methylesterification profiles, as well as extensin reorganization.

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.

Polymorphisms in plants to restrict losses to pathogens: From gene family expansions to complex network evolution

Genetic polymorphisms are the basis of the natural diversity seen in all life on earth, also in plant-pathogen interactions. Initially, studies on plant-pathogen interaction focused on reporting phenotypic variation in resistance properties and on the identification of underlying major genes. Nowadays, the field of plant-pathogen interactions is moving from focusing on families of single dominant genes involved in gene-for-gene interactions to an understanding of the plant immune system in the context of a much more complex signaling network and quantitative resistance. Simultaneously, studies on pathosystems from the wild and genome analyses advanced, revealing tremendous variation in natural plant populations. It is now imperative to place studies on genetic diversity and evolution of plant-pathogen interactions in the appropriate molecular biological, as well as evolutionary, context.

Quantitative resistance differences between and within natural populations of Solanum chilense against the oomycete pathogen Phytophthora infestans

The wild tomato species Solanum chilense is divided into geographically and genetically distinct populations that show signs of defense gene selection and differential phenotypes when challenged with several phytopathogens, including the oomycete causal agent of late blight Phytophthora infestans. To better understand the phenotypic diversity of this disease resistance in S. chilense and to assess the effect of plant genotype versus pathogen isolate, respectively, we evaluated infection frequency in a systematic approach and with large sample sizes. We studied 85 genetically distinct individuals representing nine geographically separated populations of S. chilense. This showed that differences in quantitative resistance can be observed between but also within populations at the level of individual plants. Our data also did not reveal complete immunity in any of the genotypes. We further evaluated the resistance of a subset of the plants against P. infestans isolates with diverse virulence properties. This confirmed that the relative differences in resistance phenotypes between individuals were mainly determined by the plant genotype under consideration with modest effects of pathogen isolate used in the study. Thus, our report suggests that the observed quantitative resistance against P. infestans in natural populations of a wild tomato species S. chilense is the result of basal defense responses that depend on the host genotype and are pathogen isolate-unspecific.

Population Genetic Analyses of Botrytis cinerea Isolates From Michigan Vineyards Using a High-Throughput Marker System Approach

As sequencing costs continue to decrease, new tools are being developed for assessing pathogen diversity and population structure. Traditional marker types, such as microsatellites, are often more cost effective than single-nucleotide polymorphism (SNP) panels when working with small numbers of individuals, but may not allow for fine scale evaluation of low or moderate structure in populations. Botrytis cinerea is a necrotrophic plant pathogen with high genetic variability that can infect more than 200 plant species worldwide. A panel of 52 amplicons were sequenced for 82 isolates collected from four Michigan vineyards representing 2 years of collection and varying fungicide resistance. A panel of nine microsatellite markers previously described was also tested across 74 isolates from the same population. A microsatellite and SNP marker analysis of B. cinerea populations was performed to assess the genetic diversity and population structure of Michigan vineyards, and the results from both marker types were compared. Both methods were able to detect population structure associated with resistance to the individual fungicides thiabendazole and boscalid, and multiple fungicide resistance (MFR). Microsatellites were also able to differentiate population structure associated with another fungicide, fluopyram, while SNPs were able to additionally differentiate structure based on year. For both methods, AMOVA results were similar, with microsatellite results explaining a smaller portion of the variation compared with the SNP results. The SNP-based markers presented here were able to successfully differentiate population structure similar to microsatellite results. These SNP markers represent new tools to discriminate B. cinerea isolates within closely related populations using multiple targeted sequences.

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.

Tomato Domestication Attenuated Responsiveness to a Beneficial Soil Microbe for Plant Growth Promotion and Induction of Systemic Resistance to Foliar Pathogens

Crop domestication events followed by targeted breeding practices have been pivotal for improvement of desirable traits and to adapt cultivars to local environments. Domestication also resulted in a strong reduction in genetic diversity among modern cultivars compared to their wild relatives, though the effect this could have on tripartite relationships between plants, belowground beneficial microbes and aboveground pathogens remains undetermined. We quantified plant growth performance, basal resistance and induced systemic resistance (ISR) by Trichoderma harzianum, a beneficial soil microbe against Botrytis cinerea, a necrotrophic fungus and Phytophthora infestans, a hemi-biotrophic oomycete, in 25 diverse tomato genotypes. Wild tomato related species, tomato landraces and modern commercial cultivars that were conventionally or organically bred, together, representing a domestication gradient were evaluated. Relationships between basal and ISR, plant physiological status and phenolic compounds were quantified to identify potential mechanisms. Trichoderma enhanced shoot and root biomass and ISR to both pathogens in a genotype specific manner. Moreover, improvements in plant performance in response to Trichoderma gradually decreased along the domestication gradient. Wild relatives and landraces were more responsive to Trichoderma, resulting in greater suppression of foliar pathogens than modern cultivars. Photosynthetic rate and stomatal conductance of some tomato genotypes were improved by Trichoderma treatment whereas leaf nitrogen status of the majority of tomato genotypes were not altered. There was a negative relationship between basal resistance and induced resistance for both diseases, and a positive correlation between Trichoderma-ISR to B. cinerea and enhanced total flavonoid contents. These findings suggest that domestication and breeding practices have altered plant responsiveness to beneficial soil microbes. Further studies are needed to decipher the molecular mechanisms underlying the differential promotion of plant growth and resistance among genotypes, and identify molecular markers to integrate selection for responsiveness into future breeding programs.

Proteome-Wide Analysis of Lysine 2-Hydroxyisobutyrylation in the Phytopathogenic Fungus Botrytis cinerea

Posttranslational modifications (PTMs) of the whole proteome have become a hot topic in the research field of epigenetics, and an increasing number of PTM types have been identified and shown to play significant roles in different cellular processes. Protein lysine 2-hydroxyisobutyrylation (K_hib) is a newly detected PTM, and the 2-hydroxyisobutyrylome has been identified in several species. Botrytis cinerea is recognized as one of the most destructive pathogens due to its broad host distribution and very large economic losses; thus the many aspects of its pathogenesis have been continuously studied. However, distribution and function of K_hib in this phytopathogenic fungus are not clear. In this study, a proteome-wide analysis of K_hib in B. cinerea was performed, and 5,398 K_hib sites on 1,181 proteins were identified. Bioinformatics analysis showed that the 2-hydroxyisobutyrylome in B. cinerea contains both conserved proteins and novel proteins when compared with K_hib proteins in other species. Functional classification, functional enrichment and protein interaction network analyses showed that K_hib proteins are widely distributed in cellular compartments and involved in diverse cellular processes. Significantly, 37 proteins involved in different aspects of regulating the pathogenicity of B. cinerea were detected as K_hib proteins. Our results provide a comprehensive view of the 2-hydroxyisobutyrylome and lay a foundation for further studying the regulatory mechanism of K_hib in both B. cinerea and other plant pathogens.

Enhanced resistance to fungal and bacterial diseases in tomato and Arabidopsis expressing BSR2 from rice

The overexpression of rice BSR2 would offer a simple and effective strategy to protect plants from multiple devastating diseases in tomato and Arabidopsis. Many devastating plant diseases are caused by pathogens possessing a wide host range. Fungal Botrytis cinerea and Rhizoctonia solani, as well as bacterial Pseudomonas syringae and Ralstonia pseudosolanacearum are four such pathogens that infect hundreds of plant species, including agronomically important crops, and cause serious diseases, leading to severe economic losses. However, reports of genes that can confer resistance to broad host-range pathogens via traditional breeding methods are currently limited. We previously reported that Arabidopsis plants overexpressing rice BROAD-SPECTRUM RESISTANCE2 (BSR2/CYP78A15) showed tolerance not only to bacterial P. syringae pv. tomato DC3000 but also to fungal Colletotrichum higginsianum and R. solani. Rice plants overexpressing BSR2 displayed tolerance to two R. solani anastomosis groups. In the present study, first, BSR2-overexpressing (OX) Arabidopsis plants were shown to be additionally tolerant to B. cinerea, R. solani, and R. pseudosolanacearum. Next, tomato 'Micro-Tom' was used as a model to determine whether such tolerance by BSR2 can be introduced into dicot crops to prevent infection from pathogens possessing wide host range. BSR2-OX tomato displayed broad-spectrum disease tolerance to fungal B. cinerea and R. solani, as well as to bacterial P. syringae and R. pseudosolanacearum. Additionally, undesirable traits such as morphological changes were not detected. Thus, BSR2 overexpression can offer a simple and effective strategy to protect crops from multiple destructive diseases.

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.

Phylotranscriptomics of the Pentapetalae Reveals Frequent Regulatory Variation in Plant Local Responses to the Fungal Pathogen Sclerotinia sclerotiorum

Quantitative disease resistance (QDR) is a conserved form of plant immunity that limits infections caused by a broad range of pathogens. QDR has a complex genetic determinism. The extent to which molecular components of the QDR response vary across plant species remains elusive. The fungal pathogen Sclerotinia sclerotiorum, causal agent of white mold diseases on hundreds of plant species, triggers QDR in host populations. To document the diversity of local responses to S. sclerotiorum at the molecular level, we analyzed the complete transcriptomes of six species spanning the Pentapetalae (Phaseolus vulgaris, Ricinus communis, Arabidopsis [Arabidopsis thaliana], Helianthus annuus, Solanum lycopersicum, and Beta vulgaris) inoculated with the same strain of S. sclerotiorum About one-third of plant transcriptomes responded locally to S. sclerotiorum, including a high proportion of broadly conserved genes showing frequent regulatory divergence at the interspecific level. Evolutionary inferences suggested a trend toward the acquisition of gene induction relatively recently in several lineages. Focusing on a group of ABCG transporters, we propose that exaptation by regulatory divergence contributed to the evolution of QDR. This evolutionary scenario has implications for understanding the QDR spectrum and durability. Our work provides resources for functional studies of gene regulation and QDR molecular mechanisms across the Pentapetalae.

Pathogen Genetic Control of Transcriptome Variation in the Arabidopsis thaliana - Botrytis cinerea Pathosystem

In plant-pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host-pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms' transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Genomics of sorghum local adaptation to a parasitic plant

Host-parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR-Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

The polyphagous plant pathogenic fungus Botrytis cinerea encompasses host-specialized and generalist populations

The host plant is often the main variable explaining population structure in fungal plant pathogens, because specialization contributes to reduce gene flow between populations associated with different hosts. Previous population genetic analysis revealed that French populations of the grey mould pathogen Botrytis cinerea were structured by hosts tomato and grapevine, suggesting host specialization in this highly polyphagous pathogen. However, these findings raised questions about the magnitude of this specialization and the possibility of specialization to other hosts. Here we report specialization of B. cinerea populations to tomato and grapevine hosts but not to other tested plants. Population genetic analysis revealed two pathogen clusters associated with tomato and grapevine, while the other clusters co-occurred on hydrangea, strawberry and bramble. Measurements of quantitative pathogenicity were consistent with host specialization of populations found on tomato, and to a lesser extent, populations found on grapevine. Pathogen populations from hydrangea and strawberry appeared to be generalist, while populations from bramble may be weakly specialized. Our results suggest that the polyphagous B. cinerea is more accurately described as a collection of generalist and specialist individuals in populations. This work opens new perspectives for grey mould management, while suggesting spatial optimization of crop organization within agricultural landscapes.

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.