A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Widespread premature transcription termination of Arabidopsis thaliana NLR genes by the spen protein FPA

Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Panoramic: A package for constructing eukaryotic pan-genomes

The study of intraspecific genomic variation in eukaryotic species has been the focus of numerous genome resequencing projects in recent years. One emerging approach for the analysis of intraspecific diversity uses the concept of a pan-genome, which theoretically represents the full set of genomic sequences and coding genes from all individuals of a given species. This approach has many advantages over reference-based methods and has been successfully applied to study both prokaryotic and eukaryotic species. However, the process of pan-genome construction still presents considerable scientific and technical challenges, especially for eukaryotic species with large and complex genomes. Although general approaches for the construction of pan-genomes have been devised, currently available software tools implement only certain modules of the entire computational procedure. Therefore, each pan-genome project requires the development of tailored analysis pipelines, thus complicating and prolonging the process and impairing research reproducibility and comparison across studies. Here, we present Panoramic, a software package for the automatic construction of eukaryotic pan-genomes. Panoramic takes raw sequencing reads as input and applies two alternative approaches for pan-genome construction. Panoramic makes pan-genome construction a considerably easier task by providing simple user interface and efficient data processing algorithms. We demonstrate the use of Panoramic by constructing the pan-genome of the model plant species Arabidopsis thaliana from sequencing data of 20 diverse ecotypes.

A high-quality Brassica napus genome reveals expansion of transposable elements, subgenome evolution and disease resistance

Rapeseed (Brassica napus L.) is a recent allotetraploid crop, which is well known for its high oil production. Here, we report a high-quality genome assembly of a typical semi-winter rapeseed cultivar, 'Zhongshuang11' (hereafter 'ZS11'), using a combination of single-molecule sequencing and chromosome conformation capture (Hi-C) techniques. Most of the high-confidence sequences (93.1%) were anchored to the individual chromosomes with a total of 19 centromeres identified, matching the exact chromosome count of B. napus. The repeat sequences in the A and C subgenomes in B. napus expanded significantly from 500 000 years ago, especially over the last 100 000 years. These young and recently amplified LTR-RTs showed dispersed chromosomal distribution but significantly preferentially clustered into centromeric regions. We exhaustively annotated the nucleotide-binding leucine-rich repeat (NLR) gene repertoire, yielding a total of 597 NLR genes in B. napus genome and 17.4% of which are paired (head-to-head arrangement). Based on the resequencing data of 991 B. napus accessions, we have identified 18 759 245 single nucleotide polymorphisms (SNPs) and detected a large number of genomic regions under selective sweep among the three major ecotype groups (winter, semi-winter and spring) in B. napus. We found 49 NLR genes and five NLR gene pairs colocated in selective sweep regions with different ecotypes, suggesting a rapid diversification of NLR genes during the domestication of B. napus. The high quality of our B. napus 'ZS11' genome assembly could serve as an important resource for the study of rapeseed genomics and reveal the genetic variations associated with important agronomic traits.

Hybrid Incompatibility of the Plant Immune System: An Opposite Force to Heterosis Equilibrating Hybrid Performances

Hybridization is a core element in modern rice breeding as beneficial combinations of two parental genomes often result in the expression of heterosis. On the contrary, genetic incompatibility between parents can manifest as hybrid necrosis, which leads to tissue necrosis accompanied by compromised growth and/or reduced reproductive success. Genetic and molecular studies of hybrid necrosis in numerous plant species revealed that such self-destructing symptoms in most cases are attributed to autoimmunity: plant immune responses are inadvertently activated in the absence of pathogenic invasion. Autoimmunity in hybrids predominantly occurs due to a conflict involving a member of the major plant immune receptor family, the nucleotide-binding domain and leucine-rich repeat containing protein (NLR; formerly known as NBS-LRR). NLR genes are associated with disease resistance traits, and recent population datasets reveal tremendous diversity in this class of immune receptors. Cases of hybrid necrosis involving highly polymorphic NLRs as major causes suggest that diversified R gene repertoires found in different lineages would require a compatible immune match for hybridization, which is a prerequisite to ensure increased fitness in the resulting hybrids. In this review, we overview recent genetic and molecular findings on hybrid necrosis in multiple plant species to provide an insight on how the trade-off between growth and immunity is equilibrated to affect hybrid performances. We also revisit the cases of hybrid weakness in which immune system components are found or implicated to play a causative role. Based on our understanding on the trade-off, we propose that the immune system incompatibility in plants might play an opposite force to restrict the expression of heterosis in hybrids. The antagonism is illustrated under the plant fitness equilibrium, in which the two extremes lead to either hybrid necrosis or heterosis. Practical proposition from the equilibrium model is that breeding efforts for combining enhanced disease resistance and high yield shall be achieved by balancing the two forces. Reverse breeding toward utilizing genomic data centered on immune components is proposed as a strategy to generate elite hybrids with balanced immunity and growth.

