Search

More than 170 years after causing the potato famine in Ireland, late blight is still considered one of the most devastating crop diseases. Commercial potato breeding efforts depend on natural sources of resistance to protect crops from the rapidly evolving late blight pathogen,Phytophthora infestans. We have identified and mapped a novel broad-spectrum disease resistance gene effective againstP. infestansfrom the wild, diploid potato speciesSolanum bulbocastanum. Diagnostic resistance gene enrichment sequencing (dRenSeq) was used to confirm the uniqueness of the identified resistance. RenSeq and GenSeq-based mapping of the resistance, referred to asRpi-blb4, alongside recombinant screening, positioned the locus responsible for the resistance to potato Chromosome 5. The interval spans approximately 2.3 MB and corresponds to the DM reference genome positions of 11.25 and 13.56 Mb.

Potato is the third most important food crop in the world. Diverse pathogens threaten sustainable crop production but can be controlled, in many cases, through the deployment of disease resistance genes belonging to the family of nucleotide-binding, leucine-rich-repeat (NLR) genes.To identify functional NLRs in established varieties, we have successfully established SMRT-AgRenSeq in tetraploid potatoes and have further enhanced the methodology by including dRenSeq in an approach that we term SMRT-AgRenSeq-d. The inclusion of dRenSeq enables the filtering of candidates after the association analysis by establishing a presence/absence matrix across resistant and susceptible potatoes that is translated into an F1 score. Using a SMRT-RenSeq based sequence representation of the NLRome from the cultivar Innovator, SMRT-AgRenSeq-d analyses reliably identified the late blight resistance benchmark genesR1, R2-like, R3aandR3bin a panel of 117 varieties with variable phenotype penetrations. All benchmark genes were identified with an F1 score of 1 which indicates absolute linkage in the panel.When applied to the elusive nematode disease resistance geneGpa5that controls the Potato Cyst Nematode (PCN) speciesGlobodera pallida(pathotypes Pa2/3), SMRT-AgRenSeq-d identified nine strong candidates. These map to the previously established position on potato chromosome 5 and are potential homologs of the late blight resistance geneR1.Assuming that NLRs are involved in controlling many types of resistances, SMRT-AgRenSeq-d can readily be applied to diverse crops and pathogen systems. In potato, SMRT-AgRenSeq-d lends itself, for example, to further study the elusive PCN resistancesH1orH3for which phenotypic data exist.

This work was supported by the Rural & Environment Science & Analytical Services (RESAS) Division of the Scottish Government through project JHI-B1-1, the Biotechnology and Biological Sciences Research Council (BBSRC) through award BB/S015663/1 and the Royal Society through award NAF\R1\201061. MS and LHB were supported through the East of Scotland Bioscience Doctoral Training Partnership (EASTBIO DTP), funded by the BBSRC award BB/T00875X/1. YW was supported through the CSC scholarship program, China. Background : In the ten years since the initial publication of the RenSeq protocol, the method has proved to be a powerful tool for studying disease resistance in plants and providing target genes for breeding programmes. Since the initial publication of the methodology, it has continued to be developed as new technologies have become available and the increased availability of computing power has made new bioinformatic approaches possible. Most recently, this has included the development of a k-mer based association genetics approach, the use of PacBio HiFi data, and graphical genotyping with diagnostic RenSeq. However, there is not yet a unified workflow available and researchers must instead configure approaches from various sources themselves. This makes reproducibility and version control a challenge and limits the ability to perform these analyses to those with bioinformatics expertise. Results : Here we present HISS, consisting of three workflows which take a user from raw RenSeq reads to the identification of candidates for disease resistance genes. These workflows conduct the assembly of enriched HiFi reads from an accession with the resistance phenotype of interest. A panel of accessions both possessing and lacking the resistance are then used in an association genetics approach (AgRenSeq) to identify contigs positively associated with the resistance phenotype. Candidate genes are then identified on these contigs and assessed for their presence or absence in the panel with a graphical genotyping approach that uses dRenSeq. These workflows are implemented via Snakemake, a python-based workflow manager. Software dependencies are either shipped with the release or handled with conda. All code is freely available and is distributed under the GNU GPL-3.0 license. Conclusions : HISS provides a user-friendly, portable, and easily customised approach for identifying novel disease resistance genes in plants. It is easily installed with all dependencies handled internally or shipped with the release and represents a significant improvement in the ease of use of these bioinformatics analyses. Publisher PDF

Plants are resistant to most microbial species due to nonhost resistance (NHR), providing broad-spectrum and durable immunity. However, the molecular components contributing to NHR are poorly characterised. We address the question of whether failure of pathogen effectors to manipulate nonhost plants plays a critical role in NHR. RxLR (Arg-any amino acid-Leu-Arg) effectors from two oomycete pathogens, Phytophthora infestans and Hyaloperonospora arabidopsidis , enhanced pathogen infection when expressed in host plants ( Nicotiana benthamiana and Arabidopsis, respectively) but the same effectors performed poorly in distantly related nonhost pathosystems. Putative target proteins in the host plant potato were identified for 64 P . infestans RxLR effectors using yeast 2-hybrid (Y2H) screens. Candidate orthologues of these target proteins in the distantly related non-host plant Arabidopsis were identified and screened using matrix Y2H for interaction with RxLR effectors from both P . infestans and H . arabidopsidis . Few P . infestans effector-target protein interactions were conserved from potato to candidate Arabidopsis target orthologues (cAtOrths). However, there was an enrichment of H . arabidopsidis RxLR effectors interacting with cAtOrths. We expressed the cAtOrth AtPUB33, which unlike its potato orthologue did not interact with P . infestans effector PiSFI3, in potato and Nicotiana benthamiana. Expression of AtPUB33 significantly reduced P . infestans colonization in both host plants. Our results provide evidence that failure of pathogen effectors to interact with and/or correctly manipulate target proteins in distantly related non-host plants contributes to NHR. Moreover, exploiting this breakdown in effector-nonhost target interaction, transferring effector target orthologues from non-host to host plants is a strategy to reduce disease.

