Simultaneous improvement of grain yield and grain protein concentration in durum wheat by using association tests and weighted GBLUP.

KEY MESSAGE: Simultaneous improvement for GY and GPC by using GWAS and GBLUP suggested a significant application in durum wheat breeding. Despite the importance of grain protein concentration (GPC) in determining wheat quality, its negative correlation with grain yield (GY) is still one of the major challenges for breeders. Here, a durum wheat panel of 200 genotypes was evaluated for GY, GPC, and their derived indices (GPD and GYD), under eight different agronomic conditions. The plant material was genotyped with the Illumina 25 k iSelect array, and a genome-wide association study was performed. Two statistical models revealed dozens of marker-trait associations (MTAs), each explaining up to 30%. phenotypic variance. Two markers on chromosomes 2A and 6B were consistently identified by both models and were found to be significantly associated with GY and GPC. MTAs identified for phenological traits co-mapped to well-known genes (i.e., Ppd-1, Vrn-1). The significance values (p-values) that measure the strength of the association of each single nucleotide polymorphism marker with the target traits were used to perform genomic prediction by using a weighted genomic best linear unbiased prediction model. The trained models were ultimately used to predict the agronomic performances of an independent durum wheat panel, confirming the utility of genomic prediction, although environmental conditions and genetic backgrounds may still be a challenge to overcome. The results generated through our study confirmed the utility of GPD and GYD to mitigate the inverse GY and GPC relationship in wheat, provided novel markers for marker-assisted selection and opened new ways to develop cultivars through genomic prediction approaches.

Phenolic Compounds and Capsaicinoids in Three Capsicum annuum Varieties: From Analytical Characterization to In Silico Hypotheses on Biological Activity.

The affinity of specific phenolic compounds (PCs) and capsaicinoids (CAPs) present in three Capsicum annuum varieties (Friariello, Cayenne and Dzuljunska Sipka) to the transient receptor potential vanilloid member 1 (TRPV1) was investigated by integrating an analytic approach for the simultaneous extraction and analysis through high-performance liquid chromatography coupled with ion trap mass spectrometry (HPLC/ITMS) and UV detection (HPLC-UV) of PCs and CAPs and structural bioinformatics based on the protein modelling and molecular simulations of protein-ligand docking. Overall, a total of 35 compounds were identified in the different samples and CAPs were quantified. The highest content of total polyphenols was recorded in the pungent Dzuljunska Sipka variety (8.91 ± 0.05 gGAE/Kg DW) while the lowest was found in the non-pungent variety Friariello (3.58 ± 0.02 gGAE/Kg DW). Protein modelling generated for the first time a complete model of the homotetrameric human TRPV1, and it was used for docking simulations with the compounds detected via the analytic approach, as well as with other compounds, as an inhibitor reference. The simulations indicate that different capsaicinoids can interact with the receptor, providing details on the molecular interaction, with similar predicted binding energy values. These results offer new insights into the interaction of capsaicinoids with TRPV1 and their possible actions.

Integration of high-throughput phenotyping with anatomical traits of leaves to help understanding lettuce acclimation to a changing environment.

MAIN CONCLUSION: The combination of image-based phenotyping with in-depth anatomical analysis allows for a thorough investigation of plant physiological plasticity in acclimation, which is driven by environmental conditions and mediated by anatomical traits. Understanding the ability of plants to respond to fluctuations in environmental conditions is critical to addressing climate change and unlocking the agricultural potential of crops both indoor and in the field. Recent studies have revealed that the degree of eco-physiological acclimation depends on leaf anatomical traits, which show stress-induced alterations during organogenesis. Indeed, it is still a matter of debate whether plant anatomy is the bottleneck for optimal plant physiology or vice versa. Here, we cultivated 'Salanova' lettuces in a phenotyping chamber under two different vapor pressure deficits (VPDs; low, high) and watering levels (well-watered, low-watered); then, plants underwent short-term changes in VPD. We aimed to combine high-throughput phenotyping with leaf anatomical analysis to evaluate their capability in detecting the early stress signals in lettuces and to highlight the different degrees of plants' eco-physiological acclimation to the change in VPD, as influenced by anatomical traits. The results demonstrate that well-watered plants under low VPD developed a morpho-anatomical structure in terms of mesophyll organization, stomatal and vein density, which more efficiently guided the acclimation to sudden changes in environmental conditions and which was not detected by image-based phenotyping alone. Therefore, we emphasized the need to complement high-throughput phenotyping with anatomical trait analysis to unveil crop acclimation mechanisms and predict possible physiological behaviors after sudden environmental fluctuations due to climate changes.

