Isolation of a family of resistance gene analogue sequences of the nucleotide binding site (NBS) type from Lens species

Conventional and new-breeding technologies for improving disease resistance in lentil (Lens culinaris Medik)

Lentil, an important cool season food legume, is a rich source of easily digestible protein, folic acid, bio-available iron, and zinc nutrients. Lentil grows mainly as a sole crop in the winter after harvesting rice in South Asia. However, the annual productivity is low due to its slow growth during the early phase, competitive weed infestation, and disease outbreaks during the crop growth period. Disease resistance breeding has been practiced for a long time to enhance resistance to various diseases. Often the sources of resistance are available in wild crop relatives. Thus, wide hybridization and the ovule rescue technique have helped to introgress the resistance trait into cultivated lentils. Besides hybridization, induced mutagenesis contributed immensely in creating variability for disease tolerance, and several disease-resistant mutant lines have been developed. However, to overcome the limitations of traditional breeding approaches, advancement in molecular marker technologies, and genomics has helped to develop disease-resistant and climate-resilient lentil varieties with more precision and efficiency. This review describes types of diseases, disease screening methods, the role of conventional and new breeding technologies in alleviating disease-incurred damage and progress toward making lentil varieties more resilient to disease outbreaks under the shadow of climate change.

Resistance screening and in silico characterization of cloned novel RGA from multi-race resistant lentil germplasm against Fusarium wilt (Fusarium oxysporum f. sp. lentis)

Fusarium wilt caused by Fusarium oxysporum f. sp. lentis (Fol) is the most devastating disease of lentil present worldwide. Identification of multi-race fusarium wilt resistance genes and their incorporation into existing cultivars will help to reduce yield losses. In the present study, 100 lentil germplasms belonging to seven lentil species were screened against seven prevalent races of Fol, and accessions IC201561 (Lens culinaris subsp. culinaris), EC714243 (L. c. subsp. odemensis), and EC718238 (L. nigricans) were identified as resistant. The typical R gene codes for the nucleotide-binding site and leucine-rich repeats (NBS-LRR) at the C terminal are linked to either the Toll/interleukin 1-like receptor (TIR) or coiled coil (CC) at the N terminal. In the present study, degenerate primers, designed from the NBS region amplifying the P-loop to the GLPLA motif, isolated forty-five resistance gene analogues (RGAs) from identified resistant accessions. The sequence alignment identified both classes of RGAs, TIR and non-TIR, based on the presence of aspartate (D) and tryptophan (W) at the end of the kinase motif, respectively. The phylogenetic analysis grouped the RGAs into six classes, from LRGA1 to LRGA6, which determined the diversity of the RGAs present in the host. Grouping of the RGAs identified from Lens nigricans, LnRGA 2, 9, 13 with I2 revealed the structural similarity with the fusarium resistance gene. The similarity index ranged from 27.85% to 86.98% among the RGAs and from 26.83% to 49.41% among the known R genes, I2, Gpa2, M, and L6. The active binding sites present along the conserved motifs grouped the RGAs into 13 groups. ADP/ATP, being the potential ligand, determines the ATP binding and ATP hydrolysis activity of the RGAs. The isolated RGAs can be used to develop markers linked to the functional R gene. Furthermore, expression analysis and full-length gene isolation pave the path to identifying the molecular mechanism involved in resistance.

NB-LRR Lineage-Specific Equipment Is Sorted Out by Sequence Pattern Adaptation and Domain Segment Shuffling

The nucleotide-binding and leucine-rich repeat (NB-LRR) genes, also known as resistance (R)-genes, play an important role in the activation of immune responses. In recent years, large-scale studies have been performed to highlight the diversification of plant NB-LRR repertories. It is well known that, to provide new functionalities, NB-LRR sequences are subject to duplication, domain fusions and acquisition and other kinds of mutations. Although some mechanisms that govern NB-LRR protein domain adaptations have been uncovered, to retrace the plant-lineage-specific evolution routes of R protein structure, a multi-genome comparative analysis was performed. This study allowed us to define groups of genes sharing homology relationships across different species. It is worth noting that the most populated groups contained well-characterized R proteins. The arsenal profile of such groups was investigated in five botanical families, including important crop species, to underline specific adaptation signatures. In addition, the dissection of 70 NB domains of well-characterized R-genes revealed the NB core motifs from which the three main R protein classes have been diversified. The structural remodeling of domain segments shaped the specific NB-LRR repertoires observed in each plant species. This analysis provided new evolutionary and functional insights on NB protein domain shuffling. Taken together, such findings improved our understanding of the molecular adaptive selection mechanisms occurring at plant R loci.

