Solanum americanum genome-assisted discovery of immune receptors that detect potato late blight pathogen effectors.

Potato (Solanum tuberosum) and tomato (Solanum lycopersicon) crops suffer severe losses to late blight caused by the oomycete pathogen Phytophthora infestans. Solanum americanum, a relative of potato and tomato, is globally distributed and most accessions are highly blight resistant. We generated high-quality reference genomes of four S. americanum accessions, resequenced 52 accessions, and defined a pan-NLRome of S. americanum immune receptor genes. We further screened for variation in recognition of 315P. infestans RXLR effectors in 52 S. americanum accessions. Using these genomic and phenotypic data, we cloned three NLR-encoding genes, Rpi-amr4, R02860 and R04373, that recognize cognate P. infestans RXLR effectors PITG_22825 (AVRamr4), PITG_02860 and PITG_04373. These genomic resources and methodologies will support efforts to engineer potatoes with durable late blight resistance and can be applied to diseases of other crops.

A potato late blight resistance gene protects against multiple Phytophthora species by recognizing a broadly conserved RXLR-WY effector.

Species of the genus Phytophthora, the plant killer, cause disease and reduce yields in many crop plants. Although many Resistance to Phytophthora infestans (Rpi) genes effective against potato late blight have been cloned, few have been cloned against other Phytophthora species. Most Rpi genes encode nucleotide-binding domain, leucine-rich repeat-containing (NLR) immune receptor proteins that recognize RXLR (Arg-X-Leu-Arg) effectors. However, whether NLR proteins can recognize RXLR effectors from multiple Phytophthora species has rarely been investigated. Here, we identified a new RXLR-WY effector AVRamr3 from P. infestans that is recognized by Rpi-amr3 from a wild Solanaceae species Solanum americanum. Rpi-amr3 associates with AVRamr3 in planta. AVRamr3 is broadly conserved in many different Phytophthora species, and the recognition of AVRamr3 homologs by Rpi-amr3 activates resistance against multiple Phytophthora pathogens, including the tobacco black shank disease and cacao black pod disease pathogens P. parasitica and P. palmivora. Rpi-amr3 is thus the first characterized resistance gene that acts against P. parasitica or P. palmivora. These findings suggest a novel path to redeploy known R genes against different important plant pathogens.

A complex resistance locus in Solanum americanum recognizes a conserved Phytophthora effector.

Late blight caused by Phytophthora infestans greatly constrains potato production. Many Resistance (R) genes were cloned from wild Solanum species and/or introduced into potato cultivars by breeding. However, individual R genes have been overcome by P. infestans evolution; durable resistance remains elusive. We positionally cloned a new R gene, Rpi-amr1, from Solanum americanum, that encodes an NRC helper-dependent CC-NLR protein. Rpi-amr1 confers resistance in potato to all 19 P. infestans isolates tested. Using association genomics and long-read RenSeq, we defined eight additional Rpi-amr1 alleles from different S. americanum and related species. Despite only ~90% identity between Rpi-amr1 proteins, all confer late blight resistance but differentially recognize Avramr1 orthologues and paralogues. We propose that Rpi-amr1 gene family diversity assists detection of diverse paralogues and alleles of the recognized effector, facilitating durable resistance against P. infestans.

Identification of RipAZ1 as an avirulence determinant of Ralstonia solanacearum in Solanum americanum.

