The de Novo Reference Genome and Transcriptome Assemblies of the Wild Tomato Species Solanum chilense Highlights Birth and Death of NLR Genes Between Tomato Species

High-Resolution Disease Phenotyping Reveals Distinct Resistance Mechanisms of Tomato Crop Wild Relatives against Sclerotinia sclerotiorum

Besides the well-understood qualitative disease resistance, plants possess a more complex quantitative form of resistance: quantitative disease resistance (QDR). QDR is commonly defined as a partial but more durable form of resistance and, therefore, might display a valuable target for resistance breeding. The characterization of QDR phenotypes, especially of wild crop relatives, displays a bottleneck in deciphering QDR's genomic and regulatory background. Moreover, the relationship between QDR parameters, such as infection frequency, lag-phase duration, and lesion growth rate, remains elusive. High hurdles for applying modern phenotyping technology, such as the low availability of phenotyping facilities or complex data analysis, further dampen progress in understanding QDR. Here, we applied a low-cost (<1.000 €) phenotyping system to measure lesion growth dynamics of wild tomato species (e.g., Solanum pennellii or Solanum pimpinellifolium). We provide insight into QDR diversity of wild populations and derive specific QDR mechanisms and their cross-talk. We show how temporally continuous observations are required to dissect end-point severity into functional resistance mechanisms. The results of our study show how QDR can be maintained by facilitating different defense mechanisms during host-parasite interaction and that the capacity of the QDR toolbox highly depends on the host's genetic context. We anticipate that the present findings display a valuable resource for more targeted functional characterization of the processes involved in QDR. Moreover, we show how modest phenotyping technology can be leveraged to help answer highly relevant biological questions.

Pathogen elicitor peptide (pep), systemin, and their receptors in tomato: sequence analysis sheds light on standing disagreements about biotic stress signaling components

Peps are endogenous damage-associated polypeptides that evoke defense responses in plants. Like other damage-associated molecular patterns, Pep signals are transduced by receptors. PEPRs are the receptors that transduce Pep danger signals. This paper identifies new putative Peps in the Solanaceae (including Solanum spp., Nicotiana spp., and Petunia spp.) and Coffea and explores their properties. Using these newly identified Peps we derive sequence logos that present a refinement of the current understanding of the importance of specific residues in the Pep signaling molecules in Solanaceae, including several arginines, prolines that restrict peptide's conformations, and C-terminal asparagine. We examine the degree of disorder in Pep, which is likely important to the mechanism of Pep perception. This work also calls into question some of the evolutionary relationships between Peps in Solanaceae and specific Arabidopsis Peps published in previous literature, culminating in a conclusion that SlPep should not be named SlPep6 due to the lack of conservation of protein sequences in AtPROPEP6 and SlPROPEP, and that SlPep probably does not have two receptors in tomato, based on phylogenetic analysis. Our analyses advance understanding of the Pep signaling system in Solanaceae.

An S-locus F-box protein as pollen S determinant targets non-self S-RNase underlying self-incompatibility in Citrus

Self-incompatibility (SI) is a crucial mechanism that prevents self-fertilization and inbreeding in flowering plants. Citrus exhibits SI regulated by a polymorphic S-locus containing an S-RNase gene and multiple S-locus F-box (SLF) genes. It has been documented that S-RNase functions as the pistil S determinant, but there is no direct evidence that the SLF genes closely linked with S-RNase function as pollen S determinants in Citrus. This study assembled the genomes of two pummelo (Citrus grandis) plants, obtained three novel complete and well-annotated S-haplotypes, and isolated 36 SLF or SLF-like alleles on the S-loci. Phylogenetic analysis of 138 SLFs revealed that the SLF genes were classified into 12 types, including six types with divergent or missing alleles. Furthermore, transformation experiments verified that the conserved S6-SLF7a protein can lead to the transition of SI to self-compatibility by recognizing non-self S8-RNase in 'Mini-Citrus' plants (S7S8 and S8S29, Fortunella hindsii), a model plant for citrus gene function studies. In vitro assays demonstrated interactions between SLFs of different S haplotypes and the Skp1-Cullin1-F-box subunit CgSSK1 protein. This study provides direct evidence that SLF controls the pollen function in Citrus, demonstrating its role in the 'non-self recognition' SI system.

A chromosome-level genome assembly of Solanum chilense, a tomato wild relative associated with resistance to salinity and drought

Introduction

Solanum chilense is a wild relative of tomato reported to exhibit resistance to biotic and abiotic stresses. There is potential to improve tomato cultivars via breeding with wild relatives, a process greatly accelerated by suitable genomic and genetic resources.

