Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean

Genome-Wide Association Study of Bacterial Blight Resistance in Soybean

Bacterial blight caused by Pseudomonas syringae pv. glycinea (Psg) is a widespread foliar disease. Although four Resistance to Pseudomonas syringae pv. glycinea (Rpg) 1 to 4 (Rpg1∼4) genes that have been observed to segregate in a Mendelian pattern have been reported to confer resistance to Psg in soybean, the genetic basis of quantitative resistance to bacterial blight in soybean remains unclear. In the present study, the Psg resistance of two soybean association panels consisting of 573 and 213 lines, respectively, was phenotyped in multiple environments in 2014 to 2016. Genome-wide association study was performed using two models, FarmCPU and BLINK, to identify Psg resistance loci. A total of 40 soybean varieties with high level of Psg resistance were identified, and 14 quantitative trait loci (QTLs) were detected on 12 soybean chromosomes. These QTLs were identified for the first time. The majority of the QTLs were detected only in one or the other association panels, while qRPG-18-1 was detected in both association panels for at least one growing season. A total of 46 candidate Psg resistance genes were identified from the qRpg_13_1, qRPG-15-1, and qRPG-18-1 loci based on gene function annotation. In addition, we found the genomic region covering rpg1-b and rpg1-r harbored the synteny with a genomic region on chromosome 15 and identified 16 nucleotide binding site-leucine-rich repeat (NBS-LRR) genes as the candidate Psg resistance genes from the synteny blocks. This study provides new information for dissecting the genetic control of Psg resistance in soybean.

Genome-Wide Association Study for Resistance to Phytophthora sojae in Soybean [Glycine max (L.) Merr.]

Phytophthora sojae (Kauffman and Gerdemann) is an oomycete pathogen that threatens soybean (Glycine max L.) production worldwide. The development of soybean cultivars with resistance to this pathogen is of paramount importance for the sustainable management of the disease. The objective of this study was to identify genomic regions associated with resistance to P. sojae isolate 40468 through genome-wide association analyses of 983 soybean germplasms. To elucidate the genetic basis of resistance, three statistical models were employed: the compressed mixed linear model (CMLM), Bayesian-information and linkage disequilibrium iteratively nested keyway (BLINK), and fixed and random model circulating probability unification (FarmCPU). The three models consistently identified a genomic region (3.8-5.3 Mbp) on chromosome 3, which has been previously identified as an Rps cluster. A total of 18 single nucleotide polymorphisms demonstrated high statistical significance across all three models, which were distributed in eight linkage disequilibrium (LD) blocks within the aforementioned interval. Of the eight, LD3-2 exhibited the discernible segregation of phenotypic reactions by haplotype. Specifically, over 93% of accessions with haplotypes LD3-2-F or LD3-2-G displayed resistance, whereas over 91% with LD3-2-A, LD3-2-C, or LD3-2-D exhibited susceptibility. Furthermore, the BLINK and FarmCPU models identified new genomic variations significantly associated with the resistance on several other chromosomes, indicating that the resistance observed in this panel was due to the presence of different alleles of multiple Rps genes. These findings underscore the necessity for robust statistical models to accurately detect true marker-trait associations and provide valuable insights into soybean genetics and breeding.

Evolutionary analysis of TIR- and non-TIR-NBS-LRR disease resistance genes in wild strawberries

Introduction

NBS-LRR genes (NLRs) are the most extensive category of plant resistance genes (R genes) and play a crucial role in pathogen defense. Understanding the diversity and evolutionary dynamics of NLRs in different plant species is essential for improving disease resistance. This study investigates the NLR gene family in eight diploid wild strawberry species to explore their structural characteristics, evolutionary relationships, and potential for enhancing disease resistance.

Methods

We conducted a comprehensive genome-wide identification and structural analysis of NLRs across eight diploid wild strawberry species. Phylogenetic analysis was performed to examine the relationships between TIR-NLRs (TNLs), Non-TIR-NLRs (non-TNLs), CC-NLRs (CNLs), and RPW8-NLRs (RNLs). Gene structures were compared, and gene expression was profiled across different NLR subfamilies. Additionally, in vitro leaf inoculation assays with Botrytis cinerea were performed to assess the resistance of various strawberry species.

Results

Our analysis revealed that non-TNLs constitute over 50% of the NLR gene family in all eight strawberry species, surpassing the proportion of TNLs. Phylogenetic analysis showed that TNLs diverged into two subclades: one grouping with CNLs and the other closely related to RNLs. A significantly higher number of non-TNLs were under positive selection compared to TNLs, indicating their rapid diversification. Gene structure analysis demonstrated that non-TNLs have shorter gene structures than TNLs and exhibit higher expression levels, particularly RNLs. Notably, non-TNLs showed dominant expression under both normal and infected conditions. In vitro leaf inoculation assays revealed that Fragaria pentaphylla and Fragaria nilgerrensis, which have the highest proportion of non-TNLs, exhibited significantly greater resistance to Botrytis cinerea compared to Fragaria vesca, which has the lowest proportion of non-TNLs.

Discussion

The findings of this study provide important insights into the evolutionary dynamics of NLRs in strawberries, particularly the significant role of non-TNLs in pathogen defense. The rapid diversification and higher expression levels of non-TNLs suggest their potential contribution to enhanced disease resistance. This research highlights the value of non-TNLs in strawberry breeding programs aimed at improving resistance to pathogens such as Botrytis cinerea.

The landscape of sequence variations between resistant and susceptible hot peppers to predict functional candidate genes against bacterial wilt disease

BACKGROUND

Bacterial wilt (BW), caused by Ralstonia solanacearum (Ral), results in substantial yield losses in pepper crops. Developing resistant pepper varieties through breeding is the most effective strategy for managing BW. To achieve this, a thorough understanding of the genetic information connected with resistance traits is essential. Despite identifying three major QTLs for bacterial wilt resistance in pepper, Bw1 on chromosome 8, qRRs-10.1 on chromosome 10, and pBWR-1 on chromosome 1, the genetic information of related BW pepper varieties has not been sufficiently studied. Here, we resequenced two pepper inbred lines, C. annuum 'MC4' (resistant) and C. annuum 'Subicho' (susceptible), and analyzed genomic variations through SNPs and Indels to identify candidate genes for BW resistance.

RESULTS

An average of 139.5 Gb was generated among the two cultivars, with coverage ranging from 44.94X to 46.13X. A total of 8,815,889 SNPs was obtained between 'MC4' and 'Subicho'. Among them, 31,190 (0.35%) were non-synonymous SNPs (nsSNPs) corresponding to 10,926 genes, and these genes were assigned to 142 Gene Ontology (GO) terms across the two cultivars. We focused on three known BW QTL regions by identifying genes with sequence variants through gene set enrichment analysis and securing those belonging to high significant GO terms. Additionally, we found 310 NLR genes with nsSNP variants between 'MC4' (R) and 'Subicho' (S) within these regions. Also, we performed an Indel analysis on these genes. By integrating all this data, we identified eight candidate BW resistance genes, including two NLR genes with nsSNPs variations in qRRs-10.1 on chromosome 10.

CONCLUSION

We identified genomic variations in the form of SNPs and Indels by re-sequencing two pepper cultivars with contrasting traits for bacterial wilt. Specifically, the four genes associated with pBWR-1 and Bw1 that exhibit both nsSNP and Indel variations simultaneously in 'Subicho', along with the two NLR genes linked to qRRs-10.1, which are known for their direct involvement in immune responses, are identified as most likely BW resistance genes. These variants in leading candidate genes associated with BW resistance can be used as important markers for breeding pepper varieties.

Genome-Wide Association Studies and QTL Mapping Reveal a New Locus Associated with Resistance to Bacterial Pustule Caused by Xanthomonas citri pv. glycines in Soybean

Bacterial pustule (BP), caused by Xanthomonas citri pv. glycines, is an important disease that, under favorable conditions, can drastically affect soybean production. We performed a genome-wide association study (GWAS) with a panel containing Brazilian and American cultivars, which were screened qualitatively and quantitatively against two Brazilian X. citri isolates (IBS 333 and IBS 327). The panel was genotyped using a genotyping by sequencing (GBS) approach, and we identified two main new regions in soybeans associated with X. citri resistance on chromosomes 6 (IBS 333) and 18 (IBS 327), different from the traditional rxp gene located on chromosome 17. The region on chromosome 6 was also detected by QTL mapping using a biparental cross between Williams 82 (R) and PI 416937 (S), showing that Williams 82 has another recessive resistance gene besides rxp, which was also detected in nine BP-resistant ancestors of the Brazilian cultivars (including CNS, S-100), based on haplotype analysis. Furthermore, we identified additional SNPs in strong LD (0.8) with peak SNPs by exploring variation available in WGS (whole genome sequencing) data among 31 soybean accessions. In these regions in strong LD, two candidate resistance genes were identified (Glyma.06g311000 and Glyma.18g025100) for chromosomes 6 and 18, respectively. Therefore, our results allowed the identification of new chromosomal regions in soybeans associated with BP disease, which could be useful for marker-assisted selection and will enable a reduction in time and cost for the development of resistant cultivars.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

The soybean plasma membrane GmDR1 protein conferring broad-spectrum disease and pest resistance regulates several receptor kinases and NLR proteins

Overexpression of Glycine max disease resistant 1 (GmDR1) exhibits broad-spectrum resistance against Fusarium virguliforme, Heterodera glycines (soybean cyst nematode), Tetranychus urticae (Koch) (spider mites), and Aphis glycines Matsumura (soybean aphids) in soybean. To understand the mechanisms of broad-spectrum immunity mediated by GmDR1, the transcriptomes of a strong and a weak GmDR1-overexpressor following treatment with chitin, a pathogen- and pest-associated molecular pattern (PAMP) common to these organisms, were investigated. The strong and weak GmDR1-overexpressors exhibited altered expression of 6098 and 992 genes, respectively, as compared to the nontransgenic control following chitin treatment. However, only 192 chitin- and 115 buffer-responsive genes exhibited over two-fold changes in expression levels in both strong and weak GmDR1-overexpressors as compared to the control. MapMan analysis of the 192 chitin-responsive genes revealed 64 biotic stress-related genes, of which 53 were induced and 11 repressed as compared to the control. The 53 chitin-induced genes include nine genes that encode receptor kinases, 13 encode nucleotide-binding leucine-rich repeat (NLR) receptor proteins, seven encode WRKY transcription factors, four ethylene response factors, and three MYB-like transcription factors. Investigation of a subset of these genes revealed three receptor protein kinases, seven NLR proteins, and one WRKY transcription factor genes that are induced following F. virguliforme and H. glycines infection. The integral plasma membrane GmDR1 protein most likely recognizes PAMPs including chitin and activates transcription of genes encoding receptor kinases, NLR proteins and defense-related genes. GmDR1 could be a pattern recognition receptor that regulates the expression of several NLRs for expression of PAMP-triggered immunity and/or priming the effector triggered immunity.

Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone)

MAIN CONCLUSION

The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Genome-wide association study of powdery mildew resistance in cultivated soybean from Northeast China

Powdery mildew (PMD), caused by the pathogen Microsphaera diffusa, leads to substantial yield decreases in susceptible soybean under favorable environmental conditions. Effective prevention of soybean PMD damage can be achieved by identifying resistance genes and developing resistant cultivars. In this study, we genotyped 331 soybean germplasm accessions, primarily from Northeast China, using the SoySNP50K BeadChip, and evaluated their resistance to PMD in a greenhouse setting. To identify marker-trait associations while effectively controlling for population structure, we conducted genome-wide association studies utilizing factored spectrally transformed linear mixed models, mixed linear models, efficient mixed-model association eXpedited, and compressed mixed linear models. The results revealed seven single nucleotide polymorphism (SNP) loci strongly associated with PMD resistance in soybean. Among these, one SNP was localized on chromosome (Chr) 14, and six SNPs with low linkage disequilibrium were localized near or in the region of previously mapped genes on Chr 16. In the reference genome of Williams82, we discovered 96 genes within the candidate region, including 17 resistance (R)-like genes, which were identified as potential candidate genes for PMD resistance. In addition, we performed quantitative real-time reverse transcription polymerase chain reaction analysis to evaluate the gene expression levels in highly resistant and susceptible genotypes, focusing on leaf tissues collected at different times after M. diffusa inoculation. Among the examined genes, three R-like genes, including Glyma.16G210800, Glyma.16G212300, and Glyma.16G213900, were identified as strong candidates associated with PMD resistance. This discovery can significantly enhance our understanding of soybean resistance to PMD. Furthermore, the significant SNPs strongly associated with resistance can serve as valuable markers for genetic improvement in breeding M. diffusa-resistant soybean cultivars.

Oligogenic Control of Quantitative Resistance Against Powdery Mildew Revealed in Portuguese Common Bean Germplasm

Common bean (Phaseolus vulgaris L.) is one of the most important food legumes worldwide, and its production is severely affected by fungal diseases such as powdery mildew. Portugal has a diverse germplasm, with accessions of Andean, Mesoamerican, and admixed origin, making it a valuable resource for common bean genetic studies. In this work, we evaluated the response of a Portuguese collection of 146 common bean accessions to Erysiphe diffusa infection, observing a wide range of disease severity and different levels of compatible and incompatible reactions, revealing the presence of different resistance mechanisms. We identified 11 incompletely hypersensitive resistant and 80 partially resistant accessions. We performed a genome-wide association study to clarify its genetic control, resulting in the identification of eight disease severity-associated single-nucleotide polymorphisms, spread across chromosomes Pv03, Pv09, and Pv10. Two of the associations were unique to partial resistance and one to incomplete hypersensitive resistance. The proportion of variance explained by each association varied between 15 and 86%. The absence of a major locus, together with the relatively small number of loci controlling disease severity, suggested an oligogenic inheritance of both types of resistance. Seven candidate genes were proposed, including a disease resistance protein (toll interleukin 1 receptor-nucleotide binding site-leucine-rich repeat class), an NF-Y transcription factor complex component, and an ABC-2 type transporter family protein. This work contributes with new resistance sources and genomic targets valuable to develop selection molecular tools and support powdery mildew resistance precision breeding in common bean.

Comparative analysis of soybean transcriptional profiles reveals defense mechanisms involved in resistance against Diaporthe caulivora

Soybean stem canker (SSC) caused by the fungal pathogen Diaporthe caulivora is an important disease affecting soybean production worldwide. However, limited information related to the molecular mechanisms underlying soybean resistance to Diaporthe species is available. In the present work, we analyzed the defense responses to D. caulivora in the soybean genotypes Williams and Génesis 5601. The results showed that compared to Williams, Génesis 5601 is more resistant to fungal infection evidenced by significantly smaller lesion length, reduced disease severity and pathogen biomass. Transcriptional profiling was performed in untreated plants and in D. caulivora-inoculated and control-treated tissues at 8 and 48 h post inoculation (hpi). In total, 2.322 and 1.855 genes were differentially expressed in Génesis 5601 and Williams, respectively. Interestingly, Génesis 5601 exhibited a significantly higher number of upregulated genes compared to Williams at 8 hpi, 1.028 versus 434 genes. Resistance to D. caulivora was associated with defense activation through transcriptional reprogramming mediating perception of the pathogen by receptors, biosynthesis of phenylpropanoids, hormone signaling, small heat shock proteins and pathogenesis related (PR) genes. These findings provide novel insights into soybean defense mechanisms leading to host resistance against D. caulivora, and generate a foundation for the development of resistant SSC varieties within soybean breeding programs.

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.