Ralstonia solanacearum causes bacterial wilt disease in many plant species. Type III-secreted effectors (T3Es) play crucial roles in bacterial pathogenesis. However, some T3Es are recognized by corresponding disease resistance proteins and activate plant immunity. In this study, we identified the R. solanacearum T3E protein RipAZ1 (Ralstonia injected protein AZ1) as an avirulence determinant in the black nightshade species Solanum americanum. Based on the S. americanum accession-specific avirulence phenotype of R. solanacearum strain Pe_26, 12 candidate avirulence T3Es were selected for further analysis. Among these candidates, only RipAZ1 induced a cell death response when transiently expressed in a bacterial wilt-resistant S. americanum accession. Furthermore, loss of ripAZ1 in the avirulent R. solanacearum strain Pe_26 resulted in acquired virulence. Our analysis of the natural sequence and functional variation of RipAZ1 demonstrated that the naturally occurring C-terminal truncation results in loss of RipAZ1-triggered cell death. We also show that the 213 amino acid central region of RipAZ1 is sufficient to induce cell death in S. americanum. Finally, we show that RipAZ1 may activate defence in host cell cytoplasm. Taken together, our data indicate that the nucleocytoplasmic T3E RipAZ1 confers R. solanacearum avirulence in S. americanum. Few avirulence genes are known in vascular bacterial phytopathogens and ripAZ1 is the first one in R. solanacearum that is recognized in black nightshades. This work thus opens the way for the identification of disease resistance genes responsible for the specific recognition of RipAZ1, which can be a source of resistance against the devastating bacterial wilt disease.

Identification of Avramr1 from Phytophthora infestans using long read and cDNA pathogen-enrichment sequencing (PenSeq).

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.

Evolution of NLR resistance genes with noncanonical N-terminal domains in wild tomato species.

Nucleotide-binding and leucine-rich repeat immune receptors (NLRs) provide resistance against diverse pathogens. To create comparative NLR resources, we conducted resistance gene enrichment sequencing (RenSeq) with single-molecule real-time sequencing of PacBio for 18 accessions in Solanaceae, including 15 accessions of five wild tomato species. We investigated the evolution of a class of NLRs, CNLs with extended N-terminal sequences previously named Solanaceae Domain. Through comparative genomic analysis, we revealed that the extended CNLs (exCNLs) anciently emerged in the most recent common ancestor between Asterids and Amaranthaceae, far predating the Solanaceae family. In tomatoes, the exCNLs display exceptional modes of evolution in a clade-specific manner. In the clade G3, exCNLs have substantially elongated their N-termini through tandem duplications of exon segments. In the clade G1, exCNLs have evolved through recent proliferation and sequence diversification. In the clade G6, an ancestral exCNL has lost its N-terminal domains in the course of evolution. Our study provides high-quality NLR gene models for close relatives of domesticated tomatoes that can serve as a useful resource for breeding and molecular engineering for disease resistance. Our findings regarding the exCNLs offer unique backgrounds and insights for future functional studies of the NLRs.

Extreme resistance to Potato virus Y in potato carrying the Rysto gene is mediated by a TIR-NLR immune receptor.

Potato virus Y (PVY) is a major potato (Solanum tuberosum L.) pathogen that causes severe annual crop losses worth billions of dollars worldwide. PVY is transmitted by aphids, and successful control of virus transmission requires the extensive use of environmentally damaging insecticides to reduce vector populations. Rysto , from the wild relative S. stoloniferum, confers extreme resistance (ER) to PVY and related viruses and is a valuable trait that is widely employed in potato resistance breeding programmes. Rysto was previously mapped to a region of potato chromosome XII, but the specific gene has not been identified to date. In this study, we isolated Rysto using resistance gene enrichment sequencing (RenSeq) and PacBio SMRT (Pacific Biosciences single-molecule real-time sequencing). Rysto was found to encode a nucleotide-binding leucine-rich repeat (NLR) protein with an N-terminal TIR domain and was sufficient for PVY perception and ER in transgenic potato plants. Rysto -dependent extreme resistance was temperature-independent and requires EDS1 and NRG1 proteins. Rysto may prove valuable for creating PVY-resistant cultivars of potato and other Solanaceae crops.

Gene expression polymorphism underpins evasion of host immunity in an asexual lineage of the Irish potato famine pathogen.