Complete Sequence, Multichromosomal Architecture and Transcriptome Analysis of the Solanum tuberosum Mitochondrial Genome.

Mitochondrial genomes (mitogenomes) in higher plants can induce cytoplasmic male sterility and be somehow involved in nuclear-cytoplasmic interactions affecting plant growth and agronomic performance. They are larger and more complex than in other eukaryotes, due to their recombinogenic nature. For most plants, the mitochondrial DNA (mtDNA) can be represented as a single circular chromosome, the so-called master molecule, which includes repeated sequences that recombine frequently, generating sub-genomic molecules in various proportions. Based on the relevance of the potato crop worldwide, herewith we report the complete mtDNA sequence of two S. tuberosum cultivars, namely Cicero and Désirée, and a comprehensive study of its expression, based on high-coverage RNA sequencing data. We found that the potato mitogenome has a multi-partite architecture, divided in at least three independent molecules that according to our data should behave as autonomous chromosomes. Inter-cultivar variability was null, while comparative analyses with other species of the Solanaceae family allowed the investigation of the evolutionary history of their mitogenomes. The RNA-seq data revealed peculiarities in transcriptional and post-transcriptional processing of mRNAs. These included co-transcription of genes with open reading frames that are probably expressed, methylation of an rRNA at a position that should impact translation efficiency and extensive RNA editing, with a high proportion of partial editing implying frequent mis-targeting by the editing machinery.

Transcriptome and metabolome of synthetic Solanum autotetraploids reveal key genomic stress events following polyploidization.

Polyploids are generally classified as autopolyploids, derived from a single species, and allopolyploids, arising from interspecific hybridization. The former represent ideal materials with which to study the consequences of genome doubling and ascertain whether there are molecular and functional rules operating following polyploidization events. To investigate whether the effects of autopolyploidization are common to different species, or if species-specific or stochastic events are prevalent, we performed a comprehensive transcriptomic and metabolomic characterization of diploids and autotetraploids of Solanum commersonii and Solanum bulbocastanum. Autopolyploidization remodelled the transcriptome and the metabolome of both species. In S. commersonii, differentially expressed genes (DEGs) were highly enriched in pericentromeric regions. Most changes were stochastic, suggesting a strong genotypic response. However, a set of robustly regulated transcripts and metabolites was also detected, including purine bases and nucleosides, which are likely to underlie a common response to polyploidization. We hypothesize that autopolyploidization results in nucleotide pool imbalance, which in turn triggers a genomic shock responsible for the stochastic events observed. The more extensive genomic stress and the higher number of stochastic events observed in S. commersonii with respect to S. bulbocastanum could be the result of the higher nucleoside depletion observed in this species.

Drought and flooding have distinct effects on herbivore-induced responses and resistance in Solanum dulcamara.

In the field, biotic and abiotic stresses frequently co-occur. As a consequence, common molecular signalling pathways governing adaptive responses to individual stresses can interact, resulting in compromised phenotypes. How plant signalling pathways interact under combined stresses is poorly understood. To assess this, we studied the consequence of drought and soil flooding on resistance of Solanum dulcamara to Spodoptera exigua and their effects on hormonal and transcriptomic profiles. The results showed that S. exigua larvae performed less well on drought-stressed plants than on well-watered and flooded plants. Both drought and insect feeding increased abscisic acid and jasmonic acid (JA) levels, whereas flooding did not induce JA accumulation. RNA sequencing analyses corroborated this pattern: drought and herbivory induced many biological processes that were repressed by flooding. When applied in combination, drought and herbivory had an additive effect on specific processes involved in secondary metabolism and defence responses, including protease inhibitor activity. In conclusion, drought and flooding have distinct effects on herbivore-induced responses and resistance. Especially, the interaction between abscisic acid and JA signalling may be important to optimize plant responses to combined drought and insect herbivory, making drought-stressed plants more resistant to insects than well-watered and flooded plants.

Characterization of the glutathione S-transferase gene family through ESTs and expression analyses within common and pigmented cultivars of Citrus sinensis (L.) Osbeck.

