Integrated RNA-seq and sRNA-seq revealed differences in transcriptome between susceptible and resistant tomato responding to Fusarium oxysporum

Tomato wilt disease caused by Fusarium oxysporum f. sp. lycopersici (FOL) is a worldwide destructive disease of tomato. As exploring gene expression and function approaches constitute an initial point for investigating pathogen-host interaction, we performed RNA-seq and sRNA-seq analysis to unravel regulated genes and miRNAs in tomato infected by FOL. Differentially expressed (DE) protein coding gene and miRNA gene profiles upon inoculation with FOL were presented at twenty-four hours post-inoculation including four treatments. Total of more than 182.6 million and 132.2 million high quality clean reads were obtained by RNA-seq and sRNA-seq, respectively. A large overlap was found in DE mRNAs between susceptible cultivar Moneymaker and resistant cultivar Motelle. All Gene Ontology terms were mainly classified into catalytic activity, metabolic process and binding. Combining with qRT-PCR, five disease resistance genes, Solyc01g095630, Solyc03g059080, Solyc00g174340, Solyc11g071750 and Solyc05g050350, were verified to involved in the disease resistance in the resistant cultivar Motelle treated with FOL. Northern blot analysis further confirmed the results from sRNA-Seq and demonstrated that several miRNAs including Sly-miR477-5p, sly-miR167a, novel_mir_675, novel_mir_504 and novel_mir_762 conferred FOL infection. Our data resulted that pathogen resistant genes/miRNAs may play a critical role with the benefit of a coordinated machinery in prompting the response in prompting FOL response in tomato, which offered us with a future direction and surely help in generating models of mediated resistance responses with assessment of genomic gene expression patterns.