Intragenic recombination and diversifying selection contribute to the evolution of downy mildew resistance at the RPP8 locus of Arabidopsis

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

A single amino acid at position 31 in the N-terminus of the coat protein of cucumber mosaic virus determines its avirulence function for RCY1-conferred virus resistance

The coat protein (CP) of the cucumber mosaic virus (CMV) yellow strain [CMV(Y)], but not the CMV B2 strain [CMV(B2)], serves as an avirulence determinant against the NB-LRR class RCY1 of Arabidopsis thaliana. To investigate the avirulence function, a series of binary vectors were constructed by partially exchanging the CP coding sequence between CMV(Y) and CMV(B2) or introducing nucleotide substitutions. These vectors were transiently expressed in Nicotiana benthamiana leaves transformed with modified RCY1 cDNA. Analysis of hypersensitive resistance-cell death (HCD), CP accumulation, and defense gene expression at leaf sites infiltrated with Agrobacterium indicated that a single amino acid at position 31 of the CP seems to determine the avirulence function.

Barley MLA3 recognizes the host-specificity effector Pwl2 from Magnaporthe oryzae

Plant nucleotide-binding leucine-rich repeat (NLRs) immune receptors directly or indirectly recognize pathogen-secreted effector molecules to initiate plant defense. Recognition of multiple pathogens by a single NLR is rare and usually occurs via monitoring for changes to host proteins; few characterized NLRs have been shown to recognize multiple effectors. The barley (Hordeum vulgare) NLR gene Mildew locus a (Mla) has undergone functional diversification, and the proteins encoded by different Mla alleles recognize host-adapted isolates of barley powdery mildew (Blumeria graminis f. sp. hordei [Bgh]). Here, we show that Mla3 also confers resistance to the rice blast fungus Magnaporthe oryzae in a dosage-dependent manner. Using a forward genetic screen, we discovered that the recognized effector from M. oryzae is Pathogenicity toward Weeping Lovegrass 2 (Pwl2), a host range determinant factor that prevents M. oryzae from infecting weeping lovegrass (Eragrostis curvula). Mla3 has therefore convergently evolved the capacity to recognize effectors from diverse pathogens.

Integrated Molecular and Bioinformatics Approaches for Disease-Related Genes in Plants

Modern plant pathology relies on bioinformatics approaches to create novel plant disease diagnostic tools. In recent years, a significant amount of biological data has been generated due to rapid developments in genomics and molecular biology techniques. The progress in the sequencing of agriculturally important crops has made it possible to develop a better understanding of plant-pathogen interactions and plant resistance. The availability of host-pathogen genome data offers effective assistance in retrieving, annotating, analyzing, and identifying the functional aspects for characterization at the gene and genome levels. Physical mapping facilitates the identification and isolation of several candidate resistance (R) genes from diverse plant species. A large number of genetic variations, such as disease-causing mutations in the genome, have been identified and characterized using bioinformatics tools, and these desirable mutations were exploited to develop disease resistance. Moreover, crop genome editing tools, namely the CRISPR (clustered regulatory interspaced short palindromic repeats)/Cas9 (CRISPR-associated) system, offer novel and efficient strategies for developing durable resistance. This review paper describes some aspects concerning the databases, tools, and techniques used to characterize resistance (R) genes for plant disease management.

Fungal effectors versus defense-related genes of B. juncea and the status of resistant transgenics against fungal pathogens

Oilseed brassica has become instrumental in securing global food and nutritional security. B. juncea, colloquially known as Indian mustard, is cultivated across tropics and subtropics including Indian subcontinent. The production of Indian mustard is severely hampered by fungal pathogens which necessitates human interventions. Chemicals are often resorted to as they are quick and effective, but due to their economic and ecological unsustainability, there is a need to explore their alternatives. The B. juncea-fungal pathosystem is quite diverse as it covers broad-host range necrotrophs (Sclerotinia sclerotiorum), narrow-host range necrotrophs (Alternaria brassicae and A. brassicicola) and biotrophic oomycetes (Albugo candida and Hyaloperonospora brassica). Plants ward off fungal pathogens through two-step resistance mechanism; PTI which involves recognition of elicitors and ETI where the resistance gene (R gene) interacts with the fungal effectors. The hormonal signalling is also found to play a vital role in defense as the JA/ET pathway is initiated at the time of necrotroph infection and SA pathway is induced when the biotrophs attack plants. The review discuss the prevalence of fungal pathogens of Indian mustard and the studies conducted on effectoromics. It covers both pathogenicity conferring genes and host-specific toxins (HSTs) that can be used for a variety of purposes such as identifying cognate R genes, understanding pathogenicity and virulence mechanisms, and establishing the phylogeny of fungal pathogens. It further encompasses the studies on identifying resistant sources and characterisation of R genes/quantitative trait loci and defense-related genes identified in Brassicaceae and unrelated species which, upon introgression or overexpression, confer resistance. Finally, the studies conducted on developing resistant transgenics in Brassicaceae have been covered in which chitinase and glucanase genes are mostly used. The knowledge gained from this review can further be used for imparting resistance against major fungal pathogens.

