Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

Genome-wide characterization of nitric oxide-induced NBS-LRR genes from Arabidopsis thaliana and their association in monocots and dicots

BACKGROUND

Nitric oxide (NO) is pivotal in regulating the activity of NBS-LRR specific R genes, crucial components of the plant's immune system. It is noteworthy that previous research has not included a genome-wide analysis of NO-responsive NBS-LRR genes in plants.

RESULTS

The current study examined 29 NO-induced NBS-LRR genes from Arabidopsis thaliana, along with two monocots (rice and maize) and two dicots (soybean and tomato) using genome-wide analysis tools. These NBS-LRR genes were subjected to comprehensive characterization, including analysis of their physio-chemical properties, phylogenetic relationships, domain and motif identification, exon/intron structures, cis-elements, protein-protein interactions, prediction of S-Nitrosylation sites, and comparison of transcriptomic and qRT-PCR data. Results showed the diverse distribution of NBS-LRR genes across chromosomes, and variations in amino acid number, exons/introns, molecular weight, and theoretical isoelectric point, and they were found in various cellular locations like the plasma membrane, cytoplasm, and nucleus. These genes predominantly harbor the NB-ARC superfamily, LRR, LRR_8, and TIR domains, as also confirmed by motif analysis. Additionally, they feature species-specific PLN00113 superfamily and RX-CC_like domain in dicots and monocots, respectively, both responsive to defense against pathogen attacks. The NO-induced NBS-LRR genes of Arabidopsis reveal the presence of cis-elements responsive to phytohormones, light, stress, and growth, suggesting a wide range of responses mediated by NO. Protein-protein interactions, coupled with the prediction of S-Nitrosylation sites, offer valuable insights into the regulatory role of NO at the protein level within each respective species.

CONCLUSION

These above findings aimed to provide a thorough understanding of the impact of NO on NBS-LRR genes and their relationships with key plant species.

Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.)

KEY MESSAGE

Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone)

MAIN CONCLUSION

The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Comparative genomics and bioinformatics approaches revealed the role of CC-NBS-LRR genes under multiple stresses in passion fruit

Passion fruit is widely cultivated in tropical, subtropical regions of the world. The attack of bacterial and fungal diseases, and environmental factors heavily affect the yield and productivity of the passion fruit. The CC-NBS-LRR (CNL) gene family being a subclass of R-genes protects the plant against the attack of pathogens and plays a major role in effector-triggered immunity (ETI). However, no information is available regarding this gene family in passion fruit. To address the underlying problem a total of 25 and 21 CNL genes have been identified in the genome of purple (Passiflora edulis Sims.) and yellow (Passiflora edulis f. flavicarpa) passion fruit respectively. Phylogenetic tree was divided into four groups with PeCNLs present in 3 groups only. Gene structure analysis revealed that number of exons ranged from 1 to 9 with 1 being most common. Most of the PeCNL genes were clustered at the chromosome 3 and underwent strong purifying selection, expanded through segmental (17 gene pairs) and tandem duplications (17 gene pairs). PeCNL genes contained cis-elements involved in plant growth, hormones, and stress response. Transcriptome data indicated that PeCNL3, PeCNL13, and PeCNL14 were found to be differentially expressed under Cucumber mosaic virus and cold stress. Three genes were validated to be multi-stress responsive by applying Random Forest model of machine learning. To comprehend the biological functions of PeCNL proteins, their 3D structure and gene ontology (GO) enrichment analysis were done. Our research analyzed the CNL gene family in passion fruit to understand stress regulation and improve resilience. This study lays the groundwork for future investigations aimed at enhancing the genetic composition of passion fruit to ensure robust growth and productivity in challenging environments.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection

BACKGROUND

The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.).

RESULTS

The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains.

CONCLUSION

