Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Multiplex primer prediction software for divergent targets</title><author><name sortKey="Gardner, Shea N" sort="Gardner, Shea N" uniqKey="Gardner S" first="Shea N." last="Gardner">Shea N. Gardner</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Hiddessen, Amy L" sort="Hiddessen, Amy L" uniqKey="Hiddessen A" first="Amy L." last="Hiddessen">Amy L. Hiddessen</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Williams, Peter L" sort="Williams, Peter L" uniqKey="Williams P" first="Peter L." last="Williams">Peter L. Williams</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Hara, Christine" sort="Hara, Christine" uniqKey="Hara C" first="Christine" last="Hara">Christine Hara</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Wagner, Mark C" sort="Wagner, Mark C" uniqKey="Wagner M" first="Mark C." last="Wagner">Mark C. Wagner</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Colston, Bill W" sort="Colston, Bill W" uniqKey="Colston B" first="Bill W." last="Colston">Bill W. Colston</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">19759213</idno><idno type="pmc">2770652</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2770652</idno><idno type="RBID">PMC:2770652</idno><idno type="doi">10.1093/nar/gkp659</idno><date when="2009">2009</date><idno type="wicri:Area/Pmc/Corpus">000F45</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000F45</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Multiplex primer prediction software for divergent targets</title><author><name sortKey="Gardner, Shea N" sort="Gardner, Shea N" uniqKey="Gardner S" first="Shea N." last="Gardner">Shea N. Gardner</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Hiddessen, Amy L" sort="Hiddessen, Amy L" uniqKey="Hiddessen A" first="Amy L." last="Hiddessen">Amy L. Hiddessen</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Williams, Peter L" sort="Williams, Peter L" uniqKey="Williams P" first="Peter L." last="Williams">Peter L. Williams</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Hara, Christine" sort="Hara, Christine" uniqKey="Hara C" first="Christine" last="Hara">Christine Hara</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Wagner, Mark C" sort="Wagner, Mark C" uniqKey="Wagner M" first="Mark C." last="Wagner">Mark C. Wagner</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author><author><name sortKey="Colston, Bill W" sort="Colston, Bill W" uniqKey="Colston B" first="Bill W." last="Colston">Bill W. Colston</name><affiliation><nlm:aff id="AFF1"></nlm:aff></affiliation></author></analytic><series><title level="j">Nucleic Acids Research</title><idno type="ISSN">0305-1048</idno><idno type="eISSN">1362-4962</idno><imprint><date when="2009">2009</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><title>Abstract</title><p>We describe a Multiplex Primer Prediction (MPP) algorithm to build multiplex compatible primer sets to amplify all members of large, diverse and unalignable sets of target sequences. The MPP algorithm is scalable to larger target sets than other available software, and it does not require a multiple sequence alignment. We applied it to questions in viral detection, and demonstrated that there are no universally conserved priming sequences among viruses and that it could require an unfeasibly large number of primers (∼3700 18-mers or ∼2000 10-mers) to generate amplicons from all sequenced viruses. We then designed primer sets separately for each viral family, and for several diverse species such as foot-and-mouth disease virus (FMDV), hemagglutinin (HA) and neuraminidase (NA) segments of influenza A virus, Norwalk virus, and HIV-1. We empirically demonstrated the application of the software with a multiplex set of 16 short (10 nt) primers designed to amplify the Poxviridae family to produce a specific amplicon from vaccinia virus.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Victoria, Jg" uniqKey="Victoria J">JG Victoria</name></author><author><name sortKey="Kapoor, A" uniqKey="Kapoor A">A Kapoor</name></author><author><name sortKey="Li, L" uniqKey="Li L">L Li</name></author><author><name sortKey="Blinkova, O" uniqKey="Blinkova O">O Blinkova</name></author><author><name sortKey="Slikas, B" uniqKey="Slikas B">B Slikas</name></author><author><name sortKey="Wang, C" uniqKey="Wang C">C Wang</name></author><author><name sortKey="Naeem, A" uniqKey="Naeem A">A Naeem</name></author><author><name sortKey="Zaidi, S" uniqKey="Zaidi S">S Zaidi</name></author><author><name sortKey="Delwart, E" uniqKey="Delwart E">E Delwart</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edwards, Ra" uniqKey="Edwards R">RA Edwards</name></author><author><name sortKey="Rohwer, F" uniqKey="Rohwer F">F Rohwer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nakamura, S" uniqKey="Nakamura S">S Nakamura</name></author><author><name sortKey="Yang, C S" uniqKey="Yang C">C-S Yang</name></author><author><name sortKey="Sakon, N" uniqKey="Sakon N">N Sakon</name></author><author><name sortKey="Ueda, M" uniqKey="Ueda M">M Ueda</name></author><author><name sortKey="Tougan, T" uniqKey="Tougan T">T Tougan</name></author><author><name sortKey="Yamashita, A" uniqKey="Yamashita A">A Yamashita</name></author><author><name sortKey="Goto, N" uniqKey="Goto N">N Goto</name></author><author><name sortKey="Takahashi, K" uniqKey="Takahashi K">K Takahashi</name></author><author><name sortKey="Yasunaga, T" uniqKey="Yasunaga T">T Yasunaga</name></author><author><name sortKey="Ikuta, K" uniqKey="Ikuta K">K Ikuta</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kistler, A" uniqKey="Kistler A">A Kistler</name></author><author><name sortKey="Avila, P" uniqKey="Avila P">P Avila</name></author><author><name sortKey="Rouskin, S" uniqKey="Rouskin S">S Rouskin</name></author><author><name sortKey="Wang, D" uniqKey="Wang D">D Wang</name></author><author><name sortKey="Ward, T" uniqKey="Ward T">T Ward</name></author><author><name sortKey="Yagi, S" uniqKey="Yagi S">S Yagi</name></author><author><name sortKey="Schnurr, D" uniqKey="Schnurr D">D Schnurr</name></author><author><name sortKey="Ganem, D" uniqKey="Ganem D">D Ganem</name></author><author><name sortKey="Derisi, J" uniqKey="Derisi J">J DeRisi</name></author><author><name sortKey="Boushey, H" uniqKey="Boushey H">H Boushey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, B" uniqKey="Lin B">B Lin</name></author><author><name sortKey="Blaney, Km" uniqKey="Blaney K">KM Blaney</name></author><author><name sortKey="Malanoski, Ap" uniqKey="Malanoski A">AP Malanoski</name></author><author><name sortKey="Ligler, Ag" uniqKey="Ligler A">AG Ligler</name></author><author><name sortKey="Schnur, Jm" uniqKey="Schnur J">JM Schnur</name></author><author><name sortKey="Metzgar, D" uniqKey="Metzgar D">D Metzgar</name></author><author><name sortKey="Russell, Kl" uniqKey="Russell K">KL Russell</name></author><author><name sortKey="Stenger, Da" uniqKey="Stenger D">DA Stenger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, B" uniqKey="Lin B">B Lin</name></author><author><name sortKey="Wang, Z" uniqKey="Wang Z">Z Wang</name></author><author><name sortKey="Vora, Gj" uniqKey="Vora G">GJ Vora</name></author><author><name sortKey="Thornton, Ja" uniqKey="Thornton J">JA Thornton</name></author><author><name sortKey="Schnur, Jm" uniqKey="Schnur J">JM Schnur</name></author><author><name sortKey="Thach, Dc" uniqKey="Thach D">DC Thach</name></author><author><name sortKey="Blaney, Km" uniqKey="Blaney K">KM Blaney</name></author><author><name sortKey="Ligler, Ag" uniqKey="Ligler A">AG Ligler</name></author><author><name sortKey="Malanoski, Ap" uniqKey="Malanoski A">AP Malanoski</name></author><author><name sortKey="Santiago, J" uniqKey="Santiago J">J Santiago</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, D" uniqKey="Wang D">D Wang</name></author><author><name sortKey="Coscoy, L" uniqKey="Coscoy L">L Coscoy</name></author><author><name sortKey="Zylberberg, M" uniqKey="Zylberberg M">M Zylberberg</name></author><author><name sortKey="Avila, P" uniqKey="Avila P">P Avila</name></author><author><name sortKey="Boushey, H" uniqKey="Boushey H">H Boushey</name></author><author><name sortKey="Ganem, D" uniqKey="Ganem D">D Ganem</name></author><author><name sortKey="Derisi, J" uniqKey="Derisi J">J DeRisi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Quan, Pl" uniqKey="Quan P">PL Quan</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Jabado, Oj" uniqKey="Jabado O">OJ Jabado</name></author><author><name sortKey="Conlan, S" uniqKey="Conlan S">S Conlan</name></author><author><name sortKey="Hirschberg, Dl" uniqKey="Hirschberg D">DL Hirschberg</name></author><author><name sortKey="Pozo, F" uniqKey="Pozo F">F Pozo</name></author><author><name sortKey="Jack, Pjm" uniqKey="Jack P">PJM Jack</name></author><author><name sortKey="Cisterna, D" uniqKey="Cisterna D">D Cisterna</name></author><author><name sortKey="Renwick, N" uniqKey="Renwick N">N Renwick</name></author><author><name sortKey="Hui, J" uniqKey="Hui J">J Hui</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Quan, P" uniqKey="Quan P">P Quan</name></author><author><name sortKey="Jabado, O" uniqKey="Jabado O">O Jabado</name></author><author><name sortKey="Conlan, S" uniqKey="Conlan S">S Conlan</name></author><author><name sortKey="Hirschberg, D" uniqKey="Hirschberg D">D Hirschberg</name></author><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name></author><author><name sortKey="Zhai, J" uniqKey="Zhai J">J Zhai</name></author><author><name sortKey="Renwick, N" uniqKey="Renwick N">N Renwick</name></author><author><name sortKey="Hui, J" uniqKey="Hui J">J Hui</name></author><author><name sortKey="Hegyi, H" uniqKey="Hegyi H">H Hegyi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ecker, D" uniqKey="Ecker D">D Ecker</name></author><author><name sortKey="Drader, J" uniqKey="Drader J">J Drader</name></author><author><name sortKey="Gutierrez, J" uniqKey="Gutierrez J">J Gutierrez</name></author><author><name sortKey="Gutierrez, A" uniqKey="Gutierrez A">A Gutierrez</name></author><author><name sortKey="Hannis, J" uniqKey="Hannis J">J Hannis</name></author><author><name sortKey="Schink, A" uniqKey="Schink A">A Schink</name></author><author><name sortKey="Sampath, R" uniqKey="Sampath R">R Sampath</name></author><author><name sortKey="Blyn, Lb" uniqKey="Blyn L">LB Blyn</name></author><author><name sortKey="Eschoo, Mw" uniqKey="Eschoo M">MW Eschoo</name></author><author><name sortKey="Hall, Ta" uniqKey="Hall T">TA Hall</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sampath, R" uniqKey="Sampath R">R Sampath</name></author><author><name sortKey="Hall, Ta" uniqKey="Hall T">TA Hall</name></author><author><name sortKey="Massire, C" uniqKey="Massire C">C Massire</name></author><author><name sortKey="Li, F" uniqKey="Li F">F Li</name></author><author><name sortKey="Blyn, Lb" uniqKey="Blyn L">LB Blyn</name></author><author><name sortKey="Eshoo, Mw" uniqKey="Eshoo M">MW Eshoo</name></author><author><name sortKey="Hofstadler, Sa" uniqKey="Hofstadler S">SA Hofstadler</name></author><author><name sortKey="Ecker, Dj" uniqKey="Ecker D">DJ Ecker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lamson, D" uniqKey="Lamson D">D Lamson</name></author><author><name sortKey="Renwick, N" uniqKey="Renwick N">N Renwick</name></author><author><name sortKey="Kapoor, V" uniqKey="Kapoor V">V Kapoor</name></author><author><name sortKey="Liu, Z" uniqKey="Liu Z">Z Liu</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Ju, J" uniqKey="Ju J">J Ju</name></author><author><name sortKey="Dean, A" uniqKey="Dean A">A Dean</name></author><author><name sortKey="St George, K" uniqKey="St George K">K St.George</name></author><author><name sortKey="Briese, T" uniqKey="Briese T">T Briese</name></author><author><name sortKey="Lipkin, Wi" uniqKey="Lipkin W">WI Lipkin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Briese, T" uniqKey="Briese T">T Briese</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Kokoris, M" uniqKey="Kokoris M">M Kokoris</name></author><author><name sortKey="Jabado, O" uniqKey="Jabado O">O Jabado</name></author><author><name sortKey="Zhiqiang, L" uniqKey="Zhiqiang L">L Zhiqiang</name></author><author><name sortKey="Renwick, N" uniqKey="Renwick N">N Renwick</name></author><author><name sortKey="Kapoor, V" uniqKey="Kapoor V">V Kapoor</name></author><author><name sortKey="Casas, I" uniqKey="Casas I">I Casas</name></author><author><name sortKey="Pozo, F" uniqKey="Pozo F">F Pozo</name></author><author><name sortKey="Limberger, R" uniqKey="Limberger R">R Limberger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dominguez, Sr" uniqKey="Dominguez S">SR Dominguez</name></author><author><name sortKey="Briese, T" uniqKey="Briese T">T Briese</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Hui, J" uniqKey="Hui J">J Hui</name></author><author><name sortKey="Villari, J" uniqKey="Villari J">J Villari</name></author><author><name sortKey="Kapoor, V" uniqKey="Kapoor V">V Kapoor</name></author><author><name sortKey="Tokarz, R" uniqKey="Tokarz R">R Tokarz</name></author><author><name sortKey="Glode, Mp" uniqKey="Glode M">MP Glodé</name></author><author><name sortKey="Anderson, Ms" uniqKey="Anderson M">MS Anderson</name></author><author><name sortKey="Robinson, Cc" uniqKey="Robinson C">CC Robinson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fox, Jd" uniqKey="Fox J">JD Fox</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lenhoff, Rj" uniqKey="Lenhoff R">RJ Lenhoff</name></author><author><name sortKey="Naraghi Arani, P" uniqKey="Naraghi Arani P">P Naraghi-Arani</name></author><author><name sortKey="Thissen, Jb" uniqKey="Thissen J">JB Thissen</name></author><author><name sortKey="Olivas, J" uniqKey="Olivas J">J Olivas</name></author><author><name sortKey="Carillo, Ac" uniqKey="Carillo A">AC Carillo</name></author><author><name sortKey="Chinn, C" uniqKey="Chinn C">C Chinn</name></author><author><name sortKey="Rasmussen, M" uniqKey="Rasmussen M">M Rasmussen</name></author><author><name sortKey="Messenger, Sm" uniqKey="Messenger S">SM Messenger</name></author><author><name sortKey="Suer, Ld" uniqKey="Suer L">LD Suer</name></author><author><name sortKey="Smith, Sm" uniqKey="Smith S">SM Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Griffiths, D" uniqKey="Griffiths D">D Griffiths</name></author><author><name sortKey="Kellam, P" uniqKey="Kellam P">P Kellam</name></author><author><name sortKey="Weiss, R" uniqKey="Weiss R">R Weiss</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Quan, Pl" uniqKey="Quan P">PL Quan</name></author><author><name sortKey="Briese, T" uniqKey="Briese T">T Briese</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Lipkin, Wi" uniqKey="Lipkin W">WI Lipkin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lemmon, Gh" uniqKey="Lemmon G">GH Lemmon</name></author><author><name sortKey="Gardner, Sn" uniqKey="Gardner S">SN Gardner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jabado, Oj" uniqKey="Jabado O">OJ Jabado</name></author><author><name sortKey="Palacios, G" uniqKey="Palacios G">G Palacios</name></author><author><name sortKey="Kapoor, V" uniqKey="Kapoor V">V Kapoor</name></author><author><name sortKey="Hui, J" uniqKey="Hui J">J Hui</name></author><author><name sortKey="Renwick, N" uniqKey="Renwick N">N Renwick</name></author><author><name sortKey="Zhai, J" uniqKey="Zhai J">J Zhai</name></author><author><name sortKey="Briese, T" uniqKey="Briese T">T Briese</name></author><author><name sortKey="Lipkin, Wi" uniqKey="Lipkin W">WI Lipkin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jarman, Sn" uniqKey="Jarman S">SN Jarman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Linhart, C" uniqKey="Linhart C">C Linhart</name></author><author><name sortKey="Shamir, R" uniqKey="Shamir R">R Shamir</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rose, T" uniqKey="Rose T">T Rose</name></author><author><name sortKey="Henikoff, J" uniqKey="Henikoff J">J Henikoff</name></author><author><name sortKey="Henikoff, S" uniqKey="Henikoff S">S Henikoff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gadberry, Md" uniqKey="Gadberry M">MD Gadberry</name></author><author><name sortKey="Malcomber, St" uniqKey="Malcomber S">ST Malcomber</name></author><author><name sortKey="Doust, An" uniqKey="Doust A">AN Doust</name></author><author><name sortKey="Kellogg, Ea" uniqKey="Kellogg E">EA Kellogg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pesole, G" uniqKey="Pesole G">G Pesole</name></author><author><name sortKey="Liuni, S" uniqKey="Liuni S">S Liuni</name></author><author><name sortKey="Grillo, G" uniqKey="Grillo G">G Grillo</name></author><author><name sortKey="Belichard, P" uniqKey="Belichard P">P Belichard</name></author><author><name sortKey="Trenkle, T" uniqKey="Trenkle T">T Trenkle</name></author><author><name sortKey="Welsh, J" uniqKey="Welsh J">J Welsh</name></author><author><name sortKey="Mcclelland, M" uniqKey="Mcclelland M">M McClelland</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huang, Yc" uniqKey="Huang Y">YC Huang</name></author><author><name sortKey="Chang, Cf" uniqKey="Chang C">CF Chang</name></author><author><name sortKey="Chan, Ch" uniqKey="Chan C">Ch Chan</name></author><author><name sortKey="Yeh, Tj" uniqKey="Yeh T">TJ Yeh</name></author><author><name sortKey="Chang, Yc" uniqKey="Chang Y">YC Chang</name></author><author><name sortKey="Chen, Cc" uniqKey="Chen C">CC Chen</name></author><author><name sortKey="Kao, Cy" uniqKey="Kao C">CY Kao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Markham, N" uniqKey="Markham N">N Markham</name></author><author><name sortKey="Zuker, M" uniqKey="Zuker M">M Zuker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rachlin, J" uniqKey="Rachlin J">J Rachlin</name></author><author><name sortKey="Ding, C" uniqKey="Ding C">C Ding</name></author><author><name sortKey="Cantor, C" uniqKey="Cantor C">C Cantor</name></author><author><name sortKey="Kasif, S" uniqKey="Kasif S">S Kasif</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rachlin, J" uniqKey="Rachlin J">J Rachlin</name></author><author><name sortKey="Ding, C" uniqKey="Ding C">C Ding</name></author><author><name sortKey="Cantor, C" uniqKey="Cantor C">C Cantor</name></author><author><name sortKey="Kasif, S" uniqKey="Kasif S">S Kasif</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edgar, Rc" uniqKey="Edgar R">RC Edgar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chenna, R" uniqKey="Chenna R">R Chenna</name></author><author><name sortKey="Sugawara, H" uniqKey="Sugawara H">H Sugawara</name></author><author><name sortKey="Koike, T" uniqKey="Koike T">T Koike</name></author><author><name sortKey="Lopez, R" uniqKey="Lopez R">R Lopez</name></author><author><name sortKey="Gibson, Tj" uniqKey="Gibson T">TJ Gibson</name></author><author><name sortKey="Higgins, Dg" uniqKey="Higgins D">DG Higgins</name></author><author><name sortKey="Thompson, Jd" uniqKey="Thompson J">JD Thompson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Eddy, Sr" uniqKey="Eddy S">SR Eddy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Herold, J" uniqKey="Herold J">J Herold</name></author><author><name sortKey="Kurtz, S" uniqKey="Kurtz S">S Kurtz</name></author><author><name sortKey="Giegerich, R" uniqKey="Giegerich R">R Giegerich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name></author><author><name sortKey="Madden, T" uniqKey="Madden T">T Madden</name></author><author><name sortKey="Schaffer, A" uniqKey="Schaffer A">A Schaffer</name></author><author><name sortKey="Zhang, J" uniqKey="Zhang J">J Zhang</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author><author><name sortKey="Lipman, D" uniqKey="Lipman D">D Lipman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Casabianca, A" uniqKey="Casabianca A">A Casabianca</name></author><author><name sortKey="Gori, C" uniqKey="Gori C">C Gori</name></author><author><name sortKey="Orlandi, C" uniqKey="Orlandi C">C Orlandi</name></author><author><name sortKey="Forbici, F" uniqKey="Forbici F">F Forbici</name></author><author><name sortKey="Federico Perno, C" uniqKey="Federico Perno C">C Federico Perno</name></author><author><name sortKey="Magnani, M" uniqKey="Magnani M">M Magnani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Putkuri, N" uniqKey="Putkuri N">N Putkuri</name></author><author><name sortKey="Piiparinen, H" uniqKey="Piiparinen H">H Piiparinen</name></author><author><name sortKey="Vaheri, A" uniqKey="Vaheri A">A Vaheri</name></author><author><name sortKey="Vapalahti, O" uniqKey="Vapalahti O">O Vapalahti</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Olson, Va" uniqKey="Olson V">VA Olson</name></author><author><name sortKey="Laue, T" uniqKey="Laue T">T Laue</name></author><author><name sortKey="Laker, Mt" uniqKey="Laker M">MT Laker</name></author><author><name sortKey="Babkin, Iv" uniqKey="Babkin I">IV Babkin</name></author><author><name sortKey="Drosten, C" uniqKey="Drosten C">C Drosten</name></author><author><name sortKey="Shchelkunov, Sn" uniqKey="Shchelkunov S">SN Shchelkunov</name></author><author><name sortKey="Niedrig, M" uniqKey="Niedrig M">M Niedrig</name></author><author><name sortKey="Damon, Ik" uniqKey="Damon I">IK Damon</name></author><author><name sortKey="Meyer, H" uniqKey="Meyer H">H Meyer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pearson, Wr" uniqKey="Pearson W">WR Pearson</name></author><author><name sortKey="Robbins, G" uniqKey="Robbins G">G Robbins</name></author><author><name sortKey="Wrege, De" uniqKey="Wrege D">DE Wrege</name></author><author><name sortKey="Zhang, T" uniqKey="Zhang T">T Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Palka Santini, M" uniqKey="Palka Santini M">M Palka-Santini</name></author><author><name sortKey="Cleven, B" uniqKey="Cleven B">B Cleven</name></author><author><name sortKey="Eichinger, L" uniqKey="Eichinger L">L Eichinger</name></author><author><name sortKey="Kronke, M" uniqKey="Kronke M">M Kronke</name></author><author><name sortKey="Krut, O" uniqKey="Krut O">O Krut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Beer, Nr" uniqKey="Beer N">NR Beer</name></author><author><name sortKey="Hindson, Bj" uniqKey="Hindson B">BJ Hindson</name></author><author><name sortKey="Wheeler, Ek" uniqKey="Wheeler E">EK Wheeler</name></author><author><name sortKey="Hall, Sb" uniqKey="Hall S">SB Hall</name></author><author><name sortKey="Rose, Ka" uniqKey="Rose K">KA Rose</name></author><author><name sortKey="Kennedy, Im" uniqKey="Kennedy I">IM Kennedy</name></author><author><name sortKey="Colston, Bw" uniqKey="Colston B">BW Colston</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Caetano Anolles, G" uniqKey="Caetano Anolles G">G Caetano-Anolles</name></author><author><name sortKey="Brant, Jb" uniqKey="Brant J">JB Brant</name></author><author><name sortKey="Peter, Mg" uniqKey="Peter M">MG Peter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Williams, Jgk" uniqKey="Williams J">JGK Williams</name></author><author><name sortKey="Kubelik, Ar" uniqKey="Kubelik A">AR Kubelik</name></author><author><name sortKey="Livak, Kj" uniqKey="Livak K">KJ Livak</name></author><author><name sortKey="Rafalski, Ja" uniqKey="Rafalski J">JA Rafalski</name></author><author><name sortKey="Tingey, Sv" uniqKey="Tingey S">SV Tingey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nanda, S" uniqKey="Nanda S">S Nanda</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nanda, S" uniqKey="Nanda S">S Nanda</name></author><author><name sortKey="Jayan, G" uniqKey="Jayan G">G Jayan</name></author><author><name sortKey="Voulgaropoulou, F" uniqKey="Voulgaropoulou F">F Voulgaropoulou</name></author><author><name sortKey="Sierra Honigmann, Am" uniqKey="Sierra Honigmann A">AM Sierra-Honigmann</name></author><author><name sortKey="Uhlenhaut, C" uniqKey="Uhlenhaut C">C Uhlenhaut</name></author><author><name sortKey="Mcwatters, Bjp" uniqKey="Mcwatters B">BJP McWatters</name></author><author><name sortKey="Patel, A" uniqKey="Patel A">A Patel</name></author><author><name sortKey="Krause, Pr" uniqKey="Krause P">PR Krause</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pyrc, K" uniqKey="Pyrc K">K Pyrc</name></author><author><name sortKey="Jebbink, Mf" uniqKey="Jebbink M">MF Jebbink</name></author><author><name sortKey="Berkhout, B" uniqKey="Berkhout B">B Berkhout</name></author><author><name sortKey="Hoek, L" uniqKey="Hoek L">L Hoek</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Vries, M" uniqKey="De Vries M">M de Vries</name></author><author><name sortKey="Pyrc, K" uniqKey="Pyrc K">K Pyrc</name></author><author><name sortKey="Berkhout, R" uniqKey="Berkhout R">R Berkhout</name></author><author><name sortKey="Vermeulen Oost, W" uniqKey="Vermeulen Oost W">W Vermeulen-Oost</name></author><author><name sortKey="Dijkman, R" uniqKey="Dijkman R">R Dijkman</name></author><author><name sortKey="Jebbink, Mf" uniqKey="Jebbink M">MF Jebbink</name></author><author><name sortKey="Bruisten, S" uniqKey="Bruisten S">S Bruisten</name></author><author><name sortKey="Berkhout, B" uniqKey="Berkhout B">B Berkhout</name></author><author><name sortKey="Van Der Hoek, L" uniqKey="Van Der Hoek L">L van der Hoek</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Nucleic Acids Res</journal-id><journal-id journal-id-type="iso-abbrev">Nucleic Acids Res</journal-id><journal-id journal-id-type="publisher-id">nar</journal-id><journal-id journal-id-type="hwp">nar</journal-id><journal-title-group><journal-title>Nucleic Acids Research</journal-title></journal-title-group><issn pub-type="ppub">0305-1048</issn><issn pub-type="epub">1362-4962</issn><publisher><publisher-name>Oxford University Press</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">19759213</article-id><article-id pub-id-type="pmc">2770652</article-id><article-id pub-id-type="doi">10.1093/nar/gkp659</article-id><article-id pub-id-type="publisher-id">gkp659</article-id><article-categories><subj-group subj-group-type="heading"><subject>Computational Biology</subject></subj-group></article-categories><title-group><article-title>Multiplex primer prediction software for divergent targets</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Gardner</surname><given-names>Shea N.</given-names></name><xref ref-type="aff" rid="AFF1">1</xref><xref ref-type="corresp" rid="COR1">*</xref><pmc-comment>gardner26@llnl.gov </pmc-comment>
</contrib><contrib contrib-type="author"><name><surname>Hiddessen</surname><given-names>Amy L.</given-names></name><xref ref-type="aff" rid="AFF1">1</xref><xref ref-type="corresp" rid="COR1">*</xref><pmc-comment>gardner26@llnl.gov </pmc-comment>
</contrib><contrib contrib-type="author"><name><surname>Williams</surname><given-names>Peter L.</given-names></name><xref ref-type="aff" rid="AFF1">1</xref></contrib><contrib contrib-type="author"><name><surname>Hara</surname><given-names>Christine</given-names></name><xref ref-type="aff" rid="AFF1">1</xref></contrib><contrib contrib-type="author"><name><surname>Wagner</surname><given-names>Mark C.</given-names></name><xref ref-type="aff" rid="AFF1">1</xref></contrib><contrib contrib-type="author"><name><surname>Colston</surname><given-names>Bill W.</given-names><suffix>Jr</suffix></name><xref ref-type="aff" rid="AFF1">1</xref></contrib></contrib-group><aff id="AFF1"><label><sup>1</sup></label> Computations/Global Security,<sup>2</sup> Physical & Life Sciences, Lawrence Livermore National Laboratory and<sup>3</sup> QuantaLife, Inc., Livermore, CA, USA</aff><author-notes><corresp id="COR1">*To whom correspondence should be addressed. Tel: +1 925 422 4317; Fax: +1 925 423 6437; Email:
<email>gardner26@llnl.gov</email></corresp><corresp>Correspondence may also be addressed to Amy L. Hiddessen. Tel:
<phone>+1 925 422 4787</phone>; Fax:
<fax>+1 925 423 8920</fax>; Email:
<email>hiddessen1@llnl.gov</email></corresp></author-notes><pub-date pub-type="ppub"><day>01</day><month>10</month><year>2009</year></pub-date><pub-date iso-8601-date="2009-09-16" pub-type="epub"><day>16</day><month>9</month><year>2009</year></pub-date><pub-date pub-type="pmc-release"><day>16</day><month>9</month><year>2009</year></pub-date><pmc-comment> PMC Release delay is 0 months and 0 days and was based on the
<volume>37</volume><issue>19</issue><fpage>6291</fpage><lpage>6304</lpage><history><date date-type="received"><day>18</day><month>12</month><year>2008</year></date><date date-type="rev-recd"><day>22</day><month>7</month><year>2009</year></date><date date-type="accepted"><day>24</day><month>7</month><year>2009</year></date></history><permissions><copyright-statement>© The Author(s) 2009. Published by Oxford University Press.</copyright-statement><copyright-year>2009</copyright-year><license license-type="cc-by-nc" xlink:href="http://creativecommons.org/licenses/by-nc/2.5/uk/"><license-p><pmc-comment>CREATIVE COMMONS</pmc-comment>This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/2.5/uk/">http://creativecommons.org/licenses/by-nc/2.5/uk/</ext-link>) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p></license></permissions><self-uri xlink:href="gkp659.pdf"></self-uri><abstract><title>Abstract</title><p>We describe a Multiplex Primer Prediction (MPP) algorithm to build multiplex compatible primer sets to amplify all members of large, diverse and unalignable sets of target sequences. The MPP algorithm is scalable to larger target sets than other available software, and it does not require a multiple sequence alignment. We applied it to questions in viral detection, and demonstrated that there are no universally conserved priming sequences among viruses and that it could require an unfeasibly large number of primers (∼3700 18-mers or ∼2000 10-mers) to generate amplicons from all sequenced viruses. We then designed primer sets separately for each viral family, and for several diverse species such as foot-and-mouth disease virus (FMDV), hemagglutinin (HA) and neuraminidase (NA) segments of influenza A virus, Norwalk virus, and HIV-1. We empirically demonstrated the application of the software with a multiplex set of 16 short (10 nt) primers designed to amplify the Poxviridae family to produce a specific amplicon from vaccinia virus.</p></abstract></article-meta></front><body><sec sec-type="intro"><title>INTRODUCTION</title><p>Researchers employ numerous approaches for viral detection and discovery, including metagenomic sequencing (
<xref rid="B1" ref-type="bibr">1–3</xref>), microarrays (
<xref rid="B4" ref-type="bibr">4–9</xref>) or multiplex PCR followed by other methods of characterization such as mass spectrometry (
<xref rid="B10" ref-type="bibr">10–14</xref>), suspension arrays (
<xref rid="B15" ref-type="bibr">15</xref>,
<xref rid="B16" ref-type="bibr">16</xref>) or amplicon sequencing (
<xref rid="B17" ref-type="bibr">17</xref>). Multiplex PCR followed by analysis of amplified products fills a niche for viral identification when singleplex PCR has failed or there are a few dozen likely candidates but the expense of metagenomic sequencing or high-density microarrays is unwarranted (
<xref rid="B18" ref-type="bibr">18</xref>). However, multiplex primer design for many highly divergent targets is challenging since no universally conserved primers may exist, and finding sets of primers likely to function well in multiplex (e.g. isothermal
<italic>T</italic><sub>m</sub>’s, no primer dimers) adds to the complexity of finding conserved primer candidates. Primer design software that requires a multiple sequence alignment (MSA) as input can be problematic for diverse target sets, as MSAs can be difficult to construct, exhausting memory or available time before an alignment is completed. Even if an alignment does complete for divergent target sets such as all members of a family of RNA viruses or gene homologues across species, alignments may show little nucleotide sequence conservation, and multiple primers are required to amplify all targets. Considering the challenges of primer design for targets showing sequence variation, it is not surprising that many of the PCR-based assays in the literature are predicted to fail to detect desired targets when compared against available sequence data, and this problem is worst at higher taxonomic levels like family (
<xref rid="B19" ref-type="bibr">19</xref>).
</p><p>Most currently available multiplex primer prediction tools require an MSA (
<xref rid="B20" ref-type="bibr">20–24</xref>). None of those tools build multiplex sets in which no primers in the set are predicted to form primer-dimers, although some avoid homodimers. We attempted to run a number of alternative software tools for multiplex or degenerate primer prediction for several species level target sets ranging in size from 41 to ∼6000 sequences: Primaclade (
<xref rid="B24" ref-type="bibr">24</xref>), FastPCR (
<ext-link ext-link-type="uri" xlink:href="www.biocenter.helsinki.fi/bi/programs/fastpcr.htm">www.biocenter.helsinki.fi/bi/programs/fastpcr.htm</ext-link>), GeneUp (
<xref rid="B25" ref-type="bibr">25</xref>), PDA-MS/UniQ software (
<xref rid="B26" ref-type="bibr">26</xref>), Greene SCPrimer (
<xref rid="B20" ref-type="bibr">20</xref>) and HYDEN (
<xref rid="B22" ref-type="bibr">22</xref>). Only Greene SCPrimer and HYDEN could handle the two smallest, species-level target sets (Norwalk virus and FMDV), and none of the tools completed for the larger target sets that we attempted to run.
</p><p>Therefore, we designed the MPP software to avoid the requirement of an MSA and to scale better for large and diverse target sets. MPP builds a multiplex-compatible set of primers capable of amplifying all target sequences, attempting to minimize the number primers in the set. We set out to determine a set of highly conserved ‘universal’ primers for viruses, akin to the highly conserved 16S rRNA universal primers for bacteria. Throughout, we use the term ‘detected’ to mean that there should be at least one PCR product. This is a loose definition of ‘detected’ adopted for convenience in the following discussions, and we recognize that a PCR product may prove insufficient for viral characterization. We predicted a set of ‘universal viral primers’ for all available complete genomes of all viruses, and found that the number of universal viral primers would be impractical to implement, even if short, highly conserved priming sequences were used. Then we predicted family-level primer sets for every viral family, as well as for several highly diverse species of RNA viruses, for primers of a traditional length as well as nontraditional, shorter, more highly conserved primers, which are more likely to amplify novel, unsequenced viruses and which could be an alternative to degenerate primers. While the software uses a greedy algorithm that may settle at a local minimum which is above the true minimal set, it generated primer sets for each viral family with fewer than half the number of primers that would be expected without optimization. Although family-level primer sets for some families are too large to be practical, 64% of the families had primer sets of no more than 20 primer pairs, quite feasible for multiplex PCR.</p><p>Finally, we empirically multiplexed a set of 16 short (10-mer) primers designed for the Poxviridae family. This demonstrated specific amplification of the expected viral fragment from vaccinia Lister strain. This preliminary demonstration suggests that specific amplification with family-level multiplexed sets of short primers might be feasible for viral detection and discovery, particularly as a means for selective enrichment of viral target for downstream amplicon characterization (e.g. sequencing or probe hybridization).</p></sec><sec sec-type="methods"><title>METHODS</title><sec><title>MPP algorithm for calculating multiplexed primer sets</title><p>Here, we outline a greedy algorithm used to calculate conserved sets of multiplexed primers to amplify fragments from each member of a target set of sequences, and provide a more detailed description in the
<xref ref-type="supplementary-material" rid="sup1">Supplementary Methods</xref>. First, we enumerate all candidate oligos fitting user-specified requirements for length,
<italic>T</italic><sub>m</sub>, and lack of hairpin formation. We rank pairs of these by the number of targets in which that pair occurs within a distance range
<italic>d</italic><sub>1</sub>bases to
<italic>d</italic><sub>2</sub>bases of one another, where these might be specified so as to bound a reasonable range for a downstream characterization method such electrophoretic discrimination of bands, probe hybridization or sequencing. The most frequent pair is selected as primers. The process is repeated for the remaining targets that would not have an amplicon from the first pair, with the added consideration that new primers selected be predicted not to form dimers with other primers already selected, as well as can be predicted based on nearest neighbor thermodynamic predictions (
<xref rid="B27" ref-type="bibr">27</xref>), although free energy calculations cannot predict with certainty that primer dimers will be excluded in practice. There is an option to bin primers into reaction subsets, if desired. With binning, primers are added to a bin until that bin contains b primers, at which point a new bin is begun, following the same process. Binning primers in smaller groups avoids exclusion of the most highly conserved oligos because of primer dimer free energy constraints. The universal set of primers is the set of selected primers to amplify all genomes in the target set. The primers within a bin should be multiplexed into a single PCR reaction, but each bin should be run separately. This is a simple binning strategy, and alternative strategies could be employed such as starting a new bin with any primer pair that dimerizes with other previously selected primers regardless of the number of primers in the bin, but this could have the possible disadvantage of bin explosion to numerous singleplex or small multiplex reactions. The graph-based algorithm of MuPlex (
<xref rid="B28" ref-type="bibr">28</xref>,
<xref rid="B29" ref-type="bibr">29</xref>) is another binning strategy that could be incorporated. Output on the number of primers rejected due to the various filters (
<italic>T</italic><sub>m</sub>, hairpin free energy, etc.) is printed to standard output after each round of primer pairs are selected and also cumulatively, so a user can monitor if/which parameters might be too stringent. If an alternative set of primers is desired that does not overlap with the set selected, one can replace with N’s the subsequences matching the selected primers and their reverse complements in the target sequences and rerun the software with the modified input sequences.
</p><p>The script find_amplicons.pl predicts all the amplicons that should be generated by a list of (multiplexed) primers mixed with a set of sequences. This PCR-simulating script lists amplicon sequences, their length and position, and the forward+reverse primer combination to yield each product, as well as a summary file of the fragment length distributions enabling a quick assessment of how well each target sequence can be discriminated based purely on amplicon length patterns. We used this script to check whether some viral primer sets are predicted to generate amplicons from the human genome.</p><p>For all
<italic>T</italic><sub>m</sub>and free energy predictions, we used the following Unafold settings: [Na+] = 0.2 M, [Mg+] = 0.0015, annealing temperature = 30°C, and strand concentration of each strand = 10
<sup>–7</sup>M, making the total strand concentration of both strands = 2 × 10
<sup>–7</sup>.
</p></sec><sec><title>Imposing uniqueness requirements for target-specific priming</title><p>The MPP algorithm described here focuses on finding conserved primers, and does not require that the primers be family- or species-specific. For the runs predicting family-specific primers, we designed viral family primers that were unique relative to nontarget viral families by replacing any substring of 17 nt or longer that occurs in any non-target family with a substring of ‘N’’s, using the suffix array software vmatch (
<ext-link ext-link-type="uri" xlink:href="http://www.vmatch.de/">http://www.vmatch.de/</ext-link>). MPP eliminates oligos with N’s from consideration as primers. This approach is simple, although it risks being overly strict by eliminating some potentially successful candidate primers, since it disallows those cases where a single nonunique primer pairs with a unique primer, a pair of nonunique primers are too far apart to actually generate an amplicon in a nontarget sequence, or candidate oligos partially overlap a nonunique stretch of N’s.
</p><p>In general, searching for primers from sequence that is specific to one set of targets and excluding candidate substrings present in nontarget sequence is a useful strategy to design signatures for pathogen detection. This can be done by replacing nonunique or nonspecific substrings with N’s in the sequences input to MPP using software such as vmatch.</p></sec><sec><title>Runs predicting universal viral primers</title><p>For runs predicting universal viral primers including all viruses at once in the target set, all viral complete genomes and segments downloaded from publicly available sequence databases (Genbank, Baylor, TIGR) as of 25 April 2007 were used. Draft sequences with multiple contigs were merged into a single sequence entry, with contigs separated by 1000 N’s, a stretch sufficiently long (>
<italic>d</italic><sub>2</sub>) so that primer pairs would not be designed to fall on different contigs, although there were very few draft sequences in contigs where this was necessary. Because of the large numbers of sequences in two families, only the MP segment sequences from Orthomyxoviridae and the L segment from Bunyaviridae were included, as these are the more conserved segments, reducing the number of targets by 23 017 sequences. The total number of target sequences for these all virus runs were 11 477. We predicted a multiplex set of universal 10-mers without binning the primers into separate reactions, but it was necessary to use a very low
<italic>x</italic><sub>dimer</sub>and
<italic>x</italic><sub>homodimer</sub>of −11 kcal/mol to be possible to predict multiplexed primers to amplify every target. For the next calculations, we subdivided the primers into sub-reactions with 20 primers per reaction bin, to avoid excluding primers that were predicted to form primer dimers, and raised
<italic>x</italic><sub>dimer</sub>and
<italic>x</italic><sub>homodimer</sub>to –7 kcal/mol. Primer sets of length 5–18 were predicted (
<xref ref-type="fig" rid="F1">Figure 1</xref>). To assess the effect of removing
<italic>T</italic><sub>m</sub>constraints, we generated universal primer sets with length but no
<italic>T</italic><sub>m</sub>requirements, and compared the primer counts to those with
<italic>T</italic><sub>m</sub>constraints (
<xref ref-type="fig" rid="F2">Figure 2</xref>). Primer sets are available by contacting the authors.
</p><fig id="F1" orientation="portrait" position="float"><label>Figure 1.</label><caption><p>Hundreds to thousands of primer sequences are required to amplify all viruses. Percent of viral genomes detected versus number of primers required. Calculated for primers of different sizes, and based on sequence data as of 25 April 2007. Parameters used are given in
<xref ref-type="supplementary-material" rid="sup1">Supplementary Tables S1 and S2</xref>.
</p></caption><graphic xlink:href="gkp659f1"></graphic></fig><fig id="F2" orientation="portrait" position="float"><label>Figure 2.</label><caption><p>Effect of
<italic>T</italic><sub>m</sub>constraints and sequence database size on size of universal primer set. The number of primers in the universal set for all viruses is shown, as a function of primer size. Calculations were based either on sequence data available April 2007 and imposing
<italic>T</italic><sub>m</sub>constraints in primer selection, or without
<italic>T</italic><sub>m</sub>or GC% constraints, or based on January 2004 sequence data with
<italic>T</italic><sub>m</sub>constraints. For primer size of 5 with
<italic>T</italic><sub>m</sub>constraints, a constraint of
<italic>T</italic><sub>m</sub>> 0 was used, delivering primers with GC% of 80–100%.
</p></caption><graphic xlink:href="gkp659f2"></graphic></fig><p>We evaluated how the growth of sequence availability could affect the size of the universal set of viral primers, as well as our ability to detect unknown/unsequenced viruses. We predicted a universal viral primer set for all viral genomes and segments available as of 1 January 2004, totaling 9965 sequences, for primers of length 7–15 nt (
<xref ref-type="fig" rid="F2">Figure 2</xref>). It was not necessary to exclude any Orthomyxoviridae or Bunyaviridae segments, because these segments were not so deeply sequenced at that time. We then determined how much of the 2007 sequence data would have been detectable using the 2004 primer sets (
<xref ref-type="fig" rid="F3">Figure 3</xref>).
</p><fig id="F3" orientation="portrait" position="float"><label>Figure 3.</label><caption><p>Primers from old sequence data miss many newly sequenced viruses. Fraction of viral genomes available as of 25 April 2007 that would have been detected using primer sets developed based on the sequence data available as of 1 January 2004.</p></caption><graphic xlink:href="gkp659f3"></graphic></fig><p>To predict how contamination with human DNA might affect the ability to detect specific amplification of viruses, the average number of amplicons for the human genome was also predicted based on these primer sets with 20 primers per bin (
<xref ref-type="fig" rid="F4">Figure 4</xref>). The effect of reducing the multiplex size to 10 primers per reaction on human genome amplification was calculated by dividing priming sets in half, with 10 primers per reaction.
</p><fig id="F4" orientation="portrait" position="float"><label>Figure 4.</label><caption><p>Average number of bands predicted from PCR against the human genome using the universal viral primers for the runs as described in
<xref ref-type="fig" rid="F2">Figure 2</xref>, average ± SD of the top 10 bins for each primer length. For bin size of 10, we divided the universal primer bins of 20 primers/bin into two bins each.
</p></caption><graphic xlink:href="gkp659f4"></graphic></fig></sec><sec><title>Viral family primers</title><p>Subdividing viral targets into families, we used an updated set of sequences, downloaded 5 February 2009. For family-level primers, we computed primer sets using three alternative parameter settings (
<xref rid="T1" ref-type="table">Table 1</xref>): (i) 17–21-mer primers with
<italic>T</italic><sub>m</sub>= 55–60°C, primers not checked for uniqueness; (ii) same length and
<italic>T</italic><sub>m</sub>specifications as in (i) but eliminating from consideration as primers any oligos of at least 17 nt that occur in nontarget viral families, as described earlier; and (iii) 10–15-mer primers with
<italic>T</italic><sub>m</sub>= 40–45°C. These primers are available as
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>.
</p><table-wrap id="T1" orientation="portrait" position="float"><label>Table 1.</label><caption><p>Primer counts for viral family primer sets, as described in the methods</p></caption><table frame="hsides" rules="groups"><tbody align="left"><tr><td colspan="1" rowspan="1"><inline-graphic xlink:href="gkp659i1.jpg"></inline-graphic></td></tr></tbody></table><table-wrap-foot><fn><p>Grey shaded entries indicate where calculations were not run to completion. In other cases (not shaded) where fewer than 100% of targets were predicted to be amplified, the algorithm failed to find primer pairs that met all the required specifications for the remaining targets, that is, primers in the right length,
<italic>T</italic><sub>m</sub>, and amplicon length range, with hairpin and dimer avoidance with other primers already selected to be in the set, could not be found.
</p></fn><fn><p>For 10–15-mers,
<italic>T</italic><sub>m</sub>= 40–45°C and for 17–21-mers,
<italic>T</italic><sub>m</sub>= 55–60°C.
</p></fn><fn><p>For both, the amplicon length range was 200–800 bp.</p></fn><fn><p>None of the calculations that completed required more than the 16 GB of RAM that was available.</p></fn></table-wrap-foot></table-wrap></sec><sec><title>Species level primers for several divergent species</title><p>We also generated primer sets for several species with high sequence diversity: HIV-1, FMDV, Norwalk virus and influenza A segments HA and NA (summary in
<xref rid="T2" ref-type="table">Tables 2</xref>and
<xref rid="T3" ref-type="table">3</xref>, primer sequences available as
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>). MPP parameters were the same as in
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>, with
<italic>T</italic><sub>m</sub>= 35–50°C for 10-mer primers and
<italic>T</italic><sub>m</sub>= 55–70°C for 17–18-mer primers, and sequences were downloaded in 2007 or 2008 (influenza only). To illustrate the challenges of designing primers from an alignment, we aligned these organisms using MUSCLE (
<xref rid="B30" ref-type="bibr">30</xref>) when possible. For the HIV-1 and influenza A segments HA and NA, MUSCLE ran out of memory before completing. An alternative alignment tool, ClustalW (
<xref rid="B31" ref-type="bibr">31</xref>), had completed only a small fraction of the alignment after running for days. Therefore, for these large data sets, a random selection of ∼35 sequences for each target was aligned with MUSCLE, and this alignment was used to build a profile Hidden Markov Model (HMM) (hmmbuild) using HMMer (
<ext-link ext-link-type="uri" xlink:href="http://hmmer.wustl.edu/">http://hmmer.wustl.edu/</ext-link>).(
<xref rid="B32" ref-type="bibr">32</xref>) The full sequence set was then aligned to the HMM using hmmalign. For Norwalk virus and FMDV, we designed multiplex-degenerate primer sets using GreeneSCPrimer (
<xref rid="B20" ref-type="bibr">20</xref>) and HYDEN (
<xref rid="B22" ref-type="bibr">22</xref>), but the other three targets sets were too large for those programs. None of the other primer prediction programs we tried (as indicated in the ‘Results’ section) would scale to handle even these two smallest target sets.
</p><table-wrap id="T2" orientation="portrait" position="float"><label>Table 2.</label><caption><p>Summary of results of primer set prediction for several diverse RNA virus species for primers of length 10 or 17–18 nt, the longest stretch of bases conserved among the targets in a multiple sequence alignment, illustrating the difficulty of selecting conserved primers, and the time required to perform the primer set selections</p></caption><table frame="hsides" rules="groups"><thead align="left"><tr><th colspan="1" rowspan="1">Virus species</th><th colspan="1" rowspan="1">Number of sequences</th><th colspan="1" rowspan="1">Number of 10-mers</th><th colspan="1" rowspan="1">Number of 17–18-mer primers in set</th><th colspan="1" rowspan="1">Longest conserved region from MSA (nt)</th><th colspan="1" rowspan="1">Processing time on 1 cpu of an AMD Opteron to calculate 10-mer primers (hours)</th><th colspan="1" rowspan="1">Processing time on 1 cpu of an AMD Opteron to calculate 17–18-mer primers (hours)</th></tr></thead><tbody align="left"><tr><td colspan="1" rowspan="1">NA segment of Influenza A
<xref ref-type="table-fn" rid="TF5">*</xref></td><td colspan="1" rowspan="1">6375</td><td colspan="1" rowspan="1">52</td><td colspan="1" rowspan="1">120</td><td colspan="1" rowspan="1">5</td><td colspan="1" rowspan="1">5.7</td><td colspan="1" rowspan="1">28.9</td></tr><tr><td colspan="1" rowspan="1">HA segment of Influenza A
<xref ref-type="table-fn" rid="TF5">*</xref></td><td colspan="1" rowspan="1">5440</td><td colspan="1" rowspan="1">73</td><td colspan="1" rowspan="1">153</td><td colspan="1" rowspan="1">1</td><td colspan="1" rowspan="1">6</td><td colspan="1" rowspan="1">29.5</td></tr><tr><td colspan="1" rowspan="1">HIV-1</td><td colspan="1" rowspan="1">1175</td><td colspan="1" rowspan="1">6</td><td colspan="1" rowspan="1">16</td><td colspan="1" rowspan="1">0</td><td colspan="1" rowspan="1">1.6</td><td colspan="1" rowspan="1">9.1</td></tr><tr><td colspan="1" rowspan="1">FMDV</td><td colspan="1" rowspan="1">187</td><td colspan="1" rowspan="1">4</td><td colspan="1" rowspan="1">6</td><td colspan="1" rowspan="1">9</td><td colspan="1" rowspan="1">0.18</td><td colspan="1" rowspan="1">0.7</td></tr><tr><td colspan="1" rowspan="1">Norwalk</td><td colspan="1" rowspan="1">41</td><td colspan="1" rowspan="1">7</td><td colspan="1" rowspan="1">20</td><td colspan="1" rowspan="1">6</td><td colspan="1" rowspan="1">0.1</td><td colspan="1" rowspan="1">0.5</td></tr></tbody></table><table-wrap-foot><fn id="TF5"><p>*For Influenza A HA and NA segments, all complete sequences, including lab strains, from all hosts, countries and serotypes were downloaded from the NCBI Influenza Virus Resource database (
<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/genomes/FLU/Database/multiple.cgi">http://www.ncbi.nlm.nih.gov/genomes/FLU/Database/multiple.cgi</ext-link>) on 18 January 2008.
</p></fn><fn id="TF6"><p>Primer sequences are provided as
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>.
</p></fn></table-wrap-foot></table-wrap><table-wrap id="T3" orientation="portrait" position="float"><label>Table 3.</label><caption><p>Summary of the number of 17–18-mer primers predicted to amplify several influenza A HA serotypes</p></caption><table frame="hsides" rules="groups"><thead align="left"><tr><th colspan="1" rowspan="1">Influenza A HA Serotype</th><th colspan="1" rowspan="1">Number of sequences</th><th colspan="1" rowspan="1">Number of 17–18-mer primers in set</th></tr></thead><tbody align="left"><tr><td colspan="1" rowspan="1">H1</td><td colspan="1" rowspan="1">1080</td><td colspan="1" rowspan="1">24</td></tr><tr><td colspan="1" rowspan="1">H2</td><td colspan="1" rowspan="1">108</td><td colspan="1" rowspan="1">8</td></tr><tr><td colspan="1" rowspan="1">H3</td><td colspan="1" rowspan="1">1972</td><td colspan="1" rowspan="1">15</td></tr><tr><td colspan="1" rowspan="1">H5</td><td colspan="1" rowspan="1">1325</td><td colspan="1" rowspan="1">16</td></tr><tr><td colspan="1" rowspan="1">H7</td><td colspan="1" rowspan="1">256</td><td colspan="1" rowspan="1">8</td></tr></tbody></table><table-wrap-foot><fn id="TF7"><p>Primer sequences are provided as
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>.
</p></fn></table-wrap-foot></table-wrap></sec><sec><title>Empirical testing of Poxviridae multiplex primers</title><p>For an empirical demonstration of our algorithm, we tested a multiplex set of 16 short, 10-nt primers designed for the Poxviridae viral family against commercially purified vaccinia virus extracts. The primer sequences are provided in
<xref rid="T4" ref-type="table">Table 4</xref>along with the predicted single amplicon for vaccinia Lister strain. We chose the Poxviridae family multiplex for these first empirical demonstrations because extracted viral nucleic acid was readily available for experiments. All primers were purchased from Integrated DNA Technologies (Coralville, IA, USA) and resuspended to 100 µM stock solutions in TE buffer (pH 8.0, Teknova, Hollister, CA, USA). Working solutions containing equimolar concentrations of each of the 16 primers were used in all experiments. Purified, quantitated vaccinia Lister strain DNA was purchased at a concentration of 1.3 × 10
<sup>4</sup>copies/µl in nuclease-free water from Advanced Biotechnologies, Inc. (Columbia, MD, USA). All PCR experiments were prepared using the Superscript III RT-PCR kit from Invitrogen (Carlsbad, CA, USA). We selected the RT-PCR kit in order to establish a protocol that could later be readily applied to additional multiplex viral family reactions with viral DNA and/or RNA. Each 25 µl reaction contained 1 × SSIII buffer, 1 U of SSIII RT/Taq enzyme, 4.8 mM MgSO
<sub>4</sub>, 0.1 µM each primer, and a viral template mass of 2.7 pg (∼10
<sup>4</sup>copies). Tests were performed in triplicate and corresponding negative controls were run under identical conditions except that viral template was replaced with nuclease-free water (Ambion, Austin, TX, USA). All reactions were thermocycled on the Bio-Rad DNA Engine (Hercules, CA, USA) as follows: one cycle of 2 min at 94°C; 40 cycles of 15 s at 94°C, 30 s at 43.9°C, 1 min at 68°C; one cycle of 5 min at 68°C. In line with our goal to create a protocol applicable to both DNA and/or RNA templates, we verified that inclusion of an RT step (one cycle of 30 min at 45°C) does not alter the outcome of the subsequent PCR when working with DNA template (data not shown). As such, for all Poxviridae 16-plex results reported here, an RT step was not used. Another goal was to establish reaction conditions that can be applied to any multiplex viral primer set without the need for re-optimization, assuming the primer set is designed with the same general parameter contraints (
<italic>T</italic><sub>m</sub>s, etc.) used to compute the 16-plex presented here. As such, we first optimized master mix conditions, where it was determined that a 4.8 mM MgSO
<sub>4</sub>concentration resulted in most optimal amplification. Similarly, we determined the optimal annealing temperature based on results from annealing temperature gradient experiments. A more detailed discussion on the optimization of reaction conditions for amplification with our short primer viral multiplexes is the subject of a separate paper in preparation (Hiddessen and co-workers).
</p><table-wrap id="T4" orientation="portrait" position="float"><label>Table 4.</label><caption><p>Sequence and size (nt) information for the 16 primers in the predicted family-level multiplex for Poxviridae viruses used in these experiments</p></caption><table frame="hsides" rules="groups"><thead align="left"><tr><th colspan="1" rowspan="1">Primer sequence</th><th colspan="1" rowspan="1">Primer size (nt)</th></tr></thead><tbody align="left"><tr><td colspan="1" rowspan="1"><bold>CGGAGACCAA</bold>(FP 617)
</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1"><bold>TCGTCGTCCA</bold>(RP 617)
</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TGTTGGTGTG</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TCGTCTACGA</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TGTGGTCCTT</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">CCATGTTCGC</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TCGAGGAGAA</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">CCGGCTCCAG</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TGAACCTGGT</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">GCGCACGTAC</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TCACGCATCT</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">GGGAAACAGC</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">ACCGTTGTCA</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">TACCATCGTC</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">AACCTGTGCA</td><td colspan="1" rowspan="1">10</td></tr><tr><td colspan="1" rowspan="1">GGCGGAGGTA</td><td colspan="1" rowspan="1">10</td></tr></tbody></table><table-wrap-foot><fn id="TF8"><p>The forward and reverse primers (FP, RP) in bold are predicted to amplify vaccinia Lister with the indicated product length (617-bp amplicon size).</p></fn></table-wrap-foot></table-wrap><p>This multiplex was compared against the human genome, predicting 233 amplicons between 50 and 1000 bp from the Poxviridae 16-plex. For empirical tests against background human nucleic acids, we followed the same reaction conditions as above. In these tests, vaccinia DNA was held at a constant concentration of 2.7 pg (1.3 × 10
<sup>4</sup>copies), and serial dilutions of human genomic DNA (Novagen, Madison, WI, USA), were added to the vaccinia PCR reaction mix, starting at 2.7 pg and titrating down over 4 orders of magnitude to 2.7 × 10
<sup>–4</sup>pg. These experiments were performed using mass ratios of vaccinia:human DNA at 1 : 1, 10 : 1, 100 : 1, 1000 : 1 and 10000 : 1. These correspond to copy number ratios of vaccinia:human genomes of 1.3 × 10
<sup>4</sup>: 0.82, 1.3 × 10
<sup>4</sup>: 0.082, 1.3 × 10
<sup>4</sup>: 0.0082, 1.3 × 10
<sup>4</sup>: 0.00082 and 1.3 × 10
<sup>4</sup>: 0.000082, respectively.
</p><p>All PCR experiments were analyzed on 3% agarose TBE gels containing ethidum bromide that were purchased from Bio-Rad (Hercules, CA, USA). Blue juice™ 10× loading dye was purchased from Invitrogen (Carlsbad, CA, USA) and diluted to 2× before use. A 50-bp DNA ladder was purchased from Novagen (Madison, WI, USA). For analysis, 15 μl from each separate 25 μl PCR reaction were combined with 2 μl of of loading dye and 15 μl of the loading-dye/product mixture was loaded per well and electrophoresed for 1 h 40 min at 85 V.</p><p>For confirmation of amplicon sequence, 9 μl of amplified PCR product was mixed with 4 μl of ExoSap-it (USB, Cleveland, OH, USA) and incubated at 37°C for 15 min followed by a denaturation step at 80°C for 15 min. Sanger sequencing was then performed with the BigDye V3.1 Terminator Kit (Applied Biosystems, Inc., Foster City, CA, USA). Single reactions contained 4 μl of Ready Reaction Mix, 0.2-µM primer, 2 μl 5× Sequencing Buffer, 11.5 µl of nuclease free water (Applied Biosystems, Inc.) and 2 μl of post-Exosap-it PCR product. The sequencing reaction used the following thermocycling profile: 1 cycle of 94°C for 1 min; 25 cycles of 94°C for 15 s, 38°C for 30 s and 60°C for 4 min. Two microliters of sequencing reaction product was combined with 18 μl of Hi-Di
<sup>Tm</sup>Formamide (Applied Biosystems, Inc.) and run on the ABI 3130 (Applied Biosystems, Inc.). The results were analyzed using Sequencing Analysis v5.2 (Applied Biosystems, Inc.).
</p></sec></sec><sec sec-type="results"><title>RESULTS</title><sec><title>Universal viral primers</title><p>Predicting a set of universal viral primers that are all mixed in a single large reaction, we predict that 1008 primers of length 10 nt would be required to ensure amplification of at least one fragment of length 80–620 bp from every viral genome. These are predicted to generate a mean and maximum number of amplicons per viral genome of 13.2 and 948, respectively. Binning the primers into smaller sets of 20 primers/bin doubled the number of primers required, as
<xref ref-type="fig" rid="F1">Figure 1</xref>shows the fraction of viral genomes amplified versus the number of primers, for primers ranging from 6 to 18 nt in length. About 2000 binned 10-mers are required to amplify 100% of sequenced viruses, generating a mean of 1 and maximum of 2.9 amplicons per genome, on average. Using traditional-length 18-mer primers nearly doubles the number to ∼3700 primers. The incomplete curve for 7-mers was from inappropriate settings for this oligo length, since some genomes do not contain pairs of 7-mers with the required
<italic>T</italic><sub>m</sub>, as given in the
<xref ref-type="supplementary-material" rid="sup1">Supplementary Tables S1 and S2</xref>, within the desired amplicon length range. The concave curves showing diminishing returns reflect biased sequence availability rather than any particularly highly conserved primers: the primers that amplify the largest number of sequences are all from influenza, as far more sequences are available for this species than for any other.
</p></sec><sec><title>Increasing availability of sequence data on universal primer sets</title><p>The increase in sequence data requires ∼700 more 10-mer primers to amplify all sequenced viruses in 2007 compared to 2004 (
<xref ref-type="fig" rid="F2">Figure 2</xref>). While the increase in the number of sequences used between the two dates was only ∼15%, the number of primers required increased by 48%, illustrating the substantial increase in diversity represented by the additional sequence data.
<xref ref-type="fig" rid="F2">Figure 2</xref>also shows that removing all
<italic>T</italic><sub>m</sub>constraints allows fewer primers to be used since no conserved primers are eliminated due to
<italic>T</italic><sub>m</sub>, as some AT rich subsequences tend to be fairly conserved.
<xref ref-type="fig" rid="F3">Figure 3</xref>shows that a universal primer set predicted using the 2004 sequence data would amplify only 35% of the 2007 sequences using primers of at least 10 nt in length. Shorter primers increase this fraction to over 60%, due to the higher likelihood of occurrence and conservation of shorter oligos, but even so, a multiplex of 7-mers is not guaranteed to amplify a fragment from every virus. The minor differences in the number of genomes detected between 12, 13, 14 and 15-mers can be attributed to the facts that a greedy but not necessarily optimal algorithm is used to select one solution from among many, that the primers in a particular set depend on the
<italic>T</italic><sub>m</sub>ranges we used as well as length differences, and the unpredictable nature of novel viral sequences accumulated between 2004 and 2007, rather than any real difference in the ability to detect genomes among those primer sets.
</p></sec><sec><title>Predicting effects of contamination from human DNA</title><p>If all human DNA cannot be removed from the sample, simulations indicate that on average, multiplexes of 10-mer primers are expected to produce hundreds of amplicons from the human genome, which would appear as a smear on a gel (
<xref ref-type="fig" rid="F4">Figure 4</xref>). With primers of length 11–14, there is a large variation among bins. While most short sequences do occur in the human genome (
<xref rid="B33" ref-type="bibr">33</xref>), an amplicon requires that two occur in proximity, and primers ≥11 bases make this a sporadic event for bins with only 10 primers. Nonetheless, imperfect sample purification to eliminate eukaryotic nucleic acids could be problematic for universal viral priming using primers shorter than 15 bases, particularly for multiplexes of 10 or more primers.
</p></sec><sec><title>Family identification by sequencing products from universal 10-mer primer set versus randomly amplified fragments</title><p>If universal viral primers amplified fragments from a newly emerged virus, would the product sequences show similarity to others in the same family? We predicted amplicon sequences from 10 newly emerged viruses (Australian bat lyssavirus, Crimean Congo hemorrhagic fever, Nipah, Hendra, Venezuelan equine encephalitis, West Nile, Ebola, Hanta, hepatitis C and Marburg viruses) using the universal 10-mer primer set and BLASTed (blastn) (
<xref rid="B34" ref-type="bibr">34</xref>) those amplicon sequences against either other species in the same target family or viruses in other families. On average, 87% of the amplicons had BLAST hits in the correct family, with an average of 135 hits per amplicon. In contrast, only 11% of amplicons had hits in other, nontarget families, with an average of 2.5 hits per amplicon. For comparison, 72% of randomly selected fragments of 200 bp (or 65% and 78% for 100-bp and 600-bp fragments, respectively) had hits in the correct family with an average of 59 hits per amplicon, and 2.8% had hits in other families with an average of 2.4 hits per amplicon. Thus, universal primers do amplify more conserved regions than randomly selected fragments. These rough calculations depend heavily on available sequence data, and are likely to significantly overestimate our ability to characterize a truly novel virus, extrapolating from results of viral metagenomic sequencing in which over 60% of the sequences cannot be identified based on BLAST comparisons to existing sequence databases (
<xref rid="B2" ref-type="bibr">2</xref>).
</p></sec><sec><title>Viral family primers</title><p>The number of family-level primers for each family, and the number of genomes available for generating those primer sets, is given in
<xref rid="T1" ref-type="table">Table 1</xref>, for three alternative settings: short primers of 10–15 nt with
<italic>T</italic><sub>m</sub>40–45°C, standard length primers of 17–21 nt and
<italic>T</italic><sub>m</sub>55–60°C, and the same but requiring that each primer subsequence of at least 17 nt be unique to the target viral family relative to other viral families. At least one amplicon of 200–800 bp was required from every genome. Hypothetically, the worst-case scenario to amplify a target set of N sequences would require 2N primers. MPP requires on average only 37% or 45% of this number, for primers of length 10–15 nt or 17–21 nt without the requirement for family specific primers, respectively (averaged across families, omitting those families for which the computations were not run to completion). The most diverse families, in particular Bunyaviridae, Geminiviridae, Polydnaviridae, Reoviridae, Retroviridae, Siphoviridae and Orthomyxoviridae require so many primers that actually applying family-level amplification is probably infeasible. For these, more restricted target sets may be necessary, such as limiting to a single segment for the segmented families or to subclades, and possibly the incorporation of primers with degenerate or inosine bases. Some families with many genomes can be amplified with relatively few primers, such as Coronaviridae, Hepadnaviridae, Poxviridae, Togaviridae, Microviridae and Polyomaviridae. Using primers of length 10–15-mers or 17–21-mers, 66% or 63%, respectively, of the viral families have primer sets of 40 or fewer primers (20 primer pairs), which is feasible for typical multiplexes. The sizes of the family-level primer sets show a trend for an increase with the number of available sequences in the family (
<xref ref-type="fig" rid="F5">Figure 5</xref>,
<italic>P</italic>= 0.14). There is no clear relationship between the number of primers in the set and whether the genomes are single or double stranded RNA or DNA. All family primer sequences are available as
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>.
</p><fig id="F5" orientation="portrait" position="float"><label>Figure 5.</label><caption><p>Number of primers required to amplify fragments from each complete genome or segment in the family versus the number of available sequences in that family. Primer parameters were length 17–21-mers,
<italic>T</italic><sub>m</sub>between 55°C and 60°C.
</p></caption><graphic xlink:href="gkp659f5"></graphic></fig></sec><sec><title>Species-level primers for divergent RNA viruses</title><p>Primer design with the MPP software indicates that relatively few primers are required to amplify all sequenced genomes of HIV-1, FMDV and Norwalk virus (
<xref rid="T2" ref-type="table">Table 2</xref>), and these can be calculated in minutes. Influenza A HA and NA segments demand large numbers of 10-mer or 17–18-mer primers and hours to calculate, so one could break these into subgroups, possibly by serotype, as shown for several HA serotypes in
<xref rid="T3" ref-type="table">Table 3</xref>. The percentage of genomes amplified versus the number of primers used, for primers of either 10-mers or 17–18-mers, is shown in
<xref ref-type="fig" rid="F6">Figure 6</xref>. This plot shows that a large fraction of targets are amplified with only 2 primers, and the addition of subsequent primers shows diminishing returns in amplifying fewer, more divergent targets not detected by the initial, more conserved, primer pair, although the true diminishing returns depend on the extent to which available sequence data is an unbiased representation of diversity.
</p><fig id="F6" orientation="portrait" position="float"><label>Figure 6.</label><caption><p>The percentage of genomes amplified versus the number of primers used, for primers of either 10-mers or 17-18-mers, for several highly diverse virus species, (
<bold>A</bold>) HIV-1, FMDV and Norwalk; and (
<bold>B</bold>) Influenza A HA and NA segments. The two most highly conserved primers amplify a large fraction of genomes, and additional primers show diminishing returns in detecting the remaining, more divergent sequences.
</p></caption><graphic xlink:href="gkp659f6"></graphic></fig><p>The more traditional method of attempting to find primers from a MSA would be problematic, probably requiring manually designed primer multiplexes or highly degenerate primers. For HIV-1, for example, there is not a single position with 100% conservation across all sequenced isolates. Dropping the required conservation down to 95% (58 of the 1175 genomes could disagree with a consensus base at any position), there are three conserved regions of at least 18 bases, with positions relative to the consensus: ACAGGAGCAGATGATACAGTA starting at position 3665; TATGGAAAACAGATGGCAGG starting at 7347; and CTATGGCAGGAAGAAGCG starting at 9071. These regions are too far apart to be used as primers for most polymerases used in diagnostic PCR protocols, where amplicons must typically be less than 300 bases long for efficient amplification. A recently published study (
<xref rid="B35" ref-type="bibr">35</xref>) selected primers from the 5′ LTR U5 end to the Gag-Pol start (5′-TAGCAGTGGCGCCCGA-3′ and 5′-TCTCTCTCCTTCTAGCCTCCGC-3′), but a comparison against available genomic data indicates that 487 of the 1175 genomes (41%) do not contain a sequence match for this primer pair, so may fail to be amplified.
</p><p>For influenza A segment HA, the size of the longest conserved region from the 95% consensus is only 5 bases, and for segment NA, only 6 bases, insufficient for even a single primer. For FMDV and Norwalk virus, the longest 100% conserved regions are 9 and 6 bases, respectively. The MPP software makes it straightforward for a nonexpert to predict a multiplex-compatible set of primers to amplify all targets, even for enormous and heterogeneous target sets that cannot be aligned.</p></sec><sec><title>Comparison with other software</title><p>For comparison, we considered other software options for designing primers for these heterogeneous viruses. FastPCR (
<ext-link ext-link-type="uri" xlink:href="www.biocenter.helsinki.fi/bi/programs/fastpcr.htm">www.biocenter.helsinki.fi/bi/programs/fastpcr.htm</ext-link>) and GeneUp (
<xref rid="B25" ref-type="bibr">25</xref>) were the only programs that did not require an MSA as input. The FastPCR algorithm for group-specific PCR (i.e. universal amplification) designs PCR primer pairs individually for each target sequence without regard for primer conservation among targets, and then compares each primer pair to the other targets. This is a brute force strategy that is only suitable for small target sets and short target sequences (appropriate for gene lengths, but not for viral genome-length sequences). We ran the FastPCR software on our internal servers, but it did not complete ‘group-specific PCR’ for the smallest data set, Norwalk virus, after running for 18 h, and for ‘multiplex PCR’ gave the error message ‘No compatible combination of pair primers for multiplex PCR found’.
</p><p>GeneUp simulates PCR with pairwise combinations of candidate primers which pass length,
<italic>T</italic><sub>m</sub>, GC%, and palindrome filters against all target sequences, and uses a greedy algorithm to build a primer set to amplify all of the targets. Since testing all possible pairwise oligonucleotide combinations against each target explodes in time and memory for large target sets, a cap on the maximum number of candidate primers to be tested must be imposed, presumably using the most common oligos, although this is not explicit in the paper. Then a PCR simulation (performing a text search for the primer in the target sequence) of each pair versus each target is performed. Unfortunately, the most common oligos may not occur in the correct orientation or distance to serve as primer pairs, so that multiple iterations of the entire process must be performed before a set covering all targets is obtained. MPP, in contrast, uses an efficient ranking algorithm to favor pairs of primer candidates that will produce amplicons of the right size in the most targets. Because of the MPP hashing algorithm and data structures, those targets that are amplified is easily determined without simulating PCR. A copy of the GeneUp software for testing could not be obtained from the authors.
</p><p>The PDA-MS/UniQ approach (
<xref rid="B26" ref-type="bibr">26</xref>) uses a hash index of 4-mers and scoring heuristic to identify common regions in the target sequences with the most shared tetramers. Simulating combinations of candidate primers selected from these common regions is then performed using a genetic algorithm. Their genetic algorithm is a more scalable approach than those above, and the quality of the solution may be improved with more compute power, although a potential disadvantage of genetic algorithms is that they can be slow, particularly if there are very many possible combinations. For computational tractability, they limited primer size to 12 nt, and avoidance of hairpins and primer dimers was not modeled. This software was not available for download or on a public web server.
</p><p>The other software all required MSA as input. While these tools could not be tested on larger (e.g. family level) primer prediction problems where alignments would not be possible or appropriate, we nevertheless attempted to run these tools on the alignable species-level target sets. The Primaclade (
<xref rid="B24" ref-type="bibr">24</xref>) webserver timed out for the smallest alignments we tested, Norwalk and FMDV. CODEHOP (
<xref rid="B23" ref-type="bibr">23</xref>) requires protein alignment as input so is not appropriate for whole-genome (nucleotide) alignments.
</p><p>GreeneSCPrimer (
<xref rid="B20" ref-type="bibr">20</xref>) did generate a number of degenerate primer candidates from the MSAs for Norwalk and FMDV, requiring the user to manually select a combination of forward and reverse groups from a set of options. We ran GreeneSCPrimer using length,
<italic>T</italic><sub>m</sub>, etc. settings mirroring or more lenient than those we used for
<xref rid="T3" ref-type="table">Table 3</xref>(
<italic>T</italic><sub>m</sub>= 55–65°C, GC% = 20–80%, length 17–25 bp, 100% coverage, product size 80–620 bp, allowed
<italic>T</italic><sub>m</sub>difference 10°C, others left as defaults). HYDEN (
<xref rid="B22" ref-type="bibr">22</xref>) also generated degenerate candidates from a MSA, although it does not check
<italic>T</italic><sub>m</sub>and the length is limited to a single value rather than a range. For our tests using HYDEN, we used a length of 18 rather than 17 because of the lack of
<italic>T</italic><sub>m</sub>control, and allowed 0 mismatches. The GreeneSCPrimer option requiring the fewest total primers for the Norwalk set required 18 primers, four of which had either 2-fold or 4-fold degeneracy so the actual number of priming sequences would be 26, compared to a total of 20 non-degenerate primers predicted by MPP (
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>). HYDEN generated four degenerate primers covering only 34 of 41 sequences, each with 3- or 4-fold degeneracy for Norwalk, which translates to 15 priming sequences in the reaction. One would need to find primers to amplify the remaining seven sequences. Small degenerate priming sets (e.g. four primers in this case) are less expensive to purchase, but because of dilution effects from the many sequence combinations actually present (15 priming sequences in the PCR), sensitivity may be reduced compared to nondegenerate priming. However, using a smaller set of degenerate signatures such as those from HYDEN may be preferable, and is a capability that could improve MPP in a future version.
</p><p>For FMDV, MPP predicted six nondegenerate multiplex compatible primers to amplify all targets. GreeneSCPrimer generated a number of candidates, and manual inspection identified that the best of those primer combinations would require 6 primers, one of which had 2-fold and another had 3-fold degeneracy, totaling 9 actual priming sequences in a reaction (
<xref ref-type="supplementary-material" rid="sup1">Supplementary Data</xref>). This compares with two primers each with 4-fold degeneracy using HYDEN (eight priming sequences in a reaction), to amplify 98% (183 of 187) targets. Again, the small number of signatures predicted by HYDEN is desirable for some applications, although the degeneracy is high and these must be supplemented to pick up the few outlying sequences
</p><p>The aim of the MuPlex software (
<xref rid="B28" ref-type="bibr">28</xref>,
<xref rid="B29" ref-type="bibr">29</xref>) is to partition primer pairs into multiplex-compatible bins for SNP genotyping, and it does not employ any algorithm during primer selection to minimize the number of primers to amplify all targets, yielding a one-to-one correspondence between number of primer pairs versus number of targets. However, the MuPlex graph-based algorithm to partition primers into bins is more sophisticated than the simple binning scheme of MPP. In future work, MPP could be used to identify a universal set of primers and the backend of MuPlex used to optimize how they are binned into separate reactions.
</p><p>In summary, no software except MPP was capable of predicting primer sets for the larger target sets we examined or generating a multiplexed primer set without a MSA for any of the target sets we examined, even the smallest ones containing only a single species. Only HYDEN and GreeneSCPrimer, both requiring MSA inputs, completed for the two smallest target sets, selecting smaller sets of degenerate primers than the nondegenerate MPP primers. If target sets are sufficiently conserved so that a reasonable MSA can be built, and degenerate primers are acceptable, these tools may be preferable over MPP. However, for larger and more diverse target sets, only MPP completed.</p></sec><sec sec-type="results"><title>Poxviridae experimental results</title><p>Using the experimental conditions described in the ‘Methods’ section, we tested the Poxviridae 16-plex (
<xref rid="T4" ref-type="table">Table 4</xref>) against vaccinia virus, Lister strain DNA. As assessed by agarose gel analysis, we achieved specific amplification of the predicted 617-bp amplicon for the target vaccinia genome (lane 3,
<xref ref-type="fig" rid="F7">Figure 7</xref>). The band does not appear in the no template control shown in lane 4 (
<xref ref-type="fig" rid="F7">Figure 7</xref>).
</p><fig id="F7" orientation="portrait" position="float"><label>Figure 7.</label><caption><p>PCR amplification of vaccinia DNA, Lister strain with the 16-plex Poxviridae primer set, with (lanes 6–10) or without (lanes 3 and 4) human DNA present in the reaction, for a reaction containing 4.8 mM MgSO
<sub>4</sub>, 0.1 µM each primer and 2.7 pg (∼10
<sup>4</sup>copies) of vaccinia Lister DNA, and using an annealing temperature of 43.9°C. Lane 3 shows specific amplification of vaccinia Lister DNA followed by the no template control in lane 4. Lanes 6–10 represent the experiments with mass ratios of vaccinia:human DNA of 1 : 1, 10 : 1, 100 : 1, 1000 : 1 and 10000 : 1, respectively. Lanes 12–14 represent amplification of human DNA only (no vaccinia present in reaction) for the same masses used at the 1 : 1, 100 : 1 and 10 000 : 1 ratios (or 2.7 pg, 0.027 pg and 2.7 × 10
<sup>–4</sup>pg human DNA, respectively). Lanes 5, 11 and 15 did not contain any sample. The arrow on the left points to the expected band (617 bp). The arrows on the right correspond to the 50-bp DNA ladder (lanes 1 and 16 contain the same ladder).
</p></caption><graphic xlink:href="gkp659f7"></graphic></fig><p>To confirm that the amplicon observed in gel analysis was a specifically amplified product from vaccinia virus, Lister strain, we sequenced the product according to the procedures described above. An analysis of the high-quality sequence read data taken from the electropherogram yielded a 95% identity (maximum) to vaccinia virus, Lister strain (AY678276.1) with a query coverage of 99% and an e-score of 0 (no data shown). These values indicate that the sequenced product was vaccinia DNA and not amplification from an exogenous nucleic acid source, e.g. host cell.</p><p>Notably, our multiplex reaction did not produce a smear (lane 3,
<xref ref-type="fig" rid="F7">Figure 7</xref>) that would be indicative of nonspecific priming of either the target or exogenous host cellular nucleic acids. While the exact manufacturer’s extraction protocol is proprietary, it is generally known that a standard sucrose gradient ultracentrifugation step is used to enrich for the viral capsids prior to viral nucleic acid extraction. However, to the best of our knowledge, no nuclease digestions are performed prior to viral capsid lysis. Furthermore, while the viral ‘extract’ contains a mixture of both viral DNA and cellular nucleic acids, the exact proportions of host cell and viral nucleic acids cannot be determined. Thus, our results indicate that a multiplex of 16 primers of 10 nt each amplifies only the specific predicted band from vaccinia.
</p><p>In some applications, such as clinical or biodetection applications, the exogenous nucleic acids from other eukaryotic sources, notably human sources, may be present in varying and unknown concentrations. To test the effects of background human genomic DNA on the performance (amplification specificity) of the Poxviridae 16-plex primer set, we conducted PCR reactions across a series of mass ratios of vaccinia Lister DNA:human genomic DNA at 1 : 1, 10 : 1, 100 : 1, 1000 : 1 and 10000 : 1. These correspond to approximate copy number ratios of vaccinia : human genomes of 1.3 × 10
<sup>4</sup>: 0.82, 1.3 × 10
<sup>4</sup>: 0.082, 1.3 × 10
<sup>4</sup>: 0.0082, 1.3 × 10
<sup>4</sup>: 0.00082 and 1.3 × 10
<sup>4</sup>: 0.000082, using genome sizes of 1.9 × 10
<sup>5</sup>bp and 3 × 10
<sup>9</sup>bp for vaccinia and human DNA, respectively. For context, there are 6.58 pg or two copies of the genome in one human cell. Our algorithm predicted 233 amplicons between 50 and 1000 bp from the human genome. In experiments, this would appear as a ‘smear’ on the agarose gel, which indeed was observed for the mass ratios of 1 : 1 and 10 : 1 (
<xref ref-type="fig" rid="F7">Figure 7</xref>, lanes 6 and 7, respectively). However, even at a ratio of 10 : 1 vaccinia:human DNA (lane 7), the 617-bp vaccinia amplicon is clearly visible on the gel, despite the numerous nonspecific amplicons. At ratios of 100 : 1, 1000 : 1 and 10000 : 1 (lanes 8, 9 and 10, respectively,
<xref ref-type="fig" rid="F7">Figure 7</xref>), the smear is drastically reduced to nonexistent (at the resolution of the agarose gel), and the vaccinia amplicon is readily visible. For comparison, we tested the 16-plex primers against the same mass of human DNA in the absence of vaccinia DNA at 2.7 pg (the 1 : 1 ratio mass), 0.027 pg (100 : 1 mass) and 2.7 × 10
<sup>−4</sup>pg (10000 : 1 mass) (lanes 12, 13, and 14, respectively,
<xref ref-type="fig" rid="F7">Figure 7</xref>). The data show a similar smear at 2.7 pg as was observed when vaccinia DNA was present (lane 6). However, the low-intensity 617-bp amplicon from vaccinia is visible in lane 6 (with vaccinia) while a similar amplicon is not present in lane 12 (without vaccina). These results provide evidence that amplification with these family level primer sets will depend on viral titers in the actual sample, and that there are cases where amplification will either not be possible or require additional sample purification steps. However, as discussed further in the next section, specific amplification with a short-primer multiplex could be used for selective enrichment of viral targets by generating amplicons with known sequence, which may be feasible for viral detection if combined with a probe-based amplicon detection method such as TaqMan® or Luminex bead based suspension arrays (
<ext-link ext-link-type="uri" xlink:href="http://www.luminexcorp.com/">http://www.luminexcorp.com/</ext-link>).
</p><p>We compared recently published conserved Orthopoxvirus primers (
<xref rid="B36" ref-type="bibr">36</xref>) to the 148 Poxviridae genomes, and computational predictions suggest that 67 of the available genomes might not be amplified by the two primers they designed for the Orthopox genus. These included a number of monkeypox, ectromelia, several vaccinia and a couple of variola strains. However, in permissive hybridization conditions it is possible that primers would anneal despite mismatches to target, allowing more of the targets to be amplified. Four conserved Orthopoxvirus primers from an earlier publication (
<xref rid="B37" ref-type="bibr">37</xref>), before many of the Poxviridae genomes were available, do not match 53 genomes, including a number of monkeypox, ectromelia, camelpox and one variola minor genome. The primers in (
<xref rid="B37" ref-type="bibr">37</xref>) included inosine bases, which we replaced with each possible A, T, G and C base in all possible combinations for our primer-target comparison. Thus, published Orthopoxvirus primers may fail to amplify many desired targets.
</p></sec></sec><sec sec-type="discussion"><title>DISCUSSION AND CONCLUSIONS</title><p>Primer design for amplification and detection of divergent target sequences can be challenging, and this problem will only grow as sequencing technologies improve. Some methods are limited in scalability, particularly those requiring a MSA as input. Developing a PCR multiplex is often a tedious mix-and-match process from among primers originally designed to work in singleplex. We describe the MPP algorithm based on hashing of conserved
<italic>k</italic>-mer subsequences that requires no MSA and where multiplex-compatible primer sets are built
<italic>de novo</italic>to avoid primer dimer and hairpin formation to the extent that can be predicted based on free energy calculations.
</p><p>The algorithm seeks to minimize the number of primers to amplify all targets, although because the algorithm is heuristic and not an exhaustive search of all combinations of primers, the smallest primer sets may not always be selected. This is an NP complete problem (
<xref rid="B38" ref-type="bibr">38</xref>), so an exhaustive search for the global optimum is only practical for very small target sets. It is also beyond the scope of the current software to predict priming with mismatches between target and primer, which may occur under non-stringent hybridization conditions or for mismatches at the 5′ end of a primer. While allowing for such priming would reduce the number of primers needed to amplify all targets, it would substantially slow the algorithm and increase memory requirements. Nor do we claim to have taken all steps to optimize (e.g. through parallelization, see
<xref ref-type="supplementary-material" rid="sup1">Supplementary Methods</xref>) for speed or memory, but instead present this software as a simple embodiment of one alternative to MSA for multiplex primer design.
</p><p>The software handles a large number of input sequences, although for very diverse targets the predicted primer multiplex may be too large to be empirically feasible. For many target sets, such as those including all the genomes of a species, predicting a universal primer set requires only minutes up to a few hours, although for inputs with thousands of sequences a run may take days. We used MPP to design a universal primer set for a target set of all virus genomes, and showed that even if short primers are used, thousands of priming sequences would be required to amplify all sequenced viruses. We then applied MPP to design multiplex family level primers for every viral family, as well as for some diverse species target sets too large for other available primer prediction software.</p><p>In addition to finding primers for viral detection, another application of MPP could be to design multiplex primers for homologs in a gene family. The user can specify an appropriate product length range to ensure amplification of an adequate span across the gene. MPP could also be used to design multiplex primer sets for unrelated target sequences, for example, multiple bacterial and viral species or gene families, in a single reaction. Since no sequence alignment is required, there is no need for any sequence conservation among targets. Recent work showed that a large-scale multiplex of 800 primer pairs specifically designed to detect a diverse set of genes from nine pathogens improved sensitivity on a microarray by up to 1000-fold (
<xref rid="B39" ref-type="bibr">39</xref>). The MPP software could be used to design multiplex primer sets for similar work.
</p><p>We demonstrated the application of the algorithm experimentally using purified vaccinia DNA, amplifying the expected band with a Poxviridae short primer multiplex PCR, and confirmed the band’s expected sequence. We provided an example of the effect of background nucleic acids by spiking human DNA into the PCR reaction over a series of mass ratios. While a band for vaccinia was visualized by slab gel electrophoresis for all mass ratios, background DNA smears were also clear at the higher ratios. At a mass ratio of 1 : 1 vaccinia:human DNA, there are 2.7 pg or 0.82 copies of the human genome present in the reaction. To place this into context, there are ∼6.58 pg or two copies of genomic DNA in one human cell. Thus, the results from the 1 : 1 mass ratio reaction represents the impact that ∼41% of the total DNA content from one human cell could have on the results, when analyzed by slab gel electrophoresis. In a clinical sample, there could be a much greater number of human cells and DNA. Amplification with these family level primer sets will depend on viral titers in the actual sample, and these experiments show that there are cases where amplification will either not be possible or require additional sample purification steps. While short primer PCR multiplexes may enable amplification of diverse target sets, they will not match the specificity of longer primers. Combination of these short primer multiplexes with digitized, microfluidic-based picoliter reactions (
<xref rid="B40" ref-type="bibr">40</xref>) and/or next-generation microfluidic-based sample purification technologies that have the ability to isolate target from contaminating nucleic acids may help to overcome this limitation. However, even in the absence of such technologies, a highly conserved short primer multiplex could enrich amplification products for members of a particular viral family compared to randomly amplified or unamplified sample. Then it would need to be followed by product sequencing or array hybridization, such as TaqMan® or Luminex bead-based suspension arrays (
<ext-link ext-link-type="uri" xlink:href="http://www.luminexcorp.com/">http://www.luminexcorp.com/</ext-link>), to provide more specific information about the organism(s) present, since electrophoretic banding patterns may contain unexpected bands, particularly if the sample contains a significant amount of nonviral nucleic acids.
</p><p>Previous studies have shown the utility of short primer singleplex PCR using 9-mers or 10-mers followed by gel electrophoresis for genetic fingerprinting of eukaryotes and bacteria (
<xref rid="B41" ref-type="bibr">41</xref>,
<xref rid="B42" ref-type="bibr">42</xref>). For viruses, with much smaller and more diverse genomes, the large numbers of 9-mer or 10-mer primers required to generate at least one band from every virus as predicted by our analyses implies that primer size would need to be as short as 5-mers to rely on a gel banding pattern using only one priming sequence for fingerprinting viruses (unpublished data). However, the analyses here predict that imperfect sample purification to eliminate eukaryotic nucleic acids could be problematic for universal viral priming using primers shorter than 15 bases, particularly for multiplexes of 10 or more primers. Nanda
<italic>et al.</italic>(
<xref rid="B43" ref-type="bibr">43</xref>,
<xref rid="B44" ref-type="bibr">44</xref>) were able to achieve sufficient viral isolation from cell culture samples to allow viral identification using viral PCR with priming sequences as short as pentamers, so the problem of contaminating host nucleic acids for specific, short primer PCR of viruses is not insurmountable. They found specific pentamer PCR to be several logs more sensitive than nonspecific amplification, provided that they purified encapsidated viral nucleic acids prior to PCR. Another method that has been used for virus discovery is VIDISCA (Virus discovery cDNA-AFLP) using restriction enzyme digestion, adaptor ligation, and PCR by priming with the adaptor sequence (
<xref rid="B45" ref-type="bibr">45</xref>,
<xref rid="B46" ref-type="bibr">46</xref>). This method, like the pentamer priming used by (
<xref rid="B43" ref-type="bibr">43</xref>,
<xref rid="B44" ref-type="bibr">44</xref>), requires prior separation of encapsidated viral nucleic acids, as it generates fragments from any DNA present, viral, host or otherwise. Multiplex PCR with primers 10–15 nt in length may be yet another alternative strategy lying between these nonspecific methods and PCR with standard primers of at least 18 nt, as we have shown that it can add some measure of specificity for a viral family.
</p><p>In summary, we applied the MPP software to generate multiplex-compatible primer sets for every viral family and several divergent viral species and experimentally demonstrated application of one multiplex set to show nonrandom amplification with a set of short primers.</p></sec><sec><title>Code and Primer Availability</title><p>Primer sets are available as
<xref ref-type="supplementary-material" rid="sup1">supplemental material</xref>or by request from the authors. The MPP software is freely available for academic and nonprofit use at
<ext-link ext-link-type="uri" xlink:href="http://mpp.llnl.gov">http://mpp.llnl.gov</ext-link>.
</p></sec><sec><title>FUNDING</title><p>Laboratory Directed Research and Development (LDRD) and Computations TechBase awards from the Lawrence Livermore National Laboratory. Funding for open access charge: Lawrence Livermore National Laboratory.</p><p><italic>Conflict of interest statement</italic>. None declared.
</p></sec><sec sec-type="supplementary-material"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="sup1"><label>Supplementary Data</label><media xlink:href="gkp659_Supplementary_Data.zip"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><title>ACKNOWLEDGEMENTS</title><p>The authors wish to thank anonymous reviewers for helpful comments. Auspices statement (required by LLNL): This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.</p></ack><ref-list><title>REFERENCES</title><ref id="B1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Victoria</surname><given-names>JG</given-names></name><name><surname>Kapoor</surname><given-names>A</given-names></name><name><surname>Li</surname><given-names>L</given-names></name><name><surname>Blinkova</surname><given-names>O</given-names></name><name><surname>Slikas</surname><given-names>B</given-names></name><name><surname>Wang</surname><given-names>C</given-names></name><name><surname>Naeem</surname><given-names>A</given-names></name><name><surname>Zaidi</surname><given-names>S</given-names></name><name><surname>Delwart</surname><given-names>E</given-names></name></person-group><article-title>Metagenomic analyses of viruses in the stool of children with acute flaccid paralysis</article-title><source>J. Virol.</source><year>2009</year><comment>doi:10.1128/JVI.02301-08</comment></element-citation></ref><ref id="B2"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edwards</surname><given-names>RA</given-names></name><name><surname>Rohwer</surname><given-names>F</given-names></name></person-group><article-title>Viral metagenomics</article-title><source>Nat. Rev. Micro.</source><year>2005</year><volume>3</volume><fpage>504</fpage><lpage>510</lpage></element-citation></ref><ref id="B3"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nakamura</surname><given-names>S</given-names></name><name><surname>Yang</surname><given-names>C-S</given-names></name><name><surname>Sakon</surname><given-names>N</given-names></name><name><surname>Ueda</surname><given-names>M</given-names></name><name><surname>Tougan</surname><given-names>T</given-names></name><name><surname>Yamashita</surname><given-names>A</given-names></name><name><surname>Goto</surname><given-names>N</given-names></name><name><surname>Takahashi</surname><given-names>K</given-names></name><name><surname>Yasunaga</surname><given-names>T</given-names></name><name><surname>Ikuta</surname><given-names>K</given-names></name><etal></etal></person-group><article-title>Direct metagenomic detection of viral pathogens in nasal and fecal specimens using an unbiased high-throughput sequencing approach</article-title><source>PLoS ONE</source><year>2009</year><volume>4</volume><fpage>e4219</fpage><pub-id pub-id-type="pmid">19156205</pub-id></element-citation></ref><ref id="B4"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kistler</surname><given-names>A</given-names></name><name><surname>Avila</surname><given-names>P</given-names></name><name><surname>Rouskin</surname><given-names>S</given-names></name><name><surname>Wang</surname><given-names>D</given-names></name><name><surname>Ward</surname><given-names>T</given-names></name><name><surname>Yagi</surname><given-names>S</given-names></name><name><surname>Schnurr</surname><given-names>D</given-names></name><name><surname>Ganem</surname><given-names>D</given-names></name><name><surname>DeRisi</surname><given-names>J</given-names></name><name><surname>Boushey</surname><given-names>H</given-names></name></person-group><article-title>Pan-viral screening of respiratory tract infections in adults with and without asthma reveals unexpected human coronavirus and human rhinovirus diversity</article-title><source>J. Infect. Dis.</source><year>2007</year><volume>196</volume><fpage>817</fpage><lpage>825</lpage><pub-id pub-id-type="pmid">17703411</pub-id></element-citation></ref><ref id="B5"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname><given-names>B</given-names></name><name><surname>Blaney</surname><given-names>KM</given-names></name><name><surname>Malanoski</surname><given-names>AP</given-names></name><name><surname>Ligler</surname><given-names>AG</given-names></name><name><surname>Schnur</surname><given-names>JM</given-names></name><name><surname>Metzgar</surname><given-names>D</given-names></name><name><surname>Russell</surname><given-names>KL</given-names></name><name><surname>Stenger</surname><given-names>DA</given-names></name></person-group><article-title>Using a resequencing microarray as a multiple respiratory pathogen detection assay</article-title><source>J. Clin. Microbiol.</source><year>2007</year><volume>45</volume><fpage>443</fpage><lpage>452</lpage><pub-id pub-id-type="pmid">17135438</pub-id></element-citation></ref><ref id="B6"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname><given-names>B</given-names></name><name><surname>Wang</surname><given-names>Z</given-names></name><name><surname>Vora</surname><given-names>GJ</given-names></name><name><surname>Thornton</surname><given-names>JA</given-names></name><name><surname>Schnur</surname><given-names>JM</given-names></name><name><surname>Thach</surname><given-names>DC</given-names></name><name><surname>Blaney</surname><given-names>KM</given-names></name><name><surname>Ligler</surname><given-names>AG</given-names></name><name><surname>Malanoski</surname><given-names>AP</given-names></name><name><surname>Santiago</surname><given-names>J</given-names></name><etal></etal></person-group><article-title>Broad-spectrum respiratory tract pathogen identification using resequencing DNA microarrays</article-title><source>Genome Res.</source><year>2006</year><comment>gr.4337206</comment></element-citation></ref><ref id="B7"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>D</given-names></name><name><surname>Coscoy</surname><given-names>L</given-names></name><name><surname>Zylberberg</surname><given-names>M</given-names></name><name><surname>Avila</surname><given-names>P</given-names></name><name><surname>Boushey</surname><given-names>H</given-names></name><name><surname>Ganem</surname><given-names>D</given-names></name><name><surname>DeRisi</surname><given-names>J</given-names></name></person-group><article-title>Microarray-based detection and genotyping of viral pathogens</article-title><source>PNAS</source><year>2002</year><volume>99</volume></element-citation></ref><ref id="B8"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Quan</surname><given-names>PL</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Jabado</surname><given-names>OJ</given-names></name><name><surname>Conlan</surname><given-names>S</given-names></name><name><surname>Hirschberg</surname><given-names>DL</given-names></name><name><surname>Pozo</surname><given-names>F</given-names></name><name><surname>Jack</surname><given-names>PJM</given-names></name><name><surname>Cisterna</surname><given-names>D</given-names></name><name><surname>Renwick</surname><given-names>N</given-names></name><name><surname>Hui</surname><given-names>J</given-names></name><etal></etal></person-group><article-title>Detection of respiratory viruses and subtype identification of influenza A viruses by GreeneChipResp oligonucleotide microarray</article-title><source>J. Clin. Microbiol.</source><year>2007</year><volume>45</volume><fpage>2359</fpage><lpage>2364</lpage><pub-id pub-id-type="pmid">17553978</pub-id></element-citation></ref><ref id="B9"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Quan</surname><given-names>P</given-names></name><name><surname>Jabado</surname><given-names>O</given-names></name><name><surname>Conlan</surname><given-names>S</given-names></name><name><surname>Hirschberg</surname><given-names>D</given-names></name><name><surname>Liu</surname><given-names>Y</given-names></name><name><surname>Zhai</surname><given-names>J</given-names></name><name><surname>Renwick</surname><given-names>N</given-names></name><name><surname>Hui</surname><given-names>J</given-names></name><name><surname>Hegyi</surname><given-names>H</given-names></name><etal></etal></person-group><article-title>Panmicrobial oligonucleotide array for diagnosis of infectious diseases</article-title><source>Emerg. Infect. Dis.</source><year>2007</year><volume>13</volume><fpage>73</fpage><lpage>81</lpage><pub-id pub-id-type="pmid">17370518</pub-id></element-citation></ref><ref id="B10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ecker</surname><given-names>D</given-names></name><name><surname>Drader</surname><given-names>J</given-names></name><name><surname>Gutierrez</surname><given-names>J</given-names></name><name><surname>Gutierrez</surname><given-names>A</given-names></name><name><surname>Hannis</surname><given-names>J</given-names></name><name><surname>Schink</surname><given-names>A</given-names></name><name><surname>Sampath</surname><given-names>R</given-names></name><name><surname>Blyn</surname><given-names>LB</given-names></name><name><surname>Eschoo</surname><given-names>MW</given-names></name><name><surname>Hall</surname><given-names>TA</given-names></name><etal></etal></person-group><article-title>The Ibis T5000 universal biosensor: an automated platform for pathogen identification and strain typing</article-title><source>JALA</source><year>2006</year><volume>11</volume><fpage>341</fpage><lpage>351</lpage></element-citation></ref><ref id="B11"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sampath</surname><given-names>R</given-names></name><name><surname>Hall</surname><given-names>TA</given-names></name><name><surname>Massire</surname><given-names>C</given-names></name><name><surname>Li</surname><given-names>F</given-names></name><name><surname>Blyn</surname><given-names>LB</given-names></name><name><surname>Eshoo</surname><given-names>MW</given-names></name><name><surname>Hofstadler</surname><given-names>SA</given-names></name><name><surname>Ecker</surname><given-names>DJ</given-names></name></person-group><article-title>Rapid identification of emerging infectious agents using PCR and electrospray ionization mass spectrometry</article-title><source>Ann. NY Acad. Sci.</source><year>2007</year><volume>1102</volume><fpage>109</fpage><lpage>120</lpage><pub-id pub-id-type="pmid">17470915</pub-id></element-citation></ref><ref id="B12"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lamson</surname><given-names>D</given-names></name><name><surname>Renwick</surname><given-names>N</given-names></name><name><surname>Kapoor</surname><given-names>V</given-names></name><name><surname>Liu</surname><given-names>Z</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Ju</surname><given-names>J</given-names></name><name><surname>Dean</surname><given-names>A</given-names></name><name><surname>St.George</surname><given-names>K</given-names></name><name><surname>Briese</surname><given-names>T</given-names></name><name><surname>Lipkin</surname><given-names>WI</given-names></name></person-group><article-title>Mass tag polymerase chain reaction detection of respiratory pathogens, including a new rhinovirus genotype, that caused influenza-like illness in New York State during 2004–2005</article-title><source>J. Infect. Dis.</source><year>2006</year><volume>194</volume><fpage>1398</fpage><lpage>1402</lpage><pub-id pub-id-type="pmid">17054069</pub-id></element-citation></ref><ref id="B13"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Briese</surname><given-names>T</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Kokoris</surname><given-names>M</given-names></name><name><surname>Jabado</surname><given-names>O</given-names></name><name><surname>Zhiqiang</surname><given-names>L</given-names></name><name><surname>Renwick</surname><given-names>N</given-names></name><name><surname>Kapoor</surname><given-names>V</given-names></name><name><surname>Casas</surname><given-names>I</given-names></name><name><surname>Pozo</surname><given-names>F</given-names></name><name><surname>Limberger</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>Diagnostic system for rapid and sensitive differential detection of pathogens</article-title><source>Emerg. Infect. Dis.</source><year>2005</year><volume>11</volume><fpage>310</fpage><lpage>313</lpage><pub-id pub-id-type="pmid">15752453</pub-id></element-citation></ref><ref id="B14"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dominguez</surname><given-names>SR</given-names></name><name><surname>Briese</surname><given-names>T</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Hui</surname><given-names>J</given-names></name><name><surname>Villari</surname><given-names>J</given-names></name><name><surname>Kapoor</surname><given-names>V</given-names></name><name><surname>Tokarz</surname><given-names>R</given-names></name><name><surname>Glodé</surname><given-names>MP</given-names></name><name><surname>Anderson</surname><given-names>MS</given-names></name><name><surname>Robinson</surname><given-names>CC</given-names></name><etal></etal></person-group><article-title>Multiplex MassTag-PCR for respiratory pathogens in pediatric nasopharyngeal washes negative by conventional diagnostic testing shows a high prevalence of viruses belonging to a newly recognized rhinovirus clade</article-title><source>J. Clin. Virol.</source><year>2008</year><volume>43</volume><fpage>219</fpage><lpage>222</lpage><pub-id pub-id-type="pmid">18674964</pub-id></element-citation></ref><ref id="B15"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fox</surname><given-names>JD</given-names></name></person-group><article-title>Nucleic acid amplification tests for detection of respiratory viruses</article-title><source>J. Clin. Virol.</source><year>2007</year><volume>40</volume><fpage>S15</fpage><lpage>S23</lpage><pub-id pub-id-type="pmid">18162249</pub-id></element-citation></ref><ref id="B16"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lenhoff</surname><given-names>RJ</given-names></name><name><surname>Naraghi-Arani</surname><given-names>P</given-names></name><name><surname>Thissen</surname><given-names>JB</given-names></name><name><surname>Olivas</surname><given-names>J</given-names></name><name><surname>Carillo</surname><given-names>AC</given-names></name><name><surname>Chinn</surname><given-names>C</given-names></name><name><surname>Rasmussen</surname><given-names>M</given-names></name><name><surname>Messenger</surname><given-names>SM</given-names></name><name><surname>Suer</surname><given-names>LD</given-names></name><name><surname>Smith</surname><given-names>SM</given-names></name><etal></etal></person-group><article-title>Multiplexed molecular assay for rapid exclusion of foot-and-mouth disease</article-title><source>J. Virol. Methods</source><year>2008</year><volume>153</volume><fpage>61</fpage><lpage>69</lpage><pub-id pub-id-type="pmid">18634827</pub-id></element-citation></ref><ref id="B17"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Griffiths</surname><given-names>D</given-names></name><name><surname>Kellam</surname><given-names>P</given-names></name><name><surname>Weiss</surname><given-names>R</given-names></name></person-group><article-title>International patent no</article-title><year>2002</year><comment>WO/2002/099130</comment></element-citation></ref><ref id="B18"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Quan</surname><given-names>PL</given-names></name><name><surname>Briese</surname><given-names>T</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Lipkin</surname><given-names>WI</given-names></name></person-group><article-title>Rapid sequence-based diagnosis of viral infection</article-title><source>Antivir. Res.</source><year>2008</year><volume>79</volume><fpage>1</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">18367256</pub-id></element-citation></ref><ref id="B19"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lemmon</surname><given-names>GH</given-names></name><name><surname>Gardner</surname><given-names>SN</given-names></name></person-group><article-title>Predicting the sensitivity and specificity of published real-time PCR assays</article-title><source>Ann. Clin. Microbiol. Antimicrob.</source><year>2008</year><volume>7</volume><fpage>18</fpage><comment>doi:10.1186/1476-0711-1187-1118</comment><pub-id pub-id-type="pmid">18817537</pub-id></element-citation></ref><ref id="B20"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jabado</surname><given-names>OJ</given-names></name><name><surname>Palacios</surname><given-names>G</given-names></name><name><surname>Kapoor</surname><given-names>V</given-names></name><name><surname>Hui</surname><given-names>J</given-names></name><name><surname>Renwick</surname><given-names>N</given-names></name><name><surname>Zhai</surname><given-names>J</given-names></name><name><surname>Briese</surname><given-names>T</given-names></name><name><surname>Lipkin</surname><given-names>WI</given-names></name></person-group><article-title>Greene SCPrimer: a rapid comprehensive tool for designing degenerate primers from multiple sequence alignments</article-title><source>Nucleic Acids Res.</source><year>2006</year><volume>34</volume><fpage>6605</fpage><lpage>6611</lpage><pub-id pub-id-type="pmid">17135211</pub-id></element-citation></ref><ref id="B21"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jarman</surname><given-names>SN</given-names></name></person-group><article-title>Amplicon: software for designing PCR primers on aligned DNA sequences</article-title><source>Bioinformatics</source><year>2004</year><volume>20</volume><fpage>1644</fpage><lpage>1645</lpage><pub-id pub-id-type="pmid">14962918</pub-id></element-citation></ref><ref id="B22"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Linhart</surname><given-names>C</given-names></name><name><surname>Shamir</surname><given-names>R</given-names></name></person-group><article-title>The degenerate primer design problem</article-title><source>Bioinformatics</source><year>2002</year><volume>18</volume><fpage>S172</fpage><lpage>S181</lpage></element-citation></ref><ref id="B23"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rose</surname><given-names>T</given-names></name><name><surname>Henikoff</surname><given-names>J</given-names></name><name><surname>Henikoff</surname><given-names>S</given-names></name></person-group><article-title>CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primer) PCR primer design</article-title><source>Nuclei Acids Res.</source><year>2003</year><volume>31</volume><fpage>3763</fpage><lpage>3766</lpage></element-citation></ref><ref id="B24"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gadberry</surname><given-names>MD</given-names></name><name><surname>Malcomber</surname><given-names>ST</given-names></name><name><surname>Doust</surname><given-names>AN</given-names></name><name><surname>Kellogg</surname><given-names>EA</given-names></name></person-group><article-title>Primaclade–a flexible tool to find conserved PCR primers across multiple species</article-title><source>Bioinformatics</source><year>2005</year><volume>21</volume><fpage>1263</fpage><lpage>1264</lpage><pub-id pub-id-type="pmid">15539448</pub-id></element-citation></ref><ref id="B25"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pesole</surname><given-names>G</given-names></name><name><surname>Liuni</surname><given-names>S</given-names></name><name><surname>Grillo</surname><given-names>G</given-names></name><name><surname>Belichard</surname><given-names>P</given-names></name><name><surname>Trenkle</surname><given-names>T</given-names></name><name><surname>Welsh</surname><given-names>J</given-names></name><name><surname>McClelland</surname><given-names>M</given-names></name></person-group><article-title>GeneUp: a program to select short PCR primer pairs that occur in multiple members of sequence lists</article-title><source>BioTechniques</source><year>1998</year><volume>25</volume><fpage>112</fpage><lpage>123</lpage><pub-id pub-id-type="pmid">9668985</pub-id></element-citation></ref><ref id="B26"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname><given-names>YC</given-names></name><name><surname>Chang</surname><given-names>CF</given-names></name><name><surname>Chan</surname><given-names>Ch</given-names></name><name><surname>Yeh</surname><given-names>TJ</given-names></name><name><surname>Chang</surname><given-names>YC</given-names></name><name><surname>Chen</surname><given-names>CC</given-names></name><name><surname>Kao</surname><given-names>CY</given-names></name></person-group><article-title>Integrated minimum-set primers and unique probe design algorithms for differential detection on symptom-related pathogens</article-title><source>Bioinformatics</source><year>2005</year><volume>21</volume><fpage>4330</fpage><lpage>4337</lpage><pub-id pub-id-type="pmid">16249263</pub-id></element-citation></ref><ref id="B27"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Markham</surname><given-names>N</given-names></name><name><surname>Zuker</surname><given-names>M</given-names></name></person-group><article-title>DINAMelt web server for nucleic acid melting prediction</article-title><source>Nucleic Acids Res.</source><year>2005</year><volume>33</volume><fpage>W577</fpage><lpage>W581</lpage><pub-id pub-id-type="pmid">15980540</pub-id></element-citation></ref><ref id="B28"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rachlin</surname><given-names>J</given-names></name><name><surname>Ding</surname><given-names>C</given-names></name><name><surname>Cantor</surname><given-names>C</given-names></name><name><surname>Kasif</surname><given-names>S</given-names></name></person-group><article-title>MuPlex: multi-objective multiplex PCR assay design</article-title><source>Nucleic Acids Res.</source><year>2005</year><volume>33</volume><fpage>W544</fpage><lpage>W547</lpage><pub-id pub-id-type="pmid">15980531</pub-id></element-citation></ref><ref id="B29"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rachlin</surname><given-names>J</given-names></name><name><surname>Ding</surname><given-names>C</given-names></name><name><surname>Cantor</surname><given-names>C</given-names></name><name><surname>Kasif</surname><given-names>S</given-names></name></person-group><article-title>Computational tradeoffs in multiplex PCR assay design for SNP genotyping</article-title><source>BMC Genomics</source><year>2005</year><volume>6</volume><fpage>102</fpage><pub-id pub-id-type="pmid">16042802</pub-id></element-citation></ref><ref id="B30"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname><given-names>RC</given-names></name></person-group><article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title><source>Nucleic Acids Res.</source><year>2004</year><volume>32</volume><fpage>1792</fpage><lpage>1797</lpage><pub-id pub-id-type="pmid">15034147</pub-id></element-citation></ref><ref id="B31"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chenna</surname><given-names>R</given-names></name><name><surname>Sugawara</surname><given-names>H</given-names></name><name><surname>Koike</surname><given-names>T</given-names></name><name><surname>Lopez</surname><given-names>R</given-names></name><name><surname>Gibson</surname><given-names>TJ</given-names></name><name><surname>Higgins</surname><given-names>DG</given-names></name><name><surname>Thompson</surname><given-names>JD</given-names></name></person-group><article-title>Multiple sequence alignment with the Clustal series of programs</article-title><source>Nucleic Acids Res.</source><year>2003</year><volume>31</volume><fpage>3497</fpage><lpage>3500</lpage><pub-id pub-id-type="pmid">12824352</pub-id></element-citation></ref><ref id="B32"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eddy</surname><given-names>SR</given-names></name></person-group><article-title>Profile hidden Markov models</article-title><source>Bioinformatics</source><year>1998</year><volume>14</volume><fpage>755</fpage><lpage>763</lpage><pub-id pub-id-type="pmid">9918945</pub-id></element-citation></ref><ref id="B33"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Herold</surname><given-names>J</given-names></name><name><surname>Kurtz</surname><given-names>S</given-names></name><name><surname>Giegerich</surname><given-names>R</given-names></name></person-group><article-title>Efficient computation of absent words in genomic sequences</article-title><source>BMC Bioinformatics</source><year>2008</year><volume>9</volume><fpage>167</fpage><pub-id pub-id-type="pmid">18366790</pub-id></element-citation></ref><ref id="B34"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname><given-names>SF</given-names></name><name><surname>Madden</surname><given-names>T</given-names></name><name><surname>Schaffer</surname><given-names>A</given-names></name><name><surname>Zhang</surname><given-names>J</given-names></name><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Miller</surname><given-names>W</given-names></name><name><surname>Lipman</surname><given-names>D</given-names></name></person-group><article-title>Gapped BLAST and PSI-BLAST: a new generation of protein database search programs</article-title><source>Nucleic Acids Res.</source><year>1997</year><volume>25</volume><fpage>3389</fpage><lpage>3402</lpage><pub-id pub-id-type="pmid">9254694</pub-id></element-citation></ref><ref id="B35"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Casabianca</surname><given-names>A</given-names></name><name><surname>Gori</surname><given-names>C</given-names></name><name><surname>Orlandi</surname><given-names>C</given-names></name><name><surname>Forbici</surname><given-names>F</given-names></name><name><surname>Federico Perno</surname><given-names>C</given-names></name><name><surname>Magnani</surname><given-names>M</given-names></name></person-group><article-title>Fast and sensitive quantitative detection of HIV DNA in whole blood leucocytes by SYBR green I real-time PCR assay</article-title><source>Mol. Cell Probes</source><year>2007</year><volume>21</volume><fpage>368</fpage><lpage>378</lpage><pub-id pub-id-type="pmid">17629450</pub-id></element-citation></ref><ref id="B36"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Putkuri</surname><given-names>N</given-names></name><name><surname>Piiparinen</surname><given-names>H</given-names></name><name><surname>Vaheri</surname><given-names>A</given-names></name><name><surname>Vapalahti</surname><given-names>O</given-names></name></person-group><article-title>Detection of human Orthopoxvirus infections and differentiation of smallpox virus with real-time PCR</article-title><source>J. Med. Virol.</source><year>2009</year><volume>81</volume><fpage>146</fpage><lpage>152</lpage><pub-id pub-id-type="pmid">19031452</pub-id></element-citation></ref><ref id="B37"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Olson</surname><given-names>VA</given-names></name><name><surname>Laue</surname><given-names>T</given-names></name><name><surname>Laker</surname><given-names>MT</given-names></name><name><surname>Babkin</surname><given-names>IV</given-names></name><name><surname>Drosten</surname><given-names>C</given-names></name><name><surname>Shchelkunov</surname><given-names>SN</given-names></name><name><surname>Niedrig</surname><given-names>M</given-names></name><name><surname>Damon</surname><given-names>IK</given-names></name><name><surname>Meyer</surname><given-names>H</given-names></name></person-group><article-title>Real-time PCR system for detection of Orthopoxviruses and simultaneous identification of smallpox virus</article-title><source>J. Clin. Microbiol.</source><year>2004</year><volume>42</volume><fpage>1940</fpage><lpage>1946</lpage><pub-id pub-id-type="pmid">15131152</pub-id></element-citation></ref><ref id="B38"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pearson</surname><given-names>WR</given-names></name><name><surname>Robbins</surname><given-names>G</given-names></name><name><surname>Wrege</surname><given-names>DE</given-names></name><name><surname>Zhang</surname><given-names>T</given-names></name></person-group><article-title>A new approach to primer selection in polymerase chain reaction experiments</article-title><source>Proc. Int. Conf. Intell. Syst. Mol. Biol.</source><year>1995</year><volume>3</volume><fpage>285</fpage><lpage>291</lpage><pub-id pub-id-type="pmid">7584448</pub-id></element-citation></ref><ref id="B39"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Palka-Santini</surname><given-names>M</given-names></name><name><surname>Cleven</surname><given-names>B</given-names></name><name><surname>Eichinger</surname><given-names>L</given-names></name><name><surname>Kronke</surname><given-names>M</given-names></name><name><surname>Krut</surname><given-names>O</given-names></name></person-group><article-title>Large scale multiplex PCR improves pathogen detection by DNA microarrays</article-title><source>BMC Microbiol.</source><year>2009</year><volume>9</volume><fpage>1</fpage><pub-id pub-id-type="pmid">19121223</pub-id></element-citation></ref><ref id="B40"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beer</surname><given-names>NR</given-names></name><name><surname>Hindson</surname><given-names>BJ</given-names></name><name><surname>Wheeler</surname><given-names>EK</given-names></name><name><surname>Hall</surname><given-names>SB</given-names></name><name><surname>Rose</surname><given-names>KA</given-names></name><name><surname>Kennedy</surname><given-names>IM</given-names></name><name><surname>Colston</surname><given-names>BW</given-names></name></person-group><article-title>On-chip, real-time, single-copy polymerase chain reaction in picoliter droplets</article-title><source>Anal. Chem.</source><year>2007</year><volume>79</volume><fpage>8471</fpage><lpage>8475</lpage><pub-id pub-id-type="pmid">17929880</pub-id></element-citation></ref><ref id="B41"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caetano-Anolles</surname><given-names>G</given-names></name><name><surname>Brant</surname><given-names>JB</given-names></name><name><surname>Peter</surname><given-names>MG</given-names></name></person-group><article-title>DNA amplification fingerprinting using very short arbitrary oligonucleotide primers</article-title><source>Nat. Biotech.</source><year>1991</year><volume>9</volume><fpage>553</fpage><lpage>557</lpage></element-citation></ref><ref id="B42"><label>42</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Williams</surname><given-names>JGK</given-names></name><name><surname>Kubelik</surname><given-names>AR</given-names></name><name><surname>Livak</surname><given-names>KJ</given-names></name><name><surname>Rafalski</surname><given-names>JA</given-names></name><name><surname>Tingey</surname><given-names>SV</given-names></name></person-group><article-title>DNA polymorphisms amplified by arbitrary primers are useful as genetic markers</article-title><source>Nucleic Acids Res.</source><year>1990</year><volume>18</volume><fpage>6531</fpage><lpage>6535</lpage><pub-id pub-id-type="pmid">1979162</pub-id></element-citation></ref><ref id="B43"><label>43</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Nanda</surname><given-names>S</given-names></name></person-group><source>ASM Biodefense and Emerging Diseases Research Meeting</source><year>2007</year><publisher-loc>DC</publisher-loc><publisher-name>Washington</publisher-name><fpage>53</fpage></element-citation></ref><ref id="B44"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nanda</surname><given-names>S</given-names></name><name><surname>Jayan</surname><given-names>G</given-names></name><name><surname>Voulgaropoulou</surname><given-names>F</given-names></name><name><surname>Sierra-Honigmann</surname><given-names>AM</given-names></name><name><surname>Uhlenhaut</surname><given-names>C</given-names></name><name><surname>McWatters</surname><given-names>BJP</given-names></name><name><surname>Patel</surname><given-names>A</given-names></name><name><surname>Krause</surname><given-names>PR</given-names></name></person-group><article-title>Universal virus detection by degenerate-oligonucleotide primed polymerase chain reaction of purified viral nucleic acids</article-title><source>J. Virol. Methods</source><year>2008</year><volume>152</volume><fpage>18</fpage><lpage>24</lpage><pub-id pub-id-type="pmid">18601953</pub-id></element-citation></ref><ref id="B45"><label>45</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Pyrc</surname><given-names>K</given-names></name><name><surname>Jebbink</surname><given-names>MF</given-names></name><name><surname>Berkhout</surname><given-names>B</given-names></name><name><surname>Hoek</surname><given-names>L</given-names></name></person-group><source>SARS – and Other Coronaviruses</source><year>2008</year><fpage>1</fpage><lpage>17</lpage></element-citation></ref><ref id="B46"><label>46</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>de Vries</surname><given-names>M</given-names></name><name><surname>Pyrc</surname><given-names>K</given-names></name><name><surname>Berkhout</surname><given-names>R</given-names></name><name><surname>Vermeulen-Oost</surname><given-names>W</given-names></name><name><surname>Dijkman</surname><given-names>R</given-names></name><name><surname>Jebbink</surname><given-names>MF</given-names></name><name><surname>Bruisten</surname><given-names>S</given-names></name><name><surname>Berkhout</surname><given-names>B</given-names></name><name><surname>van der Hoek</surname><given-names>L</given-names></name></person-group><article-title>Human Parechovirus Type 1, 3, 4, 5, and 6 Detection in Picornavirus Cultures</article-title><source>J. Clin. Microbiol.</source><year>2008</year><volume>46</volume><fpage>759</fpage><lpage>762</lpage><pub-id pub-id-type="pmid">18077635</pub-id></element-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000F459 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000F459 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= MersV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |