Genetic Variation in NBS Genes: Analysis in Tolerant Accessions and Implications for Disease Resistance

Abstract

This article provides a comprehensive analysis of genetic variation in Nucleotide-Binding Site (NBS) genes, a key class of disease resistance genes in plants, with a focus on tolerant accessions. It explores the foundational diversity and evolution of NBS genes across species, details methodological approaches for their identification and variation analysis, addresses challenges in sequencing and data interpretation, and presents validation strategies through functional studies. Aimed at researchers and scientists in genetics and drug development, the content synthesizes current research to highlight how genetic variation in NBS genes of tolerant genotypes underpins resistance mechanisms, offering insights for breeding and therapeutic discovery.

Unraveling NBS Gene Diversity and Evolutionary Patterns in Plant Genomes

The Nucleotide-Binding Site (NBS) gene superfamily represents one of the most extensive classes of plant resistance (R) genes, forming the cornerstone of the innate immune system against pathogens such as viruses, bacteria, fungi, and oomycetes [1] [2] [3]. These genes encode proteins characterized by a central NBS domain that facilitates nucleotide binding and hydrolysis, which is crucial for signal transduction during plant defense responses [3] [4]. The NBS domain is typically flanked by various N-terminal and C-terminal domains, leading to a standardized classification system based on domain architecture [2] [3].

Table 1: Standard Classification of NBS-LRR Genes Based on Domain Architecture

Class	Domain Architecture	Description	Prevalence Examples
TNL	TIR-NBS-LRR	N-terminal TIR domain, NBS, C-terminal LRR	64 in N. tabacum [2]
CNL	CC-NBS-LRR	N-terminal Coiled-Coil, NBS, C-terminal LRR	74 in N. tabacum [2]
NL	NBS-LRR	Only NBS and LRR domains	306 in N. tabacum [2]
TN	TIR-NBS	Only TIR and NBS domains	9 in N. tabacum [2]
CN	CC-NBS	Only CC and NBS domains	150 in N. tabacum [2]
N	NBS	Only the NBS domain	127 in N. tomentosiformis [2]

The NBS-encoding genes represent a significant portion of plant genomes, though their numbers vary dramatically between species. A recent comparative analysis across 34 plant species, ranging from mosses to monocots and dicots, identified 12,820 NBS-domain-containing genes, revealing 168 distinct classes of domain architecture [1]. This diversity encompasses both classical patterns (e.g., NBS, NBS-LRR, TIR-NBS-LRR) and species-specific structural patterns [1].

Evolution and Genomic Distribution

The expansion and evolution of the NBS gene family are primarily driven by gene duplication events, including both tandem and segmental duplications [5] [3]. In many plant genomes, such as apple (Malus x domestica), these duplication events have played a major role in the family's expansion, allowing for the generation of new pathogen recognition specificities [3]. Studies on Rosaceae species and pear genomes have demonstrated that proximal duplications are a key factor leading to differences in NBS gene numbers between closely related species [5]. Furthermore, analyses of orthologous gene pairs often reveal a Ka/Ks ratio >1, indicating that positive selection has acted upon these genes following species divergence, which is consistent with an evolutionary "arms race" with rapidly evolving pathogens [5].

Table 2: NBS Gene Family Size Across Selected Plant Species

Plant Species	Number of NBS Genes	Noteworthy Characteristics	Source
Malus x domestica (Apple)	1,015	Approximately equal distribution of TIR (TNL) and CC (CNL) domains [3]	[3]
Nicotiana tabacum	603	Allotetraploid; ~76.6% of members traceable to parental genomes [2]	[2]
Pyrus bretschneideri (Asian Pear)	338	Difference with European pear mainly due to proximal duplications [5]	[5]
Pyrus communis (European Pear)	412	Difference with Asian pear mainly due to proximal duplications [5]	[5]
Nicotiana benthamiana	156	Model plant for virology; includes 5 TNL and 25 CNL types [4]	[4]

Molecular Function and Signaling Mechanisms

NBS-LRR proteins function as intracellular immune receptors that monitor for pathogen presence. Their activation triggers a robust defense response, often culminating in the Hypersensitive Response (HR), a form of programmed cell death at the infection site that restricts pathogen spread [3]. The functional mechanism involves distinct roles for different protein domains:

• LRR Domain: The leucine-rich repeat region is primarily responsible for the specific recognition of pathogen-derived effector proteins (avirulence factors). This can occur through direct binding or by monitoring the status of host proteins that are modified by pathogens [4].
• NBS Domain: Upon effector recognition, the nucleotide-binding site domain undergoes a conformational shift from an ADP-bound state to an ATP-bound state. This shift is essential for activating downstream signaling cascades [4].
• N-terminal Domains (TIR/CC): The activated TIR or CC domains initiate signaling. TNL and CNL proteins often utilize different downstream signaling pathways, suggesting functional diversification within the superfamily [3] [4].

The subcellular localization of NBS-LRR proteins is diverse, encompassing the cytoplasm, plasma membrane, and nucleus, which aligns with their roles in detecting pathogens and transducing signals across different cellular compartments [4].

Protocol: Genome-Wide Identification and Characterization of NBS Genes

This protocol is adapted from multiple genomic studies [2] [5] [3] and provides a standard workflow for identifying and analyzing the NBS gene superfamily in a plant genome of interest.

Step 1: Identification of NBS-Encoding Genes

• HMMER Search: Use the hmmsearch program from HMMER v3.1b2 or later to scan the plant's proteome or predicted protein sequences. The standard Hidden Markov Model (HMM) profile for the NBS domain is PF00931 (NB-ARC) obtained from the Pfam database. An E-value cutoff of < 1e-04 is commonly used [3] [4]. For increased sensitivity, some studies employ a more stringent cutoff (e.g., < 1e-20) [4] or build a custom, species-specific HMM profile from initial high-confidence hits [3].
• Domain Verification: Confirm the presence of the NBS domain in all candidate sequences using tools like PfamScan or the NCBI Conserved Domain Database (CDD). Manually verify the completeness of the domain. This step is crucial for removing false positives [2] [3].

Step 2: Classification and Structural Analysis

• Amino-Terminal Analysis: Identify the N-terminal domain (TIR, CC, or other) in the candidate proteins. The TIR domain can be identified using Pfam models (e.g., PF01582, PF00560), while the CC domain is often confirmed via the NCBI CDD or prediction tools like COILS [2].
• LRR Identification: Scan for C-terminal LRR domains using Pfam models (e.g., PF07723, PF07725, PF12779, PF13516, PF13855) [2]. Based on the presence/absence of TIR, CC, and LRR domains, classify the genes into the standard classes (TNL, CNL, NL, TN, CN, N) [4].

Step 3: Phylogenetic and Evolutionary Analysis

• Multiple Sequence Alignment: Perform a multiple alignment of the full-length NBS protein sequences or the NBS domains using tools like MUSCLE or ClustalW with default parameters [2] [4].
• Phylogenetic Tree Construction: Use software such as MEGA11 to build a phylogenetic tree (e.g., using the Maximum Likelihood method) with a bootstrap test (e.g., 1000 replicates) to assess node reliability [2] [5]. This tree helps visualize the evolutionary relationships and major clades (e.g., TNL vs. non-TNL) within the gene family.
• Duplication and Selection Analysis:
  • Use MCScanX to identify tandem and segmental duplication events [2].
  • Calculate non-synonymous (Ka) and synonymous (Ks) substitution rates for orthologous gene pairs using tools like KaKs_Calculator 2.0 to infer selection pressures. A Ka/Ks > 1 indicates positive selection [2] [5].

Figure 1: A standard workflow for genome-wide identification and analysis of the NBS gene superfamily.

Application Note: Analyzing Genetic Variation in Tolerant Accessions

Understanding the genetic variation of NBS genes in disease-tolerant plant accessions is a key strategy for uncovering natural sources of resistance and guiding breeding programs.

Case Study: Cotton Leaf Curl Disease (CLCuD)

A comparative study between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions revealed substantial genetic variation in their NBS genes. The analysis identified 6,583 unique variants in the tolerant Mac7 accession compared to 5,173 in the susceptible Coker 312 [1]. This suggests that the tolerant line possesses a richer repertoire of polymorphisms, which may contribute to its enhanced disease resistance phenotype.

Protocol: Genetic Variation Analysis Pipeline

• Population Resequencing: Sequence a population of plant accessions, including both tolerant and susceptible varieties, using whole-genome resequencing.
• Variant Calling: Map reads to a reference genome and call SNPs and InDels using standard variant calling pipelines (e.g., GATK).
• Variant Annotation: Annotate identified variants based on their location (e.g., promoter, coding sequence, intron) and predicted effect (e.g., synonymous, non-synonymous) on NBS genes.
• Population Genetics Analysis: Calculate nucleotide diversity (Ï€) and fixation indices (FST) between tolerant and susceptible groups to identify regions of the genome, particularly within NBS genes, that show signatures of selection [5].
• Association Analysis: Correlate specific haplotypes or non-synonymous variants in NBS genes with the tolerance phenotype to pinpoint candidate causal polymorphisms.

Protocol: Functional Validation Using Virus-Induced Gene Silencing (VIGS)

To confirm the functional role of a candidate NBS gene in disease resistance, in planta validation is essential. VIGS is a powerful technique for transient gene silencing.

• Vector Construction: Clone a 300-500 bp fragment of the target NBS gene into a VIGS vector (e.g., based on Tobacco Rattle Virus).
• Plant Infiltration: At the 2-4 leaf stage, infiltrate tolerant and susceptible plants with the recombinant VIGS vector using Agrobacterium tumefaciens-mediated delivery. Include control plants infiltrated with an empty vector.
• Pathogen Challenge: Once silencing is established (typically 2-3 weeks post-infiltration), challenge the plants with the target pathogen.
• Phenotypic Assessment: Monitor and score disease symptoms over time. Compare the disease progression in silenced plants versus control plants.
• Molecular Confirmation:
  • Use qPCR to verify the knockdown efficiency of the target NBS gene in silenced plants.
  • Quantify pathogen biomass (e.g., through viral titer quantification) to assess the effect of silencing on resistance.

This approach was successfully used to validate the role of a candidate NBS gene (GaNBS from orthogroup OG2), where silencing in a resistant cotton background led to increased virus titers, demonstrating its putative role in virus resistance [1].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents and Tools for NBS Gene Research

Reagent/Tool	Function/Application	Example Sources/Platforms
Pfam HMM Profile PF00931	Core model for identifying the NBS domain in protein sequences.	Pfam Database
HMMER Software Suite	For performing sensitive homology searches using HMMs.	http://www.hmmer.org/ [2]
NCBI Conserved Domain Database (CDD)	For verifying and annotating conserved protein domains.	https://www.ncbi.nlm.nih.gov/cdd [2]
VIGS Vectors	For transient functional validation of NBS genes through silencing.	TRV-based vectors [1]
qPCR Reagents & Validated Primers	For gene expression validation and silencing efficiency confirmation.	SYBR Green kits; primers spanning exon-exon junctions [6]
Interaction Networks (e.g., PCNet)	For systems biology approaches and stratification analysis.	[7]
Fosimdesonide	Fosimdesonide, MF:C49H61N4O15PS, MW:1009.1 g/mol	Chemical Reagent
Tenacissoside C	Tenacissoside C, MF:C53H76O19, MW:1017.2 g/mol	Chemical Reagent

Figure 2: Simplified signaling mechanism of a typical NBS-LRR protein, showing the roles of its major domains in pathogen perception and defense activation.

The nucleotide-binding site leucine-rich repeat (NBS-LRR) genes represent one of the largest and most critical plant gene families, encoding intracellular immune receptors that perceive pathogen-derived effectors and initiate effector-triggered immunity (ETI) [8] [9]. These genes are broadly classified into two major subfamilies based on their N-terminal domains: Toll/interleukin-1 receptor (TIR) NBS-LRR (TNL) and coiled-coil (CC) NBS-LRR (CNL) genes [8] [10]. A third, smaller class contains a resistance to powdery mildew 8 (RPW8) domain (RNL) [9]. The architectural diversity of these genes—from their canonical domain structures to species-specific variations—plays a fundamental role in plant pathogen recognition and defense signaling. This diversity is not random but follows evolutionary patterns shaped by genomic duplication events and pathogen pressures, particularly evident in comparative analyses across plant families [9] [11]. Understanding these patterns is crucial for genetic variation analysis in disease-tolerant plant accessions, providing a foundation for identifying durable resistance genes for crop improvement.

Table 1: NBS-LRR Gene Classification and Characteristics

Class	N-Terminal Domain	Prevalence in Dicots/Monocots	Key Structural Features	Representative Genes
TNL	TIR (Toll/Interleukin-1 Receptor)	Abundant in dicots; rare in cereals [10] [12]	Adopts a flavodoxin-like fold; self-associates for signaling [8]	RPS4, RRS1, Ma, L6 [8] [12]
CNL	CC (Coiled-Coil)	Predominant in monocots; present in dicots [8] [10]	Largely helical structure; functional diversity in self-association and signaling [8]	RPS2, RPS5, Rp1-D21, MLA10 [8]
RNL	RPW8 (Resistance to Powdery Mildew 8)	Present in both dicots and monocots [9]	Functions downstream in signal transduction from TNLs/CNLs [9]	ADR1, NRG1 [9]

Classical Domain Architectures and Functional Mechanisms

The canonical structure of NBS-LRR genes comprises three core domains: a variable N-terminal domain (TIR, CC, or RPW8), a central nucleotide-binding site (NBS) domain, and a C-terminal leucine-rich repeat (LRR) domain [9]. The N-terminal domain is primarily involved in signal transduction, the NBS domain binds and hydrolyzes nucleotides and acts as a molecular switch, while the LRR domain is crucial for pathogen recognition specificity and is often under diversifying selection [8] [10].

TNL Signaling and Architecture

TNL proteins contain a TIR domain that adopts a conserved flavodoxin-like fold consisting of five α-helices surrounding a five-strand β-sheet [8]. Signal transduction is intimately linked to the TIR domain's ability to self-associate through specific interfaces formed by surface-exposed α-helices [8]. Upon pathogen perception, a conformational change occurs, facilitating TIR-TIR interaction and oligomerization, which triggers downstream defense signaling cascades leading to the hypersensitive response (HR) [12]. Many TNLs also feature a C-terminal region beyond the LRR domain, whose function is less characterized but may play roles in ligand binding or intramolecular interactions [12]. For instance, the Ma gene from peach possesses an unusually large C-terminal region with five duplicated post-LRR (PL) domains, longer than the rest of the protein combined [12].

CNL Diversity and Classification

CNLs are characterized by an N-terminal coiled-coil (CC) domain, a largely helical structure that exhibits significant functional diversity [8]. CC domains have been categorized into several classes based on specific motifs and functions:

• CCEDVID: Characterized by a conserved EDVID motif, suggested to be involved in intramolecular interactions with the NB domain (e.g., Sr33, MLA10, Rx) [8].
• CCR: Shares similarity with RPW8 and is often associated with broad-spectrum resistance [8].
• CCCAN: The "canonical" CC domain found in NLRs like RPS2 and RPS5 [8].
• SD-CC: Includes a large auxiliary solanaceous domain N-terminal to the CC domain (e.g., Sw-5b, Prf) [8].
• I2-like: A monophyletic clade with similarity to the tomato I2 CNL, also possessing an EDVID motif but distinct from CCEDVID class [8].

Unlike TIR domains, the overall structure of CC domains remains debated, complicating functional categorization. However, they are often involved in oligomerization and protein-protein interactions necessary for immune signaling [8].

Evolutionary Patterns and Species-Specific Expansion

Genome-wide comparative analyses reveal that NBS-LRR genes have undergone dynamic and species-specific evolutionary patterns across plant lineages, primarily driven by gene duplication and loss events [9] [11]. These patterns reflect adaptation to specific pathogen pressures and ecological niches.

Comparative Genomics in Rosaceae

A comprehensive study of five Rosaceae species—woodland strawberry (Fragaria vesca), apple (Malus × domestica), pear (Pyrus bretschneideri), peach (Prunus persica), and mei (Prunus mume)—revealed striking differences in NBS-LRR gene numbers, from 144 in strawberry to 748 in apple [11]. All species contained more non-TNLs (CNLs and XNLs) than TNLs, but the proportion of TNLs varied significantly, from 15.97% in strawberry to 47.12% in pear [11].

Table 2: NBS-LRR Gene Repertoire in Five Rosaceae Species

Species	Total NBS-LRR Genes	TNL Genes (%)	CNL Genes (%)	XNL Genes (%)	Multi-Gene Proportion (%)
F. vesca (Strawberry)	144	23 (15.97%)	89 (61.81%)	32 (22.22%)	32.64%
M. × domestica (Apple)	748	219 (29.28%)	446 (59.63%)	83 (11.10%)	68.98%
P. bretschneideri (Pear)	469	221 (47.12%)	201 (42.86%)	47 (10.02%)	63.33%
P. persica (Peach)	354	128 (36.16%)	193 (54.52%)	33 (9.32%)	65.82%
P. mume (Mei)	352	153 (43.47%)	173 (49.15%)	26 (7.39%)	53.98%

These differences reflect distinct evolutionary trajectories: woody perennial species (apple, pear, peach, mei) showed higher proportions of multi-gene families (exceeding 50% in apple, pear, and peach) compared to the herbaceous strawberry (32.64%), indicating more extensive recent duplications in woody species [11]. Phylogenetic analysis revealed 385 species-specific duplicate clades, with high percentages of NBS-LRR genes derived from species-specific duplication (e.g., 66.04% in apple, 48.61% in pear, 37.01% in peach) [11]. This suggests that recent, species-specific duplications have been the primary driver of NBS-LRR expansion in Rosaceae, rather than preservation of ancestral genes.

Evolutionary Dynamics and Selection Pressures

The evolution of NBS-LRR genes follows different patterns between gene classes and lineages. In Rosaceae, TNLs exhibited significantly higher Ks (synonymous substitution rate) and Ka/Ks (nonsynonymous/synonymous substitution rate) values compared to non-TNLs, with most NBS-LRRs having Ka/Ks ratios less than 1, indicating evolution under purifying selection [11]. However, the higher Ka/Ks for TNLs suggests they may be evolving more rapidly to adapt to different pathogens compared to non-TNLs [11].

Studies in broader plant contexts confirm these divergent patterns. In cereals, TIR-domain-containing NBS-LRR genes were not amplified during evolution, unlike in dicots, leading to a predominance of CNLs in grass species [10]. Some cereal lineages even evolved unique NBS-LRR classes, such as a 50-member group in rice that encodes proteins similar to N-termini and NBS domains but lacks LRR domains entirely [10].

Diagram 1: Evolutionary divergence of NBS-LRR genes between dicots and monocots, showing the independent expansion and contraction patterns that led to species-specific domain architectures. Created with DOT language.

Application Notes: Experimental Protocols for Genetic Variation Analysis

Genome-Wide Identification and Classification of NBS-LRR Genes

Purpose: To systematically identify and classify NBS-LRR genes in plant genomes for genetic variation studies.

Materials and Reagents:

• High-quality genomic DNA and annotation files
• HMMER software suite
• Pfam and NCBI-CDD databases
• Multiple sequence alignment tools (e.g., MAFFT, MUSCLE)
• Phylogenetic analysis software (e.g., RAxML, MrBayes)

Procedure:

• Sequence Retrieval: Download whole genome sequences and annotation files from relevant databases (e.g., Genome Database for Rosaceae for Rosaceae species) [9].
• Initial Screening: Perform BLAST and HMMER searches using the hidden Markov model of the NB-ARC domain (PF00931) as query with an E-value threshold of 1.0 for BLAST and default parameters for HMMER [9].
• Domain Validation: Submit candidate genes to online Pfam and NCBI-CDD search to confirm presence of N-terminal domains (CC/TIR/RPW8) and NBS domains using an E-value cutoff of 10⁻⁴ [9].
• Classification: Classify validated NBS-LRR genes into TNL, CNL, and RNL classes based on their N-terminal domain [9].
• Gene Structure Analysis: Use tools like GSDS2.0 to analyze intron-exon structures based on coding and genomic DNA sequences [9].
• Motif Identification: Identify conserved motifs using MEME suite with parameters set to discover 10 motifs [9].

Troubleshooting Tips:

• For fragmented genome assemblies, use iterative BLAST searches with known NBS-LRR sequences from related species.
• For domain verification, manually inspect borderline cases with E-values near the cutoff.
• For classification ambiguity, use multiple domain prediction tools and phylogenetic analysis.

Functional Analysis of CC Domains in CNLs

Purpose: To characterize the structure and function of CC domains in CNL proteins.

Materials and Reagents:

• Cloning vectors and host strains
• Site-directed mutagenesis kit
• Model plants (e.g., Nicotiana benthamiana)
• Agrobacterium tumefaciens strains for transient expression
• Cell death assay reagents
• Co-immunoprecipitation materials

Procedure:

• Construct Design: Design CC domain constructs through detailed sequence and secondary structure comparisons. Define domain boundaries based on sequence alignment with previously characterized CC domains [8].
• Cloning: Clone full-length CNL genes and CC domain fragments into appropriate expression vectors.
• Transient Expression: Express constructs in N. benthamiana leaves via Agrobacterium-mediated transformation [8].
• Cell Death Assay: Monitor hypersensitive response (HR) cell death symptoms over 2-5 days post-infiltration [8].
• Oligomerization Studies: Test self-association capability through co-immunoprecipitation and bimolecular fluorescence complementation assays [8].
• Structure-Function Analysis: Create targeted mutations in conserved motifs (e.g., EDVID) and assess functional impact [8].

Validation Methods:

• Quantitative cell death scoring using electrolyte leakage or Evans Blue staining
• Confocal microscopy for subcellular localization
• Western blotting to verify protein expression
• Yeast two-hybrid assays for protein-protein interactions

Diagram 2: Experimental workflow for genome-wide identification and classification of NBS-LRR genes, showing key bioinformatics steps from initial screening to detailed analysis. Created with DOT language.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Research Reagents and Resources for NBS-LRR Gene Analysis

Reagent/Resource	Function/Application	Example Use Cases	Key Features
HMMER Suite	Hidden Markov Model-based sequence analysis	Identification of NBS-LRR genes using PF00931 (NB-ARC) profile	Sensitive detection of divergent family members; handles domain architecture
Pfam Database	Protein family classification	Validation of N-terminal (TIR/CC/RPW8) and NBS domains	Curated multiple sequence alignments and HMM profiles; E-value statistics
MEME Suite	Motif discovery and analysis	Identification of conserved motifs in NBS domains	Discovers ungapped motifs; statistical significance assessment
Nicotiana benthamiana	Transient expression system	Functional assays for cell death signaling	Susceptible to Agrobacterium; high protein expression; rapid results
Agrobacterium tumefaciens	Plant transformation vector delivery	Transient expression of NBS-LRR constructs	Efficient DNA transfer to plant cells; compatible with binary vectors
Phylogenetic Software (RAxML)	Evolutionary relationship inference	Reconstruction of NBS-LRR gene families	Handles large datasets; maximum likelihood approach; bootstrap support
Ivermectin B1a monosaccharide	Ivermectin B1a monosaccharide, MF:C41H62O11, MW:730.9 g/mol	Chemical Reagent	Bench Chemicals
Cdki-IN-1	Cdki-IN-1, MF:C16H15ClN2O, MW:286.75 g/mol	Chemical Reagent	Bench Chemicals

The architectural diversity of NBS-LRR genes, from the classical TNL/CNL division to species-specific domain patterns, represents a sophisticated plant immune strategy shaped by evolutionary pressures. The dynamic expansion and contraction of these gene families across plant lineages, particularly the recent species-specific duplications observed in Rosaceae woody perennials, highlight the ongoing arms race between plants and their pathogens [9] [11]. The functional specialization of domain architectures, such as the PL domains in peach TNLs and the unique LRR-less variants in cereals, provides insights into the mechanistic diversity of pathogen recognition and signaling [10] [12]. For researchers investigating genetic variation in disease-tolerant accessions, the experimental frameworks and resources presented here offer practical pathways to identify and characterize novel resistance genes. This knowledge not only advances our understanding of plant immunity but also provides valuable tools for crop improvement through marker-assisted selection and genetic engineering of durable disease resistance.

This application note details the mechanisms behind the expansion of Nucleotide-Binding Site (NBS) gene families, which are the largest class of plant disease resistance (R) genes. Framed within broader research on genetic variation in tolerant plant accessions, this resource provides a comparative quantitative analysis of NBS genes across species and delineates the distinct contributions of whole-genome duplication (WGD) and tandem duplication to the evolution of this critical gene family. The protocols herein offer validated methodologies for identifying NBS genes and profiling duplication events, supported by specific reagent solutions and analytical workflows. This information is vital for researchers aiming to understand plant immunity and engineer durable, broad-spectrum disease resistance in crops.

NBS-LRR genes encode proteins that play a vital role in plant innate immunity by recognizing pathogen effector proteins and initiating defense responses [2]. In the arms race between plants and their pathogens, the expansion and diversification of the R-gene repertoire are essential for plant survival [13].