Genome-Wide Association Analysis for Resistance to Coniothyrium glycines Causing Red Leaf Blotch Disease in Soybean

Soybean is a high oil and protein-rich legume with several production constraints. Globally, several fungi, viruses, nematodes, and bacteria cause significant yield losses in soybean. Coniothyrium glycines (CG), the causal pathogen for red leaf blotch disease, is the least researched and causes severe damage to soybean. The identification of resistant soybean genotypes and mapping of genomic regions associated with resistance to CG is critical for developing improved cultivars for sustainable soybean production. This study used single nucleotide polymorphism (SNP) markers generated from a Diversity Arrays Technology (DArT) platform to conduct a genome-wide association (GWAS) analysis of resistance to CG using 279 soybean genotypes grown in three environments. A total of 6395 SNPs was used to perform the GWAS applying a multilocus model Fixed and random model Circulating Probability Unification (FarmCPU) with correction of the population structure and a statistical test p-value threshold of 5%. A total of 19 significant marker-trait associations for resistance to CG were identified on chromosomes 1, 5, 6, 9, 10, 12, 13, 15, 16, 17, 19, and 20. Approximately 113 putative genes associated with significant markers for resistance to red leaf blotch disease were identified across soybean genome. Positional candidate genes associated with significant SNP loci-encoding proteins involved in plant defense responses and that could be associated with soybean defenses against CG infection were identified. The results of this study provide valuable insight for further dissection of the genetic architecture of resistance to CG in soybean. They also highlight SNP variants and genes useful for genomics-informed selection decisions in the breeding process for improving resistance traits in soybean.

A genome-wide association study and genomic prediction for Phakopsora pachyrhizi resistance in soybean

Soybean brown rust (SBR), caused by Phakopsora pachyrhizi, is a devastating fungal disease that threatens global soybean production. This study conducted a genome-wide association study (GWAS) with seven models on a panel of 3,082 soybean accessions to identify the markers associated with SBR resistance by 30,314 high quality single nucleotide polymorphism (SNPs). Then five genomic selection (GS) models, including Ridge regression best linear unbiased predictor (rrBLUP), Genomic best linear unbiased predictor (gBLUP), Bayesian least absolute shrinkage and selection operator (Bayesian LASSO), Random Forest (RF), and Support vector machines (SVM), were used to predict breeding values of SBR resistance using whole genome SNP sets and GWAS-based marker sets. Four SNPs, namely Gm18_57,223,391 (LOD = 2.69), Gm16_29,491,946 (LOD = 3.86), Gm06_45,035,185 (LOD = 4.74), and Gm18_51,994,200 (LOD = 3.60), were located near the reported P. pachyrhizi R genes, Rpp1, Rpp2, Rpp3, and Rpp4, respectively. Other significant SNPs, including Gm02_7,235,181 (LOD = 7.91), Gm02_7234594 (LOD = 7.61), Gm03_38,913,029 (LOD = 6.85), Gm04_46,003,059 (LOD = 6.03), Gm09_1,951,644 (LOD = 10.07), Gm10_39,142,024 (LOD = 7.12), Gm12_28,136,735 (LOD = 7.03), Gm13_16,350,701(LOD = 5.63), Gm14_6,185,611 (LOD = 5.51), and Gm19_44,734,953 (LOD = 6.02), were associated with abundant disease resistance genes, such as Glyma.02G084100, Glyma.03G175300, Glyma.04g189500, Glyma.09G023800, Glyma.12G160400, Glyma.13G064500, Glyma.14g073300, and Glyma.19G190200. The annotations of these genes included but not limited to: LRR class gene, cytochrome 450, cell wall structure, RCC1, NAC, ABC transporter, F-box domain, etc. The GWAS based markers showed more accuracies in genomic prediction than the whole genome SNPs, and Bayesian LASSO model was the ideal model in SBR resistance prediction with 44.5% ~ 60.4% accuracies. This study aids breeders in predicting selection accuracy of complex traits such as disease resistance and can shorten the soybean breeding cycle by the identified markers.

Field Pea (Pisum sativum) Germplasm Screening for Seedling Ascochyta Blight Resistance and Genome-Wide Association Studies Reveal Loci Associated with Resistance to Peyronellaea pinodes and Ascochyta koolunga

Ascochyta blight is a damaging disease that affects the stems, leaves, and pods of field pea (Pisum sativum) and impacts yield and grain quality. In Australia, field pea Ascochyta blight is primarily caused by the necrotrophic fungal species Peyronellaea pinodes and Ascochyta koolunga. In this study, we screened 1,276 Pisum spp. germplasm accessions in seedling disease assays with a mix of three isolates of P. pinodes and 641 accessions with three mixed isolates of A. koolunga (513 accessions were screened with both species). A selection of three P. sativum accessions with low disease scores for either pathogen, or in some cases both, were crossed with Australian field pea varieties PBA Gunyah and PBA Oura, and recombinant inbred line populations were made. Populations at the F3:4 and F4:5 generation were phenotyped for their disease response to P. pinodes and A. koolunga, and genotypes were determined using the diversity arrays technology genotyping method. Marker-trait associations were identified using a genome-wide association study approach. Trait-associated loci were mapped to the published P. sativum genome assembly, and candidate resistance gene analogues were identified in the corresponding genomic regions. One locus on chromosome 2 (LG1) was associated with resistance to P. pinodes, and the 8 Mb genomic region contains 156 genes, two of which are serine/threonine protein kinases, putatively contributing to the resistance trait. A second locus on chromosome 5 (LG3) was associated with resistance to A. koolunga, and the 35 Mb region contains 488 genes, of which five are potential candidate resistance genes, including protein kinases, a mitogen-activated protein kinase, and an ethylene-responsive protein kinase homolog.

Comprehensive Analysis and Functional Verification of the Pinus massoniana NBS-LRR Gene Family Involved in the Resistance to Bursaphelenchus xylophilus

Pinus massoniana Lamb. is a crucial timber and resin conifer in China, but its plantation industry is threatened by outbreaks of pine wilt disease (PWD) caused by Bursaphelenchus xylophilus (pinewood nematode; PWN). However, as of yet, there is no comprehensive analysis of NBS-LRR genes in P. massoniana involved in its defense against PWN. In this study, 507 NBS genes were identified in the transcriptome of resistant and susceptible P. masoniana inoculated with the PWN. The phylogenetic analysis and expression profiles of resistant and susceptible P. massoniana revealed that the up-regulated PmNBS-LRR97 gene was involved in conferring resistance to PWN. The results of real-time quantitative PCR (qRT-PCR) showed that PmNBS-LRR97 was significantly up-regulated after PWN infection, especially in the stems. Subcellular localization indicated that PmNBS-LRR97 located to the cell membrane. PmNBS-LRR97 significantly activated the expression of reactive oxygen species (ROS)-related genes in P. massoniana. In addition, the overexpression of PmNBS-LRR97 was capable of promoting the production of ROS, aiding in plant growth and development. In summary, PmNBS-LRR97 participates in the defense response to PWN and plays an active role in conferring resistance in P. massoniana. This finding provides new insight into the regulatory mechanism of the R gene in P. massoniana.

