Comparative transcriptome analysis reveals resistant and susceptible genes in tobacco cultivars in response to infection by Phytophthora nicotianae

Phytophthora nicotianae is highly pathogenic to Solanaceous crops and is a major problem in tobacco production. The tobacco cultivar Beihart1000-1 (BH) is resistant, whereas the Xiaohuangjin 1025 (XHJ) cultivar is susceptible to infection. Here, BH and XHJ were used as models to identify resistant and susceptible genes using RNA sequencing (RNA-seq). Roots were sampled at 0, 6, 12, 24, and 60 h post infection. In total, 23,753 and 25,187 differentially expressed genes (DEGs) were identified in BH and XHJ, respectively. By mapping upregulated DEGs to the KEGG database, changes of the rich factor of “plant pathogen interaction pathway” were corresponded to the infection process. Of all the DEGs in this pathway, 38 were specifically regulated in BH. These genes included 11 disease-resistance proteins, 3 pathogenesis-related proteins, 4 RLP/RLKs, 2 CNGCs, 7 calcium-dependent protein kinases, 4 calcium-binding proteins, 1 mitogen-activated protein kinase kinase, 1 protein EDS1L, 2 WRKY transcription factors, 1 mannosyltransferase, and 1 calmodulin-like protein. By combining the analysis of reported susceptible (S) gene homologs and DEGs in XHJ, 9 S gene homologs were identified, which included 1 calmodulin-binding transcription activator, 1 cyclic nucleotide-gated ion channel, 1 protein trichome birefringence-like protein, 1 plant UBX domain-containing protein, 1 ADP-ribosylation factor GTPase-activating protein, 2 callose synthases, and 2 cellulose synthase A catalytic subunits. qRT-PCR was used to validate the RNA-seq data. The comprehensive transcriptome dataset described here, including candidate resistant and susceptible genes, will provide a valuable resource for breeding tobacco plants resistant to P. nicotianae infections.