The wild tomato species Solanum chilense shows variation in pathogen resistance between geographically distinct populations

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.

Small-spored Alternaria spp. (section Alternaria) are common pathogens on wild tomato species

The wild relatives of modern tomato crops are native to South America. These plants occur in habitats as different as the Andes and the Atacama Desert and are, to some degree, all susceptible to fungal pathogens of the genus Alternaria. Alternaria is a large genus. On tomatoes, several species cause early blight, leaf spots and other diseases. We collected Alternaria-like infection lesions from the leaves of eight wild tomato species from Chile and Peru. Using molecular barcoding markers, we characterized the pathogens. The infection lesions were caused predominantly by small-spored species of Alternaria of the section Alternaria, like A. alternata, but also by Stemphylium spp., Alternaria spp. from the section Ulocladioides and other related species. Morphological observations and an infection assay confirmed this. Comparative genetic diversity analyses show a larger diversity in this wild system than in studies of cultivated Solanum species. As A. alternata has been reported to be an increasing problem in cultivated tomatoes, investigating the evolutionary potential of this pathogen is not only interesting to scientists studying wild plant pathosystems. It could also inform crop protection and breeding programs to be aware of potential epidemics caused by species still confined to South America.

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.

Distinct Phyllosphere Microbiome of Wild Tomato Species in Central Peru upon Dysbiosis

Plants are colonized by myriads of microbes across kingdoms, which affect host development, fitness, and reproduction. Hence, plant microbiomes have been explored across a broad range of host species, including model organisms, crops, and trees under controlled and natural conditions. Tomato is one of the world's most important vegetable crops; however, little is known about the microbiota of wild tomato species. To obtain insights into the tomato microbiota occurring in natural environments, we sampled epiphytic microbes from leaves of four tomato species, Solanum habrochaites, S. corneliomulleri, S. peruvianum, and S. pimpinellifolium, from two geographical locations within the Lima region of Peru over 2 consecutive years. Here, a high-throughput sequencing approach was applied to investigate microbial compositions including bacteria, fungi, and eukaryotes across tomato species and geographical locations. The phyllosphere microbiome composition varies between hosts and location. Yet, we identified persistent microbes across tomato species that form the tomato microbial core community. In addition, we phenotypically defined healthy and dysbiotic samples and performed a downstream analysis to reveal the impact on microbial community structures. To do so, we compared microbial diversities, unique OTUs, relative abundances of core taxa, and microbial hub taxa, as well as co-occurrence network characteristics in healthy and dysbiotic tomato leaves and found that dysbiosis affects the phyllosphere microbial composition in a host species-dependent manner. Yet, overall, the present data suggests an enrichment of plant-promoting microbial taxa in healthy leaves, whereas numerous microbial taxa containing plant pathogens occurred in dysbiotic leaves.Concluding, we identify the core phyllosphere microbiome of wild tomato species, and show that the overall phyllosphere microbiome can be impacted by sampling time point, geographical location, host genotype, and plant health. Future studies in these components will help understand the microbial contribution to plant health in natural systems and can be of use in cultivated tomatoes.

Edaphoclimatic Descriptors of Wild Tomato Species (Solanum Sect. Lycopersicon) and Closely Related Species (Solanum Sect. Juglandifolia and Sect. Lycopersicoides) in South America

