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 genomesWe 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.

Signatures of natural selection in abiotic stress-responsive genes of Solanum chilense.

Environmental conditions are strong selective forces, which may influence adaptation and speciation. The wild tomato species Solanum chilense, native to South America, is exposed to a range of abiotic stress factors. To identify signatures of natural selection and local adaptation, we analysed 16 genes involved in the abiotic stress response and compared the results to a set of reference genes in 23 populations across the entire species range. The abiotic stress-responsive genes are characterized by elevated nonsynonymous nucleotide diversity and divergence. We detected signatures of positive selection in several abiotic stress-responsive genes on both the population and species levels. Local adaptation to abiotic stresses is particularly apparent at the boundary of the species distribution in populations from coastal low-altitude and mountainous high-altitude regions.

Adaptation to low temperatures in the wild tomato species Solanum chilense.

Molecular adaptation to abiotic stresses in plants is a complex process based mainly on the modifications of gene transcriptional activity and the alteration of protein-protein interactions. We used a combination of population genetic, comparative transcriptomic and plant physiology approaches to investigate the mechanisms of adaptation to low temperatures in Solanum chilense populations distributed along Andean altitudinal gradients. We found that plants from all populations have high chilling tolerance, which does not correlate with temperatures in their native habitats. In contrast, tolerance to freezing shows a significant association with altitude and temperature variables. We also observed the differences in expression patterns of cold-response genes between plants from high- and low-altitude populations. These results suggest that genetic adaptations to low temperatures evolved in high-altitude populations of S. chilense. At the transcriptional level, these adaptations may include high levels of constitutive expression of the genes encoding ICE1, the key transcription factor of the cold signalling pathway, and chloroplast ω-3 fatty acid desaturase FAD7. At the sequence level, a signature of selection associated with the adaptation to high altitudes was detected at the C-terminal part of ICE1 encoding the ACT regulatory domain.

North-South Colonization Associated with Local Adaptation of the Wild Tomato Species Solanum chilense.

After colonization population sizes may vary across the species range depending on environmental conditions and following colonizations. An interesting question is whether local adaptation occurs more frequently in large ancestral populations or in small derived populations. A higher number of new mutations and a lower effect of genetic drift should favor selection in large populations, whereas small derived populations may require an initial local adaptation event to facilitate the colonization of new habitats. Wild tomatoes are native to a broad range of different habitats characterized by variable abiotic conditions in South America, and represent an ideal system to study this interplay between demography and natural selection. Population genetic analyses and statistical inference of past demography were conducted on pooled-sequencing data from 30 genes (8,080 single nucleotide polymorphisms) from an extensive sampling of 23 Solanum chilense populations over Chile and Peru. We reveal first a north-south colonization associated with relaxed purifying selection in the south as shown by a decrease of genetic variation and an increasing proportion of nonsynonymous polymorphism from north to south, and population substructure with at least four genetic groups. Second, we uncover a dual picture of adaptation consisting of 1) a decreasing proportion of adaptive amino acid substitutions from north to south suggesting that adaptation is favored in large populations, whereas 2) signatures of local adaptation predominantly occur in the smaller populations from the marginal ranges in the south.