Molecular mechanisms of Phytophthora sojae avirulence effectors escaping host recognition

Phytophthora sojae is a well-known destructive oomycete pathogen, which causes soybean stem and root rot and poses a serious threat to global food security. Growing soybean cultivars with the appropriate resistance to P. sojae (Rps) genes are the primary management strategy to reduce losses. In most Phytophthora pathosystems, host resistance protein encoded by a specific R gene in the plant recognizes corresponding RxLR effector protein, encoded by an avirulence gene. This gene-for-gene relationship has been exploited to help breeders and agronomists deploy soybean cultivars. To date, 6 Rps genes have been incorporated into commercial soybean germplasm and trigger plant immunity in response to 8 P. sojae avirulence effectors. The incorporation of Rps genes in the soybean population creates selection pressure in favor of novel pathotypes of P. sojae. The 8 avirulence genes evolved to evade the host immune system, driven by genetic selection pressures. Understanding the evading strategies has important reference value for the prevention and control of Phytophthora stem and root rot. This investigation primarily highlights the research on the strategies of P. sojae avirulence effector evasion of host recognition, looking forward to creating durable resistance genes and thereby enabling successful disease management.

Genome-wide identification, characterization, and evolutionary analysis of NBS genes and their association with disease resistance in Musa spp

Banana is an important food crop that is susceptible to a wide range of pests and diseases that can reduce yield and quality. The primary objective of banana breeding programs is to increase disease resistance, which requires the identification of resistance (R) genes. Despite the fact that resistant sources have been identified in bananas, the genes, particularly the nucleotide-binding site (NBS) family, which play an important role in protecting plants against pathogens, have received little attention. As a result, this study included a thorough examination of the NBS disease resistance gene family's classification, phylogenetic analysis, genome organization, evolution, cis-elements, differential expression, regulation by microRNAs, and protein-protein interaction. A total of 116 and 43 putative NBS genes from M. acuminata and M. balbisiana, respectively, were identified and characterized, and were classified into seven sub-families. Structural analysis of NBS genes revealed the presence of signal peptides, their sub-cellular localization, molecular weight and pI. Eight commonly conserved motifs were found, and NBS genes were unevenly distributed across multiple chromosomes, with the majority of NBS genes being located in chr3 and chr1 of the A and B genomes, respectively. Tandem duplication occurrences have helped bananas' NBS genes spread throughout evolution. Transcriptome analysis of NBS genes revealed significant differences in expression between resistant and susceptible cultivars of fusarium wilt, eumusae leaf spot, root lesion nematode, and drought, implying that they can be used as candidate resistant genes. Ninety miRNAs were discovered to have targets in 104 NBS genes from the A genome, providing important insights into NBS gene expression regulation. Overall, this study offers a valuable genomic resource and understanding of the function and evolution of NBS genes in relation to rapidly evolving pathogens, as well as providing breeders with selection targets for fast-tracking breeding of banana varieties with more durable resistance to pathogens.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

Advances in S gene targeted genome-editing and its applicability to disease resistance breeding in selected Solanaceae crop plants

Genome-editing tools for the development of traits to tolerate abiotic and biotic adversaries are the recently devised breeding techniques revolutionizing molecular breeding by addressing the issues of rapidness and precision. To that end, disease resistance development by disrupting disease susceptibility genes (S genes) to intervene in the biological mechanism of pathogenicity has significantly improved the techniques of molecular breeding. Despite the achievements in genome-editing aimed at the intervention of the function of susceptibility determinants or gene regulatory elements, off-target effects associated with yield-related traits are still the main setbacks. The challenges are attributed to the complexity of the inheritance of traits controlled by pleiotropic genes. Therefore, a more rigorous genome-editing tool with ultra-precision and efficiency for the development of broad-spectrum and durable disease resistance applied to staple crop plants is of critical importance in molecular breeding programs. The main objective of this article is to review the most impressive progresses achieved in resistance breeding against the main diseases of three Solanaceae crops (potato, Solanum tuberosum; tomato, Solanum lycopersicum and pepper, Capsicum annuum) using genome-editing by disrupting the sequences of S genes, their promoters, or pathogen genes. In this paper, we discussed the complexity and applicability of genome-editing tools, summarized the main disease of Solanaceae crops, and compiled the recent reports on disease resistance developed by S-gene silencing and their off-target effects. Moreover, GO count and gene annotation were made for pooled S-genes from biological databases. Achievements and prospects of S-gene-based next-generation breeding technologies are also discussed.

Most S genes are membrane –anchored and are involved in infection and pre-penetration process

S gene-editing is less likely to cause an off-target effect

Gene-editing has been considered a more acceptable engineering tool

Editing S genes either from the pathogen or host ends has opened new possibilities

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

Molecular mechanisms underlying multi-level defense responses of horticultural crops to fungal pathogens

The horticultural industry helps to enrich and improve the human diet while contributing to growth of the agricultural economy. However, fungal diseases of horticultural crops frequently occur during pre- and postharvest periods, reducing yields and crop quality and causing huge economic losses and wasted food. Outcomes of fungal diseases depend on both horticultural plant defense responses and fungal pathogenicity. Plant defense responses are highly sophisticated and are generally divided into preformed and induced defense responses. Preformed defense responses include both physical barriers and phytochemicals, which are the first line of protection. Induced defense responses, which include innate immunity (pattern-triggered immunity and effector-triggered immunity), local defense responses, and systemic defense signaling, are triggered to counterstrike fungal pathogens. Therefore, to develop regulatory strategies for horticultural plant resistance, a comprehensive understanding of defense responses and their underlying mechanisms is critical. Recently, integrated multi-omics analyses, CRISPR-Cas9-based gene editing, high-throughput sequencing, and data mining have greatly contributed to identification and functional determination of novel phytochemicals, regulatory factors, and signaling molecules and their signaling pathways in plant resistance. In this review, research progress on defense responses of horticultural crops to fungal pathogens and novel regulatory strategies to regulate induction of plant resistance are summarized, and then the problems, challenges, and future research directions are examined.

Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

A cluster of atypical resistance genes in soybean confers broad-spectrum antiviral activity

Soybean mosaic virus (SMV) is a severe soybean (Glycine max) pathogen. Here we characterize a soybean SMV resistance cluster (SRC) that comprises five resistance (R) genes. SRC1 encodes a Toll/interleukin-1 receptor and nucleotide-binding site (TIR-NBS [TN]) protein, SRC4 and SRC6 encode TIR proteins with a short EF-hand domain, while SRC7 and SRC8 encode TNX proteins with a noncanonical basic secretory protein (BSP) domain at their C-termini. We mainly studied SRC7, which contains a noncanonical BSP domain and gave full resistance to SMV. SRC7 possessed broad-spectrum antiviral activity toward several plant viruses including SMV, plum pox virus, potato virus Y, and tobacco mosaic virus. The TIR domain alone was both necessary and sufficient for SRC7 immune signaling, while the NBS domain enhanced its activity. Nuclear oligomerization via the interactions of both TIR and NBS domains was essential for SRC7 function. SRC7 expression was transcriptionally inducible by SMV infection and salicylic acid (SA) treatment, and SA was required for SRC7 triggered virus resistance. SRC7 expression was posttranscriptionally regulated by miR1510a and miR2109, and the SRC7-miR1510a/miR2109 regulatory network appeared to contribute to SMV-soybean interactions in both resistant and susceptible soybean cultivars. In summary, we report a soybean R gene cluster centered by SRC7 that is regulated at both transcriptional and posttranscriptional levels, possesses a yet uncharacterized BSP domain, and has broad-spectrum antiviral activities. The SRC cluster is special as it harbors several functional R genes encoding atypical TIR-NBS-LRR (TNL) type R proteins, highlighting its importance in SMV-soybean interaction and plant immunity.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Mapping Major Disease Resistance Genes in Soybean by Genome-Wide Association Studies

