Identification and Characterization of NBS Resistance Genes in Akebia trifoliata

Arogenate dehydratase isoforms strategically deregulate phenylalanine biosynthesis in Akebia trifoliata

Arogenate dehydratase (ADT) is key for phenylalanine (Phe) biosynthesis in plants. To examine ADT components and function in Akebia trifoliata, a representative of Ranunculaceae, we first identified eight ADTs (AktADT1-8, encoding sequences varying from 1032 to 1962 bp) in the A. trifoliata reference genome and five proteins (AktADT1, AktADT4, AktADT7, AktADT8 and AktADT8s) with moonlighting prephenate dehydratase (PDT) activity and Km values varying from 0.43 to 2.17 mM. Structurally, two basic residue combinations (Val314/Ala317 and Ala314/Val317) in the PAC domain are essential for the moonlighting PDT activity of ADTs. Functionally, AktADT4 and AktADT8 successfully restored the wild-type phenotype of pha2, a knockout mutant of Saccharomyces cerevisiae. In addition, AktADTs are ubiquitously expressed, but their expression levels are tissue specific, and the half maximal inhibitory concentration (IC50) of Phe for AktADTs ranged from 49.81 to 331.17 μM. Both AktADT4 and AktADT8 and AktADT8s localized to chloroplast stromules and the cytosol, respectively, while the remaining AktADTs localized to the chloroplast stroma. These findings suggest that various strategies exist for regulating Phe biosynthesis in A. trifoliata. This provides a reasonable explanation for the high Phe content and insights for further genetic improvement of the edible fruits of A. trifoliata.

Genome-Wide Analysis of the Polygalacturonase Gene Family Sheds Light on the Characteristics, Evolutionary History, and Putative Function of Akebia trifoliata

Polygalacturonase (PG) is one of the largest families of hydrolytic enzymes in plants. It is involved in the breakdown of pectin in the plant cell wall and even contributes to peel cracks. Here, we characterize PGs and outline their expression profiles using the available reference genome and transcriptome of Akebia trifoliata. The average length and exon number of the 47 identified AktPGs, unevenly assigned on 14 chromosomes and two unassembled contigs, were 5399 bp and 7, respectively. The phylogenetic tree of 191 PGs, including 47, 57, 51, and 36 from A. trifoliata, Durio zibethinus, Actinidia chinensis, and Vitis vinifera, respectively, showed that AktPGs were distributed in all groups except group G and that 10 AktPGs in group E were older, while the remaining 37 AktPGs were younger. Evolutionarily, all AktPGs generally experienced whole-genome duplication (WGD)/segmental repeats and purifying selection. Additionally, the origin of conserved domain III was possibly associated with a histidine residue (H) substitute in motif 8. The results of both the phylogenetic tree and expression profiling indicated that five AktPGs, especially AktPG25, could be associated with the cracking process. Detailed information and data on the PG family are beneficial for further study of the postharvest biology of A. trifoliata.

Germplasm resources and genetic improvement of Akebia: A new fruit crop in China

Akebia species, belonging to Lardizabalaceae, are widespread from subtropical to temperate environments of China, Japan, and Korea. All known Akebia species have medicinal and dietary value and have been widely cultivated as a new fruit crop in many areas of China. However, compared with other crop species, the breeding improvement and commercial cultivation of Akebia remain in their infancy. This review systematically introduces the present germplasm resources, geographical distribution, biological characteristics, interspecific and intraspecific cross compatibility, molecular biology, and breeding progress in Akebia species. Akebia plants are widely distributed in Shanxi, Henan, Sichuan, Chongqing, Hunan, Hubei, Jiangxi, Zhejiang, and Fujian provinces of China, and wild Akebia plants exhibit abundant phenotypic and genetic diversity due to their wide range of geographical distribution and high adaptability in different habitats. Interspecific artificial hybridization experiments have been conducted in our Akebia germplasm resources nursery. The results showed that there was no reproductive isolation between Akebia species, and fertile progeny could be produced. The synthesis of knowledge on these species provides insights for the rational development and utilization of these germplasm resources, and can facilitate the development of new breeding lines or varieties for commercial cultivation or production. Finally, perspectives on Akebia breeding research are discussed and conclusions are provided. This review provided breeders with new insights into Akebia domestication and breeding, and we also proposed five basic steps in the domestication of new fruit crops.

Phenotypic characterization of the wheat temperature-sensitive leaf color mutant and physical mapping of mutant gene by reduced-representation sequencing

