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

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

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

Functional characterization of NBS-LRR genes reveals an NBS-LRR gene that mediates resistance against Fusarium wilt

BACKGROUND

Most disease resistance (R) genes in plants encode proteins that contain leucine-rich-repeat (LRR) and nucleotide-binding site (NBS) domains, which belong to the NBS-LRR family. The sequenced genomes of Fusarium wilt-susceptible Vernicia fordii and its resistant counterpart, Vernicia montana, offer significant resources for the functional characterization and discovery of novel NBS-LRR genes in tung tree.

RESULTS

Here, we identified 239 NBS-LRR genes across two tung tree genomes: 90 in V. fordii and 149 in V. montana. Five VmNBS-LRR paralogous were predicted in V. montana, and 43 orthologous were detected between V. fordii and V. montana. The orthologous gene pair Vf11G0978-Vm019719 exhibited distinct expression patterns in V. fordii and V. montana: Vf11G0978 showed downregulated expression in V. fordii, while its orthologous gene Vm019719 demonstrated upregulated expression in V. montana, indicating that this pair may be responsible for the resistance to Fusarium wilt in V. montana. Vm019719 from V. montana, activated by VmWRKY64, was shown to confer resistance to Fusarium wilt in V. montana by a virus-induced gene silencing (VIGS) experiment. However, in the susceptible V. fordii, its allelic counterpart, Vf11G0978, exhibited an ineffective defense response, attributed to a deletion in the promoter's W-box element.

CONCLUSIONS

This study provides the first systematic analysis of NBS-LRR genes in the tung tree and identifies a candidate gene that can be utilized for marker-assisted breeding to control Fusarium wilt in V. fordii.

Deciphering peanut complex genomes paves a way to understand its origin and domestication

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Genome wide identification and evolutionary analysis of vat like NBS-LRR genes potentially associated with resistance to aphids in cotton

Nucleotide Binding Site - Leucine Rich Repeat (NBS-LRR) genes play a significant role in plant defense against biotic stresses and are an integral part of signal transduction pathways. Vat gene has been well reported for their role in resistance to Aphis gossypii and viruses transmitted by them. Despite their importance, Vat like NBS-LRR resistance genes have not yet been identified and studied in cotton species. This study report hundreds of orthologous Vat like NBS-LRR genes from the genomes of 18 cotton species through homology searches and the distribution of those identified genes were tend to be clustered on different chromosome. Especially, in a majority of the cases, Vat like genes were located on chromosome number 13 and they all shared two conserved NBS-LRR domains, one disease resistant domain and several repeats of LRR on the investigated cotton Vat like proteins. Gene ontology study on Vat like NBS-LRR genes revealed the molecular functions viz., ADP and protein binding. Phylogenetic analysis also revealed that Vat like sequences of two diploid species, viz., G. arboreum and G. anomalum, were closely related to the sequences of the tetraploids than all other diploids. The Vat like genes of G. aridum and G. schwendimanii were distantly related among diploids and tetraploids species. Various hormones and defense related cis-acting regulatory elements were identified from the 2 kb upstream sequences of the Vat like genes implying their defensive response towards the biotic stresses. Interestingly, G. arboreum and G. trilobum were found to have more regulatory elements than larger genomes of tetraploid cotton species. Thus, the present study provides the evidence for the evolution of Vat like genes in defense mechanisms against aphids infestation in cotton genomes and allows further characterization of candidate genes for developing aphid and aphid transmitted viruses resistant crops through cotton breeding.

Deciphering evolutionary dynamics of WRKY genes in Arachis species

BACKGROUND

Cultivated peanut (Arachis hypogaea), a progeny of the cross between A. duranensis and A. ipaensis, is an important oil and protein crop from South America. To date, at least six Arachis genomes have been sequenced. WRKY transcription factors (TFs) play crucial roles in plant growth, development, and response to abiotic and biotic stresses. WRKY TFs have been identified in A. duranensis, A. ipaensis, and A. hypogaea cv. Tifrunner; however, variations in their number and evolutionary patterns across various Arachis spp. remain unclear.

RESULTS