Soybean is one of the most valuable agricultural crops in the world. Besides, this legume is constantly attacked by a wide range of pathogens (fungi, bacteria, viruses, and nematodes) compromising yield and increasing production costs. One of the major disease management strategies is the genetic resistance provided by single genes and quantitative trait loci (QTL). Identifying the genomic regions underlying the resistance against these pathogens on soybean is one of the first steps performed by molecular breeders. In the past, genetic mapping studies have been widely used to discover these genomic regions. However, over the last decade, advances in next-generation sequencing technologies and their subsequent cost decreasing led to the development of cost-effective approaches to high-throughput genotyping. Thus, genome-wide association studies applying thousands of SNPs in large sets composed of diverse soybean accessions have been successfully done. In this chapter, a comprehensive review of the majority of GWAS for soybean diseases published since this approach was developed is provided. Important diseases caused by Heterodera glycines, Phytophthora sojae, and Sclerotinia sclerotiorum have been the focus of the several GWAS. However, other bacterial and fungi diseases also have been targets of GWAS. As such, this GWAS summary can serve as a guide for future studies of these diseases. The protocol begins by describing several considerations about the pathogens and bringing different procedures of molecular characterization of them. Advice to choose the best isolate/race to maximize the discovery of multiple R genes or to directly map an effective R gene is provided. A summary of protocols, methods, and tools to phenotyping the soybean panel is given to several diseases. We also give details of options of DNA extraction protocols and genotyping methods, and we describe parameters of SNP quality to soybean data. Websites and their online tools to obtain genotypic and phenotypic data for thousands of soybean accessions are highlighted. Finally, we report several tricks and tips in Subheading 4, especially related to composing the soybean panel as well as generating and analyzing the phenotype data. We hope this protocol will be helpful to achieve GWAS success in identifying resistance genes on soybean.

Molecular interaction of charcoal rot pathogenesis in soybean: a complex interaction

Charcoal rot (CR) is a major disease of soybean, which is caused by Macrophomina phaseolina (Mp). Increasing temperatures and low rainfall in recent years have immensely benefitted the pathogen. Hence, the search for genetically acquired resistance to this pathogen is essential. The pathogen is a hemibiotroph, which germinates on the root surface and colonizes epidermal tissue. Several surface receptors initiate pathogenesis, followed by the secretion of various enzymes that provide entry to host tissue. Several enzymes and other converging cascades in the pathogen participate against host defensive responses. β-glucan of the fungal cell wall is recognized as MAMPs (microbe-associated molecular patterns) in plants, which trigger host immune responses. Kinase receptors, resistance, and pathogenesis-related genes correspond to host defense response. They work in conjunction with hormone-mediated defense pathway especially, the systemic acquired resistance, calcium-signaling, and production of phytoalexins. Due to its quantitative nature, limited QTLs have been identified in soybean for CR resistance. The present review attempts to provide a functional link between M. phaseolina pathogenicity and soybean responses. Elucidation of CR resistance responses would facilitate improved designing of breeding programs, and may help in the selection of corresponding genes to introgress CR resistant traits.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

Novel disease resistance gene paralogs created by CRISPR/Cas9 in soy

KEY MESSAGE

Novel disease resistance gene paralogues are generated by targeted chromosome cleavage of tandem duplicated NBS-LRR gene complexes and subsequent DNA repair in soybean. This study demonstrates accelerated diversification of innate immunity of plants using CRISPR. Nucleotide-binding-site-leucine-rich-repeat (NBS-LRR) gene families are key components of effector-triggered immunity. They are often arranged in tandem duplicated arrays in the genome, a configuration that is conducive to recombinations that will lead to new, chimeric genes. These rearrangements have been recognized as major sources of novel disease resistance phenotypes. Targeted chromosome cleavage by CRISPR/Cas9 can conceivably induce rearrangements and thus emergence of new resistance gene paralogues. Two NBS-LRR families of soy have been selected to demonstrate this concept: a four-copy family in the Rpp1 region (Rpp1L) and a large, complex locus, Rps1 with 22 copies. Copy-number variations suggesting large-scale, CRISPR/Cas9-mediated chromosome rearrangements in the Rpp1L and Rps1 complexes were detected in up to 58.8% of progenies of primary transformants using droplet-digital PCR. Sequencing confirmed development of novel, chimeric paralogs with intact open reading frames. These novel paralogs may confer new disease resistance specificities. This method to diversify innate immunity of plants by genome editing is readily applicable to other disease resistance genes or other repetitive loci.

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 and Analysis of NBS-LRR Genes in Actinidia chinensis Genome

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the most important disease resistance genes in plants. The genome sequence of kiwifruit (Actinidia chinensis) provides resources for the characterization of NBS-LRR genes and identification of new R-genes in kiwifruit. In the present study, we identified 100 NBS-LRR genes in the kiwifruit genome and they were grouped into six distinct classes based on their domain architecture. Of the 100 genes, 79 are truncated non-regular NBS-LRR genes. Except for 37 NBS-LRR genes with no location information, the remaining 63 genes are distributed unevenly across 18 kiwifruit chromosomes and 38.01% of them are present in clusters. Seventeen families of cis-acting elements were identified in the promoters of the NBS-LRR genes, including AP2, NAC, ERF and MYB. Pseudomonas syringae pv. actinidiae (pathogen of the kiwifruit bacterial canker) infection induced differential expressions of 16 detected NBS-LRR genes and three of them are involved in plant immunity responses. Our study provides insight of the NBS-LRR genes in kiwifruit and a resource for the identification of new R-genes in the fruit.

Identifying Plasmopara viticola resistance Loci in grapevine (Vitis amurensis) via genotyping-by-sequencing-based QTL mapping

Downy mildew, caused by Plasmopara viticola, is a major disease that affects grapevines, and a few resistance (R) genes have been identified thus far. In order to identify R genes, we investigated F₁ progeny from a cross between the downy mildew-resistant Vitis amurensis 'Shuang Hong' and the susceptible Vitis vinifera 'Cabernet Sauvignon'. The P. viticola-resistance of the progeny varied continuously and was segregated as a quantitative trait. Genotyping-by-sequencing was used to construct linkage maps. The integrated map spanned 1898.09 cM and included 5603 single nucleotide polymorphisms on 19 linkage groups (LGs). Linkage analysis identified three quantitative trait loci (QTLs) for P. viticola resistance: 22 (Rpv22) on LG 02, Rpv23 on LG15, and Rpv24 on LG18. The phenotypic variance contributed by these three QTLs ranged from 15.9 to 30.0%. qRT-PCR analysis of R-gene expression in these QTLs revealed a CC-NBS-LRR disease resistance gene RPP8, two LRR receptor-like serine/threonine-protein kinases, a serine/threonine-protein kinase BLUS1, a glutathione peroxidase 8, an ethylene-responsive transcription factor ERF038, and a transcription factor bZIP11 were induced by P. viticola, and these genes may play important role in P. viticola response.

Fine Mapping, Candidate Gene Identification and Co-segregating Marker Development for the Phytophthora Root Rot Resistance Gene RpsYD25

