The Bursaphelenchus genus (Nematoda: Parasitaphelenchidae) comprises mostly wood-inhabiting nematodes that feed on various tree-colonizing fungi. One species of the genus, B. xylophilus, has been proven as an agent causing pine wilt disease (PWD). However, involvement of other Bursaphelenchus species in the PWD remains enigmatic. In the current paper, comparative molecular analysis is performed based on nuclear ribosomal DNA (rDNA) of B. vallesianus, a species that was recently isolated from pine trees (Pinus sylvestris) exhibiting wilting and declining symptoms in the Czech Republic. Sequencing of the nuclear-encoded ITS1-5Á8S-ITS2 rDNA region confirmed previous taxonomic conclusions based on morphology. Evolutionary reconstructions resulted in a phylogenetic tree, where the Czech isolate of B. vallesianus occupied a common clade together with other species belonging to the so-called B. sexdentati group. Unexpectedly, comprehensive analysis of the sequence data revealed a genetic variation distinguishing the Czech isolate of B. vallesianus from all other species of the B. sexdentati group. This dissimilarity consists of the presence of a four nucleotide exchange found in the 5Á8S rRNA-coding gene. The newly identified genetic variation appears to affect the 5Á8S rRNA folding, as deduced from secondary structure models. Additionally, it is shown that for the first time, to the authors' knowledge, both bursaphelenchid internal transcribed spacers (ITS1 and ITS2) fold into the multibranched closed loops. While the ITS2 closed loop is formed with help of canonical 5Á8S-28S rRNA pairing, the ITS1 forms the thermodynamically stable closed loop with no support of flanking rRNA sequences. The current information on bursaphelenchid ITS rDNA sequence diversity and structure is further discussed.