WRKY TFs were identified and compared across different Arachis species, including A. duranensis, A. ipaensis, A. monticola, A. hypogaea cultivars (cv.) Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The results showed that the WRKY TFs underwent dynamic equilibrium between diploid and tetraploid peanut species, characterized by the loss of old WRKY TFs and retention of the new ones. Notably, cultivated peanuts inherited more conserved WRKY orthologs from wild tetraploid peanuts than their wild diploid donors. Analysis of the W-box elements and protein-protein interactions revealed that different domestication processes affected WRKY evolution across cultivated peanut varieties. WRKY TFs of A. hypogaea cv. Fuhuasheng and Shitouqi exhibited a similar domestication process, while those of cv. Tifrunner of the same species underwent a different domestication process based on protein-protein interaction analysis.

CONCLUSIONS

This study provides new insights into the evolution of WRKY TFs in Arachis spp.

Genome-wide identification, evolutionary and expression analysis of the cyclin-dependent kinase gene family in peanut

BACKGROUND

Cyclin-dependent kinases (CDKs) are a predominant group of serine/threonine protein kinases that have multi-faceted functions in eukaryotes. The plant CDK members have well-known roles in cell cycle progression, transcriptional regulation, DNA repair, abiotic stress and defense responses, making them promising targets for developing stress adaptable high-yielding crops. There is relatively sparse information available on the CDK family genes of cultivated oilseed crop peanut and its diploid progenitors.

RESULTS

We have identified 52 putative cyclin-dependent kinases (CDKs) and CDK-like (CDKLs) genes in Arachis hypogaea (cultivated peanut) and total 26 genes in each diploid parent of cultivated peanut (Arachis duranensis and Arachis ipaensis). Both CDK and CDKL genes were classified into eight groups based on their cyclin binding motifs and their phylogenetic relationship with Arabidopsis counterparts. Genes in the same subgroup displayed similar exon-intron structure and conserved motifs. Further, gene duplication analysis suggested that segmental duplication events played major roles in the expansion and evolution of CDK and CDKL genes in cultivated peanuts. Identification of diverse cis-acting response elements in CDK and CDKL genes promoter indicated their potential fundamental roles in multiple biological processes. Various gene expression patterns of CDKs and CDKLs in different peanut tissues suggested their involvement during growth and development. In addition, qRT-PCR analysis demonstrated that most representing CDK and CDKL gene family members were significantly down-regulated under ABA, PEG and mannitol treatments.

CONCLUSIONS

Genome-wide analysis offers a comprehensive understanding of the classification, evolution, gene structure, and gene expression profiles of CDK and CDKL genes in cultivated peanut and their diploid progenitors. Additionally, it also provides cell cycle regulatory gene resources for further functional characterization to enhance growth, development and abiotic stress tolerance.

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.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

Genome-Wide Association Study and Genomic Prediction for Bacterial Wilt Resistance in Common Bean (Phaseolus vulgaris) Core Collection

Common bean (Phaseolus vulgaris) is one of the major legume crops cultivated worldwide. Bacterial wilt (BW) of common bean (Curtobacterium flaccumfaciens pv. flaccumfaciens), being a seed-borne disease, has been a challenge in common bean producing regions. A genome-wide association study (GWAS) was conducted to identify SNP markers associated with BW resistance in the USDA common bean core collection. A total of 168 accessions were evaluated for resistance against three different isolates of BW. Our study identified a total of 14 single nucleotide polymorphism (SNP) markers associated with the resistance to BW isolates 528, 557, and 597 using mixed linear models (MLMs) in BLINK, FarmCPU, GAPIT, and TASSEL 5. These SNPs were located on chromosomes Phaseolus vulgaris [Pv]02, Pv04, Pv08, and Pv09 for isolate 528; Pv07, Pv10, and Pv11 for isolate 557; and Pv04, Pv08, and Pv10 for isolate 597. The genomic prediction accuracy was assessed by utilizing seven GP models with 1) all the 4,568 SNPs and 2) the 14 SNP markers. The overall prediction accuracy (PA) ranged from 0.30 to 0.56 for resistance against the three BW isolates. A total of 14 candidate genes were discovered for BW resistance located on chromosomes Pv02, Pv04, Pv07, Pv08, and Pv09. This study revealed vital information for developing genetic resistance against the BW pathogen in common bean. Accordingly, the identified SNP markers and candidate genes can be utilized in common bean molecular breeding programs to develop novel resistant cultivars.

The fate of drought-related genes after polyploidization in Arachis hypogaea cv. Tifrunner

