Background: Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock. Results: Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons. Conclusions: This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

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.