Unique evolutionary pattern of numbers of gramineous NBS-LRR genes

Adaptation mechanism of three Impatiens species to different habitats based on stem morphology, lignin and MYB4 gene

BACKGROUND

Impatiens is an important genus with rich species of garden plants, and its distribution is extremely extensive, which is reflected in its diverse ecological environment. However, the specific mechanisms of Impatiens' adaptation to various environments and the mechanism related to lignin remain unclear.

RESULTS

Three representative Impatiens species,Impatiens chlorosepala (wet, low degree of lignification), Impatiens uliginosa (aquatic, moderate degree of lignification) and Impatiens rubrostriata (terrestrial, high degree of lignification), were selected and analyzed for their anatomical structures, lignin content and composition, and lignin-related gene expression. There are significant differences in anatomical parameters among the stems of three Impatiens species, and the anatomical structure is consistent with the determination results of lignin content. Furthermore, the thickness of the xylem and cell walls, as well as the ratio of cell wall thickness to stem diameter have a strong correlation with lignin content. The anatomical structure and degree of lignification in Impatiens can be attributed to the plant's growth environment, morphology, and growth rate. Our analysis of lignin-related genes revealed a negative correlation between the MYB4 gene and lignin content. The MYB4 gene may control the lignin synthesis in Impatiens by controlling the structural genes involved in the lignin synthesis pathway, such as HCT, C3H, and COMT. Nonetheless, the regulation pathway differs between species of Impatiens.

CONCLUSIONS

This study demonstrated consistency between the stem anatomy of Impatiens and the results obtained from lignin content and composition analyses. It is speculated that MYB4 negatively regulates the lignin synthesis in the stems of three Impatiens species by regulating the expression of structural genes, and its regulation mechanism appears to vary across different Impatiens species. This study analyses the variations among different Impatiens plants in diverse habitats, and can guide further molecular investigations of lignin biosynthesis in Impatiens.

High-resolution mapping of Ryd4Hb, a major resistance gene to Barley yellow dwarf virus from Hordeum bulbosum

KEY MESSAGE

We mapped Ryd4Hb in a 66.5 kbp interval in barley and dissociated it from a sublethality factor. These results will enable a targeted selection of the resistance in barley breeding. Virus diseases are causing high yield losses in crops worldwide. The Barley yellow dwarf virus (BYDV) complex is responsible for one of the most widespread and economically important viral diseases of cereals. While no gene conferring complete resistance (immunity) has been uncovered in the primary gene pool of barley, sources of resistance were searched and identified in the wild relative Hordeum bulbosum, representing the secondary gene pool of barley. One such locus, Ryd4Hb, has been previously introgressed into barley, and was allocated to chromosome 3H, but is tightly linked to a sublethality factor that prevents the incorporation and utilization of Ryd4Hb in barley varieties. To solve this problem, we fine-mapped Ryd4Hb and separated it from this negative factor. We narrowed the Ryd4Hb locus to a corresponding 66.5 kbp physical interval in the barley 'Morex' reference genome. The region comprises a gene from the nucleotide-binding and leucine-rich repeat immune receptor family, typical of dominant virus resistance genes. The closest homolog to this Ryd4Hb candidate gene is the wheat Sr35 stem rust resistance gene. In addition to the fine mapping, we reduced the interval bearing the sublethality factor to 600 kbp in barley. Aphid feeding experiments demonstrated that Ryd4Hb provides a resistance to BYDV rather than to its vector. The presented results, including the high-throughput molecular markers, will permit a more targeted selection of the resistance in breeding, enabling the use of Ryd4Hb in barley varieties.

Genome-wide identification of nucleotide-binding domain leucine-rich repeat (NLR) genes and their association with green peach aphid (Myzus persicae) resistance in peach

Resistance genes (R genes) are a class of genes that are immune to a wide range of diseases and pests. In planta, NLR genes are essential components of the innate immune system. Currently, genes belonging to NLR family have been found in a number of plant species, but little is known in peach. Here, 286 NLR genes were identified on peach genome by using their homologous genes in Arabidopsis thaliana as queries. These 286 NLR genes contained at least one NBS domain and LRR domain. Phylogenetic and N-terminal domain analysis showed that these NLRs could be separated into four subfamilies (I-IV) and their promoters contained many cis-elements in response to defense and phytohormones. In addition, transcriptome analysis showed that 22 NLR genes were up-regulated after infected by Green Peach Aphid (GPA), and showed different expression patterns. This study clarified the NLR gene family and their potential functions in aphid resistance process. The candidate NLR genes might be useful in illustrating the mechanism of aphid resistance in peach.

In silico structural-functional characterization of three differentially expressed resistance gene analogs identified in Dalbergia sissoo against dieback disease reveals their role in immune response regulation

Plant immunity includes enemy recognition, signal transduction, and defensive response against pathogens. We experimented to identify the genes that contribute resistance against dieback disease to Dalbergia sissoo, an economically important timber tree. In this study, we investigated the role of three differentially expressed genes identified in the dieback-induced transcriptome in Dalbergia sissoo. The transcriptome was probed using DOP-rtPCR analysis. The identified RGAs were characterized in silico as the contributors of disease resistance that switch on under dieback stress. Their predicted fingerprints revealed involvement in stress response. Ds-DbRCaG-02-Rga.a, Ds-DbRCaG-04-Rga.b, and Ds-DbRCaG-06-Rga.c showed structural homology with the Transthyretin-52 domain, EAL associated YkuI_C domain, and Src homology-3 domain respectively, which are the attributes of signaling proteins possessing a role in regulating immune responses in plants. Based on in-silico structural and functional characterization, they were predicted to have a role in immune response regulation in D. sissoo.