The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Comparative genomics profiling revealed multi-stress responsive roles of the CC-NBS-LRR genes in three mango cultivars

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of disease resistance (R) genes in plants and is active in response to viruses, bacteria, and fungi usually involved in effector-triggered immunity (ETI). Pangenome-wide studies allow researchers to analyze the genetic diversity of multiple species or their members simultaneously, providing a comprehensive understanding of the evolutionary relationships and diversity present among them. The draft pan-genome of three Mangifera indica cultivars (Alphonso, Hong Xiang Ya, and Tommy atkins) was constructed and Presence/absence variants (PAVs) were filtered through the ppsPCP pipeline. As a result, 2823 genes and 5907 PAVs from H. Xiang Ya, and 1266 genes and 2098 PAVs from T. atkins were added to the reference genome. For the identification of CC-NBS-LRR (CNL) genes in these mango cultivars, this draft pan-genome study has successfully identified 47, 27, and 36 members in Alphonso, H. Xiang Ya, and T. atkins respectively. The phylogenetic analysis divided MiCNL proteins into four distinct subgroups. All MiCNL genes are unevenly distributed on chromosomes. Both tandem and segmental duplication events played a significant role in the expansion of the CNL gene family. These genes contain cis-elements related to light, stress, hormone, and development. The analysis of protein-protein interactions (PPI) revealed that MiCNL proteins interacted with other defense-responsive proteins. Gene Ontology (GO) analysis indicated that MiCNL genes play a role in defense mechanisms within the organism. The expression level of the identified genes in fruit peel was observed under disease and cold stress which showed that Mi_A_CNL13 and 14 were up-regulated while Mi_A_CNL15, 25, 30, 31, and 40 were down-regulated in disease stress. On the other hand, Mi_A_CNL2, 14, 41, and 45 were up-regulated and Mi_A_CNL47 is down-regulated in cold stress. Subsequently, the Random Forest (RF) classifier was used to assess the multi-stress response of MiCNLs. It was found that Mi_A_CNL14 is a gene that responds to multiple stress conditions. The CNLs have similar protein structures which show that they are involved in the same function. The above findings provide a foundation for a deeper understanding of the functional characteristics of the mango CNL gene family.

Identification, characterization, and validation of NBS-encoding genes in grass pea

Grass pea is a promising crop with the potential to provide food and fodder, but its genomics has not been adequately explored. Identifying genes for desirable traits, such as drought tolerance and disease resistance, is critical for improving the plant. Grass pea currently lacks known R-genes, including the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, which plays a key role in protecting the plant from biotic and abiotic stresses. In our study, we used the recently published grass pea genome and available transcriptomic data to identify 274 NBS-LRR genes. The evolutionary relationships between the classified genes on the reported plants and LsNBS revealed that 124 genes have TNL domains, while 150 genes have CNL domains. All genes contained exons, ranging from 1 to 7. Ten conserved motifs with lengths ranging from 16 to 30 amino acids were identified. We found TIR-domain-containing genes in 132 LsNBSs, with 63 TIR-1 and 69 TIR-2, and RX-CCLike in 84 LsNBSs. We also identified several popular motifs, including P-loop, Uup, kinase-GTPase, ABC, ChvD, CDC6, Rnase_H, Smc, CDC48, and SpoVK. According to the gene enrichment analysis, the identified genes undergo several biological processes such as plant defense, innate immunity, hydrolase activity, and DNA binding. In the upstream regions, 103 transcription factors were identified that govern the transcription of nearby genes affecting the plant excretion of salicylic acid, methyl jasmonate, ethylene, and abscisic acid. According to RNA-Seq expression analysis, 85% of the encoded genes have high expression levels. Nine LsNBS genes were selected for qPCR under salt stress conditions. The majority of the genes showed upregulation at 50 and 200 μM NaCl. However, LsNBS-D18, LsNBS-D204, and LsNBS-D180 showed reduced or drastic downregulation compared to their respective expression levels, providing further insights into the potential functions of LsNBSs under salt stress conditions. They provide valuable insights into the potential functions of LsNBSs under salt stress conditions. Our findings also shed light on the evolution and classification of NBS-LRR genes in legumes, highlighting the potential of grass pea. Further research could focus on the functional analysis of these genes, and their potential use in breeding programs to improve the salinity, drought, and disease resistance of this important crop.

Genome-wide characterization, phylogenetic and expression analysis of ABCG gene subfamily in cucumber (Cucumis sativus L.)

The ABCG is the largest subfamily of the ABC family with extensive functions, and only a few members have been identified in detail. However, more and more studies have shown that the members of this family are very important and are involved in many life processes such as plant development and response to various stresses. Cucumber is an important vegetable crops around the world. The cucumber development is essential for its production and quality. Meanwhile, various stresses have caused serious losses of cucumber. However, the ABCG genes were not well characterized and functioned in cucumber. In this study, the cucumber CsABCG gene family were identified and characterized, and their evolutionary relationship and functions were analyzed. The cis-acting elements and expression analysis showed that they played important role in development and responding to various biotic and abiotic stresses in cucumber. Phylogenetic analysis, sequence alignment and MEME (Multiple Em for Motif Elicitation) analysis indicated that the functions of ABCG proteins in different plants are evolutionarily conserved. Collinear analysis revealed that the ABCG gene family was highly conserved during the evolution. In addition, the potential binding sites of the CsABCG genes targeted by miRNA were predicted. These results will lay a foundation for further research on the function of the CsABCG genes in cucumber.

