Transcriptome analysis of a nematode resistant and susceptible upland cotton line at two critical stages of Meloidogyne incognita infection and development

Host plant resistance is the most practical approach to control the Southern root-knot nematode (Meloidogyne incognita; RKN), which has emerged as one of the most serious economic pests of Upland cotton (Gossypium hirsutum L.). Previous QTL analyses have identified a resistance locus on chromosome 11 (qMi-C11) affecting galling and another locus on chromosome-14 (qMi-C14) affecting egg production. Although these two QTL regions were fine mapped and candidate genes identified, expression profiling of genes would assist in further narrowing the list of candidate genes in the QTL regions. We applied the comparative transcriptomic approach to compare expression profiles of genes between RKN susceptible and resistance genotypes at an early stage of RKN development that coincides with the establishment of a feeding site and at the late stage of RKN development that coincides with RKN egg production. Sequencing of cDNA libraries produced over 315 million reads of which 240 million reads (76%) were mapped on to the Gossypium hirsutum genome. A total of 3,789 differentially expressed genes (DEGs) were identified which were further grouped into four clusters based on their expression profiles. A large number of DEGs were found to be down regulated in the susceptible genotype at the late stage of RKN development whereas several genes were up regulated in the resistant genotype. Key enriched categories included transcription factor activity, defense response, response to phyto-hormones, cell wall organization, and protein serine/threonine kinase activity. Our results also show that the DEGs in the resistant genotype at qMi-C11 and qMi-C14 loci displayed higher expression of defense response, detoxification and callose deposition genes, than the DEGs in the susceptible genotype.