Genome-Wide Characterization of Lectin Receptor Kinases in Saccharum spontaneum L. and Their Responses to Stagonospora tainanensis Infection

Sugarcane is an important sugar and bioenergy ethanol crop, and the hyperploidy has led to stagnant progress in sugarcane genome decipherment, which also hindered the genome-wide analyses of versatile lectin receptor kinases (LecRKs). The published genome of Saccharum spontaneum, one of the two sugarcane ancestor species, enables us to study the characterization of LecRKs and their responses to sugarcane leaf blight (SLB) triggered by Stagonospora tainanensis. A total of 429 allelic and non-allelic LecRKs, which were classified into evolved independently three types according to signal domains and phylogeny, were identified based on the genome. Regarding those closely related LecRKs in the phylogenetic tree, their motifs and exon architectures of representative L- and G-types were similar or identical. LecRKs showed an unequal distribution on chromosomes and more G-type tandem repeats may come from the gene expansion. Comparing the differentially expressed LecRKs (DELs) in response to SLB in sugarcane hybrid and ancestor species S. spontaneum, we found that the DEL number in the shared gene sets was highly variable among each sugarcane accession, which indicated that the expression dynamics of LecRKs in response to SLB were quite different between hybrids and particularly between sugarcane hybrid and S. spontaneum. In addition, C-type LecRKs may participate in metabolic processes of plant–pathogen interaction, mainly including pathogenicity and plant resistance, indicating their putative roles in sugarcane responses to SLB infection. The present study provides a basic reference and global insight into the further study and utilization of LecRKs in plants.

Diversity, structure and function of the coiled-coil domains of plant NLR immune receptors

Plant nucleotide-binding, leucine-rich repeat receptors (NLRs) perceive pathogen avirulence effectors and activate defense responses. Nucleotide-binding, leucine-rich repeat receptors are classified into coiled-coil (CC)-containing and Toll/interleukin-1 receptor (TIR)-containing NLRs. Recent advances suggest that NLR CC domains often function in signaling activation, especially for induction of cell death. In this review, we outline our current understanding of NLR CC domains, including their diversity/classification and structure, their roles in cell death induction, disease resistance, and interaction with other proteins. Furthermore, we provide possible directions for future work.

A natural diversity screen in Arabidopsis thaliana reveals determinants for HopZ1a recognition in the ZAR1-ZED1 immune complex

Pathogen pressure on hosts can lead to the evolution of genes regulating the innate immune response. By characterizing naturally occurring polymorphisms in immune receptors, we can better understand the molecular determinants of pathogen recognition. ZAR1 is an ancient Arabidopsis thaliana NLR (Nucleotide-binding [NB] Leucine-rich-repeat [LRR] Receptor) that recognizes multiple secreted effector proteins from the pathogenic bacteria Pseudomonas syringae and Xanthomonas campestris through its interaction with receptor-like cytoplasmic kinases (RLCKs). ZAR1 was first identified for its role in recognizing P. syringae effector HopZ1a, through its interaction with the RLCK ZED1. To identify additional determinants of HopZ1a recognition, we performed a computational screen for ecotypes from the 1001 Genomes project that were likely to lack HopZ1a recognition, and tested ~300 ecotypes. We identified ecotypes containing polymorphisms in ZAR1 and ZED1. Using our previously established Nicotiana benthamiana transient assay and Arabidopsis ecotypes, we tested for the effect of naturally occurring polymorphisms on ZAR1 interactions and the immune response. We identified key residues in the NB or LRR domain of ZAR1 that impact the interaction with ZED1. We demonstrate that natural diversity combined with functional assays can help define the molecular determinants and interactions necessary to regulate immune induction in response to pathogens.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

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.

A Truncated Singleton NLR Causes Hybrid Necrosis in Arabidopsis thaliana

Hybrid necrosis in plants arises from conflict between divergent alleles of immunity genes contributed by different parents, resulting in autoimmunity. We investigate a severe hybrid necrosis case in Arabidopsis thaliana, where the hybrid does not develop past the cotyledon stage and dies 3 weeks after sowing. Massive transcriptional changes take place in the hybrid, including the upregulation of most NLR (nucleotide-binding site leucine-rich repeat) disease-resistance genes. This is due to an incompatible interaction between the singleton TIR-NLR gene DANGEROUS MIX 10 (DM10), which was recently relocated from a larger NLR cluster, and an unlinked locus, DANGEROUS MIX 11 (DM11). There are multiple DM10 allelic variants in the global A. thaliana population, several of which have premature stop codons. One of these, which has a truncated LRR-PL (leucine-rich repeat [LRR]-post-LRR) region, corresponds to the DM10 risk allele. The DM10 locus and the adjacent genomic region in the risk allele carriers are highly differentiated from those in the nonrisk carriers in the global A. thaliana population, suggesting that this allele became geographically widespread only relatively recently. The DM11 risk allele is much rarer and found only in two accessions from southwestern Spain-a region from which the DM10 risk haplotype is absent-indicating that the ranges of DM10 and DM11 risk alleles may be nonoverlapping.

