Rapid change in the genetic diversity of Botrytis cinerea populations after the introduction of strains in a tomato glasshouse

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.

Aerobiological modeling I: A review of predictive models

The present work is the first of two reviews on applied modeling in the field of aerobiology. The aerobiological predictive models for pollen and fungal spores, usually defined as predictive statistical models, will, amongst other objectives, forecast airborne particles' concentration or dynamical behavior of the particles. These models can be classified into Observation Based Models (OBM), Phenological Based Models (PHM), or OTher Models (OTM). The aim of this review is to show, analyze and discuss the different predictive models used in pollen and spore aerobiological studies. The analysis was performed on published electronic scientific articles from 1998 to 2016 related to the type of model, the taxa and the modelled parameters. From a total of 503 studies, 55.5% used OBM (44.8% on pollen and 10.7% on fungal spores), 38.5% PHM (all on pollen) and 6% OTM (5.4% on pollen and 0.6% on fungal spores). OBM have been used with high frequency to forecast concentration. The most frequent model of OBM was linear regression (18.5% out of 503) on pollen and artificial neural networks (4.6%) on fungal spores. In the PHM, the principal use was to characterize the main pollen season (flowering season) based on the model of growth degree days. Finally, OTM have been used to estimate concentrations at unmonitored areas. Olea (14,5%) on pollen and Alternaria (4,8%) on fungal spores were the taxa most frequently modelled. Daily concentration was the most modelled parameter by OBM (25.2%) and season start day by PHM (35.6%). The PHM approaches include greater model diversity and use fewer independent variables than OBM. In addition, PHM show to be easier to apply than OBM; however, the wide range of criteria to define the parameters to use in PHM (e.g.: pollination start day) makes that each model is used with a lesser frequency than other models.

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.

Population structure of Venturia inaequalis, a hemibiotrophic fungus, under different host resistance specificities in the Kashmir valley

Venturia inaequalis is a notorious fungal pathogen and show classical gene for gene interaction with its apple host. Neutral markers provide clues about history, evolutionary potential, genetic diversity and population structure of V. inaequalis. The genetic diversity and population structure of fungus indicates that the pathogen is highly diverse with the capacity to breach the scab resistance genes. In the present study, we collected 108 V. inaequalis isolates from three apple cultivars differing in Rvi1 resistance gene. Based on the AMOVA, the variation was mostly distributed among the isolates, providing evidence of non-existence of subpopulation in orchards thus founder population is difficult to arise in Kashmir apple orchards. Pair wise genetic differentiation is less due to regular occurrence of gene flow between the populations residing on different orchard as infected material is transported without stringent quarantine measures. Based on principal coordinate analysis and clustering algorithm as implemented in STRUCTURE, we observed admixture between the two subpopulations, which is quite low, suggesting the existence of pre-zygotic and post-zygotic barriers to gene flow and we cannot rule out the existence of other structures shared by accessions belonging to different varieties. Due to the continuous increase in introduction and monoculture of apple varieties, mixed orchard with different host resistance specificities are more suitable for managing the apple scab in Kashmir valley.

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.

Striking Similarities Between Botrytis cinerea From Non-agricultural and From Agricultural Habitats

Investigations into life history of microorganisms that cause plant diseases have been limited mostly to contexts where they are in interaction with plants, and with cropped or otherwise managed vegetation. Therefore, knowledge about the diversity of plant pathogens, about potential reservoirs of inoculum and about the processes that contribute to their survival and adaptation is limited to these contexts. The agro-centric perspective of plant pathogen life histories is incoherent with respect to the capacity of many of them to persist as saprophytes on various substrates. In this context we have investigated the ubiquity of the broad host range necrotrophic fungus Botrytis cinerea, outside of agricultural settings and have determined if the populations in these natural habitats can be distinguished phenotypically and phylogenetically from populations isolated from diseased crops. Over a period of 5 years, we isolated B. cinerea from 235 samples of various substrates collected in France including rainfall, snowpack, river, and lake water, epilithic biofilms in mountain streams, leaf litter and plant debris, rock surfaces, bird feathers and healthy wild plants from outside of agricultural fields. All substrates except rock surfaces harbored B. cinerea leading us to establish a collection of purified strains that were compared to B. cinerea from diseased tomato, grapes and various other crops in France. Phylogenetic comparisons of 321 strains from crop plants and 100 strains from environmental substrates based on sequences of 9 microsatellite markers revealed that strains from crops and the environment could not be distinguished. Furthermore, the genetic diversity of strains outside of agriculture was just as broad as within agriculture. In tests to determine the aggressiveness of strains on tomato stems, the mean disease severity caused by strains from environmental substrates was statistically identical to the severity of disease caused by strains from tomato, but was significantly greater than the severity caused by strains from grape or other crops. Our results suggest that highly diverse populations of this plant pathogen persist outside of agriculture in association with substrates other than plants and that this part of their life history is compatible with its capacity to maintain its potential as plant 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.