Comparative transcriptome profiling of near isogenic lines PBW343 and FLW29 to unravel defense related genes and pathways contributing to stripe rust resistance in wheat

Stripe rust (Sr), caused by Puccinia striiformis f. sp. tritici (Pst), is the most devastating disease that poses serious threat to the wheat-growing nations across the globe. Developing resistant cultivars is the most challenging aspect in wheat breeding. The function of resistance genes (R genes) and the mechanisms by which they influence plant-host interactions are poorly understood. In the present investigation, comparative transcriptome analysis was carried out by involving two near-isogenic lines (NILs) PBW343 and FLW29. The seedlings of both the genotypes were inoculated with Pst pathotype 46S119. In total, 1106 differentially expressed genes (DEGs) were identified at early stage of infection (12 hpi), whereas expressions of 877 and 1737 DEGs were observed at later stages (48 and 72 hpi) in FLW29. The identified DEGs were comprised of defense-related genes including putative R genes, 7 WRKY transcriptional factors, calcium, and hormonal signaling associated genes. Moreover, pathways involved in signaling of receptor kinases, G protein, and light showed higher expression in resistant cultivar and were common across different time points. Quantitative real-time PCR was used to further confirm the transcriptional expression of eight critical genes involved in plant defense mechanism against stripe rust. The information about genes are likely to improve our knowledge of the genetic mechanism that controls the stripe rust resistance in wheat, and data on resistance response-linked genes and pathways will be a significant resource for future research.

QTL mapping of resistance to Pseudoperonospora cubensis clade 2, mating type A1, in Cucumis melo and dual-clade marker development

This is the first identification of QTLs underlying resistance in Cucumis melo to an isolate of Pseudoperonospora cubensis identified as Clade 2/mating type A1. Pseudoperonospora cubensis, causal organism of cucurbit downy mildew (CDM), causes severe necrosis and defoliation on Cucumis melo (melon). A recombinant inbred line population (N = 169) was screened against an isolate of P. cubensis (Clade 2/mating type A1) in replicated greenhouse and growth chamber experiments. SNPs (n = 5633 bins) identified in the RIL population were used for quantitative trait loci (QTL) mapping. A single major QTL on chromosome 10 (qPcub-10.3-10.4) was consistently associated with resistance across all experiments, while a second major QTL on chromosome 8 (qPcub-8.3) was identified only in greenhouse experiments. These two major QTLs were identified on the same chromosomes (8 and 10) but in different locations as two major QTLs (qPcub-8.2 and qPcub-10.1) previously identified for resistance to P. cubensis Clade 1/mating type A2. Kompetitive allele-specific PCR (KASP) markers were developed for these four major QTLs and validated in the RIL population through QTL mapping. These markers will provide melon breeders a high-throughput genotyping toolkit for development of melon cultivars with broad tolerance to CDM.

Genome-wide identification of KNOX transcription factors in cotton and the role of GhKNOX4-A and GhKNOX22-D in response to salt and drought stress

Cotton is one of the most important economic and fiber crops in the world. KNOX is one class of universal transcription factors, which plays important roles in plant growth and development as well as response to different stresses. Although there are many researches on KNOXs in other plant species, there are few reports on cotton. In this study, we systematically and comprehensively identified all KNOX genes in upland cotton and its two ancestral species; we also studied their functions by employing RNA-seq analysis and virus-induced gene silence (VIGS). A total of 89 KNOX genes were identified from three cotton species. Among them, 44 were from upland cotton, 22 and 23 were found in its ancestral species G. raimondii and G. arboreum, respectively. Plant polyploidization and domestication play a selective force driving KNOX gene evolution. Phylogenetic analysis displayed that KNOX genes were evolved into three Classes. The intron length and exon number differed in each Class. Transcriptome data showed that KNOX genes of Class II were widely expressed in multiple tissues, including fiber. The majority of KNOX genes were induced by different abiotic stresses. Additionally, we found multiple cis-elements related to stress in the promoter region of KNOX genes. VIGS silence of GhKNOX4-A and GhKNOX22-D genes showed significant growth and development effect in cotton seedlings under salt and drought treatments. Both GhKNOX4-A and GhKNOX22-D regulated plant tolerance; silencing both genes induced oxidative stresses, evidenced by reduced SOD activity and induced leave cell death, and also enhanced stomatal open and water loss. Thus, GhKNOX4-A and GhKNOX22-D may contribute to drought response by regulating stomata opening and oxidative stresses.