Few available leaf color mutants in crops have greatly limited the understanding of photosynthesis mechanisms, leading to few accomplishments in crop yield improvement via enhanced photosynthetic efficiency. Here, a noticeable albino mutant, CN19M06, was identified. A comparison between CN19M06 and the wild type CN19 at different temperatures showed that the albino mutant was temperature-sensitive and produced leaves with a decreased chlorophyll content at temperatures below 10 °C. Genetic analysis suggested that the albinism was controlled by one recessive nuclear gene named TSCA1, which was putatively assigned to the region of 718.1-729.8 Mb on chromosome 2AL using bulked-segregant analysis and double-digest restriction site-associated DNA. Finally, molecular linkage analysis physically anchored TSCA1 to a narrowed region of 718.8-725.3 Mb with a 6.5 Mb length on 2AL flanked by InDel 18 and InDel 25 with 0.7 cM genetic interval. Among the 111 annotated functional genes in the corresponding chromosomal region, only TraesCS2A01G487900 of the PAP fibrillin family was both related to chlorophyll metabolism and temperature sensitivity; therefore, it was considered the putative candidate gene of TSCA1. Overall, CN19M06 has great potential for exploring the molecular mechanism of photosynthesis and monitoring temperature changes in wheat production.

Genome-Wide Identification of Superoxide Dismutase and Expression in Response to Fruit Development and Biological Stress in Akebia trifoliata: A Bioinformatics Study

Akebia trifoliata is a newly domesticated perennial fruit tree, and the lack of molecular research on stress resistance seriously affects its genetic improvement and commercial value development. Superoxide dismutase (SOD) can effectively eliminate the accumulation of reactive oxygen species (ROS) during the rapid growth of plant organs under biotic and abiotic stresses, maintaining a steady state of physiological metabolism. In this study, 13 SODs consisting of two FeSODs (FSDs), four MnSODs (MSDs) and seven Cu/ZnSODs (CSDs) were identified in the A. trifoliata genome. Structurally, the phylogeny, intron–exon pattern and motif sequences within these three subfamilies show high conservation. Evolutionarily, segmental/wide genome duplication (WGD) and dispersed duplication form the current SOD profile of A. trifoliata. Weighted gene coexpression network analysis (WGCNA) revealed the metabolic pathways of nine (69.2%) SODs involved in fruit development, among which AktMSD4 regulates fruit development and AktCSD4 participates in the stress response. In addition, under the stress of multiple pathogens, six (46.6%) SODs were continuously upregulated in the rinds of resistant lines; of these, three SODs (AktMSD1, AktMSD2 and AktMSD3) were weakly or not expressed in susceptible lines. The results pave the way for theoretical research on SODs and afford the opportunity for genetic improvement of A. trifoliata.

Genome-wide identification, characterization, and expression profile ofNBS-LRRgene family in sweet orange (Citrussinensis)

BACKGROUND

The NBS-LRR (nucleotide-binding site-leucine-rich repeat gene) gene family, known as the plant R (resistance) gene family with the most members, plays a significant role in plant resistance to various external adversity stresses. The NBS-LRR gene family has been researched in many plant species. Citrus is one of the most vital global cash crops, the number one fruit group, and the third most traded agricultural product world wild. However, as one of the largest citrus species, a comprehensive study of the NBS-LRR gene family has not been reported on sweet oranges.

METHODS

In this study, NBS-LRR genes were identified from the Citrus sinensis genome (v3.0), with a comprehensive analysis of this gene family performed, including phylogenetic analysis, gene structure, cis-acting element of a promoter, and chromosomal localization, among others. The expression pattern of NBS-LRR genes was analyzed when sweet orange fruits were infected by Penicillium digitatum, employing experimental data from our research group. It first reported the expression patterns of NBS-LRR genes under abiotic stresses, using three transcript data from NCBI (National Center for Biotechnology Information).

RESULTS

In this study, 111 NBS-LRR genes were identified in the C. sinensis genome (v3.0) and classified into seven subfamilies according to their N-terminal and C-terminal domains. The phylogenetic tree results indicate that genes containing only the NBS structural domain are more ancient in the sweet orange NBS-LRR gene family. The chromosome localization results showed that 111 NBS-LRR genes were distributed unevenly on nine chromosomes, with the most genes distributed on chromosome 1. In addition, we identified a total of 18 tandem duplication gene pairs in the sweet orange NBS-LRR gene family, and based on the Ka/Ks ratio, all of the tandem duplication genes underwent purifying selection. Transcriptome data analysis showed a significant number of NBS-LRR genes expressed under biotic and abiotic stresses, and some reached significantly different levels of expression. It indicates that the NBS-LRR gene family is vital in resistance to biotic and abiotic stresses in sweet oranges.

