Evolution of NLR Resistance Genes in Magnoliids: Dramatic Expansions of CNLs and Multiple Losses of TNLs

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Natural selection directing molecular evolution in vertebrate viral sensors

Diseases caused by pathogens contribute to molecular adaptations in host immunity. Variety of viral pathogens challenging animal immunity can drive positive selection diversifying receptors recognising the infections. However, whether distinct virus sensing systems differ across animals in their evolutionary modes remains unclear. Our review provides a comparative overview of natural selection shaping molecular evolution in vertebrate viral-binding pattern recognition receptors (PRRs). Despite prevailing negative selection arising from the functional constraints, multiple lines of evidence now suggest diversifying selection in the Toll-like receptors (TLRs), NOD-like receptors (NLRs), RIG-I-like receptors (RLRs) and oligoadenylate synthetases (OASs). In several cases, location of the positively selected sites in the ligand-binding regions suggests effects on viral detection although experimental support is lacking. Unfortunately, in most other PRR families including the AIM2-like receptor family, C-type lectin receptors (CLRs), and cyclic GMP-AMP synthetase studies characterising their molecular evolution are rare, preventing comparative insight. We indicate shared characteristics of the viral sensor evolution and highlight priorities for future research.

Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding

Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.

Evolution of NLR genes in genus Arachis reveals asymmetric expansion of NLRome in wild and domesticated tetraploid species

Arachis hypogaea is an allotetraploid crop widely grown in the world. Wild relatives of genus Arachis are the rich source of genetic diversity and high levels of resistance to combat pathogens and climate change. The accurate identification and characterization of plant resistance gene, nucleotide binding site leucine rich repeat receptor (NLRs) substantially contribute to the repertoire of resistances and improve production. In the current study, we have studied the evolution of NLR genes in genus Arachis and performed their comparative genomics among four diploids (A. duranensis, A. ipaensis, A. cardenasii, A. stenosperma) and two tetraploid (wild: A. monticola and domesticated: A. hypogaea) species. In total 521, 354, 284, 794, 654, 290 NLR genes were identified from A. cardenasii, A. stenosperma and A. duranensis, A. hypogaea, A. monticola and A. ipaensis respectively. Phylogenetic analysis and classification of NLRs revealed that they belong to 7 subgroups and specific subgroups have expanded in each genome leading towards divergent evolution. Gene gain and loss, duplication assay reveals that wild and domesticated tetraploids species have shown asymmetric expansion of NLRome in both sub-genome (AA and BB). A-subgenome of A. monticola exhibited significant contraction of NLRome while B-subgenome shows expansion and vice versa in case of A. hypogaea probably due to distinct natural and artificial selection pressure. In addition, diploid species A. cardenasii revealed the largest repertoire of NLR genes due to higher frequency of gene duplication and selection pressure. A. cardenasii and A. monticola can be regarded as putative resistance resources for peanut breeding program for introgression of novel resistance genes. Findings of this study also emphasize the application neo-diploids and polyploids due to higher quantitative expression of NLR genes. To the best of our knowledge, this is the first study that studied the effect of domestication and polyploidy on the evolution of NLR genes in genus Arachis to identify genomic resources for improving resistance of polyploid crop with global importance on economy and food security.

Investigation of Resistance Genes in Genus Vigna Reveals Highly Variable NLRome in Parallel Domesticated Member Species

Vigna is a unique genus that consist of multiple crop species that are domesticated in parallel fashion between 7-10 thousand years ago. Here we studied the evolution of nucleotide-binding site leucine-rich repeat receptor (NLR) genes across five crop species of genus Vigna. In total identified 286, 350, 234, 250, 108 and 161 NLR genes were from Phaseolous vulgaris, Vigna. unguiculata, Vigna mungo, Vigna radiata, Vigna angularis and Vigna umbellata respectively. Comprehensive phylogenetic and clusterization analysis reveals the presence of seven subgroups of Coiled coil like NLRs (CC-NLR) genes and four distinct lineages of Toll interleukin receptor like NLRs (TIR-NLR). Subgroup CC_G10-NLR shows large scale diversification among Vigna species suggesting genus specific distinct duplication pattern in Vigna species. Mainly birth of new NLR gene families and higher rate of terminal duplication is the major determinants for expansion of NLRome in genus Vigna. Recent expansion of NLRome in V. anguiculata and V. radiata was also observed which might suggest that domestication have supported their duplication of lineage specific NLR genes. In short, large scale difference in the architecture of NLRome were observed in diploid plant species. Our findings allowed us to hypothesized that independent parallel domestication is the major drivers of highly divergent evolution of NLRome in genus Vigna.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