Isolation and characterization of fusarium wilt resistance gene analogs in radish

The resistance gene analog (RGA)-based marker strategy is an effective supplement for current marker reservoir of radish disease-resistance breeding. In this study, we identified RGAs based on the conserved nucleotide-binding site (NBS) and S-receptor-like kinase (SRLK) domains. A total of 68 NBS-RGAs and 46 SRLK-RGAs were isolated from two FW-resistant radish inbred lines, B2 and YR31, and one susceptible line, YR15. A BLASTx search revealed that the NBS-RGAs contained six conserved motifs (i.e., P loop, RNBS-A, Kinase-2, RNBS-B, RNBS-C, and GLPL) and the SRLK-RGAs, contained two conserved motifs (i.e., G-type lectin and PAN-AP). A phylogenetic analysis indicated that the NBS-RGAs could be separated into two classes (i.e., toll/interleukin receptor and coiled-coil types), with six subgroups, and the SRLK-RGAs were divided into three subgroups. Moreover, we designed RGA-specific markers from data-mining approach in radish databases. Based on marker analysis, 24 radish inbred lines were clustered into five main groups with a similarity index of 0.44 and showing genetic diversity with resistance variation in those radish inbred lines. The development of RGA-specific primers would be valuable for marker-assisted selection during the breeding of disease-resistant radish cultivars.

Genome-wide identification and tissue-specific expression analysis of nucleotide binding site-leucine rich repeat gene family in Cicer arietinum (kabuli chickpea)

The nucleotide binding site-leucine rich repeat (NBS-LRR) proteins play an important role in the defense mechanisms against pathogens. Using bioinformatics approach, we identified and annotated 104 NBS-LRR genes in chickpea. Phylogenetic analysis points to their diversification into two families namely TIR-NBS-LRR and non-TIR-NBS-LRR. Gene architecture revealed intron gain/loss events in this resistance gene family during their independent evolution into two families. Comparative genomics analysis elucidated its evolutionary relationship with other fabaceae species. Around 50% NBS-LRRs reside in macro-syntenic blocks underlining positional conservation along with sequence conservation of NBS-LRR genes in chickpea. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. Four cis-regulatory elements namely WBOX, DRE, CBF, and GCC boxes, that commonly occur in resistance genes, were present in the promoter regions of these genes. Further, the findings will provide a strong background to use candidate disease resistance NBS-encoding genes and identify their specific roles in chickpea.

Current knowledge in lentil genomics and its application for crop improvement

Most of the lentil growing countries face a certain set of abiotic and biotic stresses causing substantial reduction in crop growth, yield, and production. Until-to date, lentil breeders have used conventional plant breeding techniques of selection-recombination-selection cycle to develop improved cultivars.These techniques have been successful in mainstreaming some of the easy-to-manage monogenic traits. However, in case of complex quantitative traits, these conventional techniques are less precise. As most of the economic traits are complex, quantitative, and often influenced by environments and genotype-environment interaction, the genetic improvement of these traits becomes difficult. Genomics assisted breeding is relatively powerful and fast approach to develop high yielding varieties more suitable to adverse environmental conditions. New tools such as molecular markers and bioinformatics are expected to generate new knowledge and improve our understanding on the genetics of complex traits. In the past, the limited availability of genomic resources in lentil could not allow breeders to employ these tools in mainstream breeding program.The recent application of the next generation sequencing and genotyping by sequencing technologies has facilitated to speed up the lentil genome sequencing project and large discovery of genome-wide single nucleotide polymorphism (SNP) markers. Currently, several linkage maps have been developed in lentil through the use of expressed sequenced tag (EST) derived simple sequence repeat (SSR) and SNP markers.These maps have emerged as useful genomic resources to identify quantitative trait loci imparting tolerance to biotic and abiotic stresses in lentil. In this review, the current knowledge on available genomic resources and its application in lentil breeding program are discussed.

Identification of immunity related genes to study the Physalis peruviana--Fusarium oxysporum pathosystem