The chromosome-scale genome of Phoebe bournei reveals contrasting fates of terpene synthase (TPS)-a and TPS-b subfamilies

Terpenoids, including aromatic volatile monoterpenoids and sesquiterpenoids, function in defense against pathogens and herbivores. Phoebe trees are remarkable for their scented wood and decay resistance. Unlike other Lauraceae species investigated to date, Phoebe species predominantly accumulate sesquiterpenoids instead of monoterpenoids. Limited genomic data restrict the elucidation of terpenoid variation and functions. Here, we present a chromosome-scale genome assembly of a Lauraceae tree, Phoebe bournei, and identify 72 full-length terpene synthase (TPS) genes. Genome-level comparison shows pervasive lineage-specific duplication and contraction of TPS subfamilies, which have contributed to the extreme terpenoid variation within Lauraceae species. Although the TPS-a and TPS-b subfamilies were both expanded via tandem duplication in P. bournei, more TPS-a copies were retained and constitutively expressed, whereas more TPS-b copies were lost. The TPS-a genes on chromosome 8 functionally diverged to synthesize eight highly accumulated sesquiterpenes in P. bournei. The essential oil of P. bournei and its main component, β-caryophyllene, exhibited antifungal activities against the three most widespread canker pathogens of trees. The TPS-a and TPS-b subfamilies have experienced contrasting fates over the evolution of P. bournei. The abundant sesquiterpenoids produced by TPS-a proteins contribute to the excellent pathogen resistance of P. bournei trees. Overall, this study sheds light on the evolution and adaptation of terpenoids in Lauraceae and provides valuable resources for boosting plant immunity against pathogens in various trees and crops.

Predicting Cloned Disease Resistance Gene Homologs (CDRHs) in Radish, Underutilised Oilseeds, and Wild Brassicaceae Species

Brassicaceae crops, including Brassica, Camelina and Raphanus species, are among the most economically important crops globally; however, their production is affected by several diseases. To predict cloned disease resistance (R) gene homologs (CDRHs), we used the protein sequences of 49 cloned R genes against fungal and bacterial diseases in Brassicaceae species. In this study, using 20 Brassicaceae genomes (17 wild and 3 domesticated species), 3172 resistance gene analogs (RGAs) (2062 nucleotide binding-site leucine-rich repeats (NLRs), 497 receptor-like protein kinases (RLKs) and 613 receptor-like proteins (RLPs)) were identified. CDRH clusters were also observed in Arabis alpina, Camelina sativa and Cardamine hirsuta with assigned chromosomes, consisting of 62 homogeneous (38 NLR, 17 RLK and 7 RLP clusters) and 10 heterogeneous RGA clusters. This study highlights the prevalence of CDRHs in the wild relatives of the Brassicaceae family, which may lay the foundation for rapid identification of functional genes and genomics-assisted breeding to develop improved disease-resistant Brassicaceae crop cultivars.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Suspension cell secretome of the grain legume Lathyrus sativus (grasspea) reveals roles in plant development and defense responses

Plant secretomics has been especially important in understanding the molecular basis of plant development, stress resistance and biomarker discovery. In addition to sharing a similar role in maintaining cell metabolism and biogenesis with the animal secretome, plant-secreted proteins actively participate in signaling events crucial for cellular homeostasis during stress adaptation. However, investigation of the plant secretome remains largely overlooked, particularly in pulse crops, demanding urgent attention. To better understand the complexity of the secretome, we developed a reference map of a stress-resilient orphan legume, Lathyrus sativus (grasspea), which can be utilized as a potential proteomic resource. Secretome analysis of L. sativus led to the identification of 741 nonredundant proteins belonging to a myriad of functional classes, including antimicrobial, antioxidative and redox potential. Computational prediction of the secretome revealed that ∼29% of constituents are predicted to follow unconventional protein secretion (UPS) routes. We conducted additional in planta analysis to determine the localization of two secreted proteins, recognized as cell surface residents. Sequence-based homology comparison revealed that L. sativus shares ∼40% of the constituents reported thus far from in vitro and in planta secretome analysis in model and crop species. Significantly, we identified 571 unique proteins secreted from L. sativus involved in cell-to-cell communication, organ development, kinase-mediated signaling, and stress perception, among other critical roles. Conclusively, the grasspea secretome participates in putative crosstalk between genetic circuits that regulate developmental processes and stress resilience.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host–pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

Biotechnological approaches in management of oomycetes diseases