BACKGROUND: Outbreaks caused by asexual lineages of fungal and oomycete pathogens are a continuing threat to crops, wild animals and natural ecosystems (Fisher MC, Henk DA, Briggs CJ, Brownstein JS, Madoff LC, McCraw SL, Gurr SJ, Nature 484:186-194, 2012; Kupferschmidt K, Science 337:636-638, 2012). However, the mechanisms underlying genome evolution and phenotypic plasticity in asexual eukaryotic microbes remain poorly understood (Seidl MF, Thomma BP, BioEssays 36:335-345, 2014). Ever since the 19th century Irish famine, the oomycete Phytophthora infestans has caused recurrent outbreaks on potato and tomato crops that have been primarily caused by the successive rise and migration of pandemic asexual lineages (Goodwin SB, Cohen BA, Fry WE, Proc Natl Acad Sci USA 91:11591-11595, 1994; Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10:e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Cooke DEL, Cano LM, Raffaele S, Bain RA, Cooke LR, Etherington GJ, Deahl KL, Farrer RA, Gilroy EM, Goss EM, et al. PLoS Pathog 8:e1002940, 2012). However, the dynamics of genome evolution within these clonal lineages have not been determined. The objective of this study was to use a comparative genomics and transcriptomics approach to determine the molecular mechanisms that underpin phenotypic variation within a clonal lineage of P. infestans. RESULTS: Here, we reveal patterns of genomic and gene expression variation within a P. infestans asexual lineage by comparing strains belonging to the South American EC-1 clone that has dominated Andean populations since the 1990s (Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Delgado RA, Monteros-Altamirano AR, Li Y, Visser RGF, van der Lee TAJ, Vosman B, Plant Pathol 62:1081-1088, 2013; Forbes GA, Escobar XC, Ayala CC, Revelo J, Ordonez ME, Fry BA, Doucett K, Fry WE, Phytopathology 87:375-380, 1997; Oyarzun PJ, Pozo A, Ordonez ME, Doucett K, Forbes GA, Phytopathology 88:265-271, 1998). We detected numerous examples of structural variation, nucleotide polymorphisms and loss of heterozygosity within the EC-1 clone. Remarkably, 17 genes are not expressed in one of the two EC-1 isolates despite apparent absence of sequence polymorphisms. Among these, silencing of an effector gene was associated with evasion of disease resistance conferred by a potato immune receptor. CONCLUSIONS: Our findings highlight the molecular changes underpinning the exceptional genetic and phenotypic plasticity associated with host adaptation in a pandemic clonal lineage of a eukaryotic plant pathogen. We observed that the asexual P. infestans lineage EC-1 can exhibit phenotypic plasticity in the absence of apparent genetic mutations resulting in virulence on a potato carrying the Rpi-vnt1.1 gene. Such variant alleles may be epialleles that arose through epigenetic changes in the underlying genes.

Identification and rapid mapping of a gene conferring broad-spectrum late blight resistance in the diploid potato species Solanum verrucosum through DNA capture technologies.

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. 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.

Targeted capture and sequencing of gene-sized DNA molecules.

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.

Accelerated cloning of a potato late blight-resistance gene using RenSeq and SMRT sequencing.

Global yields of potato and tomato crops have fallen owing to potato late blight disease, which is caused by Phytophthora infestans. Although most commercial potato varieties are susceptible to blight, many wild potato relatives show variation for resistance and are therefore a potential source of Resistance to P. infestans (Rpi) genes. Resistance breeding has exploited Rpi genes from closely related tuber-bearing potato relatives, but is laborious and slow. Here we report that the wild, diploid non-tuber-bearing Solanum americanum harbors multiple Rpi genes. We combine resistance (R) gene sequence capture (RenSeq) with single-molecule real-time (SMRT) sequencing (SMRT RenSeq) to clone Rpi-amr3i. This technology should enable de novo assembly of complete nucleotide-binding, leucine-rich repeat receptor (NLR) genes, their regulatory elements and complex multi-NLR loci from uncharacterized germplasm. SMRT RenSeq can be applied to rapidly clone multiple R genes for engineering pathogen-resistant crops.

