Genomic evaluation of Bordetella spp . originating from Australia

Abstract Bordetella pertussis is the primary causative agent of pertussis, a highly infectious respiratory disease associated with prolonged coughing episodes. Pertussis infections are typically mild in adults, however in neonates, infections can be fatal. Despite successful vaccine uptake, the disease is re-emerging across the globe, therefore it is critical to determine the mechanism by which B. pertussis is escaping vaccination control. Studies have suggested that significant changes have occurred in B. pertussis genomes in response to whole cell and acellular vaccines. Continued molecular monitoring is therefore crucial for public health surveillance. High-resolution molecular surveillance of B. pertussis can be achieved through the sequencing of the whole genome. In public health laboratories, whole genome sequencing is primarily performed by short-read sequencing technologies as they are most cost-effective. However short read sequencing does not resolve the extensive genomic rearrangement evident in Bordetella genomes. This is because repeat regions present in Bordetella genomes are collapsed by downstream analysis. For example, the B. pertussis genome contains more than 200 copies of the IS 481 insertion element, hence assemblies generally consist of >200 contigs. Advancements in long-read technologies however increase the potential to circularise and close genomes by bridging the locations of the IS 481 insertion element. In this study, we aimed to contextualise the Bordetella spp. circulating in NSW, Australia and assess their relationship with global isolates utilising core genome, SNP and structural clustering analysis using long read technology. We report five closed genomes of Bordetella spp. isolated from Australian patients. Two of the three B. pertussis closed isolates, were unique with their own genomic structure, while the other structurally clustered with global isolates. We found that Australian B. holmesii and B. parapertussis strains cluster with global isolates and do not appear to be unique to Australia. Australian draft B. holmesii SNP analysis showed that between 1999 and 2007, isolates were relatively similar, however post-2012, isolates were distinct from each other. The closed isolates can also be used as high-quality reference sequences for both surveillance and other investigations into pertussis spread.