Development and application of a qPCR-based genotyping assay for Ophidiomyces ophidiicola to investigate the epidemiology of ophidiomycosis

Abstract

Ophidiomycosis (snake fungal disease) is an infectious disease caused by the fungus Ophidiomyces ophidiicola to which all snake species appear to be susceptible. Significant variation has been observed in clinical presentation, progression of disease, and response to treatment, which may be due to genetic variation in the causative agent. Recent phylogenetic analysis based on whole-genome sequencing identified that O. ophidiicola strains from the United States formed a clade distinct from European strains, and that multiple clonal lineages of the clade are present in the United States. The purpose of this study was to design a qPCR-based genotyping assay for O. ophidiicola, then apply that assay to swabextracted DNA samples to investigate whether the multiple O. ophidiicola clades and clonal lineages in the United States have specific geographic, taxonomic, or temporal predilections. To this end, six full genome sequences of O. ophidiicola representing different clades and clonal lineages were aligned to identify genomic areas shared between subsets of the isolates. Eleven hydrolysis-based Taqman primer-probe sets were designed to amplify selected gene segments and produce unique amplification patterns for each isolate, each with a limit of detection of 10 or fewer copies of the target sequence and an amplification efficiency of 90-110%. The qPCR-based approach was validated using samples from strains known to belong to specific clades and applied to swab-extracted O. ophidiicola DNA samples from multiple snake species, states, and years. When compared to full-genome sequencing, the qPCR-based genotyping assay assigned 75% of samples to the same major clade (Cohen's kappa = 0.360, 95% Confidence Interval = 0.154-0.567) with 67-77% sensitivity and 88-100% specificity, depending on clade/clonal lineage. Swab-extracted O. ophidiicola DNA samples from across the United States were assigned to six different clonal lineages, including four of the six established lineages and two newly defined groups, which likely represent recombinant strains of O. ophidiicola. Using multinomial logistic regression modeling to predict clade based on snake taxonomic group, state of origin, and year of collection, state was the most significant predictor of clonal lineage. Furthermore, clonal lineage was not associated with disease severity in the most intensely sampled species, the Lake Erie watersnake (Nerodia sipedon insularum). Overall, this assay represents a rapid, cost-effective genotyping method for O. ophidiicola that can be used to better understand the epidemiology of ophidiomycosis.