Tularemia Transmission to Humans, the Netherlands, 2011-2021.

We used national registry data on human cases of Francisella tularensis subspecies holarctica infection to assess transmission modes among all 26 autochthonous cases in the Netherlands since 2011. The results indicate predominance of terrestrial over aquatic animal transmission sources. We recommend targeting disease-risk communication toward hunters, recreationists, and outdoor professionals.

Signalling and responding to zoonotic threats using a One Health approach: a decade of the Zoonoses Structure in the Netherlands, 2011 to 2021.

In the Netherlands, the avian influenza outbreak in poultry in 2003 and the Q fever outbreak in dairy goats between 2007 and 2010 had severe consequences for public health. These outbreaks led to the establishment of an integrated human-veterinary risk analysis system for zoonoses, the Zoonoses Structure. The aim of the Zoonoses Structure is to signal, assess and control emerging zoonoses that may pose a risk to animal and/or human health in an integrated One Health approach. The Signalling Forum Zoonoses (SO-Z), the first step of the Zoonoses Structure, is a multidisciplinary committee composed of experts from the medical, veterinary, entomology and wildlife domains. The SO-Z shares relevant signals with professionals and has monthly meetings. Over the past 10 years (June 2011 to December 2021), 390 different signals of various zoonotic pathogens in animal reservoirs and humans have been assessed. Here, we describe the Zoonoses Structure with examples from signals and responses for four zoonotic events in the Netherlands (tularaemia, Brucella canis, West Nile virus, and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)). This may serve as an example for other countries on how to collaborate in a One Health approach to signal and control emerging zoonoses.

Investigation of Clostridium botulinum group III's mobilome content.

Clostridium botulinum group III is mainly responsible for botulism in animals. It could lead to high animal mortality rates and, therefore, represents a major environmental and economic concern. Strains of this group harbor the botulinum toxin locus on an unstable bacteriophage. Since the release of the first complete C. botulinum group III genome sequence (strain BKT015925), strains have been found to contain others mobile elements encoding for toxin components. In this study, seven assays targeting toxin genes present on the genetic mobile elements of C. botulinum group III were developed with the objective to better characterize C. botulinum group III strains. The investigation of 110 C. botulinum group III strains and 519 naturally contaminated samples collected during botulism outbreaks in Europe showed alpha-toxin and C2-I/C2-II markers to be systematically associated with type C/D bont-positive samples, which may indicate an important role of these elements in the pathogenicity mechanisms. On the contrary, bont type D/C strains and the related positive samples appeared to contain almost none of the markers tested. Interestingly, 31 bont-negative samples collected on farms after a botulism outbreak revealed to be positive for some of the genetic mobile elements tested. This suggests loss of the bont phage, either in farm environment after the outbreak or during laboratory handling.

Evaluation of a multiplex real-time PCR for detection of four bacterial agents commonly associated with bovine respiratory disease in bronchoalveolar lavage fluid.

BACKGROUND: Pasteurella multocida, Mannheimia haemolytica, Histophilus somni and Trueperella pyogenes are four bacterial agents commonly associated with bovine respiratory disease (BRD). In this study a bacterial multiplex real-time PCR (the RespoCheck PCR) was evaluated for the detection in bronchoalveolar lavage fluid (BALF) of these four bacterial agents. RESULTS: The analytical sensitivity of the multiplex real-time PCR assay determined on purified DNA and on bacterial cells of the four target pathogens was one to ten fg DNA/assay and 4 × 10-1 to 2 × 100 CFU/assay. The analytical specificity of the test was, as evaluated on a collection of 118 bacterial isolates, 98.3% for M. haemolytica and 100% for the other three target bacteria. A set of 160 BALF samples of calves originating from ten different herds with health problems related to BRD was examined with bacteriological methods and with the RespoCheck PCR. Using bacteriological examination as the gold standard, the diagnostic sensitivities and specificities of the four bacterial agents were respectively between 0.72 and 1.00 and between 0.70 and 0.99. Kappa values for agreement between results of bacteriological examination and PCRs were low for H. somni (0.17), moderate for P. multocida (0.52) and M. haemolytica (0.57), and good for T. pyogenes (0.79). The low and moderate kappa values seemed to be related to limitations of the bacteriological examination, this was especially the case for H. somni. CONCLUSION: It was concluded that the RespoCheck PCR assay is a valuable diagnostic tool for the simultaneous detection of the four bacterial agents in BALF of calves.