Defining the full tomato NB-LRR resistance gene repertoire using genomic and cDNA RenSeq.

BACKGROUND: The availability of draft crop plant genomes allows the prediction of the full complement of genes that encode NB-LRR resistance gene homologs, enabling a more targeted breeding for disease resistance. Recently, we developed the RenSeq method to reannotate the full NB-LRR gene complement in potato and to identify novel sequences that were not picked up by the automated gene prediction software. Here, we established RenSeq on the reference genome of tomato (Solanum lycopersicum) Heinz 1706, using 260 previously identified NB-LRR genes in an updated Solanaceae RenSeq bait library. RESULT: Using 250-bp MiSeq reads after RenSeq on genomic DNA of Heinz 1706, we identified 105 novel NB-LRR sequences. Reannotation included the splitting of gene models, combination of partial genes to a longer sequence and closing of assembly gaps. Within the draft S. pimpinellifolium LA1589 genome, RenSeq enabled the annotation of 355 NB-LRR genes. The majority of these are however fragmented, with 5'- and 3'-end located on the edges of separate contigs. Phylogenetic analyses show a high conservation of all NB-LRR classes between Heinz 1706, LA1589 and the potato clone DM, suggesting that all sub-families were already present in the last common ancestor. A phylogenetic comparison to the Arabidopsis thaliana NB-LRR complement verifies the high conservation of the more ancient CCRPW8-type NB-LRRs. Use of RenSeq on cDNA from uninfected and late blight-infected tomato leaves allows the avoidance of sequence analysis of non-expressed paralogues. CONCLUSION: RenSeq is a promising method to facilitate analysis of plant resistance gene complements. The reannotated tomato NB-LRR complements, phylogenetic relationships and chromosomal locations provided in this paper will provide breeders and scientists with a useful tool to identify novel disease resistance traits. cDNA RenSeq enables for the first time next-gen sequencing approaches targeted to this very low-expressed gene family without the need for normalization.

Elevating crop disease resistance with cloned genes.

Essentially all plant species exhibit heritable genetic variation for resistance to a variety of plant diseases caused by fungi, bacteria, oomycetes or viruses. Disease losses in crop monocultures are already significant, and would be greater but for applications of disease-controlling agrichemicals. For sustainable intensification of crop production, we argue that disease control should as far as possible be achieved using genetics rather than using costly recurrent chemical sprays. The latter imply CO2 emissions from diesel fuel and potential soil compaction from tractor journeys. Great progress has been made in the past 25 years in our understanding of the molecular basis of plant disease resistance mechanisms, and of how pathogens circumvent them. These insights can inform more sophisticated approaches to elevating disease resistance in crops that help us tip the evolutionary balance in favour of the crop and away from the pathogen. We illustrate this theme with an account of a genetically modified (GM) blight-resistant potato trial in Norwich, using the Rpi-vnt1.1 gene isolated from a wild relative of potato, Solanum venturii, and introduced by GM methods into the potato variety Desiree.

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations.

RenSeq is a NB-LRR (nucleotide binding-site leucine-rich repeat) gene-targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB-LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB-LRRs and can be accessed through a genome browser that we provide. We compared these NB-LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum 'Heinz 1706' extended the NB-LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co-segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi-ber2) and S. ruiz-ceballosii (Rpi-rzc1), we were able to apply RenSeq successfully to identify markers that co-segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy-to-adapt Galaxy pipelines.

Structures of an active-site mutant of a plant 1,3-ß-glucanase in complex with oligosaccharide products of hydrolysis.