Shaping the leaf microbiota: plant-microbe-microbe interactions

The aerial portion of a plant, namely the leaf, is inhabited by pathogenic and non-pathogenic microbes. The leaf's physical and chemical properties, combined with fluctuating and often challenging environmental factors, create surfaces that require a high degree of adaptation for microbial colonization. As a consequence, specific interactive processes have evolved to establish a plant leaf niche. Little is known about the impact of the host immune system on phyllosphere colonization by non-pathogenic microbes. These organisms can trigger plant basal defenses and benefit the host by priming for enhanced resistance to pathogens. In most disease resistance responses, microbial signals are recognized by extra- or intracellular receptors. The interactions tend to be species specific and it is unclear how they shape leaf microbial communities. In natural habitats, microbe-microbe interactions are also important for shaping leaf communities. To protect resources, plant colonizers have developed direct antagonistic or host manipulation strategies to fight competitors. Phyllosphere-colonizing microbes respond to abiotic and biotic fluctuations and are therefore an important resource for adaptive and protective traits. Understanding the complex regulatory host-microbe-microbe networks is needed to transfer current knowledge to biotechnological applications such as plant-protective probiotics.

A Meta-Analysis Reveals Opposite Effects of Biotic and Abiotic Stresses on Transcript Levels of Arabidopsis Intracellular Immune Receptor Genes

Plant intracellular immune receptor NLR (nucleotide-binding leucine-rich repeat) proteins sense the presence of pathogens and trigger strong and robust immune responses. NLR genes are known to be tightly controlled at the protein level, but little is known about their dynamics at the transcript level. In this study, we presented a meta-analysis of transcript dynamics of all 207 NLR genes in the Col-0 accession of Arabidopsis thaliana under various biotic and abiotic stresses based on 88 publicly available RNA sequencing datasets from 27 independent studies. We find that about two thirds of the NLR genes are generally induced by pathogens, immune elicitors, or salicylic acid (SA), suggesting that transcriptional induction of NLR genes might be an important mechanism in plant immunity regulation. By contrast, NLR genes induced by biotic stresses are often repressed by abscisic acid, high temperature and drought, suggesting that transcriptional regulation of NLR genes might be important for interaction between abiotic and biotic stress responses. In addition, pathogen-induced expression of some NLR genes are dependent on SA induction. Interestingly, a small group of NLR genes are repressed under certain biotic stress treatments, suggesting an unconventional function of this group of NLRs. This meta-analysis thus reveals the transcript dynamics of NLR genes under biotic and abiotic stress conditions and suggests a contribution of NLR transcript regulation to plant immunity as well as interactions between abiotic and biotic stress responses.

Haplotype-based genome-wide association increases the predictability of leaf rust (Puccinia triticina) resistance in wheat

Resistance breeding is crucial for sustainable control of wheat leaf rust and single nucleotide polymorphism (SNP)-based genome-wide association studies (GWAS) are widely used to dissect leaf rust resistance. Unfortunately, GWAS based on SNPs often explained only a small proportion of the genetic variation. We compared SNP-based GWAS with a method based on functional haplotypes (FH) considering epistasis in a comprehensive hybrid wheat mapping population composed of 133 parents plus their 1574 hybrids and characterized with 626 245 high-quality SNPs. In total, 2408 and 1 139 828 significant associations were detected in the mapping population by using SNP-based and FH-based GWAS, respectively. These associations mapped to 25 and 69 candidate regions, correspondingly. SNP-based GWAS highlighted two already-known resistance genes, Lr22a and Lr34-B, while FH-based GWAS detected associations not only on these genes but also on two additional genes, Lr10 and Lr1. As revealed by a second hybrid wheat population for independent validation, the use of detected associations from SNP-based and FH-based GWAS reached predictabilities of 11.72% and 22.86%, respectively. Therefore, FH-based GWAS is not only more powerful for detecting associations, but also improves the accuracy of marker-assisted selection compared with the SNP-based approach.

Molecular mechanism of nanochitin whisker elicits plant resistance against Phytophthora and the receptors in plants