Mycoplasma detection by triplex real-time PCR in bronchoalveolar lavage fluid from bovine respiratory disease complex cases.

BACKGROUND: In this study we evaluated the RespoCheck Mycoplasma triplex real-time PCR for the detection in bronchoalveolar lavage fluid (BALF) of Mycoplasma (M.) dispar, M. bovis and M. bovirhinis, all three associated with bovine respiratory disease (BRD). Primers and probes of the RespoCheck Mycoplasma triplex real-time PCR are based on the V3/V4 region of the 16S rRNA gene of the three Mycoplasma species. RESULTS: The analytical sensitivity of the RespoCheck triplex real-time PCR was, as determined by spiking experiments of the Mycoplasma strains in Phosphate Buffered Saline, 300 colony forming units (cfu)/mL for M. dispar, and 30 cfu/mL for M. bovis or M. bovirhinis. The analytical sensitivity of the RespoCheck Mycoplasma triplex real-time PCRwas, as determined on purified DNA, 10 fg DNA per assay for M. dispar and 100 fg fo rM. bovis and M. bovirhinis. The analytical specificity of the RespoCheck Mycoplasma triplex real-time PCR was, as determined by testing Mycoplasmas strains (n = 17) and other bacterial strains (n = 107), 100, 98.2 and 99.1% for M. bovis, M. dispar and M. bovirhinis respectively. The RespoCheck Mycoplasma triplex real-time PCR was compared with the PCR/DGGE analysis for M. bovis, M. dispar and M. bovirhinis respectively by testing 44 BALF samples from calves. CONCLUSION: In conclusion, the RespoCheck PCR assay can be a valuable tool for timely and accurate detection of three Mycoplasma species associated with in bovine respiratory disease.

Environmental surveillance during an outbreak of tularaemia in hares, the Netherlands, 2015.

Tularaemia, a disease caused by the bacterium Francisella tularensis, is a re-emerging zoonosis in the Netherlands. After sporadic human and hare cases occurred in the period 2011 to 2014, a cluster of F. tularensis-infected hares was recognised in a region in the north of the Netherlands from February to May 2015. No human cases were identified, including after active case finding. Presence of F. tularensis was investigated in potential reservoirs and transmission routes, including common voles, arthropod vectors and surface waters. F. tularensis was not detected in common voles, mosquito larvae or adults, tabanids or ticks. However, the bacterium was detected in water and sediment samples collected in a limited geographical area where infected hares had also been found. These results demonstrate that water monitoring could provide valuable information regarding F. tularensis spread and persistence, and should be used in addition to disease surveillance in wildlife.

Molecular gene profiling of Clostridium botulinum group III and its detection in naturally contaminated samples originating from various European countries.

We report the development of real-time PCR assays for genotyping Clostridium botulinum group III targeting the newly defined C. novyi sensu lato group; the nontoxic nonhemagglutinin (NTNH)-encoding gene ntnh; the botulinum neurotoxin (BoNT)-encoding genes bont/C, bont/C/D, bont/D, and bont/D/C; and the flagellin (fliC) gene. The genetic diversity of fliC among C. botulinum group III strains resulted in the definition of five major subgroups named fliC-I to fliC-V. Investigation of fliC subtypes in 560 samples, with various European origins, showed that fliC-I was predominant and found exclusively in samples contaminated by C. botulinum type C/D, fliC-II was rarely detected, no sample was recorded as fliC-III or fliC-V, and only C. botulinum type D/C samples tested positive for fliC-IV. The lack of genetic diversity of the flagellin gene of C. botulinum type C/D would support a clonal spread of type C/D strains in different geographical areas. fliC-I to fliC-III are genetically related (87% to 92% sequence identity), whereas fliC-IV from C. botulinum type D/C is more genetically distant from the other fliC types (with only 50% sequence identity). These findings suggest fliC-I to fliC-III have evolved in a common environment and support a different genetic evolution for fliC-IV. A combination of the C. novyi sensu lato, ntnh, bont, and fliC PCR assays developed in this study allowed better characterization of C. botulinum group III and showed the group to be less genetically diverse than C. botulinum groups I and II, supporting a slow genetic evolution of the strains belonging to C. botulinum group III.