Methods

In this study we generated a high-quality, chromosome-level, de novo assembly for the S. chilense accession LA1972 using a hybrid assembly strategy with ~180 Gbp of Illumina short reads and ~50 Gbp long PacBio reads. Further scaffolding was performed using Bionano optical maps and 10x Chromium reads.

Results

The resulting sequences were arranged into 12 pseudomolecules using Hi-C sequencing. This resulted in a 901 Mbp assembly, with a completeness of 95%, as determined by Benchmarking with Universal Single-Copy Orthologs (BUSCO). Sequencing of RNA from multiple tissues resulting in ~219 Gbp of reads was used to annotate the genome assembly with an RNA-Seq guided gene prediction, and for a de novo transcriptome assembly. This chromosome-level, high-quality reference genome for S. chilense accession LA1972 will support future breeding efforts for more sustainable tomato production.

Discussion

Gene sequences related to drought and salt resistance were compared between S. chilense and S. lycopersicum to identify amino acid variations with high potential for functional impact. These variants were subsequently analysed in 84 resequenced tomato lines across 12 different related species to explore the variant distributions. We identified a set of 7 putative impactful amino acid variants some of which may also impact on fruit development for example the ethylene-responsive transcription factor WIN1 and ethylene-insensitive protein 2. These variants could be tested for their ability to confer functional phenotypes to cultivars that have lost these variants.

AlphaFold-Multimer predicts cross-kingdom interactions at the plant-pathogen interface

Adapted plant pathogens from various microbial kingdoms produce hundreds of unrelated small secreted proteins (SSPs) with elusive roles. Here, we used AlphaFold-Multimer (AFM) to screen 1879 SSPs of seven tomato pathogens for interacting with six defence-related hydrolases of tomato. This screen of 11,274 protein pairs identified 15 non-annotated SSPs that are predicted to obstruct the active site of chitinases and proteases with an intrinsic fold. Four SSPs were experimentally verified to be inhibitors of pathogenesis-related subtilase P69B, including extracellular protein-36 (Ecp36) and secreted-into-xylem-15 (Six15) of the fungal pathogens Cladosporium fulvum and Fusarium oxysporum, respectively. Together with a P69B inhibitor from the bacterial pathogen Xanthomonas perforans and Kazal-like inhibitors of the oomycete pathogen Phytophthora infestans, P69B emerges as an effector hub targeted by different microbial kingdoms, consistent with a diversification of P69B orthologs and paralogs. This study demonstrates the power of artificial intelligence to predict cross-kingdom interactions at the plant-pathogen interface.

Insights into the Transcriptomics of Crop Wild Relatives to Unravel the Salinity Stress Adaptive Mechanisms

The narrow genomic diversity of modern cultivars is a major bottleneck for enhancing the crop's salinity stress tolerance. The close relatives of modern cultivated plants, crop wild relatives (CWRs), can be a promising and sustainable resource to broaden the diversity of crops. Advances in transcriptomic technologies have revealed the untapped genetic diversity of CWRs that represents a practical gene pool for improving the plant's adaptability to salt stress. Thus, the present study emphasizes the transcriptomics of CWRs for salinity stress tolerance. In this review, the impacts of salt stress on the plant's physiological processes and development are overviewed, and the transcription factors (TFs) regulation of salinity stress tolerance is investigated. In addition to the molecular regulation, a brief discussion on the phytomorphological adaptation of plants under saline environments is provided. The study further highlights the availability and use of transcriptomic resources of CWR and their contribution to pangenome construction. Moreover, the utilization of CWRs' genetic resources in the molecular breeding of crops for salinity stress tolerance is explored. Several studies have shown that cytoplasmic components such as calcium and kinases, and ion transporter genes such as Salt Overly Sensitive 1 (SOS1) and High-affinity Potassium Transporters (HKTs) are involved in the signaling of salt stress, and in mediating the distribution of excess Na⁺ ions within the plant cells. Recent comparative analyses of transcriptomic profiling through RNA sequencing (RNA-Seq) between the crops and their wild relatives have unraveled several TFs, stress-responsive genes, and regulatory proteins for generating salinity stress tolerance. This review specifies that the use of CWRs transcriptomics in combination with modern breeding experimental approaches such as genomic editing, de novo domestication, and speed breeding can accelerate the CWRs utilization in the breeding programs for enhancing the crop's adaptability to saline conditions. The transcriptomic approaches optimize the crop genomes with the accumulation of favorable alleles that will be indispensable for designing salt-resilient crops.

