Transcriptome Analysis of Wheat–Tilletia indica Interaction Provides Defense and Pathogenesis-Related Genes

Abstract

Karnal bunt (Tilletia indica Mitra) is an internationally quarantined disease of wheat. Until now, very little information has been available on the molecular basis of resistance and pathogenicity of T. indica. To investigate the molecular basis of host–pathogen interaction, the transcriptome of T. indica inoculated resistant (HD29) and susceptible (WH542) genotypes of wheat were analyzed. Approximately 58 million reads were generated using RNA sequencing by the Illumina NextSeq500 platform. These sequence reads were aligned to a reference genome of wheat to compare the expression level of genes in resistant and susceptible genotypes. The high-quality reads were deposited in the NCBI SRA database (SRP159223). More than 80,000 genes were expressed in both the resistant and susceptible wheat genotypes. Of these, 76,088 were commonly expressed genes, including 3184 significantly upregulated and 1778 downregulated genes. Four thousand one hundred thirteen and 5604 genes were exclusively expressed in susceptible and resistant genotypes, respectively. Based on the significance, 503 genes were upregulated and 387 genes were downregulated. Using gene ontology, the majority of coding sequences were associated with response to stimuli, stress, carbohydrate metabolism, developmental process, and catalytic activity. Highly differentially expressed genes (integral component of membrane, exonuclease activity, nucleic acid binding, DNA binding, metal ion binding) were validated in resistant and susceptible genotypes using qPCR analysis and similar expression levels were found in RNA-Seq. Apart from the wheat, the mapping of T. indica was 7.07% and 7.63% of resistant and susceptible hosts, respectively, upon infection, which revealed significant pathogenesis-related genes. This first study provided in-depth information and new insights into wheat–T. indica interaction for managing Karnal bunt disease of wheat.