Research across diverse plant species has consistently revealed that two primary molecular forces are responsible for the expansion of NBS gene families:

• Tandem Duplication: The localized duplication of genes in close proximity on a chromosome, often leading to clusters of related NBS genes that are hotspots for the evolution of new pathogen specificities.
• Whole-Genome Duplication (WGD): The duplication of the entire genome, which provides raw genetic material for neofunctionalization or subfunctionalization of NBS genes.

The balance between these two forces varies by species, influencing the architecture and size of the R-gene repertoire. The following sections provide a detailed, data-driven breakdown of their roles.

Quantitative Analysis of NBS Expansion Across Species

Table 1: Documented NBS Gene Counts and Dominant Expansion Mechanisms in Plants

Species	Total NBS Genes	Dominant Expansion Mechanism(s)	Key Findings	Citation
Nicotiana tabacum (Tobacco)	603	Whole-Genome Duplication (WGD)	WGD contributed significantly to NBS expansion in this allotetraploid. The total count (~603) approximates the sum of its parental species.	[2]
Akebia trifoliata	73	Tandem & Dispersed Duplication	Tandem duplications produced 33 genes; dispersed duplications produced 29 genes. 64 mapped genes were unevenly distributed, with 41 located in clusters.	[13]
Cowpea (Vigna unguiculata)	2,188 R-genes	Tandem & Dispersed Duplication	Dispersed and tandem duplication events under purifying selection were the main contributors to kinome expansion.	[14]
Multiple Species (34 species)	12,820 NBS genes	Tandem Duplication	A large-scale comparative study identified several orthogroups and observed tandem duplications as a key driver of diversity.	[1]

Experimental Protocols for Profiling NBS Genes and Duplication Events

Protocol 1: Genome-Wide Identification and Classification of NBS Genes

This protocol is adapted from methods used in recent studies of Nicotiana and Akebia trifoliata genomes [13] [2].

• Objective: To comprehensively identify and classify all NBS-LRR genes within a plant genome.
• Principle: Using hidden Markov models (HMMs) and conserved domain databases to scan a genome's proteome for the characteristic NB-ARC domain and associated structural domains.

Workflow Steps:

• Data Acquisition: Download the genome assembly and annotated protein sequences for the target species.
• HMM Search: Perform an HMM search using HMMER software with the PFAM model PF00931 (NB-ARC domain) as the query. Use a liberal E-value cutoff (e.g., 1.0) to maximize initial candidate identification.
• Domain Verification and Classification:
  • Scan all candidate sequences against the Pfam database (e.g., using PfamScan) to verify the presence of the NBS domain and identify TIR (PF01582) and LRR (PF08191, PF00560, etc.) domains.
  • Identify Coiled-coil (CC) domains using the NCBI Conserved Domain Database (CDD) and/or tools like CoiledCoil with a threshold of 0.5.
  • Classify genes into subfamilies (TNL, CNL, RNL, N, etc.) based on their domain architecture.
• Redundancy Removal: Merge results from different searches and remove duplicate entries to generate a non-redundant set of NBS genes.

Protocol 2: Analysis of Gene Duplication Events

This protocol outlines the steps to distinguish between WGD and tandem duplication events in the identified NBS gene family [13] [2].

• Objective: To determine the mode of duplication (WGD/tandem/dispersed) for NBS genes and assess evolutionary pressures.
• Principle: Syntenic block analysis identifies WGD events, while physical clustering on chromosomes indicates tandem duplication. The ratio of non-synonymous to synonymous substitutions (Ka/Ks) indicates selection pressure.

Workflow Steps:

• Self-BLAST and MCScanX:
  • Perform an all-vs-all BLASTP search of the species' protein sequences.
  • Input the BLAST results into MCScanX to identify segmental (WGD) and tandem duplication events across the genome.
• Synteny Analysis:
  • For cross-species comparison, perform reciprocal BLASTP searches between the target and a related species.
  • Use MCScanX to identify syntenic blocks. NBS genes located in corresponding syntenic blocks are likely products of ancient WGDs.
• Tandem Duplication Identification:
  • Define tandem duplicates as NBS genes of the same subfamily located within a close physical distance (e.g., ≤ 1 gene) on the same chromosome [13].
• Evolutionary Pressure Analysis:
  • For pairs of duplicated genes, align coding sequences and calculate non-synonymous (Ka) and synonymous (Ks) substitution rates using tools like KaKs_Calculator 2.0.
  • Interpret results: A Ka/Ks ratio > 1 indicates positive selection; < 1 indicates purifying selection; ≈ 1 indicates neutral evolution.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Research Reagents, Tools, and Databases for NBS Gene Analysis

Item Name	Function / Application	Specification / Note
PF00931 HMM Profile	Core model for identifying the NB-ARC domain in HMMER searches.	Sourced from the Pfam database. Critical for initial gene discovery.
Pfam & NCBI CDD	Databases for verifying protein domains (TIR, LRR, CC, etc.) to classify NBS subfamilies.	Essential for accurate structural classification of candidate genes.
MCScanX	Software tool for analyzing genome collinearity and identifying WGD and tandem duplication events.	Key for delineating the evolutionary forces behind gene family expansion.
KaKs_Calculator	Software for calculating Ka and Ks substitution rates.	Determines the selection pressure acting on duplicated gene pairs.
NBS Profiling Primers	PCR primers targeting conserved NBS motifs (P-loop, Kinase-2, GLPL) for resistance gene enrichment.	Used for targeted sequencing and diversity analysis of R-genes in germplasm [15].
Lly-283	Lly-283, MF:C17H18N4O4, MW:342.35 g/mol	Chemical Reagent
Leesggglvqpggsmk tfa	Leesggglvqpggsmk tfa, MF:C66H109F3N18O26S, MW:1659.7 g/mol	Chemical Reagent

Concluding Perspectives

Understanding the dual roles of WGD and tandem duplication is fundamental for the strategic enhancement of disease resistance in crops. WGD provides a broad-scale foundation for R-gene repertoire expansion, while tandem duplication acts as a agile mechanism for rapid, localized diversification in response to pathogen pressure [13] [2]. The protocols and data summarized in this application note provide a roadmap for researchers to dissect the evolutionary history of NBS genes in tolerant accessions. This knowledge, in turn, informs modern breeding and biotechnological approaches, such as marker-assisted selection and genome editing, to develop crops with durable and resilient immune systems.

Nucleotide-binding site (NBS) genes constitute the largest family of plant disease resistance (R) genes, playing a crucial role in plant immunity. This application note provides a comprehensive comparative genomic analysis of the NBS-LRR gene family across diverse plant species, revealing substantial variation in gene repertoire size, subfamily composition, and evolutionary dynamics. We present standardized protocols for genome-wide identification and characterization of NBS genes, along with experimental frameworks for investigating their roles in disease resistance. The findings offer valuable insights for researchers investigating genetic variation in NBS genes of tolerant accessions and facilitate the development of disease-resistant crop varieties.

Plant resistance (R) genes encoding nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains are crucial components of the plant immune system, responsible for detecting pathogen effectors and initiating defense responses [16] [17]. Approximately 80% of functionally characterized R genes belong to the NBS-LRR (or NLR) gene family, making them a major focus of plant immunity research [17]. The NBS domain functions as a molecular switch by binding and hydrolyzing ATP/GTP to activate downstream immune signaling, while the LRR domain is responsible for recognizing diverse pathogen effectors [4] [17].

With advances in sequencing technologies, genome-wide analyses of NBS-LRR genes have been performed across numerous plant species, revealing remarkable diversity in gene number, subfamily composition, and genomic distribution. This application note synthesizes current methodologies and findings from comparative genomic studies of NBS gene families, providing standardized protocols for researchers investigating genetic variation in NBS genes of tolerant accessions.

Results and Data Analysis

Comparative Analysis of NBS Gene Repertoire Across Species

Table 1: NBS-LRR Gene Family Size Variation Across Plant Species

Species	Genome Type	Total NBS Genes	CNL	TNL	RNL	Other/Partial	Reference
Nicotiana tabacum	Allotetraploid	603	224 (37.1%)	73 (12.1%)	-	306 (50.7%)	[2]
Nicotiana sylvestris	Diploid	344	130 (37.8%)	42 (12.2%)	-	172 (50.0%)	[2]
Nicotiana tomentosiformis	Diploid	279	112 (40.1%)	40 (14.3%)	-	127 (45.5%)	[2]
Akebia trifoliata	Diploid	73	50 (68.5%)	19 (26.0%)	4 (5.5%)	-	[13]
Raphanus sativus	Diploid	225	51 (22.7%)	134 (59.6%)	0 (0%)	40 (17.8%)	[18]
Salvia miltiorrhiza	Diploid	196	75 (38.3%)	2 (1.0%)	1 (0.5%)	118 (60.2%)	[17]
Nicotiana benthamiana	Diploid	156	25 (16.0%)	5 (3.2%)	-	126 (80.8%)	[4]
Arabidopsis thaliana	Diploid	164	-	-	-	-	[18]
Brassica oleracea	Diploid	244	-	-	-	-	[18]
Saccharum spontaneum	Polyploid	691	-	-	-	-	[16]

Table 2: NBS-LRR Subfamily Classification Based on Domain Architecture

Subfamily	N-Terminal Domain	Central Domain	C-Terminal Domain	Primary Function
CNL	Coiled-coil (CC)	NBS (NB-ARC)	LRR	Pathogen recognition and immunity activation
TNL	TIR (Toll/Interleukin-1 Receptor)	NBS (NB-ARC)	LRR	Pathogen recognition and immunity activation
RNL	RPW8 (Resistance to Powdery Mildew 8)	NBS (NB-ARC)	LRR	Downstream defense signaling
CN	Coiled-coil (CC)	NBS (NB-ARC)	-	Regulatory or adapter functions
TN	TIR	NBS (NB-ARC)	-	Regulatory or adapter functions
NL	-	NBS (NB-ARC)	LRR	Atypical recognition
N	-	NBS (NB-ARC)	-	Regulatory functions

The NBS gene repertoire varies substantially across plant species, influenced by factors such as genome size, ploidy, life history, and evolutionary pressure from pathogens. Comparative analysis reveals that whole genome duplication (WGD), gene expansion, and allele loss significantly impact NBS-LRR gene numbers [16]. For instance, the allotetraploid Nicotiana tabacum contains approximately the combined total of NBS genes from its diploid progenitors (N. sylvestris and N. tomentosiformis), suggesting that polyploidization contributes to NBS repertoire expansion [2].

Subfamily composition also exhibits remarkable variation. Monocot species like rice (Oryza sativa) have completely lost TNL and RNL subfamilies, while gymnosperms such as Pinus taeda show significant expansion of TNL genes (89.3% of typical NBS-LRRs) [17]. In Salvia miltiorrhiza, researchers observed a notable reduction in TNL (1.0%) and RNL (0.5%) subfamilies, with CNL genes dominating the NBS repertoire (38.3%) [17].

Genomic Distribution and Organization

NBS-LRR genes are frequently distributed unevenly across chromosomes, often clustered at chromosome ends [13]. In Akebia trifoliata, 64 mapped NBS candidates were unevenly distributed on 14 chromosomes, with most assigned to chromosome ends [13]. Similarly, in radish (Raphanus sativus), 72% of NBS-encoding genes were grouped in 48 clusters distributed in 24 crucifer blocks on chromosomes [18].

Cluster analysis reveals that NBS genes are frequently organized in complex arrays containing both complete and incomplete genes. This arrangement serves as a source of variation and a plant's reservoir for producing new functional R alleles through frameshift recombination and DNA repair processes [15].

Figure 1: Workflow for Genome-wide Identification and Analysis of NBS-LRR Genes

Experimental Protocols

Protocol 1: Genome-Wide Identification of NBS-Encoding Genes

Materials and Reagents

• High-quality genome assembly and annotation files
• High-performance computing resources
• Sequence analysis software (HMMER, BLAST, InterProScan)
• Programming environment (Python/R) for data processing

Step-by-Step Procedure

Step 1: HMM-Based Candidate Identification

• Download the NB-ARC domain HMM profile (PF00931) from the Pfam database
• Perform HMM search against the target proteome using HMMER v3.1b2 or later
• Use an expectation value (E-value) cutoff of < 1×10⁻²⁰ for initial screening [4]
• Extract candidate protein sequences for further verification

Step 2: Domain Verification and Classification

