Identification of Stress Associated microRNAs in Solanum lycopersicum by High-Throughput Sequencing

Tomato (Solanum lycopersicum) is one of the most important crops around the world and also a model plant to study response to stress. High-throughput sequencing was used to analyse the microRNA (miRNA) profile of tomato plants undergoing five biotic and abiotic stress conditions (drought, heat, P. syringae infection, B. cinerea infection, and herbivore insect attack with Leptinotarsa decemlineata larvae) and one chemical treatment with a plant defence inducer, hexanoic acid. We identified 104 conserved miRNAs belonging to 37 families and we predicted 61 novel tomato miRNAs. Among those 165 miRNAs, 41 were stress-responsive. Reverse transcription quantitative PCR (RT-qPCR) was used to validate high-throughput expression analysis data, confirming the expression profiles of 10 out of 11 randomly selected miRNAs. Most of the differentially expressed miRNAs were stress-specific, except for sly-miR167c-3p upregulated in B. cinerea and P. syringae infection, sly-newmiR26-3p upregulated in drought and Hx treatment samples, and sly-newmiR33-3p, sly-newmiR6-3p and sly-newmiR8-3p differentially expressed both in biotic and abiotic stresses. From mature miRNAs sequences of the 41 stress-responsive miRNAs 279 targets were predicted. An inverse correlation between the expression profiles of 4 selected miRNAs (sly-miR171a, sly-miR172c, sly-newmiR22-3p and sly-miR167c-3p) and their target genes (Kinesin, PPR, GRAS40, ABC transporter, GDP and RLP1) was confirmed by RT-qPCR. Altogether, our analysis of miRNAs in different biotic and abiotic stress conditions highlight the interest to understand the functional role of miRNAs in tomato stress response as well as their putative targets which could help to elucidate plants molecular and physiological adaptation to stress.