Exon shuffling potentiates a diverse repertoire of brown algal NB-ARC-TPR candidate immune receptor proteins via alternative splicing

Like other organisms, brown algae are subject to diseases caused by bacteria, fungi, and viruses. Brown algal immunity mechanisms are not well characterized; however, there is evidence suggesting that pathogen receptors exist in brown algae. One key protein family likely associated with brown algal innate immunity possesses an NB-ARC domain analogous to innate immune proteins in plants and animals. In this study, we conducted an extensive survey of NB-ARC genes in brown algae and obtained insights into the domain organization and evolutionary history of the encoded proteins. Our data show that brown algae possess an ancient NB-ARC-tetratricopeptide repeat (NB-TPR) domain architecture. We identified an N-terminal effector domain, the four-helix bundle, which was not previously found associated with NB-ARC domains. The phylogenetic tree including NB-ARC domains from all kingdoms of life suggests the three clades of brown algal NB-TPRs are likely monophyletic, whereas their TPRs seem to have distinct origins. One group of TPRs exhibit intense exon shuffling, with various alternative splicing and diversifying selection acting on them, suggesting exon shuffling is an important mechanism for evolving ligand-binding specificities. The reconciliation of gene duplication and loss events of the NB-ARC genes reveals that more independent gene gains than losses have occurred during brown algal evolution, and that tandem duplication has played a major role in the expansion of NB-ARC genes. Our results substantially enhance our understanding of the evolutionary history and exon shuffling mechanisms of the candidate innate immune repertoire of brown algae.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Comparative-genomic analysis reveals dynamic NLR gene loss and gain across Apiaceae species

Introduction: Nucleotide-binding leucine-rich repeat (NLR) genes play a crucial role in green plants' responding to various pathogens. Genome-scale evolutionary studies of NLR genes are important for discovering and applying functional NLR genes. However, little is known about the evolution of NLR genes in the Apiaceae family including agricultural and medical plants. Methods: In this study, comparative genomic analysis was performed in four Apiaceae species to trace the dynamic evolutionary patterns of NLR genes during speciation in this family. Results: The results revealed different number of NLR genes in these four Apiaceae species, namely, Angelica sinensis (95), Coriandrum sativum (183), Apium graveolens (153) and Daucus carota (149). Phylogenetic analysis demonstrated that NLR genes in these four species were derived from 183 ancestral NLR lineages and experienced different levels of gene-loss and gain events. The contraction pattern of the ancestral NLR lineages was discovered during the evolution of D. carota, whereas a different pattern of contraction after first expansion of NLR genes was observed for A. sinensis, C. sativum and A. graveolens. Discussion: Taken together, rapid and dynamic gene content variation has shaped evolutionary history of NLR genes in Apiaceae species.

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.

Genome-wide analysis of NBS-LRR genes in Rosaceae species reveals distinct evolutionary patterns

The nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes, one of the largest gene families in plants, are evolving rapidly and playing a critical role in plant resistance to pathogens. In this study, a genome-wide search in 12 Rosaceae genomes screened out 2188 NBS-LRR genes, with the gene number varied distinctively across different species. The reconciled phylogeny revealed 102 ancestral genes (7 RNLs, 26 TNLs, and 69 CNLs), which underwent independent gene duplication and loss events during the divergence of the Rosaceae. The NBS-LRR genes exhibited dynamic and distinct evolutionary patterns in the 12 Rosaceae species due to independent gene duplication/loss events, which resulted the discrepancy of NBS-LRR gene number among Rosaceae species. Specifically, Rubus occidentalis, Potentilla micrantha, Fragaria iinumae and Gillenia trifoliata, displayed a “first expansion and then contraction” evolutionary pattern; Rosa chinensis exhibited a “continuous expansion” pattern; F. vesca had a “expansion followed by contraction, then a further expansion” pattern, three Prunus species and three Maleae species shared a “early sharp expanding to abrupt shrinking” pattern. Overall, this study elucidated the dynamic and complex evolutionary patterns of NBS-LRR genes in the 12 Rosaceae species, and could assist further investigation of mechanisms driving these evolutionary patterns.

MusaRgeneDB: an online comprehensive database for disease resistance genes in Musa spp

Banana is one of the major food crops and its production is subject to many pests and diseases. Banana breeding exploits wild relatives and progenitor species for the introgression of resistant genes (R) into cultivated varieties to overcome these hurdles. With advances in sequencing technologies, whole-genome sequences are available for many Musa spp. and many of them are potential donors of disease resistance genes. Considering their potential role, R genes from these species were explored to develop an user-friendly open-access database that will be useful for studying and implementing disease resistance in bananas. MusaRgene database is complemented with complete details of 3598 R genes identified from eight Musa spp. and rice, Arabidopsis, sorghum along with its classification and separate modules on its expression under various stresses in resistant and susceptible cultivars and corresponding SSRs are also provided. This database can be regarded as the primary resource of information on R genes from bananas and their relatives. R genes from other allele mining studies are also incorporated which will enable the identification of its homolog in related Musa spp. MusaRgene database will aid in the identification of genes and markers associated, cloning of full-length R genes, and genetic transformation or gene editing of the R genes in susceptible cultivars. Multiple R genes can also be identified for pyramiding the genes to increase the level of resistance and durability. Overall, this database will facilitate the understanding of defense mechanisms in bananas against biotic or abiotic stresses leading to the development of promising disease-resistant varieties.

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.

Genomic Variations and Mutational Events Associated with Plant-Pathogen Interactions

Simple Summary

