2019-01-03

Introduction Discoveryofunconventionalviruses associatedw ithbovinerespiratory diseases– w iththehelpofhigh Non-culturab le ordifficult- Nove lorm utate d organism throughputsequencing to-culture organism

Yanyun Huang Prairie Diagnostic Services (PDS) Inc. High Throughput Sequencing (HTS)

h ttp s://www.google.ca/search ?q=ce ll+culture &rlz=1C1CHBF_e nCA758CA758&source 1 2

Introduction Introduction

First Generation Sequencing Second Generation Sequencing Why High Throughput Sequencing

Virusdiscove ry Mutation discove ry and p rofiling In 2011, Sch m allenb e rg (SBV) waside ntified and DNA Prote inInte ractions nam e d using High th rough p utseque ncing (HTS) from cattle and she e p inGe rm any and th e Ne th e rlands. Transcriptom e Analysis

Ap p lication of HTS inve te rinary m e dicine islimite d inCanada. Advantage s: 1. Se que nce a large num b e rof DNA Limitation: m oleculesinp arallel; 1. Sm allam ountsof DNA 2. Do not ne e d p riorknowledge of wh at p roce ssed p e runittime ; y ou’re seque ncing. 2. h igh costp e rb ase: firsth um an 3. Ch e ap costp e rb ase: a single h um an ge nom e $100 m illion ge nom e $1000 h ttp s://www.b iote ch articles.com /Ge ne tics-Article/Sange r-Se que ncing-Ch ain-te rmination-m e th od-of-DNA-Se que ncing-3937.h tm l 3 4

Introduction Introduction

$1 b illion, annually. Controlled with $30 p e rm orbidcalf. Th e re are sub stantialstudiesto supp ort th at viralinfe ction gre atly incre ased antimicrob ialsand susce p tibility of b e e f cattle to b acte rialinfe ction. 70-80% of th e m orbidity viralvaccination (USDA, 2013)

First60 day s Th e re are alsosom e diagnosticcases,wh e re lesionssugge stvirale tiology, h owe ve r, conve ntionalvirusescannot b e de te cte d b y PCR assay s. Use of antimicrob ials prom ote sde ve lopm e nt of antimicrob ial- re sistance (AMR). Conve ntionalviruses= traditionally associate d with BRD, include d inBRD vaccine sand availab le diagnosticassay s h ttp ://th re e rive rsve tgroup.co.uk/ve te rinary-service s/farm-animals/info-for-farme rs/b ovine -re spiratory-disease-b rd-p ne umonia/ 5 6

1 2019-01-03

Hypothesis

Sample information 240 sam p les(120 p aire d nasalswab sand trach e alwashe s)kindly p rovide d b y Dr.Tre vorW . Alex ande r (Le th b ridge Re search Ce nte r,Agriculture and Agri-Food Canada) Healthy control cattle BRD cattle Th e re are unconve ntional virusesassociate d with BRD

1. No visualBRD signs 1. One visualBRD sign, at least 2. Re ctalte m p e rature < 40 °C 2. Re ctal te m p e rature ≥ 40 °C 3. No tre atm e nt againstoth e rdisease 3. Ab norm allung sounds 4. No ab norm allung sounds 4. No tre atm e nt againstoth e rdiseases 5. If sick with in60 day s,re m ove d h ttp s://www.google.ca/search ?rlz=1C1CHBF_e nCA758CA758&tb m =isch &q=h ap p y +ste e r+cartoon&ch ips=q:h ap p y +ste e r+cartoon 7 8

Materials and Methods Materials and Methods NexteraXT library preparation and sequencing Sample preparation

5’ 3’ RNA Trach e aland Ex traction of viral 0.2µg Input DNA Nasalswab nucleicacids Re verse transcriptase + random p rime rs+ dNTPs First-strand synthesis 1 Bioanalyzer 5’ 3’ Ce ntrifuge DNA-RNA hybrid 6 3’ 5’ DNase/ RNase 2 5 Adde d RNase H and Se que nase Library Normalization and Pooling Tagmentation of Input DNA 7 PCR Amplification Hostce lls 3’ 5’ Second-strand synthesis with Index Primers 3’ 5’ 4 Bacte ria Polym e rase and p rime rs with out6 Ns

Amplification h ttp s://www.ab m good.com /m arke ting/knowledge _b ase/ne x t_ge ne ra PCR Clean-up Sample loading and Data analysis tion_seque ncing_e x p e rime ntal_de sign.p h p 9 10

Materials and Methods

Raw data Assem b led th e re st re ads