Background: Filopodia are dynamic, finger-like actin-filament bundles that overcome membrane tension by forces generated through actin polymerization at their tips to allow extension of these structures a few microns beyond the cell periphery. Actin assembly of these protrusions is regulated by accessory proteins including heterodimeric capping protein (CP) or Ena/VASP actin polymerases to either terminate or promote filament growth. Accordingly, the depletion of CP in B16-F1 melanoma cells was previously shown to cause an explosive formation of filopodia. In Ena/VASP-deficient cells, CP depletion appeared to result in ruffling instead of inducing filopodia, implying that Ena/VASP proteins are absolutely essential for filopodia formation. However, this hypothesis was not yet experimentally confirmed. Methods: Here, we used B16-F1 cells and CRISPR/Cas9 technology to eliminate CP either alone or in combination with Ena/VASP or other factors residing at filopodia tips, followed by quantifications of filopodia length and number. Results: Unexpectedly, we find massive formations of filopodia even in the absence of CP and Ena/VASP proteins. Notably, combined inactivation of Ena/VASP, unconventional myosin-X and the formin FMNL3 was required to markedly impair filopodia formation in CP-deficient cells. Conclusions: Taken together, our results reveal that, besides Ena/VASP proteins, numerous other factors contribute to filopodia formation.

Sheet-like membrane protrusions at the leading edge, termed lamellipodia, drive 2D-cell migration using active actin polymerization. Microspikes comprise actin-filament bundles embedded within lamellipodia, but the molecular mechanisms driving their formation and their potential functional relevance have remained elusive. Microspike formation requires the specific activity of clustered Ena/VASP proteins at their tips to enable processive actin assembly in the presence of capping protein, but the factors and mechanisms mediating Ena/VASP clustering are poorly understood. Systematic analyses of B16-F1 melanoma mutants lacking potential candidate proteins revealed that neither inverse BAR-domain proteins, nor lamellipodin or Abi are essential for clustering, although they differentially contribute to lamellipodial VASP accumulation. In contrast, unconventional myosin-X (MyoX) identified here as proximal to VASP was obligatory for Ena/VASP clustering and microspike formation. Interestingly, and despite the invariable distribution of other relevant marker proteins, the width of lamellipodia in MyoX-KO mutants was significantly reduced as compared to B16-F1 control, suggesting that microspikes contribute to lamellipodium stability. Consistently, MyoX removal caused marked defects in protrusion and random 2D-cell migration. Strikingly, Ena/VASP-deficiency also uncoupled MyoX cluster dynamics from actin assembly in lamellipodia, establishing their tight functional association in microspike formation.Significance StatementUnlike filopodia that protrude well beyond the cell periphery and are implicated in sensing, morphogenesis and cell-to-cell communication, the function of microspikes consisting of actin-filament bundles fully embedded within lamellipodia is less clear. Microspike formation involves specific clustering of Ena/VASP family members at filament-barbed ends to enable processive actin polymerization in the presence of capping protein, but the factors and mechanisms mediating Ena/VASP clustering have remained unknown. Here, we systematically analyzed these processes in genetic knockout mutants derived from B16-F1 cells and show that Ena/VASP clustering at microspike tips involves Lamellipodin, but not inverse BAR-domain proteins, and strictly requires unconventional Myosin-X. Complete loss of microspikes was confirmed with CRISPR/Cas9-mediated MyoX knockout in Rat2 fibroblasts, excluding cell type-specific effects.

As part of an overall programme aimed at reducing the acrylamide content of crisps, this paper explores the impact of hot-washing on potato slice sugar concentration during industrial scale manufacture. We investigated cold-washing as an alternative to hot-washing, hot-wash residence time and temperature to optimise sugar removal and therefore reduce the potential for high acrylamide levels after frying. Due to the variable nature of potatoes, an extensive variability study was performed to determine confidence boundaries of results. It was found that the cold-wash unit removed on average 21% of the initial sugar content. In the hot-wash the current operational residence time of 3.5 min at 70 °C gave a sugar reduction of 27.5%, which could be increased to 48.5% if residence time is extended to 5 min. Hot-wash temperatures of 40 °C–60 °C were found to increase glucose and fructose content and therefore the potential for acrylamide formation. A “double cold-wash” was trialled and proved to be as successful as hot-washing at 70 °C for all but the highest sugar potatoes, challenging the current operational process and offering the potential for major energy savings.

