key: cord-0762702-3im1l78y authors: Ratta, Barkha; Yadav, Brijesh Singh; Pokhriyal, Mayank; Saxena, Meeta; Sharma, Bhaskar title: Microarray Chip Based Identification of a Mixed Infection of Bovine Herpesvirus 1 and Bovine Viral Diarrhea 2 From Indian Cattle date: 2013-09-12 journal: Curr Microbiol DOI: 10.1007/s00284-013-0448-9 sha: e46fb86507e474f0ff51b5332c6a7b1b4a1ecbb5 doc_id: 762702 cord_uid: 3im1l78y Bovine herpesvirus 1 (BHV1) and bovine viral diarrhea virus 2 (BVD2) are endemic in India although no mixed infection with these viruses has been reported from India. We report first mixed infection of these viruses in cattle during routine screening with a microarray chip. 62 of the 69 probes of BHV1 and 42 of the 57 BVD2 probes in the chip gave positive signals for the virus. The virus infections were subsequently confirmed by RT-PCR. We also discuss the implications of these findings. Indian cattle and buffalo populations (200 million and 75 million, respectively) are the largest in the world, and India is emerging as one of the largest milk-producing countries of the world. The health and productivity of Indian cattle and buffalo populations are dependent on the absence or the low scale presence of many endemic diseases which are prevalent elsewhere, particularly bovine viral diarrhea (BVD) and bovine herpesvirus 1 (BHV1). BVD is a disease which is endemic in developed countries. It is only in this century [15, 17, 18] that the virus has been isolated from Indian livestock, though serological evidence of the existence of this virus was reported in 1980s [32] . Both BVD virus genotypes 1 and 2 have been reported [4, [16] [17] [18] [19] 29] . BVDV, a member of the Pestivirus genus of the family Flaviviridae, has been linked to a number of clinical conditions of cattle, viz., diarrhea, abortions, congenital abnormalities, and birth of persistently infected animals [24] . Acute BVDV infections in postnatal life are mostly inapparent [2] . BVDV generally induces severe immunosuppression in affected animals. The immunosuppression is associated with transient leukopenia, neutrophil dysfunction, and other immunological disturbances [5, 12] . The immunosuppressive effect impairs the ability of the animal to clear virus or bacteria from its blood and tissues [21, 28] . Infectious bovine rhinotracheitis/infectious pustular vulvovaginitis caused by BHV1 was first reported in the US in 1950s and in India in 1976. BHV1 is a member of genus Varicellovirus, subfamily Alphaherpesviriniae of the family Herpesviridae. BHV1 causes several diseases in cattle including rhinotracheitis [20] , vaginitis, balanoposthitis, abortion, conjunctivitis, and enteritis [8] . The World Organization for Animal Health [OIE] lists BHV1 as a List B-notifiable disease. OIEs List B diseases are transmissible diseases that are of socioeconomic and/or public health importance and that are significant in international trade. Acute BHV1 infection, regardless of whether there are clinical signs of disease, leads to latent infection. Latently infected cattle show no disease unless the latent infection is reactivated [13] . In India, the virus has been isolated from states of Orissa [23] , Karnataka [22] , and Gujarat [10] . The serological evidence for the existence of BHV1 in Indian domestic and wild livestock [34] indicates that this disease is widespread in India. BHV1 and BVD are the two diseases that are difficult to manage and eradicate because of latent infections. Mixed infections with these two viruses have been reported from different countries [1, 9] . The economic cost of mixed infection has not been determined but is apparently high. In India, the seroprevalence of both BVD and BHV1 in Indian livestock is known. So far, no report of mixed infections of these viruses from Indian livestock has been reported. We report the first mixed infection of BHV1 and BVD in Indian cattle using a microarray chip that we have developed, and discuss the implications of the findings and also the use and utility of microarray chips for identifying mixed and emerging infections. The clinical sample collected from the institute's veterinary polyclinic was obtained from a cow which had aborted fetus (5 months) and was having fever. We have been designing and testing microarray chips for identification of animal viruses. Different versions of chips have been tested and an earlier version of the chip identified Newcastle disease virus from sheep [31] . The chip used here for testing the clinical samples contained 4079 unique probes representing 124 different viruses and 3000 random probes. The probes were designed using e-array (https://earray.chem.agilent.com) or array designer (www. premierbiosoft.com) and were analyzed using BLAST program individually before including the same in the chip. The probe design criteria were (1) length of oligonucleotide probes 50-60. ( 2) The Tm (melting point) of the probes 700°C ± 5. (3) The GC content of the probes 40-60 %. (4) The maximum poly nucleotide tract allowed was six. (5) 90 % homology with target sequence with at least 20 nucleotide contiguous matches at 3 0 end and less than 85 % homology with nontarget sequence. The random probes included in the chip were taken from MDA (microbial detection array) chip [7] . The chip was fabricated commercially (Agilent Technologies, USA) as per our design. Total RNA was extracted using TRIzol Ò Reagent (Invitrogen) according to manufacturer's protocol. The quality of total RNA was checked by Nanodrop spectrophotometer and bioanalyzer. The total RNA was labeled with Cyanine 3-CTP labeling dye using Quick-Amp labeling Kit one-color (Agilent) according to manufacturer's protocol. The purified Cy3-labeled samples were hybridized on Agilent's Microarray Hybridization Chamber using Agilent Microarray Hybridization Chamber Kit (G2534A) according to manufacturer's instruction. The total RNA was amplified by Quantitect Ò Whole transcriptome amplification Kit (QIAGEN, cat no. 207043) according to the manufacturer's protocol. After hybridization, the arrays were scanned with an Agilent scanner. The microarray data from scanned images were extracted with Agilent Feature Extraction software version 9.5.3. The signal intensities, retrieved from raw data text file generated from hybridized image file using feature extraction software, were sorted by genus and species of the virus. Many methodologies have been employed for identifying a virus from microarray data. These include hybridization intensity based on expected and observed intensity [35] or on log odd ratio [7] , number of probes giving positive signal, or even percentage of probes giving positive signal [27] . We had adopted percentage of probes giving positive signals for making virus call. This method of making a virus call is easy and has 50 % as cut-off for making a virus call is stringent others have kept it at about 30 % [27] . We had kept it at 50 % based on our earlier observation that the closely related viruses give high cross reactivity. In the earlier version of the chip, we used to get 60-70 % cross hybridization of probes between BHV1 and BHV5 and CAV1 and CAV2. This has been brought down to about 30-35 % with the chip that was used in this study. The background intensity for identifying positive signal was arrived from intensity of random probes and was kept at 99th percentile of random probe signal intensity [7] . The primer sequences used for amplification of BVD2 and BHV1 genes are listed in Table 1 . The amplification for BHV1 gE gene was carried out in 50 ll reaction volume containing 19 Taq DNA Polymerase buffer, 1.5 mM MgCl 2 , 20 pmol of each primer, 200 lM of dNTPs, 4 % DMSO, 1 ll of amplified transcriptome, 2 U Taq Polymerase (Fermentas), and nuclease-free water to make up volume of 50 ll. After initial denaturation at 95°C for 4 min, the amplification was carried out for 34 cycles each of 95°C-30 s, 61°C-45 s, and 72°C-1 min with final extension of 5 min at 72°C. The amplification for BVD2 was carried out in 50 ll reaction volume containing 19 Taq DNA polymerase buffer, 1.5 mM MgCl 2 , 20 pmol of each primer, 200 lM dNTPs, 1 ll amplified transcriptome, 2 U of Taq DNA Polymerase (Fermentas), and nucleasefree water to make 50 ll. After initial denaturation at 95°C for 4 min, the amplification was carried out for 34 cycles each of 95°C-30 s, 50°C-45 s, and 72°C-1 min with final extension of 10 min at 72°C. For nested PCR of BVD2, the procedure was essentially the same except that the template was replaced by 1 ll of primary PCR amplicon and primers of set B ( Table 1 ). The annealing temperature was kept at 50°C. Of the many samples tested in one sample from cow, 62 of the 69 probes of BHV1 and 42 of the 57 probes of BVD2 gave positive signals ( Table 2) . A virus was called to be present if more than 50 % probes for that virus gave signal above cut-off. According to this criterion, both the viruses were present in this particular sample. Next virus for which we obtained some positive signal was BHV5 (13 of the 37 probes gave positive signal). BHV1 and BHV5 (formerly BHV1.3) are closely related with more than 90 % sequence similarity, and there is always some cross hybridization. We normally get this much cross hybridization with BHV5 when RNA from BHV1-infected cell culture lysate is hybridized with the chip (personal communication). We did not get any further significant signal for other viruses in this sample (Fig. 1) . The average signal intensity of BHV1 probes was consistently high as compared to BVD2 probes average signal intensity (Table 2) suggesting that this was probably not necessarily true, BHV1 was the major infection and BVD2 a minor infection. Whole transcriptome PCR results (Fig. 2 , amplified product 613 bp for BHV-1 and 198 bp and 105 bp for BVD1, and BVD2 genotype, respectively) corroborated the microarray data. The unambiguous result in both the tests and clinical history (aborted fetus and fever) of the animal makes it almost a certainty that there was a mixed infection of BHV1 and BVD2 in the animal from which the sample was drawn. BHV1, BVD, respiratory syncytial virus (RSV), bovine parainfluenza virus 3, bovine adenovirus, and coronavirus are the viruses found which are associated with the bovine respiratory disease complex [6] . In India, cattle movement is frequent; cattle are moved across long distances under highly stressed conditions; shipping fever is known to exist; and, in most cases, it is bacteria, which are isolated from these cases [30] ; the underlying virus cause of [26, 28] , except India. They increase the complexity of the symptoms and economic losses [33] . Among BVD1 and BVD2 genotypes, it is the BVD2 genotype which is more virulent [4] . Among the few clinical samples tested, one sample was positive for BHV1 and BVD2. In India, BVD is still considered as an exotic disease, and handling the virus is confined to a lab with high-security disease though serological evidence [17] indicates otherwise. Large-scale screening of clinical samples may give the actual extent of spread of mixed infections of these viruses. However, based on serological evidences [3, 17, 34] , it can be assumed that mixed infection of these viruses could be present on a much larger extent and mostly go unreported due to want of extensive search. Microarray chips for diagnosis have emerged as highly parallel diagnostic platforms, which that have the ability to detect many known, novel, and emerging pathogenic agents simultaneously. The Virochip developed originally by Palacios is one such chip [25] , microbial detection array (MDA) [7] is even more comprehensive chip and has the capability to screen almost all the viruses currently known to exist. Most of microarray chips developed and tested were for human pathogens though they contains probes for animal viruses also. The Virochip was recently used for identifying pathogens in veterinary clinical samples [14, 21] . We had also tested in-house designed chip for testing veterinary clinical samples [35] . Massive multiplexing ability of microarray as compared to other diagnostic techniques is advantageous in identifying multiple infections which may otherwise go unreported. Disadvantage of microarray chip is its high cost which makes its applicability in routine diagnosis highly improbable. Many lowdensity FDA approved microarrays are now commercially available, e.g., xTag Respiratory Viral Panel (Luminex), Filmarray (Biofire), and TessArray RPM-Flu 3.1 Array (TessArae), ResPlex II v2.0 RVP (Qiagen), they are now used regularly in many clinical microbiology labs. The use of high-density microarrays is limited for want of FDA or other regulatory agencies' approval. Their approval by regulatory agencies required testing against all the pathogens; they are capable of identifying, which is technically not feasible. Current cost of using high-density array is high (200-300 USD). The cost is expected to go down substantially as the methodology gets widely adopted. High-cost microarray can still be used as an adjunct diagnostic method for screening of multiple infections, and in those cases where routine diagnostic methods fail to identify a pathogen from clinical sample. Bovine viral diarrhoea, bovine herpesvirus and parainfluenza-3 virus infection in three cattle herds in Egypt in 2000 The causative agent of BVD: its epidemiology and pathogenesis A serological survey of antibodies against bovine herpesvirus-1 in yak (Poephagus grunniens) in Arunachal Pradesh in India Genetic and antigenic characterization of bovine viral diarrhea virus type 2 isolated from cattle in India Characterization of the major structural proteins of purified bovine viral diarrhea virus Respiratory viruses and bacteria in cattle A microbial detection array (MDA) for viral and bacterial detection Bovine herpesviruses. Part I. Bovine herpesvirus 1 Disease outbreak in calves caused by a mixed infection with infectious bovine rhinotracheitis virus and bovine virus diarrhoea virus Seroprevalence of bovine herpesvirus 1 Use of polymerase chain reaction to simultaneously detect and type bovine viral diarrhoea viruses isolated from clinical specimens Evidence for biotype-specific effects of bovine viral diarrhea virus on biological responses in acutely infected calves Latent bovineherpesvirus 1 infection in calves protected by colostral immunity Identification of a novelcoronavirus from a beluga whale by using a panviral microarray Genetic and antigenic characterization of Indian bovine viral diarrhea virus type 2 isolated from Indian goats (Capra hircus) Genetic typing of bovine viral diarrhoea virus isolates from India Identification of bovine viral diarrhea virus type 1 in Indian buffaloes and their genetic relationship with cattle strains in 5 0 UTR Prevalence of bovine viral diarrhoea virus antibodies among sheep and goats in India Molecular characterization of bovine viral diarrhea virus type 2 isolate originating from a native Indian sheep (Ovies aries) Infectious bovine rhinotracheitis virus infection and infertility in cows, heifers and bulls Pathogenesis of intrauterine infections with bovine viral diarrhea virus Isolation of infectious bovine rhinotracheitis virus in Karnataka Note on serological evidence of viral abortion in cattle in Orissa Ruminant pestiviruses Panmicrobial oligonucleotide array for diagnosis of infectious diseases Effect of bovine viral diarrhea virus infection on the distribute on of infectious bovine rhinotracheitis virus in calves Detection of respiratory viruses and subtype identification of influenza A viruses by GreenChipResp oligonucleotide microarray Association of bovine viral diarrhea virus with multiple viral infections in bovine respiratory disease outbreaks Differentiation of types 1a, 1b and 2 bovine viral diarrhoea virus (BVDV) by PCR Studies on shipping fever of cattle. 3. Comparison of Pasteurella spp. isolated from shipping fever and other infectious processes Isolation of Newcastle disease virus from nonavian host sheep and its implications Seroprevalence of bovine viral diarrhoea virus in India: a survey from 1999-2004 Prevalence of bovine viral diarrhoea virus antibodies in India Seroprevalence of infectious bovine rhinotracheitis in India E-Predict: a computational strategy for species identification based on observed DNA microarray hybridization patterns Respiratory viruses in acute respiratory tract infection in western India Acknowledgments The authors acknowledge and thank ICAR for funding and Dr. Umesh Dimri for providing clinical samples from institute polyclinic.