Ob taine d h igh quality re ads Map p e d th e (Trimm om atic) assem b led seque nce s to Re fSe q datab ase Map p e d to cow ge nom e to re m ove h ost-de rive d seque nce s

h ttp s://e n.wikipe dia.org/wiki/Illum ina_dy e _seque ncing 11 12

2 2019-01-03

Data Analysis

Sample Total No. of Reads No. of Reads unmapped to cattle Percentage

Nasal swab 33,813,931 9,618,962 28.45%

Tracheal wash 48,936,263 24,026,779 49.10%

Total 82,750,194 33,645,741 40.66%

Virus 1,812,011 2.19%

13 14

Table 1. Viruses identified by metagenomics

Largest contig Total No. of Largest contig Table 2. Virus identification rates and association between presence of viruses and the probability of an Virus Family Closest relative length (bp) sequence reads amino acid identity (%) animal showing symptoms of BRD IDV 1,587 17,297 99 (PB2) Influenza D virus No. of positive cattle/total BRD ICV Orthomyxoviridae 1,010 307 100 (PB1) Influenza C virus or control cattle (% positive) BRBV Picornaviridae 2,431 38,648 99 (polyprotein) Bovine rhinitis B virus Virus Odds ratio 95% CI for odds ratio P value BRAV Picornaviridae 1,296 1,022 100 (polyprotein) Bovine rhinitis A virus Symptomatic Asymptomatic EVE Picornaviridae 3,186 20,124 98 (polyprotein) Enterovirus E BRSV 1,169 121,005 98 (RdRp) Bovine orthopneumovirus IDV 13 (22) 3 (5) 6.145 1.435 - 26.310 0.015* BPIV3 Paramyxoviridae 279 49 99 (M) Bovine respirovirus 3 BRBV 16 (28) 6 (10) 3.836 1.245 - 11.821 0.020* BCV 7,513 197,921 99 (ORF1ab) Bovine coronavirus BRSV 10 (17) 1 (2) 13.422 1.454 - 123.885 0.022* BNV Coronaviridae 4,782 86,392 99 (PP1a/b) Bovine nidovirus TCH5 BPV2 1,149 1,427 93 (nonstructural protein) Bovine parvovirus-2 BCV 11 (19) 2 (3) 7.392 1.354 - 40.346 0.021* BAAV Parvoviridae 1,096 1,002 99 (Cap) Bovine adeno-associated virus BRAV 7 (12) 2 (3) 5.659 0.982 - 32.602 0.052 UTPV1 Parvoviridae 4,375 1,061,037 98 (NS1, VP1 and VP2) Bovine hokovirus 2 BPV2 7 (12) 3 (5) 3.289 0.682 - 15.865 0.137 UBPV6 Parvoviridae 4,518 263,902 99 (nonstructural protein) Ungulate bocaparvovirus 6 BNV 4 (7) 23 (40) 0.078 0.021 - 0.288 0.000ǂ BVDV1 602 12 97 (NS5b) Bovine viral diarrhea virus 1 ICV 0 (0) 6 (10) ------0.967 HCV Flaviviridae 528 11 97 (core protein) Bovine hepacivirus BAdV3 366 14 99 (284R) Bovine adenovirus 3 BAV 5 (9) 1 (2) 4.485 0.459 - 43.798 0.195 BAV Astroviridae 1,220 884 98 (NSP1ab) Bovine UTPV1 26 (45) 16 (28) 1.878 0.812 - 4.348 0.140 ssCDV 676 158 91 (Rep) Gila monster-associated gemycircularvirus UBPV6 8 (14) 20 (34) 0.296 0.108 - 0.814 0.019ǂ WUPyV 731 337 77 (large T antigen) WU Polyomavirus WUPyV 3 (5) 6 (10) 0.421 0.081 - 2.185 0.300 PBCV 351 288 96 (CVM1) Paramecium bursaria Chlorella virus EVE 6 (10) 2 (3) 4.258 0.704 - 25.740 0.113 HPV 763 174 100 (major capsid protein) Human papillomavirus type 40 BAdV3 1 (2) 2 (3) 0.339 0.025 - 4.619 0.414

15 16

Conclusionandimplications Table 3. Association between presence of at least one of the following four viruses and the probability of an animal showing symptoms of BRD •There are unconventional viruses in BRD samples Total No. of cattle positive •Several viruses are associated with BRD in the current data set Virus containing at least 1 of 4 viruses Odds ratio 95% CI for odds ratio PDOF value P value

Symptomatic Asymptomatic •Further studies needed to understand the roles of these viruses in IDV BRD BRBV 36 10 7.484 2.167-25.849 0.001* 0.002* BRSV •HTS is a powerful technique for virus detection BCV

17 18

3