Plant pathogenic oomycetes cause significant impact on agriculture and, therefore, their management is utmost important. Though conventional methods to combat these pathogens (resistance breeding and use of fungicides) are available but these are limited by the availability of resistant cultivars due to evolution of new pathogenic races, development of resistance in the pathogens against agrochemicals and their potential hazardous effects on the environment and human health. This has fuelled a continual search for novel and alternate strategies for management of phytopathogens. The recent advances in oomycetes genome (Phytophthora infestans, P. ramorum, P. sojae, Pythium ultimum, Albugo candida etc.) would further help in understanding host-pathogen interactions essentially needed for designing effective management strategies. In the present communication the novel and alternate strategies for the management of oomycetes diseases are discussed.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Comparative transcriptome of compatible and incompatible interaction of Erysiphe pisi and garden pea reveals putative defense and pathogenicity factors

Comparative transcriptome analysis of Erysiphe pisi-infected pea (Pisum sativum) genotypes JI-2480 (resistant) and Arkel (susceptible) at 72 hours post-inoculation (hpi) was carried to detect molecular components involved in compatible and incompatible interactions. Differential gene expression was observed in Arkel and JI-2480 genotype at 72 hpi with E. pisi isolate (Ep01) using EdgeR software. Out of 32 217 transcripts, 2755 transcripts showed significantly altered gene expression in case of plants while 530 were related to E. pisi (P < 0.05). The higher transcript number of differentially expressed genes demonstrated peak activity of pathogenicity genes in plants at 72 hpi. Glycolysis was observed to be the major pathway for energy source during fungal growth. Differential gene expression of plant transcripts revealed significant expression of putative receptor and regulatory sequences involved in defense in the resistant, JI-2480 compared to susceptible, Arkel genotype. Expression of genes involved in defense and hormonal signaling, genes related to hypersensitive response, reactive oxygen species and phenylpropanoid pathway in JI-2480 indicated their crucial role in disease resistance against E. pisi. Down-regulation of transcription factors like-WRKY-28 and up-regulation of several putative pattern recognition receptors in JI-2480 compared to Arkel also suggested activation of host-mediated defense responses against E. pisi in pea.

The rice LRR-like1 protein YELLOW AND PREMATURE DWARF 1 is involved in leaf senescence induced by high light

Senescence in plants is induced by endogenous physiological changes and exogenous stresses. In this study, we isolated two alleles of a novel rice (Oryza sativa) mutant, yellow and premature dwarf 1 (ypd1). The ypd1 mutants exhibited a yellow and dwarf phenotype from germination, and premature senescence starting at tillering. Moreover, the ypd1 mutants were sensitive to high light, which accelerated cell death and senescence. Consistent with their yellow phenotype, the ypd1 mutants had abnormal chloroplasts and lower levels of photosynthetic pigments. TUNEL assays together with histochemical staining demonstrated that ypd1 mutants showed cell death and that they accumulated reactive oxygen species. The ypd1 mutants also showed increased expression of genes associated with senescence. Map-based cloning revealed a G→A substitution in exon 6 (ypd1-1) and exon 13 (ypd1-2) of LOC_Os06g13050 that affected splicing and caused premature termination of the encoded protein. YPD1 was found to be preferentially expressed in the leaf and it encodes a LRR-like1 protein. Complementation, overexpression, and targeted deletion confirmed that the mutations in YPD1 caused the ypd1 phenotype. YPD1 was localized on the chloroplast membrane. Our results thus demonstrate that the novel rice LRR-like1 protein YPD1 affects chloroplast development and leaf senescence.

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.

Distinct Evolutionary Patterns of NBS-Encoding Genes in Three Soapberry Family (Sapindaceae) Species

Nucleotide-binding site (NBS)-type disease resistance genes (R genes) play key roles in plant immune responses and have co-evolved with pathogens over the course of plant lifecycles. Comparative genomic studies tracing the dynamic evolution of NBS-encoding genes have been conducted using many important plant lineages. However, studies on Sapindaceae species have not been performed. In this study, a discrepant number of NBS-encoding genes were identified in the genomes of Xanthoceras sorbifolium (180), Dinnocarpus longan (568), and Acer yangbiense (252). These genes were unevenly distributed and usually clustered as tandem arrays on chromosomes, with few existed as singletons. The phylogenetic analysis revealed that NBS-encoding genes formed three monophyletic clades, RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL), and CC-NBS-LRR (CNL), which were distinguished by amino acid motifs. The NBS-encoding genes of the X. sorbifolium, D. longan, and A. yangbiense genomes were derived from 181 ancestral genes (three RNL, 23 TNL, and 155 CNL), which exhibited dynamic and distinct evolutionary patterns due to independent gene duplication/loss events. Specifically, X. sorbifolium exhibited a “first expansion and then contraction” evolutionary pattern, while A. yangbiense and D. longan exhibited a “first expansion followed by contraction and further expansion” evolutionary pattern. However, further expansion in D. longan was stronger than in A. yangbiense after divergence, suggesting that D. longan gained more genes in response to various pathogens. Additionally, the ancient and recent expansion of CNL genes generated the dominance of this subclass in terms of gene numbers, while the low copy number status of RNL genes was attributed to their conserved functions.