The Cape gooseberry (Physalisperuviana L) is an Andean exotic fruit with high nutritional value and appealing medicinal properties. However, its cultivation faces important phytosanitary problems mainly due to pathogens like Fusarium oxysporum, Cercosporaphysalidis and Alternaria spp. Here we used the Cape gooseberry foliar transcriptome to search for proteins that encode conserved domains related to plant immunity including: NBS (Nucleotide Binding Site), CC (Coiled-Coil), TIR (Toll/Interleukin-1 Receptor). We identified 74 immunity related gene candidates in P. peruviana which have the typical resistance gene (R-gene) architecture, 17 Receptor like kinase (RLKs) candidates related to PAMP-Triggered Immunity (PTI), eight (TIR-NBS-LRR, or TNL) and nine (CC–NBS-LRR, or CNL) candidates related to Effector-Triggered Immunity (ETI) genes among others. These candidate genes were categorized by molecular function (98%), biological process (85%) and cellular component (79%) using gene ontology. Some of the most interesting predicted roles were those associated with binding and transferase activity. We designed 94 primers pairs from the 74 immunity-related genes (IRGs) to amplify the corresponding genomic regions on six genotypes that included resistant and susceptible materials. From these, we selected 17 single band amplicons and sequenced them in 14 F. oxysporum resistant and susceptible genotypes. Sequence polymorphisms were analyzed through preliminary candidate gene association, which allowed the detection of one SNP at the PpIRG-63 marker revealing a nonsynonymous mutation in the predicted LRR domain suggesting functional roles for resistance.

Isolation and characterization of nucleotide-binding site resistance gene homologues in common bean (Phaseolus vulgaris)

Common bean production is constrained by many fungal, viral, and bacterial pathogens. Thus, the identification of resistance (R) genes is an important focal point of common bean research. The main goal of our study was to identify resistance gene homologues (RGH) in the crop, using degenerate primers designed from conserved sequences in the nucleotide-binding site (NBS) domains of R-genes from the model legume Medicago truncatula. Total DNA of the Andean common bean genotype G19833 was used for amplification of over 500 primer combinations. Sequencing of amplicons showed that 403 cloned fragments had uninterrupted open reading frames and were considered representative of functional RGH genes. The sequences were grouped at two levels of nucleotide identity (90 and 80%) and representative sequences of each group were used for phylogenetic analyses. The RGH sequence diversity of common bean was divided into TIR and non-TIR families, each with different clusters. The TIR sequences grouped into 14 clades while non-TIR sequences grouped into seven clades. Pairwise comparisons showed purifying selection, although some sequences may have been the result of diversifying selection. Knowledge about RGH genes in common bean can allow the design of molecular markers for pyramiding of resistance genes against various pathogens.

Genome-wide identification and characterization of nucleotide binding site leucine-rich repeat genes in linseed reveal distinct patterns of gene structure

Plants employ different disease-resistance genes to detect pathogens and to induce defense responses. The largest class of these genes encodes proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains. To identify the putative NBS-LRR encoding genes from linseed, we analyzed the recently published linseed genome sequence and identified 147 NBS-LRR genes. The NBS domain was used for phylogeny construction and these genes were classified into two well-known families, non-TIR (CNL) and TIR related (TNL), and formed eight clades in the neighbor-joining bootstrap tree. Eight different gene structures were observed among these genes. An unusual domain arrangement was observed in the TNL family members, predominantly in the TNL-5 clade members belonging to class D. About 12% of the genes observed were linseed specific. The study indicated that the linseed genes probably have an ancient origin with few progenitor genes. Quantitative expression analysis of five genes showed inducible expression. The in silico expression evidence was obtained for a few of these genes, and the expression was not correlated with the presence of any particular regulatory element or with unusual domain arrangement in those genes. This study will help in understanding the evolution of these genes, the development of disease resistant varieties, and the mechanism of disease resistance in linseed.

Integration of new CAPS and dCAPS-RGA markers into a composite chickpea genetic map and their association with disease resistance