Plants, unlike animals, do not have defender cells or an adaptive immune system. Instead, plants rely on each cell’s innate immunity and systemic signals emitted from infection sites. On the other hand, not all plants, even within the same species, are genetically identical, and their genetic backgrounds determine how well they respond to stress factors. Through evolution, plants have acquired various defense mechanisms that play important roles in the never-ending fight between plants and pathogens. Genetic variation in relation to plant disease resistance can thus be contextualized to provide new insights into these defense mechanisms and evolutionary processes that lead to resistance to pathogens. By focusing on genetic variations and mutational events linked with plant–pathogen interactions, the paper explores how genome compartments facilitate plant and pathogen evolutionary processes.

Abstract

Phytopathologists are actively researching the molecular basis of plant–pathogen interactions. The mechanisms of responses to pathogens have been studied extensively in model crop plant species and natural populations. Today, with the rapid expansion of genomic technologies such as DNA sequencing, transcriptomics, proteomics, and metabolomics, as well as the development of new methods and protocols, data analysis, and bioinformatics, it is now possible to assess the role of genetic variation in plant–microbe interactions and to understand the underlying molecular mechanisms of plant defense and microbe pathogenicity with ever-greater resolution and accuracy. Genetic variation is an important force in evolution that enables organisms to survive in stressful environments. Moreover, understanding the role of genetic variation and mutational events is essential for crop breeders to produce improved cultivars. This review focuses on genetic variations and mutational events associated with plant–pathogen interactions and discusses how these genome compartments enhance plants’ and pathogens’ evolutionary processes.

Plant autoimmunity-fresh insights into an old phenomenon

The plant immune system is well equipped to ward off the attacks of different types of phytopathogens. It primarily relies on two types of immune sensors-plasma membrane-resident receptor-like kinases and intracellular nucleotide-binding domain leucine-rich repeat (NLRs) receptors that engage preferentially in pattern- and effector-triggered immunity, respectively. Delicate fine-tuning, in particular of the NLR-governed branch of immunity, is key to prevent inappropriate and deleterious activation of plant immune responses. Inadequate NLR allele constellations, such as in the case of hybrid incompatibility, and the mis-activation of NLRs or the absence or modification of proteins guarded by these NLRs can result in the spontaneous initiation of plant defense responses and cell death-a phenomenon referred to as plant autoimmunity. Here, we review recent insights augmenting our mechanistic comprehension of plant autoimmunity. The recent findings broaden our understanding regarding hybrid incompatibility, unravel candidates for proteins likely guarded by NLRs and underline the necessity for the fine-tuning of NLR expression at various levels to avoid autoimmunity. We further present recently emerged tools to study plant autoimmunity and draw a cross-kingdom comparison to the role of NLRs in animal autoimmune conditions.

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.

Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Frequent Gene Duplication/Loss Shapes Distinct Evolutionary Patterns of NLR Genes in Arecaceae Species

Nucleotide-binding leucine-rich repeat (NLR) genes play a key role in plant immune responses and have co-evolved with pathogens since the origin of green plants. Comparative genomic studies on the evolution of NLR genes have been carried out in several angiosperm lineages. However, most of these lineages come from the dicot clade. In this study, comparative analysis was performed on NLR genes from five Arecaceae species to trace the dynamic evolutionary pattern of the gene family during species speciation in this monocot lineage. The results showed that NLR genes from the genomes of Elaeis guineensis (262), Phoenix dactylifera (85), Daemonorops jenkinsiana (536), Cocos nucifera (135) and Calamus simplicifolius (399) are highly variable. Frequent domain loss and alien domain integration have occurred to shape the NLR protein structures. Phylogenetic analysis revealed that NLR genes from the five genomes were derived from dozens of ancestral genes. D. jenkinsiana and E. guineensis genomes have experienced “consistent expansion” of the ancestral NLR lineages, whereas a pattern of “first expansion and then contraction” of NLR genes was observed for P. dactylifera, C. nucifera and C. simplicifolius. The results suggest that rapid and dynamic gene content and structure variation have shaped the NLR profiles of Arecaceae species.

Identification and Characterization of NBS Resistance Genes in Akebia trifoliata

Akebia trifoliata is an important multiuse perennial plant that often suffers attacks from various pathogens due to its long growth cycle, seriously affecting its commercial value. The absence of research on the resistance (R) genes of A. trifoliata has greatly limited progress in the breeding of resistant varieties. Genes encoding proteins containing nucleotide binding sites (NBSs) and C-terminal leucine-rich repeats (LRRs), the largest family of plant resistance (R) genes, are vital for plant disease resistance. A comprehensive genome-wide analysis showed that there were only 73 NBS genes in the A. trifoliata genome, including three main subfamilies (50 coiled coil (CC)-NBS-LRR (CNL), 19 Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) and four resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) genes). Additionally, 64 mapped NBS candidates were unevenly distributed on 14 chromosomes, most of which were assigned to the chromosome ends; 41 of these genes were located in clusters, and the remaining 23 genes were singletons. Both the CNLs and TNLs were further divided into four subgroups, and the CNLs had fewer exons than the TNLs. Structurally, all eight previously reported conserved motifs were identified in the NBS domains, and both their order and their amino acid sequences exhibited high conservation. Evolutionarily, tandem and dispersed duplications were shown to be the two main forces responsible for NBS expansion, producing 33 and 29 genes, respectively. A transcriptome analysis of three fruit tissues at four developmental stages showed that NBS genes were generally expressed at low levels, while a few of these genes showed relatively high expression during later development in rind tissues. Overall, this research is the first to identify and characterize A. trifoliata NBS genes and is valuable for both the development of new resistant cultivars and the study of molecular mechanisms of resistance.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

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.

Clubroot in Brassica: recent advances in genomics, breeding, and disease management

