Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin
Identifieur interne : 000993 ( Pmc/Corpus ); précédent : 000992; suivant : 000994Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin
Auteurs : Sami Samiullah ; Juliet Roberts ; Shu-Biao WuSource :
- PLoS ONE [ 1932-6203 ] ; 2017.
Abstract
Ten reference genes were investigated for normalization of gene expression data in the shell gland of laying hens. Analyses performed with geNorm revealed that hypoxanthine phosphoribosyltransferase 1 (
Url:
DOI: 10.1371/journal.pone.0180432
PubMed: 28671969
PubMed Central: 5495395
Links to Exploration step
PMC:5495395Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin</title><author><name sortKey="Samiullah, Sami" sort="Samiullah, Sami" uniqKey="Samiullah S" first="Sami" last="Samiullah">Sami Samiullah</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author><author><name sortKey="Roberts, Juliet" sort="Roberts, Juliet" uniqKey="Roberts J" first="Juliet" last="Roberts">Juliet Roberts</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author><author><name sortKey="Wu, Shu Biao" sort="Wu, Shu Biao" uniqKey="Wu S" first="Shu-Biao" last="Wu">Shu-Biao Wu</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">28671969</idno><idno type="pmc">5495395</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5495395</idno><idno type="RBID">PMC:5495395</idno><idno type="doi">10.1371/journal.pone.0180432</idno><date when="2017">2017</date><idno type="wicri:Area/Pmc/Corpus">000993</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000993</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin</title><author><name sortKey="Samiullah, Sami" sort="Samiullah, Sami" uniqKey="Samiullah S" first="Sami" last="Samiullah">Sami Samiullah</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author><author><name sortKey="Roberts, Juliet" sort="Roberts, Juliet" uniqKey="Roberts J" first="Juliet" last="Roberts">Juliet Roberts</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author><author><name sortKey="Wu, Shu Biao" sort="Wu, Shu Biao" uniqKey="Wu S" first="Shu-Biao" last="Wu">Shu-Biao Wu</name><affiliation><nlm:aff id="aff001"></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2017">2017</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Ten reference genes were investigated for normalization of gene expression data in the shell gland of laying hens. Analyses performed with geNorm revealed that hypoxanthine phosphoribosyltransferase 1 (<italic>HPRT1</italic>) and hydroxymethylbilane synthase (<italic>HMBS</italic>) were the two most stable reference genes in response to post-oviposition time alone (POT) or with nicarbazin treatment (POT+N) of laying hens. NormFinder analyses showed that the two most stable reference genes in response to POT and POT+N were 18S ribosomal RNA (<italic>18S rRNA</italic>), ribosomal protein L4 (<italic>RPL4</italic>) and <italic>HMBS</italic>, <italic>RPL4</italic>, respectively. BestKeeper analyses showed that <italic>18S rRNA</italic>, <italic>RPL4</italic> and <italic>HPRT1</italic>, <italic>HMBS</italic> were the two most stable reference genes for POT, and POT+N, respectively. Of the ten reference genes, all except <italic>B2M</italic> showed geNorm M <0.5, suggesting that they were stably expressed in the shell gland tissue. Consensus from these three programs suggested <italic>HPRT1</italic> and <italic>HMBS</italic> could be used as the two most stable reference genes in the present study. Expression analyses of four candidate target genes with the two most and the two least stable genes showed that a combination of stable reference genes leads to more discriminable quantification of expression levels of target genes, while the least stable genes failed to do so. Therefore, <italic>HMBS</italic> and <italic>HPRT1</italic> are recommended as the two most stable reference genes for the normalization of gene expression data at different stages of eggshell formation in brown-egg laying hens. Available statistical programs for reference gene ranking should include more robust analysis capability to analyse the gene expression data generated from factorial design experiments.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Nys, Y" uniqKey="Nys Y">Y Nys</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Melek, O" uniqKey="Melek O">O Melek</name></author><author><name sortKey="Morris, T" uniqKey="Morris T">T Morris</name></author><author><name sortKey="Jennings, R" uniqKey="Jennings R">R Jennings</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Johnson, Al" uniqKey="Johnson A">AL Johnson</name></author><author><name sortKey="Cc, S" uniqKey="Cc S">S CC</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Arias, Jl" uniqKey="Arias J">JL Arias</name></author><author><name sortKey="Fink, Dj" uniqKey="Fink D">DJ Fink</name></author><author><name sortKey="Xiao, S Q" uniqKey="Xiao S">S-Q Xiao</name></author><author><name sortKey="Heuer, Ah" uniqKey="Heuer A">AH Heuer</name></author><author><name sortKey="Caplan, Ai" uniqKey="Caplan A">AI Caplan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wellman Labadie, O" uniqKey="Wellman Labadie O">O Wellman-Labadie</name></author><author><name sortKey="Picman, J" uniqKey="Picman J">J Picman</name></author><author><name sortKey="Hincke, M" uniqKey="Hincke M">M Hincke</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Z" uniqKey="Liu Z">Z Liu</name></author><author><name sortKey="Zheng, Q" uniqKey="Zheng Q">Q Zheng</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author><author><name sortKey="Lu, L" uniqKey="Lu L">L Lu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jonchere, V" uniqKey="Jonchere V">V Jonchère</name></author><author><name sortKey="Rehault Godbert, S" uniqKey="Rehault Godbert S">S Réhault-Godbert</name></author><author><name sortKey="Hennequet Antier, C" uniqKey="Hennequet Antier C">C Hennequet-Antier</name></author><author><name sortKey="Cabau, C" uniqKey="Cabau C">C Cabau</name></author><author><name sortKey="Sibut, V" uniqKey="Sibut V">V Sibut</name></author><author><name sortKey="Cogburn, La" uniqKey="Cogburn L">LA Cogburn</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dunn, Ic" uniqKey="Dunn I">IC Dunn</name></author><author><name sortKey="Wilson, Pw" uniqKey="Wilson P">PW Wilson</name></author><author><name sortKey="Lu, Z" uniqKey="Lu Z">Z Lu</name></author><author><name sortKey="Bain, Mm" uniqKey="Bain M">MM Bain</name></author><author><name sortKey="Crossan, Cl" uniqKey="Crossan C">CL Crossan</name></author><author><name sortKey="Talbot, Rt" uniqKey="Talbot R">RT Talbot</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, G" uniqKey="Song G">G Song</name></author><author><name sortKey="Seo, Hw" uniqKey="Seo H">HW Seo</name></author><author><name sortKey="Choi, Jw" uniqKey="Choi J">JW Choi</name></author><author><name sortKey="Rengaraj, D" uniqKey="Rengaraj D">D Rengaraj</name></author><author><name sortKey="Kim, Tm" uniqKey="Kim T">TM Kim</name></author><author><name sortKey="Lee, Br" uniqKey="Lee B">BR Lee</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Additives, Epo" uniqKey="Additives E">EPo Additives</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hughes, B" uniqKey="Hughes B">B Hughes</name></author><author><name sortKey="Jones, J" uniqKey="Jones J">J Jones</name></author><author><name sortKey="Toler, J" uniqKey="Toler J">J Toler</name></author><author><name sortKey="Solis, J" uniqKey="Solis J">J Solis</name></author><author><name sortKey="Castaldo, D" uniqKey="Castaldo D">D Castaldo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dorne, J" uniqKey="Dorne J">J Dorne</name></author><author><name sortKey="Fernandez Cruz, M" uniqKey="Fernandez Cruz M">M Fernández-Cruz</name></author><author><name sortKey="Bertelsen, U" uniqKey="Bertelsen U">U Bertelsen</name></author><author><name sortKey="Renshaw, D" uniqKey="Renshaw D">D Renshaw</name></author><author><name sortKey="Peltonen, K" uniqKey="Peltonen K">K Peltonen</name></author><author><name sortKey="Anadon, A" uniqKey="Anadon A">A Anadon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bilandzi, N" uniqKey="Bilandzi N">N Bilandžić</name></author><author><name sortKey="Dolenc, J" uniqKey="Dolenc J">J Dolenc</name></author><author><name sortKey="Ga Nik, Ks" uniqKey="Ga Nik K">KŠ Gačnik</name></author><author><name sortKey="Varenina, I" uniqKey="Varenina I">I Varenina</name></author><author><name sortKey="Kolanovi, Bs" uniqKey="Kolanovi B">BS Kolanović</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Leeson, S" uniqKey="Leeson S">S Leeson</name></author><author><name sortKey="Caston, L" uniqKey="Caston L">L Caston</name></author><author><name sortKey="Summers, J" uniqKey="Summers J">J Summers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Thellin, O" uniqKey="Thellin O">O Thellin</name></author><author><name sortKey="Zorzi, W" uniqKey="Zorzi W">W Zorzi</name></author><author><name sortKey="Lakaye, B" uniqKey="Lakaye B">B Lakaye</name></author><author><name sortKey="De Borman, B" uniqKey="De Borman B">B De Borman</name></author><author><name sortKey="Coumans, B" uniqKey="Coumans B">B Coumans</name></author><author><name sortKey="Hennen, G" uniqKey="Hennen G">G Hennen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lee, Pd" uniqKey="Lee P">PD Lee</name></author><author><name sortKey="Sladek, R" uniqKey="Sladek R">R Sladek</name></author><author><name sortKey="Greenwood, Cm" uniqKey="Greenwood C">CM Greenwood</name></author><author><name sortKey="Hudson, Tj" uniqKey="Hudson T">TJ Hudson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Greer, S" uniqKey="Greer S">S Greer</name></author><author><name sortKey="Honeywell, R" uniqKey="Honeywell R">R Honeywell</name></author><author><name sortKey="Geletu, M" uniqKey="Geletu M">M Geletu</name></author><author><name sortKey="Arulanandam, R" uniqKey="Arulanandam R">R Arulanandam</name></author><author><name sortKey="Raptis, L" uniqKey="Raptis L">L Raptis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Uddin, Mj" uniqKey="Uddin M">MJ Uddin</name></author><author><name sortKey="Cinar, Mu" uniqKey="Cinar M">MU Cinar</name></author><author><name sortKey="Tesfaye, D" uniqKey="Tesfaye D">D Tesfaye</name></author><author><name sortKey="Looft, C" uniqKey="Looft C">C Looft</name></author><author><name sortKey="Tholen, E" uniqKey="Tholen E">E Tholen</name></author><author><name sortKey="Schellander, K" uniqKey="Schellander K">K Schellander</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cinar, Mu" uniqKey="Cinar M">MU Cinar</name></author><author><name sortKey="Islam, Ma" uniqKey="Islam M">MA Islam</name></author><author><name sortKey="Uddin, Mj" uniqKey="Uddin M">MJ Uddin</name></author><author><name sortKey="Tholen, E" uniqKey="Tholen E">E Tholen</name></author><author><name sortKey="Tesfaye, D" uniqKey="Tesfaye D">D Tesfaye</name></author><author><name sortKey="Looft, C" uniqKey="Looft C">C Looft</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Beekman, L" uniqKey="Beekman L">L Beekman</name></author><author><name sortKey="Tohver, T" uniqKey="Tohver T">T Tohver</name></author><author><name sortKey="Dardari, R" uniqKey="Dardari R">R Dardari</name></author><author><name sortKey="Leguillette, R" uniqKey="Leguillette R">R Léguillette</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Petriccione, M" uniqKey="Petriccione M">M Petriccione</name></author><author><name sortKey="Mastrobuoni, F" uniqKey="Mastrobuoni F">F Mastrobuoni</name></author><author><name sortKey="Zampella, L" uniqKey="Zampella L">L Zampella</name></author><author><name sortKey="Scortichini, M" uniqKey="Scortichini M">M Scortichini</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vandesompele, J" uniqKey="Vandesompele J">J Vandesompele</name></author><author><name sortKey="De Preter, K" uniqKey="De Preter K">K De Preter</name></author><author><name sortKey="Pattyn, F" uniqKey="Pattyn F">F Pattyn</name></author><author><name sortKey="Poppe, B" uniqKey="Poppe B">B Poppe</name></author><author><name sortKey="Van Roy, N" uniqKey="Van Roy N">N Van Roy</name></author><author><name sortKey="De Paepe, A" uniqKey="De Paepe A">A De Paepe</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Andersen, Cl" uniqKey="Andersen C">CL Andersen</name></author><author><name sortKey="Jensen, Jl" uniqKey="Jensen J">JL Jensen</name></author><author><name sortKey=" Rntoft, Tf" uniqKey=" Rntoft T">TF Ørntoft</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pfaffl, Mw" uniqKey="Pfaffl M">MW Pfaffl</name></author><author><name sortKey="Tichopad, A" uniqKey="Tichopad A">A Tichopad</name></author><author><name sortKey="Prgomet, C" uniqKey="Prgomet C">C Prgomet</name></author><author><name sortKey="Neuvians, Tp" uniqKey="Neuvians T">TP Neuvians</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guernsey, Dl" uniqKey="Guernsey D">DL Guernsey</name></author><author><name sortKey="Jiang, H" uniqKey="Jiang H">H Jiang</name></author><author><name sortKey="Campagna, Dr" uniqKey="Campagna D">DR Campagna</name></author><author><name sortKey="Evans, Sc" uniqKey="Evans S">SC Evans</name></author><author><name sortKey="Ferguson, M" uniqKey="Ferguson M">M Ferguson</name></author><author><name sortKey="Kellogg, Md" uniqKey="Kellogg M">MD Kellogg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="May, Bk" uniqKey="May B">BK May</name></author><author><name sortKey="Dogra, Sc" uniqKey="Dogra S">SC Dogra</name></author><author><name sortKey="Sadlon, Tj" uniqKey="Sadlon T">TJ Sadlon</name></author><author><name sortKey="Bhasker, Cr" uniqKey="Bhasker C">CR Bhasker</name></author><author><name sortKey="Cox, Tc" uniqKey="Cox T">TC Cox</name></author><author><name sortKey="Bottomley, Ss" uniqKey="Bottomley S">SS Bottomley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Takahashi, S" uniqKey="Takahashi S">S Takahashi</name></author><author><name sortKey="Taketani, S" uniqKey="Taketani S">S Taketani</name></author><author><name sortKey="Akasaka, J E" uniqKey="Akasaka J">J-e Akasaka</name></author><author><name sortKey="Kobayashi, A" uniqKey="Kobayashi A">A Kobayashi</name></author><author><name sortKey="Hayashi, N" uniqKey="Hayashi N">N Hayashi</name></author><author><name sortKey="Yamamoto, M" uniqKey="Yamamoto M">M Yamamoto</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Day, Al" uniqKey="Day A">AL Day</name></author><author><name sortKey="Parsons, Bm" uniqKey="Parsons B">BM Parsons</name></author><author><name sortKey="Dailey, Ha" uniqKey="Dailey H">HA Dailey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harbin, Bm" uniqKey="Harbin B">BM Harbin</name></author><author><name sortKey="Dailey, Ha" uniqKey="Dailey H">HA Dailey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kuchipudi, Sv" uniqKey="Kuchipudi S">SV Kuchipudi</name></author><author><name sortKey="Tellabati, M" uniqKey="Tellabati M">M Tellabati</name></author><author><name sortKey="Sebastian, S" uniqKey="Sebastian S">S Sebastian</name></author><author><name sortKey="Londt, Bz" uniqKey="Londt B">BZ Londt</name></author><author><name sortKey="Jansen, C" uniqKey="Jansen C">C Jansen</name></author><author><name sortKey="Vervelde, L" uniqKey="Vervelde L">L Vervelde</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, F" uniqKey="Yang F">F Yang</name></author><author><name sortKey="Lei, X" uniqKey="Lei X">X Lei</name></author><author><name sortKey="Rodriguez Palacios, A" uniqKey="Rodriguez Palacios A">A Rodriguez-Palacios</name></author><author><name sortKey="Tang, C" uniqKey="Tang C">C Tang</name></author><author><name sortKey="Yue, H" uniqKey="Yue H">H Yue</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dalgaard, Ts" uniqKey="Dalgaard T">TS Dalgaard</name></author><author><name sortKey="Skovgaard, K" uniqKey="Skovgaard K">K Skovgaard</name></author><author><name sortKey="Norup, Lr" uniqKey="Norup L">LR Norup</name></author><author><name sortKey="Pleidrup, J" uniqKey="Pleidrup J">J Pleidrup</name></author><author><name sortKey="Permin, A" uniqKey="Permin A">A Permin</name></author><author><name sortKey="Schou, Tw" uniqKey="Schou T">TW Schou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yin, R" uniqKey="Yin R">R Yin</name></author><author><name sortKey="Liu, X" uniqKey="Liu X">X Liu</name></author><author><name sortKey="Liu, C" uniqKey="Liu C">C Liu</name></author><author><name sortKey="Ding, Z" uniqKey="Ding Z">Z Ding</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author><author><name sortKey="Tian, F" uniqKey="Tian F">F Tian</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, G" uniqKey="Li G">G Li</name></author><author><name sortKey="Chen, S" uniqKey="Chen S">S Chen</name></author><author><name sortKey="Duan, Z" uniqKey="Duan Z">Z Duan</name></author><author><name sortKey="Qu, L" uniqKey="Qu L">L Qu</name></author><author><name sortKey="Xu, G" uniqKey="Xu G">G Xu</name></author><author><name sortKey="Yang, N" uniqKey="Yang N">N Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Samiullah, S" uniqKey="Samiullah S">S Samiullah</name></author><author><name sortKey="Roberts, J" uniqKey="Roberts J">J Roberts</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Hellemans, J" uniqKey="Hellemans J">J Hellemans</name></author><author><name sortKey="Mortier, G" uniqKey="Mortier G">G Mortier</name></author><author><name sortKey="De Paepe, A" uniqKey="De Paepe A">A De Paepe</name></author><author><name sortKey="Speleman, F" uniqKey="Speleman F">F Speleman</name></author><author><name sortKey="Vandesompele, J" uniqKey="Vandesompele J">J Vandesompele</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Radoni, A" uniqKey="Radoni A">A Radonić</name></author><author><name sortKey="Thulke, S" uniqKey="Thulke S">S Thulke</name></author><author><name sortKey="Bae, H G" uniqKey="Bae H">H-G Bae</name></author><author><name sortKey="Muller, Ma" uniqKey="Muller M">MA Müller</name></author><author><name sortKey="Siegert, W" uniqKey="Siegert W">W Siegert</name></author><author><name sortKey="Nitsche, A" uniqKey="Nitsche A">A Nitsche</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pfaffl, Mw" uniqKey="Pfaffl M">MW Pfaffl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ji, H" uniqKey="Ji H">H Ji</name></author><author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name></author><author><name sortKey="Liu, J" uniqKey="Liu J">J Liu</name></author><author><name sortKey="Guo, J" uniqKey="Guo J">J Guo</name></author><author><name sortKey="Wang, Z" uniqKey="Wang Z">Z Wang</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goossens, K" uniqKey="Goossens K">K Goossens</name></author><author><name sortKey="Van Poucke, M" uniqKey="Van Poucke M">M Van Poucke</name></author><author><name sortKey="Van Soom, A" uniqKey="Van Soom A">A Van Soom</name></author><author><name sortKey="Vandesompele, J" uniqKey="Vandesompele J">J Vandesompele</name></author><author><name sortKey="Van Zeveren, A" uniqKey="Van Zeveren A">A Van Zeveren</name></author><author><name sortKey="Peelman, Lj" uniqKey="Peelman L">LJ Peelman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Svingen, T" uniqKey="Svingen T">T Svingen</name></author><author><name sortKey="Letting, H" uniqKey="Letting H">H Letting</name></author><author><name sortKey="Hadrup, N" uniqKey="Hadrup N">N Hadrup</name></author><author><name sortKey="Hass, U" uniqKey="Hass U">U Hass</name></author><author><name sortKey="Vinggaard, Am" uniqKey="Vinggaard A">AM Vinggaard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nascimento, Cs" uniqKey="Nascimento C">CS Nascimento</name></author><author><name sortKey="Barbosa, Lt" uniqKey="Barbosa L">LT Barbosa</name></author><author><name sortKey="Brito, C" uniqKey="Brito C">C Brito</name></author><author><name sortKey="Fernandes, Rp" uniqKey="Fernandes R">RP Fernandes</name></author><author><name sortKey="Mann, Rs" uniqKey="Mann R">RS Mann</name></author><author><name sortKey="Pinto, Apg" uniqKey="Pinto A">APG Pinto</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cinar, Mu" uniqKey="Cinar M">MU Cinar</name></author><author><name sortKey="Islam, Ma" uniqKey="Islam M">MA Islam</name></author><author><name sortKey="Proll, M" uniqKey="Proll M">M Pröll</name></author><author><name sortKey="Kocamis, H" uniqKey="Kocamis H">H Kocamis</name></author><author><name sortKey="Tholen, E" uniqKey="Tholen E">E Tholen</name></author><author><name sortKey="Tesfaye, D" uniqKey="Tesfaye D">D Tesfaye</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pereira Fantini, Pm" uniqKey="Pereira Fantini P">PM Pereira-Fantini</name></author><author><name sortKey="Rajapaksa, Ae" uniqKey="Rajapaksa A">AE Rajapaksa</name></author><author><name sortKey="Oakley, R" uniqKey="Oakley R">R Oakley</name></author><author><name sortKey="Tingay, Dg" uniqKey="Tingay D">DG Tingay</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shivhare, R" uniqKey="Shivhare R">R Shivhare</name></author><author><name sortKey="Lata, C" uniqKey="Lata C">C Lata</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Z J" uniqKey="Wu Z">Z-J Wu</name></author><author><name sortKey="Tian, C" uniqKey="Tian C">C Tian</name></author><author><name sortKey="Jiang, Q" uniqKey="Jiang Q">Q Jiang</name></author><author><name sortKey="Li, X H" uniqKey="Li X">X-H Li</name></author><author><name sortKey="Zhuang, J" uniqKey="Zhuang J">J Zhuang</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">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, CA USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">28671969</article-id><article-id pub-id-type="pmc">5495395</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0180432</article-id><article-id pub-id-type="publisher-id">PONE-D-16-45782</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Genetics</subject><subj-group><subject>Gene Expression</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Biochemistry</subject><subj-group><subject>Nucleic acids</subject><subj-group><subject>RNA</subject><subj-group><subject>Non-coding RNA</subject><subj-group><subject>Ribosomal RNA</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Biochemistry</subject><subj-group><subject>Ribosomes</subject><subj-group><subject>Ribosomal RNA</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Cell biology</subject><subj-group><subject>Cellular structures and organelles</subject><subj-group><subject>Ribosomes</subject><subj-group><subject>Ribosomal RNA</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Computer and Information Sciences</subject><subj-group><subject>Computer Software</subject></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Animals</subject><subj-group><subject>Vertebrates</subject><subj-group><subject>Amniotes</subject><subj-group><subject>Birds</subject><subj-group><subject>Fowl</subject><subj-group><subject>Gamefowl</subject><subj-group><subject>Chickens</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Animals</subject><subj-group><subject>Vertebrates</subject><subj-group><subject>Amniotes</subject><subj-group><subject>Birds</subject><subj-group><subject>Poultry</subject><subj-group><subject>Chickens</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Genetics</subject><subj-group><subject>Animal Genetics</subject><subj-group><subject>Bird Genetics</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Molecular Biology</subject><subj-group><subject>Molecular Biology Techniques</subject><subj-group><subject>Artificial Gene Amplification and Extension</subject><subj-group><subject>Polymerase Chain Reaction</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Research and Analysis Methods</subject><subj-group><subject>Molecular Biology Techniques</subject><subj-group><subject>Artificial Gene Amplification and Extension</subject><subj-group><subject>Polymerase Chain Reaction</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Cell Biology</subject><subj-group><subject>Cellular Structures and Organelles</subject><subj-group><subject>Cell Membranes</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Genetics</subject><subj-group><subject>Gene Amplification</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin</article-title><alt-title alt-title-type="running-head">Reference gene selection for the shell gland of laying hens</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Samiullah</surname><given-names>Sami</given-names></name><xref ref-type="aff" rid="aff001"></xref></contrib><contrib contrib-type="author"><name><surname>Roberts</surname><given-names>Juliet</given-names></name><xref ref-type="aff" rid="aff001"></xref></contrib><contrib contrib-type="author"><contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0002-1790-6015</contrib-id><name><surname>Wu</surname><given-names>Shu-Biao</given-names></name><xref ref-type="aff" rid="aff001"></xref><xref ref-type="corresp" rid="cor001">*</xref></contrib></contrib-group><aff id="aff001"><addr-line>Animal Science, School of Environmental and Rural Science, University of New England, Armidale, New South Wales, Australia</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Woloschak</surname><given-names>Gayle E.</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>Northwestern University Feinberg School of Medicine, UNITED STATES</addr-line></aff><author-notes><fn fn-type="COI-statement" id="coi001"><p><bold>Competing Interests: </bold>This study was funded by Australian Egg Corporation Limited (AECL). AECL has made no attempt to influence the reporting of the results of the present study and no aspect of the findings of the present study has any direct bearing on AECL commercial activities. This does not alter our adherence to PLOS ONE policies on sharing data and materials. AECL is one of 15 Rural Research and Development Corporations (RDCs) which receives funding from the Australian Government as well as from industry levies (see <ext-link ext-link-type="uri" xlink:href="http://www.ruralrdc.com.au/rural-research-development-corporations/">http://www.ruralrdc.com.au/rural-research-development-corporations/</ext-link>). The AECL R&D Program is administered independently from the commercial aspects of the AECL as required as a condition of Australian Government funding (see <ext-link ext-link-type="uri" xlink:href="http://www.ruralrdc.com.au/abouttherrdcs/#governance">http://www.ruralrdc.com.au/abouttherrdcs/#governance</ext-link>). SBW is currently serving on the editorial board of PLOS ONE. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.</p></fn><fn fn-type="con"><p><list list-type="simple"><list-item><p><bold>Conceptualization:</bold> SS JR SBW.</p></list-item><list-item><p><bold>Data curation:</bold> SS SBW.</p></list-item><list-item><p><bold>Formal analysis:</bold> SS SBW.</p></list-item><list-item><p><bold>Funding acquisition:</bold> JR.</p></list-item><list-item><p><bold>Investigation:</bold> SS JR SBW.</p></list-item><list-item><p><bold>Methodology:</bold> SS JR SBW.</p></list-item><list-item><p><bold>Project administration:</bold> JR SBW.</p></list-item><list-item><p><bold>Resources:</bold> JR SBW.</p></list-item><list-item><p><bold>Supervision:</bold> SBW JR.</p></list-item><list-item><p><bold>Visualization:</bold> SS.</p></list-item><list-item><p><bold>Writing – original draft:</bold> SS SBW.</p></list-item><list-item><p><bold>Writing – review & editing:</bold> SBW JR.</p></list-item></list></p></fn><corresp id="cor001">* E-mail: <email>shubiao.wu@une.edu.au</email></corresp></author-notes><pub-date pub-type="epub"><day>3</day><month>7</month><year>2017</year></pub-date><pub-date pub-type="collection"><year>2017</year></pub-date><volume>12</volume><issue>7</issue><elocation-id>e0180432</elocation-id><history><date date-type="received"><day>17</day><month>11</month><year>2016</year></date><date date-type="accepted"><day>15</day><month>6</month><year>2017</year></date></history><permissions><copyright-statement>© 2017 Samiullah et al</copyright-statement><copyright-year>2017</copyright-year><copyright-holder>Samiullah et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="pone.0180432.pdf"></self-uri><abstract><p>Ten reference genes were investigated for normalization of gene expression data in the shell gland of laying hens. Analyses performed with geNorm revealed that hypoxanthine phosphoribosyltransferase 1 (<italic>HPRT1</italic>) and hydroxymethylbilane synthase (<italic>HMBS</italic>) were the two most stable reference genes in response to post-oviposition time alone (POT) or with nicarbazin treatment (POT+N) of laying hens. NormFinder analyses showed that the two most stable reference genes in response to POT and POT+N were 18S ribosomal RNA (<italic>18S rRNA</italic>), ribosomal protein L4 (<italic>RPL4</italic>) and <italic>HMBS</italic>, <italic>RPL4</italic>, respectively. BestKeeper analyses showed that <italic>18S rRNA</italic>, <italic>RPL4</italic> and <italic>HPRT1</italic>, <italic>HMBS</italic> were the two most stable reference genes for POT, and POT+N, respectively. Of the ten reference genes, all except <italic>B2M</italic> showed geNorm M <0.5, suggesting that they were stably expressed in the shell gland tissue. Consensus from these three programs suggested <italic>HPRT1</italic> and <italic>HMBS</italic> could be used as the two most stable reference genes in the present study. Expression analyses of four candidate target genes with the two most and the two least stable genes showed that a combination of stable reference genes leads to more discriminable quantification of expression levels of target genes, while the least stable genes failed to do so. Therefore, <italic>HMBS</italic> and <italic>HPRT1</italic> are recommended as the two most stable reference genes for the normalization of gene expression data at different stages of eggshell formation in brown-egg laying hens. Available statistical programs for reference gene ranking should include more robust analysis capability to analyse the gene expression data generated from factorial design experiments.</p></abstract><funding-group><award-group id="award001"><funding-source><institution-wrap><institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100000976</institution-id><institution>Australian Egg Corporation Limited</institution></institution-wrap></funding-source><award-id>AECL1UN121</award-id><principal-award-recipient><name><surname>Roberts</surname><given-names>Juliet</given-names></name></principal-award-recipient></award-group><funding-statement>This study was financially supported by Australian Egg Corporation Limited (AECL) under grant number AECL 1UN121. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><fig-count count="5"></fig-count><table-count count="8"></table-count><page-count count="20"></page-count></counts><custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name><meta-value>All relevant data are within the paper. We do not have any supporting information in this paper.</meta-value></custom-meta></custom-meta-group></article-meta><notes><title>Data Availability</title><p>All relevant data are within the paper. We do not have any supporting information in this paper.</p></notes></front><body><sec sec-type="intro" id="sec001"><title>Introduction</title><p>The chicken reproductive tract is divided into five histologically distinct parts, the ovary, infundibulum, magnum, isthmus, and shell gland (uterus). The shell gland is an expanded pouch-like part of the oviduct where an egg remains for approximately 18–20 hours, during which shell formation takes place [<xref rid="pone.0180432.ref001" ref-type="bibr">1</xref>]. The next ovulation occurs 0.5 hour after the preceding oviposition [<xref rid="pone.0180432.ref002" ref-type="bibr">2</xref>]. Calcification of the eggshell is associated with stimuli initiated by ovulation or by neuroendocrine factors that control and coordinate both ovulation and calcium secretion [<xref rid="pone.0180432.ref003" ref-type="bibr">3</xref>]. The eggshell is a highly ordered bio-ceramic (about 5 g calcium) of fused calcite crystalline cones, formed on a protein skeleton with distinct layers and regular pores [<xref rid="pone.0180432.ref004" ref-type="bibr">4</xref>]. Like other epithelial cells, the shell gland epithelium provides antimicrobial protection for both hen and egg; thus, it is a rich source of anti-microbial proteins [<xref rid="pone.0180432.ref005" ref-type="bibr">5</xref>].</p><p>Approximately 437 peptides and ion transporters have been identified as being involved in the formation of three distinct layers of the eggshell [<xref rid="pone.0180432.ref006" ref-type="bibr">6</xref>, <xref rid="pone.0180432.ref007" ref-type="bibr">7</xref>]. However, the mechanisms of egg and eggshell biogenesis in relation to the origin and flow of the precursors at various stages of egg formation are not fully understood. The biosynthetic pathway of protoporphyrin IX (PP IX) is important in the formation of the eggshell as it contributes to the shell colour. PP IX is an immediate precursor of heme and a major component of brown eggshell pigment. To date, there is little information about the origin of its precursors and how PP IX is inhibited from converting into heme through the enzymatic activity of ferrochelatase in the shell gland of laying hens. Furthermore, it is not clear how many genes are involved in the PP IX synthesis, ultimate transportation across cell membranes and subsequent deposition into distinct eggshell layers. In the shell gland, hundreds of genes that are differentially expressed between juvenile and mature laying hens have been identified [<xref rid="pone.0180432.ref008" ref-type="bibr">8</xref>]. It is assumed that some genes express differentially in relation to the formation of distinct layers of the eggshell. Furthermore, the expression of genes associated with epithelial differentiation and tissue remodelling may vary with different levels of estrogen secretion in the presence or absence of an egg [<xref rid="pone.0180432.ref009" ref-type="bibr">9</xref>].</p><p>Nicarbazin is one of the various factors that causes lower production and/or deposition of PP IX into eggshell, when fed at recommended dosages (50–125 mg/kg of feed) to brown-egg laying hens [<xref rid="pone.0180432.ref010" ref-type="bibr">10</xref>, <xref rid="pone.0180432.ref011" ref-type="bibr">11</xref>]. It is a chemically produced drug composed of a complex equimolar amount of 4,40-dinitrocarbanilide (DNC) and 2-hydroxy-4,6-dimethylpyrimidine (HDP) and is registered for use in poultry fattening and for treating coccidiosis [<xref rid="pone.0180432.ref010" ref-type="bibr">10</xref>, <xref rid="pone.0180432.ref012" ref-type="bibr">12</xref>]. It is used either in pure or combined forms as a feed additive in poultry production [<xref rid="pone.0180432.ref013" ref-type="bibr">13</xref>]. Its effects on eggshell colour, which are dosage and time dependent, are reversible [<xref rid="pone.0180432.ref014" ref-type="bibr">14</xref>]. The pharmacodynamics of nicarbazin in the shell gland are not yet known, and the way it acts at the molecular level to alter the synthesis and/or deposition of PP IX into eggshells at different stages of shell formation needs to be further investigated. Therefore, we anticipate that the responses of the genes involved in the PP IX synthesis pathway to nicarbazin will shed some light on its actions in the alteration of eggshell colour deposition. Nicarbazin was used as a model as residues of this drug in feed are still a problem of loss in shell colour particularly in Australia. In addition, investigation of the molecular basis of nicarbazin effects on shell colour may shed light on the mechanisms of shell colour loss in commercial laying flocks from other causes.</p><p>The transcriptional profiling in the shell gland is different from other parts of the oviduct, such as the magnum and isthmus [<xref rid="pone.0180432.ref007" ref-type="bibr">7</xref>]. Thus, data on differential gene expression in the shell gland are needed and, therefore, selection of reference genes for gene expression analysis as a normalization approach is paramount. To the best of our knowledge, no studies have been performed to identify suitable reference genes for the normalization of quantitative PCR (qPCR) gene expression data in the shell gland of chickens. Traditionally, the most commonly used housekeeping genes, such as <italic>ACTB</italic>, <italic>TUBB</italic>, and <italic>GAPDH</italic> have been selected as generic reference genes. However, ample evidence has shown that the expression of these genes may not be constant across a range of experimental conditions and tissues under investigation [<xref rid="pone.0180432.ref015" ref-type="bibr">15</xref>–<xref rid="pone.0180432.ref020" ref-type="bibr">20</xref>]. Thus, it is now recommended to use these genes only as reference genes for normalization when prior analysis of their expression stability has been carried out [<xref rid="pone.0180432.ref021" ref-type="bibr">21</xref>], to ensure that cellular expression level of the reference genes is virtually identical under different conditions in the study. It is also recommended that more than one reference gene be employed to achieve more robust, accurate and reliable normalization of gene expression data [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>].</p><p>Of the available statistical software, three distinct tools have been frequently reported in the literature for ranking the overall expression stabilities of the reference genes as normalizers for gene expression studies. The geNorm module in qbase+ software version 3.0 (Biogazelle, Belgium) calculates the gene expression stability (geNorm M) as the arithmetic mean of the pairwise variation (geNorm V) between all tested genes [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>, <xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>]. The geNorm V for any given two genes is the standard deviation calculated from the log<sub>2</sub> transformed relative quantities between those two genes [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>]. Before analysis, qbase+ pre-processes the data for efficiency correction, inter-run calibration, bad replicates removal and conversion to relative quantities [<xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>]. The relative quantities are then converted to either linear or log transformed scale. The current geNorm tool does limit the minimum number of genes to 8, unlike to its previous Excel based version. qbase+ allows easy exchange of data between users, and exports tabulated data for further statistical analysis using dedicated software. The most stably expressed gene produces the lowest geNorm M value. The most stable reference gene is determined by step-wise exclusion of the least stable genes [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>]. geNorm eliminates the genes sequentially and thus a differentially expressed gene does not affect the ultimate outcomes from the analysis. Therefore, geNorm is usually less sensitive to differentially expressed genes initially included in the assay. Good reference genes have an M < 0.5 and CV (Coefficient of variance) <0.2, while M values up to 1 are acceptable for more difficult samples [<xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>]. The cut-off value for geNorm V is 0.15. geNorm does not consider treatment groups and all samples are treated as being from a single population. NormFinder (GenEx version 6.0.1) calculates the standard deviation (SD) of the genes relative to the mean expression of all the genes in the panel [<xref rid="pone.0180432.ref024" ref-type="bibr">24</xref>]. It calculates a global average expression of all the genes in all the samples, to which the individual genes are compared. Based on this comparison, SD for each reference gene is estimated. Furthermore, if the samples are from different treatment groups, NormFinder separates the variation into an intragroup and an intergroup contribution [<xref rid="pone.0180432.ref024" ref-type="bibr">24</xref>, <xref rid="pone.0180432.ref025" ref-type="bibr">25</xref>]. Hence, a low stability value reflects low inter- and intra-group variation. Similar to qbase+, GenEx does have an option to highlight bad replicates during data analysis and thus bad replicates can be excluded from the analysis easily. In the GenEx version 6.0.1, the data pre-processing is very similar to that explained in the qbase+ section. In addition, the data in GenEx, can be also converted to logarithmic scale, such as log<sub>2</sub>, log<sub>10</sub>, 1n and log(X+1). An Excel based BestKeeper (Version 1) software is used to determine the best stable reference genes based on Pearson correlation coefficient (r), coefficient of variance (CV) and standard deviation (SD). Only genes with a high r and low SD values are combined into BestKeeper index (BKI) value using the geometric mean of their Cq values [<xref rid="pone.0180432.ref026" ref-type="bibr">26</xref>]. The BKI is calculated from the geometric mean of the candidates Cq values for each specific sample [<xref rid="pone.0180432.ref026" ref-type="bibr">26</xref>]. The most stable reference genes are the ones with the lowest SD values and highest coefficients of correlation with the BKI [<xref rid="pone.0180432.ref026" ref-type="bibr">26</xref>]. BestKeeper also uses a statistical algorithm wherein the Pearson correlation coefficient for each candidate reference gene pair is calculated along with the probability of correlation significance of the pair [<xref rid="pone.0180432.ref026" ref-type="bibr">26</xref>]. Overall, a generalized opinion from the literature and scientific forums can be summarised: geNorm, NormFinder and BestKeeper basically provide similar outcomes for the overall stability of candidate reference genes.</p><p>In the present study, we aimed to select reference genes from ten housekeeping genes to be used for the analysis of gene expression levels at different stages of eggshell formation in the shell gland and in response to nicarbazin feeding of the laying hens. The most stable genes were selected on the basis of the stability of the genes across the three software. Furthermore, four candidate target genes encoding either enzymes or peptide transporter were chosen to compare the outcomes from the data by the two most and the two least stable reference genes. The solute carrier family 25, member 38 (<italic>SLC25A38</italic>), located on mitochondrial membrane, transports glycine into mitochondria for the synthesis of aminolevulinic acid, the first step in the synthesis of PP IX [<xref rid="pone.0180432.ref027" ref-type="bibr">27</xref>]. The delta-aminolevulinate synthase 1 (<italic>ALAS1</italic>) gene encodes a rate limiting non-erythroid enzyme that catalyses the reaction of succinyl co-enzyme A with glycine to form delta-aminolevulinic acid within the mitochondrial matrix [<xref rid="pone.0180432.ref028" ref-type="bibr">28</xref>]. Coproporphyrinogen oxidase (<italic>CPOX</italic>) gene encodes an enzyme in the PP IX biosynthetic pathway that converts coproporphyrinogen III into protoporhyrinogen III [<xref rid="pone.0180432.ref029" ref-type="bibr">29</xref>]. Ferrochelatase (<italic>FECH</italic>) gene encodes <italic>FECH</italic> enzyme, which converts PP IX into heme [<xref rid="pone.0180432.ref030" ref-type="bibr">30</xref>, <xref rid="pone.0180432.ref031" ref-type="bibr">31</xref>]. The outcome of the study provides a set of reference genes that are expressed in a relatively constant level for all the birds sampled at different time-points of eggshell formation and in response to treatment with nicarbazin.</p></sec><sec sec-type="materials|methods" id="sec002"><title>Materials and methods</title><p>The experimental setup was approved by the University of New England, Animal Ethics Approval Committee under Authority No. AEC15-022. The protocol was carried out in accordance with the guidelines specified in the Australian Code for the Care and Use of Animals for Scientific Purposes 8<sup>th</sup> edition 2013.</p><sec id="sec003"><title>Selection of reference genes and primer design</title><p>In the current study, ten reference genes were selected from the literature published for chickens and other animals (<xref ref-type="table" rid="pone.0180432.t001">Table 1</xref>). The primers were either sourced from previously published studies in chickens or designed using NCBI primer tool (<xref ref-type="table" rid="pone.0180432.t002">Table 2</xref>). The primer quality was checked in “Beacon Designer” software (<ext-link ext-link-type="uri" xlink:href="http://www.premierbiosoft.com/qOligo/Oligo.jsp?PID=1">http://www.premierbiosoft.com/qOligo/Oligo.jsp?PID=1</ext-link>) for the levels of secondary structures such as primer dimer, sequence repeats and palindrome. To check the sequence specificity, primers were blasted against the NCBI database using BLASTN, Ensemble Chicken Galgal4 and UCSC’s Chicken (<italic>Gallus gallus</italic>) Genome Browser Gateway. Prior to qPCR analysis, primer efficiency and specificity for each primer pair were examined with target RNA samples in 10-time serial dilutions. Only primer pairs with specific amplifications and high efficiency were used in the optimisation.</p><table-wrap id="pone.0180432.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t001</object-id><label>Table 1</label><caption><title>Functional annotations of the reference genes used in the current study.</title></caption><alternatives><graphic id="pone.0180432.t001g" xlink:href="pone.0180432.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Gene symbol</th><th align="left" rowspan="1" colspan="1">Description</th><th align="left" rowspan="1" colspan="1">Cellular localization</th><th align="left" rowspan="1" colspan="1">Biological function</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><bold><italic>18S rRNA</italic></bold></td><td align="left" rowspan="1" colspan="1">Nuclear ribosomal RNA small subunit</td><td align="left" rowspan="1" colspan="1">cytoplasm, nucleus</td><td align="left" rowspan="1" colspan="1">biogenesis and export of the 40S ribosomal subunit</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>ALB</italic></bold></td><td align="left" rowspan="1" colspan="1">Albumin</td><td align="left" rowspan="1" colspan="1">extracellular</td><td align="left" rowspan="1" colspan="1">stabilizing extracellular fluid volume</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>ACTB</italic></bold></td><td align="left" rowspan="1" colspan="1">β-actin</td><td align="left" rowspan="1" colspan="1">cytoplasm, membranes</td><td align="left" rowspan="1" colspan="1">cytoskeletal structural protein, nucleotide, and ATP binding</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>B2M</italic></bold></td><td align="left" rowspan="1" colspan="1">Beta 2-microglobulin</td><td align="left" rowspan="1" colspan="1">Golgi membrane, plasma membrane, early endosome membrane, extracellular region</td><td align="left" rowspan="1" colspan="1">cytoskeletal protein, immune response, protein binding</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>CA2</italic></bold></td><td align="left" rowspan="1" colspan="1">Carbonic anhydrase 2</td><td align="left" rowspan="1" colspan="1">cytoplasm, cell membrane</td><td align="left" rowspan="1" colspan="1">catalyses reversible hydration of carbon dioxide</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>CST3</italic></bold></td><td align="left" rowspan="1" colspan="1">Cystatin C</td><td align="left" rowspan="1" colspan="1">extracellular</td><td align="left" rowspan="1" colspan="1">as an inhibitor of cysteine proteinases</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>GAPDH</italic></bold></td><td align="left" rowspan="1" colspan="1">Glyceraldehyde 3-phosphate dehydrogenase</td><td align="left" rowspan="1" colspan="1">plasma membrane</td><td align="left" rowspan="1" colspan="1">glycolytic enzyme, oxidoreductase in glycolysis and gluconeogenesis</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>HMBS</italic></bold></td><td align="left" rowspan="1" colspan="1">Hydroxymethylbilane synthase</td><td align="left" rowspan="1" colspan="1">cytoplasm</td><td align="left" rowspan="1" colspan="1">heme synthesis, porphyrin metabolism, transferase activity</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>HPRT1</italic></bold></td><td align="left" rowspan="1" colspan="1">Hypoxanthine phosphoribosyltransferase 1</td><td align="left" rowspan="1" colspan="1">cytoplasm</td><td align="left" rowspan="1" colspan="1">purine synthesis in salvage pathway</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold><italic>RPL4</italic></bold></td><td align="left" rowspan="1" colspan="1">Ribosomal protein L4</td><td align="left" rowspan="1" colspan="1">cytoplasm</td><td align="left" rowspan="1" colspan="1">component of the 60S subunit and encodes a ribosomal protein</td></tr></tbody></table></alternatives></table-wrap><table-wrap id="pone.0180432.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t002</object-id><label>Table 2</label><caption><title>Forward (F) and reverse (R) primer sequences of the selected candidate reference and target genes.</title></caption><alternatives><graphic id="pone.0180432.t002g" xlink:href="pone.0180432.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">Primer sequence (5′-3′)</th><th align="left" rowspan="1" colspan="1">Amplicon size (bp)</th><th align="left" rowspan="1" colspan="1">Ta (°C)</th><th align="left" rowspan="1" colspan="1">Accession No.</th><th align="left" rowspan="1" colspan="1">Reference</th></tr></thead><tbody><tr><td align="left" rowspan="2" colspan="1"><bold><italic>18S rRNA</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>TGTGCCGCTAGAGGTGAAATT</monospace></td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1">AF173612.1</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref032" ref-type="bibr">32</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>TGGCAAATGCTTTCGCTTT</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>ALB</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>CCTGGACACCAAGGAAAT</monospace></td><td align="left" rowspan="2" colspan="1">197</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1">NM_205261.2</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref033" ref-type="bibr">33</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>TGTGGACGCCGATAGAAT</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>ACTB</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>CTGTGCCCATCTATGAAGGCTA</monospace></td><td align="left" rowspan="2" colspan="1">139</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi?db=nucleotide&id=45382926">NM_205518.1</ext-link></td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref033" ref-type="bibr">33</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>ATTTCTCTCTCGGCTGTGGTG</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>B2M</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>CGTCCTCAACTGCTTCGTG</monospace></td><td align="left" rowspan="2" colspan="1">194</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1">NM_001001750.1</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref033" ref-type="bibr">33</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>TCTCGTGCTCCACCTTGC</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>CA2</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>TCAAACCAAGGGGAAACAAGC</monospace></td><td align="left" rowspan="2" colspan="1">99</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/nucleotide/46048695?report=gbwithparts">NM_205317.1</ext-link></td><td align="left" rowspan="2" colspan="1">this study</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>GTAGTCAGGGAGCCAGGGTA</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>CST3</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>ATGAGAACGACGAGGGCTTG</monospace></td><td align="left" rowspan="2" colspan="1">130</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1">NM_205500.2</td><td align="left" rowspan="2" colspan="1">this study</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>ATTCCAGACACGAGCTGCC</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>GAPDH</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>GAGGGTAGTGAAGGCTGCTG</monospace></td><td align="left" rowspan="2" colspan="1">113</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1">NM_204305.1</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref034" ref-type="bibr">34</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>CATCAAAGGTGGAGGAATGG</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>HMBS</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>GGCTGGGAGAATCGCATAGG</monospace></td><td align="left" rowspan="2" colspan="1">131</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1">XM_417846.2</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref035" ref-type="bibr">35</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>TCCTGCAGGGCAGATACCAT</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>HPRT1</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>ACTGGCTGCTTCTTGTG</monospace></td><td align="left" rowspan="2" colspan="1">245</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1">NM_204848.1</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref033" ref-type="bibr">33</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>GGTTGGGTTGTGCTGTT</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>RPL4</italic></bold></td><td align="left" rowspan="1" colspan="1">F: <monospace>TTATGCCATCTGTTCTGCC</monospace></td><td align="left" rowspan="2" colspan="1">235</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi?db=nucleotide&id=55926193">NM_001007479.1</ext-link></td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref032" ref-type="bibr">32</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>GCGATTCCTCATCTTACCCT</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>SLC25A38</italic></bold><xref ref-type="table-fn" rid="t002fn001">*</xref></td><td align="left" rowspan="1" colspan="1">F: <monospace>AGACACGGTATGAGAGTGGA</monospace></td><td align="left" rowspan="2" colspan="1">139</td><td align="left" rowspan="2" colspan="1">63</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/nucleotide/513169473?report=gbwithparts">XM_418818.3</ext-link></td><td align="left" rowspan="2" colspan="1">this study</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>ATCCCAGAGAAAGGTGCGTC</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>CPOX</italic></bold><xref ref-type="table-fn" rid="t002fn001">*</xref></td><td align="left" rowspan="1" colspan="1">F: <monospace>GAGAGGACGGTATGTGGAGT</monospace></td><td align="left" rowspan="2" colspan="1">187</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/nucleotide/513161781?report=gbwithparts">XM_004938236.1</ext-link></td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref036" ref-type="bibr">36</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>TTTGGGATTGCGGAGAAC</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>ALAS1</italic></bold><xref ref-type="table-fn" rid="t002fn001">*</xref></td><td align="left" rowspan="1" colspan="1">F: <monospace>GGTGGACAGGAAAGGTAAAGA</monospace></td><td align="left" rowspan="2" colspan="1">197</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1"><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/nucleotide/65301158?report=gbwithparts">NM_001018012.1</ext-link></td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref036" ref-type="bibr">36</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>ACTGGTCATACTGGAAGGTG</monospace></td></tr><tr><td align="left" rowspan="2" colspan="1"><bold><italic>FECH</italic></bold><xref ref-type="table-fn" rid="t002fn001">*</xref></td><td align="left" rowspan="1" colspan="1">F: <monospace>TGCTTTGCCGATCACAT</monospace></td><td align="left" rowspan="2" colspan="1">112</td><td align="left" rowspan="2" colspan="1">60</td><td align="left" rowspan="2" colspan="1">U68033.1</td><td align="left" rowspan="2" colspan="1">[<xref rid="pone.0180432.ref036" ref-type="bibr">36</xref>]</td></tr><tr><td align="left" rowspan="1" colspan="1">R: <monospace>CACGGTTCACCACAGACAT</monospace></td></tr></tbody></table></alternatives><table-wrap-foot><fn id="t002fn001"><p>*Genes used as candidate target genes for validation of reference genes at three different time-points of eggshell formation</p></fn></table-wrap-foot></table-wrap></sec><sec id="sec004"><title>Laying hens and tissue sampling</title><sec id="sec005"><title>Effect of time-points on stability of reference gene expression (Experiment 1)</title><p>Based on the intensity of brown eggshell colour and uniformity in egg weight, 20 hens out of a flock of 63 Hy-Line Brown laying hens were selected. The laying production of the selected hens was 100%. The feed offered was premium top layer mash (Barastock, Australia). At the time of the experiment, hens were 36–37 weeks old. From the selected 20 hens, eggshell colour (L*) and egg weight (g) were measured using a Spectrophotometer (Konica CM-2600d Ramsey, NJ, USA) [<xref rid="pone.0180432.ref037" ref-type="bibr">37</xref>] and analytical weighing balance Quintix513-1S (Sartorius Lab Instruments GmbH & Co. KG Goettingen, Germany), respectively. The hens were divided into four groups in such a way that the average L* values and egg weight were not significantly different among the selected groups (<xref ref-type="table" rid="pone.0180432.t003">Table 3</xref>). Individual hen oviposition time was monitored using a video camera at the time of sampling. Four groups of hens were sampled based on time-points (post-oviposition time 2, 5, 15 and 23.5 hrs). The hens were euthanized by CO<sub>2</sub> and the shell gland tissue was excised within 2 minutes of the euthanization. Approximately 500 mg tissue was taken from the centre of the shell gland after opening the shell gland from the anterior-ventral position and transferred directly to RNALater (Sigma Aldrich, Australia). The samples were stored at -20°C until further processing for total RNA extraction.</p><table-wrap id="pone.0180432.t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t003</object-id><label>Table 3</label><caption><title>Eggshell colour (L*) and egg weight of the hens selected for reference gene study optimisation at different time-points.</title><p>The eggs were collected and analysed before dividing the experimental hens into various groups. On the basis of eggshell variables, hens were divided into groups in such a way that the variables were not significantly different among groups.</p></caption><alternatives><graphic id="pone.0180432.t003g" xlink:href="pone.0180432.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="center" rowspan="2" colspan="1">Variable</th><th align="center" colspan="4" rowspan="1">Time-point (hr)</th><th align="center" rowspan="2" colspan="1">P value</th></tr><tr><th align="center" rowspan="1" colspan="1">2 (5 hens)</th><th align="center" rowspan="1" colspan="1">5 (5 hens)</th><th align="center" rowspan="1" colspan="1">15 (5 hens)</th><th align="center" rowspan="1" colspan="1">23.5 (5 hens)</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">L* value</td><td align="left" rowspan="1" colspan="1">56.24±0.99</td><td align="left" rowspan="1" colspan="1">54.38±0.98</td><td align="left" rowspan="1" colspan="1">54.19±0.89</td><td align="left" rowspan="1" colspan="1">54.36±0.28</td><td align="char" char="." rowspan="1" colspan="1">0.2982</td></tr><tr><td align="left" rowspan="1" colspan="1">Egg weight (g)</td><td align="left" rowspan="1" colspan="1">62.81±3.51</td><td align="left" rowspan="1" colspan="1">63.95±2.31</td><td align="left" rowspan="1" colspan="1">62.96±1.66</td><td align="left" rowspan="1" colspan="1">65.65±1.77</td><td align="char" char="." rowspan="1" colspan="1">0.8326</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="t003fn001"><p>Values are mean ±S.E.</p></fn></table-wrap-foot></table-wrap></sec><sec id="sec006"><title>Effect of time-points and nicarbazin on stability of reference gene expression (Experiment 2)</title><p>A total of 30 hens having 100% laying efficiency were selected based on the intensity of eggshell colour (L*) and egg weight from the remaining flock of 43 Hy-Line Brown laying hens (<xref ref-type="table" rid="pone.0180432.t004">Table 4</xref>). Rearing conditions were the same as described previously. At the time of the experiment, hens were 42–45 weeks old and were divided into groups with a 2×3 factorial design (<xref ref-type="table" rid="pone.0180432.t004">Table 4</xref>). The hens were divided into groups in such a way that the average L* values and egg weight (g) were not significantly different among groups. Out of 30 hens, 15 hens were fed nicarbazin @100mg/kg of commercial layer diet while the control hens were fed only commercial layer diet. Nicarbazin was fed to each group at a time in order to allow time for the processing of the treated hens without intoxicating them. The eggshell colour (L*) and egg weight were recorded for all the hens from prior to treatment until the tissue collection. Procedures for tissue collection and handling were the same as mentioned previously.</p><table-wrap id="pone.0180432.t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t004</object-id><label>Table 4</label><caption><title>Eggshell colour (L*) and egg weight of the hens selected for reference gene optimisation with nicarbazin treatment at different time-points.</title><p>The eggs were collected and analysed before dividing the experimental hens into various groups. On the basis of eggshell variables, hens were divided into groups in such a way that the variables were not significantly different among groups.</p></caption><alternatives><graphic id="pone.0180432.t004g" xlink:href="pone.0180432.t004"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="4" colspan="1">Variable</th><th align="center" colspan="8" rowspan="1">Group</th></tr><tr><th align="left" colspan="4" rowspan="1">Control (15 hens)</th><th align="left" colspan="4" rowspan="1">Nicarbazin (15 hens)</th></tr><tr><th align="left" colspan="3" rowspan="1">Time-point (hr)</th><th align="left" rowspan="2" colspan="1">P value</th><th align="left" colspan="3" rowspan="1">Time-point (hr)</th><th align="left" rowspan="2" colspan="1">P value</th></tr><tr><th align="left" rowspan="1" colspan="1">5 (5hens)</th><th align="left" rowspan="1" colspan="1">15 (5hens)</th><th align="left" rowspan="1" colspan="1">23.5 (5hens)</th><th align="left" rowspan="1" colspan="1">5 (5hens)</th><th align="left" rowspan="1" colspan="1">15 (5hens)</th><th align="left" rowspan="1" colspan="1">23.5 (5hens)</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><bold>L* value</bold></td><td align="left" rowspan="1" colspan="1">58.69±0.59</td><td align="left" rowspan="1" colspan="1">59.54±0.89</td><td align="left" rowspan="1" colspan="1">59.39±0.63</td><td align="char" char="." rowspan="1" colspan="1">0.6827</td><td align="left" rowspan="1" colspan="1">58.46±0.67</td><td align="left" rowspan="1" colspan="1">58.89±0.52</td><td align="left" rowspan="1" colspan="1">59.06±0.51</td><td align="char" char="." rowspan="1" colspan="1">0.7548</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Egg weight (g)</bold></td><td align="left" rowspan="1" colspan="1">63.73±1.77</td><td align="left" rowspan="1" colspan="1">63.29±1.58</td><td align="left" rowspan="1" colspan="1">63.87±1.58</td><td align="char" char="." rowspan="1" colspan="1">0.6969</td><td align="left" rowspan="1" colspan="1">58.74±1.69</td><td align="left" rowspan="1" colspan="1">60.93±1.41</td><td align="left" rowspan="1" colspan="1">62.65±1.19</td><td align="char" char="." rowspan="1" colspan="1">0.2019</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="t004fn001"><p>Values are mean ±S.E.</p></fn></table-wrap-foot></table-wrap></sec></sec><sec id="sec007"><title>Total RNA extraction and purification</title><p>Total RNA was extracted using TRIsure (Bioline, Australia), according to the manufacturer's instructions. Briefly, an approximately 50 mg of tissue (wet weight) was homogenized in 1 mL of TRIsure using an IKA T10 basic Homogenizer (Wilmington, NC, USA). After the RNA pellet was washed with 75% ethanol and subsequently air-dried for 10–15 minutes, 50 μL of UltraPure<sup>™</sup> DEPC-Treated water (Ambion, USA) was used to dissolve RNA pellets. The total RNA was further purified using RNeasy Mini Kit (Qiagen, GmbH, Hilden, Germany) as per the manufacturer’s instructions. The elution of RNA from the spin column with 50 μL of RNase-free water was repeated twice and the eluted RNA solutions were mixed thoroughly. The purified RNA was analysed in a NANODROP-8000 spectrophotometer (ThermoFisher Scientific, Wilmington, DE, USA) to measure its quantity and purity. RNA integrity was examined using 1.5% agarose gel electrophoresis in 1× UltraPure<sup>™</sup> TAE Buffer (ThermoFisher Scientific, Australia) solution. The 0.1–2 Kb RNA Ladder (Ambion<sup>™</sup>, Australia) was used as a marker. A 1:1000 Lonza GelStar<sup>™</sup> Nucleic Acid Gel Stain in Orange G loading dye was used to stain RNA and the gel was viewed and photographed using a camera. Only RNA showing two distinct 28S and 18S bands without smear was regarded integrate and used in downstream qPCR assay.</p></sec><sec id="sec008"><title>Quantitative PCR</title><p>qPCR was performed with the SensiFAST SYBR<sup>®</sup> Lo-ROX One-Step RT-PCR Kit (Bioline, Australia). Master mix was prepared as per the manufacturer’s protocol and 4 μL of RNA template from 1:100 dilutions with the exception of <italic>18S rRNA</italic> (that was in 10<sup>−4</sup> dilutions) was added to the reaction wells using Corbett CAS1200 robotics (Corbett Life Science, Sydney, Australia). The reaction was run in triplicates of 20 μL in a Rotor- Gene Disc 100 (Qiagen, Sydney, Australia) with a Rotor-Gene 6000 thermocycler (Corbett Research, Sydney, Australia). No template control (NTC) and no reverse transcriptase (-RT) control were also included to detect possible contamination. Thermocycling conditions for a 2-step PCR were: reverse transcription at 45°C for 10 minutes, first denaturation at 95°C for 2 minutes, then 40 cycles of denaturation at 95°C for 5s and annealing at 60°C or 63°C for 20s. The fluorescent data were acquired at the end of each annealing step during PCR cycles. A melting step was conducted to assess the specificity of PCR amplification. The PCR products were examined on Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany) gel using DNA 1000 Kit as per the manufacturer’s instructions to estimate the size of the amplicons for specificity.</p><p>PCR amplification efficiencies and correlation coefficients (R<sup>2</sup>) were determined with the amplifications of a series of six 10-fold dilutions. The qPCR data of the genes were processed further when the PCR amplification efficiency was in a range of 90 to 105%, and linear correlation coefficient R<sup>2</sup> > 0.980 were considered of high standard [<xref rid="pone.0180432.ref038" ref-type="bibr">38</xref>].</p></sec><sec id="sec009"><title>Statistical analysis</title><p>The eggshell colour (L*) and egg weight data were analysed by Statview software (SAS Institute Inc., Version 5.0.1.0). A two-way analysis of variance (ANOVA) was conducted taking time-point and nicarbazin treatment as independent and L* and egg weight as dependent variables. Level of significance was indicated by probability of less than 5%. The Fishers LSD test was used to differentiate levels of significance between mean values.</p><p>To determine the expression stability of 10 different reference genes, the geNorm module in qbase+ software version 3.0 (Biogazelle, Belgium) was used to calculate the gene expression stability measure (geNorm M) [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>, <xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>]. The input data for qbase+ were generated using the relative quantities based on comparative quantification cycle (Cq). To be consistent, the Cq values for <italic>18S rRNA</italic>, which were from 10<sup>−4</sup> dilutions, were adjusted according to 1:100 dilutions. Any triplicate reaction with difference more than 0.5 cycle was excluded from the analysis. To select the most stable genes, geNorm re-calculates the M stability measures after removing the least stable genes and repeats the process until the one most stable gene remains [<xref rid="pone.0180432.ref022" ref-type="bibr">22</xref>, <xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>]. To test the minimum number of reference genes, geNorm calculates a pairwise variation (geNorm V) based on V<sub>n/n+1</sub> and a higher value indicates a significant effect of additional gene on data normalization. Normally, the benefit of using an extra (n+1)th reference gene is limited as soon as the V<sub>n/n+1</sub> value drops below the 0.15 threshold. Due to its sequential elimination of less stable genes so as to produce less bias on the output of analysis, this program was used as the primary base for the selection of the reference genes.</p><p>In addition, another two programs, i.e., NormFinder [<xref rid="pone.0180432.ref024" ref-type="bibr">24</xref>, <xref rid="pone.0180432.ref025" ref-type="bibr">25</xref>] and BestKeeper [<xref rid="pone.0180432.ref026" ref-type="bibr">26</xref>], were used to analyse the stability of gene expression as complementary measures to safeguard the output generated from geNorm. The raw Cq values were exported and analysed in NormFinder for reference gene expression stabilities. An Excel based BestKeeper (Version 1) software was also used to determine the best stable reference genes based on Pearson correlation coefficient (r), coefficient of variance (CV) and standard deviation (SD). The most stable genes were selected on the basis of the stability of the genes across the three software.</p><p>Data for the candidate target gene expression using the two most stable and the two least stable reference genes were analysed in qbase+ by scaling the average relative quantities across all unknown samples per target gene [<xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>, <xref rid="pone.0180432.ref039" ref-type="bibr">39</xref>]. Effect of time-points on the relative expression levels of the candidate target genes was analysed using one-way ANOVA. Tukey-Kramer method was used to correct the p-value results (corrected p value < 0.05) for the pairwise group comparisons in the ANOVA test [<xref rid="pone.0180432.ref023" ref-type="bibr">23</xref>].</p></sec></sec><sec sec-type="results" id="sec010"><title>Results</title><sec id="sec011"><title>Primers specificity and efficiency</title><p>All of the primer pairs were specific in amplifications by showing a single band on the Agilent 2100 Bioanalyzer gel (<xref ref-type="fig" rid="pone.0180432.g001">Fig 1</xref>). The melting curve analyses of all primer pairs are depicted in <xref ref-type="fig" rid="pone.0180432.g002">Fig 2</xref>. The amplification efficiency of all ten candidate reference genes was between 93% and 101%. The amplification efficiencies were 100% for <italic>18S rRNA</italic>, 98% for <italic>ACTB</italic>, 101% for <italic>ALB</italic>, 93% for <italic>B2M</italic>, 94% for <italic>CA2</italic>, 100% for each of <italic>CST3</italic> and <italic>GAPDH</italic>, 94% for <italic>HMBS</italic>, 98% for <italic>HPRT1</italic> and 94% for <italic>RPL4</italic>. The overall expression pattern (Cq values) for these ten reference genes is shown in <xref ref-type="fig" rid="pone.0180432.g003">Fig 3A</xref>. Most of the reference genes were highly expressed, with average Cq values between 12 and 22 cycles, except <italic>ALB</italic>, which showed average Cq values around 28 cycles (<xref ref-type="fig" rid="pone.0180432.g003">Fig 3A</xref>). The expression pattern of all ten reference genes was calculated in the combined dataset of four different time-points (2, 5, 15, 23.5 hr, post-oviposition times).</p><fig id="pone.0180432.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.g001</object-id><label>Fig 1</label><caption><title>Amplification of the gene fragments from the eggshell gland tissue of laying hens to assess the specificities of the primers used in the current study.</title><p>L, DNA ladder; 1. <italic>18S rRNA</italic> (63 bp); 2. <italic>ACTB</italic> (139 bp); 3. <italic>ALB</italic> (197 bp); 4. <italic>B2M</italic> (194 bp); 5. <italic>CA2</italic> (99 bp); 6. <italic>CST3</italic> (130 bp); 7. <italic>GAPDH</italic> (113 bp); 8. <italic>HMBS</italic> (131 bp); 9. <italic>HPRT1</italic> (245 bp); 10. <italic>RPL4</italic> (235); 11. <italic>SLC25A38</italic> (139 bp); 12. <italic>FECH</italic> (112 bp). All the amplified products were in accordance to the expected sizes. The amplified products were run on Agilent 2100 Bioanalyzer using Agilent DNA 1000 Kit as per the manufacturer’s instructions. The upper (purple) and lower (green) markers act as internal standards and are used to align the ladder analysis with the individual DNA sample analysis. The standard curve (plotting migration time against DNA amplicon size), in conjunction with the markers, is then used to calculate DNA fragment sizes for each well from the migration times measured (for more detail see Agilent 2100 Bioanalyzer Users Guide for Molecular Assays).</p></caption><graphic xlink:href="pone.0180432.g001"></graphic></fig><fig id="pone.0180432.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.g002</object-id><label>Fig 2</label><caption><title>Melting curves of the amplicons from 10 candidate reference genes showing that the amplifications were specific and no primer dimers were present.</title><p>All of the amplicons showed a single peak in melting curve analysis. After qPCR cycles, a melting phase at a ramp from 50°C to 99°C at 1°C increment was conducted to assess the specificity of PCR amplification.</p></caption><graphic xlink:href="pone.0180432.g002"></graphic></fig><fig id="pone.0180432.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.g003</object-id><label>Fig 3</label><caption><title>Effect of time-points and combination of time-points and nicarbazin treatment on the expression stability of reference genes in the shell gland of laying hens.</title><p>A). Mean Cq values of 10 reference genes affected by time-points (Experiment 1). Bars represent standard deviation. B). Pairwise variation (geNorm V) of the optimal number of reference genes affected by time-points (Experiment 1). Pairwise variation (V<sub>n/n+1</sub>) was analysed between the normalization factors NF<sub>n</sub> and NF<sub>n+1</sub> to determine the optimal number of reference genes. In the geNorm V graph, each bar represents change in normalization accuracy by stepwise inclusion of most to least stable gene. C). Average expression stability (geNorm M) values of reference genes excluding time-point 2 (post-oviposition time, 2 hours; Experiment 1). Most stable genes have lower M value. Time-point 2 was excluded to compare it with the control group of the time-points and nicarbazin treatment experiment (Experiment 2). D). Mean Cq values of 10 reference genes affected by time-points and nicarbazin treatment (control and nicarbazin groups combined). Bars represent standard deviation. E). Pairwise variation (geNorm V) of the optimal number of reference genes (control, nicarbazin treated and both the groups combined). F). Average expression stability (geNorm M) values of reference genes (time-points (Experiment 1) and control hens from the time-points and nicarbazin treated groups (Experiment 2)). Most stable genes have lower M value.</p></caption><graphic xlink:href="pone.0180432.g003"></graphic></fig></sec><sec id="sec012"><title>Effect of time-points on stability of reference gene expression (Experiment 1)</title><p>Based on the expression stability (geNorm M), <italic>HPRT1</italic> and <italic>HMBS</italic> were the two most stable genes. The average expression stabilities (geNorm M) of the ten reference genes were within the acceptable range (<0.50) that varied from 0.252 (<italic>HPRT1</italic>) to 0.482 (<italic>ACTB</italic>) (<xref ref-type="table" rid="pone.0180432.t005">Table 5</xref>). The pairwise variation (geNorm V) also chose <italic>HPRT1</italic> and <italic>HMBS</italic> as the best set of genes to be used for expression data analysis (<xref ref-type="fig" rid="pone.0180432.g003">Fig 3B</xref>). The geNorm V value of <italic>HPRT1</italic> and <italic>HMBS</italic>, which was 0.086, indicated that stepwise inclusion of the next most stable reference gene (<italic>18S rRNA</italic>) is not necessary for data normalization. In line with geNorm M results, all ten reference genes showed geNorm V <0.15 (a default cut-off value). The results of NormFinder and BestKeeper analyses were slightly different from those of the geNorm. NormFinder ranked <italic>18S rRNA</italic> and <italic>RPL4</italic> as the two most stable reference genes, while the two most stable reference genes in BestKeeper were <italic>18S rRNA</italic> and <italic>HPRT1</italic> (<xref ref-type="table" rid="pone.0180432.t005">Table 5</xref>). Nevertheless, all the genes analysed by NormFinder and BestKeeper had SD < 1.0, showing that these genes were overall stably expressed in the tissue under investigation. The two least stable reference genes across all the three statistical tools were <italic>CA2</italic> and <italic>ACTB</italic> (<xref ref-type="table" rid="pone.0180432.t005">Table 5</xref>). The expression data were also analysed in geNorm excluding post-oviposition time 2 hours but this had no significant effect on the ranking of genes (<xref ref-type="fig" rid="pone.0180432.g003">Fig 3C</xref>).</p><table-wrap id="pone.0180432.t005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t005</object-id><label>Table 5</label><caption><title>Stability values of reference genes affected by four different time-points performed in Experiment 1.</title><p>Means of relative expression levels of genes at four different time-points (post-oviposition times 2, 5, 15, 23.5 hrs) were used to calculate the expression stability of genes across the three statistical software.</p></caption><alternatives><graphic id="pone.0180432.t005g" xlink:href="pone.0180432.t005"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="2" colspan="1">Rank</th><th align="left" colspan="2" rowspan="1">geNorm</th><th align="left" colspan="2" rowspan="1">NormFinder</th><th align="left" colspan="2" rowspan="1">BestKeeper</th></tr><tr><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">M value</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.252</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.065</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.190</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.254</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.178</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.220</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.264</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.182</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.230</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.284</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.261</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.240</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.306</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.288</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.260</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.331</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.331</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.290</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.352</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.388</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.300</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.400</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.442</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.370</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.440</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.489</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.500</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.482</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.551</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.500</td></tr></tbody></table></alternatives></table-wrap></sec><sec id="sec013"><title>Effect of time-points and nicarbazin treatment on stability of reference gene expression (Experiment 2)</title><p>The overall expression pattern (Cq values) of all the ten reference genes is depicted in <xref ref-type="fig" rid="pone.0180432.g003">Fig 3D</xref>. Calculating the expression stability of the reference genes for both the groups combined (control and nicarbazin treated), geNorm ranked <italic>HMBS</italic> and <italic>HPRT1</italic> as the two most stable genes (<xref ref-type="table" rid="pone.0180432.t006">Table 6</xref>). The pairwise variations (geNorm V) for the control, nicarbazin treated and both the groups combined are shown in <xref ref-type="fig" rid="pone.0180432.g003">Fig 3E</xref>. The geNorm V values of the pairwise variation of the ten reference genes for the control, nicarbazin and both the groups combined were <0.15 (a default cut-off value). In order to gain insight into any difference between the expression data of the time-points experiment (Experiment 1) and the combination of time-points and nicarbazin treatment experiment (Experiment 2), the expression data of time-points and the control hens from time-points and nicarbazin treatment were analysed together with the exclusion of post-oviposition time 2 hours. The gene ranking was reshuffled and <italic>HMBS</italic> was ranked as the least stable gene (<xref ref-type="fig" rid="pone.0180432.g003">Fig 3F</xref>). NormFinder ranked <italic>HMBS</italic> and <italic>RPL4</italic>, while the BestKeeper showed <italic>HPRT1</italic> and <italic>HMBS</italic> as the two most stable genes (<xref ref-type="table" rid="pone.0180432.t006">Table 6</xref>). Categorising all the ten reference genes into most stable, middle order and least stable, the reshuffling in gene ranking was different but relatively consistent for the first four genes across all the three statistical tools. It seems that the genes falling in mid order were more variable in stability when analysed for comparison by the three statistical tools. For the overall ranking obtained by the three algorithms, the two most stable reference genes for the total dataset were <italic>HMBS</italic> and <italic>HPRT1</italic>, while the two least stable genes were <italic>B2M</italic> and <italic>CA2</italic> (<xref ref-type="table" rid="pone.0180432.t006">Table 6</xref>).</p><table-wrap id="pone.0180432.t006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t006</object-id><label>Table 6</label><caption><title>Overall stability values of reference genes affected by time-points and nicarbazin treatment.</title><p>Means of relative expression levels in respective groups at three different time-points (post-oviposition times 5, 15, 23.5 hrs) and with nicarbazin treatments (yes, no) were used to calculate the expression stability of genes in responses to the time points and nicarbazin treatment across the three statistical software.</p></caption><alternatives><graphic id="pone.0180432.t006g" xlink:href="pone.0180432.t006"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="2" colspan="1">Rank</th><th align="left" colspan="2" rowspan="1">geNorm</th><th align="left" colspan="2" rowspan="1">NormFinder</th><th align="left" colspan="2" rowspan="1">BestKeeper</th></tr><tr><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">M value</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.236</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.182</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.210</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.243</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.191</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.220</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.247</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.266</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.230</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.300</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.291</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.270</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.331</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.296</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.320</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.369</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.409</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.330</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.418</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.446</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.380</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.454</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.475</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.450</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.488</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.491</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.480</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.601</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.991</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.860</td></tr></tbody></table></alternatives></table-wrap><p>When the data were normalized in geNorm for each group separately, the control group showed that <italic>HMBS</italic> and <italic>HPRT1</italic> were the most stable reference genes (<xref ref-type="table" rid="pone.0180432.t007">Table 7</xref>). In the same group, both NormFinder and BestKeeper ranked <italic>RPL4</italic> and <italic>HMBS</italic> as the two most stable reference genes. However, the overall ranking of the genes falling in middle order was reshuffled following the analyses by all the three statistical software. In the control group, the two least stable genes across all the three applets were <italic>B2M</italic> and <italic>CA2</italic>. The only genes that showed higher M value than the cut-off value in geNorm were <italic>CA2</italic> and <italic>B2M</italic> (<xref ref-type="table" rid="pone.0180432.t007">Table 7</xref>).</p><table-wrap id="pone.0180432.t007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t007</object-id><label>Table 7</label><caption><title>Stability values of reference genes affected by time-points.</title><p>Means of relative expression levels of the genes in the control group at three different time-points (post-oviposition times 5, 15, 23.5 hrs) were used to calculate the expression stability of genes across the three statistical software.</p></caption><alternatives><graphic id="pone.0180432.t007g" xlink:href="pone.0180432.t007"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="2" colspan="1">Rank</th><th align="left" colspan="2" rowspan="1">geNorm</th><th align="left" colspan="2" rowspan="1">NormFinder</th><th align="left" colspan="2" rowspan="1">BestKeeper</th></tr><tr><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">M value</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.207</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.138</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.190</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.207</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.172</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.220</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.214</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.225</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.240</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.274</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.300</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.260</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.305</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.309</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.290</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.373</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.455</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.340</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.424</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.459</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.400</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.461</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.499</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.450</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.503</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.544</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.480</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.609</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.963</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.710</td></tr></tbody></table></alternatives></table-wrap><p>The reference genes ranking in the nicarbazin treatment group was slightly different from the control group. The geNorm ranked <italic>HMBS</italic> and <italic>GAPDH</italic> as the two most stable genes followed by <italic>HPRT1</italic> (<xref ref-type="table" rid="pone.0180432.t008">Table 8</xref>). NormFinder showed <italic>HMBS</italic> and <italic>ALB</italic>, while the BestKeeper showed <italic>HPRT1</italic> and <italic>HMBS</italic> as the two most stable reference genes (<xref ref-type="table" rid="pone.0180432.t008">Table 8</xref>). The two least stable genes across all the three applets were <italic>B2M</italic> and <italic>CST3</italic> (<xref ref-type="table" rid="pone.0180432.t008">Table 8</xref>). The stability values of <italic>B2M</italic> were slightly higher than the cut-off value (geNorm M < 0.5; SD < 1.0) when the data were analysed both in the geNorm and NormFinder.</p><table-wrap id="pone.0180432.t008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.t008</object-id><label>Table 8</label><caption><title>Stability values of reference genes affected by nicarbazin treatment (nicarbazin group only).</title><p>Means of relative expression levels of the genes in the nicarbazin group at three different time-points (post-oviposition times 5, 15, 23.5 hrs) were used to calculate the expression stability of genes across the three statistical software.</p></caption><alternatives><graphic id="pone.0180432.t008g" xlink:href="pone.0180432.t008"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="2" colspan="1">Rank</th><th align="left" colspan="2" rowspan="1">geNorm</th><th align="left" colspan="2" rowspan="1">NormFinder</th><th align="left" colspan="2" rowspan="1">BestKeeper</th></tr><tr><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">M value</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th><th align="left" rowspan="1" colspan="1">Gene</th><th align="left" rowspan="1" colspan="1">SD</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.262</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.199</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.190</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.273</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.226</td><td align="left" rowspan="1" colspan="1"><italic>HMBS</italic></td><td align="char" char="." rowspan="1" colspan="1">0.230</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.278</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.231</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.260</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.329</td><td align="left" rowspan="1" colspan="1"><italic>HPRT1</italic></td><td align="char" char="." rowspan="1" colspan="1">0.284</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.280</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"><italic>RPL4</italic></td><td align="char" char="." rowspan="1" colspan="1">0.340</td><td align="left" rowspan="1" colspan="1"><italic>GAPDH</italic></td><td align="char" char="." rowspan="1" colspan="1">0.289</td><td align="left" rowspan="1" colspan="1"><italic>ALB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.280</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.360</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.297</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.300</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.403</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.427</td><td align="left" rowspan="1" colspan="1"><italic>18S rRNA</italic></td><td align="char" char="." rowspan="1" colspan="1">0.340</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.435</td><td align="left" rowspan="1" colspan="1"><italic>ACTB</italic></td><td align="char" char="." rowspan="1" colspan="1">0.435</td><td align="left" rowspan="1" colspan="1"><italic>CA2</italic></td><td align="char" char="." rowspan="1" colspan="1">0.460</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.463</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.513</td><td align="left" rowspan="1" colspan="1"><italic>CST3</italic></td><td align="char" char="." rowspan="1" colspan="1">0.490</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.583</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">1.008</td><td align="left" rowspan="1" colspan="1"><italic>B2M</italic></td><td align="char" char="." rowspan="1" colspan="1">0.860</td></tr></tbody></table></alternatives></table-wrap><p>In order to determine the consistency of the stabilities of the reference genes analysed by geNorm, NormFinder and BestKeeper, the relative expressions of the ten genes were compared as shown in <xref ref-type="fig" rid="pone.0180432.g004">Fig 4</xref>. The results showed that the expression stability of the genes was consistent from the analyses performed by the three programs except <italic>B2M</italic> when the birds were treated with nicarbazin at three time-points of egg shell formation (birds age 42–45 weeks).</p><fig id="pone.0180432.g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.g004</object-id><label>Fig 4</label><caption><title>Comparison of the relative expression stability of 10 candidate reference genes analysed by geNorm, NormFinder and BestKeeper.</title><p>The geNorm M and SD values of 10 reference genes were graphed to assess the ranking pattern of the genes assigned by three different software. A). Four different time-points (2, 5, 15, 23.5 hr) post-oviposition in the birds at age of 36–37 weeks. B). Groups treated with or without nicarbazin at three different time-points (5, 15, 23.5 hr) post-oviposition of birds at 42–45 weeks. C). Three different time-points (5, 15, 23.5 hr) post-oviposition times of the birds without nicarbazin treatment at age of 42–45 weeks. D). Three different time-points (5, 15, 23.5 hr) post-oviposition of the birds with nicarbazin treatment at age of 42–45 weeks.</p></caption><graphic xlink:href="pone.0180432.g004"></graphic></fig></sec><sec id="sec014"><title>Expression of candidate target genes using most stable and least stable reference genes</title><p>The level of significance (p value) changed for all four candidate target genes when the relative expression data were normalized with the two most stable (<italic>HMBS</italic>, <italic>HPRT1</italic>) and the two least stable reference genes (<italic>B2M</italic>, <italic>CA2</italic>) (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5</xref>). The p value increased when the data were normalized with the two least stable reference genes. The relative expression level of <italic>SLC25A38</italic> was significantly different (p value 4.8E-10) among different time-points when the data were normalized with the two most stable reference genes (<italic>HMBS</italic> and <italic>HPRT1</italic>) (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5A</xref>). However, for the same candidate target gene, the level of significance decreased (p value = 0.0847) among different time-points when the data were normalized with the two least stable reference genes, <italic>B2M</italic> and <italic>CA2</italic> (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5A</xref>). The expression levels of <italic>ALAS1</italic>, <italic>CPOX</italic> and <italic>FECH</italic> were changed in terms of the p values when the data were normalized with the two most stable and the two least stable reference genes. The p value of <italic>ALAS1</italic> changed from 3.2E-12 when two most stable reference genes were used, to 0.0111 when the two least stable reference genes were used (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5B</xref>). The p value of <italic>CPOX</italic> changed from 0.0700 to 0.6134 when the data were normalized with the two least stable reference genes (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5C</xref>). The p value of <italic>FECH</italic> changed from 4.6E-16 to 9.9E-05 when the data were normalized with the two least stable reference genes (<xref ref-type="fig" rid="pone.0180432.g005">Fig 5D</xref>).</p><fig id="pone.0180432.g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0180432.g005</object-id><label>Fig 5</label><caption><title>Validation of the expression level of four candidate target genes affected by three different time-points of eggshell formation.</title><p>The candidate target genes were normalized with the two most stable (<italic>HMBS</italic> and <italic>HPRT1</italic>) and the two least stable reference genes (<italic>B2M</italic> and <italic>CA2</italic>). A). <italic>SLC25A38</italic> normalized with the two most stable genes (P = 4.8E-10); <italic>SLC25A38</italic> normalized with the two least stable genes (P = 0.0847). B) <italic>ALAS1</italic> normalized with the two most stable genes (P = 3.2E-12); <italic>ALAS1</italic> normalized with the two least stable genes (P = 0.0111). C). <italic>CPOX</italic> normalized with the two most stable genes (P = 0.0700); <italic>CPOX</italic> normalized with the two least stable genes (P = 0.6134). D). <italic>FECH</italic> normalized with the two most stable genes (P = 4.6E-16); <italic>FECH</italic> normalized with the two least stable genes (P = 9.9E-05). For the same gene in the same treatment, <sup>a,b</sup> across the bars indicate significant differences. B). For the four candidate target genes, normalized relative quantities were calculated in qbase+ based on (2<sup>^-ΔΔCq</sup>) [<xref rid="pone.0180432.ref040" ref-type="bibr">40</xref>] using gene specific amplification efficiencies [<xref rid="pone.0180432.ref041" ref-type="bibr">41</xref>], to show the relative expression of Cq levels in folds to the mean Cq of all samples of the genes.</p></caption><graphic xlink:href="pone.0180432.g005"></graphic></fig></sec></sec><sec sec-type="conclusions" id="sec015"><title>Discussion</title><p>We investigated the stability of ten reference genes in the shell gland region of the oviduct in relation to different stages of eggshell formation and nicarbazin treatment. The data analysed by the three different statistical software indicated that the overall stability of the reference genes was affected by different time-points (post-oviposition time) and nicarbazin treatment, and differences in the ranking of reference genes analysed using three statistical software were observed. The current study provides information on the expression stability of these candidate reference genes and most stably expressed reference genes are suggested for the normalization of gene expression data in the chicken shell gland.</p><p>The higher stability of <italic>HPRT1</italic> and <italic>HMBS</italic> across all the three software indicated that these two genes can be used as reference genes for the normalization of expression data in the shell gland of the brown-egg laying hen. In fact, most of the ten reference genes tested in the current study were in the acceptable range as reference genes with geNorm M value <0.5 and SD <1.0. <italic>B2M</italic> was the only exception that showed slightly higher geNorm M value in the time-points and nicarbazin treatment study. Taking time-points separately, or together with the nicarbazin treatment, the pairwise variation (geNorm V) showed that the variation between the first two most stable genes was under the cut-off value (<0.15). In qbase+, geNorm V indicates level of variation in the average values of reference gene stability with the sequential inclusion of the next stable reference gene to the equation Vn/n+1 (for calculation of the normalization factor). The analysis starts with the two most stably expressed genes being compared to the pair including the third (V2/3), and the process continues until the least stable gene is added (for example, V9/10). Generally, if a geNorm V (V2/3) of 0.300 is achieved using two most stable genes and a geNorm V (V3/4) of 0.14 is achieved with three most stable reference genes, then the average of the most stable three genes would be the optimal normalization factor for further data analysis. In the current study, the first two most stable reference genes were under the cut-off value (<0.15) of geNorm V and thus adding the third most stable reference gene for expression data normalization was not necessary. However, as indicated by the geNorm V, all of the genes (V2/3 to V9/10) showed pairwise variation < 0.15 and thus all could be used for accurate data normalization. Based on the geNorm V results, this demonstrates that all the genes analysed had relatively high stability in the shell gland tissue in response to nicarbazin treatment of chickens as well as to the cyclic changes in shell gland tissue during the egg lying cycle. Based on geNorm M results, <italic>B2M</italic> showed low expression stability in response to the stages of egg formation in the nicarbazin treatment experiment and therefore should be ruled out from being used in the normalisation of expression data while different stages of egg formation are involved in the study.</p><p>Furthermore, the difference in stability values of <italic>B2M</italic> during the egg formation stages may also be dependent on the age of the birds. It appears that it is more stable in response to the time-points when birds are younger (36–37 week vs 42–45 week of age). Nevertheless, consensus from the analyses performed by these three programs was that <italic>HMBS</italic> and <italic>HPRT1</italic> were the two most stable housekeeping genes and thus were chosen as reference genes in the current study and recommended for similar studies in the shell gland of laying hens.</p><p>To validate whether most stable reference genes identified in the study would result in more accurate assessment of target gene expression, the data obtained in time-points and nicarbazin treatment experiment for four candidate target genes were normalized with the two most stable (<italic>HPRT1</italic>, <italic>HMBS</italic>) and the two least stable reference genes (<italic>B2M</italic>, <italic>CA2</italic>). Results showed that normalizing candidate target gene expression data with the two least reference stable genes is not as accurate when compared with the data normalized with the two most stable reference genes. Therefore, the most stable genes can produce more robust and accurate results for gene expression data as recommended by the optimisation outcomes achieved in the current study.</p><p>To the best of our knowledge, no validated reference genes have been used for the normalization of expression data in the shell gland of avian species under various treatments. Thus, this is the first study to establish a set of stably expressed reference genes in the shell gland and can be used in chickens and possibly in other avian species. In different species, different reference genes under different treatments have been validated in other tissues of the reproductive system. For example, in geese, <italic>HPRT1</italic> and <italic>GAPDH</italic> were ranked as the two most stable reference genes in the ovary [<xref rid="pone.0180432.ref042" ref-type="bibr">42</xref>]. In bovine, the two most stable reference genes in the uterus were <italic>YWHAZ</italic> and <italic>GAPDH</italic> in relation to developmental stages of an embryo [<xref rid="pone.0180432.ref043" ref-type="bibr">43</xref>]. Similarly, under various toxicological treatments, in the ovary of mouse, the two most stable reference genes were <italic>RPL13a</italic> and <italic>GAPDH</italic> [<xref rid="pone.0180432.ref044" ref-type="bibr">44</xref>]. Based on limited studies already being performed, it appears <italic>HPRT1</italic> may be more stable in the reproductive system of avian species while <italic>GAPDH</italic> is more stable in mammalian species [<xref rid="pone.0180432.ref040" ref-type="bibr">40</xref>]. However, further studies are required to accumulate the information required to reach a more generalised conclusion.</p><p>The present study has demonstrated that the rankings of the expression stability of the 10 candidate reference genes had similar trends but discrepancies were observed among three different statistical programs, geNorm, NormFinder and BestKeeper. Similar discrepancies have been observed elsewhere with different species and treatments [<xref rid="pone.0180432.ref042" ref-type="bibr">42</xref>, <xref rid="pone.0180432.ref045" ref-type="bibr">45</xref>, <xref rid="pone.0180432.ref046" ref-type="bibr">46</xref>]. So far, there is no consensus as to which software is more powerful in the ranking of expression stability of candidate reference genes and researchers have given the same weight to all three programs. We have shown in this study a comparison in the consistency of the ranking of candidate reference genes among all the three software used in both experiments, indicating that all of them gave similar results and can be used for the analysis of expression data. As has been stated previously, the stability of all the chosen reference genes was in the acceptable range for reference genes. Therefore, the stability levels of these genes are essentially very close. With different algorithms in different programs, slight change of their stability orders can be expected by the analyses using these programs.</p><p>It is worth noting; however, that these three programs do not have an option for analysing the reference gene expression data generated from a factorial design, but can only perform analysis based on individual group as independent treatment. To the best of our knowledge, the optimisation of reference genes has not been performed in such a factorial design so far. Therefore, it is questionable whether the programs possess the capacity to generate a reliable ranking for an experiment designed in a factorial fashion. The gene expression stability analysis of reference genes has been reported in experiments exploring the roles of multiple factors; for example, geographical locations and ventilation in new born lambs [<xref rid="pone.0180432.ref047" ref-type="bibr">47</xref>], multiple stress conditions and different developmental tissues in pear millet [<xref rid="pone.0180432.ref048" ref-type="bibr">48</xref>] and different developmental stages and hormonal stimuli on leaves of tea [<xref rid="pone.0180432.ref049" ref-type="bibr">49</xref>]. However, multiple factors have not been considered as independent effects in the analysis of stability of reference gene expression or for their interactions. We suggest these programs should add such a capacity or new programs should be available for the analyses of data produced from factorial design experiments. This would permit the role of treatments in the expression stability of the reference genes to be more robustly investigated and their interactions explored. Such investigations are warranted to provide a more powerful statistical analysis protocol.</p></sec><sec sec-type="conclusions" id="sec016"><title>Conclusion</title><p>In summary, we have performed optimisation of reference genes in the samples collected at different time-points of egg/eggshell formation and with nicarbazin treatment in laying hens. All of the reference genes except <italic>B2M</italic> were stably expressed according to the cut-off values of the programs, and two most stably expressed genes, <italic>HMBS</italic> and <italic>HPRT1</italic> are recommended for the normalization of gene expression data in the shell gland of chickens while different shell formation stages and nicarbazin treatment are involved in the experiment. It is anticipated that these reference genes may be used for the study of reproductive system in other avian species.</p></sec></body><back><ack><p>This study was funded by Australian Egg Corporation Limited, Australia under grant number AECL 1UN121. We thank Kapil Chousalkar at The University of Adelaide, SA, Australia for his technical support in the design of this experiment.</p></ack><ref-list><title>References</title><ref id="pone.0180432.ref001"><label>1</label><mixed-citation publication-type="journal"><name><surname>Nys</surname><given-names>Y</given-names></name>. <article-title>Relationships between age, shell quality and individual rate and duration of shell formation in domestic hens</article-title>. <source>Br Poult Sci</source>. <year>1986</year>;<volume>27</volume>: <fpage>253</fpage>–<lpage>259</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref002"><label>2</label><mixed-citation publication-type="journal"><name><surname>Melek</surname><given-names>O</given-names></name>, <name><surname>Morris</surname><given-names>T</given-names></name>, <name><surname>Jennings</surname><given-names>R</given-names></name>. <article-title>The time factor in egg formation for hens exposed to ahemeral light-dark cycles</article-title>. <source>Br Poult Sci</source>. <year>1973</year>;<volume>14</volume>: <fpage>493</fpage>–<lpage>498</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1080/00071667308416056">10.1080/00071667308416056</ext-link></comment><pub-id pub-id-type="pmid">4748870</pub-id></mixed-citation></ref><ref id="pone.0180432.ref003"><label>3</label><mixed-citation publication-type="book"><name><surname>Johnson</surname><given-names>AL</given-names></name>. <source>Reproduction in the female</source>. <edition>6th ed</edition><name><surname>CC</surname><given-names>S</given-names></name>, editor. <publisher-loc>USA</publisher-loc>: <publisher-name>Academic Press</publisher-name>; <year>2015</year>.</mixed-citation></ref><ref id="pone.0180432.ref004"><label>4</label><mixed-citation publication-type="journal"><name><surname>Arias</surname><given-names>JL</given-names></name>, <name><surname>Fink</surname><given-names>DJ</given-names></name>, <name><surname>Xiao</surname><given-names>S-Q</given-names></name>, <name><surname>Heuer</surname><given-names>AH</given-names></name>, <name><surname>Caplan</surname><given-names>AI</given-names></name>. <article-title>Biomineralization and eggshells: cell-mediated acellular compartments of mineralized extracellular matrix</article-title>. <source>Int Rev Cytol</source>. <year>1993</year>;<volume>145</volume>: <fpage>217</fpage>–<lpage>250</lpage>. <pub-id pub-id-type="pmid">8500981</pub-id></mixed-citation></ref><ref id="pone.0180432.ref005"><label>5</label><mixed-citation publication-type="journal"><name><surname>Wellman-Labadie</surname><given-names>O</given-names></name>, <name><surname>Picman</surname><given-names>J</given-names></name>, <name><surname>Hincke</surname><given-names>M</given-names></name>. <article-title>Avian antimicrobial proteins: structure, distribution and activity</article-title>. <source>Worlds Poult Sci J</source>. <year>2007</year>;<volume>63</volume>: <fpage>421</fpage>–<lpage>438</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref006"><label>6</label><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>Z</given-names></name>, <name><surname>Zheng</surname><given-names>Q</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name>, <name><surname>Lu</surname><given-names>L</given-names></name>. <article-title>Microarray analysis of genes involved with shell strength in layer shell gland at the early stage of active calcification</article-title>. <source>Asian Austral J Anim</source>. <year>2013</year>;<volume>26</volume>: <fpage>609</fpage>–<lpage>624</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref007"><label>7</label><mixed-citation publication-type="journal"><name><surname>Jonchère</surname><given-names>V</given-names></name>, <name><surname>Réhault-Godbert</surname><given-names>S</given-names></name>, <name><surname>Hennequet-Antier</surname><given-names>C</given-names></name>, <name><surname>Cabau</surname><given-names>C</given-names></name>, <name><surname>Sibut</surname><given-names>V</given-names></name>, <name><surname>Cogburn</surname><given-names>LA</given-names></name>, <etal>et al</etal><article-title>Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg</article-title>. <source>BMC Genomics</source>. <year>2010</year>;<volume>11</volume>: <fpage>57</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1471-2164-11-57">10.1186/1471-2164-11-57</ext-link></comment><pub-id pub-id-type="pmid">20092629</pub-id></mixed-citation></ref><ref id="pone.0180432.ref008"><label>8</label><mixed-citation publication-type="journal"><name><surname>Dunn</surname><given-names>IC</given-names></name>, <name><surname>Wilson</surname><given-names>PW</given-names></name>, <name><surname>Lu</surname><given-names>Z</given-names></name>, <name><surname>Bain</surname><given-names>MM</given-names></name>, <name><surname>Crossan</surname><given-names>CL</given-names></name>, <name><surname>Talbot</surname><given-names>RT</given-names></name>, <etal>et al</etal><article-title>New hypotheses on the function of the avian shell gland derived from microarray analysis comparing tissue from juvenile and sexually mature hens</article-title>. <source>Gen Comp Endocrinol</source>. <year>2009</year>;<volume>163</volume>: <fpage>225</fpage>–<lpage>232</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.ygcen.2009.03.006">10.1016/j.ygcen.2009.03.006</ext-link></comment><pub-id pub-id-type="pmid">19303879</pub-id></mixed-citation></ref><ref id="pone.0180432.ref009"><label>9</label><mixed-citation publication-type="journal"><name><surname>Song</surname><given-names>G</given-names></name>, <name><surname>Seo</surname><given-names>HW</given-names></name>, <name><surname>Choi</surname><given-names>JW</given-names></name>, <name><surname>Rengaraj</surname><given-names>D</given-names></name>, <name><surname>Kim</surname><given-names>TM</given-names></name>, <name><surname>Lee</surname><given-names>BR</given-names></name>, <etal>et al</etal><article-title>Discovery of candidate genes and pathways regulating oviduct development in chickens</article-title>. <source>Biol Reprod</source>. <year>2011</year>;<volume>85</volume>: <fpage>306</fpage>–<lpage>314</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1095/biolreprod.110.089227">10.1095/biolreprod.110.089227</ext-link></comment><pub-id pub-id-type="pmid">21543768</pub-id></mixed-citation></ref><ref id="pone.0180432.ref010"><label>10</label><mixed-citation publication-type="journal"><name><surname>Additives</surname><given-names>EPo</given-names></name>, <collab>Products or Substances used in Animal F</collab>. <article-title>Scientific opinion on the safety and efficacy of Koffogran (nicarbazin) as a feed additive for chickens for fattening</article-title>. <source>EFSA Journal</source>. <year>2010</year>;<volume>8</volume>: <fpage>1551</fpage>.</mixed-citation></ref><ref id="pone.0180432.ref011"><label>11</label><mixed-citation publication-type="journal"><name><surname>Hughes</surname><given-names>B</given-names></name>, <name><surname>Jones</surname><given-names>J</given-names></name>, <name><surname>Toler</surname><given-names>J</given-names></name>, <name><surname>Solis</surname><given-names>J</given-names></name>, <name><surname>Castaldo</surname><given-names>D</given-names></name>. <article-title>Effects of exposing broiler breeders to nicarbazin contaminated feed</article-title>. <source>Poult Sci</source>. <year>1991</year>;<volume>70</volume>: <fpage>476</fpage>–<lpage>482</lpage>. <pub-id pub-id-type="pmid">2047341</pub-id></mixed-citation></ref><ref id="pone.0180432.ref012"><label>12</label><mixed-citation publication-type="journal"><name><surname>Dorne</surname><given-names>J</given-names></name>, <name><surname>Fernández-Cruz</surname><given-names>M</given-names></name>, <name><surname>Bertelsen</surname><given-names>U</given-names></name>, <name><surname>Renshaw</surname><given-names>D</given-names></name>, <name><surname>Peltonen</surname><given-names>K</given-names></name>, <name><surname>Anadon</surname><given-names>A</given-names></name>, <etal>et al</etal><article-title>Risk assessment of coccidostatics during feed cross-contamination: animal and human health aspects</article-title>. <source>Toxicol Appl Pharmacol</source>. <year>2013</year>;<volume>270</volume>: <fpage>196</fpage>–<lpage>208</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.taap.2010.12.014">10.1016/j.taap.2010.12.014</ext-link></comment><pub-id pub-id-type="pmid">21215766</pub-id></mixed-citation></ref><ref id="pone.0180432.ref013"><label>13</label><mixed-citation publication-type="journal"><name><surname>Bilandžić</surname><given-names>N</given-names></name>, <name><surname>Dolenc</surname><given-names>J</given-names></name>, <name><surname>Gačnik</surname><given-names>KŠ</given-names></name>, <name><surname>Varenina</surname><given-names>I</given-names></name>, <name><surname>Kolanović</surname><given-names>BS</given-names></name>. <article-title>Feed additives diclazuril and nicarbazin in egg and liver samples from Croatian farms</article-title>. <source>Food Addit Contam B</source>. <year>2013</year>;<volume>6</volume>: <fpage>90</fpage>–<lpage>97</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref014"><label>14</label><mixed-citation publication-type="journal"><name><surname>Leeson</surname><given-names>S</given-names></name>, <name><surname>Caston</surname><given-names>L</given-names></name>, <name><surname>Summers</surname><given-names>J</given-names></name>. <article-title>The effect of graded levels of nicarbazin on reproductive performance of laying hens</article-title>. <source>Can J Anim Sci</source>. <year>1989</year>;<volume>69</volume>: <fpage>757</fpage>–<lpage>764</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref015"><label>15</label><mixed-citation publication-type="journal"><name><surname>Thellin</surname><given-names>O</given-names></name>, <name><surname>Zorzi</surname><given-names>W</given-names></name>, <name><surname>Lakaye</surname><given-names>B</given-names></name>, <name><surname>De Borman</surname><given-names>B</given-names></name>, <name><surname>Coumans</surname><given-names>B</given-names></name>, <name><surname>Hennen</surname><given-names>G</given-names></name>, <etal>et al</etal><article-title>Housekeeping genes as internal standards: use and limits</article-title>. <source>J Biotechnol</source>. <year>1999</year>;<volume>75</volume>: <fpage>291</fpage>–<lpage>295</lpage>. <pub-id pub-id-type="pmid">10617337</pub-id></mixed-citation></ref><ref id="pone.0180432.ref016"><label>16</label><mixed-citation publication-type="journal"><name><surname>Lee</surname><given-names>PD</given-names></name>, <name><surname>Sladek</surname><given-names>R</given-names></name>, <name><surname>Greenwood</surname><given-names>CM</given-names></name>, <name><surname>Hudson</surname><given-names>TJ</given-names></name>. <article-title>Control genes and variability: absence of ubiquitous reference transcripts in diverse mammalian expression studies</article-title>. <source>Genome Res</source>. <year>2002</year>;<volume>12</volume>: <fpage>292</fpage>–<lpage>297</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1101/gr.217802">10.1101/gr.217802</ext-link></comment><pub-id pub-id-type="pmid">11827948</pub-id></mixed-citation></ref><ref id="pone.0180432.ref017"><label>17</label><mixed-citation publication-type="journal"><name><surname>Greer</surname><given-names>S</given-names></name>, <name><surname>Honeywell</surname><given-names>R</given-names></name>, <name><surname>Geletu</surname><given-names>M</given-names></name>, <name><surname>Arulanandam</surname><given-names>R</given-names></name>, <name><surname>Raptis</surname><given-names>L</given-names></name>. <article-title>Housekeeping genes; expression levels may change with density of cultured cells</article-title>. <source>J Immunol Methods</source>. <year>2010</year>;<volume>355</volume>: <fpage>76</fpage>–<lpage>79</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.jim.2010.02.006">10.1016/j.jim.2010.02.006</ext-link></comment><pub-id pub-id-type="pmid">20171969</pub-id></mixed-citation></ref><ref id="pone.0180432.ref018"><label>18</label><mixed-citation publication-type="journal"><name><surname>Uddin</surname><given-names>MJ</given-names></name>, <name><surname>Cinar</surname><given-names>MU</given-names></name>, <name><surname>Tesfaye</surname><given-names>D</given-names></name>, <name><surname>Looft</surname><given-names>C</given-names></name>, <name><surname>Tholen</surname><given-names>E</given-names></name>, <name><surname>Schellander</surname><given-names>K</given-names></name>. <article-title>Age-related changes in relative expression stability of commonly used housekeeping genes in selected porcine tissues</article-title>. <source>BMC Res Notes</source>. <year>2011</year>;<volume>4</volume>: <fpage>441</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1756-0500-4-441">10.1186/1756-0500-4-441</ext-link></comment><pub-id pub-id-type="pmid">22023805</pub-id></mixed-citation></ref><ref id="pone.0180432.ref019"><label>19</label><mixed-citation publication-type="journal"><name><surname>Cinar</surname><given-names>MU</given-names></name>, <name><surname>Islam</surname><given-names>MA</given-names></name>, <name><surname>Uddin</surname><given-names>MJ</given-names></name>, <name><surname>Tholen</surname><given-names>E</given-names></name>, <name><surname>Tesfaye</surname><given-names>D</given-names></name>, <name><surname>Looft</surname><given-names>C</given-names></name>, <etal>et al</etal><article-title>Evaluation of suitable reference genes for gene expression studies in porcine alveolar macrophages in response to LPS and LTA</article-title>. <source>BMC Res Notes</source>. <year>2012</year>;<volume>5</volume>: <fpage>107</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1756-0500-5-107">10.1186/1756-0500-5-107</ext-link></comment><pub-id pub-id-type="pmid">22340302</pub-id></mixed-citation></ref><ref id="pone.0180432.ref020"><label>20</label><mixed-citation publication-type="journal"><name><surname>Beekman</surname><given-names>L</given-names></name>, <name><surname>Tohver</surname><given-names>T</given-names></name>, <name><surname>Dardari</surname><given-names>R</given-names></name>, <name><surname>Léguillette</surname><given-names>R</given-names></name>. <article-title>Evaluation of suitable reference genes for gene expression studies in bronchoalveolar lavage cells from horses with inflammatory airway disease</article-title>. <source>BMC Mol Biol</source>. <year>2011</year>;<volume>12</volume>: <fpage>5</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1471-2199-12-5">10.1186/1471-2199-12-5</ext-link></comment><pub-id pub-id-type="pmid">21272375</pub-id></mixed-citation></ref><ref id="pone.0180432.ref021"><label>21</label><mixed-citation publication-type="journal"><name><surname>Petriccione</surname><given-names>M</given-names></name>, <name><surname>Mastrobuoni</surname><given-names>F</given-names></name>, <name><surname>Zampella</surname><given-names>L</given-names></name>, <name><surname>Scortichini</surname><given-names>M</given-names></name>. <article-title>Reference gene selection for normalization of RT-qPCR gene expression data from <italic>Actinidia delicios</italic>a leaves infected with <italic>Pseudomonas syringae pv</italic>. <italic>actinidiae</italic></article-title>. <source>Sci Rep</source>. <year>2015</year>;<volume>5</volume>: <fpage>16961</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/srep16961">10.1038/srep16961</ext-link></comment><pub-id pub-id-type="pmid">26581656</pub-id></mixed-citation></ref><ref id="pone.0180432.ref022"><label>22</label><mixed-citation publication-type="journal"><name><surname>Vandesompele</surname><given-names>J</given-names></name>, <name><surname>De Preter</surname><given-names>K</given-names></name>, <name><surname>Pattyn</surname><given-names>F</given-names></name>, <name><surname>Poppe</surname><given-names>B</given-names></name>, <name><surname>Van Roy</surname><given-names>N</given-names></name>, <name><surname>De Paepe</surname><given-names>A</given-names></name>, <etal>et al</etal><article-title>Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes</article-title>. <source>Genome Biol</source>. <year>2002</year>;<volume>3</volume>: <fpage>research0034. 0031</fpage>.<pub-id pub-id-type="pmid">12184808</pub-id></mixed-citation></ref><ref id="pone.0180432.ref023"><label>23</label><mixed-citation publication-type="other">qbase+ version 3.0. <ext-link ext-link-type="uri" xlink:href="https://www.qbaseplus.com/">https://www.qbaseplus.com/</ext-link>.</mixed-citation></ref><ref id="pone.0180432.ref024"><label>24</label><mixed-citation publication-type="other">GenEx version 6.0.1. <ext-link ext-link-type="uri" xlink:href="http://genex.gene-quantification.info/">http://genex.gene-quantification.info/</ext-link>.</mixed-citation></ref><ref id="pone.0180432.ref025"><label>25</label><mixed-citation publication-type="journal"><name><surname>Andersen</surname><given-names>CL</given-names></name>, <name><surname>Jensen</surname><given-names>JL</given-names></name>, <name><surname>Ørntoft</surname><given-names>TF</given-names></name>. <article-title>Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets</article-title>. <source>Cancer Res</source>. <year>2004</year>;<volume>64</volume>: <fpage>5245</fpage>–<lpage>5250</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1158/0008-5472.CAN-04-0496">10.1158/0008-5472.CAN-04-0496</ext-link></comment><pub-id pub-id-type="pmid">15289330</pub-id></mixed-citation></ref><ref id="pone.0180432.ref026"><label>26</label><mixed-citation publication-type="journal"><name><surname>Pfaffl</surname><given-names>MW</given-names></name>, <name><surname>Tichopad</surname><given-names>A</given-names></name>, <name><surname>Prgomet</surname><given-names>C</given-names></name>, <name><surname>Neuvians</surname><given-names>TP</given-names></name>. <article-title>Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations</article-title>. <source>Biotechnol Lett</source>. <year>2004</year>;<volume>26</volume>: <fpage>509</fpage>–<lpage>515</lpage>. <pub-id pub-id-type="pmid">15127793</pub-id></mixed-citation></ref><ref id="pone.0180432.ref027"><label>27</label><mixed-citation publication-type="journal"><name><surname>Guernsey</surname><given-names>DL</given-names></name>, <name><surname>Jiang</surname><given-names>H</given-names></name>, <name><surname>Campagna</surname><given-names>DR</given-names></name>, <name><surname>Evans</surname><given-names>SC</given-names></name>, <name><surname>Ferguson</surname><given-names>M</given-names></name>, <name><surname>Kellogg</surname><given-names>MD</given-names></name>, <etal>et al</etal><article-title>Mutations in mitochondrial carrier family gene SLC25A38 cause nonsyndromic autosomal recessive congenital sideroblastic anemia</article-title>. <source>Nat Genet</source>. <year>2009</year>;<volume>41</volume>: <fpage>651</fpage>–<lpage>653</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/ng.359">10.1038/ng.359</ext-link></comment><pub-id pub-id-type="pmid">19412178</pub-id></mixed-citation></ref><ref id="pone.0180432.ref028"><label>28</label><mixed-citation publication-type="journal"><name><surname>May</surname><given-names>BK</given-names></name>, <name><surname>Dogra</surname><given-names>SC</given-names></name>, <name><surname>Sadlon</surname><given-names>TJ</given-names></name>, <name><surname>Bhasker</surname><given-names>CR</given-names></name>, <name><surname>Cox</surname><given-names>TC</given-names></name>, <name><surname>Bottomley</surname><given-names>SS</given-names></name>. <article-title>Molecular regulation of heme biosynthesis in higher vertebrates</article-title>. <source>Prog Nucleic Acid Res Mol Biol</source>. <year>1995</year>;<volume>51</volume>: <fpage>1</fpage>–<lpage>51</lpage>. <pub-id pub-id-type="pmid">7659773</pub-id></mixed-citation></ref><ref id="pone.0180432.ref029"><label>29</label><mixed-citation publication-type="journal"><name><surname>Takahashi</surname><given-names>S</given-names></name>, <name><surname>Taketani</surname><given-names>S</given-names></name>, <name><surname>Akasaka</surname><given-names>J-e</given-names></name>, <name><surname>Kobayashi</surname><given-names>A</given-names></name>, <name><surname>Hayashi</surname><given-names>N</given-names></name>, <name><surname>Yamamoto</surname><given-names>M</given-names></name>, <etal>et al</etal><article-title>Differential regulation of <italic>coproporphyrinogen oxidase</italic> gene between erythroid and nonerythroid cells</article-title>. <source>Blood</source>. <year>1998</year>;<volume>92</volume>: <fpage>3436</fpage>–<lpage>3444</lpage>. <pub-id pub-id-type="pmid">9787185</pub-id></mixed-citation></ref><ref id="pone.0180432.ref030"><label>30</label><mixed-citation publication-type="journal"><name><surname>Day</surname><given-names>AL</given-names></name>, <name><surname>Parsons</surname><given-names>BM</given-names></name>, <name><surname>Dailey</surname><given-names>HA</given-names></name>. <article-title>Cloning and characterization of Gallus and Xenopus ferrochelatases: Presence of the [2Fe-2S] cluster in nonmammalian ferrochelatase</article-title>. <source>Arch Biochem Biophys</source>. <year>1998</year>;<volume>359</volume>: <fpage>160</fpage>–<lpage>169</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1006/abbi.1998.0910">10.1006/abbi.1998.0910</ext-link></comment><pub-id pub-id-type="pmid">9808757</pub-id></mixed-citation></ref><ref id="pone.0180432.ref031"><label>31</label><mixed-citation publication-type="journal"><name><surname>Harbin</surname><given-names>BM</given-names></name>, <name><surname>Dailey</surname><given-names>HA</given-names></name>. <article-title>Orientation of ferrochelatase in bovine liver mitochondria</article-title>. <source>Biochemistry (Mosc)</source>. <year>1985</year>;<volume>24</volume>: <fpage>366</fpage>–<lpage>370</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref032"><label>32</label><mixed-citation publication-type="journal"><name><surname>Kuchipudi</surname><given-names>SV</given-names></name>, <name><surname>Tellabati</surname><given-names>M</given-names></name>, <name><surname>Sebastian</surname><given-names>S</given-names></name>, <name><surname>Londt</surname><given-names>BZ</given-names></name>, <name><surname>Jansen</surname><given-names>C</given-names></name>, <name><surname>Vervelde</surname><given-names>L</given-names></name>, <etal>et al</etal><article-title>Highly pathogenic avian influenza virus infection in chickens but not ducks is associated with elevated host immune and pro-inflammatory responses</article-title>. <source>Vet Res</source>. <year>2014</year>;<volume>45</volume>: <fpage>118</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/s13567-014-0118-3">10.1186/s13567-014-0118-3</ext-link></comment><pub-id pub-id-type="pmid">25431115</pub-id></mixed-citation></ref><ref id="pone.0180432.ref033"><label>33</label><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>F</given-names></name>, <name><surname>Lei</surname><given-names>X</given-names></name>, <name><surname>Rodriguez-Palacios</surname><given-names>A</given-names></name>, <name><surname>Tang</surname><given-names>C</given-names></name>, <name><surname>Yue</surname><given-names>H</given-names></name>. <article-title>Selection of reference genes for quantitative real-time PCR analysis in chicken embryo fibroblasts infected with avian leukosis virus subgroup J</article-title>. <source>BMC Res Notes</source>. <year>2013</year>;<volume>6</volume>: <fpage>402</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1756-0500-6-402">10.1186/1756-0500-6-402</ext-link></comment><pub-id pub-id-type="pmid">24099561</pub-id></mixed-citation></ref><ref id="pone.0180432.ref034"><label>34</label><mixed-citation publication-type="journal"><name><surname>Dalgaard</surname><given-names>TS</given-names></name>, <name><surname>Skovgaard</surname><given-names>K</given-names></name>, <name><surname>Norup</surname><given-names>LR</given-names></name>, <name><surname>Pleidrup</surname><given-names>J</given-names></name>, <name><surname>Permin</surname><given-names>A</given-names></name>, <name><surname>Schou</surname><given-names>TW</given-names></name>, <etal>et al</etal><article-title>Immune gene expression in the spleen of chickens experimentally infected with <italic>Ascaridia galli</italic></article-title>. <source>Vet Immunol Immunopathol</source>. <year>2015</year>;<volume>164</volume>: <fpage>79</fpage>–<lpage>86</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.vetimm.2015.01.003">10.1016/j.vetimm.2015.01.003</ext-link></comment><pub-id pub-id-type="pmid">25649508</pub-id></mixed-citation></ref><ref id="pone.0180432.ref035"><label>35</label><mixed-citation publication-type="journal"><name><surname>Yin</surname><given-names>R</given-names></name>, <name><surname>Liu</surname><given-names>X</given-names></name>, <name><surname>Liu</surname><given-names>C</given-names></name>, <name><surname>Ding</surname><given-names>Z</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name>, <name><surname>Tian</surname><given-names>F</given-names></name>, <etal>et al</etal><article-title>Systematic selection of housekeeping genes for gene expression normalization in chicken embryo fibroblasts infected with Newcastle disease virus</article-title>. <source>Biochem Biophys Res Commun</source>. <year>2011</year>;<volume>413</volume>: <fpage>537</fpage>–<lpage>540</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.bbrc.2011.08.131">10.1016/j.bbrc.2011.08.131</ext-link></comment><pub-id pub-id-type="pmid">21925148</pub-id></mixed-citation></ref><ref id="pone.0180432.ref036"><label>36</label><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>G</given-names></name>, <name><surname>Chen</surname><given-names>S</given-names></name>, <name><surname>Duan</surname><given-names>Z</given-names></name>, <name><surname>Qu</surname><given-names>L</given-names></name>, <name><surname>Xu</surname><given-names>G</given-names></name>, <name><surname>Yang</surname><given-names>N</given-names></name>. <article-title>Comparison of protoporphyrin IX content and related gene expression in the tissues of chickens laying brown-shelled eggs</article-title>. <source>Poult Sci</source>. <year>2013</year>;<volume>92</volume>: <fpage>3120</fpage>–<lpage>3124</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3382/ps.2013-03484">10.3382/ps.2013-03484</ext-link></comment><pub-id pub-id-type="pmid">24235220</pub-id></mixed-citation></ref><ref id="pone.0180432.ref037"><label>37</label><mixed-citation publication-type="journal"><name><surname>Samiullah</surname><given-names>S</given-names></name>, <name><surname>Roberts</surname><given-names>J</given-names></name>. <article-title>The location of protoporphyrin in the eggshell of brown-shelled eggs</article-title>. <source>Poult Sci</source>. <year>2013</year>;<volume>92</volume>: <fpage>2783</fpage>–<lpage>2788</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3382/ps.2013-03051">10.3382/ps.2013-03051</ext-link></comment><pub-id pub-id-type="pmid">24046428</pub-id></mixed-citation></ref><ref id="pone.0180432.ref038"><label>38</label><mixed-citation publication-type="other">Bio-Rad Application User Guide. (2006) Available at: <ext-link ext-link-type="uri" xlink:href="http://www.bio-rad.com/webroot/web/pdf/lsr/literature/Bulletin_5279.pdf">http://www.bio-rad.com/webroot/web/pdf/lsr/literature/Bulletin_5279.pdf</ext-link>. (Accessed: 4th June 2016).</mixed-citation></ref><ref id="pone.0180432.ref039"><label>39</label><mixed-citation publication-type="journal"><name><surname>Hellemans</surname><given-names>J</given-names></name>, <name><surname>Mortier</surname><given-names>G</given-names></name>, <name><surname>De Paepe</surname><given-names>A</given-names></name>, <name><surname>Speleman</surname><given-names>F</given-names></name>, <name><surname>Vandesompele</surname><given-names>J</given-names></name>. <article-title>qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data</article-title>. <source>Genome Biol</source>. <year>2007</year>;<volume>8</volume>: <fpage>R19</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/gb-2007-8-2-r19">10.1186/gb-2007-8-2-r19</ext-link></comment><pub-id pub-id-type="pmid">17291332</pub-id></mixed-citation></ref><ref id="pone.0180432.ref040"><label>40</label><mixed-citation publication-type="journal"><name><surname>Radonić</surname><given-names>A</given-names></name>, <name><surname>Thulke</surname><given-names>S</given-names></name>, <name><surname>Bae</surname><given-names>H-G</given-names></name>, <name><surname>Müller</surname><given-names>MA</given-names></name>, <name><surname>Siegert</surname><given-names>W</given-names></name>, <name><surname>Nitsche</surname><given-names>A</given-names></name>. <article-title>Reference gene selection for quantitative real-time PCR analysis in virus infected cells: SARS corona virus, Yellow fever virus, Human Herpesvirus-6, Camelpox virus and Cytomegalovirus infections</article-title>. <source>Virol J</source>. <year>2005</year>;<volume>2</volume>: <fpage>7</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1743-422X-2-7">10.1186/1743-422X-2-7</ext-link></comment><pub-id pub-id-type="pmid">15705200</pub-id></mixed-citation></ref><ref id="pone.0180432.ref041"><label>41</label><mixed-citation publication-type="journal"><name><surname>Pfaffl</surname><given-names>MW</given-names></name>. <article-title>A new mathematical model for relative quantification in real-time RT—PCR</article-title>. <source>Nucleic Acids Res</source>. <year>2001</year>;<volume>29</volume>: <fpage>e45</fpage><pub-id pub-id-type="pmid">11328886</pub-id></mixed-citation></ref><ref id="pone.0180432.ref042"><label>42</label><mixed-citation publication-type="journal"><name><surname>Ji</surname><given-names>H</given-names></name>, <name><surname>Wang</surname><given-names>J</given-names></name>, <name><surname>Liu</surname><given-names>J</given-names></name>, <name><surname>Guo</surname><given-names>J</given-names></name>, <name><surname>Wang</surname><given-names>Z</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name>, <etal>et al</etal><article-title>Selection of reliable reference genes for real-time qRT-PCR analysis of Zi geese (<italic>Anser anser domestica</italic>) gene expression</article-title>. <source>Asian Austral J Anim</source>. <year>2013</year>;<volume>26</volume>: <fpage>423</fpage>–<lpage>432</lpage>.</mixed-citation></ref><ref id="pone.0180432.ref043"><label>43</label><mixed-citation publication-type="journal"><name><surname>Goossens</surname><given-names>K</given-names></name>, <name><surname>Van Poucke</surname><given-names>M</given-names></name>, <name><surname>Van Soom</surname><given-names>A</given-names></name>, <name><surname>Vandesompele</surname><given-names>J</given-names></name>, <name><surname>Van Zeveren</surname><given-names>A</given-names></name>, <name><surname>Peelman</surname><given-names>LJ</given-names></name>. <article-title>Selection of reference genes for quantitative real-time PCR in bovine preimplantation embryos</article-title>. <source>BMC Dev Biol</source>. <year>2005</year>;<volume>5</volume>: <fpage>27</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1471-213X-5-27">10.1186/1471-213X-5-27</ext-link></comment><pub-id pub-id-type="pmid">16324220</pub-id></mixed-citation></ref><ref id="pone.0180432.ref044"><label>44</label><mixed-citation publication-type="journal"><name><surname>Svingen</surname><given-names>T</given-names></name>, <name><surname>Letting</surname><given-names>H</given-names></name>, <name><surname>Hadrup</surname><given-names>N</given-names></name>, <name><surname>Hass</surname><given-names>U</given-names></name>, <name><surname>Vinggaard</surname><given-names>AM</given-names></name>. <article-title>Selection of reference genes for quantitative RT-PCR (RT-qPCR) analysis of rat tissues under physiological and toxicological conditions</article-title>. <source>PeerJ</source>. <year>2015</year>;<volume>3</volume>: <fpage>e855</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.7717/peerj.855">10.7717/peerj.855</ext-link></comment><pub-id pub-id-type="pmid">25825680</pub-id></mixed-citation></ref><ref id="pone.0180432.ref045"><label>45</label><mixed-citation publication-type="journal"><name><surname>Nascimento</surname><given-names>CS</given-names></name>, <name><surname>Barbosa</surname><given-names>LT</given-names></name>, <name><surname>Brito</surname><given-names>C</given-names></name>, <name><surname>Fernandes</surname><given-names>RP</given-names></name>, <name><surname>Mann</surname><given-names>RS</given-names></name>, <name><surname>Pinto</surname><given-names>APG</given-names></name>, <etal>et al</etal><article-title>Identification of suitable reference genes for real time quantitative polymerase chain reaction assays on Pectoralis major muscle in chicken (<italic>Gallus gallus</italic>)</article-title>. <source>PLoS ONE</source>. <year>2015</year>;<volume>10</volume>: <fpage>e0127935</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1371/journal.pone.0127935">10.1371/journal.pone.0127935</ext-link></comment><pub-id pub-id-type="pmid">26020643</pub-id></mixed-citation></ref><ref id="pone.0180432.ref046"><label>46</label><mixed-citation publication-type="journal"><name><surname>Cinar</surname><given-names>MU</given-names></name>, <name><surname>Islam</surname><given-names>MA</given-names></name>, <name><surname>Pröll</surname><given-names>M</given-names></name>, <name><surname>Kocamis</surname><given-names>H</given-names></name>, <name><surname>Tholen</surname><given-names>E</given-names></name>, <name><surname>Tesfaye</surname><given-names>D</given-names></name>, <etal>et al</etal><article-title>Evaluation of suitable reference genes for gene expression studies in porcine PBMCs in response to LPS and LTA</article-title>. <source>BMC Res Notes</source>. <year>2013</year>;<volume>6</volume>: <fpage>56</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1186/1756-0500-6-56">10.1186/1756-0500-6-56</ext-link></comment><pub-id pub-id-type="pmid">23394600</pub-id></mixed-citation></ref><ref id="pone.0180432.ref047"><label>47</label><mixed-citation publication-type="journal"><name><surname>Pereira-Fantini</surname><given-names>PM</given-names></name>, <name><surname>Rajapaksa</surname><given-names>AE</given-names></name>, <name><surname>Oakley</surname><given-names>R</given-names></name>, <name><surname>Tingay</surname><given-names>DG</given-names></name>. <article-title>Selection of reference genes for gene expression studies related to lung injury in a preterm lamb model</article-title>. <source>Sci Rep</source>. <year>2016</year>;<volume>6</volume>.</mixed-citation></ref><ref id="pone.0180432.ref048"><label>48</label><mixed-citation publication-type="journal"><name><surname>Shivhare</surname><given-names>R</given-names></name>, <name><surname>Lata</surname><given-names>C</given-names></name>. <article-title>Selection of suitable reference genes for assessing gene expression in pearl millet under different abiotic stresses and their combinations</article-title>. <source>Sci Rep</source>. <year>2016</year>;<volume>6</volume>: <fpage>23036</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/srep23036">10.1038/srep23036</ext-link></comment><pub-id pub-id-type="pmid">26972345</pub-id></mixed-citation></ref><ref id="pone.0180432.ref049"><label>49</label><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>Z-J</given-names></name>, <name><surname>Tian</surname><given-names>C</given-names></name>, <name><surname>Jiang</surname><given-names>Q</given-names></name>, <name><surname>Li</surname><given-names>X-H</given-names></name>, <name><surname>Zhuang</surname><given-names>J</given-names></name>. <article-title>Selection of suitable reference genes for qRT-PCR normalization during leaf development and hormonal stimuli in tea plant (<italic>Camellia sinensis</italic>)</article-title>. <source>Sci Rep</source>. <year>2016</year>;<volume>6</volume>: <fpage>19748</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/srep19748">10.1038/srep19748</ext-link></comment><pub-id pub-id-type="pmid">26813576</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/CovidV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000993 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000993 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Sante
|area= CovidV2
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:5495395
|texte= Reference gene selection for the shell gland of laying hens in response to time-points of eggshell formation and nicarbazin
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:28671969" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a CovidV2
| This area was generated with Dilib version V0.6.33. |  |