Campylobacter presence on Dutch broiler farms and associated risk factors.

Campylobacter is the most reported zoonotic pathogen in humans in the European Union. Poultry is a major source of human infection with Campylobacter. Although many studies are done on the presence of Campylobacter in broilers and theoretically effective control measures are known, their relative importance at broiler farms remains poorly understood. Therefore, the aim of this study was to investigate the presence of Campylobacter on selected broiler farms in the Netherlands, to determine the moment of introduction, and associated risk factors. A longitudinal study on 25 broiler farms was carried out between June 2017 and December 2020. Fecal samples were collected weekly from 43 broiler houses. In total 497 flocks were sampled. Putative variables on flock and farm characteristics for a risk factor analysis were gathered through questionnaires. Risk factors associated with the presence of Campylobacter in a broiler flock were determined using regression models. In total 30% of the flocks included in the study were positive for Campylobacter. Factors associated with presence of Campylobacter at slaughter age included: season, mowing lawns and presence of agricultural side activities. While summer/autumn and mowing lawns were associated with an increase in Campylobacter presence in flocks, the farmer having agricultural side activities other than poultry production was associated with a decrease. Analysis of the age at which flocks first tested Campylobacter positive revealed that slower growing breeds became positive on average 1 wk later compared to regular growers. This study revealed a delayed introduction of Campylobacter in slower grower vs. regular grower broiler flocks reared indoors. In addition, it confirmed importance of season as major risk factor. The relevance of mowing and preceding positive flocks as risk factors needs further investigation.

Validation of a real-time PCR based method for detection of Clostridium botulinum types C, D and their mosaic variants C-D and D-C in a multicenter collaborative trial.

Two real-time PCR arrays based on the GeneDisc(®) cycler platform (Pall-GeneDisc Technologies) were evaluated in a multicenter collaborative trial for their capacity to specifically detect and discriminate Clostridium botulinum types C, D and their mosaic variants C-D and D-C that are associated with avian and mammalian botulism. The GeneDisc(®) arrays developed as part of the DG Home funded European project 'AnibioThreat' were highly sensitive and specific when tested on pure isolates and naturally contaminated samples (mostly clinical specimen from avian origin). Results of the multicenter collaborative trial involving eight laboratories in five European Countries (two laboratories in France, Italy and The Netherlands, one laboratory in Denmark and Sweden), using DNA extracts issued from 33 pure isolates and 48 naturally contaminated samples associated with animal botulism cases, demonstrated the robustness of these tests. Results showed a concordance among the eight laboratories of 99.4%-100% for both arrays. The reproducibility of the tests was high with a relative standard deviation ranging from 1.1% to 7.1%. Considering the high level of agreement achieved between the laboratories these PCR arrays constitute robust and suitable tools for rapid detection of C. botulinum types C, D and mosaic types C-D and D-C. These are the first tests for C. botulinum C and D that have been evaluated in a European multicenter collaborative trial.

Environmental Sampling Methods for Detection of Salmonella Infections in Laying Hens: A Systematic Review and Meta-Analysis.

Salmonellosis is the second most commonly reported foodborne gastrointestinal infection in humans in the European Union (EU). Most outbreaks are caused by Salmonella Enteritidis, present in contaminated food products, particularly in egg and egg products. In recent years, an increase in the prevalence of Salmonella in laying hen flocks in the EU has been observed. For the effective control of infection, adequate detection is key. In laying hen flocks, the occurrence of Salmonella in the EU is monitored by the culture of environmental samples (dust, faeces, and boot swabs). The performance of sampling procedures described in the literature for the detection of Salmonella in laying hens was reviewed. In total, 924 abstracts were screened, resulting in the selection of 87 abstracts and 18 publications for qualitative and quantitative analyses, respectively. Sample sizes and sampling locations of faecal material and dust were variable and poorly described. Microbiological culture methods used to detect Salmonella were variably described in the literature and were often incomplete. Overall, the available literature indicates higher sensitivity of environmental versus individual hen matrices and points to differences in sensitivity between environmental matrices. For non-cage housing systems, boot swabs are the preferred samples, while for cage housing systems dust might be a more reliable sample.

Combining a parsimonious mathematical model with infection data from tailor-made experiments to understand environmental transmission.