Rod-like nanochitin (NC) whisker with cationic nature has a strong synergistic effect with fungicides on inhibition of tobacco root rot disease. This study we explored the activity of NC against Phytophthora and the mechanism for eliciting plant defense response and the receptors in planta. P. capsici isolates, model Nicotiana benthamiana plants and Arabidopsis thaliana were treated with 0.005% of NC suspension and 1 μM of flg22. Infection control efficacy against P. capsici isolates, biosynthetic enzyme activities and the PR genes expression were determined at different hours post treatment in plant. The infection control efficacy, ROS generation, and PTI maker gene expression were re-analyzed in A. thaliana Col-0, bak1 and cerk1 mutants. The results showed that NC did not exhibit inhibitory effect on vegetative growth of P. capsici, but enhanced the resistance against P. capsici by systemically enhanced phenylalanine ammonia-lyase activity and PR gene expression. P. capsici resistance, PTI maker gene promotion, and ROS production in A. thaliana induced by NC depended not only on chitin receptor CERK1, but also BAK1. NC and flg22 induced oomycete immunity through a mechanism of a cross-microbe protection via the BAK1-CERK1 pathway in plant, pointing to the complexity of the plant immunity system.

NOD-like receptor-mediated plant immunity: from structure to cell death

Animal and plant immune systems use intracellular nucleotide-binding oligomerization domain (NOD)-like receptors (NLRs) to detect pathogens, resulting in the activation of immune responses that are often associated with localized host cell death. Whereas vertebrate NLRs detect evolutionarily conserved molecular patterns and have undergone comparatively little copy number expansion, plant NLRs detect virulence factors that have often diversified in plant pathogen populations, and thus plant NLRs have been subject to parallel diversification. Plant NLRs sense the presence of virulence factors with enzymatic virulence activity often indirectly through their modification of host target proteins. By contrast, phytopathogenic virulence factors without enzymatic activity are usually recognized by NLRs directly by their structure. Structural and biochemical analyses have shown that both indirect and direct recognition of plant pathogens trigger the oligomerization of plant NLRs into active complexes. Assembly into three-layered ring-like structures has emerged as a common principle of NLR activation in plants and animals, but with distinct amino-terminal domains initiating different signalling pathways. Collectively, these analyses point to host cell membranes as a convergence point for activated plant NLRs and the disruption of cellular ion homeostasis as a possible major factor in NLR-triggered cell death signalling.

Frontiers in Dissecting and Managing Brassica Diseases: From Reference-Based RGA Candidate Identification to Building Pan-RGAomes

The Brassica genus contains abundant economically important vegetable and oilseed crops, which are under threat of diseases caused by fungal, bacterial and viral pathogens. Resistance gene analogues (RGAs) are associated with quantitative and qualitative disease resistance and the identification of candidate RGAs associated with disease resistance is crucial for understanding the mechanism and management of diseases through breeding. The availability of Brassica genome assemblies has greatly facilitated reference-based quantitative trait loci (QTL) mapping for disease resistance. In addition, pangenomes, which characterise both core and variable genes, have been constructed for B. rapa, B. oleracea and B. napus. Genome-wide characterisation of RGAs using conserved domains and motifs in reference genomes and pangenomes reveals their clustered arrangements and presence of structural variations. Here, we comprehensively review RGA identification in important Brassica genome and pangenome assemblies. Comparison of the RGAs in QTL between resistant and susceptible individuals allows for efficient identification of candidate disease resistance genes. However, the reference-based QTL mapping and RGA candidate identification approach is restricted by the under-represented RGA diversity characterised in the limited number of Brassica assemblies. The species-wide repertoire of RGAs make up the pan-resistance gene analogue genome (pan-RGAome). Building a pan-RGAome, through either whole genome resequencing or resistance gene enrichment sequencing, would effectively capture RGA diversity, greatly expanding breeding resources that can be utilised for crop improvement.

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.

What is the Molecular Basis of Nonhost Resistance?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.Nonhost resistance is typically considered the ability of a plant species to repel all attempts of a pathogen species to colonize it and reproduce on it. Based on this common definition, nonhost resistance is presumed to be very durable and, thus, of great interest for its potential use in agriculture. Despite considerable research efforts, the molecular basis of this type of plant immunity remains nebulous. We here stress the fact that "nonhost resistance" is a phenomenological rather than a mechanistic concept that comprises more facets than typically considered. We further argue that nonhost resistance essentially relies on the very same genes and pathways as other types of plant immunity, of which some may act as bottlenecks for particular pathogens on a given plant species or under certain conditions. Thus, in our view, the frequently used term "nonhost genes" is misleading and should be avoided. Depending on the plant-pathogen combination, nonhost resistance may involve the recognition of pathogen effectors by host immune sensor proteins, which might give rise to host shifts or host range expansions due to evolutionary-conditioned gains and losses in respective armories. Thus, the extent of nonhost resistance also defines pathogen host ranges. In some instances, immune-related genes can be transferred across plant species to boost defense, resulting in augmented disease resistance. We discuss future routes for deepening our understanding of nonhost resistance and argue that the confusing term "nonhost resistance" should be used more cautiously in the light of a holistic view of plant immunity.