Chromosome-level genome assembly and annotation of the prickly nightshade Solanum rostratum Dunal

The prickly nightshade Solanum rostratum, an annual malignant weed, is native to North America and has globally invaded 34 countries, causing serious threats to ecosystems, agriculture, animal husbandry, and human health. In this study, we constructed a chromosome-level genome assembly and annotation of S. rostratum. The contig-level genome was initially assembled in 898.42 Mb with a contig N50 of 62.00 Mb from PacBio high-fidelity reads. With Hi-C sequencing data scaffolding, 96.80% of the initially assembled sequences were anchored and orientated onto 12 pseudo-chromosomes, generating a genome of 869.69 Mb with a contig N50 of 72.15 Mb. We identified 649.92 Mb (72.26%) of repetitive sequences and 3,588 non-coding RNAs in the genome. A total of 29,694 protein-coding genes were predicted, with 28,154 (94.81%) functionally annotated genes. We found 99.5% and 91.3% complete embryophyta_odb10 genes in the pseudo-chromosomes genome and predicted gene datasets by BUSCO assessment. The present genomic resource provides essential information for subsequent research on the mechanisms of environmental adaptation of S. rostratum and host shift in Colorado potato beetles.

Laminarin-triggered defence responses are geographically dependent in natural populations of Solanum chilense

Natural plant populations are polymorphic and show intraspecific variation in resistance properties against pathogens. The activation of the underlying defence responses can depend on variation in perception of pathogen-associated molecular patterns or elicitors. To dissect such variation, we evaluated the responses induced by laminarin (a glucan, representing an elicitor from oomycetes) in the wild tomato species Solanum chilense and correlated this to observed infection frequencies of Phytophthora infestans. We measured reactive oxygen species burst and levels of diverse phytohormones upon elicitation in 83 plants originating from nine populations. We found high diversity in basal and elicitor-induced levels of each component. Further we generated linear models to explain the observed infection frequency of P. infestans. The effect of individual components differed dependent on the geographical origin of the plants. We found that the resistance in the southern coastal region, but not in the other regions, was directly correlated to ethylene responses and confirmed this positive correlation using ethylene inhibition assays. Our findings reveal high diversity in the strength of defence responses within a species and the involvement of different components with a quantitatively different contribution of individual components to resistance in geographically separated populations of a wild plant species.

Selective sweeps linked to the colonization of novel habitats and climatic changes in a wild tomato species

Positive selection is the driving force underpinning local adaptation and leaves footprints of selective sweeps on the underlying major genes. Quantifying the timing of selection and revealing the genetic bases of adaptation in plant species occurring in steep and varying environmental gradients are crucial to predict a species' ability to colonize new niches. We use whole-genome sequence data from six populations across three different habitats of the wild tomato species Solanum chilense to infer the past demographic history and search for genes under strong positive selection. We then correlate current and past climatic projections with the demographic history, allele frequencies, the age of selection events and distribution shifts. Several selective sweeps occur at regulatory networks involved in root-hair development in low altitude and response to photoperiod and vernalization in high-altitude populations. These sweeps appear to occur in a concerted fashion in a given regulatory gene network at particular periods of substantial climatic change. Using a unique combination of genome scans and modelling of past climatic data, we quantify the timing of selection at genes likely underpinning local adaptation to semiarid habitats.

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Reap the crop wild relatives for breeding future crops

Crop wild relatives (CWRs) have provided breeders with several 'game-changing' traits or genes that have boosted crop resilience and global agricultural production. Advances in breeding and genomics have accelerated the identification of valuable CWRs for use in crop improvement. The enhanced genetic diversity of breeding pools carrying optimum combinations of favorable alleles for targeted crop-growing regions is crucial to sustain genetic gain. In parallel, growing sequence information on wild genomes in combination with precise gene-editing tools provide a fast-track route to transform CWRs into ideal future crops. Data-informed germplasm collection and management strategies together with adequate policy support will be equally important to improve access to CWRs and their sustainable use to meet food and nutrition security targets.

Comparison of Drought and Heat Resistance Strategies among Six Populations of Solanum chilense and Two Cultivars of Solanum lycopersicum