Although most infections are transmitted through the environment, the processes underlying the environmental stage of transmission are still poorly understood for most systems. Improved understanding of the environmental transmission dynamics is important for effective non-pharmaceutical intervention strategies. To study the mechanisms underlying environmental transmission we formulated a parsimonious modelling framework including hypothesised mechanisms of pathogen dispersion and decay. To calibrate and validate the model, we conducted a series of experiments studying distance-dependent transmission of Campylobacter jejuni in broilers. We obtained informative simultaneous estimates for all three model parameters: the parameter of C. jejuni inactivation, the diffusion coefficient describing pathogen dispersion, and the transmission rate parameter. The time and distance dependence of transmission in the fitted model is quantitatively consistent with marked spatiotemporal patterns in the experimental observations. These results, for C. jejuni in broilers, show that the application of our modelling framework to suitable transmission data can provide mechanistic insight in environmental pathogen transmission.

Single-Domain Antibody Multimers for Detection of Botulinum Neurotoxin Serotypes C, D, and Their Mosaics in Endopep-MS.

Botulinum neurotoxins (BoNTs) are highly toxic proteins that require high-affinity immunocapture reagents for use in endopeptidase-based assays. Here, 30 novel and 2 earlier published llama single-domain antibodies (VHHs) against the veterinary-relevant BoNT serotypes C and D were yeast-produced. These VHHs recognized 10 independent antigenic sites, and many cross-reacted with the BoNT/DC and CD mosaic variants. As VHHs are highly suitable for genetically linking to increase antigen-binding affinity, 52 VHH multimers were produced and their affinity for BoNT/C, D, DC, and CD was determined. A selection of 15 multimers with high affinity (KD < 0.1 nM) was further shown to be resilient to a high salt wash that is used for samples from complex matrices and bound native BoNTs from culture supernatants as shown by Endopep-MS. High-affinity multimers suitable for further development of a highly sensitive Endopep-MS assay include four multimers that bind both BoNT/D and CD with KD of 14-99 pM, one multimer for BoNT/DC (65 pM) that also binds BoNT/C (75 pM), and seven multimers for BoNT/C (<1-19 pM), six of which also bind BoNT/DC with lower affinity (93-508 pM). In addition to application in diagnostic tests, these VHHs could be used for the development of novel therapeutics for animals or humans.

Neurotoxin gene profiling of clostridium botulinum types C and D native to different countries within Europe.

Clostridium botulinum types C and D, as well as their mosaic variants C-D and D-C, are associated with avian and mammalian botulism. This study reports on the development of low-density macroarrays based on the GeneDisc cycler platform (Pall-GeneDisc Technologies) applied to the simultaneous detection of the C. botulinum subtypes C, C-D, D, and D-C. The limit of detection of the PCR assays was 38 fg of total DNA, corresponding to 15 genome copies. Artificially contaminated samples of cecum showed a limit of detection below 50 spores/g. The tests were performed with a large variety of bacterial strains, including C. botulinum types C (n = 12), C-D (n = 29), D (n = 5), and D-C (n = 10), other botulinum neurotoxin (BoNT)-producing Clostridium strains (n = 20), non-BoNT-producing clostridia (n = 20), and other bacterial species (n = 23), and showed a high specificity. These PCR assays were compared to previously published real-time PCRs for the detection of C. botulinum in 292 samples collected from cases of botulism events in four European regions. The majority of the samples originated from wild birds (n = 108), poultry (n = 60), and bovines (n = 56). Among the 292 samples, 144 were positive for either the bont/C-D or the bont/D-C gene by using the GeneDisc arrays. The reliability of the results tallied to 97.94%. Interestingly, only BoNT mosaics, types C-D and D-C, were found in naturally contaminated samples whatever their animal origin and their geographical location. Further investigations should now be performed in order to check that mosaic types dominate in Europe and that acquisition of mosaic types helps in survival or adaptation to particular niche.

Serum protein profiles as potential biomarkers for infectious disease status in pigs.

