Soybean cyst nematode (SCN) is one of the most important diseases in soybean. Currently, the main management strategy relies on planting resistant cultivars. However, the overuse of a single resistance source has led to the selection of virulent SCN populations, although the mechanisms by which the nematode overcomes the resistance genes remain unknown. In this study, we used a nematode-adapted single-cell RNA-seq approach to identify SCN genes potentially involved in resistance breakdown in Peking and PI 88788 parental soybean lines. We established for the first time the full transcriptome of single SCN individuals allowing us to identify a list of putative virulence genes against both major SCN resistance sources. Our analysis identified 48 differentially expressed putative effectors (secreted proteins required for infection) alongside 40 effectors showing evidence of novel structural variants, and 11 effector genes containing phenotype-specific sequence polymorphisms. Additionally, a differential expression analysis revealed an interesting phenomenon of coexpressed gene regions with some containing putative effectors. The selection of virulent SCN individuals on Peking resulted in a profoundly altered transcriptome, especially for genes known to be involved in parasitism. Several sequence polymorphisms were also specific to these virulent nematodes and could potentially play a role in the acquisition of nematode virulence. On the other hand, the transcriptome of virulent individuals on PI 88788 was very similar to avirulent ones with the exception of a few genes, which suggest a distinct virulence strategy to Peking.
Keywords: Heterodera glycines; PI 88788; Peking; RNA-seq; disease resistance; effector; evolution; genomics; host−parasite interactions; nematodes; pathogen effectors; plant-parasitic nematode; population biology; single-cell sequencing; virulence.