Clubroot disease, caused by Plasmodiophora brassicae, affects Brassica oilseed and vegetable production worldwide. This review is focused on various aspects of clubroot disease and its management, including understanding the pathogen and resistance in the host plants. Advances in genetics, molecular biology techniques, and omics research have helped to identify several major loci, QTL, and genes from the Brassica genomes involved in the control of clubroot resistance. Transcriptomic studies have helped to extend our understanding of the mechanism of infection by the pathogen and the molecular basis of resistance/susceptibility in the host plants. A comprehensive understanding of the clubroot disease and host resistance would allow developing a better strategy by integrating the genetic resistance with cultural practices to manage this disease from a long-term perspective.

Adaptive evolution driving the young duplications in six Rosaceae species

Background

In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.

Results

In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.

Conclusions

The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07422-7.

Cotton CC-NBS-LRR Gene GbCNL130 Confers Resistance to Verticillium Wilt Across Different Species

Verticillium wilt (VW) is a destructive disease in cotton caused by Verticillium dahliae and has a significant impact on yield and quality. In the absence of safe and effective chemical control, VW is difficult to manage. Thus, at present, developing resistant varieties is the most economical and effective method of controlling Verticillium wilt of cotton. The CC-NBS-LRR (CNL) gene family is an important class of plant genes involved in disease resistance. This study identified 141 GbCNLs in Gossypium barbadense genome, with 37.5% (53 genes) GbCNLs enriched in 12 gene clusters (GC01-GC12) based on gene distribution in the chromosomes. Especially, seven GbCNLs from two largest clusters (GC11 and GC12) were significantly upregulated in the resistant cultivar (Hai No. 7124) and the susceptible (Giza No. 57). Virus-induced gene silencing of GbCNL130 in G. barbadense, one typical gene in the gene cluster 12 (GC12), significantly altered the response to VW, compromising plant resistance to V. dahliae. In contrast, GbCNL130 overexpression significantly increased the resistance to VW in the wild-type Arabidopsis thaliana. Based on our research findings presented here, we conclude that GbCNL130 promotes resistance to VW by activating the salicylic acid (SA)-dependent defense response pathway resulting in strong accumulation of reactive oxygen species and upregulation of pathogenesis-related (PR) genes. In conclusion, our study resulted in the discovery of a new CNL resistance gene in cotton, GbCNL130, that confers resistance to VW across different hosts.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

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.

A high-quality sponge gourd (Luffa cylindrica) genome

Sponge gourd (Luffa cylindrica) is an important cultivated vegetable and medicinal plant in the family Cucurbitaceae. In this study, a draft genome sequence of the sponge gourd inbred line P93075 was analyzed. Using Illumina, PacBio, and 10× Genomics sequencing techniques as well as new assembly techniques such as FALCON and chromatin interaction mapping (Hi-C), a chromosome-scale genome of approximately 656.19 Mb, with an N50 scaffold length of 48.76 Mb, was generated. From this assembly, 25,508 protein-coding gene loci were identified, and 63.81% of the whole-genome consisted of transposable elements, which are major contributors to the expansion of the sponge gourd genome. According to a phylogenetic analysis of conserved genes, the sponge gourd lineage diverged from the bitter gourd lineage approximately 41.6 million years ago. Additionally, many genes that respond to biotic and abiotic stresses were found to be lineage specific or expanded in the sponge gourd genome, as demonstrated by the presence of 462 NBS-LRR genes, a much greater number than are found in the genomes of other cucurbit species; these results are consistent with the high stress resistance of sponge gourd. Collectively, our study provides insights into genome evolution and serves as a valuable reference for the genetic improvement of sponge gourd.

Genome assembly and characterization of a complex zfBED-NLR gene-containing disease resistance locus in Carolina Gold Select rice with Nanopore sequencing

Long-read sequencing facilitates assembly of complex genomic regions. In plants, loci containing nucleotide-binding, leucine-rich repeat (NLR) disease resistance genes are an important example of such regions. NLR genes constitute one of the largest gene families in plants and are often clustered, evolving via duplication, contraction, and transposition. We recently mapped the Xo1 locus for resistance to bacterial blight and bacterial leaf streak, found in the American heirloom rice variety Carolina Gold Select, to a region that in the Nipponbare reference genome is NLR gene-rich. Here, toward identification of the Xo1 gene, we combined Nanopore and Illumina reads and generated a high-quality Carolina Gold Select genome assembly. We identified 529 complete or partial NLR genes and discovered, relative to Nipponbare, an expansion of NLR genes at the Xo1 locus. One of these has high sequence similarity to the cloned, functionally similar Xa1 gene. Both harbor an integrated zfBED domain, and the repeats within each protein are nearly perfect. Across diverse Oryzeae, we identified two sub-clades of NLR genes with these features, varying in the presence of the zfBED domain and the number of repeats. The Carolina Gold Select genome assembly also uncovered at the Xo1 locus a rice blast resistance gene and a gene encoding a polyphenol oxidase (PPO). PPO activity has been used as a marker for blast resistance at the locus in some varieties; however, the Carolina Gold Select sequence revealed a loss-of-function mutation in the PPO gene that breaks this association. Our results demonstrate that whole genome sequencing combining Nanopore and Illumina reads effectively resolves NLR gene loci. Our identification of an Xo1 candidate is an important step toward mechanistic characterization, including the role(s) of the zfBED domain. Finally, the Carolina Gold Select genome assembly will facilitate identification of other useful traits in this historically important variety.

