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

A two-sequence motif-based method for the inventory of gene families in fragmented and poorly annotated genome sequences

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H+-ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled.

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.

Crop Wild Relatives: A Valuable Source of Tolerance to Various Abiotic Stresses

Global climate change is one of the major constraints limiting plant growth, production, and sustainability worldwide. Moreover, breeding efforts in the past years have focused on improving certain favorable crop traits, leading to genetic bottlenecks. The use of crop wild relatives (CWRs) to expand genetic diversity and improve crop adaptability seems to be a promising and sustainable approach for crop improvement in the context of the ongoing climate challenges. In this review, we present the progress that has been achieved towards CWRs exploitation for enhanced resilience against major abiotic stressors (e.g., water deficiency, increased salinity, and extreme temperatures) in crops of high nutritional and economic value, such as tomato, legumes, and several woody perennial crops. The advances in -omics technologies have facilitated the elucidation of the molecular mechanisms that may underlie abiotic stress tolerance. Comparative analyses of whole genome sequencing (WGS) and transcriptomic profiling (RNA-seq) data between crops and their wild relative counterparts have unraveled important information with respect to the molecular basis of tolerance to abiotic stressors. These studies have uncovered genomic regions, specific stress-responsive genes, gene networks, and biochemical pathways associated with resilience to adverse conditions, such as heat, cold, drought, and salinity, and provide useful tools for the development of molecular markers to be used in breeding programs. CWRs constitute a highly valuable resource of genetic diversity, and by exploiting the full potential of this extended allele pool, new traits conferring abiotic-stress tolerance may be introgressed into cultivated varieties leading to superior and resilient genotypes. Future breeding programs may greatly benefit from CWRs utilization for overcoming crop production challenges arising from extreme environmental conditions.

Positive selection and heat-response transcriptomes reveal adaptive features of the Brassicaceae desert model, Anastatica hierochuntica

Plant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts. We generated an A. hierochuntica reference transcriptome and identified extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat, and detecting positively selected genes in A. hierochuntica. The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat 'stress-ready' state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype. We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.

A microRNA-microRNA crosstalk network inferred from genome-wide single nucleotide polymorphism variants in natural populations of Arabidopsis thaliana

To adapt to variable natural conditions, plants have evolved several strategies to respond to different environmental stresses. MicroRNA (miRNA)-mediated gene regulation is one of such strategies. Variants, e.g., single nucleotide polymorphisms (SNPs) within the mature miRNAs or their target sites may cause the alteration of regulatory networks and serious phenotype changes. In this study, we proposed a novel approach to construct a miRNA-miRNA crosstalk network in Arabidopsis thaliana based on the notion that two cooperative miRNAs toward common targets are under a strong pressure to be inherited together across ecotypes. By performing a genome-wide scan of the SNPs within the mature miRNAs and their target sites, we defined a "regulation fate profile" to describe a miRNA-target regulation being static (kept) or dynamic (gained or lost) across 1,135 ecotypes compared with the reference genome of Col-0. The cooperative miRNA pairs were identified by estimating the similarity of their regulation fate profiles toward the common targets. The reliability of the cooperative miRNA pairs was supported by solid expressional correlation, high PPImiRFS scores, and similar stress responses. Different combinations of static and dynamic miRNA-target regulations account for the cooperative miRNA pairs acting on various biological characteristics of miRNA conservation, expression, homology, and stress response. Interestingly, the targets that are co-regulated dynamically by both cooperative miRNAs are more likely to be responsive to stress. Hence, stress-related genes probably bear selective pressures in a certain group of ecotypes, in which miRNA regulations on the stress genes reprogram. Finally, three case studies showed that reprogramming miRNA-miRNA crosstalk toward the targets in specific ecotypes was associated with these ecotypes' climatic variables and geographical locations. Our study highlights the potential of miRNA-miRNA crosstalk as a genetic basis underlying environmental adaptation in natural populations.

The Halophyte Species Solanum chilense Dun. Maintains Its Reproduction despite Sodium Accumulation in Its Floral Organs