BACKGROUND: Glutathione S-transferases (GSTs) represent a ubiquitous gene family encoding detoxification enzymes able to recognize reactive electrophilic xenobiotic molecules as well as compounds of endogenous origin. Anthocyanin pigments require GSTs for their transport into the vacuole since their cytoplasmic retention is toxic to the cell. Anthocyanin accumulation in Citrus sinensis (L.) Osbeck fruit flesh determines different phenotypes affecting the typical pigmentation of Sicilian blood oranges. In this paper we describe: i) the characterization of the GST gene family in C. sinensis through a systematic EST analysis; ii) the validation of the EST assembly by exploiting the genome sequences of C. sinensis and C. clementina and their genome annotations; iii) GST gene expression profiling in six tissues/organs and in two different sweet orange cultivars, Cadenera (common) and Moro (pigmented). RESULTS: We identified 61 GST transcripts, described the full- or partial-length nature of the sequences and assigned to each sequence the GST class membership exploiting a comparative approach and the classification scheme proposed for plant species. A total of 23 full-length sequences were defined. Fifty-four of the 61 transcripts were successfully aligned to the C. sinensis and C. clementina genomes. Tissue specific expression profiling demonstrated that the expression of some GST transcripts was 'tissue-affected' and cultivar specific. A comparative analysis of C. sinensis GSTs with those from other plant species was also considered. Data from the current analysis are accessible at http://biosrv.cab.unina.it/citrusGST/, with the aim to provide a reference resource for C. sinensis GSTs. CONCLUSIONS: This study aimed at the characterization of the GST gene family in C. sinensis. Based on expression patterns from two different cultivars and on sequence-comparative analyses, we also highlighted that two sequences, a Phi class GST and a Mapeg class GST, could be involved in the conjugation of anthocyanin pigments and in their transport into the vacuole, specifically in fruit flesh of the pigmented cultivar.

Mitochondrial DNA editing in potato through mitoTALEN and mitoTALECD: molecular characterization and stability of editing events.

BACKGROUND: The aim of this study was to evaluate and characterize the mutations induced by two TALE-based approaches, double-strand break (DSB) induction by the FokI nuclease (mitoTALEN) and targeted base editing by the DddA cytidine deaminase (mitoTALECD), to edit, for the first time, the mitochondrial genome of potato, a vegetatively propagated crop. The two methods were used to knock out the same mitochondrial target sequence (orf125). RESULTS: Targeted chondriome deletions of different sizes (236-1066 bp) were induced by mitoTALEN due to DSB repair through ectopic homologous recombination of short direct repeats (11-12 bp) present in the target region. Furthermore, in one case, the induced DSB and subsequent repair resulted in the amplification of an already present substoichiometric molecule showing a 4288 bp deletion spanning the target sequence. With the mitoTALECD approach, both nonsense and missense mutations could be induced by base substitution. The deletions and single nucleotide mutations were either homoplasmic or heteroplasmic. The former were stably inherited in vegetative offspring. CONCLUSIONS: Both editing approaches allowed us to obtain plants with precisely modified mitochondrial genomes at high frequency. The use of the same plant genotype and mtDNA region allowed us to compare the two methods for efficiency, accuracy, type of modifications induced and stability after vegetative propagation.

Genomic analysis of the native European Solanum species, S. dulcamara.

BACKGROUND: Solanum dulcamara (bittersweet, climbing nightshade) is one of the few species of the Solanaceae family native to Europe. As a common weed it is adapted to a wide range of ecological niches and it has long been recognized as one of the alternative hosts for pathogens and pests responsible for many important diseases in potato, such as Phytophthora. At the same time, it may represent an alternative source of resistance genes against these diseases. Despite its unique ecology and potential as a genetic resource, genomic research tools are lacking for S. dulcamara. We have taken advantage of next-generation sequencing to speed up research on and use of this non-model species. RESULTS: In this work, we present the first large-scale characterization of the S. dulcamara transcriptome. Through comparison of RNAseq reads from two different accessions, we were able to predict transcript-based SNP and SSR markers. Using the SNP markers in combination with genomic AFLP and CAPS markers, the first genome-wide genetic linkage map of bittersweet was generated. Based on gene orthology, the markers were anchored to the genome of related Solanum species (tomato, potato and eggplant), revealing both conserved and novel chromosomal rearrangements. This allowed a better estimation of the evolutionary moment of rearrangements in a number of cases and showed that chromosomal breakpoints are regularly re-used. CONCLUSION: Knowledge and tools developed as part of this study pave the way for future genomic research and exploitation of this wild Solanum species. The transcriptome assembly represents a resource for functional analysis of genes underlying interesting biological and agronomical traits and, in the absence of the full genome, provides a reference for RNAseq gene expression profiling aimed at understanding the unique biology of S. dulcamara. Cross-species orthology-based marker selection is shown to be a powerful tool to quickly generate a comparative genetic map, which may speed up gene mapping and contribute to the understanding of genome evolution within the Solanaceae family.