Features of air masses associated with the deposition of Pseudomonas syringae and Botrytis cinerea by rain and snowfall

Clarifying the role of precipitation in microbial dissemination is essential for elucidating the processes involved in disease emergence and spread. The ecology of Pseudomonas syringae and its presence throughout the water cycle makes it an excellent model to address this issue. In this study, 90 samples of freshly fallen rain and snow collected from 2005-2011 in France were analyzed for microbiological composition. The conditions favorable for dissemination of P. syringae by this precipitation were investigated by (i) estimating the physical properties and backward trajectories of the air masses associated with each precipitation event and by (ii) characterizing precipitation chemistry, and genetic and phenotypic structures of populations. A parallel study with the fungus Botrytis cinerea was also performed for comparison. Results showed that (i) the relationship of P. syringae to precipitation as a dissemination vector is not the same for snowfall and rainfall, whereas it is the same for B. cinerea and (ii) the occurrence of P. syringae in precipitation can be linked to electrical conductivity and pH of water, the trajectory of the air mass associated with the precipitation and certain physical conditions of the air mass (i.e. temperature, solar radiation exposure, distance traveled), whereas these predictions are different for B. cinerea. These results are pertinent to understanding microbial survival, emission sources and atmospheric processes and how they influence microbial dissemination.

Remarkable predominance of a small number of genotypes in greenhouse populations of Botrytis cinerea

Although Botrytis cinerea is known for its ability to produce high amounts of spores on diseased plants, enabling it to complete rapidly numerous developmental cycles in favorable environments, population genetics studies of this fungus indicate enormous diversity and limited clonal spread. Here, we report an exception to this situation in the settings of commercial tomato greenhouses. The genotypic characterization of 712 isolates collected from the air and from diseased plants, following the development of gray mold epidemics in four greenhouses in southern France, revealed the presence of a few predominant genotypes in a background of highly diverse populations. The comparison of genotypic profiles for isolates collected in the air or on the plants was compatible with the hypothesis of an entry in the greenhouse of substantial amounts of inoculum from the outside environment but it also highlighted the importance of secondary inoculum produced within the crop. The overall results of this work suggest that sporulation could be an important target for disease management strategies in the greenhouse.

Nitrogen fertilization of the host plant influences production and pathogenicity of Botrytis cinerea secondary inoculum

The influence of nitrogen (N) nutrition on a plant's susceptibility to Botrytis spp. and other pathogens is well documented. However, little is known of possible effects on sporulation of the pathogen on diseased tissue and on the pathogenicity of resulting secondary inoculum. To address this question, sporulation by two strains of Botrytis cinerea was quantified on tomato plants produced under different N irrigation regimes with inputs of NO(3)- at 0.5 to 45 mmol liter(-1) (mM). Sporulation decreased significantly (P < 0.05) with increasing N fertilization up to NO(3)- at 15 to 30 mM. The secondary inoculum was collected and used to inoculate pruning wounds on tomato plants produced under a standard fertilization regime. Pathogenicity of the spores was significantly influenced by the nutritional status of their production substrate. Disease severity was highest with spores produced on plants with very low or very high N fertilization (NO(3)- at 0.5 or 30 mM). It was lowest for inoculum from plants with moderate levels of N fertilization. These results suggest that it may be possible to find an optimum level of N fertilization to reduce the production of secondary inoculum and its pathogenicity to tomato.

Comparison of Botrytis cinerea populations isolated from two open-field cultivated host plants

The necrotrophic fungus Botrytis cinerea is reported to infect more than 220 host plants worldwide. In phylogenetical-taxonomical terms, the pathogen is considered a complex of two cryptic species, group I and group II. We sampled populations of B. cinerea on sympatric strawberry and raspberry cultivars in the North-East of Hungary for three years during flowering and the harvest period. Four hundred and ninety group II B. cinerea isolates were analyzed for the current study. Three different data sets were generated: (i) PCR-RFLP patterns of the ADP-ATP translocase and nitrate reductase genes, (ii) MSB1 minisatellite sequence data, and (iii) the fragment sizes of five microsatellite loci. The structures of the different populations were similar as indicated by Nei's gene diversity and haplotype diversity. The F statistics (Fst, Gst), and the gene flow indicated ongoing differentiation within sympatric populations. The population genetic parameters were influenced by polymorphisms within the three data sets as assessed using Bayesian algorithms. Data Mining analysis pointed towards the five microsatellite loci as the most defining markers to study differentiation in the 490 isolates. The results suggest the occurrence of host-specific, sympatric divergence of generalist phytoparasites in perennial hosts.

Emaravirus: a novel genus of multipartite, negative strand RNA plant viruses