Late blight is a devastating potato disease worldwide, caused by Phytophthora infestans. The P. infestans strain 2013-18-306 from Yunnan is a “supervirulent race” that overcomes all 11 known late blight resistance genes (R1 to R11) from Solanum demissum. In a previous study, we identified a diploid wild-type potato JAM1-4 (S. jamesii) with high resistance to 2013-18-306. dRenSeq analysis indicated the presence of novel R genes in JAM1-4. RNA-Seq was used to analyze the late blight resistance response genes and defense regulatory mechanisms of JAM1-4 against 2013-18-306. Gene ontology enrichment and KEGG pathway analysis showed that many disease-resistant pathways were significantly enriched. Analysis of differentially expressed genes (DEGs) revealed an active disease resistance mechanism of JAM1-4, and the essential role of multiple signal transduction pathways and secondary metabolic pathways comprised of SA-JA-ET in plant immunity. We also found that photosynthesis in JAM1-4 was inhibited to promote the immune response. Our study reveals the pattern of resistance-related gene expression in response to a super race strain of potato late blight and provides a theoretical basis for further exploration of potato disease resistance mechanisms, discovery of new late blight resistance genes, and disease resistance breeding.

SummaryThe identification of immune receptors in crop plants is time-consuming but important for disease control. Previously, resistance gene enrichment sequencing (RenSeq) was developed to accelerate mapping of nucleotide-binding domain and leucine-rich repeat containing (NLR) genes. However, resistances mediated by pattern recognition receptors (PRRs) remain less utilised.Here, our pipeline shows accelerated mapping of PRRs. Effectoromics leads to precise identification of plants with target PRRs, and subsequent RLP/K enrichment sequencing (RLP/KSeq) leads to detection of informative SNPs that are linked to the trait.Using Phytophthora infestans as a model, we identified Solanum microdontum plants that recognize the apoplastic effectors INF1 or SCR74. RLP/KSeq in a segregating Solanum population confirmed the localization of the INF1 receptor on chromosome 12, and lead to the rapid mapping of the response to SCR74 to chromosome 9. By using markers obtained from RLP/KSeq in conjunction with additional markers, we fine-mapped the SCR74 receptor to a 43-kbp G-LecRK locus.Our findings show that RLP/KSeq enables rapid mapping of PRRs and is especially beneficial for crop plants with large and complex genomes. This work will enable the elucidation and characterisation of the non-NLR plant immune receptors and ultimately facilitate informed resistance breeding.

Late blight, caused by Phytophthora infestans (P. infestans), is a devastating disease in potato worldwide. Our previous study revealed that the Solanum andigena genotype 03112-233 is resistant to P. infestans isolate 90128, but susceptible to the super race isolate, CN152. In this study, we confirmed by diagnostic resistance gene enrichment sequencing (dRenSeq) that the resistance of 03112-233 toward 90128 is most likely based on a distinct new R gene(s). To gain an insight into the mechanism that governs resistance or susceptibility in 03112-223, comparative transcriptomic profiling analysis based on RNAseq was initiated. Changes in transcription at two time points (24 h and 72 h) after inoculation with isolates 90128 or CN152 were analyzed. A total of 8,881 and 7,209 genes were differentially expressed in response to 90128 and CN152, respectively, and 1,083 differentially expressed genes (DEGs) were common to both time points and isolates. A substantial number of genes were differentially expressed in an isolate-specific manner with 3,837 genes showing induction or suppression following infection with 90128 and 2,165 genes induced or suppressed after colonization by CN152. Hierarchical clustering analysis suggested that isolates with different virulence profiles can induce different defense responses at different time points. Further analysis revealed that the compatible interaction caused higher induction of susceptibility genes such as SWEET compared with the incompatible interaction. The salicylic acid, jasmonic acid, and abscisic acid mediated signaling pathways were involved in the response against both isolates, while ethylene and brassinosteroids mediated defense pathways were suppressed. Our results provide a valuable resource for understanding the interactions between P. infestans and potato.

Key MessageNovel major gene resistance againstPotato virus Yin diploid populations ofSolanum tuberosumGroups Phureja and Tuberosum was biologically and genetically characterised. Named Ry(o)phu, it mapped to chromosome 9.AbstractA new source of genetic resistance derived fromSolanum tuberosumGroupPhurejaagainstPotato virus Y(PVY) was identified and genetically characterised in three diploid biparental potato populations. Segregation data for two populations (05H1 and 08H1) suggested the presence of a single dominant gene for resistance to PVY which, following DaRT analysis of the 08H1 cross, was mapped to chromosome 9. More detailed genetic analysis of resistance utilised a well-characterised SNP-linkage map for the 06H1 population, together with newly generated marker data. In these plants, which have bothS. tuberosumGroupPhurejaandS. tuberosumGroupTuberosumin their pedigree, the resistance was shown to map to chromosome 9 at a locus not previously associated with PVY resistance, although there is evidence for at least one other genetic factor controlling PVY infection. The resistance factor location on chromosome 9 (named as Ry(o)phu) suggests a potential role of NB-LRR genes in this resistance. Phenotypic analysis using a GUS-tagged virus revealed that a small amount of PVY replication occurred in occasional groups of epidermal cells in inoculated leaves of resistant plants, without inducing any visible hypersensitive response. However, the virus did not enter the vascular system and systemic spread was completely prevented.