CONCLUSION

Our study provides the first comprehensive framework on the NBS-LRR family of genes, which provides a basis for further in-depth studies on the biological functions of NBS-LRR in growth, development, and response to abiotic stresses in sweet orange.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

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.

The chromosome-level genome of Akebia trifoliata as an important resource to study plant evolution and environmental adaptation in the Cretaceous

The environmental adaptation of eudicots is the most reasonable explanation for why they compose the largest clade of modern plants (>70% of angiosperms), which indicates that the basal eudicots would be valuable and helpful to study their survival and ability to thrive throughout evolutionary processes. Here, we detected two whole-genome duplication (WGD) events in the high-quality assembled Akebia trifoliata genome (652.73 Mb) with 24 138 protein-coding genes based on the evidence of intragenomic and intergenomic collinearity, synonymous substitution rate (K_S ) values and polyploidization and diploidization traces; these events putatively occurred at 85.15 and 146.43 million years ago (Mya). The integrated analysis of 16 species consisting of eight basal and eight core eudicots further revealed that there was a putative ancient WGD at the early stage of eudicots (temporarily designated θ) at 142.72 Mya, similar to the older WGD of Akebia trifoliata, and a putative core eudicot-specific WGD (temporarily designated ω). Functional enrichment analysis of retained duplicate genes following the θ event is suggestive of adaptation to the extreme environment change in both the carbon dioxide concentration and desiccation around the Jurassic-Cretaceous boundary, while the retained duplicate genes following the ω event is suggestive of adaptation to the extreme droughts, possibly leading to the rapid spread of eudicots in the mid-Cretaceous. Collectively, the A. trifoliata genome experienced two WGD events, and the older event may have occurred at the early stage of eudicots, which likely increased plant environmental adaptability and helped them survive in ancient extreme environments.

Evolution patterns of NBS genes in the genus Dendrobium and NBS-LRR gene expression in D. officinale by salicylic acid treatment

BACKGROUND

Dendrobium officinale Kimura et Migo, which contains rich polysaccharides, flavonoids and alkaloids, is a Traditional Chinese Medicine (TCM) with important economic benefits, while various pathogens have brought huge losses to its industrialization. NBS gene family is the largest class of plant disease resistance (R) genes, proteins of which are widely distributed in the upstream and downstream of the plant immune systems and are responsible for receiving infection signals and regulating gene expression respectively. It is of great significance for the subsequent disease resistance breeding of D. officinale to identify NBS genes by using the newly published high-quality chromosome-level D. officinale genome.

RESULTS

In this study, a total of 655 NBS genes were uncovered from the genomes of D. officinale, D. nobile, D. chrysotoxum, V. planifolia, A. shenzhenica, P. equestris and A. thaliana. The phylogenetic results of CNL-type protein sequences showed that orchid NBS-LRR genes have significantly degenerated on branches a and b. The Dendrobium NBS gene homology analysis showed that the Dendrobium NBS genes have two obvious characteristics: type changing and NB-ARC domain degeneration. Because the NBS-LRR genes have both NB-ARC and LRR domains, 22 D. officinale NBS-LRR genes were used for subsequent analyses, such as gene structures, conserved motifs, cis-elements and functional annotation analyses. All these results suggested that D. officinale NBS-LRR genes take part in the ETI system, plant hormone signal transduction pathway and Ras signaling pathway. Finally, there were 1,677 DEGs identified from the salicylic acid (SA) treatment transcriptome data of D. officinale. Among them, six NBS-LRR genes (Dof013264, Dof020566, Dof019188, Dof019191, Dof020138 and Dof020707) were significantly up-regulated. However, only Dof020138 was closely related to other pathways from the results of WGCNA, such as pathogen identification pathways, MAPK signaling pathways, plant hormone signal transduction pathways, biosynthetic pathways and energy metabolism pathways.

CONCLUSION

Our results revealed that the NBS gene degenerations are common in the genus Dendrobium, which is the main reason for the diversity of NBS genes, and the NBS-LRR genes generally take part in D. officinale ETI system and signal transduction pathways. In addition, the D. officinale NBS-LRR gene Dof020138, which may have an important breeding value, is indirectly activated by SA in the ETI system.