Plants lack adaptive immunity, and instead contain repeat-rich, disease resistance genes that evolve rapidly through duplication, recombination, and transposition. The number, variation, and often clustered arrangement of these genes make them challenging to sequence and catalog. The US heirloom rice variety Carolina Gold Select has resistance to two important bacterial diseases. Toward identifying the responsible gene(s), we combined long- and short-read sequencing technologies to assemble the whole genome and identify the resistance gene repertoire. We previously narrowed the location of the gene(s) to a region on chromosome four. The region in Carolina Gold Select is larger than in the rice reference genome (Nipponbare) and contains twice as many resistance genes. One shares unusual features with a known bacterial disease resistance gene, suggesting that it confers the resistance. Across diverse varieties and related species, we identified two widely-distributed groups of such genes. The results are an important step toward mechanistic characterization and deployment of the bacterial disease resistance. The genome assembly also identified a resistance gene for a fungal disease and predicted a marker phenotype used in breeding for resistance. Thus, the Carolina Gold Select genome assembly can be expected to aid in the identification and deployment of other valuable traits.

Genome- Wide Analysis of the Nucleotide Binding Site Leucine-Rich Repeat Genes of Four Orchids Revealed Extremely Low Numbers of Disease Resistance Genes

Orchids are one of the most diverse flowering plant families, yet possibly maintain the smallest number of the nucleotide-binding site-leucine-rich repeat (NBS-LRR) type plant resistance (R) genes among the angiosperms. In this study, a genome-wide search in four orchid taxa identified 186 NBS-LRR genes. Furthermore, 214 NBS-LRR genes were identified from seven orchid transcriptomes. A phylogenetic analysis recovered 30 ancestral lineages (29 CNL and one RNL), far fewer than other angiosperm families. From the genetics aspect, the relatively low number of ancestral R genes is unlikely to explain the low number of R genes in orchids alone, as historical gene loss and scarce gene duplication has continuously occurred, which also contributes to the low number of R genes. Due to recent sharp expansions, Phalaenopsis equestris and Dendrobium catenatum having 52 and 115 genes, respectively, and exhibited an “early shrinking to recent expanding” evolutionary pattern, while Gastrodia elata and Apostasia shenzhenica both exhibit a “consistently shrinking” evolutionary pattern and have retained only five and 14 NBS-LRR genes, respectively. RNL genes remain in extremely low numbers with only one or two copies per genome. Notably, all of the orchid RNL genes belong to the ADR1 lineage. A separate lineage, NRG1, was entirely absent and was likely lost in the common ancestor of all monocots. All of the TNL genes were absent as well, coincident with the RNL NRG1 lineage, which supports the previously proposed notion that a potential functional association between the TNL and RNL NRG1 genes.

Animal NLRs continue to inform plant NLR structure and function

Plant NLRs share many of the structural hallmarks of their animal counterparts. At a functional level, the central nucleotide-binding pocket appears to have binding and hydrolysis activities, similar to that of animal NLRs. The TIR domains of plant NLRs have been shown to self-associate, and there is emerging evidence that full-length plant NLRs may do so as well. It is therefore tempting to speculate that plant NLRs may form higher-order complexes similar to those of the mammalian inflammasome. Here we review the available knowledge on structure-function relationships in plant NLRs, focusing on how the information available on animal NLRs informs the mechanism of plant NLR function, and highlight the evidence that innate immunity signalling pathways in multicellular organisms often require the formation of higher-order protein complexes.

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.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

The disease resistance protein SNC1 represses the biogenesis of microRNAs and phased siRNAs

Plants evolved an array of disease resistance genes (R genes) to fight pathogens. In the absence of pathogen infection, NBS-LRR genes, which comprise a major subfamily of R genes, are suppressed by a small RNA cascade involving microRNAs (miRNAs) that trigger the biogenesis of phased siRNAs (phasiRNAs) from R gene transcripts. However, whether or how R genes influence small RNA biogenesis is unknown. In this study, we isolate a mutant with global defects in the biogenesis of miRNAs and phasiRNAs in Arabidopsis thaliana and trace the defects to the over accumulation and nuclear localization of an R protein SNC1. We show that nuclear SNC1 represses the transcription of miRNA and phasiRNA loci, probably through the transcriptional corepressor TPR1. Intriguingly, nuclear SNC1 reduces the accumulation of phasiRNAs from three source R genes and concomitantly, the expression of a majority of the ~170R genes is up-regulated. Taken together, this study suggests an R gene-miRNA-phasiRNA regulatory module that amplifies plant immune responses.

A small RNA-based signaling cascade prevents the induction of plant resistance genes (R-genes) in the absence of pathogen challenge. Here Cai et al. show that nuclear accumulation of the R protein SNC1 can activate immunity by suppressing small RNA production and releasing R-gene repression.

Genome-wide identification, characterization, and evolutionary analysis of NBS-encoding resistance genes in barley

In this study, a systematic analysis of Nucleotide-Binding Site (NBS) disease resistance (R) gene family in the barley, Hordeum vulgare L. cv. Bowman, genome was performed. Using multiple computational analyses, we could identify 96 regular NBS-encoding genes and characterize them on the bases of structural diversity, conserved protein signatures, genomic distribution, gene duplications, differential expression, selection pressure, codon usage, regulation by microRNAs and phylogenetic relationships. Depending on the presence or absence of CC and LRR domains; the identified NBS genes were assigned to four distinct groups; NBS-LRR (53.1%), CC-NBS-LRR (14.6%), NBS (26%), and CC-NBS (6.3%). NBS-associated domain analysis revealed the presence of signal peptides, zinc fingers, diverse kinases, and other structural features. Eighty-five of the identified NBS-encoding genes were mapped onto the seven barley chromosomes, revealing that 50% of them were located on chromosomes 7H, 2H, and 3H, with a tendency of NBS genes to be clustered in the distal telomeric regions of the barley chromosomes. Nine gene clusters, representing 22.35% of total mapped barley NBS-encoding genes, were found, suggesting that tandem duplication stands for an important mechanism in the expansion of this gene family in barley. Phylogenetic analysis determined 31 HvNBS orthologs from rice and Brachypodium. 87 out of 96 HvNBSs were supported by expression evidence, exhibiting various and quantitatively uneven expression patterns across distinct tissues, organs, and development stages. Fourteen potential miRNA-R gene target pairs were further identified, providing insight into the regulation of NBS genes expression. These findings offer candidate target genes to engineer disease-resistant barley genotypes, and promote our understanding of the evolution of NBS-encoding genes in Poaceae crops.