Deciphering the Impact of a Bacterial Infection on Meiotic Recombination in Arabidopsis with Fluorescence Tagged Lines

Plants are under strong evolutionary pressure to maintain surveillance against pathogens. One major disease resistance mechanism is based on NB-LRR (NLR) proteins that specifically recognize pathogen effectors. The cluster organization of the NLR gene family could favor sequence exchange between NLR genes via recombination, favoring their evolutionary dynamics. Increasing data, based on progeny analysis, suggest the existence of a link between the perception of biotic stress and the production of genetic diversity in the offspring. This could be driven by an increased rate of meiotic recombination in infected plants, but this has never been strictly demonstrated. In order to test if pathogen infection can increase DNA recombination in pollen meiotic cells, we infected Arabidopsis Fluorescent Tagged Lines (FTL) with the virulent bacteria Pseudomonas syringae. We measured the meiotic recombination rate in two regions of chromosome 5, containing or not an NLR gene cluster. In all tested intervals, no significant difference in genetic recombination frequency between infected and control plants was observed. Although it has been reported that pathogen exposure can sometimes increase the frequency of recombinant progeny in plants, our findings suggest that meiotic recombination rate in Arabidopsis may be resilient to at least some pathogen attack. Alternative mechanisms are discussed.

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

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

•

Species-wide NLR diversity is high but not unlimited

•

A large fraction of NLR diversity is recovered with 40–50 accessions

•

Presence/absence variation in NLRs is widespread, resulting in a mosaic population

•

A high diversity of NLR-integrated domains favor known virulence targets

Mapping and DNA sequence characterisation of the Rysto locus conferring extreme virus resistance to potato cultivar 'White Lady'

Virus resistance genes carried by wild plant species are valuable resources for plant breeding. The Ry_sto gene, conferring a broad spectrum of durable resistance, originated from Solanum stoloniferum and was introgressed into several commercial potato cultivars, including ‘White Lady’, by classical breeding. Ry_sto was mapped to chromosome XII in potato, and markers used for marker-assisted selection in breeding programmes were identified. Nevertheless, there was no information on the identity of the Ry_sto gene. To begin to reveal the identification of Ry_sto, fine-scale genetic mapping was performed which, in combination with chromosome walking, narrowed down the locus of the gene to approximately 1 Mb. DNA sequence analysis of the locus identified six full-length NBS-LRR-type (short NLR-type) putative resistance genes. Two of them, designated TMV2 and TMV3, were similar to a TMV resistance gene isolated from tobacco and to Y-1, which co-segregates with Ry_adg, the extreme virus resistance gene originated from Solanum andigena and localised to chromosome XI. Furthermore, TMV2 of ‘White Lady’ was found to be 95% identical at the genomic sequence level with the recently isolated Ry_sto gene of the potato cultivar ‘Alicja’. In addition to the markers identified earlier, this work generated five tightly linked new markers which can serve potato breeding efforts for extreme virus resistance.

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.

Genomic insights from the first chromosome-scale assemblies of oat (Avena spp.) diploid species

BACKGROUND

Cultivated hexaploid oat (Common oat; Avena sativa) has held a significant place within the global crop community for centuries; although its cultivation has decreased over the past century, its nutritional benefits have garnered increased interest for human consumption. We report the development of fully annotated, chromosome-scale assemblies for the extant progenitor species of the As- and Cp-subgenomes, Avena atlantica and Avena eriantha respectively. The diploid Avena species serve as important genetic resources for improving common oat's adaptive and food quality characteristics.

RESULTS

The A. atlantica and A. eriantha genome assemblies span 3.69 and 3.78 Gb with an N50 of 513 and 535 Mb, respectively. Annotation of the genomes, using sequenced transcriptomes, identified ~ 50,000 gene models in each species-including 2965 resistance gene analogs across both species. Analysis of these assemblies classified much of each genome as repetitive sequence (~ 83%), including species-specific, centromeric-specific, and telomeric-specific repeats. LTR retrotransposons make up most of the classified elements. Genome-wide syntenic comparisons with other members of the Pooideae revealed orthologous relationships, while comparisons with genetic maps from common oat clarified subgenome origins for each of the 21 hexaploid linkage groups. The utility of the diploid genomes was demonstrated by identifying putative candidate genes for flowering time (HD3A) and crown rust resistance (Pc91). We also investigate the phylogenetic relationships among other A- and C-genome Avena species.

CONCLUSIONS

The genomes we report here are the first chromosome-scale assemblies for the tribe Poeae, subtribe Aveninae. Our analyses provide important insight into the evolution and complexity of common hexaploid oat, including subgenome origin, homoeologous relationships, and major intra- and intergenomic rearrangements. They also provide the annotation framework needed to accelerate gene discovery and plant breeding.

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.