BACKGROUND: In veterinary medicine and animal husbandry, there is a need for tools allowing the early warning of diseases. Preferably, tests should be available that warn farmers and veterinarians during the incubation periods of disease and before the onset of clinical signs. The objective of this study was to explore the potential of serum protein profiles as an early biomarker for infectious disease status. Serum samples were obtained from an experimental pig model for porcine circovirus-associated disease (PCVAD), consisting of Porcine Circovirus type 2 (PCV2) infection in combination with either Porcine Parvovirus (PPV) or Porcine Reproductive and Respiratory Syndrome virus (PRRSV). Sera were collected before and after onset of clinical signs at day 0, 5 and 19 post infection. Serum protein profiles were evaluated against sera from non-infected control animals. RESULTS: Protein profiles were generated by SELDI-TOF mass spectrometry in combination with the Proteominer™ technology to enrich for low-abundance proteins. Based on these protein profiles, the experimentally infected pigs could be classified according to their infectious disease status. Before the onset of clinical signs 88% of the infected animals could be classified correctly, after the onset of clinical sigs 93%. The sensitivity of the classification appeared to be high. The protein profiles could distinguish between separate infection models, although specificity was moderate to low. Classification of PCV2/PRRSV infected animals was superior compared to PCV2/PPV infected animals. Limiting the number of proteins in the profiles (ranging from 568 to 10) had only minor effects on the classification performance. CONCLUSIONS: This study shows that serum protein profiles have potential for detection and identification of viral infections in pigs before clinical signs of the disease become visible.

Low seroprevalence of equine piroplasmosis in horses exported from the Netherlands between 2015 and 2021.

Equine piroplasmosis (EP) is a tick-borne disease affecting horses, donkeys, mules and zebras, caused by the intracellular apicomplexan protozoa Babesia caballi and Theileria equi. The geographical distribution of EP is closely related to the distribution of its vector tick species belonging to the genera of Dermacentor, Rhipicephalus and Hyalomma. Since the discovery of Dermacentor reticulatus ticks in 2007 and the first reported autochthonous cases in the South of the Netherlands in 2012, no data on the (sero)prevalence of EP in horses in the Netherlands have been reported and it remains unclear whether B. caballi and T. equi have been able to establish themselves in the Netherlands. This study aims to give an update on the current status of EP in horses in the Netherlands using data from serological tests performed in the context of export and screening of 12,881 horses from 2015 through 2020. Horses were categorized as "Dutch," "Foreign," or "Unknown" based on microchip number. The overall seroprevalence of EP in Dutch horses was found to be 0.5% (95% exact CI [0.4-0.7]), compared to 1.9% (95% exact CI [1.3-2.6]) in horses in the category "Foreign" and 1.7% (95% exact CI [1.2-2.3]) in horses in the category "Unknown." In addition, the seroprevalence per country in the category "Foreign" ranged from 0% (0.95% exact CI [0-2.8]) for Ireland to 6.0% (0.95% exact CI [3.5-9.3]) for Spain. In light of the reports on the seroprevalence during the outbreak of autochthonous EP reported in 2012 and on seroprevalences of EP in other countries in Northwestern Europe, the seroprevalence of EP in horses exported from the Netherlands is very low. However, the higher seroprevalence of EP in horses from abroad warrants the need for the monitoring of EP, as tick vectors are present in the Netherlands and the import of horses from endemic areas increases the chances of EP becoming more prevalent in the Netherlands.

Monitoring Wind-Borne Particle Matter Entering Poultry Farms via the Air-Inlet: Highly Pathogenic Avian Influenza Virus and Other Pathogens Risk.

Wind-supported transport of particle matter (PM) contaminated with excreta from highly pathogenic avian influenza virus (HPAIv)-infected wild birds may be a HPAIv-introduction pathway, which may explain infections in indoor-housed poultry. The primary objective of our study was therefore to measure the nature and quantity of PM entering poultry houses via air-inlets. The air-inlets of two recently HPAIv-infected poultry farms (a broiler farm and a layer farm) were equipped with mosquito-net collection bags. PM was harvested every 5 days for 25 days. Video-camera monitoring registered wild bird visits. PM was tested for avian influenza viruses (AIV), Campylobacter and Salmonella with PCR. Insects, predominantly mosquitoes, were tested for AIV, West Nile, Usutu and Schmallenberg virus. A considerable number of mosquitoes and small PM amounts entered the air-inlets, mostly cobweb and plant material, but no wild bird feathers. Substantial variation in PM entering between air-inlets existed. In stormy periods, significantly larger PM amounts may enter wind-directed air-inlets. PM samples were AIV and Salmonella negative and insect samples were negative for all viruses and bacteria, but several broiler and layer farm PM samples tested Campylobacter positive. Regular wild (water) bird visits were observed near to the poultry houses. Air-borne PM and insects-potentially contaminated with HPAIv or other pathogens-can enter poultry air-inlets. Implementation of measures limiting this potential introduction route are recommended.

