Evolution of the rice blast pathogen on spatially structured rice landraces maintains multiple generalist fungal lineages.

Traditional agrosystems, where humans, crops and microbes have coevolved over long periods, can serve as models to understand the ecoevolutionary determinants of disease dynamics and help the engineering of durably resistant agrosystems. Here, we investigated the genetic and phenotypic relationship between rice (Oryza sativa) landraces and their rice blast pathogen (Pyricularia oryzae) in the traditional Yuanyang terraces of flooded rice paddies in China, where rice landraces have been grown and bred over centuries without significant disease outbreaks. Analyses of genetic subdivision revealed that indica rice plants clustered according to landrace names. Three new diverse lineages of rice blast specific to the Yuanyang terraces coexisted with lineages previously detected at the worldwide scale. Population subdivision in the pathogen population did not mirror pattern of population subdivision in the host. Measuring the pathogenicity of rice blast isolates on landraces revealed generalist life history traits. Our results suggest that the implementation of disease control strategies based on the emergence or maintenance of a generalist lifestyle in pathogens may sustainably reduce the burden of disease in crops.

Convergent Adaptation to Quantitative Host Resistance in a Major Plant Pathogen.

Plant pathogens can adapt to quantitative resistance, eroding its effectiveness. The aim of this work was to reveal the genomic basis of adaptation to such a resistance in populations of the fungus Pseudocercospora fijiensis, a major devastating pathogen of banana, by studying convergent adaptation on different cultivars. Samples from P. fijiensis populations showing a local adaptation pattern on new banana hybrids with quantitative resistance were compared, based on a genome scan approach, with samples from traditional and more susceptible cultivars in Cuba and the Dominican Republic. Whole-genome sequencing of pools of P. fijiensis isolates (pool-seq) sampled from three locations per country was conducted according to a paired population design. The findings of different combined analyses highly supported the existence of convergent adaptation on the study cultivars between locations within but not between countries. Five to six genomic regions involved in this adaptation were detected in each country. An annotation analysis and available biological data supported the hypothesis that some genes within the detected genomic regions may play a role in quantitative pathogenicity, including gene regulation. The results suggested that the genetic basis of fungal adaptation to quantitative plant resistance is at least oligogenic, while highlighting the existence of specific host-pathogen interactions for this kind of resistance.IMPORTANCE Understanding the genetic basis of pathogen adaptation to quantitative resistance in plants has a key role to play in establishing durable strategies for resistance deployment. In this context, a population genomic approach was developed for a major plant pathogen (the fungus Pseudocercospora fijiensis causing black leaf streak disease of banana) whereby samples from new resistant banana hybrids were compared with samples from more susceptible conventional cultivars in two countries. A total of 11 genomic regions for which there was strong evidence of selection by quantitative resistance were detected. An annotation analysis and available biological data supported the hypothesis that some of the genes within these regions may play a role in quantitative pathogenicity. These results suggested a polygenic basis of quantitative pathogenicity in this fungal pathogen and complex molecular plant-pathogen interactions in quantitative disease development involving several genes on both sides.

Spatial and Temporal Analysis of Phytophthora megakarya Epidemic in Newly Established Cacao Plantations.

Studying spatial and temporal plant disease dynamics helps us to understand pathogen dispersal processes and improve disease control recommendations. In this study, three cacao plots devoid of primary inoculum of Phytophthora megakarya (causal agent of cacao black pod rot disease) upon establishment in 2006 were monitored for presence of disease on a weekly basis from 2009 to 2016. Ripley's K(r) function, join count statistics, and Fisher's Exact test were used to analyze spatial and temporal disease dynamics. Disease distribution maps showed aggregated disease patterns in all plots; however, for the years of disease onset, exogenous primary infections were mostly randomly distributed. The K(r) function confirmed these results indicating that inoculum generally disperses only over short distances. Moreover, significant positive spatial autocorrelations showed that diseased trees were often clustered up to a distance of 3 to 9 m. Temporal disease progression was low, meaning that endogenous inoculum failed to establish itself, which is partly explained by rigorous phytosanitation and partly by unfavorable microclimatic conditions for disease development. Because P. megakarya had difficulty establishing itself in the plots, proximity to already infected cacao plantations drove infection dynamics. Thus, isolation of newly established cacao plantations from infected ones and rigorous phytosanitation as a preventive strategy appears to be an effective approach to control cacao black pod rot disease for newly established cacao plantations.