In silico and in vitro approaches allow the identification of the Prosystemin molecular network.

Tomato Prosystemin (ProSys), the precursor of Systemin, a small peptidic hormone, is produced at very low concentration in unchallenged plants, while its expression greatly increases in response to several different stressors triggering an array of defence responses. The molecular mechanisms that underpin such a wide array of defence barriers are not fully understood and are likely correlated with the intrinsically disordered (ID) structure of the protein. ID proteins interact with different protein partners forming complexes involved in the modulation of different biological mechanisms. Here we describe the ProSys-protein network that shed light on the molecular mechanisms underpinning ProSys associated defence responses. Three different approaches were used. In silico prediction resulted in 98 direct interactors, most clustering in phytohormone biosynthesis, transcription factors and signal transduction gene classes. The network shows the central role of ProSys during defence responses, that reflects its role as central hub. In vitro ProSys interactors, identified by Affinity Purification-Mass Spectrometry (AP-MS), revealed over three hundred protein partners, while Bimolecular Fluorescent Complementation (BiFC) experiments validated in vivo some interactors predicted in silico and in vitro. Our results demonstrate that ProSys interacts with several proteins and reveal new key molecular events in the ProSys-dependent defence response of tomato plant.

Corrigendum to "In silico and in vitro approaches allow the identification of the Prosystemin molecular network" [Comput Struct Biotechnol J 21 (2023) 212-223].

[This corrects the article DOI: 10.1016/j.csbj.2022.12.006.].

Genetic diversity and signature of divergence in the genome of grapevine clones of Southern Italy varieties.

Sexual reproduction has contributed to a significant degree of variability in cultivated grapevine populations. However, the additional influence of spontaneous somatic mutations has played a pivotal role in shaping the diverse landscape of grapevine agrobiodiversity. These naturally occurring selections, termed 'clones,' represent a vast reservoir of potentially valuable traits and alleles that hold promise for enhancing grape quality and bolstering plant resilience against environmental and biotic challenges. Despite their potential, many of these clones remain largely untapped.In light of this context, this study aims to delve into the population structure, genetic diversity, and distinctive genetic loci within a collection of 138 clones derived from six Campanian and Apulian grapevine varieties, known for their desirable attributes in viticulture and winemaking. Employing two reduced representation sequencing methods, we extracted Single-Nucleotide Polymorphism (SNP) markers. Population structure analysis and fixation index (FST) calculations were conducted both between populations and at individual loci. Notably, varieties originating from the same geographical region exhibited pronounced genetic similarity.The resulting SNP dataset facilitated the identification of approximately two hundred loci featuring divergent markers (FST ≥ 0.80) within annotated exons. Several of these loci exhibited associations with essential traits like phenotypic adaptability and environmental responsiveness, offering compelling opportunities for grapevine breeding initiatives. By shedding light on the genetic variability inherent in these treasured traditional grapevines, our study contributes to the broader understanding of their potential. Importantly, it underscores the urgency of preserving and characterizing these valuable genetic resources to safeguard their intra-varietal diversity and foster future advancements in grapevine cultivation.

Revealing impaired pathways in the an11 mutant by high-throughput characterization of Petunia axillaris and Petunia inflata transcriptomes.

Petunia is an excellent model system, especially for genetic, physiological and molecular studies. Thus far, however, genome-wide expression analysis has been applied rarely because of the lack of sequence information. We applied next-generation sequencing to generate, through de novo read assembly, a large catalogue of transcripts for Petunia axillaris and Petunia inflata. On the basis of both transcriptomes, comprehensive microarray chips for gene expression analysis were established and used for the analysis of global- and organ-specific gene expression in Petunia axillaris and Petunia inflata and to explore the molecular basis of the seed coat defects in a Petunia hybrida mutant, anthocyanin 11 (an11), lacking a WD40-repeat (WDR) transcription regulator. Among the transcripts differentially expressed in an11 seeds compared with wild type, many expected targets of AN11 were found but also several interesting new candidates that might play a role in morphogenesis of the seed coat. Our results validate the combination of next-generation sequencing with microarray analyses strategies to identify the transcriptome of two petunia species without previous knowledge of their genome, and to develop comprehensive chips as useful tools for the analysis of gene expression in P. axillaris, P. inflata and P. hybrida.