(1) Background: Although resistance to pathogens and pests has been researched in many potato cultivars and breeding lines with DNA markers, there is scarce evidence as to the efficiency of the marker-assisted selection (MAS) for these traits when applied at the early stages of breeding. A goal of this study was to estimate the potential of affordable DNA markers to track resistance genes that are effective against the pathogen Phytophthora infestans (Rpi genes), as a practical breeding tool on a progeny of 68 clones derived from a cross between the cultivar Sudarynya and the hybrid 13/11-09. (2) Methods: this population was studied for four years to elucidate the distribution of late blight (LB) resistance and other agronomical desirable or simple to phenotype traits such as tuber and flower pigmentation, yield capacity and structure. LB resistance was phenotypically evaluated following natural and artificial infection and the presence/absence of nine Rpi genes was assessed with 11 sequence-characterized amplified region (SCAR) markers. To validate this analysis, the profile of Rpi genes in the 13/11-09 parent was established using diagnostic resistance gene enrichment sequencing (dRenSeq) as a gold standard. (3) Results: at the early stages of a breeding program, when screening the segregation of F1 offspring, MAS can halve the workload and selected SCAR markers for Rpi genes provide useful tools.

Three Solanaceae hosts (TSHs), S. tuberosum, N. benthamiana and S. lycopersicum, represent the three major phylogenetic clades of Solanaceae plants infected by Phytophthora infestans, which causes late blight, one of the most devastating diseases seriously affecting crop production. However, details regarding how different Solanaceae hosts respond to P. infestans are lacking. Here, we conducted RNA-seq to analyze the transcriptomic data from the TSHs at 12 and 24 h post P. infestans inoculation to capture early expression effects. Macroscopic and microscopic observations showed faster infection processes in S. tuberosum than in N. benthamiana and S. lycopersicum under the same conditions. Analysis of the number of genes and their level of expression indicated that distinct response models were adopted by the TSHs in response to P. infestans. The host-specific infection process led to overlapping but distinct in GO terms and KEGG pathways enriched for differentially expressed genes, many were tightly linked to the immune response in the TSHs. S. tuberosum showed the fastest response and strongest accumulation of reactive oxygen species compared with N. benthamiana and S. lycopersicum, which also had similarities and differences in hormone regulation. Collectively, our study provides an important reference for a better understanding of late blight response mechanisms of different Solanaceae host interactions.

Diagnostic Resistance gene enrichment Sequencing (dRenSeq) is a method for tracking and verifying the integrity of multiple functional NLR genes in complex genetic backgrounds. Traditionally, these tasks have been addressed by laborious and time-consuming PCR-based approaches or, where cognate avirulence genes are available, effector recognition studies. Both these approaches are superseded by dRenSeq in terms of cost-effectiveness, efficiency, and accuracy.

This work was supported by the Rural & Environment Science & Analytical Services Division of the Scottish Government, the BBSRC, through the joint projects CRF/2009/SCRI/SOP 0929, BB/L008025/1 and BB/K018299/1. Additional funding was obtained through the James Hutton Institute SEEDCORN initiative, AHDB Potato, the Perry Foundation and The Felix Cobbold Trust. Amanpreet Kaur was supported by the Commonwealth Scholarship Commission through a Commonwealth split-site Ph.D. grant. Key message: The nematode resistance gene H2 was mapped to the distal end of chromosome 5 in tetraploid potato. The H2 resistance gene, introduced into cultivated potatoes from the wild diploid species Solanum multidissectum, confers a high level of resistance to the Pa1 pathotype of the potato cyst nematode Globodera pallida. A cross between tetraploid H2-containing breeding clone P55/7 and susceptible potato variety Picasso yielded an F1 population that segregated approximately 1:1 for the resistance phenotype, which is consistent with a single dominant gene in a simplex configuration. Using genome reduction methodologies RenSeq and GenSeq, the segregating F1 population enabled the genetic characterisation of the resistance through a bulked segregant analysis. A diagnostic RenSeq analysis of the parents confirmed that the resistance in P55/7 cannot be explained by previously characterised resistance genes. Only the variety Picasso contained functionally characterised disease resistance genes Rpi-R1, Rpi-R3a, Rpi-R3b variant, Gpa2 and Rx, which was independently confirmed through effector vacuum infiltration assays. RenSeq and GenSeq independently identified sequence polymorphisms linked to the H2 resistance on the top end of potato chromosome 5. Allele-specific KASP markers further defined the locus containing the H2 gene to a 4.7 Mb interval on the distal short arm of potato chromosome 5 and to positions that correspond to 1.4 MB and 6.1 MB in the potato reference genome. Publisher PDF

