Richard Vipond

Molecular Epidemiologic Investigation of an Anthrax Outbreak among Heroin Users, Europe

Heroin may have been accidentally contaminated by an animal-derived source along a major drug trafficking route.

Genome sequencing shows that European isolates of Francisella tularensis subspecies tularensis are almost identical to US laboratory strain Schu S4.

Background Francisella tularensis causes tularaemia, a life-threatening zoonosis, and has potential as a biowarfare agent. F. tularensis subsp. tularensis, which causes the most severe form of tularaemia, is usually confined to North America. However, a handful of isolates from this subspecies was obtained in the 1980s from ticks and mites from Slovakia and Austria. Our aim was to uncover the origins of these enigmatic European isolates. Methodology/Principal Findings We determined the complete genome sequence of FSC198, a European isolate of F. tularensis subsp. tularensis, by whole-genome shotgun sequencing and compared it to that of the North American laboratory strain Schu S4. Apparent differences between the two genomes were resolved by re-sequencing discrepant loci in both strains. We found that the genome of FSC198 is almost identical to that of Schu S4, with only eight SNPs and three VNTR differences between the two sequences. Sequencing of these loci in two other European isolates of F. tularensis subsp. tularensis confirmed that all three European isolates are also closely related to, but distinct from Schu S4. Conclusions/Significance The data presented here suggest that the Schu S4 laboratory strain is the most likely source of the European isolates of F. tularensis subsp. tularensis and indicate that anthropogenic activities, such as movement of strains or animal vectors, account for the presence of these isolates in Europe. Given the highly pathogenic nature of this subspecies, the possibility that it has become established wild in the heartland of Europe carries significant public health implications.

Description of a Coxiella burnetii abortion outbreak in a dairy goat herd, and associated serology, PCR and genotyping results

Using PCR on aborted foetal material, Coxiella burnetii infection was confirmed as the cause of abortions in a dairy goat herd with over 1000 adults. Ninety-five (22%) abortions and 355 normal births were recorded from 440 goats over 2 months. The herd was sampled three times in 6 months to look at the within-herd seroprevalence, with the 1st visit done 24 days after the last recorded abortion. The true seroprevalence in the herd was 79.2%, 66.5% and 45.7% on each of these visits, but introduction of a group of young goats prior to the 3d visit influenced these results. Using PCR, widespread environmental contamination was demonstrated in surface dust, bedding, muck heaps, milk, bird droppings and drinking water in the goat shed. MST and MLVA analysis showed the C. burnetii from this outbreak to be of a genotype previously observed in the UK and different from those recorded in the Netherlands outbreak of 2007-2011.

Interlaboratory Comparison of Real-Time Polymerase Chain Reaction Methods to Detect Coxiella Burnetii, The Causative Agent of Q Fever

The bacterium Coxiella burnetii, which has a wide host range, causes Q fever. Infection with C. burnetii can cause abortions, stillbirth, and the delivery of weak offspring in ruminants. Coxiella burnetii infection is zoonotic, and in human beings it can cause chronic, potentially fatal disease. Real-time polymerase chain reaction (PCR) is increasingly being used to detect the organism and to aid in diagnosis both in human and animal cases. Many different real-time PCR methods, which target different genes, have been described. To assess the comparability of the C. burnetii real-time PCR assays in use in different European laboratories, a panel of nucleic acid extracts was dispatched to 7 separate testing centers. The testing centers included laboratories from both human and animal health agencies. Each laboratory tested the samples using their in-house real-time PCR methods. The results of this comparison show that the most common target gene for real-time PCR assays is the IS1111 repeat element that is present in multiple copies in the C. burnetii genome. Many laboratories also use additional real-time PCR tests that target single-copy genes. The results of the current study demonstrate that the assays in use in the different laboratories are comparable, with general agreement of results for the panel of samples.

Evaluation of the Biofire FilmArray BioThreat-E Test (v2.5) for Rapid Identification of Ebola Virus Disease in Heat-Treated Blood Samples Obtained in Sierra Leone and the United Kingdom