Within the tomato clade, Solanum chilense is considered one of the most promising sources of genes for tomato (S. lycopersicum) selection to biotic and abiotic stresses. In this study, we compared the effects of drought, high temperature, and their combination in two cultivars of S. lycopersicum and six populations of S. chilense, differing in their local habitat. Plants were grown at 21/19 °C or 28/26 °C under well-watered and water-stressed conditions. Plant growth, physiological responses, and expression of stress-responsive genes were investigated. Our results demonstrated strong variability among accessions. Differences in plant growth parameters were even higher among S. chilense populations than between species. The effects of water stress, high temperature, and their combination also differed according to the accession, suggesting differences in stress resistance between species and populations. Overall, water stress affected plants more negatively than temperature from a morpho-physiological point of view, while the expression of stress-responsive genes was more affected by temperature than by water stress. Accessions clustered in two groups regarding resistance to water stress and high temperature. The sensitive group included the S. lycopersicum cultivars and the S. chilense populations LA2931 and LA1930, and the resistant group included the S. chilense populations LA1958, LA2880, LA2765, and LA4107. Our results suggested that resistance traits were not particularly related to the environmental conditions in the natural habitat of the populations. The expression of stress-responsive genes was more stable in resistant accessions than in sensitive ones in response to water stress and high temperature. Altogether, our results suggest that water stress and high temperature resistance in S. chilense did not depend on single traits but on a combination of morphological, physiological, and genetic traits.

De Novo Genome Assembly Of Solanum Sitiens Reveals Structural Variation Associated With Drought And Salinity Tolerance

MOTIVATION

Solanum sitiens is a self-incompatible wild relative of tomato, characterised by salt and drought resistance traits, with the potential to contribute through breeding programmes to crop improvement in cultivated tomato. This species has a distinct morphology, classification and ecotype compared to other stress resistant wild tomato relatives such as S. pennellii and S. chilense. Therefore, the availability of a reference genome for S. sitiens will facilitate the genetic and molecular understanding of salt and drought resistance.

RESULTS

A high-quality de novo genome and transcriptome assembly for S. sitiens (Accession LA1974) has been developed. A hybrid assembly strategy was followed using Illumina short reads (∼159X coverage) and PacBio long reads (∼44X coverage), generating a total of ∼262 Gbp of DNA sequence. A reference genome of 1,245 Mbp, arranged in 1,483 scaffolds with a N50 of 1.826 Mbp was generated. Genome completeness was estimated at 95% using the Benchmarking Universal Single-Copy Orthologs (BUSCO) and the K-mer Analysis Tool (KAT). In addition, ∼63 Gbp of RNA-Seq were generated to support the prediction of 31,164 genes from the assembly, and to perform a de novo transcriptome. Lastly, we identified three large inversions compared to S. lycopersicum, containing several drought resistance related genes, such as beta-amylase 1 and YUCCA7.

AVAILABILITY

S. sitiens (LA1974) raw sequencing, transcriptome and genome assembly have been deposited at the NCBI's Sequence Read Archive, under the BioProject number "PRJNA633104".All the commands and scripts necessary to generate the assembly are available at the following github repository: https://github.com/MCorentin/Solanum_sitiens_assembly.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multi-omics approach highlights differences between RLP classes in Arabidopsis thaliana

BACKGROUND

The Leucine rich-repeat (LRR) receptor-like protein (RLP) family is a complex gene family with 57 members in Arabidopsis thaliana. Some members of the RLP family are known to be involved in basal developmental processes, whereas others are involved in defence responses. However, functional data is currently only available for a small subset of RLPs, leaving the remaining ones classified as RLPs of unknown function.

RESULTS

Using publicly available datasets, we annotated RLPs of unknown function as either likely defence-related or likely fulfilling a more basal function in plants. Then, using these categories, we can identify important characteristics that differ between the RLP subclasses. We found that the two classes differ in abundance on both transcriptome and proteome level, physical clustering in the genome and putative interaction partners. However, the classes do not differ in the genetic di versity of their individual members in accessible pan-genome data.

CONCLUSIONS

Our work has several implications for work related to functional studies on RLPs as well as for the understanding of RLP gene family evolution. Using our annotations, we can make suggestions on which RLPs can be identified as potential immune receptors using genetics tools and thereby complement disease studies. The lack of differences in nucleotide diversity between the two RLP subclasses further suggests that non-synonymous diversity of gene sequences alone cannot distinguish defence from developmental genes. By contrast, differences in transcript and protein abundance or clustering at genomic loci might also allow for functional annotations and characterisation in other plant species.

Quantitative resistance differences between and within natural populations of Solanum chilense against the oomycete pathogen Phytophthora infestans