Population Genetics of the Highly Polymorphic RPP8 Gene Family

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

Oomycetes Used in Arabidopsis Research

Arabidopsis plants in their natural environment are susceptible to infection by oomycete pathogens, in particular to downy mildew and white rust diseases. These naturally occurring infectious agents have imposed evolutionary pressures on Arabidopsis populations and are therefore highly relevant for the study of host-pathogen co-evolution. In addition, the study of oomycete diseases, including infections caused by several Phytophthora species, has led to many scientific discoveries on Arabidopsis immunity and disease. Herein, we describe the major oomycete species used for experiments on Arabidopsis, and how these pathosystems have been used to provide significant insights into mechanistic and evolutionary aspects of plant-oomycete interactions. We also highlight understudied aspects of plant-oomycete interactions, as well as translational approaches, that can be productively addressed using the reference pathosystems described in this article.

Surviving in a Hostile World: Plant Strategies to Resist Pests and Diseases

As primary producers, plants are under constant pressure to defend themselves against potentially deadly pathogens and herbivores. In this review, we describe short- and long-term strategies that enable plants to cope with these stresses. Apart from internal immunological strategies that involve physiological and (epi)genetic modifications at the cellular level, plants also employ external strategies that rely on recruitment of beneficial organisms. We discuss these strategies along a gradient of increasing timescales, ranging from rapid immune responses that are initiated within seconds to (epi)genetic adaptations that occur over multiple plant generations. We cover the latest insights into the mechanistic and evolutionary underpinnings of these strategies and present explanatory models. Finally, we discuss how knowledge from short-lived model species can be translated to economically and ecologically important perennials to exploit adaptive plant strategies and mitigate future impacts of pests and diseases in an increasingly interconnected and changing world.

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

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Long-term balancing selection drives evolution of immunity genes in Capsella

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

The Arabidopsis RRM domain protein EDM3 mediates race-specific disease resistance by controlling H3K9me2-dependent alternative polyadenylation of RPP7 immune receptor transcripts

The NLR-receptor RPP7 mediates race-specific immunity in Arabidopsis. Previous screens for enhanced downy mildew (edm) mutants identified the co-chaperone SGT1b (EDM1) and the PHD-finger protein EDM2 as critical regulators of RPP7. Here, we describe a third edm mutant compromised in RPP7 immunity, edm3. EDM3 encodes a nuclear-localized protein featuring an RNA-recognition motif. Like EDM2, EDM3 promotes histone H3 lysine 9 dimethylation (H3K9me2) at RPP7. Global profiling of H3K9me2 showed EDM3 to affect this silencing mark at a large set of loci. Importantly, both EDM3 and EDM2 co-associate in vivo with H3K9me2-marked chromatin and transcripts at a critical proximal polyadenylation site of RPP7, where they suppress proximal transcript polyadeylation/termination. Our results highlight the complexity of plant NLR gene regulation, and establish a functional and physical link between a histone mark and NLR-transcript processing.

Plant NLRs: From discovery to application

Plants require a complex immune system to defend themselves against a wide range of pathogens which threaten their growth and development. The nucleotide-binding leucine-rich repeat proteins (NLRs) are immune sensors that recognize effectors delivered by pathogens. The first NLR was cloned more than twenty years ago. Since this initial discovery, NLRs have been described as key components of plant immunity responsible for pathogen recognition and triggering defense responses. They have now been described in most of the well-studied mulitcellular plant species, with most having large NLR repertoires. As research has progressed so has the understanding of how NLRs interact with their recognition substrates and how they in turn activate downstream signalling. It has also become apparent that NLR regulation occurs at the transcriptional, post-transcriptional, translational, and post-translational levels. Even before the first NLR was cloned, breeders were utilising such genes to increase crop performance. Increased understanding of the mechanistic details of the plant immune system enable the generation of plants resistant against devastating pathogens. This review aims to give an updated summary of the NLR field.

Intragenic Meiotic Crossovers Generate Novel Alleles with Transgressive Expression Levels