Ringspot symptoms in European mountain ash (Sorbus aucuparia L.), fig mosaic, rose rosette, raspberry leaf blotch, pigeonpea sterility mosaic (Cajanus cajan) and High Plains disease of maize and wheat were found to be associated with viruses that share several characteristics. They all have single-stranded multipartite RNA genomes of negative orientation. In some cases, double membrane-bound virus-like particles of 80 to 200 nm in diameter were found in infected tissue. Furthermore, at least five of these viruses were shown to be vectored by eriophyid mites. Sequences of European mountain ash ringspot-associated virus (EMARaV), Fig mosaic virus (FMV), rose rosette virus (RRV), raspberry leaf blotch virus (RLBV), pigeonpea sterility mosaic virus and High Plains virus strongly support their potential phylogenetic relationship. Therefore, after characterization of EMARaV, the novel genus Emaravirus was established, and FMV was the second virus species assigned to this genus. The recently sequenced RRV and RLBV are supposed to be additional members of this new group of plant RNA viruses.

Differences in Frequency of Transposable Elements Presence in Botrytis cinerea Populations from Several Hosts in Greece

This study was conducted primarily to investigate the presence and frequency distribution of the transposable elements Boty and Flipper in populations of the necrotroph plant pathogen Botrytis cinerea in Greece. In total, 334 isolates were collected from diseased grape, strawberry, tomato, cucumber, kiwifruit, and apple fruit during 2009. The presence of the two transposable elements was based on polymerase chain reaction detection. Results showed that all the sampled hosts occurred in sympatry, with four possible different genotypes (transposa type carrying both transposable elements, Boty type carrying only the Boty element, Flipper type carrying only the Flipper element, and vacuma type carrying neither transposable element). Marked differences in genotype frequencies among populations were observed. In tomato, cucumber, grape, and strawberry, transposa isolates carrying both elements were predominant in the populations whereas, in kiwifruit and apple fruit populations, the vacuma isolates were prevailing. Furthermore, in kiwi and apple fruit populations, high frequencies of Flipper-type isolates were observed. In an attempt to explain the observed predominance of vacuma isolates in kiwifruit populations, the mycelial growth rate of a set of vacuma isolates was compared with the mycelial growth rate of a set of transposa isolates at three different temperatures (0, 10, and 20°C). The same set of isolates was used to compare pathogenicity of isolates on wound-inoculated kiwifruit incubated at two different temperatures (0 and 20°C), in terms of disease incidence and disease severity. In addition, the selected isolates were used to compare their ability in causing latent infections on kiwifruit in the field. The results showed that vacuma and transposa isolates had similar mycelial growth rates at the limiting temperatures of 0 and 10°C, while vacuma isolates grew faster at the optimum temperature of 20°C. Similarly, there was no significant difference regarding pathogenicity on kiwifruit between transposa and vacuma isolates. However, artificial inoculations conducted on blossoms in the field showed that vacuma isolates caused significantly higher incidence of latent infections.

Functional and structural comparison of pyrrolnitrin- and iprodione-induced modifications in the class III histidine-kinase Bos1 of Botrytis cinerea

Dicarboximides and phenylpyrroles are commonly used fungicides against plant pathogenic ascomycetes. Although their effect on fungal osmosensing systems has been shown in many studies, their modes-of-action still remain unclear. Laboratory- or field-mutants of fungi resistant to either or both fungicide categories generally harbour point mutations in the sensor histidine kinase of the osmotic signal transduction cascade.In the present study we compared the mechanisms of resistance to the dicarboximide iprodione and to pyrrolnitrin, a structural analogue of phenylpyrrole fungicides, in Botrytis cinerea. Pyrrolnitrin-induced mutants and iprodione-induced mutants of B. cinerea were produced in vitro. For the pyrrolnitrin-induced mutants, a high level of resistance to pyrrolnitrin was associated with a high level of resistance to iprodione. For the iprodione-induced mutants, the high level of resistance to iprodione generated variable levels of resistance to pyrrolnitrin and phenylpyrroles. All selected mutants showed hypersensitivity to high osmolarity and regardless of their resistance levels to phenylpyrroles, they showed strongly reduced fitness parameters (sporulation, mycelial growth, aggressiveness on plants) compared to the parental phenotypes. Most of the mutants presented modifications in the osmosensing class III histidine kinase affecting the HAMP domains. Site directed mutagenesis of the bos1 gene was applied to validate eight of the identified mutations. Structure modelling of the HAMP domains revealed that the replacements of hydrophobic residues within the HAMP domains generally affected their helical structure, probably abolishing signal transduction. Comparing mutant phenotypes to the HAMP structures, our study suggests that mutations perturbing helical structures of HAMP2-4 abolish signal-transduction leading to loss-of-function phenotype. The mutation of residues E529, M427, and T581, without consequences on HAMP structure, highlighted their involvement in signal transduction. E529 and M427 seem to be principally involved in osmotic signal transduction.