Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

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.

An angiosperm NLR atlas reveals that NLR gene reduction is associated with ecological specialization and signal transduction component deletion

Nucleotide-binding site-leucine-rich repeat receptor (NLR) genes comprise the largest family of plant disease resistance genes. NLR genes are phylogenetically divided into the TNL, CNL, and RNL subclasses. NLR copy numbers and subclass composition vary tremendously across angiosperm genomes. However, the evolutionary associations between genomic NLR content and plant lifestyle, or between NLR content and signal transduction components, are poorly characterized due to limited genome availability. Here, we established an angiosperm NLRatlas (ANNA, http://compbio.nju.edu.cn/app/ANNA/), which includes NLR genes from over 300 angiosperm genomes. Using ANNA, we revealed that NLR copy numbers differ up to 66-fold among closely related species due to rapid gene loss and gain. Interestingly, NLR contraction was associated with adaptations to aquatic, parasitic, and carnivorous lifestyles. The convergent NLR reduction in aquatic plants resembles the long-term evolutionary silence of NLR genes in green algae before the colonization of land. A co-evolutionary pattern between NLR subclasses and plant immune-pathway components was also identified, suggesting that immune pathway deficiencies may drive TNL loss. Finally, we recovered a conserved TNL lineage that may function independently of the RNL pathway. Our findings provide new insights into the evolution of NLR genes in the context of plant lifestyles and genome content variation.

The genome of Chinese flowering cherry (Cerasus serrulata) provides new insights into Cerasus species

Cerasus serrulata is a flowering cherry germplasm resource for ornamental purposes. In this work, we present a de novo chromosome-scale genome assembly of C. serrulata by the use of Nanopore and Hi-C sequencing technologies. The assembled C. serrulata genome is 265.40 Mb across 304 contigs and 67 scaffolds, with a contig N50 of 1.56 Mb and a scaffold N50 of 31.12 Mb. It contains 29,094 coding genes, 27,611 (94.90%) of which are annotated in at least one functional database. Synteny analysis indicated that C. serrulata and C. avium have 333 syntenic blocks composed of 14,072 genes. Blocks on chromosome 01 of C. serrulata are distributed on all chromosomes of C. avium, implying that chromosome 01 is the most ancient or active of the chromosomes. The comparative genomic analysis confirmed that C. serrulata has 740 expanded gene families, 1031 contracted gene families, and 228 rapidly evolving gene families. By the use of 656 single-copy orthologs, a phylogenetic tree composed of 10 species was constructed. The present C. serrulata species diverged from Prunus yedoensis ~17.34 million years ago (Mya), while the divergence of C. serrulata and C. avium was estimated to have occurred ∼21.44 Mya. In addition, a total of 148 MADS-box family gene members were identified in C. serrulata, accompanying the loss of the AGL32 subfamily and the expansion of the SVP subfamily. The MYB and WRKY gene families comprising 372 and 66 genes could be divided into seven and eight subfamilies in C. serrulata, respectively, based on clustering analysis. Nine hundred forty-one plant disease-resistance genes (R-genes) were detected by searching C. serrulata within the PRGdb. This research provides high-quality genomic information about C. serrulata as well as insights into the evolutionary history of Cerasus species.

Interhomolog polymorphism shapes meiotic crossover within RAC1 and RPP13 disease resistance genes

During meiosis chromosomes undergo DNA double-strand breaks (DSBs), which can produce crossovers via interhomolog repair. Meiotic recombination frequency is variable along chromosomes and concentrates in narrow hotspots. We mapped crossovers within Arabidopsis thaliana hotspots located within the RAC1 and RPP13 disease resistance genes, using varying haplotypic combinations. We observed a negative non-linear relationship between interhomolog divergence and crossover frequency, consistent with polymorphism suppressing crossover repair of DSBs. Anti-recombinase mutants fancm, recq4a recq4b, figl1 and msh2, or lines with increased HEI10 dosage, are known to show increased crossovers. Surprisingly, RAC1 crossovers were either unchanged or decreased in these genetic backgrounds. We employed deep-sequencing of crossovers to examine recombination topology within RAC1, in wild type, fancm and recq4a recq4b mutant backgrounds. The RAC1 recombination landscape was broadly conserved in anti-recombinase mutants and showed a negative relationship with interhomolog divergence. However, crossovers at the RAC1 5’-end were relatively suppressed in recq4a recq4b backgrounds, indicating that local context influences recombination outcomes. Our results demonstrate the importance of interhomolog divergence in shaping recombination within plant disease resistance genes and crossover hotspots.

Genome-wide identification and comparative analysis of NBS-LRR resistance genes in Brassica napus

Plant disease-resistance genes play a critical role in providing resistance against pathogens. The largest family of resistance genes are the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) genes. They are classified into two major subfamilies, toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) and coiled-coil (CC)-NBS-LRR (CNL) proteins. We have identified and characterised 641 NBS-LRR genes in Brassica napus, 249 in B. rapa and 443 in B. oleracea. A ratio of 1 : 2 of CNL : TNL genes was found in the three species. Domain structure analysis revealed that 57% of the NBS-LRR genes are typical resistance genes and contain all three domains (TIR/CC, NBS, LRR), whereas the remaining genes are partially deleted or truncated. Of the NBS-LRR genes, 59% were found to be physically clustered, and individual genes involved in clusters were more polymorphic than those not clustered. Of the NBS-LRR genes in B. napus, 50% were identified as duplicates, reflecting a high level of genomic duplication and rearrangement. Comparative analysis between B. napus and its progenitor species indicated that >60% of NBS-LRR genes are conserved in B. napus. This study provides a valuable resource for the identification and characterisation of candidate NBS-LRR genes.