Genome-Wide Analysis of NLR Disease Resistance Genes in an Updated Reference Genome of Barley

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

Immune gene variation associated with chromosome-scale differences among individual zebrafish genomes

Immune genes have evolved to maintain exceptional diversity, offering robust defense against pathogens. We performed genomic assembly to examine immune gene variation in zebrafish. Gene pathway analysis identified immune genes as significantly enriched among genes with evidence of positive selection. A large subset of genes was absent from analysis of coding sequences due to apparent lack of reads, prompting us to examine genes overlapping zero coverage regions (ZCRs), defined as 2 kb stretches without mapped reads. Immune genes were identified as highly enriched within ZCRs, including over 60% of major histocompatibility complex (MHC) genes and NOD-like receptor (NLR) genes, mediators of direct and indirect pathogen recognition. This variation was most highly concentrated throughout one arm of chromosome 4 carrying a large cluster of NLR genes, associated with large-scale structural variation covering more than half of the chromosome. Our genomic assemblies uncovered alternative haplotypes and distinct complements of immune genes among individual zebrafish, including the MHC Class II locus on chromosome 8 and the NLR gene cluster on chromosome 4. While previous studies have shown marked variation in NLR genes between vertebrate species, our study highlights extensive variation in NLR gene regions between individuals of the same species. Taken together, these findings provide evidence of immune gene variation on a scale previously unknown in other vertebrate species and raise questions about potential impact on immune function.

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.

Comparative-genomic analysis reveals dynamic NLR gene loss and gain across Apiaceae species

Introduction: Nucleotide-binding leucine-rich repeat (NLR) genes play a crucial role in green plants' responding to various pathogens. Genome-scale evolutionary studies of NLR genes are important for discovering and applying functional NLR genes. However, little is known about the evolution of NLR genes in the Apiaceae family including agricultural and medical plants. Methods: In this study, comparative genomic analysis was performed in four Apiaceae species to trace the dynamic evolutionary patterns of NLR genes during speciation in this family. Results: The results revealed different number of NLR genes in these four Apiaceae species, namely, Angelica sinensis (95), Coriandrum sativum (183), Apium graveolens (153) and Daucus carota (149). Phylogenetic analysis demonstrated that NLR genes in these four species were derived from 183 ancestral NLR lineages and experienced different levels of gene-loss and gain events. The contraction pattern of the ancestral NLR lineages was discovered during the evolution of D. carota, whereas a different pattern of contraction after first expansion of NLR genes was observed for A. sinensis, C. sativum and A. graveolens. Discussion: Taken together, rapid and dynamic gene content variation has shaped evolutionary history of NLR genes in Apiaceae species.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes

BACKGROUND

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state.

RESULTS

We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains.

CONCLUSIONS

The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

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.

Optimization of immune receptor-related hypersensitive cell death response assay using agrobacterium-mediated transient expression in tobacco plants

BACKGROUND

The study of the regulatory mechanisms of evolutionarily conserved Nucleotide-binding leucine-rich repeat (NLR) resistance (R) proteins in animals and plants is of increasing importance due to understanding basic immunity and the value of various crop engineering applications of NLR immune receptors. The importance of temperature is also emerging when applying NLR to crops responding to global climate change. In particular, studies of pathogen effector recognition and autoimmune activity of NLRs in plants can quickly and easily determine their function in tobacco using agro-mediated transient assay. However, there are conditions that should not be overlooked in these cell death-related assays in tobacco.

RESULTS

Environmental conditions play an important role in the immune response of plants. The system used in this study was to establish conditions for optimal hypertensive response (HR) cell death analysis by using the paired NLR RPS4/RRS1 autoimmune and AvrRps4 effector recognition system. The most suitable greenhouse temperature for growing plants was fixed at 22 °C. In this study, RPS4/RRS1-mediated autoimmune activity, RPS4 TIR domain-dependent cell death, and RPS4/RRS1-mediated HR cell death upon AvrRps4 perception significantly inhibited under conditions of 65% humidity. The HR is strongly activated when the humidity is below 10%. Besides, the leaf position of tobacco is important for HR cell death. Position #4 of the leaf from the top in 4-5 weeks old tobacco plants showed the most effective HR cell death.

CONCLUSIONS

As whole genome sequencing (WGS) or resistance gene enrichment sequencing (RenSeq) of various crops continues, different types of NLRs and their functions will be studied. At this time, if we optimize the conditions for evaluating NLR-mediated HR cell death, it will help to more accurately identify the function of NLRs. In addition, it will be possible to contribute to crop development in response to global climate change through NLR engineering.

Genome-wide Identification and Evolutionary Analysis of NBS-LRR Genes From Secale cereale

Secale cereale is an important crop in the Triticeae tribe of the Poaceae family, and it has unique agronomic characteristics and genome properties. It possesses resistance to many diseases and serves as an important resource for the breeding of other Triticeae crops. We performed a genome-wide study on S. cereale to identify the largest group of plant disease resistance genes (R genes), the nucleotide-binding site-leucine-rich repeat receptor (NBS-LRR) genes. In its genome, 582 NBS-LRR genes were identified, including one from the RNL subclass and 581 from the CNL subclass. The NBS-LRR gene number in the S. cereale genome is greater than that in barley and the diploid wheat genomes. S. cereale chromosome 4 contains the largest number of NBS-LRR genes among the seven chromosomes, which is different from the pattern in barley and the genomes B and D of wheat but similar to that in the genome A of wheat. Further synteny analysis suggests that more NBS-LRR genes on chromosome 4 have been inherited from a common ancestor by S. cereale and the wheat genome A than the wheat genomes B and D. Phylogenetic analysis revealed that at least 740 NBS-LRR lineages are present in the common ancestor of S. cereale, Hordeum vulgare and Triticum urartu. However, most of them have only been inherited by one or two species, with only 65 of them preserved in all three species. The S. cereale genome inherited 382 of these ancestral NBS-LRR lineages, but 120 of them have been lost in both H. vulgare and T. urartu. This study provides the full NBS-LRR profile of the S. cereale genome, which is a resource for S. cereale breeding and indicates that S. cereale can be an important material for the molecular breeding of other Triticeae crops.

Dynamic Diversity of NLR Genes in Triticum and Mining of Promising NLR Alleles for Disease Resistance

Bread wheat is an essential crop with the second-highest global production after maize. Currently, wheat diseases are a serious threat to wheat production. Therefore, efficient breeding for disease resistance is extremely urgent in modern wheat. Here, we identified 2012 NLR genes from hexaploid wheat, and Ks values of paired syntenic NLRs showed a significant peak at 3.1–6.3 MYA, which exactly coincided with the first hybridization event between A and B genome lineages at ~5.5 MYA. We provided a landscape of dynamic diversity of NLRs from Triticum and Aegilops and found that NLR genes have higher diversity in wild progenitors and relatives. Further, most NLRs had opposite diversity patterns between genic and 2 Kb-promoter regions, which might respectively link sub/neofunctionalization and loss of duplicated NLR genes. Additionally, we identified an alien introgression of chromosome 4A in tetraploid emmer wheat, which was similar to that in hexaploid wheat. Transcriptome data from four experiments of wheat disease resistance helped to profile the expression pattern of NLR genes and identified promising NLRs involved in broad-spectrum disease resistance. Our study provided insights into the diversity evolution of NLR genes and identified beneficial NLRs to deploy into modern wheat in future wheat disease-resistance breeding.