There is no universally appropriate basis for delimiting species in protists, including parasites. Many molecular markers used for species delimitation are part of the ribosomal RNA (rRNA) array, with different regions of the array being used for different parasitic protist taxa. However, little is known about sequence variability across the rRNA in most organisms, and there is no standard threshold at which divergence in the sequence of a particular gene can be used as a basis for species delimitation.
Here we demonstrate a method to generate the full rRNA array of parasitic protists by amplification of the array in two long, overlapping fragments followed by Illumina and Nanopore sequencing to produce high quality assemblies, to determine variations in sequence variability across the array. We apply this approach to two pairs of closely related ascetosporean parasites of crustaceans and molluscs [respectively
Phylogenetic analyses of the transcribed regions of the rRNA array demonstrate maximal support for, and separation of, all four parasite species. Sliding window global alignment analysis determined the regions of the rRNA array that had the most consistent nucleotide differences between the closely related parasites in a 1 kb region of the array. For
Our method could be adapted to quickly generate sequence data and determine regions more suitable for diagnostic assays for a wide diversity of parasite groups. It also allows the generation of sequence data for regions of the rRNA not commonly studied (e.g. regions of the intergenic spacer), thus enabling research into their suitability as marker regions.