key: cord-0838472-oujdl459
authors: Lung, O.; Ohene‐Adjei, S.; Buchanan, C.; Joseph, T.; King, R.; Erickson, A.; Detmer, S.; Ambagala, A.
title: Multiplex PCR and Microarray for Detection of Swine Respiratory Pathogens
date: 2015-12-12
journal: Transbound Emerg Dis
DOI: 10.1111/tbed.12449
sha: 1eb731dcdcd8565f3865b47be6b7f722b596a416
doc_id: 838472
cord_uid: oujdl459

Porcine respiratory disease complex (PRDC) is one of the most important health concerns for pig producers and can involve multiple viral and bacterial pathogens. No simple, single‐reaction diagnostic test currently exists for the simultaneous detection of major pathogens commonly associated with PRDC. Furthermore, the detection of most of the bacterial pathogens implicated in PRDC currently requires time‐consuming culture‐based methods that can take several days to obtain results. In this study, a novel prototype automated microarray that integrates and automates all steps of post‐PCR microarray processing for the simultaneous detection and typing of eight bacteria and viruses commonly associated with PRDC is described along with associated multiplex reverse transcriptase PCR. The user‐friendly assay detected and differentiated between four viruses [porcine reproductive and respiratory syndrome virus (PRRSV), influenza A virus, porcine circovirus type 2, porcine respiratory corona virus], four bacteria (Mycoplasma hyopneumoniae, Pasteurella multocida, Salmonella enterica serovar Choleraesuis, Streptococcus suis), and further differentiated between type 1 and type 2 PRRSV as well as toxigenic and non‐toxigenic P. multocida. The assay accurately identified and typed a panel of 34 strains representing the eight targeted pathogens and was negative when tested with 34 relevant and/or closely related non‐target bacterial and viral species. All targets were also identified singly or in combination in a panel of clinical lung samples and/or experimentally inoculated biological material.