During fungal infections, plant cells secrete chitinases, which digest chitin in the fungal cell walls. The recognition of released chitin oligomers via lysin motif (LysM)-containing immune host receptors results in the activation of defense signaling pathways. We report here that Verticillium nonalfalfae, a hemibiotrophic xylem-invading fungus, prevents these digestion and recognition processes by secreting a carbohydrate-binding motif 18 (CBM18)-chitin-binding protein, VnaChtBP, which is transcriptionally activated specifically during the parasitic life stages. VnaChtBP is encoded by the Vna8.213 gene, which is highly conserved within the species, suggesting high evolutionary stability and importance for the fungal lifestyle. In a pathogenicity assay, however, Vna8.213 knockout mutants exhibited wilting symptoms similar to the wild-type fungus, suggesting that Vna8.213 activity is functionally redundant during fungal infection of hop. In a binding assay, recombinant VnaChtBP bound chitin and chitin oligomers in vitro with submicromolar affinity and protected fungal hyphae from degradation by plant chitinases. Moreover, the chitin-triggered production of reactive oxygen species from hop suspension cells was abolished in the presence of VnaChtBP, indicating that VnaChtBP also acts as a suppressor of chitin-triggered immunity. Using a yeast-two-hybrid assay, circular dichroism, homology modeling, and molecular docking, we demonstrated that VnaChtBP forms dimers in the absence of ligands and that this interaction is stabilized by the binding of chitin hexamers with a similar preference in the two binding sites. Our data suggest that, in addition to chitin-binding LysM (CBM50) and Avr4 (CBM14) fungal effectors, structurally unrelated CBM18 effectors have convergently evolved to prevent hydrolysis of the fungal cell wall against plant chitinases and to interfere with chitin-triggered host immunity.

Knowledge of the evolutionary processes which govern pathogen recognition is critical to understanding durable disease resistance. We determined how Phytophthora infestans effector PiAVR2 is recognised by evolutionarily distinct resistance proteins R2 and Rpi-mcq1. We employed yeast two-hybrid, co-immunoprecipitation, virus-induced gene silencing, transient overexpression, and phosphatase activity assays to investigate the contributions of BSL phosphatases to R2- and Rpi-mcq1-mediated hypersensitive response (R2 HR and Rpi-mcq1 HR, respectively). Silencing PiAVR2 target BSL1 compromises R2 HR. Rpi-mcq1 HR is compromised only when BSL2 and BSL3 are silenced. BSL1 overexpression increases R2 HR and compromises Rpi-mcq1. However, overexpression of BSL2 or BSL3 enhances Rpi-mcq1 and compromises R2 HR. Okadaic acid, which inhibits BSL phosphatase activity, suppresses both recognition events. Moreover, expression of a BSL1 phosphatase-dead (PD) mutant suppresses R2 HR, whereas BSL2-PD and BSL3-PD mutants suppress Rpi-mcq1 HR. R2 interacts with BSL1 in the presence of PiAVR2, but not with BSL2 and BSL3, whereas no interactions were detected between Rpi-mcq1 and BSLs. Thus, BSL1 activity and association with R2 determine recognition of PiAVR2 by R2, whereas BSL2 and BSL3 mediate Rpi-mcq1 perception of PiAVR2. R2 and Rpi-mcq1 utilise distinct mechanisms to detect PiAVR2 based on association with different BSLs, highlighting central roles of these effector targets for both disease and disease resistance.

Diagnostic Resistance gene enrichment Sequencing (dRenSeq) is a method for tracking and verifying the integrity of multiple functional NLR genes in complex genetic backgrounds. Traditionally, these tasks have been addressed by laborious and time-consuming PCR-based approaches or, where cognate avirulence genes are available, effector recognition studies. Both these approaches are superseded by dRenSeq in terms of cost-effectiveness, efficiency, and accuracy.

This work was funded through an AHDB PhD award, the USDA GLOBAL Project and by the Rural and Environmental Science and Analytical Services division of the Scottish Government. Bioinformatics and computational analyses for the generation of the new G. pallida genome assembly were supported by the University of St Andrews Bioinformatics Unit, which is funded by a Wellcome Trust ISSF award (grant 105621/Z/14/Z). Although the use of natural resistance is the most effective management approach against the potato cyst nematode (PCN) Globodera pallida, the existence of pathotypes with different virulence characteristics constitutes a constraint towards this goal. Two resistance sources, GpaV (from Solanum vernei) and H3 from S. tuberosum ssp. andigena CPC2802 (from the Commonwealth Potato Collection) are widely used in potato breeding programmes in European potato industry. However, the use of resistant cultivars may drive strong selection towards virulence, which allows the increase in frequency of virulent alleles in the population and therefore, the emergence of highly virulent nematode lineages. This study aimed to identify Avirulence (Avr) genes in G. pallida populations selected for virulence on the above resistance sources, and the genomic impact of selection processes on the nematode. The selection drive in the populations was found to be specific to their genetic background. At the genomic level, 11 genes were found that represent candidate Avr genes. Most of the variant calls determining selection were associated with H3-selected populations, while many of them seem to be organised in genomic islands facilitating selection evolution. These phenotypic and genomic findings combined with histological studies performed revealed potential mechanisms underlying selection in G. pallida. Publisher PDF

SummaryPotato late blight, caused by the oomycete pathogenPhytophthora infestans, significantly hampers potato production. Recently, a newResistance to Phytophthora infestans(Rpi) gene,Rpi-amr1, was cloned from a wildSolanumspecies,Solanum americanum. Identification of the corresponding recognized effector (Avirulence, orAvr)genes fromP. infestansis key to elucidating their naturally occurring sequence variation, which in turn informs the potential durability of the cognate late blight resistance.To identify theP. infestanseffector recognized byRpi-amr1, we screened available effector libraries and used long read and cDNA pathogen-enrichment sequencing (PenSeq) on fourP. infestansisolates to explore the untested effectors.By using SMRT and cDNA PenSeq, we identified 47 highly expressed effectors fromP. infestans, including PITG_07569 which triggers a highly specific cell death response when transiently co-expressed withRpi-amr1inNicotiana benthamiana, suggesting thatPITG_07569isAvramr1.Here we demonstrate that long read and cDNA PenSeq enables the identification of full-length RxLR effector families, and their expression profile. This study has revealed key insights into the evolution and polymorphism of a complex RxLR effector family that is associated with the recognition byRpi-amr1.