Pattern of local adaptation to quantitative host resistance in a major pathogen of a perennial crop.

Understanding the mechanisms involved in pathogen adaptation to quantitative resistance in plants has a key role to play in establishing durable strategies for resistance deployment, especially in perennial crops. The erosion of quantitative resistance has been recently suspected in Cuba and the Dominican Republic for a major fungal pathogen of such a crop: Pseudocercospora fijiensis, causing black leaf streak disease on banana. This study set out to test whether such erosion has resulted from an adaptation of P. fijiensis populations, and to determine whether or not the adaptation is local. Almost 600 P. fijiensis isolates from Cuba and the Dominican Republic were sampled using a paired-population sampling design on resistant and susceptible banana varieties. A low genetic structure of the P. fijiensis populations was detected in each country using 16 microsatellite markers. Cross-inoculation experiments using isolates from susceptible and resistant cultivars were carried out, measuring a quantitative trait (the diseased leaf area) related to pathogen fitness on three varieties. A further analysis based on those data suggested the existence of a local pattern of adaptation to resistant cultivars in both of the study countries, due to the existence of specific (or genotype by genotype) host-pathogen interactions. However, neither cost nor benefit effects for adapted populations were found on the widely used "Cavendish" banana group. These results highlight the need to study specific host-pathogen interactions and pathogen adaptation on a wide range of quantitative resistance phenotypes in banana, in order to develop durable strategies for resistance deployment.

Multi-scale spatial genetic structure of the vector-borne pathogen 'Candidatus Phytoplasma prunorum' in orchards and in wild habitats.

Inferring the dispersal processes of vector-borne plant pathogens is a great challenge because the plausible epidemiological scenarios often involve complex spread patterns at multiple scales. The spatial genetic structure of 'Candidatus Phytoplasma prunorum', responsible for European stone fruit yellows disease, was investigated by the application of a combination of statistical approaches to genotype data of the pathogen sampled from cultivated and wild compartments in three French Prunus-growing regions. This work revealed that the prevalence of the different genotypes is highly uneven both between regions and compartments. In addition, we identified a significant clustering of similar genotypes within a radius of 50 km or less, but not between nearby wild and cultivated Prunus. We also provide evidence that infected plants are transferred between production areas, and that both species of the Cacopsylla pruni complex can spread the pathogen. Altogether, this work supports a main epidemiological scenario where 'Ca. P. prunorum' is endemic in - and generally acquired from - wild Prunus by its immature psyllid vectors. The latter then migrate to shelter plants that epidemiologically connect sites less than 50 km apart by later providing infectious mature psyllids to their "migration basins". Such multi-scale studies could be useful for other pathosystems.

Specific recognition of two MAX effectors by integrated HMA domains in plant immune receptors involves distinct binding surfaces.

The structurally conserved but sequence-unrelated MAX (Magnaporthe oryzae avirulence and ToxB-like) effectors AVR1-CO39 and AVR-PikD from the blast fungus M. oryzae are recognized by the rice nucleotide-binding domain and leucine-rich repeat proteins (NLRs) RGA5 and Pikp-1, respectively. This involves, in both cases, direct interaction of the effector with a heavy metal-associated (HMA) integrated domain (ID) in the NLR. Here, we solved the crystal structures of a C-terminal fragment of RGA5 carrying the HMA ID (RGA5_S), alone, and in complex with AVR1-CO39 and compared it to the structure of the Pikp1HMA/AVR-PikD complex. In both complexes, HMA ID/MAX effector interactions involve antiparallel alignment of ß-sheets from each partner. However, effector-binding occurs at different surfaces in Pikp1HMA and RGA5HMA, indicating that these interactions evolved independently by convergence of these two MAX effectors to the same type of plant target proteins. Interestingly, the effector-binding surface in RGA5HMA overlaps with the surface that mediates RGA5HMA self-interaction. Mutations in the HMA-binding interface of AVR1-CO39 perturb RGA5HMA-binding, in vitro and in vivo, and affect the recognition of M. oryzae in a rice cultivar containing Pi-CO39 Our study provides detailed insight into the mechanisms of effector recognition by NLRs, which has substantial implications for future engineering of NLRs to expand their recognition specificities. In addition, we propose, as a hypothesis for the understanding of effector diversity, that in the structurally conserved MAX effectors the molecular mechanism of host target protein-binding is conserved rather than the host target proteins themselves.