Ectopic Expression of HbRPW8-a from Hevea brasiliensis Improves Arabidopsis thaliana Resistance to Powdery Mildew Fungi (Erysiphe cichoracearum UCSC1)

The RPW8s (Resistance to Powdery Mildew 8) are atypical broad-spectrum resistance genes that provide resistance to the powdery mildew fungi. Powdery mildew of rubber tree is one of the serious fungal diseases that affect tree growth and latex production. However, the RPW8 homologs in rubber tree and their role of resistance to powdery mildew remain unclear. In this study, four RPW8 genes, HbRPW8-a, b, c, d, were identified in rubber tree, and phylogenetic analysis showed that HbRPW8-a was clustered with AtRPW8.1 and AtRPW8.2 of Arabidopsis. The HbRPW8-a protein was localized on the plasma membrane and its expression in rubber tree was significantly induced upon powdery mildew infection. Transient expression of HbRPW8-a in tobacco leaves induced plant immune responses, including the accumulation of reactive oxygen species and the deposition of callose in plant cells, which was similar to that induced by AtRPW8.2. Consistently, overexpression of HbRPW8-a in Arabidopsis thaliana enhanced plant resistance to Erysiphe cichoracearum UCSC1 and Pseudomonas syringae pv. tomato DC30000 (PstDC3000). Moreover, such HbRPW8-a mediated resistance to powdery mildew was in a salicylic acid (SA) dependent manner. Taken together, we demonstrated a new RPW8 member in rubber tree, HbRPW8-a, which could potentially contribute the resistance to powdery mildew.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Akebiatrifoliata: A Bioinformatics Study

WRKY transcription factors have been found in most plants and play an important role in regulating organ growth and disease response. Outlining the profile of WRKY genes is a very useful project for studying morphogenesis and resistance formation. In the present study, a total of 63 WRKY genes consisting of 13 class I, 41 class II, and 9 class III genes were identified from the newly published A. trifoliata genome, of which 62 were physically distributed on all 16 chromosomes. Structurally, two AkWRKY genes (AkWRKY6 and AkWRKY52) contained four domains, and AkWRKY17 lacked the typical heptapeptide structure. Evolutionarily, 42, 16, and 5 AkWRKY genes experienced whole genome duplication (WGD) or fragmentation, dispersed duplication, and tandem duplication, respectively; 28 Ka/Ks values of 30 pairs of homologous genes were far lower than 1, while those of orthologous gene pairs between AkWRKY41 and AkWRKY52 reached up to 2.07. Transcriptome analysis showed that many of the genes were generally expressed at a low level in 12 fruit samples consisting of three tissues, including rind, flesh, and seeds, at four developmental stages, and interaction analysis between AkWRKY and AkNBS genes containing W-boxes suggested that AkWRKY24 could play a role in plant disease resistance by positively regulating AkNBS18. In summary, the WRKY gene family of A. trifoliata was systemically characterized for the first time, and the data and information obtained regarding AkWRKY could be very useful in further theoretically elucidating the molecular mechanisms of plant development and response to pathogens and practically improving favorable traits such as disease resistance.

Expression Profiles of Microsatellites in Fruit Tissues of Akebia trifoliata and Development of Efficient EST-SSR Markers

Akebia trifoliata, a member of the family Lardizabalaceae, has high exploitation potential for multiple economic purposes, so genetic improvements to meet requirements for commercial demand are needed. However, this progress is largely impeded by a lack of effective selection markers. In this study, we obtained 271.49 Gb of clean transcriptomic data from 12 samples (three tissues at four developmental stages) of A. trifoliata fruit. We identified 175,604, 194,370, and 207,906 SSRs from the de novo assembled 416,363, 463,756, and 491,680 unigene sequences obtained from the flesh, seed, and rind tissues, respectively. The profile and proportion of SSR motifs expressed in each fruit tissue and developmental stage were remarkably similar, but many trinucleotide repeats had differential expression levels among different tissues or at different developmental stages. In addition, we successfully designed 16,869 functional EST-SSR primers according to the annotated unigenes. Finally, 94 and 72 primer pairs out of 100 randomly selected primer pairs produced clear bands and polymorphic bands, respectively. These results were also used to elucidate the expression profiles of different tissues at various stages. Additionally, we provided a set of effective, polymorphic, and reliable EST-SSR markers sufficient for accelerating the discovery of metabolic and pathway-specific functional genes for genetic improvement and increased commercial productivity.

Characterization of Microsatellites in the Akebia trifoliata Genome and Their Transferability and Development of a Whole Set of Effective, Polymorphic, and Physically Mapped Simple Sequence Repeat Markers