Late blight is a devastating potato disease worldwide, caused by

Potato late blight, caused by the oomycete pathogen Phytophthora infestans, significantly hampers potato production. Recently, a new Resistance to Phytophthora infestans (Rpi) gene, Rpi‐amr1, was cloned from a wild Solanum species, Solanum americanum. Identification of the corresponding recognized effector (Avirulence or Avr) genes from P. infestans is key to elucidating their naturally occurring sequence variation, which in turn informs the potential durability of the cognate late blight resistance. To identify the P. infestans effector recognized by Rpi‐amr1, we screened available RXLR effector libraries and used long read and cDNA pathogen‐enrichment sequencing (PenSeq) on four P. infestans isolates to explore the untested effectors. Using single‐molecule real‐time sequencing (SMRT) and cDNA PenSeq, we identified 47 highly expressed effectors from P. infestans, including PITG_07569, which triggers a highly specific cell death response when transiently coexpressed with Rpi‐amr1 in Nicotiana benthamiana, suggesting that PITG_07569 is Avramr1. Here we demonstrate that long read and cDNA PenSeq enables the identification of full‐length RXLR effector families and their expression profile. This study has revealed key insights into the evolution and polymorphism of a complex RXLR effector family that is associated with the recognition by Rpi‐amr1., Long reads and cDNA pathogen enrichment sequencing (PenSeq) are developed to reveal the RXLR effector repertoires and expression pattern of four Phytophthora infestans isolates and enable to identify a new avirulence gene Avramr1.

Fifty-two infection-induced Phytophthora infestans RXLR effectors have diverse subcellular localizations, and the majority significantly enhance infection when expressed in leaves. Combining effectors targeting different defence pathways additively enhances colonization., Oomycetes such as the potato blight pathogen Phytophthora infestans deliver RXLR effectors into plant cells to manipulate host processes and promote disease. Knowledge of where they localize inside host cells is important in understanding their function. Fifty-two P. infestans RXLR effectors (PiRXLRs) up-regulated during early stages of infection were expressed as fluorescent protein (FP) fusions inside cells of the model host Nicotiana benthamiana. FP–PiRXLR fusions were predominantly nucleo-cytoplasmic, nuclear, or plasma membrane-associated. Some also localized to the endoplasmic reticulum, mitochondria, peroxisomes, or microtubules, suggesting diverse sites of subcellular activity. Seven of the 25 PiRXLRs examined during infection accumulated at sites of haustorium penetration, probably due to co-localization with host target processes; Pi16663 (Avr1), for example, localized to Sec5-associated mobile bodies which showed perihaustorial accumulation. Forty-five FP–RXLR fusions enhanced pathogen leaf colonization when expressed in Nicotiana benthamiana, revealing that their presence was beneficial to infection. Co-expression of PiRXLRs that target and suppress different immune pathways resulted in an additive enhancement of colonization, indicating the potential to study effector combinations using transient expression assays. We provide a broad platform of high confidence P. infestans effector candidates from which to investigate the mechanisms, singly and in combination, by which this pathogen causes disease.

Summary The oomycete pathogens Phytophthora infestans and P. capsici cause significant crop losses world‐wide, threatening food security. In each case, pathogenicity factors, called RXLR effectors, contribute to virulence. Some RXLRs are perceived by resistance proteins to trigger host immunity, but our understanding of the demographic processes and adaptive evolution of pathogen virulence remains poor.Here, we describe PenSeq, a highly efficient enrichment sequencing approach for genes encoding pathogenicity determinants which, as shown for the infamous potato blight pathogen Phytophthora infestans, make up, See also the Commentary on this article by Kale, 221: 1177–1179.

Late blight, caused by Phytophthora infestans (P. infestans), is a devastating disease in potato worldwide. Our previous study revealed that the Solanum andigena genotype 03112-233 is resistant to P. infestans isolate 90128, but susceptible to the super race isolate, CN152. In this study, we confirmed by diagnostic resistance gene enrichment sequencing (dRenSeq) that the resistance of 03112-233 toward 90128 is most likely based on a distinct new R gene(s). To gain an insight into the mechanism that governs resistance or susceptibility in 03112-223, comparative transcriptomic profiling analysis based on RNAseq was initiated. Changes in transcription at two time points (24 h and 72 h) after inoculation with isolates 90128 or CN152 were analyzed. A total of 8,881 and 7,209 genes were differentially expressed in response to 90128 and CN152, respectively, and 1,083 differentially expressed genes (DEGs) were common to both time points and isolates. A substantial number of genes were differentially expressed in an isolate-specific manner with 3,837 genes showing induction or suppression following infection with 90128 and 2,165 genes induced or suppressed after colonization by CN152. Hierarchical clustering analysis suggested that isolates with different virulence profiles can induce different defense responses at different time points. Further analysis revealed that the compatible interaction caused higher induction of susceptibility genes such as SWEET compared with the incompatible interaction. The salicylic acid, jasmonic acid, and abscisic acid mediated signaling pathways were involved in the response against both isolates, while ethylene and brassinosteroids mediated defense pathways were suppressed. Our results provide a valuable resource for understanding the interactions between P. infestans and potato.