Whole Genome Re-sequencing of Soybean Accession EC241780 Providing Genomic Landscape of Candidate Genes Involved in Rust Resistance

BACKGROUND

In this study, whole genome re-sequencing of rust resistant soybean genotype EC241780 was performed to understand the genomic landscape involved in the resistance mechanism.

METHODS

A total of 374 million raw reads were obtained with paired-end sequencing performed with Illumina HiSeq 2500 instrument, out of which 287.3 million high quality reads were mapped to Williams 82 reference genome. Comparative sequence analysis of EC241780 with rust susceptible cultivars Williams 82 and JS 335 was performed to identify sequence variation and to prioritise the candidate genes.

RESULTS

Comparative analysis indicates that genotype EC241780 has high sequence similarity with rust resistant genotype PI 200492 and the resistance in EC241780 is conferred by the Rpp1 locus. Based on the sequence variations and functional annotations, three genes Glyma18G51715, Glyma18G51741 and Glyma18G51765 encoding for NBS-LRR family protein were identified as the most prominent candidate for Rpp1 locus.

CONCLUSION

The study provides insights of genome-wide sequence variation more particularly at Rpp1 loci which will help to develop rust resistant soybean cultivars through efficient exploration of the genomic resource.

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.

A molecular roadmap to the plant immune system

Plant diseases caused by pathogens and pests are a constant threat to global food security. Direct crop losses and the measures used to control disease (e.g. application of pesticides) have significant agricultural, economic, and societal impacts. Therefore, it is essential that we understand the molecular mechanisms of the plant immune system, a system that allows plants to resist attack from a wide variety of organisms ranging from viruses to insects. Here, we provide a roadmap to plant immunity, with a focus on cell-surface and intracellular immune receptors. We describe how these receptors perceive signatures of pathogens and pests and initiate immune pathways. We merge existing concepts with new insights gained from recent breakthroughs on the structure and function of plant immune receptors, which have generated a shift in our understanding of cell-surface and intracellular immunity and the interplay between the two. Finally, we use our current understanding of plant immunity as context to discuss the potential of engineering the plant immune system with the aim of bolstering plant defenses against disease.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

Molecular Basis of Disease Resistance and Perspectives on Breeding Strategies for Resistance Improvement in Crops

Crop diseases are major factors responsible for substantial yield losses worldwide, which affects global food security. The use of resistance (R) genes is an effective and sustainable approach to controlling crop diseases. Here, we review recent advances on R gene studies in the major crops and related wild species. Current understanding of the molecular mechanisms underlying R gene activation and signaling, and susceptibility (S) gene-mediated resistance in crops are summarized and discussed. Furthermore, we propose some new strategies for R gene discovery, how to balance resistance and yield, and how to generate crops with broad-spectrum disease resistance. With the rapid development of new genome-editing technologies and the availability of increasing crop genome sequences, the goal of breeding next-generation crops with durable resistance to pathogens is achievable, and will be a key step toward increasing crop production in a sustainable way.

CRISPR-Cas9 enrichment and long read sequencing for fine mapping in plants

BACKGROUND

Genomic methods for identifying causative variants for trait loci applicable to a wide range of germplasm are required for plant biologists and breeders to understand the genetic control of trait variation.

RESULTS

We implemented Cas9-targeted sequencing for fine-mapping in apple, a method combining CRISPR-Cas9 targeted cleavage of a region of interest, followed by enrichment and long-read sequencing using the Oxford Nanopore Technology (ONT). We demonstrated the capability of this methodology to specifically cleave and enrich a plant genomic locus spanning 8 kb. The repeated mini-satellite motif located upstream of the Malus × domestica (apple) MYB10 transcription factor gene, causing red fruit colouration when present in a heterozygous state, was our exemplar to demonstrate the efficiency of this method: it contains a genomic region with a long structural variant normally ignored by short-read sequencing technologiesCleavage specificity of the guide RNAs was demonstrated using polymerase chain reaction products, before using them to specify cleavage of high molecular weight apple DNA. An enriched library was subsequently prepared and sequenced using an ONT MinION flow cell (R.9.4.1). Of the 7,056 ONT reads base-called using both Albacore2 (v2.3.4) and Guppy (v3.2.4), with a median length of 9.78 and 9.89 kb, respectively, 85.35 and 91.38%, aligned to the reference apple genome. Of the aligned reads, 2.98 and 3.04% were on-target with read depths of 180 × and 196 × for Albacore2 and Guppy, respectively, and only five genomic loci were off-target with read depth greater than 25 × , which demonstrated the efficiency of the enrichment method and specificity of the CRISPR-Cas9 cleavage.