Qualitative Evaluation of Causes for Routine Salmonella Monitoring False-Positive Test Results in Dutch Poultry Breeding Flocks.

The Salmonella monitoring program, as outlined in the EU Commission regulation 200/2010, asks for repeated sampling in order to ascertain progress in achievement of the EU target. According to Article 2.2.2.2.c of this regulation, the competent authority may decide to do a resample and retest when it has reasons to question the results of initial testing. In the Netherlands, the competent authorities have been resampling and retesting all initial positive samplings for several years because of doubts about false positive initial test results. An analysis of population data in the period 2015-2019 indicates that 48% of initial samplings at the farm were classified as false positive after resampling and retesting by the competent authorities. A qualitative analysis, assessing factors that could be associated with the occurrence of false positives, indicates that cross-contamination during the sampling process by the poultry farmer is probably the most likely source. Cross-contamination of samples during transport from the farm to the laboratory and/or cross-contamination at the laboratory are also considered possible sources. Given the slightly non-optimal system-specificity of the Salmonella monitoring program, there is good reason to make, or consider, standard resampling and retesting of initial positive results by the competent veterinary authorities possible within the EU.

Challenging the "gold standard" of colony-forming units - Validation of a multiplex real-time PCR for quantification of viable Campylobacter spp. in meat rinses.

Campylobacter jejuni is the leading bacterial food-borne pathogen in Europe. Despite the accepted limits of cultural detection of the fastidious bacterium, the "gold standard" in food microbiology is still the determination of colony-forming units (CFU). As an alternative, a live/dead differentiating qPCR has been established, using propidium monoazide (PMA) as DNA-intercalating crosslink agent for inactivating DNA from dead, membrane-compromised cells. The PMA treatment was combined with the addition of an internal sample process control (ISPC), i.e. a known number of dead C. sputorum cells to the samples. The ISPC enables i), monitoring the effective reduction of dead cell signal by the light-activated DNA-intercalating dye PMA, and ii), compensation for potential DNA losses during processing. Here, we optimized the method for routine application and performed a full validation of the method according to ISO 16140-2:2016(E) for the quantification of live thermophilic Campylobacter spp. in meat rinses against the classical enumeration method ISO 10272-2:2017. In order to render the method applicable and cost-effective for practical application, the ISPC was lyophilized to be distributable to routine laboratories. In addition, a triplex qPCR was established to simultaneously quantify thermophilic Campylobacter, the ISPC and an internal amplification control (IAC). Its performance was statistically similar to the two duplex qPCRs up to a contamination level of 4.7 log10Campylobacter per ml of meat rinse. The limit of quantification (LOQ) of the alternative method was around 20 genomic equivalents per PCR reaction, i.e. 2.3 log10 live Campylobacter per ml of sample. The alternative method passed a relative trueness study, confirming the robustness against different meat rinses, and displayed sufficient accuracy within the limits set in ISO 16140-2:2016(E). Finally, the method was validated in an interlaboratory ring trial, confirming that the alternative method was fit for purpose with a tendency of improved repeatability and reproducibility compared to the reference method for CFU determination. Campylobacter served as a model organism, challenging CFU as "gold standard" and could help in guidance to the general acceptance of live/dead differentiating qPCR methods for the detection of food-borne pathogens.

Complete genome for Actinobacillus pleuropneumoniae serovar 8 reference strain 405: comparative analysis with draft genomes for different laboratory stock cultures indicates little genetic variation.

We report here the complete genome sequence of the widely studied Actinobacillus pleuropneumoniae serovar 8 reference strain 405, generated using the Pacific Biosciences (PacBio) RS II platform. Furthermore, we compared draft sequences generated by Illumina sequencing of six stocks of this strain, including the same original stock used to generate the PacBio sequence, held in different countries and found little genetic variation, with only three SNPs identified, all within the degS gene. However, sequences of two small plasmids, pARD3079 and p405tetH, detected by Illumina sequencing of the draft genomes were not identified in the PacBio sequence of the reference strain.