Key message A broad-spectrum late blight disease-resistance gene from Solanum verrucosum has been mapped to potato chromosome 9. The gene is distinct from previously identified-resistance genes. Abstract We have identified and characterised a broad-spectrum resistance to Phytophthora infestans from the wild Mexican species Solanum verrucosum. Diagnostic resistance gene enrichment (dRenSeq) revealed that the resistance is not conferred by previously identified nucleotide-binding, leucine-rich repeat genes. Utilising the sequenced potato genome as a reference, two complementary enrichment strategies that target resistance genes (RenSeq) and single/low-copy number genes (Generic-mapping enrichment Sequencing; GenSeq), respectively, were deployed for the rapid, SNP-based mapping of the resistance through bulked-segregant analysis. Both approaches independently positioned the resistance, referred to as Rpi-ver1, to the distal end of potato chromosome 9. Stringent post-enrichment read filtering identified a total of 64 informative SNPs that corresponded to the expected ratio for significant polymorphisms in the parents as well as the bulks. Of these, 61 SNPs are located on potato chromosome 9 and reside within 27 individual genes, which in the sequenced potato clone DM locate to positions 45.9 to 60.9 Mb. RenSeq- and GenSeq-derived SNPs within the target region were converted into allele-specific PCR-based KASP markers and further defined the position of the resistance to a 4.3 Mb interval at the bottom end of chromosome 9 between positions 52.62–56.98 Mb. Electronic supplementary material The online version of this article (10.1007/s00122-018-3078-6) contains supplementary material, which is available to authorized users.

Targeted capture provides an efficient and sensitive means for sequencing specific genomic regions in a high-throughput manner. To date, this method has mostly been used to capture exons from the genome (the exome) using short insert libraries and short-read sequencing technology, enabling the identification of genetic variants or new members of large gene families. Sequencing larger molecules results in the capture of whole genes, including intronic and intergenic sequences that are typically more polymorphic and allow the resolution of the gene structure of homologous genes, which are often clustered together on the chromosome. Here, we describe an improved method for the capture and single-molecule sequencing of DNA molecules as large as 7 kb by means of size selection and optimized PCR conditions. Our approach can be used to capture, sequence, and distinguish between similar members of the NB-LRR gene family—key genes in plant immune systems.

• Physiological races of the oomycete Albugo candida are biotrophic pathogens of diverse plant species, primarily the Brassicaceae, and cause infections that suppress host immunity to other pathogens. However, A. candida race diversity and the consequences of host immunosuppression are poorly understood in the field. \ud \ud • We report a method that enables sequencing of DNA of plant pathogens and plantassociated microbes directly from field samples (Pathogen Enrichment Sequencing: PenSeq). We apply this method to explore race diversity in A. candida and to detect A. candidaassociated microbes in the field (91 A. candida-infected plants). \ud \ud • We show with unprecedented resolution that each host plant species supports colonization by one of 17 distinct phylogenetic lineages, each with an unique repertoire of effector candidate alleles. These data reveal the crucial role of sexual and asexual reproduction, polyploidy and host domestication in A. candida specialization on distinct plant species. Our bait design also enabled phylogenetic assignment of DNA sequences from bacteria and fungi from plants in the field. \ud \ud • This paper shows that targeted sequencing has a great potential for the study of pathogen populations while they are colonizing their hosts. This method could be applied to other microbes, especially to those that cannot be cultured.

SummaryFollowing the molecular characterisation of functional disease resistance genes in recent years, methods to track and verify the integrity of multiple genes in varieties are needed for crop improvement through resistance stacking. Diagnostic resistance gene enrichment sequencing (dRenSeq) enables the high-confidence identification and complete sequence validation of known functional resistance genes in crops. As demonstrated for tetraploid potato varieties, the methodology is more robust and cost-effective in monitoring resistances than whole-genome sequencing and can be used to appraise (trans)gene integrity efficiently. All currently known NB-LRRs effective against viruses, nematodes and the late blight pathogen Phytophthora infestans can be tracked with dRenSeq in potato and hitherto unknown polymorphisms have been identified. The methodology provides a means to improve the speed and efficiency of future disease resistance breeding in crops by directing parental and progeny selection towards effective combinations of resistance genes.

Plant pathogens deliver effectors to alter host processes. Knowledge of how effectors target and manipulate host proteins is critical to understand crop disease. Here, we show that in planta expression of the RXLR effector Pi04314 enhances leaf colonization by Phytophthora infestans via activity in the host nucleus and attenuates induction of jasmonic and salicylic acid-responsive genes. Pi04314 interacts with three host protein phosphatase 1 catalytic (PP1c) isoforms, causing their re-localization from the nucleolus to the nucleoplasm. Re-localization of PP1c-1 also occurs during infection and is dependent on an R/KVxF motif in the effector. Silencing the PP1c isoforms or overexpression of a phosphatase-dead PP1c-1 mutant attenuates infection, demonstrating that host PP1c activity is required for disease. Moreover, expression of PP1c–1mut abolishes enhanced leaf colonization mediated by in planta Pi04314 expression. We argue that PP1c isoforms are susceptibility factors forming holoenzymes with Pi04314 to promote late blight disease., Phytophtora infestans is the causative agent of late blight disease in plants. Here Boevink et al. provide evidence that host protein phosphatases are targeted by a P. infestans effector protein and are required for disease development.