Plant endo-1,3-ß-glucanases are involved in important physiological processes such as defence mechanisms, cell division and flowering. They hydrolyze (1→3)-ß-glucans, with very limited activity towards mixed (1→3,1→4)-ß-glucans and branched (1→3,1→6)-ß-glucans. Here, crystal structures of the potato (Solanum tuberosum) endo-1,3-ß-glucanase GLUB20-2 with the nucleophilic Glu259 residue substituted by alanine (E259A) are reported. Despite this active-site mutation, the protein retained residual endoglucanase activity and when incubated in the crystallization buffer with a linear hexameric substrate derived from (1→3)-ß-glucan (laminarahexose) cleaved it in two different ways, generating trisaccharides and tetrasaccharides, as confirmed by mass spectrometry. The trisaccharide (laminaratriose) shows higher binding affinity and was found to fully occupy the -1, -2 and -3 sites of the active-site cleft, even at a low molar excess of the substrate. At elevated substrate concentration the tetrasaccharide molecule (laminaratetrose) also occupies the active site, spanning the opposite sites +1, +2, +3 and +4 of the cleft. These are the first crystal structures of a plant glycoside hydrolase family 17 (GH17) member to reveal the protein-saccharide interactions and were determined at resolutions of 1.68 and 1.55 Å, respectively. The geometry of the active-site cleft clearly precludes any (1→4)-ß-glucan topology at the subsites from -3 to +4 and could possibly accommodate ß-1,6-branching only at subsites +1 and +2. The glucose units at subsites -1 and -2 interact with highly conserved protein residues. In contrast, subsites -3, +3 and +4 are variable, suggesting that the mode of glucose binding at these sites may vary between different plant endo-1,3-ß-glucanases. Low substrate affinity is observed at subsites +1 and +2, as manifested by disorder of the glycosyl units there.

Two high-resolution structures of potato endo-1,3-ß-glucanase reveal subdomain flexibility with implications for substrate binding.

Endo-1,3-ß-glucanases are widely distributed among bacteria, fungi and higher plants. They are responsible for hydrolysis of the glycosidic bond in specific polysaccharides with tracts of unsubstituted ß-1,3-linked glucosyl residues. The plant enzymes belong to glycoside hydrolase family 17 (GH17) and are also members of class 2 of pathogenesis-related (PR) proteins. X-ray diffraction data were collected to 1.40 and 1.26 Å resolution from two crystals of endo-1,3-ß-glucanase from Solanum tuberosum (potato, cultivar Désirée) which, despite having a similar packing framework, represented two separate crystal forms. In particular, they differed in the Matthews coefficient and are consequently referred to as higher density (HD; 1.40 Å resolution) and lower density (LD; 1.26 Å resolution) forms. The general fold of the protein resembles that of other known plant endo-1,3-ß-glucanases and is defined by a (ß/α)(8)-barrel with an additional subdomain built around the C-terminal half of the barrel. The structures revealed high flexibility of the subdomain, which forms part of the catalytic cleft. Comparison with structures of other GH17 endo-1,3-ß-glucanases revealed differences in the arrangement of the secondary-structure elements in this region, which can be correlated with sequence variability and may suggest distinct substrate-binding patterns. The crystal structures revealed an unusual packing mode, clearly visible in the LD structure, caused by the presence of the C-terminal His(6) tag, which extends from the compact fold of the enzyme molecule and docks in the catalytic cleft of a neighbouring molecule. In this way, an infinite chain of His-tag-linked protein molecules is formed along the c direction.

Conserved Cys residue influences catalytic properties of potato endo-(1-->3)-beta-glucanase GLUB20-2.

The synthesis and degradation of (1-->3)-beta-glycosidic bonds between glucose moieties are essential metabolic processes in plant cell architecture and function. We have found that a unique, conserved cysteine residue, positioned outside the catalytic centre of potato endo-(1-->3)-beta-glucanase - product of the gluB20-2 gene, participates in determining the substrate specificity of the enzyme. The same residue is largely responsible for endo-(1-->3)-beta-glucanase inhibition by mercury ions. Our results confirm that the spatial adjustment between an enzyme and its substrate is one of the essential factors contributing to the specificity and accuracy of enzymatic reactions.