The global pig industry produced approximately 963 million pigs and 109 million metric tons of pork in 2011 (http://faostat.fao.org). Respiratory diseases are considered to be one of the main contributors to economic losses in the swine industry (Opriessnig et al., 2011) . The 2006 United States Department of Agriculture (U.S.D.A.) National Animal Health Monitoring System (NAHMS) study of 435 swine production sites with 100 or more pigs from 17 major pork-producing States showed that respiratory problems are the main cause of nursery deaths (53.7%) and grower/finisher pig mortality (60.1%) (United States Department of Agriculture, 2008) . In Canada, 37-78% of pigs going for slaughter have cranioventral bronchopneumonia (Hansen et al., 2010a) . Porcine respiratory disease complex (PRDC) is multifactorial, with both infectious and non-infectious factors contributing to respiratory disease and predominantly seen in pigs between the ages of 3 and 6 months (Opriessnig et al., 2011) . The interaction of viral and bacterial pathogens, environmental factors, pig-specific factors and management conditions all contribute to the development and impact the severity of PRDC (Opriessnig et al., 2011) .

The type of pathogens involved in PRDC is specific to the regions and countries where production occurs (Opriessnig et al., 2011) . However, viruses most commonly associated with PRDC include porcine reproductive and respiratory syndrome virus (PRRSV) (Rammohan et al., 2012) , porcine circovirus type 2 (PCV2) (Ellis et al., 2004; Genzow et al., 2009) , influenza A virus (IAV) and porcine respiratory corona virus (PRCV) (Pensaert et al., 1986; Jung et al., 2009; Renukaradhya et al., 2010) . Bacteria such as Mycoplasma hyopneumoniae (M. hyopneumoniae) (Hansen et al., 2010b) , Pasteurella multocida (P. multocida) (Davies, 2004) , Salmonella enterica serovar Choleraesuis (S.e Choleraesuis) (Reed et al., 1986; Asai et al., 2010) and Streptococcus suis (S. suis) (Done and Paton, 1995; Silva et al., 2006; Baums et al., 2007) are also commonly associated with PRDC.

Porcine reproductive and respiratory syndrome virus is a major cause of swine production losses worldwide, and in the United States, reproductive and growing pig losses are an estimated $560 million per year (Neumann et al., 2005) . PRRSV (genus Arterivirus, family Arteriviridae) is an enveloped virus with a single-stranded, positive-sense RNA genome of approximately 15 kb. PRRSV is classified into two types with type 1 predominating in Europe and type 2 predominating in North America and Asia. PCV2 (genus Circovirus, family Ciroviridae) is a non-enveloped virus with a circular, covalently closed single-stranded DNA genome of 1767-1768 nucleotides (Meehan et al., 1997; Hamel et al., 1998) . IAV (genus Influenzavirus A, family Orthomyxoviridae) are enveloped viruses with a genome composed of eight single-stranded negative-sense RNA segments. PRCV (genus Alphacoronavirus, family Coronaviridae) are enveloped viruses with a single-stranded positivesense RNA genome of approximately 28.5 kb. It was first isolated in Belgium in 1984, and it is a natural variant of transmissible gastroenteritis virus (TGEV) (Pensaert et al., 1986 ) that contains a 5 0 deletion (621-681 nt in size) in the S gene which is used to differentiate PRCV from TGEV in PCR-based assays (Kim et al., 2000; Costantini et al., 2004) . PRCV often causes subclinical infections. M. hyopneumoniae (Hansen et al., 2010b) causes mycoplasmal pneumonia of swine and is known to be one of the most common and economically important diseases found in pig farms worldwide, having low mortality, but high morbidity (Otagiri et al., 2005) . P. multocida is a commensal found in the upper respiratory tract of pigs, but can also act as a primary or secondary pathogen responsible for pneumonia (Davies, 2004) and atrophic rhinitis in pigs . S. e. Choleraesuis is a host adapted Salmonella serovar responsible for almost all types of salmonellosis in pigs in North America and Europe (Kingsley and Baumler, 2000) . S. suis is a Gram-positive, zoonotic bacterial pathogen important in polyserositis, septicaemia, arthritis, pneumonia and endocarditis in pigs (Done and Paton, 1995; Silva et al., 2006; Baums et al., 2007) . Some of these bacteria are of low pathogenicity and exist as commensals in the upper respiratory tract of healthy pigs. They can cause severe respiratory disease by invading tissues already damaged due to a primary pathogen(s) (virus or bacteria), general immunosuppression and other factors such as poor environmental conditions and poor management practices.

Routine diagnostic methods for detection of viruses implicated in PRDC include virus isolation in cell culture, antigen detection by direct fluorescent antibody staining, and enzyme immunoassay and culture-based methods for bacteria. Such methods (Grau-Roma and Segales, 2007) are time-consuming and require independent tests for each pathogen. Furthermore, the detection of bacterial pathogens typically depends on time-consuming culture-based methods that can take several days to obtain results. Due to their high sensitivity and ease of use, PCR and real-time PCR tests have been developed for several agents implicated in the PRDC; however, these tests typically target single pathogens (Lierz et al., 2008; Lomonaco et al., 2009; Tang et al., 2009) . A diagnostic test capable of simultaneous detection of multiple pathogens involved in PRDC (Atashpaz et al., 2009; Wernike et al., 2012; Xu et al., 2012) would save time, labour and cost by providing more information with each test performed. A multiplex PCR assay capable of detecting five porcine viruses including two porcine respiratory viruses was developed (Giammarioli et al., 2008) . Duplex and triplex real-time PCR for porcine respiratory viruses have also been recently described (Chang et al., 2014; Wu et al., 2014) . However, to date, there are no diagnostic tests capable of simultaneous detection of multiple major viral and bacterial porcine respiratory pathogens in a single reaction. Microarray technology, with its capacity to incorporate a large number of capture probes, is a potentially useful tool for multiplexed detection and typing of pathogens. Here, a microarray assay with associated multiplex RT-PCRs for detection and differentiation of four viruses and four bacteria involved in PRDC using a novel user-friendly electronic microarray in which capture probe printing, hybridization, washing and reporting are fully integrated and automated is described. The electronic microarray contains 400 test sites which can be activated independently via electrodes and utilizes electrophoretically driven hybridization that can be completed instantaneously.

Databases containing all available full and partial sequences of the genomic regions encoding the matrix proteins of IAV (n = 4373) and PRRSV (n = 1829), PRCV spike protein (n = 24), and PCV2 capsid protein (n = 2048) were compiled from the National Centre for Biotechnology Information (NCBI). Similarly, databases were compiled for the kmt1 (n = 16) and toxA (n = 50) genes of P. multocida, sly (n = 106) and orfB (n = 48) genes of S. suis, a 640 bp portion of intergenic region ig-070/071 (Gardner and Minion, 2010) of M. hyopneumoniae (n = 6) and a 708 bp region between open reading frames SC4343 and SC4353 previously identified as a metabolic island for S. e. Choleraesuis (n = 1). Sequences were retrieved by searching NCBI's 'Nucleotide' database using the gene names as keywords, as well as performing BLAST homology searches (Altschul et al., 1990 ) with a representative sequence for each target, and downloading the aligned portion of all BLAST hits. Redundant sequences were removed based on accession numbers. Multiple sequence alignments for each genetic target were generated with ClustalX v. 2.0 ( Thompson et al., 1997; Larkin et al., 2007) or MAFFT v. 7.0 (Katoh and Standley, 2013) using default settings. Databases were maintained with either Mega4 (Tamura et al., 2007) or BioEdit Sequence Alignment Editor v. 7.1.9 (Hall, 1999) to ensure that all sequences were correctly oriented and aligned. Similarly, representative whole-genome sequences, as well as full and partial sequences of homologous genes from related and unrelated non-targets such as TGEV, porcine circovirus type 1 (PCV1), as well as other Salmonella enterica serovars, and Mycoplasma species were downloaded for in silico analysis of probe specificity.

Several published PCR primers suitable for the assay were adopted from the literature (Table 1) . Additional PCR primers (Table 1 ) and all target-specific capture probes (Table 2) were designed using either AlleleID (Premier Biosoft International, Palo Alto, CA, USA) or BioEdit Sequence Alignment Editor v. 7.1.9 (Hall, 1999) based on the databases described above. Primers and probes were designed to be 18-25 bp in length, with the melting temperatures between 54 and 65°C and minimal secondary structures (DG ≥ À8.0 kcal/mol). Primers and probes identified by each software were compiled and examined in silico for specificity by BLAST (Altschul et al., 1990) analysis using custom inhouse databases containing representative whole-genome sequences, as well as full or partial sequences of homologous genes from related and relevant non-targets. Primers or probes that showed significant homology to closely related or unrelated non-target species were excluded from further investigation. 

A list of the viruses, bacteria and clinical samples used in this study is presented in Table 3 . The yield of the RNeasy Mini Kit (Qiagen, Mississauga, ON, Canada), QIAamp Viral RNA Mini Kit (Qiagen) and MagMAX Total Nucleic Acid Kit (Ambion, Austin, TX, USA) was evaluated with PCV2 (a non-enveloped DNA virus) and PRRSV (an enveloped RNA virus) as per the manufacturers' recommendation. For these viruses, 100 ll of a quantified stock was serially diluted (10 À2 to 10 À6 ) in swab material from 2-dayold piglets (Prairie Swine Centre, Saskatoon, SK, Canada), previously tested to be negative for the targets. Similarly, the UltraClean Tissue and Cells DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, USA), DNeasy Blood and Tissue Kit (Qiagen) and MagMAX Total Nucleic Acid Kit were tested to evaluate the nucleic acid extraction efficiency of a Gram-positive bacteria (S. suis) and a Gram-negative bacteria (P. multocida). For these bacteria, 100 ll of culture of known CFU/ml was serially diluted in the same swab material as above and extracted using each kit in parallel according to the manufacturers' instructions. The efficiency of the extraction kits was evaluated based on the limits of detection observed after RT-PCR amplification of the extracted material. The DNeasy Blood and Tissue Kit and the Viral RNA Mini extraction kit were the most efficient for the tested bacteria and virus targets, respectively (data not shown). Therefore, all subsequent nucleic acid extractions of laboratory amplified strains were performed using these kits. Following preparation of nucleic acid extractions, the samples were subjected to PCR and microarray analysis. For the determination of the analytical sensitivity of the assay for viral targets, selected genes were amplified using the SuperScript TM III One-Step RT-PCR System with Platinum â Taq DNA polymerase (Life Technologies, Carlsbad, CA, USA) and cloned into the pJET1.2 cloning vector using the CloneJET PCR Cloning Kit (Fisher Scientific, Ottawa, ON, Canada) according to the manufacturers' specifications. Plasmids were extracted from successfully transformed bacteria using the QIAPrep MiniPrep kit (Qiagen) according to the manufacturer's specifications and were confirmed by sequencing (Eurofins Genomics, Huntsville, AL, USA). Vectors containing the target genes for each virus, except for PCV2, were linearized with the Hind III restriction enzyme (Fisher Scientific) and subjected to in vitro transcription using the MEGAScript T7 Transcription Kit (Life Technologies) according to the manufacturer's specifications. Template DNA was eliminated using successive treatments with TURBO DNase (Life Technologies,) before quantifying the RNA using the RNA BR Assay Kit and the Qubit 2.0 fluorometer (Life Technologies) according to the manufacturers' specifications. The RNA was then serially diluted 1 : 10 in UltraPure Distilled water (Life Technologies). The copy number was inferred using an online tool (http://endmemo.com/bio/dnacopynum. php) taking into account the nucleic acid concentration and nucleotide composition of the amplified region of each target. The copy number for the plasmid containing PCV2 capsid protein coding region was inferred based on the nucleic acid concentration and nucleotide composition over the entire plasmid.

For the determination of the analytical sensitivity of the assay for bacterial targets (excluding M. hyopneumoniae), frozen cultures were streaked for single colonies onto 5% sheep blood agar plates (BBL TM Blood Agar Base Infusion Agar; BD Diagnostics, Sparks, MD, USA) and incubated at 37°C overnight. A single colony was inoculated into 5 ml of Miller's LB broth (Life Technologies) and grown overnight at 37°C on a shaking incubator (150 rpm). The overnight LB broths were serially diluted 1 : 10 in PBS, and 100 ll of material was spread onto blood agar plates in triplicate and grown overnight at 37°C for enumeration using the viable plate count method. The cultures were standardized to 3.33 9 10 8 CFU/ml (S. e Choleraesuis and P. multocida) and 3.33 9 10 6 CFU/ml (S. suis), so a 30 ll aliquot from each serial dilution in the series yielded CFUs to the nearest power of base 10 (i.e. 1 9 10 6 , 1 9 10 5 , etc.). For M. hyopneumoniae, an aliquot of genomic DNA was quantified on the Qubit 2.0 fluorometer, and a genomic copy number was inferred based on the nucleic acid continued) concentration and nucleotide composition over the entire genome.

The analytical specificity of the viral and bacterial multiplex PCR assays was assessed by amplifying panels of 14 non-target viruses and 21 bacteria, respectively (Table 3) .

The forward primers were modified with 5 0 -phosphorylation (IDT, Coralville, IA, USA). All reverse primers were modified with either 5 0 -TYE665 â fluorophore using SpC3 â attachment chemistry for the slide microarray or were synthesized with the reverse complementary sequence of the reporter probe at the 5 0 end for the electronic microarray as described in Lung et al. (2012) . RT-PCRs were performed using the SuperScript TM III One-Step RT-PCR System with Platinum â Taq DNA polymerase. A multiplex RT-PCR with 10 primers targeting the genomic regions encoding the IAV and PRRSV matrix proteins, PRCV spike protein, PCV2 capsid protein, as well as an internal control, was developed (Table 1) . A multiplex PCR with 14 primers targeting the kmt1 and toxA genes of P. multocida, sly and orfB genes of S. suis, intergenic space of M. hyopneumoniae, metabolic island of S. e. Choleraesuis and an internal control was developed (Table 1) . A plasmid containing a fragment of the dengue virus genome was used as an internal PCR control for both the bacterial and viral RT-PCRs. Both assays were optimized for buffer and magnesium concentration, annealing temperature, cycle number and internal control concentration. The finalized RT-PCR mixtures consisted of 1 ll of nucleic acid, 0.01 pg internal control, 1 ll of enzyme mix, 1 lM of each primer in 19 reaction buffer in a final volume of 25 ll. Reverse transcription was carried out for 15 min at 60°C, followed by 94°C for 2 min. PCR was carried out for 35 cycles of 94°C for 30 s, 50°C for laboratory samples and 58°C for clinical samples for 15 s, 68°C for 45 s, with a final extension step of 68°C for 7 min. Following PCR, unpurified material was assayed on the QIAxcel Capillary Gel Electrophoresis System (Qiagen) for visualization of amplicons. Analytical sensitivity of each multiplex assay was determined using serial dilutions of quantified DNA or reverse-transcribed RNA as appropriate for each pathogen. The serial dilutions were amplified using the multiplex PCR assays as well as the singleplex PCR for each target. All PCR amplifications were carried out on the Veriti thermocycler (Life Technologies) and visualized on the QIAxcel using the BioCalculator v. 3.0 software (Qiagen). The limit of detection was considered to be the last dilution where amplification was greater than the default threshold on the electropherogram and described in terms of the approximate total number RNA or DNA copies in the sample.

A total of 30 probes for the detection of four target viruses, 25 probes for the detection of four bacterial targets and three control probes were initially screened by passive hybridization on low-cost conventional Epoxy glass slide microarrays (Corning, Corning, NY, USA) that were printed and processed in-house according to protocols described previously (Lung et al., 2011) . The probes were screened against a panel of five isolates representing the four target viruses, and three non-target viruses, and 11 strains representing the four target bacteria, and 11 non-target bacteria (Table 3) . Microarray data were represented using the mean pixel intensity for each probe reaction. Probe reactivity was calculated using the mean pixel fluorescent intensity (MFI) of all probes as a ratio of the non-template control. Probe reactions above 29 the ratio of the non-template control were considered positive. Probes that showed good reactivity and specificity were selected for testing and validation on a novel automated electronic microarray in which capture probes are printed on streptavidin-containing acrylamide hydrogels and hybridization, washing and reporting are automated and computer controlled. Capture probes used on the electronic microarray were modified with 5 0 -biotin group to allow attachment to streptavidin-containing test sites (IDT, Coralville, IA, USA). A selected set of 12 probes targeting the viruses, 14 probes targeting the bacteria and three control probes (a negative probe and two probes targeting the internal control), which exhibited high reactivity and specificity on glass slide microarrays, were selected for validation on the electronic microarray. The viral probes were tested against an expanded panel of 14 strains or isolates of the four target viruses and 14 non-target viruses (Table 3) . Similarly, the bacterial probes were tested against a panel of 20 strains or isolates representing the four target bacteria and 20 related or unrelated non-target bacteria ( Table 3 ). The electronic microarray assays were run using a protocol previously described (Lung et al., 2012) with modifications. The modifications included the replacement of the 'touch down' washing protocol with a 'touch up' protocol in which washing steps were carried out using Low Salt Buffer (Nexogen, Inc., San Diego, CA, USA) with incremental increases rather than decreases in temperature. Images were captured at each temperature increment. The Red Universal Reporter Probe was replaced with a 5 0 -Alexa Fluor 647 modified locked nucleic acid (LNA) variant (5 0 -/5Alex647N/TGT+CA+AGC-G+AT+AT+ACT+GC-3 0 ) (IDT, Coralville, IA, USA) to increase its thermal stability over a more robust range of wash temperatures. All electronic microarray hybridizations were performed in duplicate, and a non-template PCR control (NTC) was included in all experiments. Raw fluorescent intensity (FI) data from all utilized electrodes at each temperature increment were obtained and analysed using Microsoft Excel. For each probe, positive-to-non-template control (P : N) ratios were calculated by dividing the FI value from each analyte by the FI value produced by the NTC. For each assay, samples that produced P : N ratios of 2.0 or greater were considered positive. Average P : N data derived from the microarrays were visualized with a heat map generated using TreeView v. 1.60 (Eisen et al., 1998) .

Oral and nasal swabs were obtained from specific pathogen free pigs at the CFIA Ottawa Laboratory Fallowfield, Ontario, Canada. Oral and nasal swabs were screened for the target viral and bacterial pathogens, and pools of oral and nasal material that were negative for the target bacteria and viruses were used for spiking with target pathogens. Bacterial strains were grown and quantified as described above, and supernatants containing virus from cell culture were used. For M. hyopneumoniae, culture was not performed and a freeze-dried cell pellet purchased from ATCC was used after re-suspension in PBS and 60% glycerol. For inoculations with single agents, 120 ll aliquots of oral and nasal samples were experimentally inoculated with 20 ll of each live virus or bacteria. Samples inoculated with multiple agents were prepared by pooling 20 ll of each pathogen together, and adding 20 ll of the pooled pathogens into 120 ll of oral and nasal material (Table 4 ). Nucleic acid from the full 140 ll volume of the samples was extracted, PCR amplified and assayed on the electronic microarray as described in previous sections. Approximately 30-40 mg of a panel of lung tissue submitted in 2015 to Manitoba Agriculture, Food and Rural Initiatives (MAFRI) for diagnosis of porcine respiratory pathogens was ground in 1.5 ml of RLT buffer (Qiagen), transferred to 400 ll of MagMax lysis buffer and processed for nucleic acid extraction in the KingFisher 96 instrument (Ambion). Five to 8 ll of extracted RNA from 90 ll of elution buffer was tested by real-time RT-PCR at MAFRI, and 5 ll was tested at CFIA by electronic microarray.

Multiplex PCR/RT-PCR Two separate multiplex PCRs were developed for amplification of selected genes of four viruses and four bacteria involved in PRDC, respectively. The multiplex PCR for bacteria consisted of a total of 12 primers for the detection of the four target bacteria, typing of P. multocida and a pair of primers for an internal control (Table 1) . Two genes from S. suis and P. multocida were each targeted for amplification by the multiplex PCR. Amplification and detection of the toxA gene of P. multocida allowed for differentiation of virulent or pathogenic strains from non-pathogenic strains. Primers for the orfB gene were added to a previously designed multiplex PCR with primers for the sly gene to allow detection of S. suis strains that lack the sly gene. The multiplex PCR generated products of the expected sizes ranging from approximately 205-708 bp (Table 1) when a panel of 23 isolates representing the target species, including different serotypes of S. suis were tested (Fig. 1) .

Similarly, a multiplex RT-PCR with 10 primers was developed and used to detect the four viruses and an internal control. The multiplex RT-PCR successfully amplified a panel of 14 targeted viruses and generated amplicons of the expected size of approximately 229-534 bp ( Table 1) . The samples represented both genotypes of PRRSV, as well as different subtypes of IAV (Fig. 1) .

The internal control variably amplified in both the bacterial and viral multiplex PCRs as a result of competitive PCR. In instances where targets were strongly amplified, amplification of the internal control was either weak or absent.

Conventional glass slide microarrays were processed manually as an initial low-cost screening tool to assess the specificity of the probe set (n = 30) designed to detect the four target viruses, distinguish between genotypes 1 and 2 of PRRSV, as well as differentiate PRCV and TGEV, and the probe set (n = 25) designed to detect the four target bacteria and differentiate toxigenic and non-toxigenic strains of P. multocida. A selected set of 12 probes targeting the viruses, 14 probes targeting the bacteria and three control probes (a negative probe and two probes targeting the internal control), which exhibited the highest reactivity and specificity on slide microarrays, were selected for validation on the user-friendly automated electronic microarray. All target viruses and bacteria were accurately detected, and PRRSV and P. multocida were accurately typed using the viral and bacterial probe set on the electronic microarray platform (Fig. 2) . The assay also successfully detected targeted pathogens in clinical lung specimens, as well as porcine oral and nasal swab material experimentally inoculated with single or multiple targets (Table 4 ). The results obtained were consistent with those obtained by singleplex assays with the exception of LT-7 which was positive for P. multocida based on the electronic microarray assay, but negative by bacterial culture.

For both the RT-PCR and microarray assays, the analytical sensitivity varied with the different targets. For some targets, the multiplex assay had comparable sensitivity with the respective singleplex assay, while for other targets, the singleplex assays were more sensitive ( Table 5 ). The multiplex assay was most sensitive for detection of IAV and S. suis for viral and bacterial targets, respectively (Table 5) .

Non-target bacteria samples did not react with the probes on either the conventional slide microarray or electronic microarray. Other than TGEV, a natural variant and highly related virus to PRCV, no other non-target viruses showed amplification in the virus multiplex RT-PCR (data not shown). Due to the strong amplification of TGEV by the multiplex RT-PCR, the internal control failed to amplify (Fig. 2b, sample 21) . However, PRCV and TGEV were differentiated based on amplicon size, and in subsequent microarray characterization (Fig. 2 ).

A user-friendly microarray for the simultaneous detection and differentiation of four viruses and four bacteria associated with PRDC was developed as a test case for a novel automated 'amplicon-to-answer' electronic microarray technology. The electronic microarray integrates and automates all post-PCR steps required for microarray analysis (capture probe printing, hybridization, washing, reporting) and allows for simultaneous identification of eight pathogens, differentiation of the two PRRSV genotypes and pathogenic versus non-pathogenic P. multocida strains.

Although the amplification of bacterial DNA did not require a reverse transcriptase phase, an RT step was included in the PCR protocol as the amplicon yield was better than without the RT step. Using the same protocol would also allow potential combination of the bacterial and viral multiplex PCR into a single multiplex PCR targeting all eight pathogens. A likely explanation for the increased amplification yield observed with RT-PCRs for bacterial targets is RT-PCR could utilize not just genomic DNA, but also RNA transcripts of target genes as template. In addition, the proprietary quantity of Taq polymerase in the RT-PCR kit that was used in the RT-PCRs may be higher than that used in the PCR.

Initial screening of capture probes was performed using a conventional slide microarray platform due to the lower cost of screening large number of probes that were printed in conjunction with other projects. Subsequently, the assay was adapted to a novel, rapid and user-friendly microarray platform that automates and integrates capture probe printing with all post-PCR steps of the assay, including electrophoretically driven hybridization, washing and reporting. The automated electronic microarray assay reduces the labour, time and number of instruments needed to acquire microarray results when compared with conventional microarrays which typically require substantial manual processing, slower passive hybridization of amplicons with capture probes and multiple pieces of equipment. The electronic microarray has an open platform format with 400 test sites that can be individually activated and utilizes a single integrated instrument that automates capture probe printing and microarray processing using parameters set by the user. Thus, the novel technology also reduces the skill level required to perform microarray assays and allows immediate, on-site modification of assays depending on the needs of the enduser. These unique features eliminate the need for anticipation of future needs, and the procurement and storage of manufactured arrays that are designed for a specific set of predetermined pathogens. For example, the system will allow immediate switching between assays that detects/ types all eight pathogens simultaneously to one that detects/types a subset of the pathogens, or to assays for detection and typing of other pathogens. As the user is able to control hybridization, wash and reporting temperature the assay can have excellent specificity and can be used for differentiating variants that are genetically very similar. However, the automated electronic microarray assay requires specialized arrays and investment in instru-

Probes P. multocida (Toxigenic) Fig. 2 . Summary of microarray results from the electronic microarray representing the four bacterial targets (a), the four viral targets (b) on the electronic microarray. The reactivity of specific reactions between targets and each pathogen-specific probes is outlined in yellow. NSBP = non-specific binding probe negative control. P : N ratios ≥ 2.0 are shown in red, and P : N ratios < 2.0 are in black. NTC = no template control. IC = internal control. Amplification and detection of the internal control are not always observed when a target is present in high amounts.

mentation. The cost of consumables for testing a sample with an electronic microarray is expected to be between one and a few real-time PCR assays and will depend on the assay design (i.e. reducing the number of probes will decrease the cost by allowing more samples to be tested on each array with 400 test sites). In addition to glass and hydrogel (acrylamide)-based microarray matrices used in this study, the final probe sets have been tested successfully on reverse dot blots performed on nylon membranes (unpublished results). Thus, the probe set described should have broad applications in hybridization-based assays that use a variety of different matrices. Assays for simultaneous detection of multiple bacteria and virus implicated in PRDC have not been described previously. The limit of detection of a LAMP assay for PCV2 has been shown to be approximately 1 copy of DNA plasmid, much more sensitive than conventional PCR whose detection limit was 1 9 10 4 copies (Zhou et al., 2011) . However, LAMP assays typically detect a single pathogen and are difficult to multiplex. In multiplex real-time PCR assays, a duplex PCR assay had a detection limit of 1 TCID 50 /ml for PCV2 and 6.3 TCID 50 /ml for PRRSV (Chang et al., 2014) , while a triplex real-time PCR assay had a detection sensitivity of 10 copies/ll for PCV2 and PRSV and 100 copies/ll for PPV (Wu et al., 2014) . The multiplex PCR described here did not reduce the limits of detection for five targets although the detection limit of the PCR for three pathogens was an order of magnitude lower in the multiplex format. Thus, further improvements in sensitivity are desirable and may be partly achieved by increasing the number of PCR cycles, reducing amplicon lengths or reducing the amount of internal control used in the assay. The limit of detection for PCV2, PRRSV, and M. hyoneumoniae was lower with the electronic microarray in comparison with the multiplex RT-PCR indicating the use of the electronic microarray platform reduces sensitivity, and unpublished results show that only amplicons that were visible after agarose gel electrophoresis can be detected by the electronic microarray. Despite the reduced analytical sensitivity of the electronic microarray assay for some of the targets, the bacterial and viral targets were detected singly or in combination in clinical samples submitted for laboratory diagnosis. Furthermore, the assay detected P. multocida in a clinical lung sample that was not detected by traditional culture methods. This discrepancy may be due to the higher sensitivity of the microarray assay, or the lack of active infection (e.g. the use of antibiotics).

While primers and probes were evaluated using all sequences available on NCBI at the time of assay design and samples representative of both target and non-target bacteria and viruses were tested in this study, regular reevaluation of the coverage of the primers and probes and additional validation with clinical samples is desired. The four swine respiratory viruses targeted in this study have also been successfully detected with the Virochip panviral array which consists of probes for detection of all known viruses at the time of design (Nicholson et al., 2011) . The high-density Virochip is a useful tool for virus discovery, but needs approximately 24 h to obtain results and requires high viral titre for positive detection. In contrast, the electronic microarray assay described here can be completed in less than 4 hours with little user handling plus approximately 1.5 h for the RT-PCR described. New instrumentation that further simplifies the workflow by integrating the PCR and array processes is now commercially available. To our knowledge, the automated microarray assay described here is the first one that simultaneously 

Lipman, 1990: Basic local alignment search tool

Molecular typing and antimicrobial resistance of Salmonella enterica subspecies enterica serovar Choleraesuis isolates from diseased pigs in Japan

Rapid virulence typing of Pasteurella multocida by multiplex PCR

Prevalence of Streptococcus suis genotypes in wild boars of Northwestern Germany

The application of a duplex reverse transcription real-time PCR for the surveillance of porcine reproductive and respiratory syndrome virus and porcine circovirus type 2

Respiratory and fecal shedding of porcine respiratory coronavirus (PRCV) in sentinel weaned pigs and sequence of the partial S-gene of the PRCV isolates

Genetic diversity among Pasteurella multocida strains of avian, bovine, ovine and porcine origin from England and Wales by comparative sequence analysis of the 16S rRNA gene

Porcine reproductive and respiratory syndrome: clinical disease, pathology and immunosuppression

Cluster analysis and display of genome-wide expression patterns

Porcine circovirus-2 and concurrent infections in the field

Detection and quantification of intergenic transcription in Mycoplasma hyopneumoniae

The effect of vaccination against Porcine reproductive and respiratory syndrome virus (PRRSV) on the Porcine circovirus-2 (PCV-2) load in porcine circovirus associated disease (PCVAD) affected pigs

Development of a novel hot-start multiplex PCR for simultaneous detection of classical swine fever virus, African swine fever virus, porcine circovirus type 2, porcine reproductive and respiratory syndrome virus and porcine parvovirus

Detection of porcine reproductive and respiratory syndrome virus, porcine circovirus type 2, swine influenza virus and Aujeszky's disease virus in cases of porcine proliferative and necrotizing pneumonia (PNP) in Spain

BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT

Nucleotide sequence of porcine circovirus associated with postweaning multisystemic wasting syndrome in pigs

Occurrence and tissue distribution of porcine circovirus type 2 identified by immunohistochemistry in Danish finishing pigs at slaughter

An investigation of the pathology and pathogens associated with porcine respiratory disease complex in Denmark

Diagnostic evaluation of multiplexed reverse transcription-PCR microsphere array assay for detection of foot-and-mouth and lookalike disease viruses

Porcine reproductive and respiratory syndrome virus modifies innate immunity and alters disease outcome in pigs subsequently infected with porcine respiratory coronavirus: implications for respiratory viral co-infections

MAFFT multiple sequence alignment software version 7: improvements in performance and usability

Molecular characterization and pathogenesis of transmissible gastroenteritis coronavirus (TGEV) and porcine respiratory coronavirus (PRCV) field isolates co-circulating in a swine herd

Host adaptation and the emergence of infectious disease: the Salmonella paradigm

Higgins, 2007: Clustal W and Clustal X version 2.0

Direct PCR analysis for toxigenic Pasteurella multocida

Occurrence of mycoplasmas in free-ranging birds of prey in Germany

Detection of Salmonella in finishing pigs on farm

Multiplex RT-PCR detection and microarray typing of vesicular disease viruses

Electronic microarray assays for avian influenza and Newcastle disease virus

Sequence of porcine circovirus DNA: affinities with plant circoviruses

Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States

Utility of a panviral microarray for detection of swine respiratory viruses in clinical samples

Polymicrobial respiratory disease in pigs

Detection of Mycoplasma hyopneumoniae in lung and nasal swab samples from pigs by nested PCR and culture methods

Isolation of a porcine respiratory, non-enteric coronavirus related to transmissible gastroenteritis

Increased prevalence of torque teno viruses in porcine respiratory disease complex affected pigs

Studies on the pathogenesis of Salmonella typhimurium and Salmonella choleraesuis var kunzendorf infection in weanling pigs

Porcine reproductive and respiratory syndrome virusinduced immunosuppression exacerbates the inflammatory response to porcine respiratory coronavirus in pigs

Virulence-associated gene profiling of Streptococcus suis isolates by PCR

Rapid diagnosis of dengue viremia by reverse transcriptase-polymerase chain reaction using 3'-noncoding region universal primers

MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0

Isolation, antimicrobial resistance, and virulence genes of Pasteurella multocida strains from swine in China

The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools

Part III: Reference of Swine Health and Health Management in the United States

Detection and typing of highly pathogenic porcine reproductive and respiratory syndrome virus by multiplex real-time rt-PCR

Rapid multiplex PCR and real-time TaqMan PCR assays for detection of Salmonella enterica and the highly virulent serovars Choleraesuis and Paratyphi C

A sensitive multiplex real-time PCR panel for rapid diagnosis of viruses associated with porcine respiratory and reproductive disorders

2012: Development of multiplex PCR for simultaneous detection of six swine DNA and RNA viruses

Loop-mediated isothermal amplification for detection of porcine circovirus type 2

1.4 9 10 3 1.4 9 10 3 1.4 9 10 3 1.4 9 10 4 PRCV 62 620 620 620 PRRSV 3.8 9 10 3 3.8 9 10 3 3.8 9 10 4 3.8 9 10 4 IAV 160 160 160 160 M. hyopneumoniae 480 480 480 4.8 9 10 3 P. multocida 1.0 9 10 3 1.0 9 10 3 1.0 9 10 3 1.0 9 10 3 S. e. Choleraesuis 1.0 9 10 3 1.0 9 10 4 1.0 9 10 3 1.0 9 10 4 S. suis 

None of the authors of this paper has a financial or personal relationship with other people or organizations that could inappropriately influence or bias the content of the paper.