CONCLUSIONS

We demonstrated that this method can isolate and resolve single-nucleotide and structural variants at the haplotype level in plant genomic regions. The combination of CRISPR-Cas9 target enrichment and ONT sequencing provides a more efficient technology for fine-mapping loci than genome-walking approaches.

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.

Origins and Immunity Networking Functions of EDS1 Family Proteins

The EDS1 family of structurally unique lipase-like proteins EDS1, SAG101, and PAD4 evolved in seed plants, on top of existing phytohormone and nucleotide-binding-leucine-rich-repeat (NLR) networks, to regulate immunity pathways against host-adapted biotrophic pathogens. Exclusive heterodimers between EDS1 and SAG101 or PAD4 create essential surfaces for resistance signaling. Phylogenomic information, together with functional studies in Arabidopsis and tobacco, identify a coevolved module between the EDS1-SAG101 heterodimer and coiled-coil (CC) HET-S and LOP-B (CC_HELO) domain helper NLRs that is recruited by intracellular Toll-interleukin1-receptor (TIR) domain NLR receptors to confer host cell death and pathogen immunity. EDS1-PAD4 heterodimers have a different and broader activity in basal immunity that transcriptionally reinforces local and systemic defenses triggered by various NLRs. Here, we consider EDS1 family protein functions across seed plant lineages in the context of networking with receptor and helper NLRs and downstream resistance machineries. The different modes of action and pathway connectivities of EDS1 family members go some way to explaining their central role in biotic stress resilience.

What natural variation can teach us about resistance durability

Breeding a crop variety to be resistant to a pathogen usually takes years. This is problematic because pathogens, with short generation times and fluid genomes, adapt quickly to overcome resistance. The triumph of the pathogen is not inevitable, however, as there are numerous examples of durable resistance, particularly in wild plants. Which factors then contribute to such resistance stability over millennia? We review current knowledge of wild and agricultural pathosystems, detailing the importance of genetic, species and spatial heterogeneity in the prevention of pathogen outbreaks. We also highlight challenges associated with increasing resistance diversity in crops, both in light of pathogen (co-)evolution and breeding practices. Historically it has been difficult to incorporate heterogeneity into agriculture due to reduced efficiency in harvesting. Recent advances implementing computer vision and automation in agricultural production may improve our ability to harvest mixed genotype and mixed species plantings, thereby increasing resistance durability.

Plant NLRs get by with a little help from their friends

Many plant NLR (nucleotide-binding, leucine-rich repeat) immune receptors require other NLRs for their function. In pairs of chromosomally adjacent sensor/helper NLRs, the sensor typically carries an integrated domain (ID) that mimics the authentic target of a pathogen effector. The RPW8-NLR clade supports the function of many diverse plant NLRs, particularly those with a TIR N-terminal domain, in concert with a family of EP-domain containing signalling partners. The NRC clade of NLRs are required for the function of many unlinked sensor NLRs in Solanaceous plants. We evaluate recent advances in paired NLR biology in the context of the structure and possible mechanisms of the first defined plant inflammasome containing ZAR1.

Exploring the application of wild species for crop improvement in a changing climate

Modern agriculture is currently facing challenges from a burgeoning population and changing climate, which requires improved crops with adaptation to climate and elite yield and quality traits. While there is a breeding bottleneck caused by intensive selection, gene banks containing conserved wild relatives and landraces can be used as breeding resources. However, with limited genetic information available on these wild relatives, the application has been hindered. With the development of both genomics and bioinformatics techniques, it is now easier to identify the genetic variation in wild species, which can be utilized for the introgression of elite traits. These wild species can therefore play an important role in food security and breeding sustainability.

Plant pan-genomes are the new reference

Recent years have seen a surge in plant genome sequencing projects and the comparison of multiple related individuals. The high degree of genomic variation observed led to the realization that single reference genomes do not represent the diversity within a species, and led to the expansion of the pan-genome concept. Pan-genomes represent the genomic diversity of a species and includes core genes, found in all individuals, as well as variable genes, which are absent in some individuals. Variable gene annotations often show similarities across plant species, with genes for biotic and abiotic stress commonly enriched within variable gene groups. Here we review the growth of pan-genomics in plants, explore the origins of gene presence and absence variation, and show how pan-genomes can support plant breeding and evolution studies.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