A composite linkage map was constructed based on two interspecific recombinant inbred line populations derived from crosses between Cicer arietinum (ILC72 and ICCL81001) and Cicer reticulatum (Cr5-10 or Cr5-9). These mapping populations segregate for resistance to ascochyta blight (caused by Ascochyta rabiei), fusarium wilt (caused by Fusarium oxysporum f. sp. ciceris) and rust (caused by Uromyces ciceris-arietini). The presence of single nucleotide polymorphisms in ten resistance gene analogs (RGAs) previously isolated and characterized was exploited. Six out of the ten RGAs were novel sequences. In addition, classes RGA05, RGA06, RGA07, RGA08, RGA09 and RGA10 were considerate putatively functional since they matched with several legume expressed sequences tags (ESTs) obtained under infection conditions. Seven RGA PCR-based markers (5 CAPS and 2 dCAPS) were developed and successfully genotyped in the two progenies. Six of them have been mapped in different linkage groups where major quantitative trait loci conferring resistance to ascochyta blight and fusarium wilt have been reported. Genomic locations of RGAs were compared with those of known Cicer R-genes and previously mapped RGAs. Association was detected between RGA05 and genes controlling resistance to fusarium wilt caused by races 0 and 5.

Characterization of resistance gene analogs with a nucleotide binding site isolated from a triploid white poplar

The majority of cloned plant disease resistance genes (R genes) encode a nucleotide binding site (NBS) and a leucine-rich repeat (LRR) domain. In this study, to better understand the R genes in white poplar, 59 resistance gene analogues (RGAs) were identified from a triploid white poplar [(Populus tomentosa x Populus bolleana) x P. tomentosa], based on conserved NBS regions. The 59 RGAs were phylogenetically classified into 10 subfamilies, and 54 RGAs with open-reading frames (ORFs) were further grouped into two classes, toll and interleukin-1 receptor (TIR) and non-TIR. BLAST searches with reference to the genomic sequence of Populus trichocarpa found 96 highly homologous regions distributed in 37 loci, suggesting the abundance and divergence of NBS-encoding genes in the triploid poplar genome. Within subfamilies 1-3, the average non-synonymous/synonymous substitution (omega) rates were < 1, indicating purifying selection on these RGAs, but some sites were clearly under diversifying selection with omega > 1. Many intergenic exchanges were also detected among these RGAs, indicating a probable role in homogenising NBS domains. Quantitative real-time PCR analysis revealed dramatic variations in the transcript level of 18 RGAs in the mature leaves, bark and roots of the triploid poplar, and identified two RGAs that had significantly higher level of transcripts in bark, four RGAs in mature leaves, and 14 in the above-ground portion of poplars, suggesting their probable roles in resistance against diseases attacking the organs. Our results shed light on genetic resources of poplar resistance and will be useful for further resistance gene isolation and exploitation.

Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula

The nucleotide-binding site (NBS)-Leucine-rich repeat (LRR) gene family accounts for the largest number of known disease resistance genes, and is one of the largest gene families in plant genomes. We have identified 333 nonredundant NBS-LRRs in the current Medicago truncatula draft genome (Mt1.0), likely representing 400 to 500 NBS-LRRs in the full genome, or roughly 3 times the number present in Arabidopsis (Arabidopsis thaliana). Although many characteristics of the gene family are similar to those described on other plant genomes, several evolutionary features are particularly pronounced in M. truncatula, including a high degree of clustering, evidence of significant numbers of ectopic translocations from clusters to other parts of the genome, a small number of more evolutionarily stable NBS-LRRs, and numerous truncations and fusions leading to novel domain compositions. The gene family clearly has had a large impact on the structure of the genome, both through ectopic translocations (potentially, a means of seeding new NBS-LRR clusters), and through two extraordinarily large superclusters. Chromosome 6 encodes approximately 34% of all TIR-NBS-LRRs, while chromosome 3 encodes approximately 40% of all coiled-coil-NBS-LRRs. Almost all atypical domain combinations are in the TIR-NBS-LRR subfamily, with many occurring within one genomic cluster. This analysis shows the gene family not only is important functionally and agronomically, but also plays a structural role in the genome.

Diversity and evolutionary relationship of nucleotide binding site-encoding disease-resistance gene analogues in sweet potato (Ipomoea batatas Lam.)