Salinity is a growing global concern that affects the yield of crop species, including tomato (Solanum lycopersicum). Its wild relative Solanum chilense was reported to have halophyte properties. We compared salt resistance of both species during the reproductive phase, with a special focus on sodium localization in the flowers. Plants were exposed to NaCl from the seedling stage. Salinity decreased the number of inflorescences in both species but the number of flowers per inflorescence and sepal length only in S. lycopersicum. External salt supply decreased the stamen length in S. chilense, and it was associated with a decrease in pollen production and an increase in pollen viability. Although the fruit set was not affected by salinity, fruit weight and size decreased in S. lycopersicum. Concentrations and localization of Na, K, Mg, and Ca differed in reproductive structures of both species. Inflorescences and fruits of S. chilense accumulated more Na than S. lycopersicum. Sodium was mainly located in male floral organs of S. chilense but in non-reproductive floral organs in S. lycopersicum. The expression of Na transporter genes differed in flowers of both species. Overall, our results indicated that S. chilense was more salt-resistant than S. lycopersicum during the reproductive phase and that differences could be partly related to dissimilarities in element distribution and transport in flowers.

Evolutionary footprints of a cold relic in a rapidly warming world

With accelerating global warming, understanding the evolutionary dynamics of plant adaptation to environmental change is increasingly urgent. Here, we reveal the enigmatic history of the genus Cochlearia (Brassicaceae), a Pleistocene relic that originated from a drought-adapted Mediterranean sister genus during the Miocene. Cochlearia rapidly diversified and adapted to circum-Arctic regions and other cold-characterized habitat types during the Pleistocene. This sudden change in ecological preferences was accompanied by a highly complex, reticulate polyploid evolution, which was apparently triggered by the impact of repeated Pleistocene glaciation cycles. Our results illustrate that two early diversified Arctic-alpine diploid gene pools contributed differently to the evolution of this young polyploid genus now captured in a cold-adapted niche. Metabolomics revealed central carbon metabolism responses to cold in diverse species and ecotypes, likely due to continuous connections to cold habitats that may have facilitated widespread adaptation to alpine and subalpine habitats, and which we speculate were coopted from existing drought adaptations. Given the growing scientific interest in the adaptive evolution of temperature-related traits, our results provide much-needed taxonomic and phylogenomic resolution of a model system as well as first insights into the origins of its adaptation to cold.

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.

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.

Fungal and Bacterial Microbiome Associated with the Rhizosphere of Native Plants from the Atacama Desert

The rhizosphere microbiome is key in survival, development, and stress tolerance in plants. Salinity, drought, and extreme temperatures are frequent events in the Atacama Desert, considered the driest in the world. However, little information of the rhizosphere microbiome and its possible contribution to the adaptation and tolerance of plants that inhabit the desert is available. We used a high-throughput Illumina MiSeq sequencing approach to explore the composition, diversity, and functions of fungal and bacterial communities of the rhizosphere of Baccharis scandens and Solanum chilense native plants from the Atacama Desert. Our results showed that the fungal phyla Ascomycota and Basidiomycota and the bacterial phyla Actinobacteria and Proteobacteria were the dominant taxa in the rhizosphere of both plants. The linear discriminant analysis (LDA) effect size (LefSe) of the rhizosphere communities associated with B. scandens showed the genera Penicillium and Arthrobacter were the preferential taxa, whereas the genera Oidiodendron and Nitrospirae was the preferential taxa in S. chilense. Both plant showed similar diversity, richness, and abundance according to Shannon index, observed OTUs, and evenness. Our results indicate that there are no significant differences (p = 0.1) between the fungal and bacterial communities of both plants, however through LefSe, we find taxa associated with each plant species and the PCoA shows a separation between the samples of each species. This study provides knowledge to relate the assembly of the microbiome to the adaptability to drought stress in desert plants.

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.

Effectiveness of Genomic Selection by Response to Selection for Winter Wheat Variety Improvement