The evolving landscape of cell surface pattern recognition across plant immune networks

To recognize diverse threats, plants monitor extracellular molecular patterns and transduce intracellular immune signaling through receptor complexes at the plasma membrane. Pattern recognition occurs through a prototypical network of interacting proteins, comprising A) receptors that recognize inputs associated with a growing number of pest and pathogen classes (bacteria, fungi, oomycetes, caterpillars), B) co-receptor kinases that participate in binding and signaling, and C) cytoplasmic kinases that mediate first stages of immune output. While this framework has been elucidated in reference accessions of model organisms, network components are part of gene families with widespread variation, potentially tuning immunocompetence for specific contexts. Most dramatically, variation in receptor repertoires determines the range of ligands acting as immunogenic inputs for a given plant. Diversification of receptor kinase (RK) and related receptor-like protein (RLP) repertoires may tune responses even within a species. Comparative genomics at pangenome scale will reveal patterns and features of immune network variation.

Variation Patterns of NLR Clusters in Arabidopsis thaliana Genomes

The nucleotide-binding domain and leucine-rich repeat (NLR) gene family is highly expanded in the plant lineage with extensive sequence and structure polymorphisms. To survey the landscape of NLR expansion, we mined the published long-read data generated by the resistance gene enrichment sequencing of 64 diverse Arabidopsis thaliana accessions. We found that the hot spots of massive multi-gene NLR cluster expansion did not typically span the whole cluster; instead, they were restricted to a handful of, or only one, dominant radiation(s). All sequences in such a radiation were distinct from other genes in the cluster but not from each other in the clade, making it difficult to assign trustworthy reference-based orthologies when multiple reference genes were present in the radiation. Consequently, NLR genes can be broadly divided into two types: radiating or high-fidelity, where high-fidelity genes are well conserved and well separated from other clades. A similar distinction could be made for NLR clusters, depending on whether cluster size was determined primarily by extensive radiation or the presence of numerous high-fidelity genes. We also identified groups of well-conserved NLR clades that were missing from the Columbia-0 reference genome. This suggests that the classification of NLRs using gene IDs from a single reference accession can rarely capture all major paralogs in a cluster accurately and representatively and that a reference-agnostic perspective is required to properly characterize these additional variations. Finally, we present a quantitative visualization method for differentiating these situations in a given clade of interest.

Current status of the multinational Arabidopsis community

The multinational Arabidopsis research community is highly collaborative and over the past thirty years these activities have been documented by the Multinational Arabidopsis Steering Committee (MASC). Here, we (a) highlight recent research advances made with the reference plant Arabidopsis thaliana; (b) provide summaries from recent reports submitted by MASC subcommittees, projects and resources associated with MASC and from MASC country representatives; and (c) initiate a call for ideas and foci for the "fourth decadal roadmap," which will advise and coordinate the global activities of the Arabidopsis research community.

High-resolution expression profiling of selected gene sets during plant immune activation

The plant immune system involves detection of pathogens via both cell-surface and intracellular receptors. Both receptor classes can induce transcriptional reprogramming that elevates disease resistance. To assess differential gene expression during plant immunity, we developed and deployed quantitative sequence capture (CAP-I). We designed and synthesized biotinylated single-strand RNA bait libraries targeted to a subset of defense genes, and generated sequence capture data from 99 RNA-seq libraries. We built a data processing pipeline to quantify the RNA-CAP-I-seq data, and visualize differential gene expression. Sequence capture in combination with quantitative RNA-seq enabled cost-effective assessment of the expression profile of a specified subset of genes. Quantitative sequence capture is not limited to RNA-seq or any specific organism and can potentially be incorporated into automated platforms for high-throughput sequencing.

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

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

Evolution of Plant NLRs: From Natural History to Precise Modifications

Nucleotide-binding leucine-rich repeat receptors (NLRs) monitor the plant intracellular environment for signs of pathogen infection. Several mechanisms of NLR-mediated immunity arose independently across multiple species. These include the functional specialization of NLRs into sensors and helpers, the independent emergence of direct and indirect recognition within NLR subfamilies, the regulation of NLRs by small RNAs, and the formation of NLR networks. Understanding the evolutionary history of NLRs can shed light on both the origin of pathogen recognition and the common constraints on the plant immune system. Attempts to engineer disease resistance have been sparse and rarely informed by evolutionary knowledge. In this review, we discuss the evolution of NLRs, give an overview of previous engineering attempts, and propose how to use evolutionary knowledge to advance future research in the generation of novel disease-recognition capabilities.

Pathogen-Associated Molecular Pattern-Triggered Immunity Involves Proteolytic Degradation of Core Nonsense-Mediated mRNA Decay Factors During the Early Defense Response