Whole-exome sequencing of selected bread wheat recombinant inbred lines as a useful resource for allele mining and bulked segregant analysis.

Although wheat (Triticum aestivum L.) is the main staple crop in the world and a major source of carbohydrates and proteins, functional genomics and allele mining are still big challenges. Given the advances in next-generation sequencing (NGS) technologies, the identification of causal variants associated with a target phenotype has become feasible. For these reasons, here, by combining sequence capture and target-enrichment methods with high-throughput NGS re-sequencing, we were able to scan at exome-wide level 46 randomly selected bread wheat individuals from a recombinant inbred line population and to identify and classify a large number of single nucleotide polymorphisms (SNPs). For technical validation of results, eight randomly selected SNPs were converted into Kompetitive Allele-Specific PCR (KASP) markers. This resource was established as an accessible and reusable molecular toolkit for allele data mining. The dataset we are making available could be exploited for novel studies on bread wheat genetics and as a foundation for starting breeding programs aimed at improving different key agronomic traits.

Sweet Cherry Diversity and Relationships in Modern and Local Varieties Based on SNP Markers.

The sweet cherry is an important fruit species that is widespread globally. In addition to the well-known traditional and modern varieties, a myriad of landraces is present in Europe, as well as in southern Italy. This study aims to evaluate the population structure, genetic relationships, and cases of duplicate samples in a collection of 143 accessions using GBS-derived SNP markers. The genetic material under investigation includes modern commercial varieties, ancient European and American varieties, landraces, and individuals retrieved from small orchards. Some of the known varieties were genetically analyzed here for the first time. In addition, several genotypes were collected from the Basilicata region (southern Italy), an area largely unexplored for sweet cherry genetic resources. The relationships among genotypes were assessed using four different methods: allele frequency and ancestry estimation, principal component analysis, Neighbor-Joining tree, and identity-by-state estimation. The analyses returned quite congruent results and highlighted the presence of four main genetic groups, namely: (i) American varieties, (ii) the 'Germersdorfer-Ferrovia' cluster, (iii) the 'Burlat' group, and (iv) the group of Italian landraces. The main drivers of clustering were ancestry, geographical distribution, and some important traits such as self-compatibility. The sweet cherries from Basilicata, herewith examined for the first time, were mostly distributed within the group of Italian landraces, being particularly linked to the autochthonous varieties of the Campania region. However, some genotypes were outside this group, thus suggesting the introduction of genetic material from other Italian regions or from European countries. The considerable amount of American and European modern varieties analyzed are genetically very closely related, suggesting a reduced genetic basis. In addition, we highlighted the discriminating ability of SNP markers to distinguish between an original variety and its mutant. Overall, our results may be useful in defining conservation strategies for sweet cherry germplasm and developing future breeding programs to enlarge the genetic basis of commercial varieties.

Phyto-Friendly Soil Bacteria and Fungi Provide Beneficial Outcomes in the Host Plant by Differently Modulating Its Responses through (In)Direct Mechanisms.

Sustainable agricultural systems based on the application of phyto-friendly bacteria and fungi are increasingly needed to preserve soil fertility and microbial biodiversity, as well as to reduce the use of chemical fertilizers and pesticides. Although there is considerable attention on the potential applications of microbial consortia as biofertilizers and biocontrol agents for crop management, knowledge on the molecular responses modulated in host plants because of these beneficial associations is still incomplete. This review provides an up-to-date overview of the different mechanisms of action triggered by plant-growth-promoting microorganisms (PGPMs) to promote host-plant growth and improve its defense system. In addition, we combined available gene-expression profiling data from tomato roots sampled in the early stages of interaction with Pseudomonas or Trichoderma strains to develop an integrated model that describes the common processes activated by both PGPMs and highlights the host's different responses to the two microorganisms. All the information gathered will help define new strategies for the selection of crop varieties with a better ability to benefit from the elicitation of microbial inoculants.

Merging genotyping-by-sequencing data from two ex situ collections provides insights on the pea evolutionary history.