Prediction performance for winter wheat grain yield and end-use quality traits. Prediction accuracies evaluated by cross-validations are significantly overestimated. Nonparametric algorithms outperform the parametric alternatives in cross-year predictions. Strategically designing training population improves response to selection. Response to selection varies across growing seasons and environments. Considering the practicality of applying genomic selection (GS) in the line development stage of a hard red winter (HRW) wheat (Triticum aestivum L.) variety development program (VDP), the effectiveness of GS was evaluated by prediction accuracy and by the response to selection across field seasons that demonstrated challenges for crop improvement under significant climate variability. Important breeding targets for wheat improvement in the southern Great Plains of the United States, including grain yield, kernel weight, wheat protein content, and sodium dodecyl sulfate (SDS) sedimentation volume as a rapid test for predicting bread-making quality, were used to estimate the effectiveness of GS across harvest years from 2014 (drought) to 2016 (normal). In general, nonparametric algorithms reproducing kernel Hilbert space (RKHS) and random forest (RF) produced higher accuracies in both same-year cross-validations (CVs) and cross-year predictions for the purpose of line selection. Further, the stability of GS performance was greatest for SDS sedimentation volume but least for wheat protein content. To ensure long-term genetic gain, our study on selection response suggested that across this sample of environmental variability, and though there are cases where phenotypic selection (PS) might be still preferred, training conducted under drought or in suboptimal conditions could provide an encouraging prediction outcome when selection decisions were made in normal conditions. However, it is not advisable to use training information collected from a normal season to predict trait performance under drought conditions. Finally, the superiority of response to selection was most evident if the training population (TP) can be optimized.

Identification and transcriptional analysis of SNARE vesicle fusion regulators in tomato (Solanum lycopersicum) during plant development and a comparative analysis of the response to salt stress with wild relatives

Intracellular vesicular trafficking ensures the exchange of lipids and proteins between the membranous compartments. Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNARE) play a central role in membrane fusion and they are key factors for vesicular trafficking in plants, including crops economically important such as tomato (Solanum lycopersicum). Taking advantage of the complete genome sequence available of S. lycopersicum, we identified 63 genes that encode putative SNARE proteins. Then, phylogenetic analysis allowed the classification of SNAREs in five main groups and recognizing their possible functions. A structure analysis of the genes, their syntenic relationships and their location in the chromosomes were also carried out for their characterization. In addition, the expression profiles of SNARE genes in different tissues were investigated using microarray-based analysis. The results indicated that specific SNAREs had a higher induction in leaf, root, flower and mature green fruit. S. lycopersicum is characterized for being a crop sensitive to saline stress unlike its wild relatives, such as Solanum pennellii, Solanum pimpinellifolium, Solanum habrochaites or Solanum chilense, which are tolerant. In this context, we analyzed different microarrays and evaluated and validated the transcript levels through qRT-PCR experiments. The results showed that SlGOS12.2, SlVAMP727 and SlSYP51.2 could have a positive relationship with salt stress and probably an important role in their tolerance. All these data increase our knowledge and can also be utilized to identify potential molecular targets for conferring tolerance to various stresses in tomato.

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.

Transgenic tomatoes for abiotic stress tolerance: status and way ahead

Tomato (Solanum lycopersicum) is one of the most important vegetable crops; its production, productivity and quality are adversely affected by abiotic stresses. Abiotic stresses such as drought, extreme temperature and high salinity affect almost every stage of tomato life cycle. Depending upon the plant stage and duration of the stress, abiotic stress causes about 70% yield loss. Several wild tomato species have the stress tolerance genes; however, it is very difficult to transfer them into cultivars due to high genetic distance and crossing barriers. Transgenic technology is an alternative potential tool for the improvement of tomato crop to cope with abiotic stress, as it allows gene transfer across species. In recent decades, many transgenic tomatoes have been developed, and many more are under progress against abiotic stress using transgenes such as DREBs, Osmotin, ZAT12 and BADH2. The altered expression of these transgenes under abiotic stresses are involved in every step of stress responses, such as signaling, control of transcription, proteins and membrane protection, compatible solute (betaines, sugars, polyols, and amino acids) synthesis, and free-radical and toxic-compound scavenging. The stress-tolerant transgenic tomato development is based on introgression of a gene with known function in stress response and putative tolerance. Transgenic tomato plants have been developed against drought, heat and salt stress with the help of various transgenes, expression of which manages the stress at the cellular level by modulating the expression of downstream genes to ultimately improve growth and yield of tomato plants and help in sustainable agricultural production. The transgenic technology could be a faster way towards tomato improvement against abiotic stress. This review provides comprehensive information about transgenic tomato development against abiotic stress such as drought, heat and salinity for researcher attention and a better understanding of transgenic technology used in tomato improvement and sustainable agricultural production.