Drought stress affects plant growth and development. Cultivated peanut (Arachis hypogaea) was formed by a cross between A. duranensis and A. ipaensis. The drought tolerance of A. duranensis and A. ipaensis is reportedly stronger than that of cultivated peanut. However, there has been little study of drought tolerance genes in Arachis. In this study, we compared drought tolerance genes between A. hypogaea cv. Tifrunner and its diploid donors. We have observed that polyploidization does not generate more drought tolerance genes in A. hypogaea cv. Tifrunner but promotes the loss of many ancient drought tolerance genes. Although putative drought tolerance genes occurred on gene duplication events in A. hypogaea cv. Tifrunner, most copies lacked drought tolerance. These findings suggest that the loss of drought tolerance genes in A. hypogaea cv. Tifrunner could possibly result in weaker drought tolerance. In addition, we have observed that the three Arachis species stochastically lost putative drought tolerance genes. The evolution of drought tolerance genes could possibly have correlated with environmental changes. Our results enhance the current understanding of drought tolerance and polyploidy evolution in Arachis species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01198-0.

Dissection of valine-glutamine genes and their responses to drought stress in Arachis hypogaea cv. Tifrunner

Valine-glutamine sequences (VQs) interact with WRKY transcription factors (TFs), forming VQ-WRKY protein complexes crucial for plant development and response to environmental changes. Cultivated peanut (Arachis hypogaea) is a tetraploid from A. duranensis and A. ipaensis cross. The Arachis spp. WRKY TFs have been identified, but Arachis VQs are largely unknown. This study identified VQs in A. duranensis, A. ipaensis, A. monticola, A. hypogaea cv. Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The study analyzed the homologous relationships between VQs in these Arachis spp. The VQ drought-tolerant genes were detected and VQ-WRKY interactions were determined in A. hypogaea cv. Tifrunner. The results showed that tetraploid Arachis spp. retained duplicated VQs, but lost ancestral VQs after allopolyploidization. The number of VQs in A. monticola, A. hypogaea cv. Fuhuasheng, and A. hypogaea cv. Shitouqi increased relative to their diploid ancestors. RNA-seq and quantitative real-time PCR experiments confirmed that three AhTVQs tolerate drought stress in A. hypogaea cv. Tifrunner. However, evidence of VQ-WRKY interaction for drought stress response is lacking in A. hypogaea cv. Tifrunner. Nevertheless, this study identified VQ-WRKY interactions, which possibly have multiple functions in A. hypogaea cv. Tifrunner. Altogether, this study dissected Arachis VQs, providing insights into Arachis VQ evolution and drought function.

Genome-wide identification of auxin response factor (ARF) gene family and the miR160-ARF18-mediated response to salt stress in peanut (Arachis hypogaea L.)

Auxin response factors (ARFs) are transcription factors that regulate the transcription of auxin-responsive genes during plant growth and development. In this study, 29 and 30 ARF members were identified from the two wild peanut species, A. duranensis and A. ipaensis, respectively. The ARFs, including their classifications, conserved domains and evolutionary relationships were characterized. RNA-seq analyses revealed that some of the ARF genes were responsive to abiotic stress, particularly high salinity. In addition to abiotic stress, the expression of 2 ARF members was also regulated by biotic stress, specifically Bradyrhizobium infection in A. duranensis. The ARF gene Arahy.7DXUOK was predicted to be a potential target of miR160. Overexpression of miR160 could cause degradation of the Arahy.7DXUOK target gene transcript and increased salt tolerance in miR160OX transgenic plants. Therefore, these molecular characterization and expression profile analyses provide comprehensive information on ARF family members and will help to elucidate their functions to facilitate further research on peanuts.

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.

Muscle synergy differences between voluntary and reactive backward stepping

Reactive stepping responses are essential to prevent falls after a loss of balance. It has previously been well described that both voluntary and reactive step training could improve the efficacy of reactive stepping in different populations. However, the effect of aging on neuromuscular control during voluntary and reactive stepping remains unclear. Electromyography (EMG) signals during both backward voluntary stepping in response to an auditory cue and backward reactive stepping elicited by a forward slip-like treadmill perturbation during stance were recorded in ten healthy young adults and ten healthy older adults. Using muscle synergy analysis, we extracted the muscle synergies for both voluntary and reactive stepping. Our results showed that fewer muscle synergies were used during reactive stepping than during voluntary stepping in both young and older adults. Minor differences in the synergy structure were observed for both voluntary and reactive stepping between age groups. Our results indicate that there is a low similarity of muscle synergies between voluntary stepping and reactive stepping and that aging had a limited effect on the structure of muscle synergies. This study enhances our understanding of the neuromuscular basis of both voluntary and reactive stepping as well as the potential effect of aging on neuromuscular control during balance tasks.