Pathogen effectors and plant immunity determine specialization of the blast fungus to rice subspecies.

Understanding how fungi specialize on their plant host is crucial for developing sustainable disease control. A traditional, centuries-old rice agro-system of the Yuanyang terraces was used as a model to show that virulence effectors of the rice blast fungus Magnaporthe oryzaeh play a key role in its specialization on locally grown indica or japonica local rice subspecies. Our results have indicated that major differences in several components of basal immunity and effector-triggered immunity of the japonica and indica rice varieties are associated with specialization of M. oryzae. These differences thus play a key role in determining M. oryzae host specificity and may limit the spread of the pathogen within the Yuanyang agro-system. Specifically, the AVR-Pia effector has been identified as a possible determinant of the specialization of M. oryzae to local japonica rice.

Evolution of Compatibility Range in the Rice-Magnaporthe oryzae System: An Uneven Distribution of R Genes Between Rice Subspecies.

Efficient strategies for limiting the impact of pathogens on crops require a good understanding of the factors underlying the evolution of compatibility range for the pathogens and host plants, i.e., the set of host genotypes that a particular pathogen genotype can infect and the set of pathogen genotypes that can infect a particular host genotype. Until now, little is known about the evolutionary and ecological factors driving compatibility ranges in systems implicating crop plants. We studied the evolution of host and pathogen compatibility ranges for rice blast disease, which is caused by the ascomycete Magnaporthe oryzae. We challenged 61 rice varieties from three rice subspecies with 31 strains of M. oryzae collected worldwide from all major known genetic groups. We determined the compatibility range of each plant variety and pathogen genotype and the severity of each plant-pathogen interaction. Compatibility ranges differed between rice subspecies, with the most resistant subspecies selecting for pathogens with broader compatibility ranges and the least resistant subspecies selecting for pathogens with narrower compatibility ranges. These results are consistent with a nested distribution of R genes between rice subspecies.

Surface polysaccharides and quorum sensing are involved in the attachment and survival of Xanthomonas albilineans on sugarcane leaves.

Xanthomonas albilineans, the causal agent of sugarcane leaf scald, is a bacterial plant pathogen that is mainly spread by infected cuttings and contaminated harvesting tools. However, some strains of this pathogen are known to be spread by aerial means and are able to colonize the phyllosphere of sugarcane before entering the host plant and causing disease. The objective of this study was to identify the molecular factors involved in the survival or growth of X. albilineans on sugarcane leaves. We developed a bioassay to test for the attachment of X. albilineans on sugarcane leaves using tissue-cultured plantlets grown in vitro. Six mutants of strain XaFL07-1 affected in surface polysaccharide production completely lost their capacity to survive on the sugarcane leaf surface. These mutants produced more biofilm in vitro and accumulated more cellular poly-ß-hydroxybutyrate than the wild-type strain. A mutant affected in the production of small molecules (including potential biosurfactants) synthesized by non-ribosomal peptide synthetases (NRPSs) attached to the sugarcane leaves as well as the wild-type strain. Surprisingly, the attachment of bacteria on sugarcane leaves varied among mutants of the rpf gene cluster involved in bacterial quorum sensing. Therefore, quorum sensing may affect polysaccharide production, or both polysaccharides and quorum sensing may be involved in the survival or growth of X. albilineans on sugarcane leaves.

Ethylene Response Factors Are Controlled by Multiple Harvesting Stresses in Hevea brasiliensis.

Tolerance of recurrent mechanical wounding and exogenous ethylene is a feature of the rubber tree. Latex harvesting involves tapping of the tree bark and ethephon is applied to increase latex flow. Ethylene is an essential element in controlling latex production. The ethylene signalling pathway leads to the activation of Ethylene Response Factor (ERF) transcription factors. This family has been identified in Hevea brasiliensis. This study set out to understand the regulation of ERF genes during latex harvesting in relation to abiotic stress and hormonal treatments. Analyses of the relative transcript abundance were carried out for 35 HbERF genes in latex, in bark from mature trees and in leaves from juvenile plants under multiple abiotic stresses. Twenty-one HbERF genes were regulated by harvesting stress in laticifers, revealing an overrepresentation of genes in group IX. Transcripts of three HbERF-IX genes from HbERF-IXc4, HbERF-IXc5 and HbERF-IXc6 were dramatically accumulated by combining wounding, methyl jasmonate and ethylene treatments. When an ethylene inhibitor was used, the transcript accumulation for these three genes was halted, showing ethylene-dependent induction. Subcellular localization and transactivation experiments confirmed that several members of HbERF-IX are activator-type transcription factors. This study suggested that latex harvesting induces mechanisms developed for the response to abiotic stress. These mechanisms probably depend on various hormonal signalling pathways. Several members of HbERF-IX could be essential integrators of complex hormonal signalling pathways in Hevea.