Wild species related to cultivated tomato are essential genetic resources in breeding programs focused on food security to face future challenges. The ecogeographic analysis allows identifying the species adaptive ranges and most relevant environmental variables explaining their patterns of actual distribution. The objective of this research was to identify the diversity, ecological descriptors, and statistical relationship of 35 edaphoclimatic variables (20 climatic, 1 geographic and 14 edaphic variables) from 4,649 accessions of 12 wild tomato species and 4 closely related species classified in Solanum sect. Lycopersicon and clustered into four phylogenetic groups, namely “Lycopersicon group” (S. pimpinellifolium, S. cheesmaniae, and S. galapagense), “Arcanum group” (S. arcanum, S. chmielewskii, and S. neorickii), “Eriopersicon group” (S. habrochaites, S. huaylasense, S. corneliomulleri, S. peruvianum, and S. chilense), “Neolycopersicon group” (S. pennellii); and two phylogenetically related groups in Solanum sect. Juglandifolia (S. juglandifolium and S. ochranthum), and section Lycopersicoides (S. lycopersicoides and S. sitiens). The relationship between the climate and edaphic variables were determined by the canonical correlation analysis, reaching 89.2% of variation with the first three canonical correlations. The most significant climatic variables were related to humidity (annual evapotranspiration, annual precipitation, and precipitation of driest month) and physicochemical soil characteristics (bulk density, pH, and base saturation percentage). In all groups, ecological descriptors and diversity patterns were consistent with previous reports. Regarding edaphoclimatic diversity, 12 climate types and 17 soil units were identified among all species. This approach has promissory applications for biodiversity conservation and uses valuable genetic resources related to a leading crop.

Review: Holistic pest management against early blight disease towards sustainable agriculture

Alternaria species are well-known aggressive pathogens that are widespread globally and warmer temperatures caused by climate change might increase their abundance more drastically. Early blight (EB) disease, caused mainly by Alternaria solani, and brown spot, caused by Alternaria alternata, are major concerns in potato, tomato and eggplant production. The development of EB is strongly linked to varieties, crop development stages, environmental factors, cultivation and field management. Several forecasting models for pesticide application to control EB were created in the last century and more recent scientific advances have included modern breeding technology to detect resistant genes and precision agriculture with hyperspectral sensors to pinpoint damage locations on plants. This paper presents an overview of the EB disease and provides an evaluation of recent scientific advances to control the disease. First of all, we describe the outline of this disease, encompassing biological cycles of the Alternaria genus, favorite climate and soil conditions as well as resistant plant species. Second, versatile management practices to minimize the effect of this pathogen at field level are discussed, covering their limitations and pitfalls. A better understanding of the underlying factors of this disease and the potential of novel research can contribute to implementing integrated pest management systems for an ecofriendly farming system. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Polymorphisms in plants to restrict losses to pathogens: From gene family expansions to complex network evolution

Genetic polymorphisms are the basis of the natural diversity seen in all life on earth, also in plant-pathogen interactions. Initially, studies on plant-pathogen interaction focused on reporting phenotypic variation in resistance properties and on the identification of underlying major genes. Nowadays, the field of plant-pathogen interactions is moving from focusing on families of single dominant genes involved in gene-for-gene interactions to an understanding of the plant immune system in the context of a much more complex signaling network and quantitative resistance. Simultaneously, studies on pathosystems from the wild and genome analyses advanced, revealing tremendous variation in natural plant populations. It is now imperative to place studies on genetic diversity and evolution of plant-pathogen interactions in the appropriate molecular biological, as well as evolutionary, context.

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.

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.

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

Wild tomato species, like Solanum chilense, are important germplasm resources for enhanced biotic and abiotic stress resistance in tomato breeding. S. chilense also serves as a model to study adaptation of plants to drought and the evolution of seed banks. The absence of a well-annotated reference genome in this compulsory outcrossing, very diverse species limits in-depth studies on the genes involved.

We generated ∼134 Gb of DNA and 157 Gb of RNA sequence data for S chilense, which yielded a draft genome with an estimated length of 914 Mb, encoding 25,885 high-confidence predicted gene models, which show homology to known protein-coding genes of other tomato species. Approximately 71% of these gene models are supported by RNA-seq data derived from leaf tissue samples. Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis of predicted gene models retrieved 93.3% of BUSCO genes. To further verify the genome annotation completeness and accuracy, we manually inspected the NLR resistance gene family and assessed its assembly quality. We find subfamilies of NLRs unique to S. chilense. Synteny analysis suggests significant degree of the gene order conservation between the S. chilense, S. lycopersicum and S. pennellii genomes.