Association of differentially expressed R-gene candidates with leaf spot resistance in peanut (Arachis hypogaea L.)

Early leaf spot (ELS) and late leaf spot (LLS) are major fungal diseases of peanut that can severely reduce yield and quality. Development of acceptable genetic resistance has been difficult due to a strong environmental component and many major and minor QTLs. Resistance genes (R-genes) are an important component of plant immune system and have been identified in peanut. Association of specific R-genes to leaf spot resistance will provide molecular targets for marker-assisted breeding strategies. In this study, advanced breeding lines from different pedigrees were evaluated for leaf spot resistance and 76 candidate R-genes expression study was applied to susceptible and resistant lines. Thirty-six R-genes were differentially expressed and significantly correlated with resistant lines, of which a majority are receptor like kinases (RLKs) and receptor like proteins (RLPs) that sense the presence of pathogen at the cell surface and initiate protection response. The largest group was receptor-like cytoplasmic kinases (RLCKs) VII that are involved in pattern-triggered kinase signaling resulting in the production reactive oxygen species (ROS). Four R-genes were homologous to TMV resistant protein N which has shown to confer resistance against tobacco mosaic virus (TMV). When mapped to peanut genomes, 36 R-genes were represented in most chromosomes except for A09 and B09. Low levels of gene-expression in resistant lines suggest expression is tightly controlled to balance the cost of R-gene expression to plant productively. Identification and association of R-genes involved in leaf spot resistance will facilitate genetic selection of leaf spot resistant lines with good agronomic traits.

QTL mapping of web blotch resistance in peanut by high-throughput genome-wide sequencing

BACKGROUND

Web blotch is one of the most important foliar diseases worldwide in peanut (Arachis hypogaea L.). The identification of quantitative trait loci (QTLs) for peanut web blotch resistance represents the basis for gene mining and the application of molecular breeding technologies.

RESULTS

In this study, a peanut recombinant inbred line (RIL) population was used to map QTLs for web blotch resistance based on high-throughput genome-wide sequencing. Frequency distributions of disease grade and disease index in five environments indicated wide phenotypic variations in response to web blotch among RILs. A high-density genetic map was constructed, containing 3634 bin markers distributed on 20 peanut linkage groups (LGs) with an average genetic distance of 0.5 cM. In total, eight QTLs were detected for peanut web blotch resistance in at least two environments, explaining from 2.8 to 15.1% of phenotypic variance. Two major QTLs qWBRA04 and qWBRA14 were detected in all five environments and were linked to 40 candidate genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) or other proteins related to disease resistances.

CONCLUSIONS

The results of this study provide a basis for breeding peanut cultivars with web blotch resistance.

Genome-Wide Identification and Evolutionary Analysis of NBS-LRR Genes From Dioscorea rotundata

Dioscorea rotundata is an important food crop that is mainly cultivated in subtropical regions of the world. D. rotundata is frequently infected by various pathogens during its lifespan, which results in a substantial economic loss in terms of yield and quality. The disease resistance gene (R gene) profile of D. rotundata is largely unknown, which has greatly hampered molecular study of disease resistance in this species. Nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes are the largest group of plant R genes, and they play important roles in plant defense responses to various pathogens. In this study, 167 NBS-LRR genes were identified from the D. rotundata genome. Subsequently, one gene was assigned to the resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) subclass and the other 166 genes to the coiled coil (CC)-NBS-LRR (CNL) subclass. None of the Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes were detected in the genome. Among them, 124 genes are located in 25 multigene clusters and 43 genes are singletons. Tandem duplication serves as the major force for the cluster arrangement of NBS-LRR genes. Segmental duplication was detected for 18 NBS-LRR genes, although no whole-genome duplication has been documented for the species. Phylogenetic analysis revealed that D. rotundata NBS-LRR genes share 15 ancestral lineages with Arabidopsis thaliana genes. The NBS-LRR gene number increased by more than a factor of 10 during D. rotundata evolution. A conservatively evolved ancestral lineage was identified from D. rotundata, which is orthologs to the Arabidopsis RPM1 gene. Transcriptome analysis for four different tissues of D. rotundata revealed a low expression of most NBS-LRR genes, with the tuber and leaf displaying a relatively high NBS-LRR gene expression than the stem and flower. Overall, this study provides a complete set of NBS-LRR genes for D. rotundata, which may serve as a fundamental resource for mining functional NBS-LRR genes against various pathogens.