Genome-Wide Identification, Molecular Evolution, and Expression Divergence of Aluminum-Activated Malate Transporters in Apples

Aluminum-activated malate transporters (ALMTs) play an important role in aluminum tolerance, stomatal opening, and fruit acidity in plants. However, the evolutionary pattern of the ALMT gene family in apples remains relatively unknown. In this study, a total of 25 MdALMT genes were identified from the apple reference genome of the "Golden Delicious" doubled-haploid tree (GDDH13). The physiological and biochemical properties, gene structure, and conserved motifs of MdALMT genes were examined. Chromosome location and gene-duplication analysis indicated that whole-genome duplication/segmental duplication played an important role in the expansion of the MdALMT gene family. The Ka/Ks ratio of duplicated MdALMT genes showed that members of this family have undergone strong purifying selection. Through exploration of the phylogenetic relationships, seven subgroups were classified, and higher old gene duplication frequency and significantly different evolutionary rates of the ALMT gene families were detected. In addition, the functional divergence of ALMT genes occurred during the evolutionary process of Rosaceae species. Furthermore, the functional divergence of MdALMT genes was confirmed by expression discrepancy and different subcellular localizations. This study provides the foundation to better understand the molecular evolution of MdALMT genes and further facilitate functional analysis to unravel their exact role in apples.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

The coming of age of EvoMPMI: evolutionary molecular plant-microbe interactions across multiple timescales

Plant-microbe interactions are great model systems to study co-evolutionary dynamics across multiple timescales. However, mechanistic research on plant-microbe interactions has often been conducted with little consideration of evolutionary concepts and methods. Conversely, evolutionary research has rarely integrated the range of mechanisms and models from the molecular plant-microbe interactions field. In recent years, the incipient field of evolutionary molecular plant-microbe interactions (EvoMPMI) has emerged to bridge this gap. Here, we report on some of the recent advances in EvoMPMI. In particular, we highlight new systems to study microbe interactions with early diverging land plants, and new findings from studies of adaptive evolution in pathogens and plants. By linking mechanistic and evolutionary research, EvoMPMI promises to expand our understanding of plant-microbe interactions.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

Evolution, functional differentiation, and co-expression of the RLK gene family revealed in Jilin ginseng, Panax ginseng C.A. Meyer

Most genes in a genome exist in the form of a gene family; therefore, it is necessary to have knowledge of how a gene family functions to comprehensively understand organismal biology. The receptor-like kinase (RLK)-encoding gene family is one of the most important gene families in plants. It plays important roles in biotic and abiotic stress tolerances, and growth and development. However, little is known about the functional differentiation and relationships among the gene members within a gene family in plants. This study has isolated 563 RLK genes (designated as PgRLK genes) expressed in Jilin ginseng (Panax ginseng C.A. Meyer), investigated their evolution, and deciphered their functional diversification and relationships. The PgRLK gene family is highly diverged and formed into eight types. The LRR type is the earliest and most prevalent, while only the Lec type originated after P. ginseng evolved. Furthermore, although the members of the PgRLK gene family all encode receptor-like protein kinases and share conservative domains, they are functionally very diverse, participating in numerous biological processes. The expressions of different members of the PgRLK gene family are extremely variable within a tissue, at a developmental stage and in the same cultivar, but most of the genes tend to express correlatively, forming a co-expression network. These results not only provide a deeper and comprehensive understanding of the evolution, functional differentiation and correlation of a gene family in plants, but also an RLK genic resource useful for enhanced ginseng genetic improvement.

Lineage-specific duplications of NBS-LRR genes occurring before the divergence of six Fragaria species

Background

Plant disease resistance (R) genes are evolving rapidly and play a critical role in the innate immune system of plants. The nucleotide binding sites-leucine rich repeat (NBS-LRR) genes are one of the largest classes in plant R genes. Previous studies have focused on the NBS-LRR genes from one or several species of different genera, and the sequenced genomes of the genus Fragaria offer the opportunity to study the evolutionary processes of these R genes among the closely related species.

Results

In this study, 325, 155, 190, 187, and 133 NBS-LRRs were discovered from F. x ananassa, F. iinumae, F. nipponica, F. nubicola, and F. orientalis, respectively. Together with the 144 NBS-LRR genes from F. vesca, a total of 1134 NBS-LRRs containing 866 multi-genes comprised 184 gene families across the six Fragaria genomes. Extremely short branch lengths and shallow nodes were widely present in the phylogenetic tree constructed with all of the NBS-LRR genes of the six strawberry species. The identities of the orthologous genes were highly significantly greater than those of the paralogous genes, while the Ks ratios of the former were very significantly lower than those of the latter in all of the NBS-LRR gene families. In addition, the Ks and Ka/Ks values of the TIR-NBS-LRR genes (TNLs) were significantly greater than those of the non-TIR-NBS-LRR genes (non-TNLs). Furthermore, the expression patterns of the NBS-LRR genes revealed that the same gene expressed differently under different genetic backgrounds in response to pathogens.

Conclusions