Pea (Pisum sativum L. subsp. sativum) is one of the oldest domesticated species and a widely cultivated legume. In this study, we combined next generation sequencing (NGS) data referring to two genotyping-by-sequencing (GBS) libraries, each one prepared from a different Pisum germplasm collection. The selection of single nucleotide polymorphism (SNP) loci called in both germplasm collections caused some loss of information; however, this did not prevent the obtainment of one of the largest datasets ever used to explore pea biodiversity, consisting of 652 accessions and 22 127 markers. The analysis of population structure reflected genetic variation based on geographic patterns and allowed the definition of a model for the expansion of pea cultivation from the domestication centre to other regions of the world. In genetically distinct populations, the average decay of linkage disequilibrium (LD) ranged from a few bases to hundreds of kilobases, thus indicating different evolutionary histories leading to their diversification. Genome-wide scans resulted in the identification of putative selective sweeps associated with domestication and breeding, including genes known to regulate shoot branching, cotyledon colour and resistance to lodging, and the correct mapping of two Mendelian genes. In addition to providing information of major interest for fundamental and applied research on pea, our work describes the first successful example of integration of different GBS datasets generated from ex situ collections - a process of potential interest for a variety of purposes, including conservation genetics, genome-wide association studies, and breeding.

Multi-omics data integration provides insights into the post-harvest biology of a long shelf-life tomato landrace.

In this study we investigated the transcriptome and epigenome dynamics of the tomato fruit during post-harvest in a landrace belonging to a group of tomatoes (Solanum lycopersicum L.) collectively known as "Piennolo del Vesuvio", all characterized by a long shelf-life. Expression of protein-coding genes and microRNAs as well as DNA methylation patterns and histone modifications were analysed in distinct post-harvest phases. Multi-omics data integration contributed to the elucidation of the molecular mechanisms underlying processes leading to long shelf-life. We unveiled global changes in transcriptome and epigenome. DNA methylation increased and the repressive histone mark H3K27me3 was lost as the fruit progressed from red ripe to 150 days post-harvest. Thousands of genes were differentially expressed, about half of which were potentially epi-regulated as they were engaged in at least one epi-mark change in addition to being microRNA targets in ~5% of cases. Down-regulation of the ripening regulator MADS-RIN and of genes involved in ethylene response and cell wall degradation was consistent with the delayed fruit softening. Large-scale epigenome reprogramming that occurred in the fruit during post-harvest likely contributed to delayed fruit senescence.

Guidelines for Setting Up a mRNA Sequencing Experiment and Best Practices for Bioinformatic Data Analysis.

RNA-sequencing, commonly referred to as RNA-seq, is the most recently developed method for the analysis of transcriptomes. It uses high-throughput next-generation sequencing technologies and has revolutionized our understanding of the complexity and dynamics of whole transcriptomes.In this chapter, we recall the key developments in transcriptome analysis and dissect the different steps of the general workflow that can be run by users to design and perform a mRNA-seq experiment as well as to process mRNA-seq data obtained by the Illumina technology. The chapter proposes guidelines for completing a mRNA-seq study properly and makes available recommendations for best practices based on recent literature and on the latest developments in technology and algorithms. We also remark the large number of choices available (especially for bioinformatic data analysis) in front of which the scientist may be in trouble.In the last part of the chapter we discuss the new frontiers of single-cell RNA-seq and isoform sequencing by long read technology.

The Tomato Interspecific NB-LRR Gene Arsenal and Its Impact on Breeding Strategies.

Tomato (Solanum lycopersicum L.) is a model system for studying the molecular basis of resistance in plants. The investigation of evolutionary dynamics of tomato resistance (R)-loci provides unique opportunities for identifying factors that promote or constrain genome evolution. Nucleotide-binding domain and leucine-rich repeat (NB-LRR) receptors belong to one of the most plastic and diversified families. The vast amount of genomic data available for Solanaceae and wild tomato relatives provides unprecedented insights into the patterns and mechanisms of evolution of NB-LRR genes. Comparative analysis remarked a reshuffling of R-islands on chromosomes and a high degree of adaptive diversification in key R-loci induced by species-specific pathogen pressure. Unveiling NB-LRR natural variation in tomato and in other Solanaceae species offers the opportunity to effectively exploit genetic diversity in genomic-driven breeding programs with the aim of identifying and introducing new resistances in tomato cultivars. Within this motivating context, we reviewed the repertoire of NB-LRR genes available for tomato improvement with a special focus on signatures of adaptive processes. This issue is still relevant and not thoroughly investigated. We believe that the discovery of mechanisms involved in the generation of a gene with new resistance functions will bring great benefits to future breeding strategies.