ABSTRACT Rapid Ebola virus (EBOV) detection is crucial for appropriate patient management and care. The performance of the FilmArray BioThreat-E test (v2.5) using whole-blood samples was evaluated in Sierra Leone and the United Kingdom and was compared with results generated by a real-time Ebola Zaire PCR reference method. Samples were tested in diagnostic laboratories upon availability, included successive samples from individual patients, and were heat treated to facilitate EBOV inactivation prior to PCR. The BioThreat-E test had a sensitivity of 84% (confidence interval [CI], 64% to 95%) and a specificity of 89% (CI, 73% to 97%) in Sierra Leone (n = 60; 44 patients) and a sensitivity of 75% (CI, 19% to 99%) and a specificity of 100% (CI, 97% to 100%) in the United Kingdom (n = 108; 70 patients) compared to the reference real-time PCR. Statistical analysis (Fishers exact test) indicated there was no significant difference between the methods at the 99% confidence level in either country. In 9 discrepant results (5 real-time PCR positives and BioThreat-E test negatives and 4 real-time PCR negatives and BioThreat-E test positives), the majority (n = 8) were obtained from samples with an observed or probable low viral load. The FilmArray BioThreat-E test (v2.5) therefore provides an attractive option for laboratories (either in austere field settings or in countries with an advanced technological infrastructure) which do not routinely offer an EBOV diagnostic capability.

Whole Genome Analysis of Injectional Anthrax Identifies Two Disease Clusters Spanning More Than 13 Years

Background Anthrax is a rare disease in humans but elicits great public fear because of its past use as an agent of bioterrorism. Injectional anthrax has been occurring sporadically for more than ten years in heroin consumers across multiple European countries and this outbreak has been difficult to trace back to a source. Methods We took a molecular epidemiological approach in understanding this disease outbreak, including whole genome sequencing of Bacillus anthracis isolates from the anthrax victims. We also screened two large strain repositories for closely related strains to provide context to the outbreak. Findings Analyzing 60 Bacillus anthracis isolates associated with injectional anthrax cases and closely related reference strains, we identified 1071 Single Nucleotide Polymorphisms (SNPs). The synapomorphic SNPs (350) were used to reconstruct phylogenetic relationships, infer likely epidemiological sources and explore the dynamics of evolving pathogen populations. Injectional anthrax genomes separated into two tight clusters: one group was exclusively associated with the 2009–10 outbreak and located primarily in Scotland, whereas the second comprised more recent (2012–13) cases but also a single Norwegian case from 2000. Interpretation Genome-based differentiation of injectional anthrax isolates argues for at least two separate disease events spanning > 12 years. The genomic similarity of the two clusters makes it likely that they are caused by separate contamination events originating from the same geographic region and perhaps the same site of drug manufacturing or processing. Pathogen diversity within single patients challenges assumptions concerning population dynamics of infecting B. anthracis and host defensive barriers for injectional anthrax. Funding This work was supported by the United States Department of Homeland Security grant no. HSHQDC-10-C-00,139 and via a binational cooperative agreement between the United States Government and the Government of Germany. This work was supported by funds from the German Ministry of Defense (Sonderforschungsprojekt 25Z1-S-431,214). Support for sequencing was also obtained from Illumina, Inc. These sources had no role in the data generation or interpretation, and had not role in the manuscript preparation. Panel 1: Research in Context Systematic Review We searched PubMed for any article published before Jun. 17, 2015, with the terms “Bacillus anthracis” and “heroin”, or “injectional anthrax”. Other than our previously published work (Price et al., 2012), we found no other relevant studies on elucidating the global phylogenetic relationships of B. anthracis strains associated with injectional anthrax caused by recreational heroin consumption of spore-contaminated drug. There were, however, publically available genome sequences of two strains involved (Price et al., 2012, Grunow et al., 2013) and the draft genome sequence of Bacillus anthracis UR-1, isolated from a German heroin user (Ruckert et al., 2012) with only limited information on the genotyping of closely related strains (Price et al., 2012, Grunow et al., 2013). Lay Person Interpretation Injectional anthrax has been plaguing heroin drug users across Europe for more than 10 years. In order to better understand this outbreak, we assessed genomic relationships of all available injectional anthrax strains from four countries spanning a > 12 year period. Very few differences were identified using genome-based analysis, but these differentiated the isolates into two distinct clusters. This strongly supports a hypothesis of at least two separate anthrax spore contamination events perhaps during the drug production processes. Identification of two events would not have been possible from standard epidemiological analysis. These comprehensive data will be invaluable for classifying future injectional anthrax isolates and for future geographic attribution.

Genome sequence analysis of Ebola virus in clinical samples from three British healthcare workers, August 2014 to March 2015.

We determined complete viral genome sequences from three British healthcare workers infected with Ebola virus (EBOV) in Sierra Leone, directly from clinical samples. These sequences closely resemble those previously observed in the current Ebola virus disease outbreak in West Africa, with glycoprotein and polymerase genes showing the most sequence variation. Our data indicate that current PCR diagnostic assays remain suitable for detection of EBOV in this epidemic and provide confidence for their continued use in diagnosis.