Meiotic recombination is an evolutionary force that generates new genetic diversity upon which selection can act. Whereas multiple studies have assessed genome-wide patterns of recombination and specific cases of intragenic recombination, few studies have assessed intragenic recombination genome-wide in higher eukaryotes. We identified recombination events within or near genes in a population of maize recombinant inbred lines (RILs) using RNA-sequencing data. Our results are consistent with case studies that have shown that intragenic crossovers cluster at the 5′ ends of some genes. Further, we identified cases of intragenic crossovers that generate transgressive transcript accumulation patterns, that is, recombinant alleles displayed higher or lower levels of expression than did nonrecombinant alleles in any of ∼100 RILs, implicating intragenic recombination in the generation of new variants upon which selection can act. Thousands of apparent gene conversion events were identified, allowing us to estimate the genome-wide rate of gene conversion at SNP sites (4.9 × 10⁻⁵). The density of syntenic genes (i.e., those conserved at the same genomic locations since the divergence of maize and sorghum) exhibits a substantial correlation with crossover frequency, whereas the density of nonsyntenic genes (i.e., those which have transposed or been lost subsequent to the divergence of maize and sorghum) shows little correlation, suggesting that crossovers occur at higher rates in syntenic genes than in nonsyntenic genes. Increased rates of crossovers in syntenic genes could be either a consequence of the evolutionary conservation of synteny or a biological process that helps to maintain synteny.

Structural, Functional, and Genomic Diversity of Plant NLR Proteins: An Evolved Resource for Rational Engineering of Plant Immunity

Plants employ a diverse intracellular system of NLR (nucleotide binding-leucine-rich repeat) innate immune receptors to detect pathogens of all types. These receptors represent valuable agronomic traits that plant breeders rely on to maximize yield in the face of devastating pathogens. Despite their importance, the mechanistic underpinnings of NLR-based disease resistance remain obscure. The rapidly increasing numbers of plant genomes are revealing a diverse array of NLR-type immune receptors. In parallel, mechanistic studies are describing diverse functions for NLR immune receptors. In this review, we intend to broadly describe how the structural, functional, and genomic diversity of plant immune receptors can provide a valuable resource for rational engineering of plant immunity.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4806-7) contains supplementary material, which is available to authorized users.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

NtTPN1: a RPP8-like R gene required for Potato virus Y-induced veinal necrosis in tobacco

Potato virus Y (PVY) is one of the most damaging viruses of tobacco. In particular, aggressive necrotic strains (PVYN ) lead to considerable losses in yield. The main source of resistance against PVY is linked to the va locus. However, va-overcoming PVY isolates inducing necrotic symptoms were observed in several countries. In this context, it is important to find va-independent protection strategies. In a previous study, the phenotyping of 162 tobacco varieties revealed 10 accessions that do not carry the va allele and do not exhibit typical PVYN -induced veinal necrosis. Despite the absence of necrotic symptoms, normal viral accumulation in these plants suggests a va-independent mechanism of tolerance to PVYN -induced systemic veinal necrosis. Fine mapping of the genetic determinant(s) was performed in a segregating F2 population. The tolerance trait is inherited as a single recessive gene, and allelism tests demonstrated that eight of the 10 tolerant varieties carry the same determinant. Anchoring the linkage map to the tobacco genome physical map allowed the identification of a RPP8-like R gene, called NtTPN1 (for Nicotiana tabacum Tolerance to PVY-induced Necrosis1), with the same single-nucleotide polymorphism in the eight tolerant accessions. Functional assays using homozygous NtTPN1 EMS mutants confirmed the role of NtTPN1 in the tolerance phenotype. PVYN -induced systemic veinal necrosis in tobacco likely represents an inefficient defense response with hypersensitive response-like characteristics. The identification of NtTPN1 opens breeding options to minimize the impact of emerging and so far uncontrolled va-breaking necrotic PVY isolates.

Functional diversification of structurally alike NLR proteins in plants

In due course of evolution many pathogens alter their effector molecules to modulate the host plants' metabolism and immune responses triggered upon proper recognition by the intracellular nucleotide-binding oligomerization domain containing leucine-rich repeat (NLR) proteins. Likewise, host plants have also evolved with diversified NLR proteins as a survival strategy to win the battle against pathogen invasion. NLR protein indeed detects pathogen derived effector proteins leading to the activation of defense responses associated with programmed cell death (PCD). In this interactive process, genome structure and plasticity play pivotal role in the development of innate immunity. Despite being quite conserved with similar biological functions in all eukaryotes, the intracellular NLR immune receptor proteins happen to be structurally distinct. Recent studies have made progress in identifying transcriptional regulatory complexes activated by NLR proteins. In this review, we attempt to decipher the intracellular NLR proteins mediated surveillance across the evolutionarily diverse taxa, highlighting some of the recent updates on NLR protein compartmentalization, molecular interactions before and after activation along with insights into the finer role of these receptor proteins to combat invading pathogens upon their recognition. Latest information on NLR sensors, helpers and NLR proteins with integrated domains in the context of plant pathogen interactions are also discussed.

A Suppressor/Avirulence Gene Combination in Hyaloperonospora arabidopsidis Determines Race Specificity in Arabidopsis thaliana