Long-distance wind-dispersal of spores in a fungal plant pathogen: estimation of anisotropic dispersal kernels from an extensive field experiment.

Given its biological significance, determining the dispersal kernel (i.e., the distribution of dispersal distances) of spore-producing pathogens is essential. Here, we report two field experiments designed to measure disease gradients caused by sexually- and asexually-produced spores of the wind-dispersed banana plant fungus Mycosphaerella fijiensis. Gradients were measured during a single generation and over 272 traps installed up to 1000 m along eight directions radiating from a traceable source of inoculum composed of fungicide-resistant strains. We adjusted several kernels differing in the shape of their tail and tested for two types of anisotropy. Contrasting dispersal kernels were observed between the two types of spores. For sexual spores (ascospores), we characterized both a steep gradient in the first few metres in all directions and rare long-distance dispersal (LDD) events up to 1000 m from the source in two directions. A heavy-tailed kernel best fitted the disease gradient. Although ascospores distributed evenly in all directions, average dispersal distance was greater in two different directions without obvious correlation with wind patterns. For asexual spores (conidia), few dispersal events occurred outside of the source plot. A gradient up to 12.5 m from the source was observed in one direction only. Accordingly, a thin-tailed kernel best fitted the disease gradient, and anisotropy in both density and distance was correlated with averaged daily wind gust. We discuss the validity of our results as well as their implications in terms of disease diffusion and management strategy.

The variety mixture strategy assessed in a G × G experiment with rice and the blast fungus Magnaporthe oryzae.

Frequent and devastating epidemics of parasites are one of the major issues encountered by modern agriculture. To manage the impact of pathogens, resistant plant varieties have been selected. However, resistances are overcome by parasites requiring the use of pesticides and causing new economical and food safety issues. A promising strategy to maintain the epidemic at a low level and hamper pathogen's adaptation to varietal resistance is the use of mixtures of varieties such that the mix will form a heterogeneous environment for the parasite. A way to find the good combination of varieties that will actually constitute a heterogeneous environment for pathogens is to look for genotype × genotype (G × G) interactions between pathogens and plant varieties. A pattern in which pathogens have a high fitness on one variety and a poor fitness on other varieties guarantees the efficiency of the mixture strategy. In the present article, we inoculated 18 different genotypes of the fungus Magnaporthe oryzae on three rice plant varieties showing different levels of partial resistance in order to find a variety combination compatible with the requirements of the variety mixture strategy, i.e., showing appropriate G × G interactions. We estimated the success of each plant-fungus interaction by measuring fungal fitness and three fungal life history traits: infection success, within-host growth, sporulation capacity. Our results show the existence of G × G interactions between the two varieties Ariete and CO39 on all measured traits and fungal fitness. We also observed that these varieties have different resistance mechanisms; Ariete is good at controlling infection success of the parasite but is not able to control its growth when inside the leaf, while CO39 shows the opposite pattern. We also found that Maratelli's resistance has been eroded. Finally, correlation analyses demonstrated that not all infectious traits are positively correlated.

Characterization of a cassiicolin-encoding gene from Corynespora cassiicola, pathogen of rubber tree (Hevea brasiliensis).