We generated the first genome and transcriptome sequence assemblies for the wild tomato species Solanum chilense and demonstrated their value in comparative genomics analyses. These data are an important resource for studies on adaptation to biotic and abiotic stress in Solanaceae, on evolution of self-incompatibility and for tomato breeding.

Subsets of NLR genes show differential signatures of adaptation during colonization of new habitats

Nucleotide binding site, leucine-rich repeat receptors (NLRs) are canonical resistance (R) genes in plants, fungi and animals, functioning as central (helper) and peripheral (sensor) genes in a signalling network. We investigate NLR evolution during the colonization of novel habitats in a model tomato species, Solanum chilense. We used R-gene enrichment sequencing to obtain polymorphism data at NLRs of 140 plants sampled across 14 populations covering the whole species range. We inferred the past demographic history of habitat colonization by resequencing whole genomes from three S. chilense plants from three key populations and performing approximate Bayesian computation using data from the 14 populations. Using these parameters, we simulated the genetic differentiation statistics distribution expected under neutral NLR evolution and identified small subsets of outlier NLRs exhibiting signatures of selection across populations. NLRs under selection between habitats are more often helper genes, whereas those showing signatures of adaptation in single populations are more often sensor-NLRs. Thus, centrality in the NLR network does not constrain NLR evolvability, and new mutations in central genes in the network are key for R-gene adaptation during colonization of different habitats.

A Molecular View of Plant Local Adaptation: Incorporating Stress-Response Networks

Ecological specialization in plants occurs primarily through local adaptation to different environments. Local adaptation is widely thought to result in costly fitness trade-offs that result in maladaptation to alternative environments. However, recent studies suggest that such trade-offs are not universal. Further, there is currently a limited understanding of the molecular mechanisms responsible for fitness trade-offs associated with adaptation. Here, we review the literature on stress responses in plants to identify potential mechanisms underlying local adaptation and ecological specialization. We focus on drought, high and low temperature, flooding, herbivore, and pathogen stresses. We then synthesize our findings with recent advances in the local adaptation and plant molecular biology literature. In the process, we identify mechanisms that could cause fitness trade-offs and outline scenarios where trade-offs are not a necessary consequence of adaptation. Future studies should aim to explicitly integrate molecular mechanisms into studies of local adaptation.

Gene Presence-Absence Polymorphism in Castrating Anther-Smut Fungi: Recent Gene Gains and Phylogeographic Structure

Gene presence-absence polymorphisms segregating within species are a significant source of genetic variation but have been little investigated to date in natural populations. In plant pathogens, the gain or loss of genes encoding proteins interacting directly with the host, such as secreted proteins, probably plays an important role in coevolution and local adaptation. We investigated gene presence-absence polymorphism in populations of two closely related species of castrating anther-smut fungi, Microbotryum lychnidis-dioicae (MvSl) and M. silenes-dioicae (MvSd), from across Europe, on the basis of Illumina genome sequencing data and high-quality genome references. We observed presence-absence polymorphism for 186 autosomal genes (2% of all genes) in MvSl, and only 51 autosomal genes in MvSd. Distinct genes displayed presence-absence polymorphism in the two species. Genes displaying presence-absence polymorphism were frequently located in subtelomeric and centromeric regions and close to repetitive elements, and comparison with outgroups indicated that most were present in a single species, being recently acquired through duplications in multiple-gene families. Gene presence-absence polymorphism in MvSl showed a phylogeographic structure corresponding to clusters detected based on SNPs. In addition, gene absence alleles were rare within species and skewed toward low-frequency variants. These findings are consistent with a deleterious or neutral effect for most gene presence-absence polymorphism. Some of the observed gene loss and gain events may however be adaptive, as suggested by the putative functions of the corresponding encoded proteins (e.g., secreted proteins) or their localization within previously identified selective sweeps. The adaptive roles in plant and anther-smut fungi interactions of candidate genes however need to be experimentally tested in future studies.