Development and Application of Full Genome Sequencing to Support Epidemiological Investigations During FMD Outbreaks
Total Page:16
File Type:pdf, Size:1020Kb
© Daniel Dalet / d-maps.com 500 km 300 mi O/BUL/1/2010 (Burgas 30/12/2010) O/TUR/926/2010* (Bursa 26/07/2010) Bulgaria O/TUR/1086/2010* (Antalya 16/08/2010) O/TUR/840/2010* (Agri 15/07/2010) Gümüşhane O/TUR/1094/2010* (Giresun 11/08/2010) Kastamonu O/TUR/868/2010* (Eskisehir 20/07/2010) O/TUR/1003/2010* (Sivas 10/08/2010) Bursa Eskişehir Sivas Ağrı O/TUR/18/2010 (Gumushane 09/07/2010) Turkey O/TUR/883/2010* (Kastamonu 23/07/2010) Antalya Iran O/TUR/35/2010 (Erzincan 13/08/2010) O/TUR/153/2010* (Gaziantep) O/TUR/154/2010* (Gaziantep) O/IRN/94/2010 (West Azerbaijan 15/04/2010) O/IRN/92/2010 (West Azerbaijan 14/04/2010) 100 O/AFG/59/2010 (Samangan 14/05/2010) 97 O/PAK/36/2010 (Gilgit-Baltistan 29/07/2010) O/TUR/36/2010 (Kocaeli 13/08/2010) O/UKG/12/2001 (AJ311724) O1/Manisa/TUR/69 (AJ251477) 0.01 Development and application of full genome sequencing to support epidemiological investigations during FMD outbreaks Begoña Valdazo-Gonzalez, Nick J. Knowles, Donald P. King Molecular Characterisation and Diagnostics Group FMD programme Institute for Animal Health GU24 0NF Ash Road, Pirbright, Surrey UNITED KINGDOM Summary and Conclusion Aim • To develop FG sequencing protocols for FMDV • To apply FG sequencing data to discover the transmission pathways of FMDV in 2011 Bulgarian outbreak M&M • 7 epithelium from Bulgaria,5 from Middle East • Full genome sequencing based on 21 overlapping fragments • Statistical parsimony methods (TCS) Preliminary results • Same ancestor for FMDVs collected from: • Kosti; • Rezovo; • Cluster of outbreaks in April • Closer sequenced FMDV: • Wild boar near Kosti, 4nt changes (minimun) • Bursa, Turkey, 26/07/10, 31 nt changes (minimum) Discussion and conclusions • Independent spread in Bulgarian livestock and wildlife from single introduction (from Turkey, most likely) • Cluster of outbreaks in April are related • Existence of intermediary hosts (wild boar most likely) • Careful interpretation of VP1 results FMDV: VP1 molecular characterization * putative functions Membrane-binding Genome-linked 5’UTR Carboxy-terminal (VPg) 3’UTR 2x VP2 Protease VP3 Capsid self-cleaving NTP binding* Protease Polymerase 5x VP1 5x VPG VP1 L 1A 1B 1C 1D 2B 2C 3A 3B 3C 3D VP4 VP2 VP3 VP1 2A VP2 VP3 VP1 AAA (n) 3x VP1 VP3 VP2 Poly(C) VP2 2x VP2 VP3 VP3 2x Primary cleavages L 2A 3C VP3 VP2 VP1 VP1 Secondary cleavages 1B/RNA? 3C 3C 3C 3C 3C 3C 5x Kilobases 0 1 2 3 4 5 6 7 8 VP1 region 640 nt • FMDV is a rapidly evolving RNA virus • VP1 phylogenetic analysis widely used for broad scale FMDV strain characterisation • Useful for regional and country-level epidemiology (transboundary movements) © Daniel Dalet / d-maps.com 500 km 300 mi © Daniel Dalet / d-maps.com 500 km 300 mi O/BUL/1/2010 (Burgas 30/12/2010) O/TUR/926/2010* (Bursa 26/07/2010) Bulgaria O/TUR/1086/2010* (Antalya 16/08/2010) O/TUR/840/2010* (Agri 15/07/2010) Gümüşhane O/TUR/1094/2010* (Giresun 11/08/2010) Kastamonu O/TUR/868/2010* (Eskisehir 20/07/2010) O/TUR/1003/2010* (Sivas 10/08/2010) Bursa Eskişehir Sivas Ağrı O/TUR/18/2010 (Gumushane 09/07/2010) Turkey O/TUR/883/2010* (Kastamonu 23/07/2010) Antalya Iran O/TUR/35/2010 (Erzincan 13/08/2010) O/TUR/153/2010* (Gaziantep) O/TUR/154/2010* (Gaziantep) O/IRN/94/2010 (West Azerbaijan 15/04/2010) Bulgaria 2011 outbreak: O/IRN/92/2010 (West Azerbaijan 14/04/2010) 100 O/AFG/59/2010 (Samangan 14/05/2010) 97 O/PAK/36/2010 (Gilgit-Baltistan 29/07/2010) O/TUR/36/2010 (Kocaeli 13/08/2010) O/UKG/12/2001 (AJ311724) O1/Manisa/TUR/69 (AJ251477) 0.01 VP1 molecular characterization Bulgaria Kastamonu Gümüşhane Bursa Eskişehir Sivas Ağrı Turkey Iran Antalya • FMD viruses from Bulgaria characterized as belonging to ME-SA/PanAsia-2ANT-10 lineage • Closely related PanAsia-2 viruses have been recovered from the Middle East during 2010 • VP1 sequences cannot discrimminate between closely related FMDV field strains Use of complete genome sequences for high- resolution tracing: • Uncertainty about precise source of infection for the majority of IPs during field outbreaks since FMDV is a rapidly evolving virus: • Nucleotide changes accrue linearly with time and are inherited • Can we use full genome sequence data to trace the spread of FMDV during an outbreak? Background: 2007 UK outbreak IAH2 IP6b AY593815 IP1b(2) IP3c MAH IAH1 IP1b(1) IP3b IP8 IP4b IP2b IP7 Sampled virus IP2c IP5 Putative ancestor virus IP5 • Detected by sero-surveillance M4 Nt change • After IP3 and IP4 Aa change • Seropositive cattle and sheep WINDSOR HEATHROW His to Arg • No acute clinical signs Asp to Gly • Evidence of healed lesions 8 8 6 6 EGHAM 3b3b nd 3c • 2 phase of outbreaks (IP3 – IP8) 7 7 4 4 3c X 5 5 X shares all the unique changes X X common to 1st phase M3 Therefore outbreaks are linked WOKING X M25 X and not due to independent Pirbright 2b 2b sources 1b1b X 2c GUILDFORD • IP5 (farm with FMD serology ALDERSHOT 10 km positive cattle and sheep) bridges gap GODALMING between two phases of the outbreak X 1c KEY FMD confirmed Preclinical (lab only) X No evidence of infection Goals of complete genome sequence analysis [1] to identify the most likely source of the FMD outbreaks in Bulgaria [2] reconstruct the relationship between outbreaks within Bulgaria (build a transmission tree) [3] assess whether there are gaps in the sampling of disease outbreaks (evidence for undisclosed infection) Full genome sequencing: material & methods Original clinical sample One operator: 48 hours • Epithelium Two operators: 24 hours - suspension BULGARIA – One per herd, higher CT value RNA extraction Wild boar – BUL1/2011 • RNeasy Mini Kit (QIAGEN) Starting point: Firat-Saraç, M. SAP, Turkey FAO/EuFMD (Project PR41764) 12LP1 – Kosti DTU Veterinary. National 12LPN3 – Rezovo Veterinary Institute Reverse Transcription BUL11/11 - Kirovo • Oligo-dT primer (Rev 6) (Ryan et al. 2008) BUL20/2011 – Golyamo Bukovo Negative controls~700b BUL26/2011 -Granichar cDNA clean up BUL30/2011 – Fakia ~700b PCR MIDDLE EAST – VP1 5UTR+P1 NSP+PoliA Backup ISR2/2011 cc DNA clean up 4 epithelium from Turkey Improved protocol Cycle sequencing reaction Ethanol Precipitation Negative controls ABI PRISM 3730 DNA Analyzer ~700b Data analysis:Statistical parsimony methods (TCS) 2-4 coverage /site 5UTR+P1 NSP+PoliA Backup Origins of the outbreaks in Bulgaria: Preliminary results: 7087 nt (most L-fragment ) ISR/02/2011 TUR/926/2010 2/2011 (Bursa) 26/07/10 10 nt changes 57 nt changes 30/12/2010 21 nt changes BUL/1/2010 (wild boar) 3 nt changes 21 nt changes 9 nt changes 12 nt changes 41 nt changes TUR/18/2010 (Gümüşhane) TUR/840/2010 09/07/2010 (Agri) TUR/36/2010 15/07/10 (Kocaeli) 13/08/2010 Putative common ancestor of the Turkish and Israeli viruses Putative common ancestor of Bulgarian viruses Relationship between Bulgarian field strains: TCS analysis - FG (First phase/wave) Kosti village 11 nt changes 12LPN1 (IP 1) 4 nt changes TUR/18/2010 BUL/1/2010 35 nt changes (Gümüşhane) (wild boar) 09/07/2010 30/12/2010 12LPN3 (IP 4) 14 nt changes Rezovo village putative common ancestor of the Bulgarian viruses All nucleotide substitutions are unique includes data from DTU, Lindholm Interpretation of data from first wave/phase • The FMD outbreaks in Kosti and Rezovo villages are unlikely to be directly linked. • Both have a common ancestor close to the virus from wild boar. • The long branches between the putative common ancestor and each of the two outbreaks examined (IP1 and IP4) suggest possible intermediate hosts (either wildlife or domesticated animals). Full genome sequencing: TCS analysis Preliminary results (26/4/11): 7087 nt (most L-fragment ) 30/12/2010 “TREE-LIKE” ≠ Linear pattern (UK2007) BUL/1/2010 02/2011 Wild boar 8 nt changes 12LPN1 Kosti 4 nt changes 02/2011 14 nt changes 12LPN3 Rezovo 6 nt changes BUL/11/11 Kirovo changes (Bursa) 19/03/2011 26/07/2010 4 nt changes 31 nt 31 TUR/926/2010 28/03/2011 13 nt changes BUL/26/11 Granichar 28/03/2011 11 nt changes BUL/20/11 BUL/30/11 Golyamo Fakia Bukovo 01/04/2011 Putative common ancestor of Bulgarian wild board and first phase of the outbreaks Putative common ancestor of the second phase of the Bulgarian outbreaks Full genome sequencing: TCS analisys Preliminary results: 7087 nt (most L-fragment ) Secon Phase/wave March -April Goliamo Black Sea Bukovo Fakia Granichar Kirovo Kosti Wild boar First phase/wave January/February Rezovo Turkey Future work • Keep sequencing Gaps TUR8/2011 Momina Tsarkva • Link genetic and epidemiological data • Further analysis: Comparison with previous studies: 2001 and 2007 UK outbreaks Molecular clock: BEAST Acknowledgements Institute for Animal Health Diagnosis of Vesicular disease group:Nigel Ferris and Geoff Hutchings Molecular diagnostic and diagnostic group:Jemma Wadsworth, Miki Madi and Valerie Mioulet SAP Foot-and-Mouth Disease Institute, Ulus, Ankara, TURKEY Müge Firat-Saraç, Ünal Parlak and Fuat Ozyoruk National Diagnostic and Research Veterinary Medical Institute Lilyana Polyhronova and Georgi Kirilov Georgiev National Veterinary Institute (DTU) Graham Belsham Food and Agriculture Organization of the United Nations/European Commission for the Control of Foot-and- Mouth Disease (Project PR41764) Department for Environment, Food and Rural Affairs (SE2938) Biotechnology and Biological Sciences Research Council Item 11 REPORT OF THE SESSION OF THE RESEARCH GROUP OF THE STANDING TECHNICAL COMMITTEE OF THE EUROPEAN COMMISSION FOR THE CONTROL OF FOOT-AND-MOUTH DISEASE (EUFMD) HELD AT Vienna (Austria) 29 September-1 OCTOBER 2010 (OPEN SESSION) 1 OCTOBER 2010 (CLOSED SESSION) FOOD AND AGRICULTURE ORGANIZATION OF THE UNITED NATIONS ROME, 2010 Table of Contents Acknowledgements ......................................................................................................................................