Phytophthora root rot (PRR) caused by Phytophthora sojae is a serious disease of soybean. The most effective disease-control strategy is to deploy resistant cultivars carrying Rps genes. Soybean cultivar Yudou25 can effectively resist pathotypes of P. sojae in China. Previous studies have mapped the Rps gene in Yudou25, RpsYD25, on chromosome 3. In this study, at first RpsYD25 was located between SSR markers Satt1k3 (2.2 cM) and BARCSOYSSR_03_0253 (4.5 cM) by using an F_2:3 population containing 165 families derived from Zaoshu18 and Yudou25. Then the recombination sites were identified in 1127 F_3:4 families derived from Zaoshu18 and Yudou25 using the developed PCR-based SNP, InDel and SSR markers, and RpsYD25 was finely mapped in the a 101.3 kb genomic region. In this region, a zinc ion binding and nucleic acid binding gene Glyma.03g034700 and two NBS-LRR genes Glyma.03g034800 and Glyma.03g034900 were predicted as candidate genes of RpsYD25, and five co-segregated SSR markers with RpsYD25 were identified and validated to be diagnostic markers. Combined with the resistance reaction to multiple P. sojae isolates, seven of 178 soybean genotypes were detected to contain RpsYD25 using the five co-segregated SSR markers. The soybean genotypes carrying RpsYD25 and the developed co-segregated markers can be effectively applied in the breeding for P. sojae resistance in China.

In silico guided structural and functional analysis of genes with potential involvement in resistance to coffee leaf rust: A functional marker based approach

Physiology-based differentiation of S_H genes and Hemileia vastatrix races is the principal method employed for the characterization of coffee leaf rust resistance. Based on the gene-for-gene theory, nine major rust resistance genes (S_H1-9) have been proposed. However, these genes have not been characterized at the molecular level. Consequently, the lack of molecular data regarding rust resistance genes or candidates is a major bottleneck in coffee breeding. To address this issue, we screened a BAC library with resistance gene analogs (RGAs), identified RGAs, characterized and explored for any S_H related candidate genes. Herein, we report the identification and characterization of a gene (gene 11), which shares conserved sequences with other S_H genes and displays a characteristic polymorphic allele conferring different resistance phenotypes. Furthermore, comparative analysis of the two RGAs belonging to CC-NBS-LRR revealed more intense diversifying selection in tomato and grape genomes than in coffee. For the first time, the present study has unveiled novel insights into the molecular nature of the S_H genes, thereby opening new avenues for coffee rust resistance molecular breeding. The characterized candidate RGA is of particular importance for further biological function analysis in coffee.

Fine mapping of a Phytophthora-resistance locus RpsGZ in soybean using genotyping-by-sequencing

BACKGROUND

Phytophthora root rot (PRR) caused by Phytophthora sojae (P. sojae) is one of the most serious limitations to soybean production worldwide. The identification of resistance gene(s) and their incorporation into elite varieties is an effective approach for breeding to prevent soybean from being harmed by this disease. A valuable mapping population of 228 F8:11 recombinant inbred lines (RILs) derived from a cross of the resistant cultivar Guizao1 and the susceptible cultivar BRSMG68 and a high-density genetic linkage map with an average distance of 0.81 centimorgans (cM) between adjacent bin markers in this population were used to map and explore candidate gene(s).

RESULTS

PRR resistance in Guizao1 was found to be controlled by a single Mendelian locus and was finely mapped to a 367.371-kb genomic region on chromosome 3 harbouring 19 genes, including 7 disease resistance (R)-like genes, in the reference Willliams 82 genome. Quantitative real-time PCR assays of possible candidate genes revealed that Glyma.03 g05300 was likely involved in PRR resistance.

CONCLUSIONS

These findings from the fine mapping of a novel Rps locus will serve as a basis for the cloning and transfer of resistance genes in soybean and the breeding of P. sojae-resistant soybean cultivars through marker-assisted selection.

Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita

Background

The Mi-1 gene was the first identified and cloned gene that provides resistance to root-knot nematodes (RKNs) in cultivated tomato. However, owing to its temperature sensitivity, this gene does not meet the need for breeding disease-resistant plants that grow under high temperature. In this study, Mi-3 was isolated from the wild species PI 126443 (LA3858) and was shown to display heat-stable resistance to RKNs. However, the mechanism that regulates this resistance remains unknown.

Results

In this study, 4760, 1024 and 137 differentially expressed genes (DEGs) were enriched on the basis of pairwise comparisons (34 °C vs. 25 °C) at 0 (before inoculation), 3 and 6 days post-inoculation (dpi), respectively. A total of 7035 DEGs were identified from line LA3858 in the respective groups under the different soil temperature treatments. At 3 dpi, most DEGs were enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to plant biotic responses, such as “plant-pathogen interaction” and “plant hormone signal transduction”. Significantly enriched DEGs were found to encode key proteins such as R proteins and heat-shock proteins (HSPs). Moreover, other DEGs were found to participate in Ca²⁺ signal transduction; the production of ROS; DEGs encoding transcription factors (TFs) from the bHLH, TGA, ERF, heat-shock transcription factor (HSF) and WRKY families were highly expressed, which contribute to be involved into the formation of phytohormones, such as salicylic acid (SA), jasmonic acid (JA) and ethylene (ET), the expression of most was upregulated at 3 dpi at the 25 °C soil temperature compared with the 34 °C soil temperature.

Conclusion

Taken together, the results of our study revealed reliable candidate genes from wild materials LA3858, that are related to Mi-3-mediate resistance to Meloidogyne incognita. A large number of vital pathways and DEGs were expressed specifically in accession LA3858 grown at 34 °C and 25 °C soil temperatures at 3 dpi. Upon infection by RKNs, pattern-recognition receptors (PRRs) specifically recognized conserved pathogen-associated molecular patterns (PAMPs) as a result of pathogen-triggered immunity (PTI), and the downstream defensive signal transduction pathway was likely activated through Ca²⁺ signal channels. The expression of various TFs was induced to synthesize phytohormones and activate R proteins related to resistance, resulting in the development of effector-triggered immunity (ETI). Last, a hypersensitive response in the roots occurred, which was probably induced by the accumulation of ROS.

Gene presence-absence variation associates with quantitative Verticillium longisporum disease resistance in Brassica napus

Although copy number variation (CNV) and presence-absence variation (PAV) have been discovered in selected gene families in most crop species, the global prevalence of these polymorphisms in most complex genomes is still unclear and their influence on quantitatively inherited agronomic traits is still largely unknown. Here we analyze the association of gene PAV with resistance of oilseed rape (Brassica napus) against the important fungal pathogen Verticillium longisporum, as an example for a complex, quantitative disease resistance in the strongly rearranged genome of a recent allopolyploid crop species. Using Single Nucleotide absence Polymorphism (SNaP) markers to efficiently trace PAV in breeding populations, we significantly increased the resolution of loci influencing V. longisporum resistance in biparental and multi-parental mapping populations. Gene PAV, assayed by resequencing mapping parents, was observed in 23-51% of the genes within confidence intervals of quantitative trait loci (QTL) for V. longisporum resistance, and high-priority candidate genes identified within QTL were all affected by PAV. The results demonstrate the prominent role of gene PAV in determining agronomic traits, suggesting that this important class of polymorphism should be exploited more systematically in future plant breeding.

Development of Molecular Marker Linked with Bacterial Fruit Blotch Resistance in Melon (Cucumis melo L.)