These results, combined with the shared hotspot regions of the duplicated NBS-LRRs on the chromosomes, indicated that the lineage-specific duplication of the NBS-LRR genes occurred before the divergence of the six Fragaria species. The Ks and Ka/Ks ratios suggested that the TNLs are more rapidly evolving and driven by stronger diversifying selective pressures than the non-TNLs.

Electronic supplementary material

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

Characterization and functional assessment of the NLRC3-like molecule of the goldfish (Carassius auratus L.)

The NLRC3-like (NLRC3L) molecule from the goldfish transcriptome database was identified and characterized. Quantitative gene expression analysis revealed the highest mRNA levels of NLRC3L were in the spleen and intestine, with lower mRNA levels observed in muscle and liver. Goldfish NLRC3L was differentially expressed in goldfish immune cell populations with highest mRNA levels measured in PBLs and macrophages. We generated a recombinant form of the molecule (rgfNLRC3L) and an anti-CT-NLRC3L IgG. Treatment of goldfish primary kidney macrophages in vitro with ATP, LPS and heat-killed Aeromonas salmonicida up-regulated the NLRC3L mRNA and protein. Confocal microscopy and co-immunoprecipitation assays indicated that goldfish rgfNLRC3L interacted with apoptosis-associated spec-like protein (ASC) in eukaryotic cells, indicating that NLRC3L may participate in the regulation of the inflammasome responses. The dual-luciferase reporter assay showed that NLRC3L over-expression did not cause the activation of NF-κB, but that it cooperated with RIP2 to down-regulate NF-κB activation. Our results indicate that the NLRC3L may function as a regulator of NLR pathways in teleosts.

Diversification of defensins and NLRs in Arabidopsis species by different evolutionary mechanisms

BACKGROUND

Genes encoding proteins underlying host-pathogen co-evolution and which are selected for new resistance specificities frequently are under positive selection, a process that maintains diversity. Here, we tested the contribution of natural selection, recombination and transcriptional divergence to the evolutionary diversification of the plant defensins superfamily in three Arabidopsis species. The intracellular NOD-like receptor (NLR) family was used for comparison because positive selection has been well documented in its members. Similar to defensins, NLRs are encoded by a large and polymorphic gene family and many of their members are involved in the immune response.

RESULTS

Gene trees of Arabidopsis defensins (DEFLs) show a high prevalence of clades containing orthologs. This indicates that their diversity dates back to a common ancestor and species-specific duplications did not significantly contribute to gene family expansion. DEFLs are characterized by a pervasive pattern of neutral evolution with infrequent positive and negative selection as well as recombination. In comparison, most NLR alignment groups are characterized by frequent occurrence of positive selection and recombination in their leucine-rich repeat (LRR) domain as well negative selection in their nucleotide-binding (NB-ARC) domain. While major NLR subgroups are expressed in pistils and leaves both in presence or absence of pathogen infection, the members of DEFL alignment groups are predominantly transcribed in pistils. Furthermore, conserved groups of NLRs and DEFLs are differentially expressed in response to Fusarium graminearum regardless of whether these genes are under positive selection or not.

CONCLUSIONS

The present analyses of NLRs expands previous studies in Arabidopsis thaliana and highlights contrasting patterns of purifying and diversifying selection affecting different gene regions. DEFL genes show a different evolutionary trend, with fewer recombination events and significantly fewer instances of natural selection. Their heterogeneous expression pattern suggests that transcriptional divergence probably made the major contribution to functional diversification. In comparison to smaller families encoding pathogenesis-related (PR) proteins under positive selection, DEFLs are involved in a wide variety of processes that altogether might pose structural and functional trade-offs to their family-wide pattern of evolution.

In situ conservation-harnessing natural and human-derived evolutionary forces to ensure future crop adaptation

Ensuring the availability of the broadest possible germplasm base for agriculture in the face of increasingly uncertain and variable patterns of biotic and abiotic change is fundamental for the world's future food supply. While ex situ conservation plays a major role in the conservation and availability of crop germplasm, it may be insufficient to ensure this. In situ conservation aims to maintain target species and the collective genotypes they represent under evolution. A major rationale for this view is based on the likelihood that continued exposure to changing selective forces will generate and favor new genetic variation and an increased likelihood that rare alleles that may be of value to future agriculture are maintained. However, the evidence that underpins this key rationale remains fragmented and has not been examined systematically, thereby decreasing the perceived value and support for in situ conservation for agriculture and food systems and limiting the conservation options available. This study reviews evidence regarding the likelihood and rate of evolutionary change in both biotic and abiotic traits for crops and their wild relatives, placing these processes in a realistic context in which smallholder farming operates and crop wild relatives continue to exist. It identifies areas of research that would contribute to a deeper understanding of these processes as the basis for making them more useful for future crop adaptation.