The pathosystem of Arabidopsis thaliana and diploid biotrophic oomycete Hyaloperonospora arabidopsidis (Hpa) has been a model for investigating the molecular basis of Flor's gene-for-gene hypothesis. The isolates Hpa-Noks1 and Hpa-Cala2 are virulent on Arabidopsis accession RMX-A02 whilst an F₁ generated from a cross between these two isolates was avirulent. The F₂ progeny segregated 3,1 (avirulent, virulent), indicating a single major effect AVR locus in this pathogen. SNP-based linkage mapping confirmed a single AVR locus within a 14 kb map interval containing two genes encoding putative effectors. The Hpa-Cala2 allele of one gene, designated H. arabidopsidiscryptic1 (HAC1), encodes a protein with a signal peptide and an RxLR/dEER motif, and triggers a defense response in RMX-A02. The second gene is heterozygous in Hpa-Cala2. One allele, designated Suppressor ofHAC1^Cala2 (S-HAC1^Cala2 ) encodes a protein with a signal peptide and a dKEE motif with no RxLR motif; the other allele (s-hac1^Cala2 ) encodes a protein with a signal peptide, a dEEE motif and is divergent in sequence from the S-HAC1^Cala2 allele. In selfed progeny from Hpa-Cala2, dominant S-HAC1^Cala2 allele carrying progeny correlates with virulence in RMX-A02, whereas homozygous recessive s-hac1^Cala2 carrying progeny were avirulent. Genetic investigations suggested other heterozygous suppressor loci might exist in the Hpa-Cala2 genome.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.

Modification of meiotic recombination by natural variation in plants

Meiosis is a specialized cell division that produces haploid gametes required for sexual reproduction. During the first meiotic division, homologous chromosomes pair and undergo reciprocal crossing over, which recombines linked sequence variation. Meiotic recombination frequency varies extensively both within and between species. In this review, we will examine the molecular basis of meiotic recombination rate variation, with an emphasis on plant genomes. We first consider cis modification caused by polymorphisms at the site of recombination, or elsewhere on the same chromosome. We review cis effects caused by mismatches within recombining joint molecules, the effect of structural hemizygosity, and the role of specific DNA sequence motifs. In contrast, trans modification of recombination is exerted by polymorphic loci encoding diffusible molecules, which are able to modulate recombination on the same and/or other chromosomes. We consider trans modifiers that act to change total recombination levels, hotspot locations, or interactions between homologous and homeologous chromosomes in polyploid species. Finally, we consider the significance of genetic variation that modifies meiotic recombination for adaptation and evolution of plant species.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Genetics of mycorrhizal symbiosis in winter wheat (Triticum aestivum)

Bread wheat (Triticum aestivum) is a major staple food and therefore of prime importance for feeding the Earth's growing population. Mycorrhiza is known to improve plant growth, but although extensive knowledge concerning the interaction between mycorrhizal fungi and plants is available, genotypic differences concerning the ability of wheat to form mycorrhizal symbiosis and quantitative trait loci (QTLs) involved in mycorrhization are largely unknown. Therefore, a diverse set of 94 bread wheat genotypes was evaluated with regard to root colonization by arbuscular mycorrhizal fungi. In order to identify genomic regions involved in mycorrhization, these genotypes were analyzed using the wheat 90k iSelect chip, resulting in 17 823 polymorphic mapped markers, which were used in a genome-wide association study. Significant genotypic differences (P < 0.0001) were detected in the ability to form symbiosis and 30 significant markers associated with root colonization, representing six QTL regions, were detected on chromosomes 3A, 4A and 7A, and candidate genes located in these QTL regions were proposed. The results reported here provide key insights into the genetics of root colonization by mycorrhizal fungi in wheat.

Comparative Transcriptome Analyses of Resistant and Susceptible Near-Isogenic Wheat Lines following Inoculation with Blumeria graminis f. sp. tritici

Powdery mildew is one of the most important diseases of wheat. In this study, the leaf RNA samples of wheat NILs carrying powdery mildew resistant and susceptible Pm2 alleles (L031 and Chancellor) and its F1 hybrid at two time points (16 h and 96 h postinoculation) were used for RNA-seq analysis. We carry comparison between similar materials at different times and between different materials at same times. The overlapping DEGs between the dominant phenotypes (L031 and F1 hybrid) and the recessive phenotype (Chancellor) were 1028 and 2214 DEGs, which were clearly lower than those between the dominant and recessive parents and thus could provide relatively accurate and valuable information. GO and KEGG enrichment analysis of DEGs revealed that other than the expected defense-related genes, differential up- and downregulation of genes from many other signaling networks were also involved. Comparative transcriptome analysis also revealed that early-stage postinoculation is important and suitable time points to study expression profiles and signaling pathways of resistance-related genes following fungal inoculation. qRT-PCR analyses showed highly consistent expression patterns of genes with RNA-seq data. The results will aid in the identification of genes and signaling pathways involved in powdery mildew response in wheat.