Akebia trifoliata is a perennial climbing woody liana plant with a high potential for commercial exploitation and theoretical research. Similarly, microsatellites (simple sequence repeats, SSRs) also have dual roles: as critical markers and as essential elements of the eukaryotic genome. To characterize the profile of SSRs and develop molecular markers, the high-quality assembled genome of A. trifoliata was used. Additionally, to determine the potential transferability of SSR loci, the genomes of Amborella trichopoda, Oryza sativa, Vitis vinifera, Arabidopsis thaliana, Papaver somniferum, and Aquilegia coerulea were also used. We identified 434,293 SSRs with abundant short repeats, such as 290,868 (66.98%) single-nucleotide repeats (SNRs) and 113,299 (26.09%) dinucleotide repeats (DNRs) in the A. trifoliata genome. 398,728 (91.81%) SSRs on 344,283 loci were physically mapped on the chromosomes, and a positive correlation (r = 0.98) was found between the number of SSRs and chromosomal length. Additionally, 342,916 (99.60%) potential SSR markers could be designed from the 344,283 physically mapped loci, while only 36,160 could be viewed as high-polymorphism-potential (HPP) markers, findings that were validated by PCR. Finally, SSR loci exhibited broad potential transferability, particularly DNRs such as the "AT/AT" and "AG/CT" loci, among all angiosperms, a finding that was not related to the genetic divergence distance. Practically, we developed a whole set of effective, polymorphic, and physically anchored molecular markers and found that, evolutionarily, DNRs could be responsible for microsatellite origin and protecting gene function.

Assessment of the Breeding Potential of a Set of Genotypes Selected from a Natural Population of Akebia trifoliata (Three–Leaf Akebia)

Akebia trifoliata (three-leaf akebia) has long been used as a medicinal herb and has the potential to be used in diverse ways, especially as a fruit crop. However, efforts to domesticate and cultivate new varieties for commercial use are only in their infancy. Here, we evaluated the genetic diversity of 29 genotypes, which were previously selected from a natural population consisting of 1447 genotypes and exhibiting high resistance to fungal diseases and a smooth peel of A. trifoliata using 85 genome-specific single sequence repeat (SSR) markers. We also characterized variation in 19 phenotypic traits and nutritional components. Large variation in phenotypic traits and nutritional components was observed, especially in vitamin C, seed/pulp, and fruit color. Correlation analyses revealed that many phenotypic traits and nutritional components were significantly correlated. A principal component analysis identified five principal components, which explained 83.2% of the total variation in the data. The results of the SSR analysis revealed that 80 of the 85 SSR markers were polymorphic; the total number of alleles amplified was 532. The expected heterozygosity was 0.672, and Shannon’s information index was 1.328. A Ward dendrogram and unweighted pair group method with arithmetic mean dendrogram revealed high diversity among the 29 genotypes and suggested that the measured morphological and nutritional traits were genetically independent of disease resistance and texture traits, as well as SSR marker loci. Finally, our results suggest that additional rounds of selection from the selected population, despite its small size, could be effective for the development of new A. trifoliata fruit cultivars.

NLRexpress-A bundle of machine learning motif predictors-Reveals motif stability underlying plant Nod-like receptors diversity

Examination of a collection of over 80,000 Plant Nod-like receptors (NLRs) revealed an overwhelming sequence diversity underlying functional specificity of pathogen detection, signaling and cooperativity. The NLR canonical building blocks-CC/TIR/RPW8, NBS and LRR-contain, however, a number of conserved sequence motifs showing a significant degree of invariance amongst different NLR groups. To identify these motifs we developed NLRexpress-a bundle of 17 machine learning (ML)-based predictors, able to swiftly and precisely detect CC, TIR, NBS, and LRR motifs while minimizing computing time without accuracy losses-aimed as an instrument scalable for screening overall proteomes, transcriptomes or genomes for identifying integral NLRs and discriminating them against incomplete sequences lacking key motifs. These predictors were further used to screen a subset of ∼34,000 regular plant NLR sequences. Motifs were analyzed using unsupervised ML techniques to assess the structural correlations hidden underneath pattern variabilities. Both the NB-ARC switch domain which admittedly is the most conserved region of NLRs and the highly diverse LRR domain with its vastly variable lengths and repeat irregularities-show well-defined relations between motif subclasses, highlighting the importance of structural invariance in shaping NLR sequence diversity. The online NLRexpress webserver can be accessed at https://nlrexpress.biochim.ro.