Genome-wide characterization of the NLR gene family in tomato (Solanum lycopersicum) and their relatedness to disease resistance

Nucleotide-binding leucine-rich-repeat receptors (NLR), the largest group of genes associated with plant disease resistance (R), have attracted attention due to their crucial role in protecting plants from pathogens. Genome-wide studies of NLRs have revealed conserved domains in the annotated tomato genome. The 321 NLR genes identified in the tomato genome have been randomly mapped to 12 chromosomes. Phylogenetic analysis and classification of NLRs have revealed that 211 genes share full-length domains categorized into three major clades (CNL, TNL, and RNL); the remaining 110 NLRs share partial domains and are classified in CN, TN, and N according to their motifs and gene structures. The cis-regulatory elements of NLRs exhibit the maximum number of these elements and are involved in response to biotic and abiotic stresses, pathogen recognition, and resistance. Analysis of the phylogenetic relationship between tomato NLRs and orthologs in other species has shown conservation among Solanaceae members and variation with A. thaliana. Synteny and Ka/Ks analyses of Solanum lycopersicum and Solanum tuberosum orthologs have underscored the importance of NLR conservation and diversification from ancestral species millions of years ago. RNA-seq data and qPCR analysis of early and late blight diseases in tomatoes revealed consistent NLR expression patterns, including upregulation in infected compared to control plants (with some exceptions), suggesting the role of NLRs as key regulators in early blight resistance. Moreover, the expression levels of NLRs associated with late blight resistance (Solyc04g007060 [NRC4] and Solyc10g008240 [RIB12]) suggested that they regulate S. lycopersicum resistance to P. infestans. These findings provide important fundamental knowledge for understanding NLR evolution and diversity and will empower the broader characterization of disease resistance genes for pyramiding through speed cloning to develop disease-tolerant varieties.

NLRscape: an atlas of plant NLR proteins

NLRscape is a webserver that curates a collection of over 80 000 plant protein sequences identified in UniProtKB to contain NOD-like receptor signatures, and hosts in addition a number of tools aimed at the exploration of the complex sequence landscape of this class of plant proteins. Each entry gathers sequence information, domain and motif annotations from multiple third-party sources but also in-house advanced annotations aimed at addressing caveats of the existing broad-based annotations. NLRscape provides a top-down perspective of the NLR sequence landscape but also services for assisting a bottom-up approach starting from a given input sequence. Sequences are clustered by their domain organization layout, global homology and taxonomic spread-in order to allow analysis of how particular traits of an NLR family are scattered within the plant kingdom. Tools are provided for users to locate their own protein of interest in the overall NLR landscape, generate custom clusters centered around it and perform a large number of sequence and structural analyses using included interactive online instruments. Amongst these, we mention: taxonomy distribution plots, homology cluster graphs, identity matrices and interactive MSA synchronizing secondary structure and motif predictions. NLRscape can be found at: https://nlrscape.biochim.ro/.

Molecular innovations in plant TIR-based immunity signaling

A protein domain (Toll and Interleukin-1 receptor [TIR]-like) with homology to animal TIRs mediates immune signaling in prokaryotes and eukaryotes. Here, we present an overview of TIR evolution and the molecular versatility of TIR domains in different protein architectures for host protection against microbial attack. Plant TIR-based signaling emerges as being central to the potentiation and effectiveness of host defenses triggered by intracellular and cell-surface immune receptors. Equally relevant for plant fitness are mechanisms that limit potent TIR signaling in healthy tissues but maintain preparedness for infection. We propose that seed plants evolved a specialized protein module to selectively translate TIR enzymatic activities to defense outputs, overlaying a more general function of TIRs.