Genome-wide identification and molecular characterization of CRK gene family in cucumber (Cucumis sativus L.) under cold stress and sclerotium rolfsii infection

BACKGROUND

The plant cysteine-rich receptor-like kinases (CRKs) are a large family having multiple roles, including defense responses under both biotic and abiotic stress. However, the CRK family in cucumbers (Cucumis sativus L.) has been explored to a limited extent. In this study, a genome-wide characterization of the CRK family has been performed to investigate the structural and functional attributes of the cucumber CRKs under cold and fungal pathogen stress.

RESULTS

A total of 15 C. sativus CRKs (CsCRKs) have been characterized in the cucumber genome. Chromosome mapping of the CsCRKs revealed that 15 genes are distributed in cucumber chromosomes. Additionally, the gene duplication analysis of the CsCRKs yielded information on their divergence and expansion in cucumbers. Phylogenetic analysis divided the CsCRKs into two clades along with other plant CRKs. Functional predictions of the CsCRKs suggested their role in signaling and defense response in cucumbers. The expression analysis of the CsCRKs by using transcriptome data and via qRT-PCR indicated their involvement in both biotic and abiotic stress responses. Under the cucumber neck rot pathogen, Sclerotium rolfsii infection, multiple CsCRKs exhibited induced expressions at early, late, and both stages. Finally, the protein interaction network prediction results identified some key possible interacting partners of the CsCRKs in regulating cucumber physiological processes.

CONCLUSIONS

The results of this study identified and characterized the CRK gene family in cucumbers. Functional predictions and validation via expression analysis confirmed the involvement of the CsCRKs in cucumber defense response, especially against S. rolfsii. Moreover, current findings provide better insights into the cucumber CRKs and their involvement in defense responses.

Genome-Wide Identification and Expression Analysis of the TIR-NBS-LRR Gene Family and Its Response to Fungal Disease in Rose (Rosa chinensis)

Roses, which are one of the world's most important ornamental plants, are often damaged by pathogens, resulting in serious economic losses. As a subclass of the disease resistance gene family of plant nucleotide-binding oligomerization domain (NOD)-like receptors, TIR-NBS-LRR (TNL) genes play a vital role in identifying pathogen effectors and activating defense responses. However, a systematic analysis of the TNL gene family is rarely reported in roses. Herein, 96 intact TNL genes were identified in Rosa chinensis. Their phylogenies, physicochemical characteristics, gene structures, conserved domains and motifs, promoter cis-elements, microRNA binding sites, and intra- and interspecific collinearity relationships were analyzed. An expression analysis using transcriptome data revealed that RcTNL genes were dominantly expressed in leaves. Some RcTNL genes responded to gibberellin, jasmonic acid, salicylic acid, Botrytis cinerea, Podosphaera pannosa, and Marssonina rosae (M. rosae); the RcTNL23 gene responded significantly to three hormones and three pathogens, and exhibited an upregulated expression. Furthermore, the black spot pathogen was identified as M. rosae. After inoculating rose leaves, an expression pattern analysis of the RcTNL genes suggested that they act during different periods of pathogen infection. The present study lays the foundations for an in-depth investigation of the TNL gene function and the mining of disease resistance genes in roses.

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

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

PlantNLRatlas: a comprehensive dataset of full- and partial-length NLR resistance genes across 100 chromosome-level plant genomes

Plants have evolved two layers of protection against biotic stress: PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). The primary mechanism of ETI involves nucleotide-binding leucine-rich repeat immune receptors (NLRs). Although NLR genes have been studied in several plant species, a comprehensive database of NLRs across a diverse array of species is still lacking. Here, we present a thorough analysis of NLR genes across 100 high-quality plant genomes (PlantNLRatlas). The PlantNLRatlas includes a total of 68,452 NLRs, of which 3,689 are full-length and 64,763 are partial-length NLRs. The majority of NLR groups were phyletically clustered. In addition, the domain sequences were found to be highly conserved within each NLR group. Our PlantNLRatlas dataset is complementary to RefPlantNLR, a collection of NLR genes which have been experimentally confirmed. The PlantNLRatlas should prove helpful for comparative investigations of NLRs across a range of plant groups, including understudied taxa. Finally, the PlantNLRatlas resource is intended to help the field move past a monolithic understanding of NLR structure and function.