Nonsense-mediated mRNA decay (NMD), an mRNA quality control process, is thought to function in plant immunity. A subset of fully spliced (FS) transcripts of Arabidopsis (Arabidopsis thaliana) resistance (R) genes are upregulated during bacterial infection. Here, we report that 81.2% and 65.1% of FS natural TIR-NBS-LRR (TNL) and CC-NBS-LRR transcripts, respectively, retain characteristics of NMD regulation, as their transcript levels could be controlled posttranscriptionally. Both bacterial infection and the perception of bacteria by pattern recognition receptors initiated the destruction of core NMD factors UP-FRAMESHIFT1 (UPF1), UPF2, and UPF3 in Arabidopsis within 30 min of inoculation via the independent ubiquitination of UPF1 and UPF3 and their degradation via the 26S proteasome pathway. The induction of UPF1 and UPF3 ubiquitination was delayed in mitogen-activated protein kinase3 (mpk3) and mpk6, but not in salicylic acid-signaling mutants, during the early immune response. Finally, previously uncharacterized TNL-type R transcripts accumulated in upf mutants and conferred disease resistance to infection with a virulent Pseudomonas strain in plants. Our findings demonstrate that NMD is one of the main regulatory processes through which PRRs fine-tune R transcript levels to reduce fitness costs and achieve effective immunity.

Atypical Resistance Protein RPW8/HR Triggers Oligomerization of the NLR Immune Receptor RPP7 and Autoimmunity

In certain plant hybrids, immunity signaling is initiated when immune components interact in the absence of a pathogen trigger. In Arabidopsis thaliana, such autoimmunity and cell death are linked to variants of the NLR RPP7 and the RPW8 proteins involved in broad-spectrum resistance. We uncover the molecular basis for this autoimmunity and demonstrate that a homolog of RPW8, HR4Fei-0, can trigger the assembly of a higher-order RPP7 complex, with autoimmunity signaling as a consequence. HR4Fei-0-mediated RPP7 oligomerization occurs via the RPP7 C-terminal leucine-rich repeat (LRR) domain and ATP-binding P-loop. RPP7 forms a higher-order complex only in the presence of HR4Fei-0 and not with the standard HR4 variant, which is distinguished from HR4Fei-0 by length variation in C-terminal repeats. Additionally, HR4Fei-0 can independently form self-oligomers, which directly kill cells in an RPP7-independent manner. Our work provides evidence for a plant resistosome complex and the mechanisms by which RPW8/HR proteins trigger cell death.

Chromosome-level assemblies of multiple Arabidopsis genomes reveal hotspots of rearrangements with altered evolutionary dynamics

Despite hundreds of sequenced Arabidopsis genomes, very little is known about the degree of genomic collinearity within single species, due to the low number of chromosome-level assemblies. Here, we report chromosome-level reference-quality assemblies of seven Arabidopsis thaliana accessions selected across its global range. Each genome reveals between 13–17 Mb rearranged, and 5–6 Mb non-reference sequences introducing copy-number changes in ~5000 genes, including ~1900 non-reference genes. Quantifying the collinearity between the genomes reveals ~350 euchromatic regions, where accession-specific tandem duplications destroy the collinearity between the genomes. These hotspots of rearrangements are characterized by reduced meiotic recombination in hybrids and genes implicated in biotic stress response. This suggests that hotspots of rearrangements undergo altered evolutionary dynamics, as compared to the rest of the genome, which are mostly based on the accumulation of new mutations and not on the recombination of existing variation, and thereby enable a quick response to the biotic stress.

Despite tremendous genomic resources in the Arabidopsis community, only a few whole genome de novo assemblies are available. Here, the authors report chromosome-level reference-quality assemblies of seven A. thaliana accessions and reveal hotspots of rearrangements with altered evolutionary dynamics.

Plant NLRs: The Whistleblowers of Plant Immunity

The study of plant diseases is almost as old as agriculture itself. Advancements in molecular biology have given us much more insight into the plant immune system and how it detects the many pathogens plants may encounter. Members of the primary family of plant resistance (R) proteins, NLRs, contain three distinct domains, and appear to use several different mechanisms to recognize pathogen effectors and trigger immunity. Understanding the molecular process of NLR recognition and activation has been greatly aided by advancements in structural studies, with ZAR1 recently becoming the first full-length NLR to be visualized. Genetic and biochemical analysis identified many critical components for NLR activation and homeostasis control. The increased study of helper NLRs has also provided insights into the downstream signaling pathways of NLRs. This review summarizes the progress in the last decades on plant NLR research, focusing on the mechanistic understanding that has been achieved.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.