• Verify the presence of NBS domains using Pfam database (E-value < 0.01) [4]
• Identify additional domains (TIR, CC, RPW8, LRR) using:
  • NCBI Conserved Domain Database (CDD)
  • SMART tool (http://smart.embl-heidelberg.de/)
  • Coiled-coil domains predicted with Coiledcoil (threshold 0.5) [13]
• Classify genes into subfamilies (CNL, TNL, RNL, CN, TN, NL, N) based on domain architecture

Step 3: Manual Curation and Non-redundant Set Generation

• Remove redundant genes and pseudogenes
• Verify domain organization through multiple databases
• Generate a non-redundant set of NBS-encoding genes

Protocol 2: Evolutionary and Phylogenetic Analysis

Materials and Reagents

• Multiple sequence alignment software (Clustal W, MUSCLE)
• Phylogenetic analysis package (MEGA11, IQ-TREE)
• Synteny analysis tool (MCScanX)
• Selection pressure analysis software (KaKs_Calculator)

Step-by-Step Procedure

Step 1: Sequence Alignment and Phylogenetic Reconstruction

• Perform multiple sequence alignment of NBS protein sequences using MUSCLE v3.8.31 with default parameters [2]
• Construct phylogenetic trees using maximum likelihood method in MEGA11 [2]
• Use bootstrap analysis with 1000 replicates to assess node support [4]
• Classify genes into clades based on phylogenetic relationships

Step 2: Synteny and Duplication Analysis

• Perform self-BLASTP based on protein sequences to identify paralogs [2]
• Identify segmental and tandem duplications using MCScanX under default configurations [2]
• Calculate non-synonymous (Ka) and synonymous (Ks) substitution rates using KaKs_Calculator 2.0 [2]
• Interpret evolutionary dynamics based on Ka/Ks ratios

Step 3: Selection Pressure Analysis

• Identify orthologous gene pairs across related species
• Calculate Ka/Ks ratios for each orthologous pair
• Classify genes under positive (Ka/Ks > 1), purifying (Ka/Ks < 1), or neutral (Ka/Ks = 1) selection
• Correlate selection patterns with gene function and expression

Protocol 3: Expression Analysis Under Pathogen Stress

Materials and Reagents

• Plant materials with different disease resistance phenotypes
• Pathogen isolates (e.g., Fusarium oxysporum for radish) [18]
• RNA extraction kit
• RNA-seq library preparation kit
• Quantitative PCR reagents

Step-by-Step Procedure

Step 1: Experimental Design and Sample Collection

• Establish disease resistance assays using appropriate pathogens
• Collect tissue samples at multiple time points post-inoculation
• Include biological replicates (minimum three per treatment)
• Preserve samples in RNA stabilization reagent

Step 2: Transcriptome Sequencing and Analysis

• Extract total RNA using standardized protocols
• Prepare RNA-seq libraries following manufacturer's instructions
• Sequence libraries on Illumina platform (minimum 20 million reads per sample)
• Map cleaned reads to reference genome using HISAT2 [2]
• Calculate expression levels (FPKM) using Cufflinks v2.2.1 [2]
• Identify differentially expressed genes using Cuffdiff [2]

Step 3: Validation and Functional Correlation

• Select candidate NBS genes based on expression patterns
• Design gene-specific primers for qRT-PCR validation
• Perform qRT-PCR using SYBR Green chemistry
• Correlate expression patterns with disease resistance phenotypes
• Identify key regulatory NBS genes for further functional studies

Figure 2: NBS-LRR Mediated Disease Resistance Signaling Pathway

The Scientist's Toolkit

Table 3: Essential Research Reagents and Computational Tools for NBS Gene Analysis

Category	Item/Software	Specific Function	Application Notes
Domain Databases	Pfam (PF00931)	NB-ARC domain identification	Primary HMM profile for initial screening
	NCBI CDD	Multiple domain verification	Confirms NBS and other associated domains
	SMART tool	Domain architecture analysis	Identifies complete domain organization
Analysis Software	HMMER v3.1b2	Hidden Markov Model searches	Core tool for identifying NBS domain-containing proteins
	MEME Suite	Conserved motif detection	Identifies conserved motifs within NBS domains
	MCScanX	Synteny and duplication analysis	Detects segmental and tandem duplications
	MEGA11	Phylogenetic reconstruction	Constructs evolutionary relationships
	KaKs_Calculator	Selection pressure analysis	Calculates Ka/Ks ratios for evolutionary analysis
Experimental Tools	HISAT2	RNA-seq read alignment	Maps transcriptomic data to reference genomes
	Cufflinks/Cuffdiff	Expression quantification	Identifies differentially expressed genes
	PlantCARE	Cis-element prediction	Analyzes promoter regions for regulatory elements
Validation Reagents	SYBR Green	qRT-PCR quantification	Validates RNA-seq expression patterns
	Pathogen isolates	Disease resistance assays	Tests functional relevance of candidate NBS genes
Semustine	Semustine, CAS:33073-59-5, MF:C10H18ClN3O2, MW:247.72 g/mol	Chemical Reagent	Bench Chemicals
Necroptosis-IN-5	Necroptosis-IN-5, MF:C30H31FN6O3, MW:542.6 g/mol	Chemical Reagent	Bench Chemicals

Discussion and Application

The variation in NBS gene repertoire across species reflects diverse evolutionary paths and adaptation to distinct pathogen pressures. Studies in sugarcane revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern cultivars, with the proportion significantly higher than expected, indicating S. spontaneum's greater contribution to disease resistance [16]. This finding demonstrates how comparative genomics can identify valuable genetic resources for breeding programs.

Evolutionary analyses indicate that both tandem and dispersed duplications are major forces responsible for NBS expansion. In Akebia trifoliata, tandem and dispersed duplications produced 33 and 29 NBS genes respectively [13]. This expansion mechanism allows plants to rapidly diversify their recognition capabilities against evolving pathogens.

Expression analyses across multiple species reveal that NBS genes are generally expressed at low levels, with a subset showing induced expression during pathogen challenge. In radish, 75 NBS-encoding genes contributed to resistance against Fusarium wilt, with specific genes like RsTNL03 and RsTNL09 showing positive regulation of resistance, while RsTNL06 appeared to function as a negative regulator [18]. This highlights the complex regulatory networks governing NBS-mediated immunity.

The protocols and analyses presented here provide a framework for investigating NBS gene variation in tolerant accessions, enabling researchers to identify key genetic determinants of disease resistance across crop species.

This application note provides comprehensive protocols for comparative analysis of NBS gene repertoires across plant species, highlighting the substantial variation in gene number, subfamily composition, and evolutionary dynamics. The standardized methodologies for genome-wide identification, evolutionary analysis, and expression profiling enable systematic investigation of NBS genes in the context of disease resistance. These approaches facilitate the identification of candidate R genes for molecular breeding and provide insights into the evolutionary mechanisms shaping plant immune systems. The integration of genomic and transcriptomic analyses outlined in this document offers a powerful strategy for elucidating the genetic basis of disease resistance in tolerant accessions.

Within the broader context of analyzing genetic variation in Nucleotide-Binding Site Leucine-Rich Repeat (NBS-LRR) genes from tolerant plant accessions, orthogroup analysis serves as a powerful computational framework for identifying sets of orthologous genes across multiple species. This protocol details the application of orthogroup analysis to delineate core genes conserved across lineages from species-specific genes that may underlie unique resistance capabilities. The NBS-LRR gene family constitutes the largest class of plant resistance (R) genes, with approximately 80% of characterized R genes encoding proteins containing these domains, which are essential for pathogen recognition and defense activation [19]. In pepper (Capsicum annuum L.), for instance, 252 NBS-LRR resistance genes have been identified, distributed unevenly across all chromosomes, with 54% forming 47 gene clusters [19]. The methodology outlined herein enables researchers to systematically characterize the genetic architecture of disease resistance, providing insights into evolutionary patterns and identifying potential genetic targets for crop improvement, particularly for enhancing resistance to devastating pathogens such as Ralstonia solanacearum, which causes bacterial wilt disease [20].

Key Concepts and Definitions

• Orthogroup: A set of genes descended from a single gene in the last common ancestor of the species being compared. Orthogroups represent gene families and can contain both orthologs and paralogs.
• Core Resistance Gene Set: Orthogroups present in all investigated species, representing conserved resistance mechanisms.
• Species-Specific Resistance Gene Set: Orthogroups found only in a single species or a subset of species, potentially conferring specialized resistance traits.
• NBS-LRR Genes: Genes encoding nucleotide-binding site and leucine-rich repeat proteins, which are central components of the plant immune system [19].
• TNL and nTNL Subfamilies: Two major NBS-LRR subclasses differentiated by their N-terminal domains: TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (nTNL), which includes CC-NBS-LRR (CNL) genes [19].

Experimental Protocol

Protein Sequence Collection and Curation

Purpose: To compile a comprehensive dataset of NBS-LRR protein sequences from genomes of interest for subsequent orthogroup inference.

Materials:

• Genome assemblies and corresponding annotation files (GFF/GTF format) for all species in the analysis
• High-performance computing resources for large-scale sequence analysis
• Bioinformatics tools: OrthoFinder, OrthoMCL, or similar orthogroup inference software

Procedure:

• Retrieve Sequences: Extract all protein-coding sequences from the annotated genomes of your target species. For example, in a study of Solanaceae resistance, this would include species such as Solanum dulcamara, Solanum nigrum, Capsicum annuum, and other relevant taxa [20].
• Identify NBS-LRR Candidates: a. Perform HMMER searches against the Pfam database using NBS-LRR-associated hidden Markov models (e.g., PF00931 for NBS domain) [19]. b. Conduct BLASTP searches using known NBS-LRR protein sequences as queries with an E-value cutoff of 1e-5. c. Combine results and remove redundant sequences using CD-HIT with 90% identity threshold.
• Validate Domain Architecture: a. Confirm the presence of characteristic NBS and LRR domains using Pfam scan [19]. b. Identify N-terminal domains (CC or TIR) using the COILS program and Pfam TIR models [19]. c. Classify sequences into nTNL (including CNL) and TNL subfamilies based on domain composition.

Table 1: Example NBS-LRR Gene Distribution Across Species

Species	Total NBS-LRR Genes	nTNL Genes	TNL Genes	Genes in Clusters	Reference
Capsicum annuum (pepper)	252	248	4	136 (54%)	[19]
Solanum dulcamara (bittersweet nightshade)	Data not specified in results; resistant to R. solanacearum	-	-	-	[20]
Solanum nigrum (black nightshade)	Data not specified in results; susceptible to R. solanacearum	-	-	-	[20]

Orthogroup Inference and Analysis

Purpose: To cluster NBS-LRR protein sequences into orthogroups and identify core and species-specific sets.

Procedure:

• Run Orthogroup Inference: a. Execute OrthoFinder with default parameters using the curated NBS-LRR protein sequences from all species. b. Alternatively, use OrthoMCL with an inflation parameter of 1.5 for more conservative clustering.
• Analyze Results: a. Core Orthogroups: Extract orthogroups present in all analyzed species. b. Species-Specific Orthogroups: Identify orthogroups unique to individual species. c. Lineage-Specific Orthogroups: Identify orthogroups shared only among related species (e.g., within a genus).
• Functional Annotation: a. Annotate orthogroups with gene ontology terms using InterProScan. b. Identify enriched functions in species-specific orthogroups using Fisher's exact test with FDR correction.

Table 2: Orthogroup Analysis Results from a Hypothetical Solanaceae Study

Orthogroup Category	Number of Orthogroups	Representative Genes	Potential Functional Significance
Core Orthogroups	15	Capana03g004459, Sdul02467, Snig12345	Conserved pathogen recognition mechanisms
S. dulcamara-Specific	8	Sdul13842, Sdul20916	Potential bacterial wilt resistance factors
C. annuum-Specific	12	Capana12g002134, Capana03g003876	Specialized recognition of pepper pathogens
S. nigrum-Specific	5	Snig08763, Snig15432	Unknown susceptibility factors
TNL-Enriched	6	Capana09g001234, Sdul_18765	TIR-mediated signaling pathways

Evolutionary and Structural Analyses

Purpose: To investigate evolutionary relationships and structural features within and between orthogroups.

Procedure:

• Phylogenetic Analysis: a. Select representative sequences from significant orthogroups. b. Perform multiple sequence alignment using MAFFT with L-INS-i strategy. c. Construct maximum-likelihood trees using IQ-TREE with 1000 ultrafast bootstraps.
• Motif Analysis: a. Identify conserved protein motifs using MEME suite with maximum 15 motifs. b. Compare motif compositions between core and species-specific orthogroups.
• Genomic Distribution: a. Map orthogroup members to their genomic positions using BEDTools. b. Identify gene clusters as genomic regions with ≥3 NBS-LRR genes within 200 kb. c. Visualize distribution using MapDraw or similar tools [19].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Materials and Reagents for Orthogroup Analysis in NBS-LRR Research

Item	Function/Application	Example/Specifications
Oxford Nanopore Technologies (ONT)	Long-read sequencing for genome assembly	Enables high-contiguity genome assemblies for accurate gene annotation [20]
Illumina Sequencing	Short-read sequencing for error correction	Paired-end reads (PE150) provide accurate base-level resolution [20]
OrthoFinder Software	Orthogroup inference from protein sequences	Computationally efficient method for identifying orthogroups across multiple species
HMMER Suite	Protein domain identification	Pfam hidden Markov models for NBS, LRR, TIR, and CC domains [19]
MEME Suite	Protein motif discovery	Identifies conserved motifs within NBS-LRR proteins (e.g., P-loop, RNBS-A, kinase-2) [19]
Ralstonia solanacearum Strains	Pathogen challenge experiments	Race 3 Biovar 2 strain UW551 for assessing resistance phenotypes [20]
Bay-707	Bay-707, MF:C15H20N4O2, MW:288.34 g/mol	Chemical Reagent
Biotin-16-dCTP	Biotin-16-dCTP, MF:C32H51N8O17P3S, MW:944.8 g/mol	Chemical Reagent

Workflow and Pathway Visualizations

Orthogroup Analysis Workflow for NBS-LRR Genes

NBS-LRR Gene Domain Architecture

Applications and Implications

Orthogroup analysis of NBS-LRR genes enables the identification of evolutionary patterns and functional specialization within plant immune systems. Studies in pepper have demonstrated a striking dominance of the nTNL subfamily (248 genes) over the TNL subfamily (only 4 genes), reflecting lineage-specific adaptations and evolutionary pressures [19]. The clustering of these genes in specific genomic regions—54% of pepper NBS-LRR genes form 47 physical clusters—highlights the dynamic evolution of resistance genes through mechanisms such as tandem duplications and genomic rearrangements [19].

Comparative analyses across Solanaceae species, including resistant accessions like Solanum dulcamara and susceptible ones like Solanum nigrum, can reveal orthogroups associated with bacterial wilt resistance [20]. The unique genes identified in resistant S. dulcamara with functions related to auxin transport, along with specific pattern recognition receptors (PRRs) present only in resistant species, provide promising candidates for further functional characterization [20]. These findings can directly inform breeding strategies for disease-resistant crops by identifying key genetic components that could be transferred to susceptible crop varieties.

Orthogroup analysis represents a powerful systematic approach for deciphering the complex landscape of plant resistance genes. By identifying both conserved and lineage-specific NBS-LRR genes, researchers can prioritize candidates for functional validation and potential integration into crop breeding programs. The protocols and methodologies outlined here provide a framework for conducting such analyses specifically within the context of genetic variation studies in tolerant plant accessions. As genomic technologies continue to advance, enabling more comprehensive genome assemblies and annotations [20], orthogroup analysis will become increasingly refined, offering deeper insights into the evolution of plant immunity and accelerating the development of durable disease resistance in agricultural crops.

Methodologies for Identifying NBS Genes and Profiling Genetic Variation

The genome-wide identification of specific gene families is a cornerstone of modern genomics, enabling researchers to understand the genetic basis of traits, including disease resistance. This process is particularly crucial for studying Nucleotide-Binding Site (NBS) genes, which constitute the largest family of plant disease resistance (R) genes. The identification and characterization of these genes in tolerant accessions provides invaluable insights into plant immunity mechanisms and offers potential genetic targets for crop improvement. This Application Note details standardized protocols for employing HMMER-based profile hidden Markov models (HMMs) and domain analysis pipelines for comprehensive NBS gene identification, with a specific focus on analyzing genetic variation between susceptible and tolerant plant accessions.

The following table catalogues the key computational tools and databases essential for executing a successful genome-wide identification project.

Table 1: Key Research Reagent Solutions for HMMER and Domain Analysis

Item Name	Function/Application	Specifications/Examples
HMMER Software Suite	Performs sequence similarity searches using probabilistic methods. Core programs include hmmscan, hmmsearch, and phmmer. [21]	Available from http://eddylab.org/software/hmmer/. Critical for building and searching with profile HMMs. [22]
Pfam Database	A curated collection of protein families and domains, each represented by multiple sequence alignments and profile HMMs. [23]	The NB-ARC domain (Pfam accession: PF00931) is the definitive model for identifying NBS genes. [13]
NCBI Conserved Domain Database (CDD)	Used for the functional classification and validation of domains present in identified candidate genes. [13]	Identifies TIR (PF01582), RPW8 (PF05659), and LRR (PF08191) domains for NBS gene sub-classification.
MEME Suite	Discovers conserved motifs within protein sequences, aiding in the structural characterization of gene families. [13]	Typically configured to identify ~10 motifs with widths of 6-50 amino acids.
OrthoFinder	Infers orthogroups and gene families across multiple species, providing insights into evolutionary relationships. [24]	Uses DIAMOND for fast sequence similarity searches and MCL for clustering.

Core Methodological Framework

Theoretical Foundation: HMMER and Domain Analysis

Profile Hidden Markov Models (HMMs) are probabilistic models that capture the consensus of a multiple sequence alignment, including position-specific conservation and variation. The HMMER software suite implements this technology for sensitive homology detection, functioning as a core engine in many bioinformatics pipelines [21].

For genome-wide studies, the typical workflow involves searching a whole proteome against a domain database like Pfam using hmmscan to identify proteins containing a domain of interest. The NB-ARC (NBS) domain is the universal signature for plant NBS-LRR resistance genes [13]. Subsequent domain analysis classifies identified genes into subfamilies (e.g., CNL, TNL, RNL) based on their N-terminal and C-terminal domains, providing the initial structural framework for understanding their potential function [13] [25].

Detailed Protocol: Genome-Wide Identification of NBS Genes

This protocol is adapted from established methodologies used in recent studies of NBS genes across various plant species [24] [13] [25].

Step 1: Data Acquisition

• Download the complete proteome file (all predicted protein sequences in FASTA format) for the target organism(s) from a public database like NCBI, Phytozome, or Plaza.

Step 2: HMMER-based Identification of NBS-Encoding Genes

• Obtain the HMM Profile: Download the NB-ARC domain HMM profile (PF00931) from the Pfam database.
• Perform Domain Scan: Use the hmmscan program from the HMMER suite to search the proteome against the NB-ARC profile.

• Process Results: Parse the hmmscan output to extract sequences with significant hits to the NB-ARC domain. Remove redundant sequences.

Step 3: Validation and Classification with Domain Analysis

• Verify NBS Domain: Re-analyze the candidate sequences against the full Pfam database to confirm the presence of the NB-ARC domain and identify other associated domains.
• Classify into Subfamilies:
  • Use NCBI's CDD search to identify TIR, RPW8, and LRR domains.
  • Use a tool like Coiledcoil (with a threshold of 0.5) to predict Coiled-Coil (CC) domains, as they are not always detected by Pfam [13].
  • Classify genes based on domain architecture:
    • CNL: CC-NBS-LRR
    • TNL: TIR-NBS-LRR
    • RNL: RPW8-NBS-LRR
    • NL: NBS-LRR (no clear N-terminal domain)
    • CN/TN: CC-NBS / TIR-NBS (lacking LRR) [25]

Step 4: Advanced Structural and Evolutionary Analysis

• Motif Discovery: Use the MEME suite to identify conserved motifs within the NBS domains of the identified genes. A standard setup uses 10 motifs with widths between 6-50 amino acids [13].
• Orthogroup Analysis: Use OrthoFinder with identified NBS proteins from multiple species to infer orthogroups. This helps identify evolutionarily conserved "core" NBS genes and species-specific expansions [24].
• Genetic Variation Analysis (For Tolerant vs. Susceptible Accessions): Compare the sequences, expression profiles, and haplotypes of candidate NBS genes between tolerant and susceptible accessions. Look for unique variants (SNPs, indels) and differential expression under stress conditions to pinpoint candidates responsible for tolerance [24].

Workflow Visualization

The following diagram illustrates the logical flow of the genome-wide identification and analysis pipeline.

Diagram Title: Genome-Wide NBS Gene Identification Pipeline

Applications in Genetic Variation Analysis of NBS Genes

The pipeline described above is instrumental in framing research on the genetic variation underlying disease tolerance. A comparative study of Gossypium hirsutum accessions—Mac7 (tolerant) and Coker312 (susceptible) to cotton leaf curl disease—identified 6,583 and 5,173 unique variants, respectively, in their NBS genes. This suggests that specific haplotypes in the tolerant accession may be key to its resistance [24]. Similarly, research on tung tree Fusarium wilt resistance compared the susceptible Vernicia fordii (90 NBS genes) and the resistant Vernicia montana (149 NBS genes). This analysis revealed that a specific NBS-LRR gene, Vm019719, was activated in the resistant species and confirmed via functional validation to confer resistance [25].

Table 2: Key Findings from NBS Gene Variation Studies in Tolerant Accessions

Study System	Tolerant Accession	Key Finding	Methodological Insight
Cotton Leaf Curl Disease [24]	G. hirsutum (Mac7)	6,583 unique NBS gene variants in the tolerant Mac7 vs. 5,173 in susceptible Coker312.	Haplotype analysis of NBS genes can pinpoint superior alleles for disease tolerance.
Tung Tree Fusarium Wilt [25]	Vernicia montana	149 NBS genes identified; Vm019719 (ortholog of Vf11G0978) confers resistance.	Comparative genomics between resistant and susceptible genotypes identifies causal R-genes.
Rice Low-Nitrogen Tolerance [26]	MAGIC Lines & Germplasm	Superior haplotype of LOC_Os06g06440 associated with higher relative grain yield under low nitrogen.	GWAS combined with haplotype analysis detects advantaged alleles for abiotic stress tolerance.

Troubleshooting and Best Practices

• Handling Frameshifts in Sequencing Data: When working with data from sequencing platforms prone to indel errors (e.g., pyrosequencing), standard HMMER may produce marginal alignments. For small-scale projects, consider using HMM-FRAME, a tool that incorporates sequencing error models to correct frameshifts during domain classification [23].
• Improving Specificity with Score Dissection: To better distinguish true homologous hits from spurious alignments driven by non-globular sequence segments (e.g., low-complexity regions), the dissectHMMER framework can be applied. It dissects the HMMER alignment score into fold-critical and remnant contributions, providing a more statistically robust E-value for homology inference [27].
• Ensuring Comprehensive Domain Annotation: Relying on a single tool for domain prediction, especially for CC domains, can lead to misclassification. Always use a combination of Pfam, CDD, and a dedicated coiled-coil prediction tool like Coiledcoil for accurate subfamily classification of NBS genes [13].

Leveraging Genome-Wide Association Studies (GWAS) to Link NBS Loci to Tolerance

Genome-wide association studies (GWAS) have emerged as a powerful method for identifying genetic variants associated with complex traits, including tolerance to various abiotic and biotic stresses in crops. A key class of genes involved in plant stress response and disease resistance is the Nucleotide-Binding Site Leucine-Rich Repeat (NBS-LRR) gene family. These genes represent a major component of the plant immune system, encoding proteins that recognize pathogen effectors and trigger defense responses [28].

This protocol details how to leverage GWAS to connect specific NBS loci with tolerance mechanisms, enabling researchers to identify candidate NBS genes underlying stress resilience. The methodology is framed within the broader context of investigating genetic variation in NBS genes across tolerant accessions, providing a robust framework for dissecting the genetic architecture of tolerance traits.

Key Principles and Background

The Role of NBS-LRR Genes in Plant Tolerance

NBS-LRR genes are part of the plant's innate immune system, responsible for detecting pathogens and initiating defense signaling cascades. They are characterized by a conserved nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) domain. The LRR domain is involved in pathogen recognition, while the NBS domain is responsible for ATP/GTP binding and activation of downstream signaling [28]. These genes often reside in clusters within plant genomes and can exhibit significant copy number variation, which contributes to evolutionary innovation in pathogen recognition.

GWAS Fundamentals for Gene Discovery

GWAS leverages historical recombination events in natural populations to identify statistical associations between genetic markers (typically SNPs) and phenotypic variation. Compared to traditional biparental QTL mapping, GWAS offers higher mapping resolution and enables the examination of a broader range of genetic diversity [29]. For tolerance traits, which are often quantitatively inherited, GWAS can detect loci with both major and minor effects, providing a more comprehensive view of the genetic architecture.

Experimental Design and Workflow

The following diagram illustrates the comprehensive workflow for linking NBS loci to tolerance traits using GWAS:

Figure 1: Comprehensive workflow for linking NBS loci to tolerance traits using GWAS.

Population Selection

The first critical step involves assembling a diverse germplasm panel with sufficient genetic variation to ensure adequate statistical power for association mapping. For crop species, this typically involves:

• Panel Size: 200-1000 accessions, depending on genetic architecture and marker density [30] [31].
• Diversity Sources: Landraces, wild relatives, and improved cultivars to capture a wide allele spectrum [28] [30].
• Population Structure Assessment: Use principal component analysis (PCA) and software like ADMIXTURE to quantify and account for population stratification [30] [32].

For example, a study on Korean landrace soybeans utilized 1,693 accessions to ensure sufficient genetic diversity for detecting associations with agronomic traits [30].

High-Throughput Phenotyping

Precise phenotyping for tolerance traits is essential for robust association mapping. Key considerations include:

• Trait Standardization: Define clear, measurable tolerance indices (e.g., relative root length, survival rate, disease severity score).
• Replication: Implement completely randomized designs with 2-3 biological replications to account for environmental variation [33] [34].
• Multi-Environment Trials: When possible, conduct phenotyping across multiple locations and seasons to account for genotype × environment interactions [31].

In rice salinity tolerance studies, researchers commonly use relative values (stress/control conditions) for root length, root dry weight, and other growth parameters to quantify tolerance [34].

Genotyping and Quality Control

High-density genotyping forms the foundation of GWAS. The protocol includes:

• Platform Selection: Utilize high-density SNP arrays (e.g., 7K SNP array in rice [33]) or sequencing-based approaches (GBS, RAD-seq, WGS).
• Quality Control: Apply filters for minor allele frequency (MAF > 0.05), missing data (<20%), and Hardy-Weinberg equilibrium [34] [31].
• Imputation: Use reference panels to impute missing genotypes and increase marker density.

For instance, a soybean flooding tolerance study utilized 34,718 high-quality SNPs after rigorous filtering of a 42,449 SNP marker set [31].

GWAS Statistical Analysis

Appropriate statistical models are crucial for reliable association detection:

• Model Selection: Implement mixed linear models (MLM) that account for population structure (Q matrix) and familial relatedness (K matrix) to reduce false positives [35] [34].
• Multiple Testing Correction: Apply false discovery rate (FDR) control or Bonferroni correction based on the number of independent tests.
• Multi-Locus Methods: Consider using multi-locus models (FarmCPU, BLINK) that can improve power for detecting small-effect loci [30] [31].

In a large-scale GWAS of risk tolerance in humans, researchers used LD Score regression to distinguish true polygenic signals from confounding biases such as population stratification [35].

NBS Loci Identification and Integration

The specific workflow for NBS gene discovery involves:

• Genome Annotation: Use tools like InterProScan or Pfam to identify NBS-domain containing genes in the reference genome.
• Physical Position Mapping: Cross-reference significant GWAS hits with the physical positions of annotated NBS genes.
• LD-Based Candidate Region Definition: Define candidate intervals based on linkage disequilibrium decay patterns (e.g., 309 kb in soybean where r² drops to 0.2 [30]).

A sorghum rust resistance study successfully identified a 57 kbp genomic region on chromosome 8 containing a cluster of five homologous NBS-LRR genes using this approach [28].

Data Analysis and Interpretation

Significance Thresholds and Interpretation

Establish appropriate significance thresholds is critical for credible association mapping:

Table 1: Statistical Parameters for GWAS Significance Interpretation

Parameter	Typical Threshold	Purpose	Example Application
Genome-wide Significance	P < 5×10⁻⁸ (for ~1M SNPs)	Control type I error rate	Human behavioral genetics [35]
Suggestive Significance	P < 1×10⁻⁵	Identify potential loci for validation	Plant stress tolerance studies [33]
Linkage Disequilibrium Decay	r² = 0.1-0.2	Define candidate intervals	~309 kb in soybean [30]
FDR Threshold	Q < 0.05	Control false discoveries	Multi-trait studies [31]
Phenotypic Variance Explained (R²)	5-30% per locus	Estimate effect size	Rice salt tolerance QTLs [34]

Candidate Gene Prioritization

Once significant associations are detected near NBS loci, prioritize candidates using:

• Gene Expression Data: Check if candidate genes are expressed in relevant tissues under stress conditions.
• Haplotype Analysis: Group significant SNPs into haplotypes and test for association with tolerance phenotypes [30] [34].
• Non-synonymous Mutations: Identify coding variants that alter protein function in tolerant vs. susceptible accessions.
• Co-localization with Known QTLs: Check if associations overlap with previously mapped tolerance QTLs [29] [28].

In the sorghum rust resistance study, researchers identified 61 point mutations in the first exon of a candidate NBS-LRR gene (Sobic.008G178200) that defined seven haplotypes with differing resistance levels [28].

Validation and Functional Characterization

Experimental Validation

GWAS hits require experimental validation to establish causality:

• Haplotype-Based Phenotypic Validation: Test for significant phenotypic differences among haplotypes of the candidate NBS gene [30].
• Gene Expression Analysis: Perform qRT-PCR to examine expression patterns under stress conditions in tolerant vs. susceptible accessions [34].
• Mutagenesis: Use CRISPR/Cas9 to create knockout mutants and validate gene function [34].

For example, in a rice salt tolerance study, CRISPR/Cas9 mutants of LOC_Os02g36880 showed altered salinity tolerance, confirming its role in salt stress response [34].

Pathway and Network Analysis

Place candidate NBS genes in biological context:

• Gene Ontology Enrichment: Identify over-represented biological processes among candidate genes [30] [31].
• Protein-Protein Interaction Networks: Map candidate genes onto known defense signaling networks.
• Cis-Element Analysis: Identify stress-responsive promoter elements in candidate NBS genes.

A soybean flooding tolerance study predicted eight candidate genes and three hub genes through protein-protein interaction network analysis, providing insights into molecular mechanisms [31].

Research Reagent Solutions

Table 2: Essential Research Reagents and Resources for GWAS of NBS-Mediated Tolerance

Category	Specific Product/Resource	Application	Key Considerations
Genotyping Platforms	7K SNP array (rice) [33], 180K Axiom Soya SNP array [30], RAD-seq [36]	High-density genotyping	Balance between density, cost, and sample throughput
Bioinformatics Tools	TASSEL (MLM) [34], PLINK [34], ADMIXTURE [32], LDBlockShow [34]	Population genetics & association analysis	Compatibility with species and sample size
NBS Annotation Resources	InterProScan, Pfam (NBS domain models), RGAugury	Identify NBS-LRR genes in genomes	Sensitivity for divergent NBS domains
Validation Reagents	CRISPR/Cas9 vectors [34], qPCR reagents [34], antibodies for protein detection	Functional characterization	Species-specific optimization required
Germplasm Resources	NPGS Sudan sorghum core [28], Korean landrace soybeans [30], IRRI rice accessions [36]	Source of natural variation	Photoperiod sensitivity, adaptation

Troubleshooting and Optimization

Common Challenges and Solutions

• Population Structure: If structure is confounding associations, increase the number of principal components in the model or use a relatedness matrix.
• Low Heritability: If trait heritability is low (<0.3), increase replication and improve phenotyping precision.
• Missing Heritability: If significant associations explain little phenotypic variance, consider rare variants, structural variations, or epistatic interactions.
• NBS Gene Family Complexity: If distinguishing between tandemly duplicated NBS genes, use long-read sequencing to resolve complex regions.

In peanut research, pangenome approaches have been developed to capture structural variations missing from single reference genomes, revealing SVs associated with seed size and weight traits [32].

This protocol provides a comprehensive framework for linking NBS loci to tolerance traits using GWAS. The integrated approach—combining diverse germplasm, high-quality phenotyping and genotyping, robust statistical models, and functional validation—enables researchers to dissect the genetic architecture of stress tolerance and identify causal NBS genes. These candidate genes can subsequently be deployed in marker-assisted breeding or gene editing programs to develop improved cultivars with enhanced tolerance to biotic and abiotic stresses.

This Application Note provides a detailed protocol for using RNA sequencing (RNA-seq) to identify and characterize responsive NBS-LRR genes within the context of research on genetic variation in stress-tolerant plant accessions. The NBS-LRR gene family is the largest class of plant disease resistance (R) proteins, serving as key determinants of the plant immune system by recognizing pathogen effectors and triggering defense responses [2] [37] [38]. This document outlines a comprehensive workflow—from experimental design to data analysis—tailored for researchers and scientists aiming to link genetic variation in NBS genes to stress tolerance phenotypes.

Background: NBS-LRR Genes in Plant Immunity

NBS-LRR proteins are intracellular receptors central to effector-triggered immunity (ETI). They typically contain a conserved nucleotide-binding site (NBS) domain and a C-terminal leucine-rich repeat (LRR) domain [37] [38]. The NBS domain is responsible for nucleotide binding and signal transduction, while the LRR domain is often involved in specific pathogen recognition [2] [38]. Based on their N-terminal domains, they are classified into major subfamilies:

• TNL: Contains a Toll/Interleukin-1 Receptor (TIR) domain.
• CNL: Contains a Coiled-Coil (CC) domain.
• RNL: Contains a Resistance to Powdery Mildew 8 (RPW8) domain [2] [37].

The number of NBS-LRR genes varies dramatically across plant species, from 73 in Akebia trifoliata to over 2,000 in wheat (Triticum aestivum), reflecting their diverse and specialized roles in adaptation [2] [38]. Their expression is often modulated in response to both biotic stresses (e.g., fungal, bacterial, or viral infections) and abiotic stresses (e.g., drought, salinity), making them prime candidates for research on stress-tolerant accessions [37] [38].

Experimental Design and Workflow

A robust experimental design is crucial for generating meaningful and reproducible RNA-seq data. The overarching goal is to identify differentially expressed NBS-LRR genes between stress-treated and control samples from tolerant and susceptible accessions.

The diagram below outlines the complete experimental and computational workflow for pinpointing stress-responsive NBS genes.

Key Considerations for Experimental Design

• Biological Replicates: A minimum of three biological replicates per condition is essential to account for natural variation and ensure statistical power in downstream differential expression analysis [39].
• Controlling Batch Effects: Technical variations introduced during different stages of the experiment (e.g., RNA isolation, library preparation, sequencing runs) can confound results. To mitigate this, process samples from different experimental groups randomly and in parallel whenever possible [39].
• Phenotypic Data: Correlate RNA-seq findings with measurable phenotypic traits (e.g., lesion size, biomass, photosynthetic efficiency) to ensure differential gene expression links to actual tolerance mechanisms [40] [41].

Detailed Protocols for Key Experiments

Protocol: RNA Extraction, Library Preparation, and Sequencing

This protocol is adapted from methods used in recent plant transcriptome studies [40].

1. RNA Extraction from Plant Tissues

• Material: Collect tissue (e.g., leaves, roots) from control and stress-treated plants. Immediately freeze in liquid nitrogen and store at -80°C.
• Grinding: Grind frozen tissue to a fine powder in liquid nitrogen using a mortar and pestle.
• Extraction: Use a CTAB (cetyltrimethylammonium bromide)-based method or a commercial plant RNA extraction kit. The CTAB method is effective for polysaccharide- and polyphenol-rich plant tissues [40].
  • Add pre-warmed extraction buffer to the powder and incubate at 65°C.
  • Extract with chloroform:isoamyl alcohol (24:1) and centrifuge.
  • Precipitate RNA from the aqueous phase using LiCl or isopropanol.
• DNase Treatment: Treat the purified RNA with DNase to remove genomic DNA contamination [40].
• Quality Control: Assess RNA concentration and purity using a spectrophotometer (e.g., NanoDrop). Evaluate RNA integrity via agarose gel electrophoresis or an instrument such as the Agilent TapeStation (RIN > 7.0 is recommended) [40] [39].

2. Library Preparation and Sequencing

• Enrichment: Use oligo(dT) beads to enrich for messenger RNA (mRNA) from total RNA.
• Library Construction: Prepare sequencing libraries using a stranded library preparation kit (e.g., Illumina TruSeq Stranded mRNA). This involves cDNA synthesis, end-repair, adapter ligation, and PCR amplification [40].
• Sequencing: Sequence the libraries on an Illumina platform (e.g., NextSeq 500) to generate a minimum of 20-30 million single-end or paired-end reads per sample [39].

Protocol: Bioinformatics Analysis of RNA-seq Data

The following pipeline is optimized for accuracy, as demonstrated in benchmarking studies [42].

1. Read Trimming and Quality Control

• Tool: Use fastp [42] or Trim Galore.
• Action: Remove adapter sequences and low-quality bases. Generate QC reports before and after trimming.
• Parameters: Set trimming parameters based on the quality profile of your data. For example, fastp can be run with parameters --cut_front --cut_tail --cut_window_size 4 --cut_mean_quality 20 to perform sliding window trimming [42].

2. Read Alignment

• Tool: Use a splice-aware aligner such as HISAT2 [2].
• Reference Genome: Align cleaned reads to the appropriate reference genome for your species (e.g., Solanium lycopersicum for tomato).
• Output: A BAM file containing aligned reads.

3. Read Quantification

• Tool: Use HTSeq [39] or featureCounts.
• Action: Count the number of reads mapped to each gene based on the provided genome annotation file (GTF/GFF).
• Output: A raw count table for all genes in each sample.

4. Differential Expression Analysis

• Tool: Use edgeR [39] or DESeq2, which are based on a negative binomial distribution model suitable for RNA-seq count data.
• Action: Compare read counts between stress and control groups to identify Differentially Expressed Genes (DEGs). A common significance threshold is an adjusted p-value (FDR) < 0.05 and an absolute log2 fold change > 1.

5. Genome-Wide Identification of NBS-LRR Genes

• Tool: Use HMMER (v3.1b2 or later) [2] [37].
• Hidden Markov Model (HMM): Search the proteome of your species of interest using the PFAM model PF00931 (NB-ARC domain) [2] [38].
• Domain Validation: Confirm the presence of associated domains (TIR, CC, LRR) in the candidate proteins using the NCBI Conserved Domain Database (CDD) [2] or Pfam.
• Classification: Classify identified genes into subfamilies (CNL, TNL, RNL, etc.) based on their domain architecture [37].

6. Data Integration

• Action: Intersect the list of DEGs from Step 4 with the list of NBS-LRR genes from Step 5 to pinpoint responsive NBS genes.
• Downstream Analysis: Perform gene ontology (GO) enrichment analysis on the responsive NBS genes to identify overrepresented biological processes. Construct phylogenetic trees to understand evolutionary relationships.

The Scientist's Toolkit: Essential Research Reagents

Table 1: Key research reagents and resources for studying NBS gene expression.

Item	Function/Description	Example Products/Tools
RNA Extraction Kit	Isolate high-quality, intact total RNA from plant tissues, often challenging due to secondary metabolites.	CTAB method [40], PicoPure RNA Isolation Kit [39]
Library Prep Kit	Prepare stranded cDNA libraries for Illumina sequencing from mRNA.	Illumina TruSeq Stranded mRNA Kit, NEBNext Ultra DNA Library Prep Kit [39]
HMMER Suite	Identify protein sequences containing conserved domains (e.g., NBS) using profile hidden Markov models.	HMMER v3.1b2 with PF00931 model [2] [37]
NCBI CDD	Database for confirming and visualizing conserved protein domains in candidate NBS-LRR genes.	CD-Search Tool [2]
Alignment Software	Map RNA-seq reads to a reference genome, accounting for intron splicing.	HISAT2 [2], TopHat2 [39]
Differential Expression Tool	Statistically identify genes with significant expression changes between conditions.	edgeR [39], DESeq2, Cufflinks/cuffdiff [2]
WIZ degrader 4	WIZ degrader 4, MF:C20H25N3O4, MW:371.4 g/mol	Chemical Reagent
Wilfornine A	Wilfornine A, MF:C45H51NO20, MW:925.9 g/mol	Chemical Reagent

Case Study: NBS Gene Identification inNicotianaspp.

A 2025 study on three Nicotiana species (N. tabacum, N. sylvestris, N. tomentosiformis) provides an excellent example of this workflow in action [2].

• Identification: Researchers used HMMER with the PF00931 model to identify 1,226 NBS genes across the three genomes. N. tabacum, an allotetraploid, contained 603 members, roughly the sum of its parental species [2].
• Classification: Genes were classified by domain composition, revealing a majority of NBS-only proteins (45.5%), followed by CC-NBS (23.3%) [2].
• Expression Analysis: RNA-seq data from disease resistance experiments (e.g., against black shank and bacterial wilt) were analyzed. Reads were aligned with HISAT2, quantified, and DEGs were identified using Cuffdiff, leading to the discovery of key NBS genes involved in disease resistance, including a multi-disease resistance gene [2].

Advanced Application: Integrating Structural Variation

For a more comprehensive analysis, consider integrating structural variations (SVs). A 2024 study in Brassica napus demonstrated that SVs (deletions, insertions, inversions) can have a profound regulatory impact on gene expression [41].

• SV-eQTL Mapping: By coupling population-scale transcriptome data with SV maps, researchers identified nearly 286,000 SV-expression quantitative trait loci (SV-eQTLs). These SVs were associated with the expression changes of 73,580 genes, revealing a pervasive effect of SVs on the transcriptome [41].
• Implication: In your research on tolerant accessions, significant expression differences in key NBS-LRR genes could be linked to underlying SVs (e.g., a promoter insertion) rather than single nucleotide polymorphisms (SNPs). Integrating SV analysis can provide a more complete mechanistic explanation for observed expression variation [41].

Pathway Diagram: NBS-LRR Gene Function in Plant Immunity

The following diagram illustrates the central role of NBS-LRR proteins in plant immune signaling, connecting pathogen recognition to defense activation.

The integrated workflow of RNA-seq expression profiling and genome-wide NBS-LRR gene identification is a powerful approach for deciphering the genetic basis of stress tolerance in plants. By following the detailed protocols and leveraging the tools outlined in this document, researchers can systematically pinpoint key responsive NBS genes, thereby identifying prime candidates for future functional validation and crop improvement strategies.

Haplotype Analysis of NBS Genes in Tolerant versus Susceptible Accessions

Nucleotide-binding site-leucine rich repeat (NBS-LRR) genes constitute the largest family of plant disease resistance (R) genes and play a critical role in effector-triggered immunity [43]. These genes are characterized by a conserved NBS domain and variable LRR regions that determine pathogen recognition specificity [24]. Haplotype analysis of NBS genes provides powerful insights into the genetic basis of disease resistance, enabling the identification of superior alleles associated with tolerant accessions [5]. This approach has become fundamental for understanding plant-pathogen co-evolution and informing molecular breeding strategies [5] [43].

The NBS gene family exhibits extraordinary diversity across plant species, with significant structural variation observed between tolerant and susceptible genotypes [24] [5]. Studies in chickpea identified 121 NBS-LRR genes, while research across 34 plant species revealed 12,820 NBS-domain-containing genes with 168 distinct domain architecture patterns [24] [43]. This natural variation provides the genetic substrate for haplotype-based studies aimed at elucidating the molecular mechanisms of disease resistance.

Key Concepts and Genetic Principles

NBS Gene Classification and Domain Architecture

NBS-encoding genes are classified based on their domain composition at the N-terminus. The two major subclasses are TIR-NBS-LRR (TNL) and CC-NBS-LRR (CNL), distinguished by the presence of Toll/interleukin-1 receptor (TIR) or coiled-coil (CC) domains, respectively [24] [43]. A comprehensive comparative analysis across land plants has identified both classical and species-specific structural patterns, including TIR-NBS-TIR-Cupin1-Cupin1 and Sugar_tr-NBS architectures [24].

Table 1: NBS Gene Classification and Distribution in Various Plant Species

Plant Species	Total NBS Genes	TNL Genes	CNL Genes	Truncated/Other	Genome Distribution
Chickpea (Cicer arietinum)	121	35	63	23	Uneven across 8 chromosomes, ~50% in clusters [43]
Asian Pear (P. bretschneideri)	338	36 (10.95%)	90 (26.6%)	212 (62.45%)	-
European Pear (P. communis)	412	-	38 (9.22%)	374 (90.78%)	-
Various Land Plants	12,820	18,707 (across 304 angiosperms)	70,737 (across 304 angiosperms)	1,847 RNL genes	603 orthogroups identified [24]

Population Genetic Variation in NBS Genes

Comparative studies between Asian and European pears revealed that 15.79% of orthologous NBS gene pairs exhibited Ka/Ks ratios >1, indicating strong positive selection following species divergence [5]. Domestication has differentially affected nucleotide diversity, with Asian pear cultivars showing decreased diversity (Ï€=6.23E-03 in cultivated vs. 6.47E-03 in wild), while European pears displayed the opposite trend (Ï€=6.48E-03 in cultivated vs. 5.91E-03 in wild) [5].

Genetic variation analysis between tolerant and susceptible cotton accessions identified substantially more unique variants in the tolerant genotype (Mac7: 6583 variants) compared to the susceptible variety (Coker312: 5173 variants) [24]. This pattern highlights the role of genetic diversity in disease resistance mechanisms.

Experimental Protocols

Genome-Wide Identification of NBS-Encoding Genes

Protocol: Identification and Classification of NBS Genes

• Step 1: Sequence Retrieval - Obtain latest genome assemblies from publicly available databases (NCBI, Phytozome, Plaza) [24].
• Step 2: Domain Screening - Use PfamScan.pl HMM search script with default e-value (1.1e-50) and Pfam-A_hmm model to identify genes containing NB-ARC domains [24].
• Step 3: Architecture Classification - Classify genes based on domain architecture using established systems grouping similar patterns into classes [24].
• Step 4: Phylogenetic Analysis - Perform multiple sequence alignment using MAFFT 7.0 and construct phylogenetic trees with FastTreeMP (1000 bootstrap replicates) [24].

Haplotype Analysis of NBS Genes in Tolerant and Susceptible Accessions

Protocol: Comparative Haplotype Analysis

• Step 1: Population Selection - Assemble diverse germplasm including tolerant and susceptible accessions. For example, in wheat salt tolerance studies, a panel of 228 hexaploid spring wheat accessions was used [44].
• Step 2: Genotyping - Conduct whole-genome resequencing or high-density SNP array analysis. Map sequencing reads using BWA software with default parameters [45].
• Step 3: Variant Calling - Identify genomic variations using GATK following best practices. Annotate variant sites using SnpEff with parameter "-upDownStreamLen 2000" [45].
• Step 4: Haplotype Block Definition - Use linkage disequilibrium (LD)-based grouping to identify haplotype blocks within NBS gene regions [44] [46].
• Step 5: Phenotype-Genotype Association - Perform genome-wide association studies (GWAS) to identify significant associations between haplotypes and disease resistance traits [5] [43].

Diagram 1: Experimental workflow for haplotype analysis of NBS genes in tolerant versus susceptible accessions.

Expression Analysis of Candidate NBS Genes

Protocol: Expression Profiling of NBS Genes

• Step 1: Experimental Design - Collect tissue samples from resistant and susceptible genotypes at multiple time points after pathogen inoculation [43].
• Step 2: RNA Extraction and Sequencing - Extract total RNA using standard protocols (e.g., CTAB method). Prepare libraries and sequence using appropriate platforms [45].
• Step 3: Transcript Quantification - Map RNA-seq reads to reference genomes using STAR with parameters "--outFilterMultimapNmax 1 --twopassMode Basic". Calculate expression values using RSEM software [45].
• Step 4: Differential Expression Analysis - Identify significantly differentially expressed genes using criteria such as |log2FC| > 1 and p < 0.05 [47].
• Step 5: Validation - Confirm expression patterns using real-time quantitative PCR on independent biological replicates [43].

Data Analysis and Interpretation

Case Study: NBS Genes in Chickpea-Ascochyta Blight Interaction

A comprehensive study in chickpea identified 121 NBS-LRR genes, 30 of which co-localized with previously reported ascochyta blight QTLs [43]. Expression profiling in resistant (CDC Corinne, CDC Luna) and susceptible (ICCV 96029) genotypes revealed that 27 NBS-LRR genes showed differential expression in response to Ascochyta rabiei infection [43]. Five NBS-LRR genes demonstrated genotype-specific expression patterns, highlighting their potential role in differential disease resistance [43].

Table 2: Haplotype Effects on Agronomic Traits in Stress Conditions

Crop Species	Stress Condition	Gene / Haplotype Block	Trait Association	Effect Size
Wheat [46]	High night-time temperature (HNT)	Hap1 TraesCS1A02G305700	Higher biomass and spike number under HNT	Significant (p<0.05)
Wheat [46]	High night-time temperature (HNT)	Hap1 TraesCS2B02G599800	Higher biomass under HNT	Significant (p<0.05)
Wheat [46]	High night-time temperature (HNT)	Hap2 TraesCS4B02G264300	Higher biomass under both control and HNT	Significant (p<0.05)
Rice [48]	Low nitrogen	LOC_Os06g06440 (AA haplotype)	Higher relative grain yield (RGY)	0.95 and 1.53 in MAGIC and germplasm varieties
Pear [5]	Domestication	Pbr025269.1 and Pbr019876.1	>5x upregulation in wild vs cultivated; >2x after A. alternata inoculation	Expression differences significant

Identification of Superior Haplotypes

Superior haplotypes can be identified through association analysis between haplotype variations and phenotypic performance. In wheat salt tolerance research, haplotype analysis identified superior haplotypes of genes encoding a sodium symporter (TraesCS1B02G413800) and peptide transporter (TraesCS5A02G004400) [44]. These superior haplotypes were predominantly present in landraces but often lost in modern cultivars due to artificial selection during breeding programs [44].

The Scientist's Toolkit

Table 3: Essential Research Reagents and Computational Tools for NBS Gene Haplotype Analysis

Category	Tool/Reagent	Specific Function	Application Notes
Sequencing & Genotyping	Illumina HiSeq Platform	Whole genome sequencing	Provides high-quality sequencing data for variant discovery [45]
	PacBio de novo sequencing	Genome assembly	Enables telomere-to-telomere genome assembly [49]
Bioinformatics Tools	BWA (0.7.17-r1188)	Read alignment to reference genome	"BWA-MEM" algorithm with default parameters recommended [45]
	GATK (Genome Analysis Toolkit)	Variant calling	Follow GATK best practices for optimal results [45]
	SnpEff (4.3s)	Variant annotation	Use parameter "-upDownStreamLen 2000" for comprehensive annotation [45]
	OrthoFinder v2.5.1	Orthogroup analysis	Identifies evolutionary relationships among NBS genes [24]
	STAR (2.7.1a)	RNA-seq read alignment	Parameters: --outFilterMultimapNmax 1 --twopassMode Basic [45]
	RSEM	Expression quantification	Calculates FPKM values for gene expression analysis [45]
Experimental Materials	CTAB method	DNA extraction from leaves	Standard protocol for high-quality DNA [45]
	Field-based heat tents	Phenotyping under stress conditions	Custom-designed for HNT stress studies in wheat [46]

Functional Validation Strategies

Gene Expression Manipulation

Virus-induced gene silencing (VIGS) provides an efficient approach for functional validation. Silencing of GaNBS (OG2) in resistant cotton demonstrated its putative role in virus tittering, confirming its importance in disease resistance [24]. This method allows for rapid assessment of gene function without the need for stable transformation.

Protein Interaction Studies

Protein-ligand and protein-protein interaction analyses can reveal molecular mechanisms. Studies in cotton showed strong interaction of putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus, providing insights into resistance mechanisms [24].

Diagram 2: NBS-LRR protein mediated defense signaling pathway in plants.

Haplotype analysis of NBS genes in tolerant versus susceptible accessions provides powerful insights into the genetic basis of disease resistance in plants. The protocols outlined in this application note enable comprehensive identification, characterization, and validation of NBS gene haplotypes associated with disease resistance. The integration of genomic, transcriptomic, and functional validation approaches facilitates the discovery of superior haplotypes that can be deployed in marker-assisted breeding programs to enhance crop disease resistance.

The Ka/Ks ratio, also known as ω or dN/dS ratio, is a fundamental metric in molecular evolution that estimates the balance between neutral mutations, purifying selection, and beneficial mutations acting on protein-coding genes [50]. This ratio compares the rate of non-synonymous substitutions (Ka), which alter the amino acid sequence, against the rate of synonymous substitutions (Ks), which do not change the encoded protein [50].

Synonymous substitutions are typically considered neutral evolutionary events, as they do not affect protein function. Therefore, the Ka/Ks ratio serves as a powerful indicator of selective pressure:

• Ka/Ks > 1: Suggests positive selection, where non-synonymous mutations are favored, often driving adaptive evolution.
• Ka/Ks = 1: Indicates neutral evolution, where sequence changes are neither favored nor selected against.
• Ka/Ks < 1: Reflects purifying selection, which removes deleterious mutations to conserve protein function and sequence [50].

This analytical framework is particularly valuable for investigating plant disease resistance genes, such as the Nucleotide-Binding Site Leucine-Rich Repeat (NBS-LRR) family. These genes frequently exhibit signatures of positive selection in pathogen-interaction domains, revealing an evolutionary arms race between plants and their pathogens [51].

Article 2: Key Interpretations and Computational Considerations of Ka/Ks

Table 1: Interpretation of Ka/Ks Ratio Values

Ka/Ks Value	Type of Selection	Evolutionary Interpretation	Example in NBS Genes
> 1	Positive/Darwinian Selection	Diversifying changes are advantageous; driven by adaptive evolution.	LRR domains in pathogen recognition [51].
= 1	Neutral Evolution	No selective constraint; mutations are fixed randomly.	Rare in functional R genes; may indicate pseudogenization.
< 1	Purifying/Stabilizing Selection	Conservation of amino acid sequence; deleterious changes are removed.	NBS (NB-ARC) domain for signal transduction [24].

Computational Methods and Limitations

Several methods exist for calculating Ka/Ks ratios, each with specific strengths. Approximate methods (e.g., Nei & Gojobori) involve counting sites and substitutions but may oversimplify. Maximum-likelihood methods (e.g., implemented in codeml/PAML) use probability theory to simultaneously estimate parameters like divergence and transition/transversion bias, offering greater robustness [50].

Key limitations must be considered:

• Selective Heterogeneity: A gene average Ka/Ks = 1 can mask a mixture of strongly purifying and positively selected sites [50].
• Beyond Coding Regions: Ka/Ks detects selection only in protein-coding regions, missing evolutionary changes in regulatory sequences [50].
• Synonymous Mutation Neutrality: Recent evidence suggests some synonymous mutations are not neutral, potentially challenging the core assumption [50].

For NBS-LRR genes, site-specific models in maximum-likelihood frameworks are crucial as they can identify individual codons under positive selection, often found in the solvent-exposed residues of the LRR domain [51].

Article 3: Experimental Protocol for Ka/Ks Analysis on NBS Gene Families

This protocol provides a step-by-step workflow for conducting a Ka/Ks analysis to detect selection pressure on NBS-LRR genes using data from resistant and susceptible plant accessions.

Workflow Diagram

Step-by-Step Protocol

Step 1: Data Collection and Identification of NBS Family Members

• Retrieve Genome Sequences: Obtain annotated genome assemblies for target species (e.g., Gossypium hirsutum accessions with differential disease tolerance) from databases such as Phytozome, CottonMD, or NCBI [52] [53] [2].
• Identify NBS-LRR Genes: Perform HMMER searches against the proteome using the NB-ARC (PF00931) hidden Markov model from Pfam [2] [24].
• Validate Domain Architecture: Confirm the presence of associated domains (TIR, CC, LRR) using the NCBI Conserved Domain Database (CDD) and Pfam [2]. Classify genes into subfamilies (TNL, CNL, NL, etc.).

Step 2: Sequence Alignment and Phylogenetic Analysis

• Multiple Sequence Alignment: Use tools like MUSCLE or MAFFT to align the protein sequences of the identified NBS genes [2] [24].
• Construct Phylogenetic Tree: Build a maximum-likelihood phylogenetic tree using MEGA11 or FastTreeMP with 1000 bootstrap replicates to assess evolutionary relationships and group sequences into orthologous clusters [52] [24].

Step 3: Ortholog Identification and Synteny Analysis

• Identify Orthologous Pairs: For comparative analysis between accessions or species, use BLASTP followed by synteny analysis with MCScanX to identify true orthologs [53] [2].
• Extract Coding Sequences: For each orthologous gene pair, extract the corresponding CDS sequences.

Step 4: Ka/Ks Calculation

• Prepare Sequence Pairs: Align the CDS sequences of each ortholog pair using codon-aware alignment in PAL2NAL.
• Calculate Substitution Rates: Use KaKs_Calculator 2.0 with the Nei-Gojobori (NG) or other models to compute Ka, Ks, and the Ka/Ks ratio for each gene pair [2].

Step 5: Statistical Testing and Identification of Positive Selection

• Site-Specific Models: To detect positive selection acting on specific codons, use the codeml program in the PAML package. Compare a null model that does not allow sites with ω > 1 (e.g., M7) against an alternative model that does (e.g., M8) using a likelihood ratio test [51].
• Bayesian Empirical Bayes: Identify specific codons with a high posterior probability of belonging to the class under positive selection.

Article 4: The Scientist's Toolkit

Table 2: Essential Research Reagents and Tools for Ka/Ks Analysis

Tool/Reagent	Function/Description	Application in NBS Gene Studies
HMMER Suite	Identifies protein domains using hidden Markov models.	Genome-wide identification of NBS-LRR genes using PF00931 (NB-ARC domain) [52] [24].
NCBI CDD	Database of conserved domain alignments.	Validation of NBS, TIR, CC, and LRR domain architecture in candidate genes [52] [2].
MCScanX	Analyzes genomic collinearity and gene duplication events.	Differentiates between orthologs and paralogs; identifies segmental/tandem duplications in NBS clusters [52] [2].
KaKs_Calculator	Computes Ka and Ks values from codon-aligned CDS.	Calculating substitution rates for orthologous NBS gene pairs [2].
PAML (codeml)	Uses maximum likelihood for phylogenetic analysis.	Detecting site-specific positive selection in NBS-LRR domains [51].
MEME Suite	Discovers conserved motifs in nucleotide or protein sequences.	Analyzing motif conservation and diversity among NBS-LRR subfamilies [52].

Article 5: Analysis and Interpretation of Results in the Context of NBS Genes

Interpreting Ka/Ks Output in Disease Resistance Studies

In NBS-LRR gene research, different protein domains exhibit distinct evolutionary pressures. The LRR domain often shows the strongest signal of positive selection (Ka/Ks > 1) because its solvent-exposed residues directly interact with rapidly evolving pathogen effectors [51] [54]. In contrast, the NBS (NB-ARC) domain, crucial for nucleotide binding and activation signaling, is typically under strong purifying selection (Ka/Ks < 1) to maintain functional integrity [24].

Case Study: NBS Gene Evolution in Tung Trees

A 2024 study on tung trees (Vernicia fordii and V. montana) identified 239 NBS-LRR genes across the two genomes [53] [25]. Researchers identified orthologous gene pairs and analyzed their evolutionary rates. One key finding was the orthologous pair Vf11G0978-Vm019719, which showed differential expression and potentially different evolutionary constraints linked to Fusarium wilt resistance [53] [25].

Structural Mapping of Positive Selection

A landmark study in Arabidopsis thaliana mapped positively selected sites in NBS-LRR genes onto protein secondary structure [51]. The analysis found that positively selected positions were disproportionately located in the LRR domain, particularly in a nine–amino acid β-strand submotif that is likely solvent-exposed [51]. This provides mechanistic insight into how plant R genes adapt to recognize changing pathogen targets.

Article 6: Advanced Applications and Integration with Functional Studies

Correlating Evolutionary Signatures with Expression Data

Modern evolutionary genomics integrates Ka/Ks analysis with functional data. For example:

• RNA-seq Analysis: Following Ka/Ks calculations, researchers can perform expression profiling of NBS genes under stress conditions (e.g., pathogen infection, drought) in tolerant and susceptible accessions [2] [24].
• Differential Expression: Identify NBS genes under positive selection that are also upregulated during defense responses. These represent high-priority candidates for functional validation.

Functional Validation through Gene Silencing

The role of candidate NBS genes identified via evolutionary analysis can be tested experimentally. Virus-Induced Gene Silencing (VIGS) is a powerful technique for this purpose [53] [24]. For instance, silencing of a specific NBS gene in a resistant tung tree accession (Vm019719) demonstrated its critical role in conferring resistance to Fusarium wilt [53] [25].

Evolutionary Dynamics and Duplication Events

Understanding the expansion mechanisms of the NBS family is crucial. Whole-genome duplication (WGD) and tandem duplication are primary drivers. Ka/Ks analysis of duplicated genes can reveal their evolutionary fates. Most duplicates are under purifying selection, but some can undergo neofunctionalization, acquiring new resistance specificities [52] [2].

Navigating Technical Challenges in NBS Gene Sequencing and Analysis

Addressing Mapping Issues in High-Homology NBS Regions with Short-Read Sequencing

Regions of high sequence homology present one of the most significant technical challenges for short-read next-generation sequencing (NGS) in genetic variation analysis. This is particularly problematic when studying Nucleotide-Binding Site (NBS) genes, which constitute the largest class of disease resistance (R) genes in plants and play critical roles in pathogen defense mechanisms [5]. The inherent characteristics of NBS genes—including tandem duplications, conserved protein domains, and extensive paralogous families—create genomic contexts where short sequencing reads cannot be uniquely mapped to their correct genomic origin [55] [56]. This mapping ambiguity leads to both false-positive and false-negative variant calls, potentially compromising the identification of genetic variations associated with disease tolerance in plant accessions. For researchers investigating the genetic basis of stress tolerance, these technical limitations can obstruct the discovery of meaningful associations between genotype and phenotype. This application note details the specific challenges and provides validated protocols to overcome mapping issues in high-homology NBS regions, enabling more reliable genetic variation analysis in tolerance research.

Quantitative Assessment of Mapping Issues

Impact of Read Length and Homology on Mapping Accuracy

Experimental simulations on genomic data have quantified the relationship between read length, degree of homology, and mapping accuracy. One study systematically evaluating 158 genes relevant to newborn screening (including many with high homology) demonstrated that while longer reads significantly improve mapping performance, certain highly homologous regions remain problematic even at 250 bp read lengths [55].

Table 1: Effect of Read Length on Mapping Accuracy in High-Homology Regions

Read Length	Correctly Mapped Reads	Incorrectly Mapped Reads	Unmapped Reads	Average Depth of Coverage
75 bp	>99%	<0.5%	<0.5%	78.2X ± 48.3
100 bp	>99%	<0.3%	<0.3%	85.7X ± 45.1
150 bp	>99%	<0.2%	<0.2%	92.3X ± 41.7
250 bp	>99%	<0.1%	<0.1%	95.8X ± 39.2

Despite high overall mapping rates, the critical finding was that four genes (SMN1, SMN2, CBS, and CORO1A) exhibited low-coverage regions within exons even at the longest read length (250 bp) [55]. These problem genes shared a common characteristic: nearly zero mismatches and minimal differences in alignment length compared to their homologous counterparts. This indicates that the degree and length of homology are the primary factors affecting mapping success rather than general read length considerations.

NBS Gene Classification and Variation Patterns

In plant genomes, NBS-encoding genes display distinct evolutionary patterns that contribute to mapping challenges. A comparative analysis of Asian (P. bretschneideri) and European (P. communis) pear genomes identified 338 and 412 NBS-encoding genes respectively, with different distributions across structural classes [5].

Table 2: NBS-Encoding Gene Classification in Asian and European Pear Genomes

Gene Class	P. bretschneideri (Asian)	P. communis (European)	Primary Difference Driver
NBS-LRR	36.4% (123 genes)	25.7% (106 genes)	Proximal duplication
CC-NBS-LRR	26.6% (90 genes)	9.22% (38 genes)	Proximal duplication
TIR-NBS-LRR	10.95% (37 genes)	19.9% (82 genes)	Species-specific expansion
NBS (truncated)	10.36% (35 genes)	24.0% (99 genes)	Differential evolution
CC-NBS	9.46% (32 genes)	7.04% (29 genes)	Relative conservation
TIR-NBS	6.21% (21 genes)	13.35% (55 genes)	Species-specific expansion

The study further revealed that approximately 15.79% of orthologous NBS gene pairs exhibited Ka/Ks ratios greater than one, indicating strong positive selection following the divergence of Asian and European pears [5]. This rapid evolution under positive selection creates additional challenges for read mapping, as reference genomes may not capture the full diversity of NBS genes across different accessions.

Bioinformatic Solutions for High-Homology Regions

Multi-Region Joint Detection (MRJD) Algorithm

To address the fundamental limitations of conventional mapping approaches, Illumina developed the Multi-Region Joint Detection (MRJD) algorithm as part of their DRAGEN v4.3 platform [57]. Unlike traditional variant callers that rely on aligners to uniquely map reads before variant detection, MRJD employs a fundamentally different approach:

• Paralogous Region Identification: MRJD first identifies all genomic locations with high sequence similarity to the target region.

• Ambiguous Read Retention: Instead of discarding reads with multiple possible mappings, MRJD retains all potentially informative reads regardless of mapping quality.

• Joint Haplotype Construction: The algorithm builds haplotypes simultaneously across all paralogous regions using both uniquely and ambiguously mapped reads.

• Variant Placement: Probabilistic models determine the most likely origin of variants across the homologous regions.

When applied to the challenging PMS2/PMS2CL region (associated with Lynch syndrome), MRJD demonstrated significant improvements over conventional methods, achieving 99.7% recall for SNVs and 97.1% recall for indels in high-homology regions [57]. The algorithm operates in two modes: a balanced default mode and a high-sensitivity mode recommended for applications where detection of all potential pathogenic variants is paramount.

Practical Implementation of MRJD

The MRJD workflow can be implemented through the following protocol:

Protocol 1: MRJD Implementation for NBS Gene Analysis

• Platform Requirements: Illumina DRAGEN v4.3 or later with MRJD capability enabled.

• Input Data Preparation:

  • Whole-genome sequencing data (PCR-free preferred)
  • Minimum coverage: 30X across target regions
  • Read length: 150 bp recommended

• Parameter Configuration:

  • Enable "MRJD High Sensitivity" mode for maximal recall
  • Specify target NBS genes and known paralogous regions
  • Set minimum mapping quality threshold to 0 to retain ambiguous reads

• Output Interpretation:

  • Review variants flagged with ambiguous placement
  • Cross-reference with orthogonal validation for critical findings
  • Utilize population frequency data to filter likely spurious calls

In performance validation using 22 non-cell-line samples, MRJD successfully detected all expected clinically relevant variants in high-homology regions, demonstrating its utility in real-world settings [57].

Experimental Design Considerations

Sequencing Strategy Optimization

Based on empirical studies, the following sequencing strategies significantly improve mapping in high-homology NBS regions:

Longer Read Lengths: While not sufficient alone, increasing read length from 75 bp to 250 bp reduces incorrectly mapped reads by 80% and improves coverage consistency in moderately homologous regions [55].

Paired-End Sequencing: Minimum 2×100 bp paired-end reads provide critical mapping information from both ends of DNA fragments, helping to anchor reads in correct genomic locations.

Increased Coverage Depth: Target 50-100X coverage for NBS regions of interest to ensure sufficient unique reads for variant calling after filtering ambiguous mappings.

Reference Genome Selection and Preparation

Pan-Genome References: For species with available pan-genome resources (e.g., maize with 26 high-quality genomes), mapping to multiple references can rescue variants missing from single reference genomes [58]. In maize annexin genes, 9 of 12 genes were "core" (present in all 26 lines), while 3 were "near-core" (present in 24-25 lines), highlighting the limitation of single-reference approaches [58].

Custom Reference Modification: For critical projects, consider creating a comprehensive reference by incorporating known NBS haplotypes from tolerant accessions into a modified reference genome.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Research Reagents and Platforms for NBS Region Analysis

Reagent/Platform	Function	Application Note
Illumina DRAGEN v4.3+ with MRJD	Bioinformatic pipeline for variant calling in high-homology regions	Enables joint analysis of paralogous regions; supports PMS2, SMN1, SMN2, NEB, STRC, IKBKG, and TTN genes [57]
Twist Bioscience Target Enrichment	Custom capture panel for NBS gene regions	Used in BabyDetect project; covers 1.5-1.6 Mb target regions with high specificity [59]
QIAsymphony SP with DNA Investigator Kit	Automated DNA extraction from dried blood spots or plant tissues	Ensures high molecular weight DNA suitable for long-range PCR validation [59]
PCR-Free WGS Library Prep	Minimizes amplification bias in homologous regions	Reduces false structural variant calls in GC-rich NBS regions [57]
Humanomics v3.15 Pipeline	Open-source variant calling pipeline	Incorporates BWA-MEM, elPrep, HaplotypeCaller; customizable for non-model organisms [59]

Integrated Workflow for NBS Gene Variation Analysis

Comprehensive Validation Protocol

Protocol 2: Orthogonal Validation of NBS Region Variants

For variants identified in high-homology NBS regions, especially those with potential functional significance in tolerant accessions, orthogonal validation is essential:

• Long-Range PCR Design:

  • Design primers in unique flanking sequences ≥500 bp from homology boundaries
  • Amplify entire homologous region in a single amplicon when possible
  • Use high-fidelity polymerase with proofreading capability

• Sanger Sequencing:

  • Sequence long-range PCR products with internal primers
  • Ensure complete coverage of the variant position
  • Compare with NGS-derived variant calls

• Segmentation Analysis:

  • For structural variants, employ digital droplet PCR or MLPA
  • Establish copy number baseline using control genes

This comprehensive approach addresses the critical challenge of accurately identifying genetic variations in high-homology NBS regions, enabling more reliable association studies in tolerance research. By implementing these specialized protocols and bioinformatic solutions, researchers can overcome the technical limitations of short-read sequencing and advance our understanding of the genetic basis of disease resistance in plants and beyond.

Optimizing Variant Calling in Paralog-Rich Gene Families

The accurate identification of genetic variation within paralog-rich gene families presents a significant challenge in genomic analysis. These regions, characterized by sequences with high similarity due to gene duplication events, are prevalent in plant genomes, particularly within Nucleotide-Binding Site-Leucine Rich Repeat (NBS-LRR) gene families that mediate disease resistance [24] [5]. Standard variant calling pipelines frequently misalign short sequencing reads across paralogous sequences, leading to inaccurate variant representation and potentially missing clinically or functionally important mutations [60]. This application note details optimized experimental and computational protocols for reliable variant detection in these complicated genomic regions, with specific application to studying genetic variation in NBS genes of tolerant plant accessions.

The Challenge of Paralogs in Variant Calling

Paralogs, genes arising from duplication events, can exhibit high sequence similarity, complicating their individual analysis. In short-read sequencing data, reads originating from one paralogous region often align equally well to other homologous regions, resulting in low mapping quality scores and ambiguous alignments [60]. Consequently, variant callers may fail to identify true variants or may generate false positives.

The NBS-LRR gene family, one of the largest classes of plant disease resistance (R) genes, is particularly affected by these challenges due to its tendency to expand through tandem and segmental duplications [24] [5]. Studies in pear genomes have revealed that proximal duplications primarily drive differences in NBS-encoding gene numbers between species, creating regions where standard variant calling approaches are insufficient [5]. Furthermore, ectopic gene conversion events, where a sequence is copied from one genomic region to a distant paralogous region, can introduce new genetic variation that remains undetectable with conventional methods [60].

Optimized Experimental Design

Sequencing Technology Selection

Choosing appropriate sequencing technology is fundamental to addressing paralog-related challenges. While short-read sequencing (SRS) platforms like Illumina NovaSeq provide high accuracy for unique genomic regions, long-read sequencing (LRS) technologies from PacBio or Oxford Nanopore generate reads that can span repetitive regions and paralog-specific variants.

Table 1: Sequencing Platform Comparison for Paralog-Rich Regions

Technology	Read Length	Advantages for Paralogs	Limitations
Illumina SRS	75-300 bp	High per-base accuracy, low cost	Limited ability to resolve paralogs
PacBio HiFi	15-20 kb	High accuracy (Q30+), long reads	Higher DNA input requirements
Oxford Nanopore	>10 kb	Very long reads, direct epigenetics	Higher error rate requires correction

For critical applications involving paralogous regions, a hybrid approach using both SRS and LRS provides optimal results. LRS enables phasing of haplotypes and resolution of structural variants, while SRS offers cost-effective high coverage for validation [61]. Recent advances have made LRS accessible for low-input samples, including dried blood spots and low-yield DNA extracts, expanding its application to diverse sample types [61].

Targeted Enrichment Strategies

For focused studies on specific paralogous gene families like NBS genes, targeted enrichment approaches can be employed:

• Custom capture panels designed to target unique regions (SUNs - singly unique nucleotides) that differentiate paralogs
• Long-range PCR amplification of specific paralogous regions followed by sequencing
• Hybridization-based enrichment using probes designed against divergent regions of paralogous genes

When designing custom capture panels for NBS genes, it is crucial to identify regions with sufficient sequence divergence to ensure specific capture. Tools such as Paralog Explorer can help identify such regions by providing information on paralogous gene pairs and their sequence similarities [62] [63].

Computational Methods for Variant Detection in Paralogs

Specialized Variant Calling Pipelines

Standard variant callers like GATK HaplotypeCaller and DeepVariant demonstrate excellent performance in unique genomic regions but show zero sensitivity in perfectly identical paralogous regions [60]. Specialized approaches are required to overcome this limitation.

The Chameleolyser method specifically addresses variant calling in paralogous regions by:

• Extracting reads from homologous regions (approximately 3.5% of the exome)
• Re-aligning them to a reference genome where all but one paralog per set are masked
• Performing sensitive variant calling on the uniquely aligned reads [60]

This approach enables identification of single nucleotide variants (SNVs), small insertions/deletions (Indels), copy number variants (CNVs), and ectopic gene conversion events in paralogous regions that remain undetectable by standard analysis. Application to 41,755 exome samples revealed 2,529,791 rare SNVs/Indels and 338,084 variants resulting from gene conversion events, none detectable by regular analysis techniques [60].

Read Alignment and Processing

Optimized alignment strategies are critical for paralog-rich regions:

• Reference genome masking: Mask all but one paralog in each set to force unique alignment of reads [60] [24]
• Alternative aligners: Use aligners with sensitive settings for repetitive regions
• Duplicate marking: Identify and mark PCR duplicates using tools like Picard or Sambamba to prevent artificial inflation of coverage [64] [65]

The following workflow diagram illustrates the optimized computational pipeline for variant calling in paralog-rich regions:

Validation and Quality Control

Rigorous validation is essential for variants called in paralogous regions:

• Orthogonal validation using long-read sequencing technologies [61] [60]
• Trio analysis to verify inheritance patterns (99.0% of SNVs/Indels in offspring should be present in parents) [60]
• Comparison with benchmark resources such as Genome in a Bottle (GIAB) for human studies [64] [65]

Table 2: Validation Metrics for Paralogous Variant Calls

Validation Method	Expected Performance	Application
PacBio HiFi LRS	>88% concordance for SNVs/Indels	Gold-standard validation
Inheritance Checking	99.0% transmission in trios	Quality control in family studies
GIAB Benchmarking	>83% concordance in high-confidence regions	Pipeline performance assessment

Application to NBS Gene Analysis in Tolerant Accessions

Identifying NBS Genes and Paralogs

The first step in analyzing variation in NBS genes is comprehensive identification of all NBS-encoding genes in the target genome. This involves:

• HMM-based identification using PfamScan with the NB-ARC domain (PF00931) as query [24]
• Classification based on domain architecture (TIR-NBS-LRR, CC-NBS-LRR, NBS-LRR, etc.) [24] [5]
• Paralog group definition using tools like OrthoFinder and DIOPT to identify paralogous relationships [24] [62]

In pear genomes, this approach identified 338 and 412 NBS-encoding genes in Asian and European pears, respectively, with differential expansion primarily driven by proximal duplications [5]. Similarly, analysis across 34 plant species identified 12,820 NBS-domain-containing genes with 168 distinct domain architecture classes [24].

Genetic Variation Analysis

Once NBS genes and their paralogous relationships are defined, optimized variant calling can be applied to identify genetic variation between tolerant and susceptible accessions:

• Population sequencing of wild and domesticated accessions to identify selection signatures [5]
• Comparative analysis of NBS gene content and variation between resistant and susceptible varieties [45]
• Expression correlation of genetic variants with gene expression changes under pathogen stress [5]

In Asian and European pears, domestication caused contrasting patterns of genetic diversity in NBS genes, with decreased nucleotide diversity in Asian cultivars but increased diversity in European cultivars, suggesting independent domestication histories [5]. Many NBS-encoding genes showed Ka/Ks ratios >1, indicating positive selection has shaped their diversity.

Functional Validation

Candidate variants identified through optimized calling require functional validation:

• Virus-Induced Gene Silencing (VIGS) to test gene function, as demonstrated with GaNBS in cotton, which showed a putative role in virus tittering [24]
• Expression analysis after pathogen inoculation to identify responsive NBS genes [5]
• Protein-ligand interaction studies to validate the functional impact of variants on protein function [24]

Research Reagent Solutions

Table 3: Essential Research Reagents and Tools for Paralog-Rich Variant Analysis

Category	Specific Tools/Reagents	Application	Key Features
Variant Callers	Chameleolyser, GATK, DeepVariant	Specialized detection in paralogs	Masking-based approach for paralogs
Paralog Identification	DIOPT, Paralog Explorer, OrthoFinder	Defining paralogous relationships	Integrates multiple prediction algorithms
Sequencing Kits	PacBio SMRTbell prep kit 3.0, Illumina cfDNA/FFPE kit	Library preparation	Optimized for low-input or long-read data
Alignment Tools	BWA-MEM, minimap2, STAR	Read alignment to reference	Sensitivity for repetitive regions
Validation Tools	Picard, SAMtools, BCFtools	File processing and manipulation	Standardized format handling
Expression Analysis	RSEM, STAR, pheatmap	Correlating variants with expression	Quantification and visualization

Accurate variant calling in paralog-rich gene families like NBS disease resistance genes requires specialized approaches throughout the sequencing and analysis workflow. By implementing the optimized experimental and computational methods described here, researchers can overcome the limitations of standard variant calling and uncover previously undetectable genetic variation. This enables more comprehensive studies of genetic diversity in NBS genes of tolerant accessions, ultimately facilitating the identification of causal variants contributing to disease resistance in plants. The continued development of long-read sequencing technologies and specialized algorithms will further enhance our ability to resolve complex genomic regions and advance plant genetic research.

Strategies for Resolving Complex NBS Loci in Polyploid Genomes

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes constitute the largest family of plant disease resistance (R) genes and play a critical role in effector-triggered immunity [53] [24]. However, characterizing these genes in polyploid genomes presents significant challenges due to high sequence similarity between duplicated homoeologous chromosomes and the clustered arrangement of NBS-LRR genes [66] [67]. Polyploidization events, common in plant evolution, generate multiple genomic copies that undergo asymmetric gene loss, fractionation, and rearrangement [68]. This technical complexity obstructs the identification of specific NBS alleles associated with disease resistance in polyploid crops, limiting the application of molecular breeding strategies. This Application Note details integrated experimental and computational strategies to overcome these hurdles, with particular emphasis on their application in identifying genetic variations in tolerant plant accessions.

Key Challenges in Polyploid NBS Genomic Analysis

Structural Complexity and Sequence Divergence

Polyploid genomes contain multiple sub-genomes with high sequence conservation between homoeologous chromosomes. NBS-LRR genes are often arranged in complex clusters with numerous paralogs and truncated genes, complicating their assembly and annotation [66] [15]. In wheat, for example, the three sub-genomes exhibit such high similarity that sequencing reads cannot be uniquely mapped to their specific sub-genome of origin without specialized approaches [66].

Reference Genome Limitations

Most available reference genomes for polyploid species are incomplete or represent collapsed assemblies that fail to distinguish between sub-genomes [66]. This is problematic for NBS gene analysis because reference inaccuracies can lead to misinterpretation of gene copy number variation (CNV) and failure to identify rare resistance alleles. In sugarcane, the high ploidy level and aneuploidy have made traditional NLR gene prediction particularly challenging [67].

Differential Gene Retention and Selection Pressure

Following polyploidization events, NBS genes undergo asymmetric evolution between sub-genomes. Studies in Asian and European pears have revealed differential retention of NBS gene types, with 15.79% of orthologous gene pairs showing evidence of positive selection (Ka/Ks >1) [5]. This divergent evolution creates challenges for comprehensive NBS gene cataloging across sub-genomes and requires specialized analytical approaches.

Table 1: NBS-LRR Gene Diversity in Selected Polyploid and Diploid Plants

Plant Species	Ploidy	Total NBS Genes	NBS with LRR Domain	Key Features	Reference
Triticum aestivum (Bread Wheat)	Hexaploid	~2,012	Not specified	Genes distributed across A, B, D sub-genomes	[24]
Pyrus bretschneideri (Asian Pear)	Diploid	338	74.0%	NBS-LRR class most frequent (36.4%)	[5]
Pyrus communis (European Pear)	Diploid	412	55.6%	Differential expansion due to proximal duplications	[5]
Vernicia montana (Tung Tree)	Diploid	149	16.1%	Contains TIR-NBS-LRR genes; Resistant to Fusarium wilt	[53]
Vernicia fordii (Tung Tree)	Diploid	90	26.7%	Lacks TIR-NBS-LRR genes; Susceptible to Fusarium wilt	[53]
Solanum tuberosum (Potato)	Tetraploid	587 domains	Not specified	NBS domains organized in complex clusters	[15]

Integrated Strategy for NBS Loci Resolution

Experimental Design and Sequencing Approaches

A multi-faceted sequencing strategy is essential for comprehensive NBS gene resolution in polyploids. The following workflow illustrates the integrated experimental and computational approach:

Advanced Sequencing Technologies

Combining multiple long-read technologies significantly improves assembly continuity across complex NBS regions. The Human Genome Structural Variation Consortium demonstrated that integrating PacBio HiFi reads (~18 kb length with high accuracy) and ultra-long Oxford Nanopore Technologies (ONT) reads (>100 kb length) enabled the closure of 92% of assembly gaps and achieved telomere-to-telomere status for 39% of chromosomes [69]. This approach is directly applicable to complex plant genomes, particularly for resolving repetitive NBS-LRR clusters.

Targeted Enrichment Strategies

When whole-genome sequencing is cost-prohibitive for large polyploid genomes, targeted enrichment of NBS domains provides a cost-effective alternative. A method successfully applied in wheat uses a 110 MB NimbleGen SeqCap EZ gene capture probe set to enrich genic regions prior to sequencing, effectively reducing genomic complexity [66]. For specialized NBS profiling, a PCR-based approach using primers targeting conserved NBS motifs (P-loop, Kinase-2, and GLPL) can generate "NBS tags" for high-throughput sequencing [15]. This method has efficiently captured 587 NBS domains from potato cultivars.

Computational and Bioinformatics Strategies

Specialized Pipelines for Polyploid Genomes

The DaapNLRSeek (Diploidy-Assisted Annotation of Polyploid NLRs) pipeline has been specifically developed to address NLR gene prediction challenges in complex polyploid genomes like sugarcane [67]. This approach leverages syntenic relationships with diploid progenitors to improve annotation accuracy in polyploids, enabling identification of complete NLR genes, including paired NLRs and TIR-only genes that are often missed by standard annotation pipelines.

Pseudo-Chromosome Construction for Mapping

For species without complete reference genomes, constructing synteny-based pseudo-chromosomes enables effective mapping of NBS loci. This approach, successfully implemented in wheat, involves using syntenic relationships with related species (e.g., Brachypodium) to infer long-range gene order [66]. These pseudo-chromosomes serve as reference sequences for sliding window mapping-by-sequencing analyses, effectively bypassing the challenges of non-uniquely mapping reads between highly similar sub-genomes.

Orthogroup Analysis and Evolutionary Classification

Comparative phylogenomics across multiple accessions helps identify core and lineage-specific NBS genes. A pan-genome analysis of 34 plant species identified 12,820 NBS-domain-containing genes classified into 168 domain architecture classes [24]. This analysis revealed 603 orthogroups (OGs), with some core OGs (e.g., OG0, OG1, OG2) conserved across species and others specific to particular lineages. Such classification enables researchers to prioritize NBS genes with potential functional significance in disease resistance.

Table 2: Research Reagent Solutions for NBS Loci Analysis

Reagent/Resource	Function	Application Example	Considerations
NimbleGen SeqCap EZ Gene Capture Probes (110 MB)	Enrichment of genic regions	Reducing complexity for wheat genome sequencing [66]	Effectively eliminates repetitive sequences
NBS-Targeting PCR Primers (P-loop, Kinase-2, GLPL)	Amplification of NBS domains	NBS profiling in potato germplasm [15]	Targets conserved motifs; designs require multiple sequence alignment
DaapNLRSeek Pipeline	NLR gene annotation in polyploids	Sugarcane NLR annotation [67]	Leverages synteny with diploid progenitors
Verkko Assembler	Haplotype-resolved genome assembly	Human genome assembly [69]	Combines HiFi and ultra-long reads for continuity
OrthoFinder with MCL Algorithm	Orthogroup clustering	Evolutionary analysis of NBS genes across species [24]	Identifies core and lineage-specific NBS orthologs

Detailed Experimental Protocols

Protocol 1: NBS Domain Enrichment and Sequencing

Primer Design for NBS Domain Amplification

• Reference Sequence Collection: Compile NBS domain sequences from annotated R genes in a related reference genome. In potato, this involved extracting 435 distinct NBS sequences from DM genome v3.4 [15].
• Multiple Sequence Alignment: Perform ClustalW or MAFFT alignment of NBS sequences to identify conserved motif regions.
• Degenerate Primer Design: Design primers complementary to P-loop, Kinase-2, and GLPL motifs with degeneracy at polymorphic positions to maximize coverage of NBS diversity.
• Primer Validation: Test primer functionality using PCR on genomic DNA templates; select primers yielding abundant, specific amplicons.

Library Preparation and Sequencing

• DNA Extraction: Isolate high-molecular-weight genomic DNA using CTAB method with RNAse treatment [70].
• NBS Tag Amplification: Perform PCR using validated NBS-targeting primers with Illumina adapter sequences.
• Library Quality Control: Verify amplicon size distribution (200-480 bp expected) using Bioanalyzer or TapeStation.
• High-Throughput Sequencing: Pool libraries and sequence on Illumina platform (e.g., HiSeq 2500 or HiSeq X Ten) with 150 bp paired-end reads [15].

Protocol 2: Mapping-by-Sequencing in Polyploids

Reference Preparation

• Pseudo-Chromosome Construction: If complete chromosome-level reference is unavailable, create pseudo-chromosomes using syntenic relationships with related species (e.g., using Brachypodium for wheat) [66].
• Synteny Mapping: Use MCscan or similar tools to establish collinear blocks between target species and reference.

Read Mapping and Variant Calling

• Quality Control: Process raw reads with FastQC and trim adapters using Trimmomatic.
• Read Mapping: Align sequences to reference using BWA-MEM or Bowtie2 with sensitive parameters.
• Variant Calling: Identify SNPs and indels using GATK HaplotypeCaller or SAMtools mpileup.
• Sliding Window Analysis: Implement custom pipeline (e.g., in iPlant environment) to scan across reference sequence and identify regions of homozygosity in mapping populations [66].

Protocol 3: Functional Validation of Candidate NBS Genes

Virus-Induced Gene Silencing (VIGS)

• Candidate Gene Selection: Select NBS genes showing differential expression or variation between resistant and susceptible accessions.
• VIGS Construct Design: Clone 200-300 bp fragment of target NBS gene into TRV2 vector.
• Agroinfiltration: Infiltrate plant tissues with Agrobacterium tumefaciens containing TRV1 and TRV2-NBS constructs.
• Phenotypic Assessment: Challenge silenced plants with pathogen and assess disease symptoms compared to controls [53] [24].

Expression Analysis

• RNA Extraction: Isolate total RNA from pathogen-infected and mock-treated tissues at multiple time points.
• qRT-PCR: Design gene-specific primers outside the VIGS target region and perform quantitative PCR with reference genes.
• Pathogen Quantification: Measure pathogen biomass using species-specific qPCR assays to correlate NBS expression with resistance [53].

Applications in Tolerant Accession Research

Identification of Resistance-Associated Variants

Comparative analysis of NBS genes between resistant and susceptible accessions can reveal structural variants associated with disease resistance. In tung trees, comparison between Fusarium wilt-resistant Vernicia montana and susceptible V. fordii identified 149 and 90 NBS-LRR genes, respectively, with specific TIR-NBS-LRR genes present only in the resistant species [53]. Functional validation confirmed that one candidate gene (Vm019719) conferred resistance to Fusarium wilt, while its allelic counterpart in susceptible V. fordii contained a promoter deletion that rendered it ineffective.

Domestication History and Selection Patterns

Analysis of NBS gene diversity in wild and domesticated accessions can reveal selection signatures during domestication. In Asian and European pears, domestication caused contrasting patterns: Asian cultivars showed decreased nucleotide diversity in NBS genes (6.23E-03 in cultivated vs. 6.47E-03 in wild), while European pears showed increased diversity (6.48E-03 in cultivated vs. 5.91E-03 in wild) [5]. This suggests independent domestication histories with different impacts on the NBS gene repertoire.

Pan-NLRome Construction for Breeding

The integration of multiple genome assemblies enables the construction of a pan-NLRome representing the full diversity of NBS-LRR genes within a species. In potato, screening of 91 cultivars identified an average of 26 nucleotide polymorphisms per NBS domain, providing a rich resource for marker development [15]. Similar approaches in cotton have identified orthogroups (e.g., OG2, OG6, OG15) with upregulated expression in tolerant accessions under biotic stress, highlighting potential targets for marker-assisted selection [24].

Resolving complex NBS loci in polyploid genomes requires an integrated approach combining advanced sequencing technologies, specialized computational pipelines, and systematic functional validation. The strategies outlined in this Application Note provide a roadmap for characterizing the extensive genetic variation in NBS genes of tolerant accessions, ultimately accelerating the development of disease-resistant crop varieties through molecular breeding. As genomic technologies continue to advance, the resolution of complex NBS loci will become increasingly routine, enabling more comprehensive understanding of plant immunity systems in even the most challenging polyploid genomes.

Overcoming Reference Bias in Diverse Germplasm and Wild Accessions

Reference bias is a significant drawback in genomic analyses where alignment tools systematically miss or incorrectly report alignments for sequencing reads containing non-reference alleles [71]. This bias confounds measurements and leads to incorrect results, particularly affecting studies on hypervariable regions, allele-specific effects, and epigenomic signals [71]. In the context of analyzing Nucleotide-Binding Site-Leucine Rich Repeat (NBS-LRR) genes—the primary disease resistance genes in plants—reference bias poses a substantial challenge when working with diverse germplasm collections and wild accessions that contain extensive genetic variation not present in reference genomes [24] [72].

The issue is particularly acute for plant genetic resources (PGR) maintained in genebanks worldwide, which include landraces, crop wild relatives, and mutants harboring valuable alleles for breeding [73]. Over 7 million accessions are held in more than 1,700 genebanks globally, with collections for agriculturally important species like wheat, rice, and barley comprising hundreds of thousands of accessions each [73]. Traditional alignment to single reference genomes impedes the comprehensive analysis of this diversity, especially for rapidly evolving gene families like NBS-LRRs that show remarkable structural variation across accessions [24] [72].

Quantitative Assessment of Reference Bias

Metrics for Measuring Bias

Biastools provides a framework for categorizing and measuring reference bias through several balance metrics calculated at heterozygous variant sites [71]. The key metrics include:

Table 1: Metrics for Quantifying Reference Bias

Metric	Definition	Interpretation
Simulation Balance (SB)	Proportion of simulated reads overlapping a heterozygous site that originated from the reference-carrying haplotype	Baseline expectation from simulation
Mapping Balance (MB)	Allelic balance considering only reads that truly originated from and overlapped the heterozygous site after alignment	Measures bias introduced during read mapping
Assignment Balance (AB)	Allelic balance after using an algorithm to determine haplotype of origin for each read overlapping the site	Measures bias introduced during haplotype assignment
Normalized Mapping Balance (NMB)	NMB = MB - SB	Values >0 indicate mapping bias toward reference allele
Normalized Assignment Balance (NAB)	NAB = AB - SB	Values >0 indicate assignment bias toward reference allele

Categorizing Bias Events

Using these metrics, Biastools categorizes reference bias into distinct types [71]:

Table 2: Categories of Reference Bias

Bias Type	NMB-NAB Signature	Primary Cause
Loss Bias	Points along diagonal in upper-right quadrant	Systematic failure to align reads containing alternate alleles
Flux Bias	Vertically displaced from origin with near-zero NMB	Reads with low mapping quality placed incorrectly across repeat regions
Local Bias	Vertically displaced from origin with high mapping quality reads	Assignment step errors, often at short tandem repeats

Experimental data shows that aligners favoring local alignments with soft clipping exhibit more bias around gaps, while end-to-end alignment modes reduce bias at insertions and deletions [71]. Approximately 79% of local bias events occur at sites annotated by Repeatmasker, with 1,012 sites in simple repeats [71].

Computational Strategies to Overcome Reference Bias

Pangenome Graph Approaches

Pangenome graphs that incorporate known genetic variants from multiple assemblies significantly reduce reference bias by removing alignment penalties for known alternate alleles [71] [32]. The frequency of structural variations (SVs) varies between subgenomes, as demonstrated in peanut where SV frequency in subgenome A is higher than in subgenome B [32].

Table 3: Pangenome Construction Statistics from Recent Studies

Species	Number of Genomes	Core Gene Families	Distributed Gene Families	Private Gene Families
Peanut (Arachis hypogaea)	8	17,137	22,232	5,643 [32]
Angiosperms (NLR genes only)	304	Not specified	Not specified	~90,000 NLR genes total [24]
34 Plant Species (NBS genes only)	34	168 domain architecture classes	Species-specific structural patterns	603 orthogroups [24]

Bias-Aware Alignment Workflow

The following workflow illustrates a comprehensive approach to minimizing reference bias in NBS gene analysis:

Figure 1: Comprehensive workflow for overcoming reference bias in NBS gene analysis of diverse germplasm.

Experimental Protocol for NBS Gene Analysis in Diverse Germplasm

Genome Assembly and Pangenome Construction

Materials:

• Diverse germplasm accessions representing target species
• High-molecular-weight DNA (>50 kb)
• PacBio HiFi, Oxford Nanopore, and Hi-C sequencing reagents

Methodology:

• Sequencing: Perform deep sequencing using PacBio HiFi long reads (N50 > 19 kb), Oxford Nanopore ultra-long reads (N50 > 70 kb), and generate Hi-C data for chromosome scaffolding [74] [32].
• Assembly: Conduct de novo assembly using NextDenovo or similar tools, followed by chromosome-level scaffolding with HiC-Pro or comparable software [32].
• Annotation: Identify repetitive elements using Extensive de novo repeat annotation followed by protein-coding gene prediction integrating RNA-seq evidence, ab initio predictions, and homology-based searches [74].
• Pangenome Construction: Cluster genes into orthogroups using OrthoFinder with DIAMOND for sequence similarity searches and MCL for clustering [24]. Classify gene families as core (present in all accessions), soft-core (missing in few), distributed (variable presence), or private (accession-specific) [32].

NBS-LRR Gene Identification and Classification

Materials:

• Annotated genome assemblies from multiple accessions
• High-performance computing resources
• HMMER suite and custom hidden Markov models

Methodology:

• Domain Identification: Scan all predicted proteins using HMMER with the NB-ARC domain (PF00931) hidden Markov model (E-value < 1×10⁻²⁰) [72]. Confirm hits through manual curation and verification of intact NBS domains.
• Architecture Classification: Identify associated domains (TIR, CC, LRR, RPW8) using hmmpfam against Pfam databases. Use Paircoil2 (P-score cutoff 0.03) for coiled-coil domain prediction [72].
• Orthogroup Analysis: Perform phylogenetic analysis of NB-ARC domains using Maximum Likelihood method in MEGA6 or similar software. Cluster sequences into orthogroups using OrthoFinder with 1000 bootstrap replicates [24].
• Genetic Variation Mapping: Call variants using GATK or similar pipelines with strict filtering (MAF > 0.05, missing rate < 0.5). Identify structural variations using assembly-based methods and graph-based detection [32].

Bias Measurement and Validation

Materials:

• Simulated and empirical sequencing data from heterozygous accessions
• Biastools software package
• Multiple reference genomes and pangenome graphs

Methodology:

• Simulation Experiments: Use Biastools simulate mode with known variant files (VCF) from high-quality projects like GIAB to generate synthetic reads from both haplotypes (~15× coverage each) [71].
• Alignment Comparison: Align simulated reads to both linear references and pangenome graphs using Bowtie 2, BWA-MEM, Minimap2, and VG Giraffe with consistent parameters [71].
• Bias Quantification: Calculate Simulation Balance (SB), Mapping Balance (MB), and Assignment Balance (AB) for all heterozygous sites. Compute Normalized Mapping Balance (NMB) and Normalized Assignment Balance (NAB) to categorize bias events [71].
• Functional Validation: Select NBS genes showing significant diversity for experimental validation using virus-induced gene silencing (VIGS) in resistant accessions and transgenic complementation in susceptible backgrounds [24].

Research Reagent Solutions

Table 4: Essential Research Reagents and Tools for Bias-Aware NBS Gene Analysis

Category	Specific Tools/Reagents	Function/Application
Sequencing Technologies	PacBio HiFi, Oxford Nanopore Ultra-long, Hi-C	De novo genome assembly and structural variant detection [74] [32]
Alignment Tools	VG Giraffe, HISAT2, Bowtie 2 (end-to-end mode), BWA-MEM	Read mapping to linear references and pangenome graphs [71]
Bias Measurement	Biastools (simulate, predict, scan modes)	Quantifying and categorizing reference bias events [71]
Variant Callers	GATK, bcftools, graph-based variant callers	SNP and structural variation discovery [32]
NBS Gene Annotation	HMMER (Pfam NB-ARC domain), OrthoFinder, MEME	Identification and classification of resistance gene families [24] [72]
Functional Validation	VIGS vectors, Agrobacterium-mediated transformation, CRISPR-Cas9	Functional characterization of NBS gene alleles [24]

Case Study: NBS Gene Analysis in Cotton Accessions

A recent study on Gossypium hirsutum accessions provides a practical example of this approach [24]. Researchers identified 12,820 NBS-domain-containing genes across 34 plant species, classifying them into 168 domain architecture classes. The analysis revealed several species-specific structural patterns including TIR-NBS-TIR-Cupin1-Cupin1 and TIR-NBS-Prenyltransf architectures [24].

Expression profiling demonstrated upregulation of specific orthogroups (OG2, OG6, OG15) in different tissues under various biotic and abiotic stresses in cotton accessions with contrasting responses to cotton leaf curl disease (CLCuD) [24]. Genetic variation analysis between susceptible (Coker 312) and tolerant (Mac7) accessions identified 6,583 unique variants in Mac7 and 5,173 in Coker312 within NBS genes [24]. Functional validation through silencing of GaNBS (OG2) in resistant cotton demonstrated its role in virus tolerance, confirming the value of comprehensive bias-aware analysis [24].

Overcoming reference bias is essential for fully leveraging the genetic diversity present in germplasm collections and wild accessions. Pangenome graph approaches combined with systematic bias measurement using tools like Biastools enable more comprehensive characterization of NBS gene diversity, revealing valuable alleles for crop improvement. As reference bias affects all genomic analyses to varying degrees, incorporating these strategies into standard practice will accelerate the discovery of genetic elements underlying stress tolerance and disease resistance in plants.

Best Practices for Functional Annotation of Novel NBS Gene Variants

Nucleotide-binding site (NBS) genes constitute one of the largest and most critical gene families in plant innate immunity, encoding proteins that function as major immune receptors for effector-triged immunity (ETI) [24]. These genes, often characterized by a conserved NBS domain alongside leucine-rich repeats (LRRs) and variable N-terminal domains (TIR, CC, or RPW8), are notoriously polymorphic at the population level and exhibit complex evolutionary dynamics driven by duplication events and selective pressures [75] [24]. In the context of researching tolerant plant accessions, the precise functional annotation of novel NBS gene variants is a critical step for unraveling the genetic basis of disease resistance. This process translates raw sequence variations into testable hypotheses about gene function, enabling researchers to connect genetic differences to resilient phenotypes. This protocol outlines a comprehensive framework for the functional annotation of novel NBS variants, integrating current bioinformatic tools and evolutionary principles to guide research in plant immunity and resistance breeding.

The Biological Significance of NBS Genes in Plant Immunity

NBS-LRR genes (NLRs) are the predominant class of plant resistance (R) genes and play an indispensable role in pathogen surveillance [24] [76]. They function as modular intracellular immune receptors that recognize specific pathogen effector proteins, triggering a robust defense response often accompanied by programmed cell death [75]. This gene family is divided into major subclasses based on the N-terminal domain:

• TNLs: Characterized by a Toll/Interleukin-1 Receptor (TIR) domain.
• CNLs: Characterized by a Coiled-Coil (CC) domain.
• RNLs: Characterized by an RPW8 domain and often acting as helpers in signaling cascades [24] [76].

A key feature of the NBS gene family is its remarkable diversification. The number of NLR genes can vary substantially across plant species, from fewer than 100 to over 1,000, and they are often organized in complex genomic clusters [75] [24]. This diversification is driven by various evolutionary mechanisms, with tandem and segmental duplications playing a major role, followed by positive selection that shapes the antigen recognition surfaces of the proteins [24] [76]. Furthermore, due to the significant fitness costs associated with constitutive expression, plants have evolved sophisticated regulatory mechanisms, including post-transcriptional control by diverse microRNA (miRNA) families that target conserved NBS-LRR motifs [75]. This intricate evolutionary and regulatory landscape makes the functional annotation of NBS variants a specialized and critical endeavor.

Table 1: Key Characteristics of Major NBS-LRR Subfamilies

Subfamily	N-Terminal Domain	Presence in Monocots/Dicots	General Function	Evolutionary Notes
TNL (TIR-NBS-LRR)	TIR (Toll/Interleukin-1 Receptor)	Dicots only	Sensor; effector recognition & signaling	Often evolves under strong positive selection [76].
CNL (CC-NBS-LRR)	CC (Coiled-Coil)	Both monocots and dicots	Sensor; effector recognition & signaling	The most prevalent subclass in many plants [24].
RNL (RPW8-NBS-LRR)	RPW8 (Resistance to Powdery Mildew 8)	Both monocots and dicots	Helper; amplifies defense signals	More conserved; can be activated by other NLRs [76].

A Step-by-Step Annotation Protocol for Novel NBS Variants

Stage 1: Variant Calling and Quality Control

Objective: To identify high-confidence sequence variants from whole-genome sequencing data of tolerant and susceptible accessions.

Materials & Protocols:

• Input Data: Whole-genome sequencing (WGS) data (Illumina short-read or PacBio/Oxford Nanopore long-read) from paired plant accessions (e.g., tolerant vs. susceptible) aligned to a reference genome. For human genetics, the recommended reference is GRCh38/hg38 [77]. For plants, use the most recent and complete assembly available.
• Variant Calling: Use a standard germline variant calling pipeline, such as the GATK Best Practices workflow. The output is a Variant Calling Format (VCF) file.
• Quality Control (QC):
  • Filtering: Apply hard filters (e.g., QD < 2.0 || FS > 60.0 || MQ < 40.0) or use variant quality score recalibration (VQSR) to remove low-quality calls.
  • VCF Annotation (Preliminary): Use tools like bcftools csq to perform basic consequence prediction (e.g., synonymous, missense, frameshift) [78].

Deliverable: A high-confidence VCF file containing annotated NBS gene variants for downstream analysis.

Stage 2: Comprehensive Variant Annotation and Prioritization

Objective: To annotate variants with functional information and prioritize those most likely to impact NBS gene function.

Materials & Protocols:

• Functional Annotation Pipeline: Utilize a comprehensive annotation tool like NIRVANA [79] or Ensembl VEP. These tools integrate multiple data sources.
  • Input: High-confidence VCF file from Stage 1.
  • Configuration: Specify the appropriate reference genome and transcript database (e.g., Ensembl Plants). The tool will generate a detailed annotation table.
• Key Annotations to Extract [79]:
  • Transcript Consequence: e.g., missense_variant, stop_gained, splice_region_variant.
  • Protein Change: In HGVS nomenclature (e.g., p.(Ser534Pro)). Note that descriptions should be based on an accepted reference sequence, with a letter prefix such as 'c.' for coding DNA and 'p.' for protein [77].
  • Population Frequency: Leverage fields like gnomad_all_af to filter out common polymorphisms unlikely to cause rare, large-effect phenotypes.
  • In Silico Predictions: Integrate scores from tools like SpliceAI (for splice effects) and REVEL (for missense pathogenicity) [78] [79].
• NBS-Specific Prioritization:
  • Focus on variants that alter conserved protein motifs, particularly the P-loop and other nucleotide-binding pockets within the NBS domain, as these are frequently targeted by regulatory miRNAs and are critical for function [75].
  • Prioritize non-synonymous variants in the solvent-exposed residues of the LRR domain, as these are often involved in direct effector recognition and display high diversity [75] [24].
  • For tolerant accessions, specifically look for gain-of-function variants (e.g., altered autoinhibition, enhanced ATPase activity) or variants that might modify gene regulation.

Table 2: A Tiered System for Prioritizing NBS Gene Variants

Tier	Variant Type	Rationale for Prioritization in Tolerant Accessions	Follow-up Analysis
Tier 1 (High)	Stop-gain, Frameshift, Canonical splice-site	Likely complete loss-of-function; may be relevant if susceptibility is dominant.	Check for haploinsufficiency; study in heterozygous state.
Tier 2 (Medium)	Missense in NBS domain, LRR domain	May directly affect protein function, ATP binding, or effector recognition.	Molecular dynamics simulation; co-segregation analysis.
Tier 3 (Medium)	Non-coding variants in miRNA binding sites	May disrupt miRNA-mediated repression (e.g., miR482/2118), leading to over-expression and autoimmunity [75] [78].	Dual-luciferase reporter assay to confirm disrupted miRNA interaction.
Tier 4 (Low)	Synonymous, deep intronic	Unlikely to affect protein function or splicing.	Usually deprioritized unless strong evidence links them to a regulatory function.

Stage 3: Functional Validation and Experimental Design

Objective: To experimentally verify the functional impact of prioritized NBS gene variants.

Materials & Protocols:

• Gene Expression Analysis:
  • Method: Quantitative RT-PCR (qRT-PCR) or RNA-seq on RNA extracted from pathogen-infected and mock-treated tissues of tolerant and susceptible accessions.
  • Analysis: Assess if the variant correlates with differential expression of the NBS gene or downstream defense response genes.
• Functional Assays for Immune Response:
  • Hypersensitive Response (HR) Assay: Transiently express the wild-type and mutant NBS gene variants in a model system like Nicotiana benthamiana. An accelerated or constitutive HR can indicate autoactivation and a gain-of-function phenotype [24].
  • Virus-Induced Gene Silencing (VIGS): Silencing the candidate NBS gene in a resistant/tolerant accession (e.g., GaNBS in cotton) and challenging with a pathogen can demonstrate its requirement for resistance, as shown in functional studies [24].
• Protein-Ligand Interaction Studies:
  • Method: Conduct computational protein-ligand docking simulations to assess how a missense variant affects the binding affinity for nucleotides (ADP/ATP) or, potentially, for pathogen effector proteins [24].
  • Validation: Follow up with experimental methods like Surface Plasmon Resonance (SPR) or Isothermal Titration Calorimetry (ITC) for biophysical confirmation.

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents and Resources for NBS Gene Functional Annotation

Item/Tool Name	Category	Function & Application in NBS Annotation
EggNOG-mapper [80]	Functional Annotation Tool	Assigns functional terms (GO, KEGG) by comparing protein sequences to ortholog groups. Provides initial functional hypotheses.
InterProScan [80]	Functional Annotation Tool	Identifies protein domains, motifs, and families by scanning against multiple databases. Crucial for confirming NBS, LRR, TIR, CC domains.
SpliceAI [78] [79]	In Silico Prediction	Predicts the effect of sequence variants on splicing, critical for interpreting non-coding and intronic variants.
OrthoDB [81]	Protein Sequence Database	A catalog of orthologs used by BRAKER and other tools to provide taxon-specific evidence for gene prediction and annotation.
VIGS Vector System	Functional Validation	Used for transient, post-transcriptional gene silencing to knock down candidate NBS genes and assess their role in disease resistance phenotypes [24].
Gateway Cloning System	Molecular Biology	Facilitates the rapid transfer of NBS gene coding sequences into various expression vectors for transient assays or stable transformation.

Workflow and Data Visualization

The following diagram summarizes the core bioinformatic workflow for the annotation and prioritization of novel NBS gene variants, from raw sequencing data to a shortlist of high-priority candidates.

NBS Variant Annotation and Prioritization Workflow

The systematic functional annotation of novel NBS gene variants is a powerful approach to deciphering the molecular mechanisms of disease resistance in plants. By integrating robust bioinformatic pipelines, a deep understanding of NBS gene family evolution, and targeted experimental validation, researchers can efficiently sift through genomic variation to identify causal polymorphisms. This protocol provides a concrete roadmap for such analyses, directly supporting the broader objective of identifying key genetic elements in tolerant accessions. The application of these best practices will accelerate the discovery of functional resistance genes and inform strategies for developing durable disease resistance in crops, a cornerstone of sustainable agriculture.

Validating NBS Gene Function and Comparative Analysis of Tolerant Accessions

Functional validation of candidate genes is a critical step in modern plant genomics, directly linking genetic sequences to biological traits. Within the context of researching genetic variation in NBS (Nucleotide-Binding Site-Leucine-Rich Repeat) genes from stress-tolerant accessions, robust validation methods are indispensable for confirming gene function. Virus-Induced Gene Silencing (VIGS) has emerged as a powerful, rapid technique for assessing gene function in planta, bypassing the need for stable transformation. This protocol details the application of a highly efficient Tobacco Rattle Virus (TRV)-based VIGS system for the functional analysis of NBS-LRR genes in soybean, a method that can be adapted for other dicot species. When combined with insights from stable transgenic approaches, these methods provide a comprehensive toolkit for elucidating the role of disease resistance genes in plant immunity, thereby accelerating the development of resilient crop varieties [82] [2] [83].

Application Notes: The Role of VIGS in Functional Genomics of NBS-LRR Genes

The NBS-LRR gene family constitutes a major class of plant disease resistance (R) genes, which are pivotal in triggering immune responses upon pathogen recognition [2]. In a recent study, 603 NBS-LRR members were identified in the Nicotiana tabacum genome alone, highlighting the vast genetic repertoire available for exploration [2]. Functional studies have demonstrated that silencing specific NBS-LRR genes can significantly compromise disease resistance. For instance, silencing an NBS-LRR gene in cotton reduced resistance to Verticillium dahliae, while in wheat, an allelic mutation in an NBS-LRR gene was the causal factor for premature senescence [2] [83]. These findings underscore the critical function of NBS-LRR genes and the importance of precise functional validation tools like VIGS.

Advantages of VIGS over Stable Transformation:

• Speed: Allows for rapid functional assessment within weeks, as opposed to the months or years required for developing stable transgenic lines [82].
• Flexibility: Enables the screening of multiple candidate genes identified from genome-wide association studies (GWAS) or transcriptomic analyses of tolerant accessions [82] [84].
• Bypasses Complex Genetics: Ideal for functional analysis in polyploid species or plants with recalcitrant transformation systems, such as soybean [82].

Protocols for TRV-Mediated VIGS in Soybean

This protocol, optimized from a recent 2025 study, describes an Agrobacterium-mediated VIGS method for soybean, achieving a high silencing efficiency of 65% to 95% [82].

Research Reagent Solutions

The following table lists the essential materials required for implementing the TRV-VIGS system.

Table 1: Key Research Reagents for TRV-VIGS

Reagent/Material	Function/Description	Source/Example
pTRV1 and pTRV2 Vectors	Binary TRV vectors; pTRV1 contains replication and movement proteins, pTRV2 carries the target gene insert for silencing.	[82]
Agrobacterium tumefaciens GV3101	Strain used for delivering the TRV vectors into plant cells.	[82]
Gene-Specific Primers	Designed to amplify a ~200-500 bp fragment of the target NBS-LRR gene for cloning into pTRV2.	[82]
Restriction Enzymes (EcoRI, XhoI)	Used for directional cloning of the target gene fragment into the pTRV2 vector.	[82]
Soybean Seeds	Plant material for transformation. The cultivar 'Tianlong 1' was used in the original study.	[82]
Antibiotics (Kanamycin, Rifampicin)	For selection of transformed Agrobacterium and plasmid maintenance.	[82]

Step-by-Step Experimental Workflow

The following diagram illustrates the complete experimental workflow for TRV-mediated VIGS, from vector construction to phenotypic analysis.

Detailed Methodology

Step 1: TRV2 Vector Construction

• Amplify Target Fragment: Isolate a 200-500 bp specific fragment of the target NBS-LRR gene (e.g., GmRpp6907 or GmRPT4) from soybean cDNA using gene-specific primers with engineered EcoRI and XhoI restriction sites [82].
• Clone into pTRV2: Digest the pTRV2-GFP vector and the PCR product with EcoRI and XhoI. Ligate the target fragment into the linearized pTRV2 vector [82].
• Transform and Sequence: Transform the ligation product into E. coli DH5α competent cells. Select positive clones, confirm the insert sequence via Sanger sequencing, and then introduce the validated recombinant plasmid into Agrobacterium tumefaciens strain GV3101 [82].

Step 2: Agrobacterium Culture Preparation

• Inoculate single colonies of Agrobacterium containing pTRV1 and the recombinant pTRV2 vectors separately in liquid LB media with appropriate antibiotics (e.g., Kanamycin, Rifampicin).
• Grow cultures at 28°C with shaking (200 rpm) until they reach an OD₆₀₀ of 0.8-1.0.
• Centrifuge the cultures and resuspend the pellets in an induction medium (e.g., containing 10 mM MES, 10 mM MgClâ‚‚, and 200 µM acetosyringone).
• Incubate the suspensions for 3-4 hours at room temperature.
• Mix the pTRV1 and pTRV2 cultures in a 1:1 ratio before infection [82].

Step 3: Plant Material Preparation and Agroinfiltration

• Critical Optimization: Conventional infiltration methods (e.g., leaf injection) are inefficient in soybean due to a thick cuticle and dense trichomes. The optimized protocol uses cotyledon nodes for high-efficiency transformation [82].
• Surface-sterilize soybean seeds and imbibe them in sterile water until swollen.
• Longitudinally bisect the swollen seeds to create half-seed explants, ensuring the cotyledonary node is exposed [82].
• Immerse the fresh explants in the prepared Agrobacterium suspension for 20-30 minutes, with gentle agitation [82].
• Co-cultivate the infected explants on sterile filter paper in the dark for 2-3 days.

Step 4: Plant Growth and Monitoring

• Transfer the explants to a standard tissue culture medium or soil and maintain them in a growth chamber with a 16/8 h light/dark cycle at 22-25°C.
• Monitor for systemic silencing symptoms. For a positive control like GmPDS (phytoene desaturase), photobleaching of leaves typically appears within 14-21 days post-inoculation (dpi) [82].

Validation and Phenotypic Analysis

Efficiency Assessment:

• Fluorescence Check: At 4 dpi, examine the cotyledonary nodes under a fluorescence microscope for GFP signals to confirm successful infection. The optimized protocol achieves an effective infectivity efficiency of over 80% [82].
• Molecular Validation:
  • Quantitative PCR (qPCR): Quantify the transcript levels of the silenced target NBS-LRR gene in newly emerged systemic leaves. Compare the levels to plants inoculated with an empty TRV vector (pTRV:empty) [82].
  • Calculate silencing efficiency using the 2^(-ΔΔCt) method. A significant reduction (e.g., >70%) in transcript abundance confirms successful silencing.

Phenotypic Evaluation:

• For NBS-LRR genes involved in disease resistance, challenge the silenced plants with the corresponding pathogen (e.g., soybean rust for GmRpp6907).
• Document the loss-of-resistance phenotype, such as increased disease susceptibility, compared to control plants [82]. Correlate the molecular data (reduced transcript levels) with the observed phenotypic changes to confirm gene function.

Table 2: Quantitative Results from an Optimized TRV-VIGS Study in Soybean

Target Gene	Gene Function	Silencing Efficiency	Observed Phenotype	Key Application
GmPDS	Carotenoid biosynthesis	65% - 95%	Photobleaching (white leaves)	Positive control for silencing [82]
GmRpp6907	Rust resistance	High	Compromised rust immunity	Validated disease resistance function [82]
GmRPT4	Defense-related	High	Induced significant phenotypic changes	Confirmed role in plant defense [82]

Integration with Broader Research on Genetic Variation in NBS Genes

The VIGS protocol is exceptionally powerful when deployed within a research pipeline that begins with the identification of genetic variation. Genome-wide association studies (GWAS) in diverse populations, such as the Aegilops tauschii-derived wheat population, can pinpoint specific marker-trait associations (MTAs) for stress resilience [84]. Candidate genes located near these MTAs, particularly NBS-LRR genes, become prime targets for functional validation using the VIGS system described above.

The structure of a typical NBS-LRR gene and its functional domains, which are the target of such analyses, is shown below.

Diagram 2: Key domains of an NBS-LRR disease resistance gene. The TIR/CC domain is often involved in signaling, the NBS domain is responsible for nucleotide binding and activation, and the LRR domain is responsible for specific pathogen recognition [2].

This combined approach—from population-level genetic analysis to rapid gene-level functional validation—creates a robust framework for discovering and characterizing key genetic players in stress tolerance, ultimately contributing to the development of climate-resilient crops [82] [84].

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes represent the largest class of disease resistance (R) genes in plants, playing a crucial role in effector-triggered immunity (ETI) against various pathogens, including viruses [24] [85]. In cotton (Gossypium spp.), Cotton Leaf Curl Disease (CLCuD), caused by begomoviruses in the Geminiviridae family, is a devastating constraint to production, transmitted by the whitefly (Bemisia tabaci) vector [24] [86]. This case study, situated within a broader thesis investigating genetic variation in NBS genes of tolerant accessions, explores the functional validation of a specific NBS gene via virus-induced gene silencing (VIGS) and its direct impact on viral titre in a resistant cotton host.

Background and Key Findings

A comprehensive genome-wide study identified 12,820 NBS-domain-containing genes across 34 plant species, revealing significant diversity and several species-specific architectural patterns [24]. Orthogroup (OG) analysis clustered these genes into 603 groups, with expression profiling indicating the putative upregulation of OG2, OG6, and OG15 under various biotic stresses in cotton accessions with differing susceptibilities to CLCuD [24].

Genetic variation analysis between a susceptible (Coker 312) and a tolerant (Mac7) accession of Gossypium hirsutum identified 6,583 unique variants in the NBS genes of the tolerant Mac7 and 5,173 in the susceptible Coker 312 [24]. This genetic divergence underscores the potential for functional differences in NBS-mediated resistance mechanisms. The core experimental finding demonstrated that silencing of GaNBS (a candidate gene from OG2) in a resistant cotton background through VIGS compromised plant defense, confirming its critical role in limiting virus accumulation [24].

Table 1: Key Quantitative Findings from Comparative NBS Gene Analysis

Analysis Parameter	Finding	Significance
Total NBS Genes Identified	12,820 genes across 34 species [24]	Highlights the extensive diversity of this gene family.
Orthogroups (OGs) with Tandem Duplications	603 OGs identified [24]	Indicates duplication as a key evolutionary mechanism.
Putative Upregulated OGs under Stress	OG2, OG6, OG15 [24]	Pinpoints candidate OGs involved in stress response.
Unique NBS Gene Variants in Tolerant Mac7	6,583 variants [24]	Suggests a genetic basis for tolerance in specific accessions.
Unique NBS Gene Variants in Susceptible Coker 312	5,173 variants [24]	Provides a comparative baseline for genetic variation.

Experimental Protocol: VIGS for Functional Validation of NBS Genes

This protocol details the methodology for silencing a candidate NBS gene (e.g., GaNBS from OG2) in resistant cotton to assess its impact on virus titre, based on established VIGS procedures [24].

Principle

VIGS is a reverse genetics technique that uses a recombinant viral vector to deliver a host gene fragment, triggering post-transcriptional gene silencing (PTGS) of the corresponding endogenous mRNA. This allows for the functional analysis of a gene by observing the phenotypic consequences of its knockdown [86].

Materials

• Plant Material: Seeds of resistant cotton accession (e.g., Gossypium arboreum or Mac7 G. hirsutum).
• VIGS Vector: A Tobacco rattle virus (TRV)-based vector (e.g., TRV1 and TRV2).
• Target Gene Fragment: A ~200-300 bp fragment specific to the candidate NBS gene (e.g., GaNBS), cloned into the TRV2 vector in sense/antisense orientation to form a hairpin [86].
• Agrobacterium tumefaciens: Strain GV3101.
• Growth Media: Luria-Bertani (LB) broth and agar with appropriate antibiotics (kanamycin, rifampicin).
• Infiltration Buffer: 10 mM MES, 10 mM MgClâ‚‚, 200 µM acetosyringone, pH 5.6.

Procedure

• Vector Preparation: Clone the specific fragment of the target NBS gene into the TRV2 vector. Verify the construct by sequencing and restriction digestion [86].
• Agrobacterium Transformation: Introduce the recombinant TRV2 and the helper TRV1 plasmids into Agrobacterium separately. Select positive colonies on LB plates with antibiotics and confirm via colony PCR [86].
• Agrobacterium Culture Preparation:
  • Inoculate single colonies of Agrobacterium containing TRV1 and the recombinant TRV2 into 5 mL of LB broth with antibiotics. Incubate at 28°C overnight with shaking.
  • Sub-culture the bacteria into 50 mL of fresh LB broth with antibiotics and 10 mM MES. Grow until OD₆₀₀ reaches ~1.0.
  • Pellet the cells by centrifugation and resuspend in infiltration buffer to a final OD₆₀₀ of 1.5.
  • Incubate the suspensions at room temperature for 3-4 hours.
• Inoculum Mixing: Combine the TRV1 and recombinant TRV2 suspensions in a 1:1 ratio.
• Plant Infiltration:
  • Use cotton plants at the two-leaf stage.
  • Using a needle-less syringe, infiltrate the mixed Agrobacterium culture into the abaxial side of the cotyledons or true leaves.
  • Include control plants infiltrated with a TRV2 vector containing a non-functional insert (e.g., empty vector or a fragment from an unrelated gene).
• Post-Infiltration Care: Maintain infiltrated plants in a greenhouse or growth chamber at 22-25°C with a 16/8-hour light/dark cycle. High humidity should be maintained for 24-48 hours post-infiltration to facilitate infection.
• Silencing Validation and Phenotyping:
  • After 2-3 weeks, observe the emergence of a silencing phenotype in positive control plants (e.g., TRV2::PDS causing photobleaching).
  • Harvest leaf tissue from the silenced areas of experimental and control plants.
  • Molecular Confirmation: Use quantitative RT-PCR (qRT-PCR) to quantify the transcript level of the target NBS gene, confirming significant knockdown in plants infiltrated with TRV2::GaNBS compared to controls.
• Virus Challenge and Titre Quantification:
  • Challenge the silenced and control plants with CLCuD via viruliferous whiteflies [87] or agro-inoculation.
  • After a set period (e.g., 14-21 days), collect leaf tissue showing symptoms.
  • Virus Titre Quantification: Extract total DNA and use quantitative PCR (qPCR) with primers specific to a conserved region of the cotton leaf curl virus (e.g., coat protein gene) to measure viral DNA accumulation.

Data Analysis

Compare the viral titre between the GaNBS-silenced plants and the control plants. A statistically significant increase in virus titre in the silenced plants demonstrates that the NBS gene is necessary for limiting viral replication, thereby confirming its functional role in resistance [24].

Signaling Pathways in NBS-Mediated Disease Resistance

The following diagram illustrates the generalized signaling pathway activated by CNL-type R proteins, which is implicated in resistance against pathogens like Verticillium dahliae and viruses, and how its disruption via VIGS leads to susceptibility.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Research Reagents for Functional NBS Gene Analysis

Reagent / Solution	Function / Application	Key Considerations
TRV-based VIGS Vectors	To induce targeted silencing of candidate NBS genes in planta.	Ensures efficient and transient knockdown; requires careful fragment selection to ensure specificity [24] [86].
Agrobacterium tumefaciens (GV3101)	A biological vector for delivering the VIGS construct into plant cells.	Preferred for high transformation efficiency in cotton; requires acetosyringone in the infiltration buffer [86].
Virus-Specific Primers for qPCR	To accurately quantify viral load (e.g., CLCuV titre) in challenged plants.	Targets a conserved viral region (e.g., coat protein); essential for measuring the functional outcome of silencing [24] [88].
Gene-Specific Primers for qRT-PCR	To validate the knockdown efficiency of the target NBS gene post-VIGS.	Confirms the molecular efficacy of the silencing protocol before phenotyping [24].
NBS Domain HMM Profile (PF00931)	For the genome-wide identification and annotation of NBS-encoding genes.	Found in the Pfam database; used with HMMER software for initial gene discovery [85] [89].

This application note demonstrates a validated protocol for functionally characterizing NBS genes in cotton resistance to CLCuD. The key finding that silencing of GaNBS (OG2) increases virus titre provides direct evidence of its indispensable role in the plant's defense system [24]. Integrating this functional validation with data on genetic variation in tolerant accessions like Mac7 offers a powerful, comprehensive strategy for identifying and deploying robust R genes in cotton breeding programs, ultimately contributing to the sustainable management of devastating viral diseases.

Application Note: Bitopic Nanobody-Ligand Conjugates for Targeted GPCR Signaling

G protein-coupled receptors (GPCRs) represent the largest family of plasma membrane-embedded signaling proteins and are involved in a wide array of physiological processes, making them attractive targets for drug development [90]. A significant challenge in GPCR pharmacology is achieving receptor subtype specificity and tissue-selective activation, as orthosteric binding sites are often highly conserved across receptor subtypes. This application note details a novel chemical biology approach using bitopic nanobody-ligand conjugates to achieve logic-gated signaling, wherein activation occurs only when two distinct receptors are co-expressed in the same cell population [90]. This methodology offers a path towards tissue-specific pharmacology with implications for GPCR drug discovery and therapies with reduced side effects.

Key Experimental Findings

Recent research has demonstrated that conjugating small-molecule GPCR agonists to nanobodies (Nbs) creates bitopic ligands that span orthosteric and allosteric binding sites [90]. These constructs enable:

• High-potency activation of engineered A2A adenosine receptor (A2AR) variants
• Strong and enduring signaling responses through tethering mechanisms
• Logic-gated activity dependent on co-expression of two distinct receptor targets
• Selective targeting of receptor pairs over individual receptors

The bitopic conjugates induce signaling responses that diverge from those induced by monovalent ligands, providing new opportunities for cell type-selective signaling [90].

Table 1: Logic-Gated Signaling Outcomes with Bitopic Nb-Ligand Conjugates

Receptor Pair Targeted	Signaling Response with Single Receptor	Signaling Response with Receptor Pair	Fold Increase in cAMP	Specificity Ratio
A2AR + BC2-tagged A2AR	Minimal	Robust	8.7 ± 0.9	>20:1
A2AR + ALFA-tagged A2AR	Minimal	Robust	9.2 ± 1.1	>25:1
A2AR + PTHR1	Minimal	Moderate	5.3 ± 0.7	>15:1
Control (A2AR only)	Minimal	Minimal	1.0 ± 0.2	1:1

Table 2: Binding Affinity and Functional Potency of CGS21680-Based Conjugates

Compound	A2AR Binding Kd (nM)	cAMP EC50 (nM)	Signal Duration (t1/2, min)	Nb Epitope Specificity
CGS21680 (parent)	25 ± 3	42 ± 5	18 ± 2	N/A
CGS-NbBC2 conjugate	18 ± 2	15 ± 2	75 ± 8	BC2 epitope
CGS-NbALFA conjugate	16 ± 3	12 ± 3	82 ± 7	ALFA epitope
CGS-Nb6E conjugate	22 ± 4	28 ± 4	68 ± 6	6E epitope

Experimental Protocols

Protocol: Synthesis of Bitopic Nanobody-Ligand Conjugates

Materials and Reagents

• CGS21680 pharmacophore with PEG3-azide linker
• DBCO-functionalized nanobodies (NbBC2, NbALFA, Nb6E)
• Phosphate Buffered Saline (PBS), pH 7.4
• Dimethyl sulfoxide (DMSO), anhydrous
• Size exclusion chromatography columns (Superdex 75)
• LC-MS system for reaction monitoring

Conjugation Procedure

• Prepare ligand solution: Dissolve CGS21680-PEG3-azide in anhydrous DMSO to 10 mM stock concentration.
• Prepare nanobody solution: Dialyze DBCO-functionalized nanobodies into PBS, pH 7.4, and concentrate to 50-100 μM.
• Click chemistry conjugation: Mix ligand and nanobody solutions at 1.5:1 molar ratio (ligand:nanobody) in PBS containing 5% DMSO.
• Reaction conditions: Incubate reaction mixture for 12-16 hours at 4°C with gentle rotation.
• Purification: Purify conjugates by size exclusion chromatography using Superdex 75 column equilibrated with PBS.
• Characterization: Verify conjugation success by LC-MS and determine final concentration by UV-Vis spectroscopy.
• Quality control: Assess binding functionality by ELISA or surface plasmon resonance against target epitopes.

Reaction time may be optimized based on specific nanobody stability. Monitor conversion hourly after 8 hours to determine optimal reaction duration.

Protocol: Evaluation of Logic-Gated Signaling in Cellular Systems

Materials and Reagents

• HEK293T cell line
• Engineered A2AR constructs with epitope tags (BC2, ALFA, 6E) in mammalian expression vectors
• Transfection reagent (e.g., polyethyleneimine or lipofectamine)
• cAMP BRET biosensor (e.g., CAMYEL)
• Coelenterazine h substrate
• assay buffer (HBSS with 5 mM HEPES, pH 7.4)
• Bitopic nanobody-ligand conjugates (from Protocol 3.1)
• Control ligands (CGS21680, NECA)

Cellular Assay Procedure

• Cell culture: Maintain HEK293T cells in DMEM supplemented with 10% FBS at 37°C, 5% CO2.
• Transfection: Co-transfect cells with A2AR constructs and cAMP BRET biosensor at 90% confluence.
  • For logic-gating experiments: Transfect with combinations of tagged and untagged receptors
  • Include single-receptor controls for specificity determination
• Assay preparation: 48 hours post-transfection, harvest cells and resuspend in assay buffer at 0.5×10^6 cells/mL.
• BRET measurements:
  • Distribute cell suspension in white 96-well plates (90 μL/well)
  • Add 10 μL of coelenterazine h (5 μM final concentration)
  • Incubate 10 minutes in dark
  • Add bitopic conjugates or control ligands (10 μL in serial dilutions)
  • Measure BRET signal immediately using compatible plate reader
• Data collection: Record donor emission (475 nm) and acceptor emission (535 nm) every 30 seconds for 60 minutes.
• Data analysis: Calculate BRET ratio (acceptor/donor) and normalize to basal levels.

Include control wells with unconjugated nanobody plus ligand to verify conjugate-dependent signaling. Perform each condition in technical triplicates with at least three biological replicates.

Signaling Pathway and Experimental Workflow Diagrams

Logic-Gated GPCR Activation by Bitopic Nb-Ligand Conjugates

Experimental Workflow for Bitopic Conjugate Evaluation

Research Reagent Solutions

Table 3: Essential Research Reagents for Bitopic Ligand Studies

Reagent/Category	Specific Examples	Function and Application
Nanobodies	NbBC2, NbALFA, Nb6E	Target engineered epitope tags on GPCRs; provide high-affinity allosteric binding [90]
Small Molecule Ligands	CGS21680-PEG3-azide	Orthosteric GPCR agonist; activates A2A receptor signaling cascade [90]
Chemical Linkers	DBCO-PEG4-NHS ester, Azide-PEG3-NHS	Facilitate conjugation via click chemistry; provide spatial flexibility [90]
Engineered GPCRs	A2AR-BC2, A2AR-ALFA, A2AR-6E	Receptors with extracellular epitope tags for nanobody recognition [90]
Signaling Reporters	CAMYEL BRET biosensor	Measure cAMP production as downstream signaling readout [90]
Cell Lines	HEK293T	Heterologous expression system for receptor characterization [90]

Application Notes & Protocols

Nucleotide-binding site (NBS) encoding genes represent one of the largest families of plant disease resistance (R) genes and play a crucial role in innate immunity by recognizing diverse pathogens and initiating defense responses [24] [5]. These genes are characterized by the presence of an NBS domain and are frequently accompanied by leucine-rich repeat (LRR) domains and either a coiled-coil (CC) or Toll/Interleukin-1 receptor (TIR) domain at the N-terminus, leading to classifications such as CNL (CC-NBS-LRR) and TNL (TIR-NBS-LRR) [91] [25]. In the context of a broader thesis on genetic variation in tolerant accessions, understanding the population genomics of NBS genes is paramount. Domestication and subsequent breeding have profoundly shaped the genetic diversity of these genes, often leading to selective sweeps where favorable alleles rapidly increase in frequency. Comparative population genomics provides powerful tools to identify these sweeps, contrast variation patterns between wild and cultivated populations, and uncover genes responsible for resistance in tolerant accessions [5] [92].

Core Experimental Protocols in Population Genomics

The following section details the fundamental methodologies for conducting a comparative population genomic analysis of NBS genes.

Protocol: Identification of NBS-Encoding Genes from Genome Assemblies

Objective: To comprehensively identify and classify NBS-encoding genes from sequenced genomes of wild and domesticated accessions.

• Step 1: Data Collection

  • Obtain high-quality genome assemblies for your target species, including both wild and domesticated accessions. A chromosome-level assembly is ideal for studying genomic distribution [92].

• Step 2: HMMER-based Domain Scanning

  • Use HMMER software (e.g., PfamScan.pl) to scan all predicted protein sequences against the Pfam database.
  • The core Hidden Markov Model (HMM) to use is PF00931 (NB-ARC domain), which is characteristic of NBS genes [24] [25].
  • Filtering parameters: Use a default e-value cutoff (e.g., 1.1e-50) to ensure high-confidence hits. All genes containing the NB-ARC domain are considered NBS genes for downstream analysis [24].

• Step 3: Architectural Classification

  • Further scan the identified NBS proteins for additional domains using the same HMMER/Pfam approach.
  • Key domains to identify:
    • LRR (Leucine-Rich Repeat): PF00560, PF07723, PF07725, PF12799, PF13516, PF13855, PF14580
    • TIR (Toll/Interleukin-1 Receptor): PF01582
    • CC (Coiled-Coil): Predicted using tools like COILS or Marcoil.
  • Classify genes into structural subgroups (e.g., CNL, TNL, NL, CN, TN, N) based on the presence or absence of these domains [25].

• Step 4: Chromosomal Mapping and Synteny Analysis

  • Map the physical positions of all identified NBS genes onto the chromosomes.
  • Perform synteny analysis between the genomes of different accessions or related species to identify orthologous genomic blocks and visualize the conservation or rearrangement of NBS gene loci [25].

The workflow below illustrates the bioinformatic pipeline for identifying and classifying NBS genes.

Protocol: Population Genomic Analysis for Detecting Selective Sweeps

Objective: To identify genomic regions, particularly around NBS genes, that have undergone selection during domestication.

• Step 1: Population Sequencing and SNP Calling

  • Resequencing: Conduct whole-genome or exome sequencing of a panel of accessions, including wild populations and domesticated varieties/cultivars. A sample size of >50 individuals per group is recommended for robust statistics [93] [94].
  • Variant Calling: Map sequencing reads to a reference genome using BWA-MEM [93]. Call SNPs using GATK following its best practices pipeline [93] [45]. Apply strict filtering (e.g., QD < 2.0, FS > 60.0, MQ < 40.0, MAF > 0.01) to obtain a high-quality SNP set [93].

• Step 2: Population Structure Analysis

  • Principal Component Analysis (PCA): Perform PCA using software like PLINK to visualize genetic clustering and confirm the differentiation between wild and domesticated groups [93] [45].
  • Population Tree: Construct a Neighbour-joining tree with tools like MEGA to infer phylogenetic relationships [93].
  • Admixture Analysis: Use ADMIXTURE to estimate individual ancestries and identify potential gene flow between populations [93].

• Step 3: Scan for Selective Sweeps

  • Calculate the following statistics in sliding windows (e.g., 10-kb windows with 5-kb steps) across the genome:
    • Population Differentiation (FST): High FST values between wild and domesticated populations indicate regions under divergent selection. Calculate using VCFtools [93] [5].
    • Nucleotide Diversity (Ï€): A reduction in diversity (Ï€) in domesticated groups compared to wild groups (Ï€ domest./Ï€ wild) in a specific region is a classic signature of a selective sweep [93] [5].
    • Extended Haplotype Homozygosity (XP-EHH): This statistic identifies haplotypes that have reached high frequency in one population (domesticated) relative to another (wild), indicating recent positive selection. Calculate using Selscan [93].
  • Define Significantly Selected Regions: Merge significant outlier windows (e.g., top 1% or 5% of values) that are within 200 kb of each other into selective sweep regions [93].

• Step 4: Intersection with NBS Genes and Convergence Analysis

  • Overlap the identified selective sweep regions with the catalog of NBS genes to find candidate NBS genes under selection.
  • For comparative studies across species (e.g., wheat and barley), identify "orthoSweeps"—selective sweeps affecting orthologous genes—to find evidence of convergent selection [94].

The following diagram summarizes the key steps in the selective sweep analysis pipeline.

Protocol: Functional Validation of Candidate NBS Genes

Objective: To experimentally confirm the role of candidate NBS genes identified from genomic analyses in disease resistance.

• Step 1: Expression Analysis via RNA-seq

  • Collect RNA from tissues of interest (e.g., roots, leaves) from resistant and susceptible accessions, both under control conditions and after pathogen challenge.
  • Perform RNA sequencing and map reads to the reference genome using STAR or HISAT2.
  • Calculate expression levels (e.g., FPKM or TPM) and identify differentially expressed genes (DEGs). Candidate NBS genes showing significant upregulation in resistant accessions post-infection are strong targets for validation [45] [25].

• Step 2: Virus-Induced Gene Silencing (VIGS)

  • Cloning: Clone a ~200-300 bp fragment of the candidate NBS gene into a VIGS vector (e.g., TRV-based vector).
  • Transformation: Introduce the vector into Agrobacterium tumefaciens.
  • Infiltration: Infiltrate the Agrobacterium suspension into leaves of a resistant accession.
  • Phenotyping: After successful gene silencing, challenge the plants with the target pathogen. A loss of resistance (increased disease symptoms) in silenced plants compared to controls (e.g., empty vector) confirms the gene's functional role in immunity [24] [25].

Data Presentation and Key Findings

Table 1: Quantitative Patterns of NBS Gene Diversity in Domesticated vs. Wild Populations from Published Studies

Species / Study	Wild π (NBS genes)	Domesticated π (NBS genes)	Key Finding	Citation
Asian Pear (P. bretschneideri)	6.47E-03	6.23E-03	Decreased nucleotide diversity in cultivars, suggesting a domestication bottleneck.	[5]
European Pear (P. communis)	5.91E-03	6.48E-03	Increased nucleotide diversity in cultivars, potentially due to divergent breeding.	[5]
Tung Tree (V. montana vs V. fordii)	N/A	N/A	Resistant V. montana has 149 NBS-LRRs; susceptible V. fordii has only 90.	[25]
Lima Bean (P. lunatus)	N/A	N/A	1,917 genes with conserved disease resistance domains annotated.	[92]

Table 2: Essential Research Reagent Solutions for NBS Gene Studies

Reagent / Tool	Category	Function in Workflow	Example Use Case
HMMER (PfamScan)	Bioinformatics	Identifies protein domains (e.g., NB-ARC) from sequence data.	Initial identification of NBS-encoding genes in a genome [24] [25].
GATK	Bioinformatics	Calls and filters SNPs from resequencing data.	Generating high-quality variant sets for population analysis [93] [45].
Selscan	Bioinformatics	Calculates XP-EHH and other selection statistics.	Detecting signatures of recent positive selection [93].
OrthoFinder	Bioinformatics	Infers orthogroups across multiple species.	Identifying conserved orthologs for convergent selection analysis [24] [94].
TRV VIGS Vector	Molecular Biology	Silences target genes in plants via RNA interference.	Functional validation of candidate NBS genes in resistant accessions [24] [25].
RNA-seq Data	Data Resource	Provides gene expression profiles across tissues/conditions.	Correlating NBS gene expression with resistance phenotypes [45] [25].

The Scientist's Toolkit: Research Reagent Solutions

The table above (Table 2) details key reagents and tools. For a laboratory setting, the following materials are essential:

• Sequencing Platforms: Illumina NovaSeq for population resequencing; PacBio or Oxford Nanopore for high-quality genome assemblies.
• Computational Resources: High-performance computing (HPC) cluster with sufficient RAM (>64 GB) and storage for large-scale genomic analyses.
• Plant Materials: Well-characterized accessions, including tolerant/resistant and susceptible lines, as well as wild relatives.
• Pathogen Strains: Authenticated strains of the target pathogen for functional inoculation assays.
• Cloning Kits: Gateway or restriction enzyme-based kits for constructing VIGS vectors.

Concluding Remarks

Integrating comparative population genomics with functional validation provides a powerful framework for unraveling the evolution and mechanism of disease resistance in plants. The protocols outlined here—from in silico gene identification and sweep scanning to experimental VIGS validation—enable researchers to pinpoint key NBS genes that have been targeted during domestication and are responsible for tolerance in specific accessions. These candidate genes serve as direct targets for marker-assisted breeding to improve crop resilience.

Integrating Genomic and Transcriptomic Data to Propose a Resistance Model

Application Note

This application note details a comprehensive bioinformatic workflow for integrating genomic and transcriptomic data to model disease resistance in plants, with a specific focus on the Nucleotide-Binding Site Leucine-Rich Repeat (NBS-LRR) gene family. This protocol is designed for researchers investigating the genetic basis of disease tolerance in plant accessions, providing a robust framework for identifying key resistance genes and understanding their regulatory mechanisms.

Plant resistance to pathogens is a complex trait often governed by a sophisticated immune system. A key component of this system is the NBS-LRR gene family, which constitutes the largest class of plant disease resistance (R) genes [2] [16]. These genes encode proteins that recognize pathogen effectors and trigger a robust immune response, known as Effector-Triggered Immunity (ETI) [16]. The development of high-throughput sequencing technologies has enabled the genome-wide identification of these genes and the analysis of their expression dynamics in response to pathogen challenge.

Integrating genomic data (which reveals the potential genetic repertoire) with transcriptomic data (which reveals active gene expression under stress) provides a powerful approach to move from mere correlation to causation in resistance modeling. This integrated strategy allows for the pinpointing of specific NBS-LRR genes that are not only present in tolerant accessions but are also functionally deployed during infection. This protocol outlines a standardized pipeline for this integration, facilitating the discovery of candidate resistance genes for functional validation and breeding applications [2] [95] [16].

Materials

Research Reagent Solutions

Table 1: Essential research reagents, software tools, and databases for integrated genomic and transcriptomic analysis of resistance genes.

Item Name	Function/Application	Specifications/Examples
HMMER Suite	Identification of NBS-LRR genes using hidden Markov models.	Uses PFAM model PF00931 (NB-ARC domain) for initial search [2].
NCBI CDD	Validation of protein domains (e.g., CC, TIR, LRR).	Confirms domain completeness after HMM search [2].
Muscle	Multiple sequence alignment of NBS-LRR protein sequences.	Used for phylogenetic analysis [2].
MCScanX	Analysis of gene duplication events (tandem and segmental).	Identifies evolutionary mechanisms behind NBS-LRR family expansion [2] [16].
KaKs_Calculator	Calculation of non-synonymous (Ka) and synonymous (Ks) substitution rates.	Measures evolutionary selection pressure on NBS-LRR genes [2].
Salmon/DESeq2	Transcript quantification and differential expression analysis from RNA-seq data.	Identifies NBS-LRR genes significantly upregulated upon pathogen infection [95].
RegTools	Identification of splice-associated variants from integrated genomic and transcriptomic data.	Detects genetic variants that may cause aberrant splicing of resistance genes [96].
eQTL Analysis	Mapping genomic loci that control transcript expression levels.	Links genetic variation to expression differences in NBS-LRR genes, explaining resistance phenotypes [97].

Methods

Genome-Wide Identification and Classification of NBS-LRR Genes

• Data Acquisition: Obtain the genome assembly and annotated protein sequences for your plant species of interest from public databases (e.g., Phytozome, EnsemblPlants) [16].

• HMMER Search: Perform a hidden Markov model (HMM) search against the proteome using HMMER v3.1b2 or later and the PF00931 (NB-ARC) model from the PFAM database. Use a significance threshold (e-value) of 0.01 or lower.

• Domain Validation: Subject the initial candidate proteins to domain analysis using the NCBI Conserved Domain Database (CDD) and PFAM to confirm the presence of associated domains (e.g., TIR: PF01582, PF07725; LRR: PF00560, PF07723, PF12799; CC: confirmed via CDD). Retain only genes containing the NB-ARC domain in combination with other relevant R-gene domains [2].

• Gene Classification: Classify the identified NBS-LRR genes into subfamilies based on their N-terminal and C-terminal domain architecture:
  • CNL: CC-NBS-LRR
  • TNL: TIR-NBS-LRR
  • NL: NBS-LRR (no clear N-terminal domain)
  • RNL: RPW8-NBS-LRR [2] [16].

Phylogenetic and Evolutionary Analysis

• Sequence Alignment: Perform multiple sequence alignment of the full-length NBS-LRR protein sequences using MUSCLE or MAFFT with default parameters [2] [16].
• Phylogenetic Tree Construction: Construct a phylogenetic tree using a tool such as IQ-TREE or MEGA11. Use the maximum likelihood method with 1000 bootstrap replicates to assess branch support.
• Selection Pressure Analysis: Calculate the non-synonymous (Ka) and synonymous (Ks) substitution rates for paired NBS-LRR genes from syntenic blocks using KaKs_Calculator 2.0 under the Nei-Gojobori (NG) model. A Ka/Ks ratio > 1 indicates positive selection, which is often observed in resistance genes [2] [16].

Transcriptomic Analysis of Disease Response

• RNA-seq Data Processing: Download RNA-seq data from public repositories (e.g., NCBI SRA) for resistant and susceptible accessions, both mock-treated and pathogen-inoculated. Convert SRA files to FASTQ format and perform quality control and adapter trimming using tools like Trimmomatic [2] [95].
• Read Alignment and Quantification: Align the cleaned reads to the reference genome using a splice-aware aligner such as HISAT2. Then, quantify transcript abundances using alignment-free methods like Salmon or by counting reads per gene with featureCounts [2] [95].
• Differential Expression Analysis: Identify Differentially Expressed Genes (DEGs) using R/Bioconductor packages such as DESeq2 or edgeR. Apply a threshold of |log2 fold change| > 1 and an adjusted p-value (Benjamini-Hochberg) ≤ 0.05 to define significance. Cross-reference DEGs with the list of identified NBS-LRR genes to find those responsive to infection [95].

Integrated Genomic-Transcriptomic Analysis

• Variant Calling and Annotation: Call genetic variants (SNPs, Indels) from whole-genome sequencing data of your accessions. Annotate the resulting VCF file using RegTools' variants annotate command, flagging variants in splicing regions (e.g., within 2 bp of an exon-intron junction) [96].
• Junction Analysis: From the RNA-seq BAM files, extract splice junctions using RegTools' junctions extract command. Annotate these junctions against a reference transcriptome to identify novel or altered splicing events [96].
• Identify Splice-Associated Variants: Use RegTools' cis-splice-effects identify module to associate genomic variants with aberrant splicing events observed in the transcriptome. This pinpoints genetic mutations that directly disrupt the splicing of key NBS-LRR genes [96].
• Expression Quantitative Trait Loci (eQTL) Mapping: For a population of genetically diverse accessions, perform eQTL analysis to identify genomic loci that control the expression levels of NBS-LRR genes. This reveals trans-acting regulatory mechanisms that may underpin resistance [97].

Results and Data Analysis

Quantitative Profiling of NBS-LRR Genes

The following table summarizes the typical output from the genome-wide identification protocol, illustrating the variation in NBS-LRR family size and composition across species or accessions.

Table 2: Example identification and classification of NBS-LRR genes in three Nicotiana species [2].

Species	NBS	TIR-NBS	CC-NBS	TIR-NBS-LRR	CC-NBS-LRR	Total
N. tomentosiformis	127	7	65	33	47	279
N. sylvestris	172	5	82	37	48	344
N. tabacum	306	9	150	64	74	603

Differential Expression of NBS-LRR Genes

Analysis of RNA-seq data from pathogen-challenged samples reveals which NBS-LRR genes are functionally engaged in the resistance response.

Table 3: Example categories of disease-responsive NBS-LRR genes identified through transcriptomic analysis, as evidenced in sugarcane and banana studies [95] [16].

Gene Category	Description	Implication for Resistance Model
Early-Upregulated NBS-LRRs	Genes showing significant upregulation within hours (e.g., 12-24 h) post-inoculation.	Potential key sensors in Effector-Triggered Immunity (ETI); candidates for broad-spectrum resistance [95].
Allele-Specific Expressors	In allopolyploids, genes where expression is biased towards one progenitor's genome (e.g., from S. spontaneum in sugarcane) [16].	Explains the differential contribution of subgenomes to resistance; guides selection in breeding.
Multi-Disease Responsive NBS-LRRs	Genes differentially expressed in response to multiple distinct pathogens.	Prime candidates for engineering durable resistance against multiple diseases [16].

Workflow and Data Integration Visualization

The following diagram illustrates the complete integrated workflow for proposing a resistance model, from data input to final model validation.

Diagram 1: Integrated workflow for resistance gene discovery and modeling.

Candidate Gene Prioritization Logic

The core of the resistance model lies in the integration step, where candidates are prioritized based on converging evidence from genomic and transcriptomic analyses. The logic for this prioritization is detailed below.

Diagram 2: Logic for prioritizing candidate resistance genes.

Conclusion

The analysis of genetic variation in NBS genes of tolerant accessions reveals a complex landscape shaped by duplication, positive selection, and functional diversification. Core methodologies from genomics and transcriptomics are essential for identifying key variants, while addressing technical challenges is critical for accurate data interpretation. Functional studies confirm the direct role of specific NBS genes in pathogen resistance. Future research should leverage long-read sequencing to resolve complex loci, employ genome editing for functional characterization, and integrate population-scale data to deploy these natural genetic variations in breeding programs and the development of novel resistance strategies, ultimately enhancing crop resilience and informing biomedical discovery.

References

• 1. Comparative analysis, diversification, and functional ... [https://pubmed.ncbi.nlm.nih.gov/38789717/]
• 2. Genome-Wide Identification of NBS-LRR Family in Three ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12192410/]
• 3. Genome-Wide Identification and Expression Analysis of NBS ... [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0107987]
• 4. Genome-wide identification and characterization of NBS ... [https://www.nature.com/articles/s41598-025-03507-5]
• 5. Contrasting genetic variation and positive selection followed ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-020-07226-1]
• 6. Evaluation of reference genes for gene expression ... [https://www.nature.com/articles/s41598-019-53544-0]
• 7. Integrating genetic and gene expression data in network ... [https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-025-06143-y]
• 8. Structure/Function Approaches to Understanding CC Domain ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6290485/]
• 9. Genome-wide analysis of NBS-LRR genes in Rosaceae ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9685399/]
• 10. Diversity in Nucleotide Binding Site–Leucine-Rich Repeat ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC187567/]
• 11. Species-specific duplications driving the recent expansion of ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-015-1291-0]
• 12. TNL genes in peach: insights into the post-LRR domain [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-016-2635-0]
• 13. Identification and Characterization of NBS Resistance ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8585750/]
• 14. Whole genome sequencing reveals transcriptional and ... [https://www.nature.com/articles/s41598-025-18334-x]
• 15. A compendium of genome-wide sequence reads from NBS ... [https://www.nature.com/articles/s41598-020-67848-z]
• 16. Genome-wide analysis of NBS-LRR genes revealed ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2023.1091567/full]
• 17. Genome-Wide Characterization and Expression Pattern ... [https://www.mdpi.com/1422-0067/26/18/9063]
• 18. Genome-wide identification and characterization of NBS ... [https://bmcplantbiol.biomedcentral.com/articles/10.1186/s12870-020-02803-8]
• 19. Phylogenetic, Structural, and Evolutionary Insights into ... [https://www.mdpi.com/1422-0067/26/5/1828]
• 20. Assembly and annotation of Solanum dulcamara and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12239606/]
• 21. HMMER web server: interactive sequence similarity ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3125773/]
• 22. HMMER: biological sequence analysis using profile HMMs [https://github.com/EddyRivasLab/hmmer]
• 23. HMM-FRAME: accurate protein domain classification for ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3115854/]
• 24. Comparative analysis, diversification, and functional ... [https://www.nature.com/articles/s41598-024-62876-5]
• 25. Functional characterization of NBS-LRR genes reveals an ... [https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-024-01836-x]
• 26. Identification of Advantaged Genes for Low-Nitrogen ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12193684/]
• 27. dissectHMMER: a HMMER-based score dissection framework ... [https://biologydirect.biomedcentral.com/articles/10.1186/s13062-015-0068-3]
• 28. The NPGS Sudan sorghum germplasm collection reveals a ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12438027/]
• 29. Genome-Wide Association Study (GWAS) to Identify Salt ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7797017/]
• 30. Integrative genome-wide association and haplotype-based ... [https://bmcplantbiol.biomedcentral.com/articles/10.1186/s12870-025-07479-6]
• 31. Genome-wide association study uncovers major genetic ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8555181/]
• 32. Pangenome analysis reveals structural variation ... [https://www.nature.com/articles/s41588-025-02170-w]
• 33. Genome-Wide Association Mapping to Identify Genetic Loci ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7047875/]
• 34. GWAS combined with linkage analysis reveals major QTLs ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2024.1462856/full]
• 35. Genome-wide association analyses of risk tolerance and risky ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6713272/]
• 36. Genetic Mapping by Sequencing More Precisely Detects ... [https://www.mdpi.com/2223-7747/10/4/705]
• 37. Genome-Wide Characterization and Expression Pattern ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12469828/]
• 38. Genome-wide identification, characterization, and ... [https://www.sciencedirect.com/science/article/abs/pii/S0378111922009374]
• 39. A Beginner's Guide to Analysis of RNA Sequencing Data [https://pmc.ncbi.nlm.nih.gov/articles/PMC6096346/]
• 40. Conserved stress-responsive genes involved in the early ... [https://www.nature.com/articles/s41598-025-01436-x]
• 41. Structural variation reshapes population gene expression ... [https://www.nature.com/articles/s41588-024-01957-7]
• 42. A comprehensive workflow for optimizing RNA-seq data analysis [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-024-10414-y]
• 43. Genetic Analysis of NBS-LRR Gene Family in Chickpea ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2017.00838/full]
• 44. Association analysis identified superior haplotypes for ... [https://www.sciencedirect.com/science/article/pii/S2667064X2500168X]
• 45. Genome-wide and pan-genomic analysis reveals rich ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11033514/]
• 46. Genetic mapping and haplotype analysis identify novel ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12256648/]
• 47. Mapping of low-nitrogen tolerance genes in Japonica rice ... [https://www.nature.com/articles/s41598-025-17400-8]
• 48. Identification of Advantaged Genes for Low-Nitrogen ... [https://www.mdpi.com/1422-0067/26/12/5749]
• 49. Functional annotation of a divergent genome using sequence ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-023-09924-y]
• 50. Ka/Ks ratio [https://en.wikipedia.org/wiki/Ka/Ks_ratio]
• 51. Patterns of Positive Selection in the Complete NBS-LRR Gene ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC186657/]
• 52. Genome-wide identification of GhEDS1 gene family members ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12481744/]
• 53. Functional characterization of NBS-LRR genes reveals an ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10898138/]
• 54. Organization and molecular evolution of a disease-resistance ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-12-240]
• 55. Next-generation sequencing of newborn screening genes : the... [https://www.nature.com/articles/s41525-020-00142-z?error=cookies_not_supported&code=4cf7f599-c3da-499c-8e86-9b914ef151ef]
• 56. Navigating highly homologous genes in a molecular ... [https://www.sciencedirect.com/science/article/pii/S1098360021014258]
• 57. Overcoming the challenges in PMS2 high-homology ... [https://www.illumina.com/science/genomics-research/articles/PMS2-small-variant-detection.html]
• 58. The evolution, variation and expression patterns of ... [https://www.nature.com/articles/s41598-025-89119-5]
• 59. Analytical Validation of a Genomic Newborn Screening ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12551001/]
• 60. Systematic analysis of paralogous regions in 41755 ... [https://www.nature.com/articles/s41467-023-42531-9]
• 61. CFTR haplotype phasing using long-read genome ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11869909/]
• 62. Paralog Explorer: A resource for mining information about ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9712503/]
• 63. Paralog Explorer [https://www.flyrnai.org/tools/paralogs/web/about]
• 64. Best practices for variant calling in clinical sequencing - PMC [https://pmc.ncbi.nlm.nih.gov/articles/PMC7586657/]
• 65. Best practices for variant calling in clinical sequencing [https://genomemedicine.biomedcentral.com/articles/10.1186/s13073-020-00791-w]
• 66. Mapping‐by‐sequencing in complex polyploid genomes ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5026171/]
• 67. DaapNLRSeek, prediction, and evolution of resistance ... [https://www.sciencedirect.com/science/article/pii/S2589004225011307]
• 68. The Brassica oleracea genome reveals the asymmetrical ... [https://www.nature.com/articles/ncomms4930]
• 69. Complex genetic variation in nearly complete human ... [https://www.nature.com/articles/s41586-025-09140-6]
• 70. Genome-Wide Analysis of Genetic Variations and the ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6182389/]
• 71. Measuring, visualizing, and diagnosing reference bias with ... [https://genomebiology.biomedcentral.com/articles/10.1186/s13059-024-03240-8]
• 72. Identification and distribution of the NBS-LRR gene family in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC4422547/]
• 73. NGS technologies for analyzing germplasm diversity in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3281264/]
• 74. Reference genome and population ... [https://www.nature.com/articles/s41467-025-64825-w]
• 75. The Diversification of Plant NBS-LRR Defense Genes Directs ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5026261/]
• 76. Evolutionary analysis of TIR- and non-TIR-NBS-LRR ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2024.1452251/full]
• 77. General recommendations [https://hgvs-nomenclature.org/recommendations/general/]
• 78. Recommendations for clinical interpretation of variants ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9295495/]
• 79. Variant Annotation Table [https://support.researchallofus.org/hc/en-us/articles/4615256690836-Variant-Annotation-Table]
• 80. Functional annotation of protein sequences [https://training.galaxyproject.org/training-material/topics/genome-annotation/tutorials/functional/tutorial.html]
• 81. [Tutorial] Genome annotation [https://informatics.fas.harvard.edu/resources/tutorials/how-to-annotate-a-genome/]
• 82. Efficient Virus-Induced Gene Silencing (VIGS) Method for ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12115190/]
• 83. The allelic mutation of NBS-LRR gene causes premature ... [https://www.sciencedirect.com/science/article/pii/S0168945225000123]
• 84. Genetic variation and genome-wide association analysis of ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2025.1621916/full]
• 85. Cotton CC-NBS-LRR Gene GbCNL130 Confers Resistance to ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8456104/]
• 86. Amplicon-Based RNA Interference Targeting V2 Gene of ... [https://link.springer.com/article/10.1007/s12033-016-9980-8]
• 87. Suppression of cotton leaf curl disease symptoms in ... [https://www.sciencedirect.com/science/article/abs/pii/S0168165617316887]
• 88. Rapid and Sensitive Detection of Tobacco Streak Virus in ... [https://microbiologyjournal.org/rapid-and-sensitive-detection-of-tobacco-streak-virus-in-cotton-using-rt-lamp-and-molecular-profiling/]
• 89. Identification and Characterization of NBS Resistance ... [https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2021.758559/full]
• 90. Development of bitopic nanobody-ligand conjugates targeting ... [https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3003424]
• 91. Genome-wide identification of NBS-encoding resistance ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC2777221/]
• 92. Comprehensive genomic resources related to ... [https://www.nature.com/articles/s41467-021-20921-1]
• 93. Genome loci with selection signatures revealed in a ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-025-12145-0]
• 94. Striking convergent selection history of wheat and barley ... [https://www.nature.com/articles/s41477-025-02128-0]
• 95. Transcriptome analysis reveals key defense genes associated ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-025-12056-0]
• 96. Integrated analysis of genomic and transcriptomic data for ... [https://www.nature.com/articles/s41467-023-37266-6]
• 97. eQTLs are key players in the integration of genomic and ... [https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-022-08690-7]