Comparative Analysis of Multiple Consensus Genomes of the Same Strain of Marek’s Disease Virus Reveals Intrastrain Variation

Abstract

Current strategies to understand the molecular basis of Marek’s disease virus (MDV) virulence primarily consist of cataloguing divergent nucleotides between strains with different phenotypes. However, most comparative genomic studies of MDV rely on previously published consensus genomes despite the confirmed existence of MDV strains as mixed viral populations. To assess the reliability of interstrain genomic comparisons relying on published consensus genomes of MDV, we obtained two additional consensus genomes of vaccine strain CVI988 (Rispens) and two additional consensus genomes of the very virulent strain Md5 by sequencing viral stocks and cultured field isolates. In conjunction with the published genomes of CVI988 and Md5, this allowed us to perform 3-way comparisons between multiple consensus genomes of the same strain. We found that consensus genomes of CVI988 can vary in as many as 236 positions involving 13 open reading frames (ORFs). In contrast, we found that Md5 genomes varied only in 11 positions involving a single ORF. Notably, we were able to identify 3 SNPs in the Unique Long region and 16 SNPs in the Unique Short (US) region of CVI988GenBank.BAC that were not present in either CVI988Pirbright.lab or CVI988USDA.PA.field. Recombination analyses of field strains previously described as natural recombinants of CVI988 yielded no evidence of crossover events in the US region when either CVI988Pirbright.lab or CVI988USDA.PA.field were used to represent CVI988 instead of CVI988GenBank.BAC. We were also able to confirm that both CVI988 and Md5 populations were mixed, exhibiting a total of 29 and 27 high-confidence minor variant positions, respectively. However, we did not find any evidence of minor variants in the positions corresponding to the 19 SNPs in the unique regions of CVI988GenBank.BAC. Taken together, our findings suggest that continued reliance on the same published consensus genome of CVI988 may have led to an overestimation of genomic divergence between CVI988 and virulent strains, and that multiple consensus genomes per strain may be necessary to ensure the accuracy of interstrain genomic comparisons.