Temporal Expression of Peripheral Blood Leukocyte Biomarkers in a Macaca fascicularis Infection Model of Tuberculosis; Comparison with Human Datasets and Analysis with Parametric/Non-parametric Tools for Improved Diagnostic Biomarker Identification

A temporal study of gene expression in peripheral blood leukocytes (PBLs) from a Mycobacterium tuberculosis primary, pulmonary challenge model Macaca fascicularis has been conducted. PBL samples were taken prior to challenge and at one, two, four and six weeks post-challenge and labelled, purified RNAs hybridised to Operon Human Genome AROS V4.0 slides. Data analyses revealed a large number of differentially regulated gene entities, which exhibited temporal profiles of expression across the time course study. Further data refinements identified groups of key markers showing group-specific expression patterns, with a substantial reprogramming event evident at the four to six week interval. Selected statistically-significant gene entities from this study and other immune and apoptotic markers were validated using qPCR, which confirmed many of the results obtained using microarray hybridisation. These showed evidence of a step-change in gene expression from an ‘early’ FOS-associated response, to a ‘late’ predominantly type I interferon-driven response, with coincident reduction of expression of other markers. Loss of T-cell-associate marker expression was observed in responsive animals, with concordant elevation of markers which may be associated with a myeloid suppressor cell phenotype e.g. CD163. The animals in the study were of different lineages and these Chinese and Mauritian cynomolgous macaque lines showed clear evidence of differing susceptibilities to Tuberculosis challenge. We determined a number of key differences in response profiles between the groups, particularly in expression of T-cell and apoptotic makers, amongst others. These have provided interesting insights into innate susceptibility related to different host `phenotypes. Using a combination of parametric and non-parametric artificial neural network analyses we have identified key genes and regulatory pathways which may be important in early and adaptive responses to TB. Using comparisons between data outputs of each analytical pipeline and comparisons with previously published Human TB datasets, we have delineated a subset of gene entities which may be of use for biomarker diagnostic test development.

Development and Operation of Ebola Diagnostic Laboratories Led By Public Health England in Sierra Leone during the West African Ebola Outbreak 2013-2015

The Ebola outbreak 2013-2015 created an urgent need for humanitarian response. Public Health England (PHE), in partnership with the Ministry of Defence (MoD) and Defence Science and Technology Laboratory (DSTL), were tasked by the UK Government (through the Department for International Development (DfID)) to provide Ebola Virus Disease (EVD) diagnostic laboratories. These diagnostic laboratories supported the Ebola Treatment Units (ETU) being established in Sierra Leone. PHE operated arguably the largest diagnostic facilities in Sierra Leone by one unilateral donor: operating 3 laboratories (co-located with ETUs), each processing up to 200 samples a day. During the time of operation (October 2014 to December 2015) over 400 civilian UK staff on rotation were deployed in these laboratories, and between them processed greater than 40,000 samples (~6% positivity rate). Here we summarise the laboratory set-up, design rationale, scope and processes deployed. This information can inform planning in response to future outbreaks. Correspondence to: Daniel Bailey, Public Health England, Porton Down, Salisbury SP4 0JG, UK, Tel: +44 (0)1980 619913; E-mail: Daniel.Bailey@phe.gov.uk

Whole genome sequencing investigation of animal-skin-drum-associated UK anthrax cases reveals evidence of mixed populations and relatedness to a US case.

There have been two anthrax cases affecting people that played and/or made animal-skin drums in the UK during the last 10 years, with single fatal occurrences in Scotland in 2006 and London in 2008. Investigations by the Health Protection Agency (now Public Health England) employing multi-locus-variable number tandem repeat analysis had previously linked the clinical cases to spores associated with animal skins and drums the patients had been in contact with. In this study, whole-genome sequencing of 23 Bacillus anthracis isolates harvested during the investigations was performed. High-quality draft assemblies of these genomes provided greater characterization of the B. anthracis strains present and placed them all upon a new branch of the global phylogeny. Although closely related, the clinical isolates from the two events, and another isolated from a drum-skin-associated case in New York in 2006, were distinct from each other. Multiple distinct genotypes were found during both investigations, implying either multiple contamination events or a single heterogeneous contamination. One environmental isolate from the Scottish incident was more closely related to London isolates than to the other Scottish isolates. As B. anthracis of this subgroup was present at both geographically and temporally distinct events, it may be more widespread at the source of contamination. All isolates were distinct from currently characterized West African strains, despite this being the likely origin of the drums and hides, therefore adding to our knowledge of B. anthracis diversity in the region.