Promoter variants of Xa23 alleles affect bacterial blight resistance and evolutionary pattern

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is the most important bacterial disease in rice (Oryza sativa L.). Our previous studies have revealed that the bacterial blight resistance gene Xa23 from wild rice O. rufipogon Griff. confers the broadest-spectrum resistance against all the naturally occurring Xoo races. As a novel executor R gene, Xa23 is transcriptionally activated by the bacterial avirulence (Avr) protein AvrXa23 via binding to a 28-bp DNA element (EBE_AvrXa23) in the promoter region. So far, the evolutionary mechanism of Xa23 remains to be illustrated. Here, a rice germplasm collection of 97 accessions, including 29 rice cultivars (indica and japonica) and 68 wild relatives, was used to analyze the evolution, phylogeographic relationship and association of Xa23 alleles with bacterial blight resistance. All the ~ 473 bp DNA fragments consisting of promoter and coding regions of Xa23 alleles in the germplasm accessions were PCR-amplified and sequenced, and nine single nucleotide polymorphisms (SNPs) were detected in the promoter regions (~131 bp sequence upstream from the start codon ATG) of Xa23/xa23 alleles while only two SNPs were found in the coding regions. The SNPs in the promoter regions formed 5 haplotypes (Pro-A, B, C, D, E) which showed no significant difference in geographic distribution among these 97 rice accessions. However, haplotype association analysis indicated that Pro-A is the most favored haplotype for bacterial blight resistance. Moreover, SNP changes among the 5 haplotypes mostly located in the EBE/ebe regions (EBE_AvrXa23 and corresponding ebes located in promoters of xa23 alleles), confirming that the EBE region is the key factor to confer bacterial blight resistance by altering gene expression. Polymorphism analysis and neutral test implied that Xa23 had undergone a bottleneck effect, and selection process of Xa23 was not detected in cultivated rice. In addition, the Xa23 coding region was found highly conserved in the Oryza genus but absent in other plant species by searching the plant database, suggesting that Xa23 originated along with the diversification of the Oryza genus from the grass family during evolution. This research offers a potential for flexible use of novel Xa23 alleles in rice breeding programs and provide a model for evolution analysis of other executor R genes.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Genome reconstruction in Cynara cardunculus taxa gains access to chromosome-scale DNA variation

The genome sequence of globe artichoke (Cynara cardunculus L. var. scolymus, 2n = 2x = 34) is now available for use. A survey of C. cardunculus genetic resources is essential for understanding the evolution of the species, carrying out genetic studies and for application of breeding strategies. We report on the resequencing analyses (~35×) of four globe artichoke genotypes, representative of the core varietal types, as well as a genotype of the related taxa cultivated cardoon. The genomes were reconstructed at a chromosomal scale and structurally/functionally annotated. Gene prediction indicated a similar number of genes, while distinctive variations in miRNAs and resistance gene analogues (RGAs) were detected. Overall, 23,5 M SNP/indel were discovered (range 6,34 M –14,50 M). The impact of some missense SNPs on the biological functions of genes involved in the biosynthesis of phenylpropanoid and sesquiterpene lactone secondary metabolites was predicted. The identified variants contribute to infer on globe artichoke domestication of the different varietal types, and represent key tools for dissecting the path from sequence variation to phenotype. The new genomic sequences are fully searchable through independent Jbrowse interfaces (www.artichokegenome.unito.it), which allow the analysis of collinearity and the discovery of genomic variants, thus representing a one-stop resource for C. cardunculus genomics.

Distinct Patterns of Gene Gain and Loss: Diverse Evolutionary Modes of NBS-Encoding Genes in Three Solanaceae Crop Species

Plant resistance conferred by nucleotide binding site (NBS)-encoding resistance genes plays a key role in the defense against various pathogens throughout the entire plant life cycle. However, comparative analyses for the systematic evaluation and determination of the evolutionary modes of NBS-encoding genes among Solanaceae species are rare. In this study, 447, 255, and 306 NBS-encoding genes were identified from the genomes of potato, tomato, and pepper, respectively. These genes usually clustered as tandem arrays on chromosomes; few existed as singletons. Phylogenetic analysis indicated that three subclasses [TNLs (TIR-NBS-LRR), CNLs (CC-NBS-LRR), and RNLs (RPW8-NBS-LRR)] each formed a monophyletic clade and were distinguished by unique exon/intron structures and amino acid motif sequences. By comparing phylogenetic and systematic relationships, we inferred that the NBS-encoding genes in the present genomes of potato, tomato, and pepper were derived from 150 CNL, 22 TNL, and 4 RNL ancestral genes, and underwent independent gene loss and duplication events after speciation. The NBS-encoding genes therefore exhibit diverse and dynamic evolutionary patterns in the three Solanaceae species, giving rise to the discrepant gene numbers observed today. Potato shows a “consistent expansion” pattern, tomato exhibits a pattern of “first expansion and then contraction,” and pepper presents a “shrinking” pattern. The earlier expansion of CNLs in the common ancestor led to the dominance of this subclass in gene numbers. However, RNLs remained at low copy numbers due to their specific functions. Along the evolutionary process of NBS-encoding genes in Solanaceae, species-specific tandem duplications contributed the most to gene expansions.

Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species

Background

Nucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.

Results

Based on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.

Conclusions

The structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3682-x) contains supplementary material, which is available to authorized users.

Quantitative Trait Loci and Resistance Gene Analogs Associated with Maize White Spot Resistance

Maize white spot (MWS), caused by the bacterium Pantoea ananatis, is one of the most important maize foliar diseases in tropical and subtropical regions, causing significant yield losses. Despite its economic importance, genetic studies of MWS are scarce. The aim of this study was to map quantitative trait loci (QTL) associated with MWS resistance and to identify resistance gene analogs (RGA) underlying these QTL. QTL mapping was performed in a tropical maize F_2:3 population, which was genotyped with simple-sequence repeat and RGA-tagged markers and phenotyped for the response to MWS in two Brazilian southeastern locations. Nine QTL explained approximately 70% of the phenotypic variance for MWS resistance at each location, with two of them consistently detected in both environments. Data mining using 112 resistance genes cloned from different plant species revealed 1,697 RGA distributed in clusters within the maize genome. The RGA Pto19, Pto20, Pto99, and Xa26.151.4 were genetically mapped within MWS resistance QTL on chromosomes 4 and 8 and were preferentially expressed in the resistant parental line at locations where their respective QTL occurred. The consistency of QTL mapping, in silico prediction, and gene expression analyses revealed RGA and genomic regions suitable for marker-assisted selection to improve MWS resistance.