Characterization of microRNAs in the cyst nematode Heterodera glycines identifies possible candidates involved in cross-kingdom interactions with its host Glycine max

The soybean cyst nematode (SCN – Heterodera glycines) is one of the most damaging pests to the cultivated soybean worldwide. Using a wide array of stylet-secreted effector proteins, this nematode can restructure its host cells into a complex and highly active feeding structure called the syncytium. Tight regulation of these proteins is thought to be essential to the successful formation of this syncytium. To date, multiple mechanisms have been proposed to regulate the expression of these proteins including through post-transcriptional regulation. MicroRNAs (miRNAs) are a class of small, roughly 22-nucleotide-long, non-coding RNA shown to regulate gene expression through its interaction with the 3’ untranslated region of genes. These same small RNAs have also been hypothesized to be able to cross over kingdom barriers and regulate genes in other species in a process called cross-kingdom interactions. In this study, we characterized the miRNome of the SCN via sequencing of small-RNAs isolated from whole nematodes and exosomes representing all developmental stages. We identified 121 miRNA loci encoding 96 distinct miRNA families including multiple lineage- and species-specific candidates. Using a combination of plant- and animal-specific miRNA target predictors, we generated a unique repertoire of miRNA:mRNA interacting partners in the nematode and its host plant leading to the identification of a set of nine probable cross-kingdom miRNA candidates.