Bacterial fruit blotch (BFB) causes losses in melon marketable yield. However, until now, there has been no information about the genetic loci responsible for resistance to the disease or their pattern of inheritance. We determined the inheritance pattern of BFB resistance from a segregating population of 491 F₂ individuals raised by crossing BFB-resistant (PI 353814) and susceptible (PI 614596) parental accessions. All F₁ plants were resistant to Acidovorax citrulli strain KACC18782, and F₂ plants segregated with a 3:1 ratio for resistant and susceptible phenotypes, respectively, in a seedling bioassay experiment, indicating that BFB resistance is controlled by a monogenic dominant gene. In an investigation of 57 putative disease-resistance related genes across the melon genome, only the MELO3C022157 gene (encoding TIR-NBS-LRR domain), showing polymorphism between resistant and susceptible parents, revealed as a good candidate for further investigation. Cloning, sequencing and quantitative RT-PCR expression of the polymorphic gene MELO3C022157 located on chromosome 9 revealed multiple insertion/deletions (InDels) and single nucleotide polymorphisms (SNPs), of which the SNP A²⁰³⁵T in the second exon of the gene caused loss of the LRR domain and truncated protein in the susceptible accession. The InDel marker MB157-2, based on the large (504 bp) insertion in the first intron of the susceptible accession, was able to distinguish resistant and susceptible accessions among 491 F₂ and 22 landraces/inbred accessions with 98.17% and 100% detection accuracy, respectively. This novel PCR-based, co-dominant InDel marker represents a practical tool for marker-assisted breeding aimed at developing BFB-resistant melon accessions.

Genome-Wide Mining of Disease Resistance Gene Analogs Using Conserved Domains

The production of legume crop species is severely affected by disease, imposing a significant yield loss annually worldwide. Plant resistance gene analogs (RGAs) play specific roles in plant resistance responses, and their identification and subsequent application in breeding programs help to reduce this yield loss. RGAs contain conserved domains and motifs, which can be used for their identification and classification. Nucleotide-binding site-leucine-rich repeat (NLR), receptor like kinase (RLK), and receptor like protein (RLP) genes are the main types of RGAs. Computational identification and characterization of RGAs has been performed successfully among different plant species. Here, we explain the computational workflow for genome-wide RGA identification in legumes.

Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera

NBS-LRR comprises a large class of disease resistance (R) proteins that play a widespread role in plant protection against pathogens. In grapevine, powdery mildew cause significant losses in its productivity and efforts are being directed towards finding of resistance loci or genes imparting resistance/tolerance against such fungal diseases. In the present study, we performed genome-wide analysis of NBS-LRR genes during PM infection in grapevine. We identified 18, 23, 12, 16, 10, 10, 9, 20 and 14 differentially expressed NBS-LRR genes in response to PM infection in seven partially PM-resistant (DVIT3351.27, Husseine, Karadzhandal, Khalchili, Late vavilov, O34-16, Sochal) and 2 PM-susceptible (Carignan and Thompson seedless) V. vinifera accessions. Further, the identified sequences were characterized based on chromosomal locations, physicochemical properties, gene structure and motif analysis, and functional annotation by Gene Ontology (GO) mapping. The NBS-LRR genes responsive to powdery mildew could potentially be exploited to improve resistance in grapes.

Colonisation of Oncidium orchid roots by the endophyte Piriformospora indica restricts Erwinia chrysanthemi infection, stimulates accumulation of NBS-LRR resistance gene transcripts and represses their targeting micro-RNAs in leaves

BACKGROUND

Erwinia chrysanthemi (Ec) is a destructive pathogen which causes soft-rot diseases in diverse plant species including orchids. We investigated whether colonization of Oncidium roots by the endophytic fungus Piriformospora indica (Pi) restricts Ec-induced disease development in leaves, and whether this might be related to the regulation of nucleotide binding site-leucine rich repeat (NBS-LRR) Resistance (R) genes.

RESULTS

Root colonization of Oncidium stackings by Pi restricts progression of Ec-induced disease development in the leaves. Since Pi does not inhibit Ec growth on agar plates, we tested whether NBS-LRR R gene transcripts and the levels of their potential target miRNAs in Oncidium leaves might be regulated by Pi. Using bioinformatic tools, we first identified NBS-LRR R gene sequences from Oncidium, which are predicted to be targets of miRNAs. Among them, the expression of two R genes was repressed and the accumulation of several regulatory miRNA stimulated by Ec in the leaves of Oncidium plants. This correlated with the progression of disease development, jasmonic and salicylic acid accumulation, ethylene synthesis and H₂O₂ production after Ec infection of Oncidium leaves. Interestingly, root colonization by Pi restricted disease development in the leaves, and this was accompanied by higher expression levels of several defense-related R genes and lower expression level of their target miRNA.

CONCLUSION

Based on these data we propose that Pi controls the levels of NBS-LRR R mRNAs and their target miRNAs in leaves. This regulatory circuit correlates with the protection of Oncidium plants against Ec infection, and molecular and biochemical investigations will demonstrate in the future whether, and if so, to what extent these two observations are related to each other.

Ectopic Expression of Grapevine Gene VaRGA1 in Arabidopsis Improves Resistance to Downy Mildew and Pseudomonas syringae pv. tomato DC3000 But Increases Susceptibility to Botrytis cinerea

The protein family with nucleotide binding sites and leucine-rich repeat (NBS-LRR) in plants stimulates immune responses caused by effectors and can mediate resistance to hemi-biotrophs and biotrophs. In our previous study, a Toll-interleukin-1(TIR)-NBS-LRR gene cloned from Vitis amurensis "Shuanghong", VaRGA1, was induced by Plasmopara viticola and could improve the resistance of tobacco to Phytophthora capsici. In this study, VaRGA1 in "Shuanghong" was also induced by salicylic acid (SA), but inhibited by jasmonic acid (JA). To investigate whether VaRGA1 confers broad-spectrum resistance to pathogens, we transferred this gene into Arabidopsis and then treated with Hyaloperonospora arabidopsidis (Hpa), Botrytis cinerea (B. cinerea), and Pseudomonas syringae pv. tomato DC3000 (PstDC3000). Results showed that VaRGA1 improved transgenic Arabidopsis thaliana resistance to the biotrophic Hpa and hemi-biotrophic PstDC3000, but decreased resistance to the necrotrophic B. cinerea. Additionally, qPCR assays showed that VaRGA1 plays an important role in disease resistance by activating SA and inhibiting JA signaling pathways. A 1104 bp promoter fragment of VaRGA1 was cloned and analyzed to further elucidate the mechanism of induction of the gene at the transcriptional level. These results preliminarily confirmed the disease resistance function and signal regulation pathway of VaRGA1, and contributed to the identification of R-genes with broad-spectrum resistance function.

Identification of a Major QTL (qRRs-10.1) That Confers Resistance to Ralstonia solanacearum in Pepper (Capsicum annuum) Using SLAF-BSA and QTL Mapping

The soilborne pathogen Ralstonia solanacearum is the causal agent of bacterial wilt (BW), a major disease of pepper (Capsicum annuum). The genetic basis of resistance to this disease in pepper is not well known. This study aimed to identify BW resistance markers in pepper. Analysis of the dynamics of bioluminescent R. solanacearum colonization in reciprocal grafts of a resistant (BVRC 1) line and a susceptible (BVRC 25) line revealed that the resistant rootstock effectively suppressed the spreading of bacteria into the scion. The two clear-cut phenotypic distributions of the disease severity index in 440 F2 plants derived from BVRC 25 × BVRC 1 indicated that a major genetic factor as well as a few minor factors that control BW resistance. By specific-locus amplified fragment sequencing combined with bulked segregant analysis, two adjacent resistance-associated regions on chromosome 10 were identified. Quantitative trait (QTL) mapping revealed that these two regions belong to a single QTL, qRRs-10.1. The marker ID10-194305124, which reached a maximum log-likelihood value at 9.79 and accounted for 19.01% of the phenotypic variation, was located the closest to the QTL peak. A cluster of five predicted R genes and three defense-related genes, which are located in close proximity to the significant markers ID10-194305124 or ID10-196208712, are important candidate genes that may confer BW resistance in pepper.