Corynespora Leaf Fall (CLF) is a major disease of rubber tree (Hevea brasiliensis) caused by the Ascomycota Corynespora cassiicola. Here we describe the cloning and characterization of a gene encoding cassiicolin (Cas), a glycosylated cystein-rich small secreted protein (SSP) identified as a potential CLF disease effector in rubber tree. Three isolates with contrasted levels of aggressiveness were analyzed comparatively. The cassiicolin gene was detected - and the toxin successfully purified - from the isolates with high and medium aggressiveness (CCP and CCAM3 respectively) but not from the isolate with the lowest aggressiveness (CCAM1), suggesting the existence of a different disease effector in the later. CCP and CCAM3 carried strictly identical cassiicolin genes and produced toxins of identical mass, as evidence by mass spectrometry analysis, thus suggesting conserved post-translational modifications in addition to sequence identity. The differences in aggressiveness between CCP and CCAM3 may be attributed to differences in cassiicolin transcript levels rather than qualitative variations in cassiicolin structure. Cassiicolin may play an important role in the early phase of infection since a peak of cassiicolin transcripts occurred in 1 or 2 days after inoculation (before the occurrence of the first symptoms), in both the tolerant and the susceptible cultivars.

Gene expression pattern in response to wounding, methyl jasmonate and ethylene in the bark of Hevea brasiliensis.

Natural rubber production in Hevea brasiliensis is determined by both tapping and ethephon frequencies. It is affected by a complex physiological disorder called tapping panel dryness. This syndrome is likely to be induced by environmental and latex harvesting stresses. Defence responses, including rubber biosynthesis, are dramatically mediated by wounding, jasmonate and ethylene (ET), among other factors. Using real-time RT-PCR, the effects of wounding, methyl jasmonate (MeJA) and ET on the relative transcript abundance of a set of 25 genes involved in their signalling and metabolic pathways were studied in the bark of 3-month-old epicormic shoots. Temporal regulation was found for 9 out of 25 genes. Wounding treatment regulated the transcript abundance of 10 genes. Wounding-specific regulation was noted for the HbMAPK, HbBTF3b, HbCAS1, HbLTPP and HbPLD genes. MeJA treatment regulated the transcript abundance of nine genes. Of these, the HbMYB, HbCAS2, HbCIPK and HbChi genes were shown to be specifically MeJA inducible. ET response was accompanied by regulation of the transcript abundance of eight genes, and six genes, HbETR2, HbEIN2, HbEIN3, HbCaM, HbPIP1 and HbQM, were specifically regulated by ET treatment. Additionally, the transcript level of the HbGP and HbACR genes was enhanced by all three treatments simultaneously. Overall, a large number of genes were found to be regulated 4 h after the treatments were applied. This study nevertheless revealed some jasmonic acid-independent wound signalling pathways in H. brasiliensis, provided a general characterization of signalling pathways and will serve as a new base from which to launch advanced studies of the network of pathways operating in H. brasiliensis.

A novel two T-DNA binary vector allows efficient generation of marker-free transgenic plants in three elite cultivars of rice (Oryza sativa L.).

A pilot binary vector was constructed to assess the potential of the 2 T-DNA system for generating selectable marker-free progeny plants in three elite rice cultivars (ZhongZuo321, Ariete and Khao Dawk Mali 105) known to exhibit contrasting amenabilities to transformation. The first T-DNA of the vector, delimited by Agrobacterium tumefaciens borders, contains the hygromycin phosphotransferase (hpt) selectable gene and the green fluorescent protein (gfp) reporter gene while the second T-DNA, delimited by Agrobacterium rhizogenes borders, bears the phosphinothricin acetyl transferase (bar) gene, featuring the gene of interest. 82-90% of the hygromycin-resistant primary transformants exhibited tolerance to ammonium glufosinate mediated by the bar gene suggesting very high co-transformation frequency in the three cultivars. All of the regenerated plants were analyzed by Southern blot which confirmed co-integration of the T-DNAs at frequencies consistent with those of co-expression and allowed determination of copy number for each gene as well as detection of two different vector backbone fragments extending between the two T-DNAs. Hygromycin susceptible, ammonium glufosinate tolerant phenotypes represented 14.4, 17.4 and 14.3% of the plants in T1 progenies of ZZ321, Ariete and KDML105 primary transformants, respectively. We developed a statistical model for deducing from the observed copy number of each T-DNA in T0 plants and phenotypic segregations in T1 progenies the most likely constitution and linkage of the T-DNA integration locus. Statistical analysis identified in 40 out of 42 lines a most likely linkage configuration theoretically allowing genetic separation of the two T-DNA types and out segregation of the T-DNA bearing the bar gene. Overall, though improvements of the technology would be beneficial, the 2 T-DNA system appeared to be a useful approach to generate selectable marker-free rice plants with a consistent frequency among cultivars.