During fungal infections, plant cells secrete chitinases that digest chitin in the fungal cell walls. The recognition of released chitin oligomers via lysin motif (LysM)-containing immune receptors results in the activation of defence signalling pathways. We report here thatVerticillium nonalfalfae, a hemibiotrophic xylem-invading fungus, prevents this recognition process by secreting a CBM18 (carbohydrate binding motif 18)-chitin binding protein, VnaChtBP, which is transcriptionally activated specifically during the parasitic life stages. VnaChtBP is encoded by theVna8.213gene which is highly conserved within the species, suggesting high evolutionary stability and importance for the fungal lifestyle. In a pathogenicity assay, however,Vna8.213knockout mutants exhibit wilting symptoms similar to the wild type fungus, suggesting thatVna8.213activity is functionally redundant during fungal infection of hop. In binding assay, recombinant VnaChtBP binds chitin and chitin oligomersin vitrowith submicromolar affinity and protects fungal hyphae from degradation by plant chitinases. Using a yeast-two-hybrid assay, homology modelling and molecular docking, we demonstrated that VnaChtBP forms dimers in the absence of ligands and that this interaction is stabilized by the binding of chitin hexamers with a similar preference in the two binding sites. Our data suggest that, in addition to chitin binding LysM (CBM50) and Avr4 (CBM14) fungal effectors, structurally unrelated CBM18 effectors have convergently evolved to prevent hydrolysis of the fungal cell wall against plant chitinases.

Numerous genes that determine the outcome of plant-pathogen interactions are currently being discovered and include, for example, immune receptors, susceptibility factors and pathogen effectors and their host targets. Target enrichment sequencing provides a means to preferentially resequence these genes of interest without the need to first generate a genotype-specific genome assembly. The Basic Local Alignment Search Tool (BLAST), in combination with the here developed BLASTmap, can be used to design probes that specifically target such gene(s), either by using the target species or the closest related genus as a reference. BLAST is a ubiquitous tool in biological sequence analysis and a multitude of programs are available for the visualization of BLAST alignments. However, there are currently no dedicated programs for visual comparison of large-scale BLAST output attributes such as bit score. The need to quickly and efficiently compare many thousands of BLAST results led to the development of BLASTmap, an interactive web application created using the Shiny R package, customized for clustering and viewing BLAST results as an interactive heat map. Here we show an example of how BLASTmap was successfully applied to analyze custom DNA/RNA probe sequences and to visually determine that four probes are sufficient for the specific yet inclusive enrichment of the potato R2 disease resistance gene family.

Numerous genes that determine the outcome of plant-pathogen interactions are currently being discovered and include, for example, immune receptors, susceptibility factors and pathogen effectors and their host targets. Target enrichment sequencing provides a means to preferentially resequence these genes of interest without the need to first generate a genotype-specific genome assembly. The Basic Local Alignment Search Tool (BLAST), in combination with the here developed BLASTmap, can be used to design probes that specifically target such gene(s), either by using the target species or the closest related genus as a reference. BLAST is a ubiquitous tool in biological sequence analysis and a multitude of programs are available for the visualization of BLAST alignments. However, there are currently no dedicated programs for visual comparison of large-scale BLAST output attributes such as bit score. The need to quickly and efficiently compare many thousands of BLAST results led to the development of BLASTmap, an interactive web application created using the Shiny R package, customized for clustering and viewing BLAST results as an interactive heat map. Here we show an example of how BLASTmap was successfully applied to analyze custom DNA/RNA probe sequences and to visually determine that four probes are sufficient for the specific yet inclusive enrichment of the potato R2 disease resistance gene family.

Using genetics, genomics, and diagnostic RenSeq, we demonstrate that the major NB-LRR gene R8 explains the field resistance against potato late blight in the QTL dPI09c, Following the often short-lived protection that major nucleotide binding, leucine-rich-repeat (NB-LRR) resistance genes offer against the potato pathogen Phytophthora infestans, field resistance was thought to provide a more durable alternative to prevent late blight disease. We previously identified the QTL dPI09c on potato chromosome 9 as a more durable field resistance source against late blight. Here, the resistance QTL was fine-mapped to a 186 kb region. The interval corresponds to a larger, 389 kb, genomic region in the potato reference genome of Solanum tuberosum Group Phureja doubled monoploid clone DM1-3 (DM) and from which functional NB-LRRs R8, R9a, Rpi-moc1, and Rpi_vnt1 have arisen independently in wild species. dRenSeq analysis of parental clones alongside resistant and susceptible bulks of the segregating population B3C1HP showed full sequence representation of R8. This was independently validated using long-range PCR and screening of a bespoke bacterial artificial chromosome library. The latter enabled a comparative analysis of the sequence variation in this locus in diverse Solanaceae. We reveal for the first time that broad spectrum and durable field resistance against P. infestans is conferred by the NB-LRR gene R8, which is thought to provide narrow spectrum race-specific resistance.