Evolutionary balance between LRR domain loss and young NBS-LRR genes production governs disease resistance in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea L.) is an important oil and protein crop, but it has low disease resistance; therefore, it is important to reveal the number, sequence features, function, and evolution of genes that confer resistance. Nucleotide-binding site-leucine-rich repeats (NBS-LRRs) are resistance genes that are involved in response to various pathogens.

RESULTS

We identified 713 full-length NBS-LRRs in A. hypogaea cv. Tifrunner. Genetic exchange events occurred on NBS-LRRs in A. hypogaea cv. Tifrunner, which were detected in the same subgenomes and also found in different subgenomes. Relaxed selection acted on NBS-LRR proteins and LRR domains in A. hypogaea cv. Tifrunner. Using quantitative trait loci (QTL), we found that NBS-LRRs were involved in response to late leaf spot, tomato spotted wilt virus, and bacterial wilt in A. duranensis (2 NBS-LRRs), A. ipaensis (39 NBS-LRRs), and A. hypogaea cv. Tifrunner (113 NBS-LRRs). In A. hypogaea cv. Tifrunner, 113 NBS-LRRs were classified as 75 young and 38 old NBS-LRRs, indicating that young NBS-LRRs were involved in response to disease after tetraploidization. However, compared to A. duranensis and A. ipaensis, fewer LRR domains were found in A. hypogaea cv. Tifrunner NBS-LRR proteins, partly explaining the lower disease resistance of the cultivated peanut.

CONCLUSIONS

Although relaxed selection acted on NBS-LRR proteins and LRR domains, LRR domains were preferentially lost in A. hypogaea cv. Tifrunner compared to A. duranensis and A. ipaensis. The QTL results suggested that young NBS-LRRs were important for resistance against diseases in A. hypogaea cv. Tifrunner. Our results provid insight into the greater susceptibility of A. hypogaea cv. Tifrunner to disease compared to A. duranensis and A. ipaensis.

Resistance mechanisms involved in complex immunity of wheat against rust diseases

The review is devoted to the disclosure of the modern concept of plant immunity as a hierarchical system of plant host protection, controlled by combinations of major and minor resistance genes (loci). The “zigzag” model is described in detail for discussing the molecular bases of plant immunity with key concepts: pathogen-associated molecular patterns triggering innate immunity, ambivalent effectors causing susceptibility, but when interacting with resistance genes, a hypersensitive reaction or alternative defense mechanisms. There are three types of resistance in cereals: (1) basal resistance provided by plasma membrane-localized receptors proteins; (2) racespecific resistance provided by intracellular immune R-receptors; (3) partial resistance conferred by quantitative gene loci. The system ‘wheat (Triticum aestivum) – the fungus causing leaf rust (Puccinia triticina)’ is an interesting model for observing all the resistance mechanisms listed above, since the strategy of this pathogen is aimed at the constitutive use of host resources. The review focuses on known wheat genes responsible for various types of resistance to leaf rust: race-specific genes Lr1, Lr10, Lr19, and Lr21; adult resistance genes which are hypersensitive Lr12, Lr13, Lr22a, Lr22b, Lr35, Lr48, and Lr49; nonhypersensitive genes conferring partial resistance Lr34, Lr46, Lr67, and Lr77. The involvement of some wheat R-genes in pre-haustorial resistance to leaf rust has been discovered recently: Lr1, Lr3a, Lr9, LrB, Lr19, Lr21, Lr38. The presence of these genes in the genotype ensures the interruption of early pathogenesis through the following mechanisms: disorientation and branching of the germ tube; formation of aberrant fungal penetration structures (appressorium, substomatal vesicle); accumulation of callose in mesophyll cell walls. Breeding for immunity is accelerated by implementation of data on various mechanisms of wheat resistance to rust diseases, which are summarized in this review.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Genotypes and Genomic Regions Associated With Rhizoctonia solani Resistance in Common Bean

Rhizoctonia solani Kühn (teleomorph Thanatephorus cucumeris) is an important root rot pathogen of common bean (Phaseolus vulgaris L.). To uncover genetic factors associated with resistance to the pathogen, the Andean (ADP; n = 273) and Middle American (MDP; n = 279) diversity panels, which represent much of the genetic diversity known in cultivated common bean, were screened in the greenhouse using R. solani anastomosis group 2-2. Repeatability of the assay was confirmed by the response of five control genotypes. The phenotypic data for both panels were normally distributed. The resistance responses of ∼10% of the ADP (n = 28) and ∼6% of the MDP (n = 18) genotypes were similar or higher than that of the resistant control line VAX 3. A genome-wide association study (GWAS) was performed using ∼200k single nucleotide polymorphisms to discover genomic regions associated with resistance in each panel, For GWAS, the raw phenotypic score, and polynomial and binary transformation of the scores, were individually used as the input data. A major QTL peak was observed on Pv02 in the ADP, while a major QTL was observed on Pv01 with the MDP. These regions were associated with clusters of TIR-NB_ARC-LRR (TNL) gene models encoding proteins similar to known disease resistance genes. Other QTL, unique to each panel, were mapped within or adjacent to a gene model or cluster of related genes associated with disease resistance. This is a first case study that provides evidence for major as well as minor genes involved in resistance to R. solani in common bean. This information will be useful to integrate more durable root rot resistance in common bean breeding programs and to study the genetic mechanisms associated with root diseases in this important societal legume.

Genetic mapping of powdery mildew resistance genes in soybean by high-throughput genome-wide sequencing

KEY MESSAGE

The Mendelian locus conferring resistance to powdery mildew in soybean was precisely mapped using a combination of phenotypic screening, genetic analyses, and high-throughput genome-wide sequencing. Powdery mildew (PMD), caused by the fungus Microsphaera diffusa Cooke & Peck, leads to considerable yield losses in soybean [Glycine max (L.) Merr.] under favourable environmental conditions and can be controlled by identifying germplasm resources with resistance genes. In this study, resistance to M. diffusa among resistant varieties B3, Fudou234, and B13 is mapped as a single Mendelian locus using three mapping populations derived from crossing susceptible with resistant cultivars. The position of the PMD resistance locus in B3 is located between simple sequence repeat (SSR) markers GMES6959 and Satt_393 on chromosome 16, at genetic distances of 7.1 cM and 4.6 cM, respectively. To more finely map the PMD resistance gene, a high-density genetic map was constructed using 248 F8 recombinant inbred lines derived from a cross of Guizao1 × B13. The final map includes 3748 bins and is 3031.9 cM in length, with an average distance of 0.81 cM between adjacent markers. This genotypic analysis resulted in the precise delineation of the B13 PMD resistance locus to a 188.06-kb genomic region on chromosome 16 that harbours 28 genes, including 17 disease resistance (R)-like genes in the reference Williams 82 genome. Quantitative real-time PCR assays of possible candidate genes revealed differences in the expression levels of 9 R-like genes between the resistant and susceptible parents. These results provide useful information for marker-assisted breeding and gene cloning for PMD resistance.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.