Most plant disease-resistance genes (R-genes) isolated so far encode proteins with a nucleotide binding site (NBS) domain and belong to a superfamily. NBS domains related to R-genes show a highly conserved backbone of an amino acid motif, which makes it possible to isolate resistance gene analogues (RGAs) by degenerate primers. Degenerate primers based on the conserved motif (P-loop and GLPL) of the NBS domain from R -genes were used to isolate RGAs from the genomic DNA of sweet potato cultivar Qingnong no.2. Five distinct clusters of RGAs (22 sequences) with the characteristic NBS representing a highly diverse sample were identified in sweet potato genomic DNA. Sequence identity among the 22 RGA nucleotide sequences ranged from 41.2% to 99.4%, while the deduced amino acid sequence identity from the 22 RGAs ranged from 20.6%to 100%. The analysis of sweet potato RGA sequences suggested mutation as the primary source of diversity. The phylogenetic analyses for RGA nucleotide sequences and deduced amino acids showed that RGAs from sweet potato were classified into two distinct groups--toll and interleukin receptor-1 (TIR)-NBS-LRR and non-TIR-NBS-LRR. The high degree of similarity between sweet potato RGAs and NBS sequences derived from R-genes cloned from tomato, tobacco, flax and potato suggest an ancestral relationship. Further studies showed that the ratio of non-synonymous to synonymous substitution within families was low. These data obtained from sweet potato suggest that the evolution of NBS-encoding sequences in sweet potato occur by the gradual accumulation of mutations leading to purifying selection and slow rates of divergence within distinct R-gene families.

Identification and characterization of NBS-LRR class resistance gene analogs in faba bean (Vicia faba L.) and chickpea (Cicer arietinum L.)

A PCR approach with degenerate primers designed from conserved NBS-LRR (nucleotide binding site-leucine-rich repeat) regions of known disease-resistance (R) genes was used to amplify and clone homologous sequences from 5 faba bean (Vicia faba) lines and 2 chickpea (Cicer arietinum) accessions. Sixty-nine sequenced clones showed homologies to various R genes deposited in the GenBank database. The presence of internal kinase-2 and kinase-3a motifs in all the sequences isolated confirm that these clones correspond to NBS-containing genes. Using an amino-acid sequence identity of 70% as a threshold value, the clones were grouped into 10 classes of resistance-gene analogs (RGA01 to RGA10). The number of clones per class varied from 1 to 30. RGA classes 1, 6, 8, and 9 were comprised solely of clones isolated from faba bean, whereas classes 2, 3, 4, 5, and 7 included only chickpea clones. RGA10, showing a within-class identity of 99%, was the only class consisting of both faba bean and chickpea clones. A phylogenetic tree, based on the deduced amino-acid sequences of 12 representative clones from the 10 RGA classes and the NBS domains of 6 known R genes (I2 and Prf from tomato, RPP13 from Arabidopsis, Gro1-4 from potato, N from tobacco, L6 from flax), clearly indicated the separation between TIR (Toll/interleukin-1 receptor homology: Gro1-4, L6, N, RGA05 to RGA10)- and non-TIR (I2, Prf, RPP13, RGA01 to RGA04)-type NBS-LRR sequences. The development of suitable polymorphic markers based on cloned RGA sequences to be used in genetic mapping will facilitate the assessment of their potential linkage relationships with disease-resistance genes in faba bean and chickpea. This work is the first to report on faba bean RGAs.

Isolation and characterization of a set of disease resistance-gene analogs (RGAs) from wild rice, Zizania latifolia Griseb. I. Introgression, copy number lability, sequence change, and DNA methylation alteration in several rice-Zizania introgression lines

Eight resistance-gene analogs (RGAs) were isolated from wild rice, Zizania latifolia Griseb., by degenerate primers designed according to conserved motifs at or around the nucleotide-binding site (NBS) of known NBS-containing plant resistance genes. The 8 RGAs were classified into 6 distinct groups based on their deduced amino acid sequence similarity of 60% or greater. Gel-blot hybridization of each of the RGAs to 4 rice - Z. latifolia intro gression lines indicated an array of changes at either introgressed Zizania RGAs or, more likely, their rice homologs. The changes included dramatic increase in copy number, modification at the primary DNA sequence, and alteration in DNA methylation patterns.

Plant NBS-LRR proteins: adaptable guards

Plant nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins are a large family involved in disease resistance; they may monitor the status of proteins targeted by pathogens.

The majority of disease resistance genes in plants encode nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins. This large family is encoded by hundreds of diverse genes per genome and can be subdivided into the functionally distinct TIR-domain-containing (TNL) and CC-domain-containing (CNL) subfamilies. Their precise role in recognition is unknown; however, they are thought to monitor the status of plant proteins that are targeted by pathogen effectors.