Transcriptome sequencing reveals novel Citrus bark cracking viroid (CBCVd) variants from citrus and their molecular characterization

Citrus bark cracking viroid (CBCVd), previously called Citrus viroid IV, belongs to the genus Cocadviroid within the family Pospiviroidae. CBCVd has been identified as an important causative agent in citrus and hops. In this study, we obtained the full-length genomes of different variants of all detected citrus viroids from Pakistan through transcriptome sequencing. Different CBCVd variants were first found in Pakistan. These newly discovered Pakistani CBCVd variants were provisionally called “CBCVd-LSS” for their low sequence similarity (80.9%–88.9%) with the CBCVd RefSeq sequence (NC_003539). The two most predominant CBCVd sequences from Pakistan had the closest identity, 90.6% and 87.9%, with two CBCVd sequences isolated from hops. Identification and molecular characterization of CBCVd from citrus in Pakistan and China were also reported. The length of CBCVd from China ranged from 282 to 286 nucleotides, while that of the one from Pakistan ranged from 273 to 277 nucleotides. Based on genetic diversity and phylogenetic analysis, two main CBCVd clades were identified. CBCVd sequences from Pakistan, China, and other countries were further divided into six sub-clades. Sequence alignment revealed some nucleotide changes between these sub-clades, and analysis indicated that several mutations could significantly affect the primary and secondary structure of the viroid. Our results indicated that the CBCVd sequences from Pakistan and China were significantly different with respect to genome and secondary structure and Pakistan might be one of the independent geographical origins of CBCVd worldwide.