The wild tomato species Solanum chilense is divided into geographically and genetically distinct populations that show signs of defense gene selection and differential phenotypes when challenged with several phytopathogens, including the oomycete causal agent of late blight Phytophthora infestans. To better understand the phenotypic diversity of this disease resistance in S. chilense and to assess the effect of plant genotype versus pathogen isolate, respectively, we evaluated infection frequency in a systematic approach and with large sample sizes. We studied 85 genetically distinct individuals representing nine geographically separated populations of S. chilense. This showed that differences in quantitative resistance can be observed between but also within populations at the level of individual plants. Our data also did not reveal complete immunity in any of the genotypes. We further evaluated the resistance of a subset of the plants against P. infestans isolates with diverse virulence properties. This confirmed that the relative differences in resistance phenotypes between individuals were mainly determined by the plant genotype under consideration with modest effects of pathogen isolate used in the study. Thus, our report suggests that the observed quantitative resistance against P. infestans in natural populations of a wild tomato species S. chilense is the result of basal defense responses that depend on the host genotype and are pathogen isolate-unspecific.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

Limits and convergence properties of the sequentially Markovian coalescent

Several methods based on the sequentially Markovian coalescent (SMC) make use of full genome sequence data from samples to infer population demographic history including past changes in population size, admixture, migration events and population structure. More recently, the original theoretical framework has been extended to allow the simultaneous estimation of population size changes along with other life history traits such as selfing or seed banking. The latter developments enhance the applicability of SMC methods to nonmodel species. Although convergence proofs have been given using simulated data in a few specific cases, an in-depth investigation of the limitations of SMC methods is lacking. In order to explore such limits, we first develop a tool inferring the best case convergence of SMC methods assuming the true underlying coalescent genealogies are known. This tool can be used to quantify the amount and type of information that can be confidently retrieved from given data sets prior to the analysis of the real data. Second, we assess the inference accuracy when the assumptions of SMC approaches are violated due to departures from the model, namely the presence of transposable elements, variable recombination and mutation rates along the sequence, and SNP calling errors. Third, we deliver a new interpretation of SMC methods by highlighting the importance of the transition matrix, which we argue can be used as a set of summary statistics in other statistical inference methods, uncoupling the SMC from hidden Markov models (HMMs). We finally offer recommendations to better apply SMC methods and build adequate data sets under budget constraints.

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.

Population studies of the wild tomato species Solanum chilense reveal geographically structured major gene-mediated pathogen resistance

Natural plant populations encounter strong pathogen pressure and defence-associated genes are known to be under selection dependent on the pressure by the pathogens. Here, we use populations of the wild tomato Solanum chilense to investigate natural resistance against Cladosporium fulvum, a well-known ascomycete pathogen of domesticated tomatoes. Host populations used are from distinct geographical origins and share a defined evolutionary history. We show that distinct populations of S. chilense differ in resistance against the pathogen. Screening for major resistance gene-mediated pathogen recognition throughout the whole species showed clear geographical differences between populations and complete loss of pathogen recognition in the south of the species range. In addition, we observed high complexity in a homologues of Cladosporium resistance (Hcr) locus, underlying the recognition of C. fulvum, in central and northern populations. Our findings show that major gene-mediated recognition specificity is diverse in a natural plant-pathosystem. We place major gene resistance in a geographical context that also defined the evolutionary history of that species. Data suggest that the underlying loci are more complex than previously anticipated, with small-scale gene recombination being possibly responsible for maintaining balanced polymorphisms in the populations that experience pathogen pressure.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1 Mb of the S₂ -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S_6a -haplotype of P. inflata and the S_N -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

Alternative Splicing in the Anthocyanin Fruit Gene Encoding an R2R3 MYB Transcription Factor Affects Anthocyanin Biosynthesis in Tomato Fruits

Tomato (Solanum lycopersicum) fruits are typically red at ripening, with high levels of carotenoids and a low content in flavonoids. Considerable work has been done to enrich the spectrum of their health-beneficial phytochemicals, and interspecific crosses with wild species have successfully led to purple anthocyanin-colored fruits. The Aft (Anthocyanin fruit) tomato accession inherited from Solanum chilense the ability to accumulate anthocyanins in fruit peel through the introgression of loci controlling anthocyanin pigmentation, including four R2R3 MYB transcription factor-encoding genes. Here, we carried out a comparative functional analysis of these transcription factors in wild-type and Aft plants, and tested their ability to take part in the transcriptional complexes that regulate the biosynthetic pathway and their efficiency in inducing anthocyanin pigmentation. Significant differences emerged for SlAN2like, both in the expression level and protein functionality, with splicing mutations determining a complete loss of function of the wild-type protein. This transcription factor thus appears to play a key role in the anthocyanin fruit pigmentation. Our data provide new clues to the long-awaited genetic basis of the Aft phenotype and contribute to understand why domesticated tomato fruits display a homogeneous red coloration without the typical purple streaks observed in wild tomato species.