Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Cryptic sequence features in the active postmortem transcriptome</title><author><name sortKey="Noble, Peter A" sort="Noble, Peter A" uniqKey="Noble P" first="Peter A." last="Noble">Peter A. Noble</name><affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000122986657</institution-id><institution-id institution-id-type="GRID">grid.34477.33</institution-id><institution>Department of Periodontics,</institution><institution>University of Washington,</institution></institution-wrap>Box 357444, Seattle, WA 98195 USA</nlm:aff></affiliation></author><author><name sortKey="Pozhitkov, Alexander E" sort="Pozhitkov, Alexander E" uniqKey="Pozhitkov A" first="Alexander E." last="Pozhitkov">Alexander E. Pozhitkov</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0421 8357</institution-id><institution-id institution-id-type="GRID">grid.410425.6</institution-id><institution>City of Hope, Information Sciences - Beckman Research Institute,</institution></institution-wrap>4920 Rivergrade Rd., Irwindale, CA 91706 USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">30217147</idno><idno type="pmc">6137749</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC6137749</idno><idno type="RBID">PMC:6137749</idno><idno type="doi">10.1186/s12864-018-5042-x</idno><date when="2018">2018</date><idno type="wicri:Area/Pmc/Corpus">000297</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000297</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Cryptic sequence features in the active postmortem transcriptome</title><author><name sortKey="Noble, Peter A" sort="Noble, Peter A" uniqKey="Noble P" first="Peter A." last="Noble">Peter A. Noble</name><affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000122986657</institution-id><institution-id institution-id-type="GRID">grid.34477.33</institution-id><institution>Department of Periodontics,</institution><institution>University of Washington,</institution></institution-wrap>Box 357444, Seattle, WA 98195 USA</nlm:aff></affiliation></author><author><name sortKey="Pozhitkov, Alexander E" sort="Pozhitkov, Alexander E" uniqKey="Pozhitkov A" first="Alexander E." last="Pozhitkov">Alexander E. Pozhitkov</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0421 8357</institution-id><institution-id institution-id-type="GRID">grid.410425.6</institution-id><institution>City of Hope, Information Sciences - Beckman Research Institute,</institution></institution-wrap>4920 Rivergrade Rd., Irwindale, CA 91706 USA</nlm:aff></affiliation></author></analytic><series><title level="j">BMC Genomics</title><idno type="eISSN">1471-2164</idno><imprint><date when="2018">2018</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p id="Par1">Our previous study found that more than 500 transcripts significantly increased in abundance in the zebrafish and mouse several hours to days postmortem relative to live controls. The current literature suggests that most mRNAs are post-transcriptionally regulated in stressful conditions. We rationalized that the postmortem transcripts must contain sequence features (3- to 9- mers) that are unique from those in the rest of the transcriptome and that these features putatively serve as binding sites for proteins and/or non-coding RNAs involved in post-transcriptional regulation.</p></sec><sec><title>Results</title><p id="Par2">We identified 5117 and 2245 over-represented sequence features in the mouse and zebrafish, respectively, which represents less than 1.5% of all possible features. Some of these features were disproportionately distributed along the transcripts with high densities in the 3′ untranslated regions of the zebrafish (0.3 mers/nt) and the open reading frames of the mouse (0.6 mers/nt). Yet, the highest density (2.3 mers/nt) occurred in the open reading frames of 11 mouse transcripts that lacked 3′ or 5′ untranslated regions. These results suggest the transcripts with high density of features might serve as ‘molecular sponges’ that sequester RNA binding proteins and/or microRNAs, and thus indirectly increase the stability and gene expression of other transcripts. In addition, some of the features were identified as binding sites for <italic>Rbfox</italic> and <italic>Hud</italic> proteins that are also involved in increasing transcript stability and gene expression.</p></sec><sec><title>Conclusions</title><p id="Par3">Our results are consistent with the hypothesis that transcripts involved in responding to extreme stress, such as organismal death, have sequence features that make them different from the rest of the transcriptome. Some of these features serve as putative binding sites for proteins and non-coding RNAs that determine transcript stability and fate. A small number of the transcripts have high density sequence features, which are presumably involved in sequestering RNA binding proteins and microRNAs and thus preventing regulatory interactions among other transcripts. Our results provide baseline data on post-transcriptional regulation in stressful conditions that has implications for regulation in disease, starvation, and cancer.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (10.1186/s12864-018-5042-x) contains supplementary material, which is available to authorized users.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Harvey, R" uniqKey="Harvey R">R Harvey</name></author><author><name sortKey="Dezi, V" uniqKey="Dezi V">V Dezi</name></author><author><name sortKey="Pizzinga, M" uniqKey="Pizzinga M">M Pizzinga</name></author><author><name sortKey="Willis, Ae" uniqKey="Willis A">AE Willis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hogg, Jr" uniqKey="Hogg J">JR Hogg</name></author><author><name sortKey="Collins, K" uniqKey="Collins K">K Collins</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ferreira, Pg" uniqKey="Ferreira P">PG Ferreira</name></author><author><name sortKey="Mu Oz Aguirre, M" uniqKey="Mu Oz Aguirre M">M Muñoz-Aguirre</name></author><author><name sortKey="Reverter, F" uniqKey="Reverter F">F Reverter</name></author><author><name sortKey="Sa Godinho, Cp" uniqKey="Sa Godinho C">CP Sá Godinho</name></author><author><name sortKey="Sousa, A" uniqKey="Sousa A">A Sousa</name></author><author><name sortKey="Amadoz, A" uniqKey="Amadoz A">A Amadoz</name></author><author><name sortKey="Sodaei, R" uniqKey="Sodaei R">R Sodaei</name></author><author><name sortKey="Hidalgo, Mr" uniqKey="Hidalgo M">MR Hidalgo</name></author><author><name sortKey="Pervouchine" uniqKey="Pervouchine">Pervouchine</name></author><author><name sortKey="Dmitri" uniqKey="Dmitri">Dmitri</name></author><author><name sortKey="Carbonell Caballero, J" uniqKey="Carbonell Caballero J">J Carbonell-Caballero</name></author><author><name sortKey="Nurtdinov, R" uniqKey="Nurtdinov R">R Nurtdinov</name></author><author><name sortKey="Breschi, A" uniqKey="Breschi A">A Breschi</name></author><author><name sortKey="Amador, R" uniqKey="Amador R">R Amador</name></author><author><name sortKey="Oliveira, P" uniqKey="Oliveira P">P Oliveira</name></author><author><name sortKey="Cubuk, C" uniqKey="Cubuk C">C Çubuk</name></author><author><name sortKey="Curado, J" uniqKey="Curado J">J Curado</name></author><author><name sortKey="Aguet, F" uniqKey="Aguet F">F Aguet</name></author><author><name sortKey="Oliveira, C" uniqKey="Oliveira C">C Oliveira</name></author><author><name sortKey="Dopazo, J" uniqKey="Dopazo J">J Dopazo</name></author><author><name sortKey="Sammeth, M" uniqKey="Sammeth M">M Sammeth</name></author><author><name sortKey="Ardlie, Kg" uniqKey="Ardlie K">KG Ardlie</name></author><author><name sortKey="Guig, R" uniqKey="Guig R">R Guigó</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pozhitkov, Ae" uniqKey="Pozhitkov A">AE Pozhitkov</name></author><author><name sortKey="Neme, R" uniqKey="Neme R">R Neme</name></author><author><name sortKey="Domazet Loo, T" uniqKey="Domazet Loo T">T Domazet-Loo</name></author><author><name sortKey="Leroux, Bg" uniqKey="Leroux B">BG Leroux</name></author><author><name sortKey="Soni, S" uniqKey="Soni S">S Soni</name></author><author><name sortKey="Tautz, D" uniqKey="Tautz D">D Tautz</name></author><author><name sortKey="Noble, Pa" uniqKey="Noble P">PA Noble</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dykes, Im" uniqKey="Dykes I">IM Dykes</name></author><author><name sortKey="Emanueli, C" uniqKey="Emanueli C">C Emanueli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="He, Jh" uniqKey="He J">JH He</name></author><author><name sortKey="Han, Zp" uniqKey="Han Z">ZP Han</name></author><author><name sortKey="Liu, Jm" uniqKey="Liu J">JM Liu</name></author><author><name sortKey="Zhou, Jb" uniqKey="Zhou J">JB Zhou</name></author><author><name sortKey="Zou, Mx" uniqKey="Zou M">MX Zou</name></author><author><name sortKey="Lv, Yb" uniqKey="Lv Y">YB Lv</name></author><author><name sortKey="Li, Yg" uniqKey="Li Y">YG Li</name></author><author><name sortKey="Cao, Mr" uniqKey="Cao M">MR Cao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gerstberger, S" uniqKey="Gerstberger S">S Gerstberger</name></author><author><name sortKey="Hafner, M" uniqKey="Hafner M">M Hafner</name></author><author><name sortKey="Tuschl, T" uniqKey="Tuschl T">T Tuschl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maraia, Richard J" uniqKey="Maraia R">Richard J. Maraia</name></author><author><name sortKey="Mattijssen, Sandy" uniqKey="Mattijssen S">Sandy Mattijssen</name></author><author><name sortKey="Cruz Gallardo, Isabel" uniqKey="Cruz Gallardo I">Isabel Cruz-Gallardo</name></author><author><name sortKey="Conte, Maria R" uniqKey="Conte M">Maria R. Conte</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Castella, S" uniqKey="Castella S">S Castella</name></author><author><name sortKey="Bernard, R" uniqKey="Bernard R">R Bernard</name></author><author><name sortKey="Corno, M" uniqKey="Corno M">M Corno</name></author><author><name sortKey="Fradin, A" uniqKey="Fradin A">A Fradin</name></author><author><name sortKey="Larcher, Jc" uniqKey="Larcher J">JC Larcher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Conboy, John G" uniqKey="Conboy J">John G. Conboy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Soengas, Ms" uniqKey="Soengas M">MS Soengas</name></author><author><name sortKey="Hernando, E" uniqKey="Hernando E">E Hernando</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gilot, D" uniqKey="Gilot D">D Gilot</name></author><author><name sortKey="Migault, M" uniqKey="Migault M">M Migault</name></author><author><name sortKey="Bachelot, L" uniqKey="Bachelot L">L Bachelot</name></author><author><name sortKey="Journe, F" uniqKey="Journe F">F Journé</name></author><author><name sortKey="Rogiers, A" uniqKey="Rogiers A">A Rogiers</name></author><author><name sortKey="Donnou Fournet, E" uniqKey="Donnou Fournet E">E Donnou-Fournet</name></author><author><name sortKey="Mogha, A" uniqKey="Mogha A">A Mogha</name></author><author><name sortKey="Mouchet, N" uniqKey="Mouchet N">N Mouchet</name></author><author><name sortKey="Pinel Marie, Ml" uniqKey="Pinel Marie M">ML Pinel-Marie</name></author><author><name sortKey="Mari, B" uniqKey="Mari B">B Mari</name></author><author><name sortKey="Montier, T" uniqKey="Montier T">T Montier</name></author><author><name sortKey="Corre, S" uniqKey="Corre S">S Corre</name></author><author><name sortKey="Gautron, A" uniqKey="Gautron A">A Gautron</name></author><author><name sortKey="Rambow, F" uniqKey="Rambow F">F Rambow</name></author><author><name sortKey="El Hajj, P" uniqKey="El Hajj P">P El Hajj</name></author><author><name sortKey="Ben Jouira, R" uniqKey="Ben Jouira R">R Ben Jouira</name></author><author><name sortKey="Tartare Deckert, S" uniqKey="Tartare Deckert S">S Tartare-Deckert</name></author><author><name sortKey="Marine, Jc" uniqKey="Marine J">JC Marine</name></author><author><name sortKey="Felden, B" uniqKey="Felden B">B Felden</name></author><author><name sortKey="Ghanem, G" uniqKey="Ghanem G">G Ghanem</name></author><author><name sortKey="Galibert, Md" uniqKey="Galibert M">MD Galibert</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Labno, A" uniqKey="Labno A">A Łabno</name></author><author><name sortKey="Tomecki, R" uniqKey="Tomecki R">R Tomecki</name></author><author><name sortKey="Dziembowski, A" uniqKey="Dziembowski A">A Dziembowski</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wilusz, Je" uniqKey="Wilusz J">JE Wilusz</name></author><author><name sortKey="Sunwoo, H" uniqKey="Sunwoo H">H Sunwoo</name></author><author><name sortKey="Spector, Dl" uniqKey="Spector D">DL Spector</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Eddy, Sr" uniqKey="Eddy S">SR Eddy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jeffrey, Hj" uniqKey="Jeffrey H">HJ Jeffrey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Noble, Pa" uniqKey="Noble P">PA Noble</name></author><author><name sortKey="Citek, Rw" uniqKey="Citek R">RW Citek</name></author><author><name sortKey="Ogunseitan, Oa" uniqKey="Ogunseitan O">OA Ogunseitan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Almeida, Js" uniqKey="Almeida J">JS Almeida</name></author><author><name sortKey="Carrico, Ja" uniqKey="Carrico J">JA Carriço</name></author><author><name sortKey="Maretzek, A" uniqKey="Maretzek A">A Maretzek</name></author><author><name sortKey="Noble, Pa" uniqKey="Noble P">PA Noble</name></author><author><name sortKey="Fletcher, M" uniqKey="Fletcher M">M Fletcher</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Bolognani, F" uniqKey="Bolognani F">F Bolognani</name></author><author><name sortKey="Contente Cuomo, T" uniqKey="Contente Cuomo T">T Contente-Cuomo</name></author><author><name sortKey="Perrone Bizzozero, Ni" uniqKey="Perrone Bizzozero N">NI Perrone-Bizzozero</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jacobsen, A" uniqKey="Jacobsen A">A Jacobsen</name></author><author><name sortKey="Wen, J" uniqKey="Wen J">J Wen</name></author><author><name sortKey="Marks, Ds" uniqKey="Marks D">DS Marks</name></author><author><name sortKey="Krogh, A" uniqKey="Krogh A">A Krogh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Barreau, C" uniqKey="Barreau C">C Barreau</name></author><author><name sortKey="Paillard, L" uniqKey="Paillard L">L Paillard</name></author><author><name sortKey="Osborne, Hb" uniqKey="Osborne H">HB Osborne</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ding, Cl" uniqKey="Ding C">CL Ding</name></author><author><name sortKey="Xu, G" uniqKey="Xu G">G Xu</name></author><author><name sortKey="Tang, Hl" uniqKey="Tang H">HL Tang</name></author><author><name sortKey="Zhu, Sy" uniqKey="Zhu S">SY Zhu</name></author><author><name sortKey="Zhao, Lj" uniqKey="Zhao L">LJ Zhao</name></author><author><name sortKey="Ren, H" uniqKey="Ren H">H Ren</name></author><author><name sortKey="Zhao, P" uniqKey="Zhao P">P Zhao</name></author><author><name sortKey="Qi, Zt" uniqKey="Qi Z">ZT Qi</name></author><author><name sortKey="Wang, W" uniqKey="Wang W">W Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miller, Lc" uniqKey="Miller L">LC Miller</name></author><author><name sortKey="Jiang, Z" uniqKey="Jiang Z">Z Jiang</name></author><author><name sortKey="Sang, Y" uniqKey="Sang Y">Y Sang</name></author><author><name sortKey="Harhay, Gp" uniqKey="Harhay G">GP Harhay</name></author><author><name sortKey="Lager, Km" uniqKey="Lager K">KM Lager</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cruse, G" uniqKey="Cruse G">G Cruse</name></author><author><name sortKey="Beaven, Ma" uniqKey="Beaven M">MA Beaven</name></author><author><name sortKey="Music, Sc" uniqKey="Music S">SC Music</name></author><author><name sortKey="Bradding, P" uniqKey="Bradding P">P Bradding</name></author><author><name sortKey="Gilfillan, Am" uniqKey="Gilfillan A">AM Gilfillan</name></author><author><name sortKey="Metcalfe, Dd" uniqKey="Metcalfe D">DD Metcalfe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bens, S" uniqKey="Bens S">S Bens</name></author><author><name sortKey="Zichner, T" uniqKey="Zichner T">T Zichner</name></author><author><name sortKey="Stutz, Am" uniqKey="Stutz A">AM Stütz</name></author><author><name sortKey="Caliebe, A" uniqKey="Caliebe A">A Caliebe</name></author><author><name sortKey="Wagener, R" uniqKey="Wagener R">R Wagener</name></author><author><name sortKey="Hoff, K" uniqKey="Hoff K">K Hoff</name></author><author><name sortKey="Korbel, Jo" uniqKey="Korbel J">JO Korbel</name></author><author><name sortKey="Von Bismarck, P" uniqKey="Von Bismarck P">P von Bismarck</name></author><author><name sortKey="Siebert, R" uniqKey="Siebert R">R Siebert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fernebro, J" uniqKey="Fernebro J">J Fernebro</name></author><author><name sortKey="Francis, P" uniqKey="Francis P">P Francis</name></author><author><name sortKey="Eden, P" uniqKey="Eden P">P Edén</name></author><author><name sortKey="Borg, A" uniqKey="Borg A">A Borg</name></author><author><name sortKey="Panagopoulos, I" uniqKey="Panagopoulos I">I Panagopoulos</name></author><author><name sortKey="Mertens, F" uniqKey="Mertens F">F Mertens</name></author><author><name sortKey="Vallon Christersson, J" uniqKey="Vallon Christersson J">J Vallon-Christersson</name></author><author><name sortKey="Akerman, M" uniqKey="Akerman M">M Akerman</name></author><author><name sortKey="Rydholm, A" uniqKey="Rydholm A">A Rydholm</name></author><author><name sortKey="Bauer, Hc" uniqKey="Bauer H">HC Bauer</name></author><author><name sortKey="Mandahl, N" uniqKey="Mandahl N">N Mandahl</name></author><author><name sortKey="Nilbert, M" uniqKey="Nilbert M">M Nilbert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kerr, Jr" uniqKey="Kerr J">JR Kerr</name></author><author><name sortKey="Kaushik, N" uniqKey="Kaushik N">N Kaushik</name></author><author><name sortKey="Fear, D" uniqKey="Fear D">D Fear</name></author><author><name sortKey="Baldwin, Da" uniqKey="Baldwin D">DA Baldwin</name></author><author><name sortKey="Nuwaysir, Ef" uniqKey="Nuwaysir E">EF Nuwaysir</name></author><author><name sortKey="Adcock, Im" uniqKey="Adcock I">IM Adcock</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, L" uniqKey="Chen L">L Chen</name></author><author><name sortKey="Ma, C" uniqKey="Ma C">C Ma</name></author><author><name sortKey="Bian, Y" uniqKey="Bian Y">Y Bian</name></author><author><name sortKey="Shao, C" uniqKey="Shao C">C Shao</name></author><author><name sortKey="Wang, T" uniqKey="Wang T">T Wang</name></author><author><name sortKey="Li, J" uniqKey="Li J">J Li</name></author><author><name sortKey="Chong, X" uniqKey="Chong X">X Chong</name></author><author><name sortKey="Su, L" uniqKey="Su L">L Su</name></author><author><name sortKey="Lu, J" uniqKey="Lu J">J Lu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, Yunfei" uniqKey="Zhao Y">Yunfei Zhao</name></author><author><name sortKey="Yang, Ling" uniqKey="Yang L">Ling Yang</name></author><author><name sortKey="He, Jiao" uniqKey="He J">Jiao He</name></author><author><name sortKey="Yang, Huai" uniqKey="Yang H">Huai Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, J" uniqKey="Zhou J">J Zhou</name></author><author><name sortKey="Wang, F" uniqKey="Wang F">F Wang</name></author><author><name sortKey="Liu, B" uniqKey="Liu B">B Liu</name></author><author><name sortKey="Yang, L" uniqKey="Yang L">L Yang</name></author><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name></author><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Arch, Rh" uniqKey="Arch R">RH Arch</name></author><author><name sortKey="Thompson, Cb" uniqKey="Thompson C">CB Thompson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mignone, Flavio" uniqKey="Mignone F">Flavio Mignone</name></author><author><name sortKey="Gissi, Carmela" uniqKey="Gissi C">Carmela Gissi</name></author><author><name sortKey="Liuni, Sabino" uniqKey="Liuni S">Sabino Liuni</name></author><author><name sortKey="Pesole, Graziano" uniqKey="Pesole G">Graziano Pesole</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tay, Y" uniqKey="Tay Y">Y Tay</name></author><author><name sortKey="Zhang, J" uniqKey="Zhang J">J Zhang</name></author><author><name sortKey="Thomson, Am" uniqKey="Thomson A">AM Thomson</name></author><author><name sortKey="Lim, B" uniqKey="Lim B">B Lim</name></author><author><name sortKey="Rigoutsos, I" uniqKey="Rigoutsos I">I Rigoutsos</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ray, D" uniqKey="Ray D">D Ray</name></author><author><name sortKey="Kazan, H" uniqKey="Kazan H">H Kazan</name></author><author><name sortKey="Cook, Kb" uniqKey="Cook K">KB Cook</name></author><author><name sortKey="Weirauch, Mt" uniqKey="Weirauch M">MT Weirauch</name></author><author><name sortKey="Najafabadi, Hs" uniqKey="Najafabadi H">HS Najafabadi</name></author><author><name sortKey="Li, X" uniqKey="Li X">X Li</name></author><author><name sortKey="Gueroussov, S" uniqKey="Gueroussov S">S Gueroussov</name></author><author><name sortKey="Albu, M" uniqKey="Albu M">M Albu</name></author><author><name sortKey="Zheng, H" uniqKey="Zheng H">H Zheng</name></author><author><name sortKey="Yang, A" uniqKey="Yang A">A Yang</name></author><author><name sortKey="Na, H" uniqKey="Na H">H Na</name></author><author><name sortKey="Irimia, M" uniqKey="Irimia M">M Irimia</name></author><author><name sortKey="Matzat, Lh" uniqKey="Matzat L">LH Matzat</name></author><author><name sortKey="Dale, Rk" uniqKey="Dale R">RK Dale</name></author><author><name sortKey="Smith, Sa" uniqKey="Smith S">SA Smith</name></author><author><name sortKey="Yarosh, Ca" uniqKey="Yarosh C">CA Yarosh</name></author><author><name sortKey="Kelly, Sm" uniqKey="Kelly S">SM Kelly</name></author><author><name sortKey="Nabet, B" uniqKey="Nabet B">B Nabet</name></author><author><name sortKey="Mecenas, D" uniqKey="Mecenas D">D Mecenas</name></author><author><name sortKey="Li, W" uniqKey="Li W">W Li</name></author><author><name sortKey="Laishram, Rs" uniqKey="Laishram R">RS Laishram</name></author><author><name sortKey="Qiao, M" uniqKey="Qiao M">M Qiao</name></author><author><name sortKey="Lipshitz, Hd" uniqKey="Lipshitz H">HD Lipshitz</name></author><author><name sortKey="Piano, F" uniqKey="Piano F">F Piano</name></author><author><name sortKey="Corbett, Ah" uniqKey="Corbett A">AH Corbett</name></author><author><name sortKey="Carstens, Rp" uniqKey="Carstens R">RP Carstens</name></author><author><name sortKey="Frey, Bj" uniqKey="Frey B">BJ Frey</name></author><author><name sortKey="Anderson, Ra" uniqKey="Anderson R">RA Anderson</name></author><author><name sortKey="Lynch, Kw" uniqKey="Lynch K">KW Lynch</name></author><author><name sortKey="Penalva, Lo" uniqKey="Penalva L">LO Penalva</name></author><author><name sortKey="Lei, Ep" uniqKey="Lei E">EP Lei</name></author><author><name sortKey="Fraser, Ag" uniqKey="Fraser A">AG Fraser</name></author><author><name sortKey="Blencowe, Bj" uniqKey="Blencowe B">BJ Blencowe</name></author><author><name sortKey="Morris, Qd" uniqKey="Morris Q">QD Morris</name></author><author><name sortKey="Hughes, Tr" uniqKey="Hughes T">TR Hughes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Geisberg, Jv" uniqKey="Geisberg J">JV Geisberg</name></author><author><name sortKey="Moqtaderi, Z" uniqKey="Moqtaderi Z">Z Moqtaderi</name></author><author><name sortKey="Fan, X" uniqKey="Fan X">X Fan</name></author><author><name sortKey="Ozsolak, F" uniqKey="Ozsolak F">F Ozsolak</name></author><author><name sortKey="Struhl, K" uniqKey="Struhl K">K Struhl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Frankel, Lb" uniqKey="Frankel L">LB Frankel</name></author><author><name sortKey="Christoffersen, Nr" uniqKey="Christoffersen N">NR Christoffersen</name></author><author><name sortKey="Jacobsen, A" uniqKey="Jacobsen A">A Jacobsen</name></author><author><name sortKey="Lindow, M" uniqKey="Lindow M">M Lindow</name></author><author><name sortKey="Krogh, A" uniqKey="Krogh A">A Krogh</name></author><author><name sortKey="Lund, Ah" uniqKey="Lund A">AH Lund</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Salmena, L" uniqKey="Salmena L">L Salmena</name></author><author><name sortKey="Poliseno, L" uniqKey="Poliseno L">L Poliseno</name></author><author><name sortKey="Tay, Y" uniqKey="Tay Y">Y Tay</name></author><author><name sortKey="Kats, L" uniqKey="Kats L">L Kats</name></author><author><name sortKey="Pandolfi, Pp" uniqKey="Pandolfi P">PP Pandolfi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Plass, M" uniqKey="Plass M">M Plass</name></author><author><name sortKey="Rasmussen, Sh" uniqKey="Rasmussen S">SH Rasmussen</name></author><author><name sortKey="Krogh, A" uniqKey="Krogh A">A Krogh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pozhitkov, Ae" uniqKey="Pozhitkov A">AE Pozhitkov</name></author><author><name sortKey="Noble, Pa" uniqKey="Noble P">PA Noble</name></author><author><name sortKey="Bryk, J" uniqKey="Bryk J">J Bryk</name></author><author><name sortKey="Tautz, D" uniqKey="Tautz D">D Tautz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pozhitkov, Ae" uniqKey="Pozhitkov A">AE Pozhitkov</name></author><author><name sortKey="Noble, Pa" uniqKey="Noble P">PA Noble</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">BMC Genomics</journal-id><journal-id journal-id-type="iso-abbrev">BMC Genomics</journal-id><journal-title-group><journal-title>BMC Genomics</journal-title></journal-title-group><issn pub-type="epub">1471-2164</issn><publisher><publisher-name>BioMed Central</publisher-name><publisher-loc>London</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">30217147</article-id><article-id pub-id-type="pmc">6137749</article-id><article-id pub-id-type="publisher-id">5042</article-id><article-id pub-id-type="doi">10.1186/s12864-018-5042-x</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Cryptic sequence features in the active postmortem transcriptome</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-6013-2588</contrib-id><name><surname>Noble</surname><given-names>Peter A.</given-names></name><address><phone>206-409-6664</phone><email>panoble2017@gmail.com</email></address><xref ref-type="aff" rid="Aff1">1</xref></contrib><contrib contrib-type="author"><name><surname>Pozhitkov</surname><given-names>Alexander E.</given-names></name><address><email>apozhitkov@coh.org</email></address><xref ref-type="aff" rid="Aff2">2</xref></contrib><aff id="Aff1"><label>1</label><institution-wrap><institution-id institution-id-type="ISNI">0000000122986657</institution-id><institution-id institution-id-type="GRID">grid.34477.33</institution-id><institution>Department of Periodontics,</institution><institution>University of Washington,</institution></institution-wrap>Box 357444, Seattle, WA 98195 USA</aff><aff id="Aff2"><label>2</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0421 8357</institution-id><institution-id institution-id-type="GRID">grid.410425.6</institution-id><institution>City of Hope, Information Sciences - Beckman Research Institute,</institution></institution-wrap>4920 Rivergrade Rd., Irwindale, CA 91706 USA</aff></contrib-group><pub-date pub-type="epub"><day>14</day><month>9</month><year>2018</year></pub-date><pub-date pub-type="pmc-release"><day>14</day><month>9</month><year>2018</year></pub-date><pub-date pub-type="collection"><year>2018</year></pub-date><volume>19</volume><elocation-id>675</elocation-id><history><date date-type="received"><day>4</day><month>4</month><year>2018</year></date><date date-type="accepted"><day>27</day><month>8</month><year>2018</year></date></history><permissions><copyright-statement>© The Author(s). 2018</copyright-statement><license license-type="OpenAccess"><license-p><bold>Open Access</bold>This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>) applies to the data made available in this article, unless otherwise stated.</license-p></license></permissions><abstract id="Abs1"><sec><title>Background</title><p id="Par1">Our previous study found that more than 500 transcripts significantly increased in abundance in the zebrafish and mouse several hours to days postmortem relative to live controls. The current literature suggests that most mRNAs are post-transcriptionally regulated in stressful conditions. We rationalized that the postmortem transcripts must contain sequence features (3- to 9- mers) that are unique from those in the rest of the transcriptome and that these features putatively serve as binding sites for proteins and/or non-coding RNAs involved in post-transcriptional regulation.</p></sec><sec><title>Results</title><p id="Par2">We identified 5117 and 2245 over-represented sequence features in the mouse and zebrafish, respectively, which represents less than 1.5% of all possible features. Some of these features were disproportionately distributed along the transcripts with high densities in the 3′ untranslated regions of the zebrafish (0.3 mers/nt) and the open reading frames of the mouse (0.6 mers/nt). Yet, the highest density (2.3 mers/nt) occurred in the open reading frames of 11 mouse transcripts that lacked 3′ or 5′ untranslated regions. These results suggest the transcripts with high density of features might serve as ‘molecular sponges’ that sequester RNA binding proteins and/or microRNAs, and thus indirectly increase the stability and gene expression of other transcripts. In addition, some of the features were identified as binding sites for <italic>Rbfox</italic> and <italic>Hud</italic> proteins that are also involved in increasing transcript stability and gene expression.</p></sec><sec><title>Conclusions</title><p id="Par3">Our results are consistent with the hypothesis that transcripts involved in responding to extreme stress, such as organismal death, have sequence features that make them different from the rest of the transcriptome. Some of these features serve as putative binding sites for proteins and non-coding RNAs that determine transcript stability and fate. A small number of the transcripts have high density sequence features, which are presumably involved in sequestering RNA binding proteins and microRNAs and thus preventing regulatory interactions among other transcripts. Our results provide baseline data on post-transcriptional regulation in stressful conditions that has implications for regulation in disease, starvation, and cancer.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (10.1186/s12864-018-5042-x) contains supplementary material, which is available to authorized users.</p></sec></abstract><kwd-group xml:lang="en"><title>Keywords</title><kwd>Motifs</kwd><kwd>Sequence features</kwd><kwd>Post-transcriptional regulation</kwd><kwd>Stress response</kwd><kwd>Postmortem gene expression</kwd><kwd>Chaos game representation</kwd><kwd>Zebrafish</kwd><kwd>Mouse</kwd><kwd>5’UTR</kwd><kwd>3’UTR</kwd><kwd>ORFs</kwd><kwd>Molecular sponge</kwd></kwd-group><custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name><meta-value>© The Author(s) 2018</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="Sec1"><title>Background</title><p id="Par12">Understanding regulatory circuits and how they influence transcriptional dynamics are important for comprehending the response of biological systems to stress such as starvation, disease, cancer and even death. Under stressful conditions, most (90%) mRNAs are regulated post-transcriptionally [<xref ref-type="bibr" rid="CR1">1</xref>] -- presumably because it is more energetically favorable than regulation at the transcriptional level [<xref ref-type="bibr" rid="CR2">2</xref>].</p><p id="Par13">Two studies have recently shown that hundreds of transcripts increase in abundance in vertebrate organs/tissues in response to organismal death [<xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR4">4</xref>]. These increases could be due to active transcription and/or post-transcriptional regulation of the nascent transcripts. Post-transcriptional regulation involves RNA binding proteins (RBPs) and non-coding RNAs (ncRNAs) [<xref ref-type="bibr" rid="CR5">5</xref>, <xref ref-type="bibr" rid="CR6">6</xref>] that form complexes with RNA motifs and regions of secondary structure within the RNAs [<xref ref-type="bibr" rid="CR7">7</xref>]. While the binding of RBPs to specific motifs in a transcript is well documented [<xref ref-type="bibr" rid="CR8">8</xref>–<xref ref-type="bibr" rid="CR10">10</xref>], the binding of ncRNA, in the form of microRNA (miRNA), circular RNA, or long ncRNA (lncRNA) to specific motifs within transcripts is less understood. Apparently, some mRNAs and ncRNAs act as “molecular sponges” that bind miRNAs preventing them from performing their functions. For example, miRNA-16 is sequestered by mRNAs encoded by the Tyrosinase-related Protein 1 (<italic>Tyrp1</italic>) gene [<xref ref-type="bibr" rid="CR11">11</xref>]. As a consequence, miRNA-16 tumour-suppressor functions are lost and cell proliferation occurs [<xref ref-type="bibr" rid="CR12">12</xref>]. Another “sponge” example is lncRNA encoded by the <italic>Meg3</italic> gene that binds miRNA-664 counteracting its inhibitory effect on production of alcohol dehydrogenase [<xref ref-type="bibr" rid="CR6">6</xref>]. These are examples of two RNAs acting as molecular sponges -- yet, not all of the functions of ncRNAs are known at this time [<xref ref-type="bibr" rid="CR13">13</xref>] -- other roles have been suggested [<xref ref-type="bibr" rid="CR14">14</xref>–<xref ref-type="bibr" rid="CR16">16</xref>].</p><p id="Par14">Our previous study revealed that some transcripts increase in abundance with postmortem time [<xref ref-type="bibr" rid="CR4">4</xref>]. As a step forward towards better understanding of possible mechanisms responsible for these increases, our present study examined sequence features (i.e., short mers) that are over- or under- represented within these transcripts. We recognize that short mers are not the only sequence features responsible for these increases – we begin with short mers because they are easily identified. That said other more complex features are probably yet to be discovered.</p><p id="Par15">We rationalized that some mers are over- or under- represented in these transcripts because they serve as binding sites for RBPs or ncRNAs involved in post-transcriptional regulation. To investigate this phenomenon, we examined the presence/absence/ frequencies of mers up to 9 nt in length and compared them to controls, which consisted of random draws of transcripts from the rest of the transcriptome (i.e., those not increasing in abundance in response to stress). The results show that several thousand mers are over-represented in the postmortem transcripts of the zebrafish and mouse. Further examination of the frequencies of the mers show that some transcripts have more unique mers than others, and that the density of the unique mers varies by transcript and region (i.e., 5’UTR, ORF, 3’UTR).</p></sec><sec id="Sec2"><title>Methods</title><p id="Par16">A schematic overview of the experimental design for the study is provided in Fig. <xref rid="Fig1" ref-type="fig">1</xref><italic>.</italic><fig id="Fig1"><label>Fig. 1</label><caption><p>Schematic representation of the study experimental design</p></caption><graphic xlink:href="12864_2018_5042_Fig1_HTML" id="MO1"></graphic></fig></p><sec id="Sec3"><title>Dataset assembly</title><p id="Par17">Messenger RNA transcripts of <italic>Danio rerio</italic> (GRC210.89) and <italic>Mus musculus</italic> (GRCm38.p5) with annotations were downloaded from NCBI. Transcript sequences containing ambiguous nucleotides (i.e., ‘N’s) and those less than 100 nt in length were removed. The final “clean” data sets were used for bioinformatic analyses.</p></sec><sec id="Sec4"><title>Extracting 2- to 9-mers from transcript sequences</title><p id="Par18">An alignment-free sequence comparison method called ‘Chaos Genome Representation’ (CGR) [<xref ref-type="bibr" rid="CR17">17</xref>–<xref ref-type="bibr" rid="CR19">19</xref>] was used to extract mers from the transcript sequences because it was more practical (computationally efficient) than string-based search algorithms (see Proof in Additional file <xref rid="MOESM1" ref-type="media">1</xref>). CGR is an iterative mapping technique that processes nucleotides in a sequence to find the <italic>x-, y-</italic> coordinates for their position in a continuous space. The <italic>x-</italic> and <italic>y-</italic> coordinates can then be used to recover sequence, which in this study were oligomers. Once the coordinates of a sequence are known, the presence/absence/ frequency of a mer of any size in a transcript sequence can easily be determined, as demonstrated below.</p></sec><sec id="Sec5"><title>Reading a sequence into CGR space</title><p id="Par19">The processing of a transcript sequence involves converting each nucleotide into <italic>x-</italic> and <italic>y-</italic> coordinates and assembling the coordinates into a CGR database. For example, the sequence ‘AAACC’ is represented by the <italic>x</italic>- and <italic>y-</italic> coordinates of + 0.53125 and − 0.53125, respectively. The coordinates are determined by reading the sequence into CGR space. The space is confined by the four possible nucleotides as vertices of a binary square with <italic>x, y</italic> position (− 1, + 1) being the vertex A, (+ 1, + 1) being the vertex T, (− 1, + 1) being the vertex G and (− 1, − 1) being the vertex C. The position of a nucleotide in the fragment is calculated by moving a pointer to half the distance between the previous position and the current binary representation.</p><sec id="Sec6"><title>An example</title><p id="Par20">Starting at point <italic>x, y</italic> (0, 0), the first nucleotide ‘A’ is plotted at half way to the vertex of A (− 1, + 1), which is coordinate (− 0.5, + 0.5). The next nucleotide is also ‘A’, therefore half way from the coordinate (− 0.5, + 0.5) to vertex of A (− 1, + 1) is (− 0.75, + 0.75). The next nucleotide is also ‘A’ so half way from the coordinate (− 0.75, + 0.75) to the vertex of A (− 1, + 1) is the coordinate (− 0.875, + 0.875). The next nucleotide ‘C’, so half-way from the coordinate (− 0.875, + 0.875) to the vertex of C (− 1, − 1) is the coordinate (+ 0.0625, − 0.0625) and so on up to the last nucleotide of the sequence with the last coordinates of <italic>x</italic> = + 0.53125 and <italic>y</italic> = − 0.53125. A depiction of reading a sequence into CGR space is shown in Fig. <xref rid="Fig1" ref-type="fig">1a</xref> of the Almeida et al. [<xref ref-type="bibr" rid="CR19">19</xref>] study.</p><p id="Par21">Once all the sequences have been read into CGR space and their coordinates stored in a database, it is possible to determine the presence/absence/frequency of mers by their coordinates and mer length (i.e., 1/resolution), which is outlined in the Mer analysis section below.</p><p id="Par22">The software for the processing of the nucleotide sequences into coordinates and recovering the sequences from the coordinates is available here: <ext-link ext-link-type="uri" xlink:href="http://peteranoble.com/software.html">http://peteranoble.com/software.html</ext-link>. Details on the mathematics of iterative mapping of nucleotide sequences have been previously published [<xref ref-type="bibr" rid="CR19">19</xref>].</p></sec></sec><sec id="Sec7"><title>Mer analysis</title><p id="Par23">Mer analysis determines the presence/absence/frequency of a mer of length <italic>z</italic> (where <italic>z</italic> is 2 to 9) in a gene transcript.</p><sec id="Sec8"><title>Finding a specific mer in a transcript</title><p id="Par24">Let us assume that a database of the <italic>x-, y-</italic> coordinates of the target sequence has been assembled and we want to determine the presence/absence of the mer ‘AAACAA’ in a target sequence. There are three steps.</p><p id="Par25">First, we process the mer AAACAA into CGR space to find it <italic>x-, y-</italic> coordinates, which are − 0.734375 and 0.734375, respectively.</p><p id="Par26">Second, we determine the resolution of the search, which depends on mer length (i.e., resolution = 2<sup>(mer length)</sup>). A 6-mer requires a resolution of 64. The inverse of the resolution (1/resolution) is the CGR space around the coordinates that contain the specific mer. The CGR space around the coordinates is expressed by the following equation:</p><p id="Par27"><disp-formula id="Equa"><alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {x}^{\prime }=x\pm \frac{1}{r},{y}^{\prime }=y\pm \frac{1}{r}, $$\end{document}</tex-math><mml:math id="M2" display="block"><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi class="MJX-variant" mathvariant="normal">′</mml:mi></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mi>x</mml:mi><mml:mo>±</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>r</mml:mi></mml:mfrac><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi class="MJX-variant" mathvariant="normal">′</mml:mi></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mi>y</mml:mi><mml:mo>±</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>r</mml:mi></mml:mfrac><mml:mo>,</mml:mo></mml:math><graphic xlink:href="12864_2018_5042_Article_Equa.gif" position="anchor"></graphic></alternatives></disp-formula>where <italic>r</italic> is 2<sup>mer_length</sup>.</p><p id="Par28">For the 6-mer AAACAA<disp-formula id="Equb"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {x}^{\prime }=-0.734375\pm \frac{1}{64},{y}^{\prime }=0.734375\pm \frac{1}{64} $$\end{document}</tex-math><mml:math id="M4" display="block"><mml:msup><mml:mi>x</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mn>0.734375</mml:mn><mml:mo>±</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>64</mml:mn></mml:mfrac><mml:mo>,</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mn>0.734375</mml:mn><mml:mo>±</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>64</mml:mn></mml:mfrac></mml:math><graphic xlink:href="12864_2018_5042_Article_Equb.gif" position="anchor"></graphic></alternatives></disp-formula></p><p id="Par29">Third, the coordinates and CGR space of the mer is then used to search the CGR space of the target transcript sequence in the database. Any transcript that have coordinates within the box of x’ and y’ of the mer represents the sequence ‘AAACAA’. Furthermore, one can tally the number of hits within the box, which represents the frequency of the mer in a target sequence. We verified the presence of the mers in the identified target sequences by textual comparisons.</p></sec></sec><sec id="Sec9"><title>Statistical and bioinformatics analyses</title><p id="Par30">Analyses were conducted using SAS/JMP (version 7.0.2), R (version 3.4.0) and Microsoft Excel (versions 14.3.0 and 11.6). Hierarchical two-way cluster analysis was conducted on the binary matrices using Wards linkage method in SAS/JMP with default settings for cluster assignments. The resulting binary matrices were collapsed by their corresponding cluster assignments using a custom-designed program in C++. The resulting files were scaled to an average of zero and standard deviation of one in MS Excel and transferred to R to produce the heatmaps with no scaling. Network analysis was conducted using Gephi 0.9.2.</p></sec><sec id="Sec10"><title>Identification of 5’UTR, ORFs and 3’UTR and RNA motifs in transcripts</title><p id="Par31">RegRNA 2.0 was used to identify functional RNA motifs and sites in the gene transcripts [<xref ref-type="bibr" rid="CR20">20</xref>]. The server identifies splicing sites, splicing regulatory motifs, polyadenylation sites, transcriptional motifs, translational motifs, UTR motifs, mRNA degradation elements, RNA editing sites, riboswitches, RNA cis-regulatory elements, RNA-RNA interaction regions, and open reading frames using a integrated software package consisting of ~ 20 programs.</p><p id="Par32">Nucleotide sequences of the transcripts were individually submitted to the server, default search parameters specified, and tab-delimited results downloaded to a computer. The downloaded file contained global and local functions of the motifs and sites, their location in the transcript sequence, motif length and the sequence of the motif. Sequences of the unique mers in a transcript were matched to the sequence information of the motifs in the transcripts.</p></sec></sec><sec id="Sec11"><title>Results</title><p id="Par33">Our previous study on postmortem gene expression dynamics [<xref ref-type="bibr" rid="CR4">4</xref>] used a 60-mer oligonucleotide microarray to measure transcript levels. These perfectly matched probes were used in the present study to identify gene transcripts in the assembled datasets of the mouse and zebrafish (Additional file <xref rid="MOESM2" ref-type="media">2</xref>). A certain portion of the transcripts has been shown to significantly increase in abundance after organismal death relative to live controls [<xref ref-type="bibr" rid="CR4">4</xref>]. Henceforth, these transcripts are referred to as the over-abundant pool (OP), and transcripts not in this category are referred to as transcripts of the control pool (CP). Additional files <xref rid="MOESM3" ref-type="media">3</xref>, <xref rid="MOESM4" ref-type="media">4</xref>, <xref rid="MOESM5" ref-type="media">5</xref> and <xref rid="MOESM6" ref-type="media">6</xref> contain probes and their corresponding transcripts. In total, the OP of the mouse and zebrafish consisted of 333 and 230 gene transcripts, respectively, and the CP consisted of 32,611 and 27,433 transcripts, respectively.</p><p id="Par34">To determine if transcript length was a contributing factor when comparing different transcripts in the OP to those in the CP, we randomly selected two sets of transcripts from the CP (each set consisting of 333 gene transcripts for the mouse and 230 gene transcripts for the zebrafish) and compared the lengths of each set to those from the OP. No significant differences were found (two-tailed T-tests with unequal variance; alpha = 0.05) in either the mouse or the zebrafish, which rules out transcript length as a factor affecting Mer analyses (Additional file <xref rid="MOESM7" ref-type="media">7</xref>).</p><p id="Par35">Analysis of the composition of transcripts according to Esembl standards revealed that 86% of the transcripts in the OP of the mouse were protein coding, 9% were non-coding RNA, 3% were pseudogenes and 1% were unknown. In the OP of the zebrafish, 96% of the transcripts were protein coding, and 1% were long non-coding RNA, 3% were pseudogenes.</p><p id="Par36">Screening for sequence motifs (e.g., AU-rich elements (ARE), cis-regulatory elements of ERPIN, cis-regulatory elements of Rfam, exon splicing enhancer (ESE), functional RNA sequences, intron splicing enhancer (ISE), intron splicing silencer (ISS), long stems, microRNA target sites, ncRNA hybridization regions, polyadenylation sites (PAS), rho-independent terminator, transcriptional regulatory motif, untranslated region (UTR) motifs) was conducted in all of the OP transcripts of the mouse and 3 sets of random draws from the CP of the mouse. That is, each transcript sequence was submitted individually to RegRNA2 [<xref ref-type="bibr" rid="CR20">20</xref>] and the outputs were assembled into a dataset. Statistical analyses of the data using one-tailed T-tests revealed no significant differences (alpha = 0.05) in the number of instances by motif categories between the OP and CP. Of note, two categories approached significance, namely “functional RNA sequences” and “ncRNA hybridization regions” – however, the categories have ambiguous meanings.</p><p id="Par37">In summary, most of the transcripts showing postmortem abundance increase in the zebrafish and mice are protein coding.</p><sec id="Sec12"><title>Mer analyses</title><p id="Par38">The occurrences of 2- to 9-mers in gene transcripts of the OP were compared to those of the controls (i.e., CP). In the zebrafish, the controls consisted of 2- to 9- mers found in 30 sets of 230 transcripts that were randomly drawn (with replacement) from the CP of the zebrafish (Additional file <xref rid="MOESM8" ref-type="media">8</xref>). In the mouse, the controls consisted of 2- to 9- mers found in 30 sets of 333 gene transcripts that were randomly drawn (with replacement) from the CP of the mouse (Additional file <xref rid="MOESM9" ref-type="media">9</xref>).</p><p id="Par39">To test the assumption that the 30 sets of random draws sufficiently represented the diversity of transcripts found in each organism, we classified an additional 3 sets of 333 and 230 transcripts from the CPs of the mouse and zebrafish, respectively (without replacement) (Additional files <xref rid="MOESM10" ref-type="media">10</xref> and <xref rid="MOESM11" ref-type="media">11</xref>). Only transcripts not previously drawn were used in this test and the results are presented in the section below.</p><p id="Par40">The average count of individual mers from the random draws of the CPs were tabulated into a spreadsheet and compared to the counts of individual mers in the OPs of each organism. The arbitrary criterion used to identify ‘unique’ mers as either under- or over-represented was: a mer in the OP having less than or greater than 5 times the standard deviation of the average abundance of a corresponding mer in the CPs.</p></sec><sec id="Sec13"><title>Mer counts</title><p id="Par41">Given that 2-mers have 16 possible nucleotide combinations (i.e., AA, AT, AC, … TT) and 3-mers have 64 combinations (i.e., AAA, AAT… TTT), all short mers (2 to 3 nt) were anticipated to be present in transcripts of the OP and CPs, and therefore, no differences between the pools should be observed. Differences between the pools however, should change with increasing mer length presumably due to real differences or random chance (i.e., false-positives; FP).</p><p id="Par42">A maximum difference between the OP and CP pools was 6-mers (<italic>n</italic> = 74 transcripts) for the mouse and 5-mers (<italic>n</italic> = 18 transcripts) for the zebrafish (Table <xref rid="Tab1" ref-type="table">1</xref>, Fig. <xref rid="Fig2" ref-type="fig">2a</xref>). When normalized to the number of possible mer combinations, the maximum difference was 7-mers for the mouse and 5-mers for the zebrafish (Fig. <xref rid="Fig2" ref-type="fig">2c</xref>). Hence, mers of 5 to 7 nt in length are optimal for distinguishing between the pools.<table-wrap id="Tab1"><label>Table 1</label><caption><p>Average ± standard deviation of 333 gene transcripts in the mouse and 230 transcripts in the zebrafish that contained unique mers by group (OP vs. CP). The absolute difference in unique mer counts by group and mer length is shown</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Animal</th><th>Mer length</th><th>Num transcripts (OP)</th><th>Num transcripts (CP)</th><th>Absolute Difference</th></tr></thead><tbody><tr><td rowspan="8">Mouse</td><td>2</td><td>333 ± 0.25</td><td>333 ± 0.4</td><td>0 ± 0.1</td></tr><tr><td>3</td><td>330 ± 6.9</td><td>329 ± 7.9</td><td>2 ± 3.0</td></tr><tr><td>4</td><td>304 ± 30.3</td><td>304 ± 30.4</td><td>7 ± 7.2</td></tr><tr><td>5</td><td>227 ± 56.0</td><td>226 ± 55.9</td><td>13 ± 9.5</td></tr><tr><td>6</td><td>119 ± 52.5</td><td>117 ± 50.7</td><td>74 ± 45.2</td></tr><tr><td>7</td><td>43 ± 27.8</td><td>42 ± 25.5</td><td>6 ± 7.3</td></tr><tr><td>8</td><td>13 ± 10.8</td><td>12 ± 9.1</td><td>3 ± 4.2</td></tr><tr><td>9</td><td>3 ± 4.0</td><td>2 ± 1.9</td><td>3 ± 3.3</td></tr><tr><td rowspan="8">Zebrafish</td><td>2</td><td>230 ± 0.0</td><td>230 ± 0.1</td><td>0 ± 0.0</td></tr><tr><td>3</td><td>230 ± 0.6</td><td>229 ± 1.8</td><td>1 ± 1.5</td></tr><tr><td>4</td><td>220 ± 11.8</td><td>211 ± 15.5</td><td>9 ± 6.1</td></tr><tr><td>5</td><td>166 ± 35.7</td><td>148 ± 33.5</td><td>18 ± 8.1</td></tr><tr><td>6</td><td>81 ± 34.9</td><td>70 ± 29.4</td><td>12 ± 8.9</td></tr><tr><td>7</td><td>27 ± 17.5</td><td>23 ± 14.2</td><td>5 ± 5.2</td></tr><tr><td>8</td><td>8 ± 6.6</td><td>7 ± 5.0</td><td>2 ± 2.5</td></tr><tr><td>9</td><td>2 ± 2.4</td><td>2 ± 1.5</td><td>1 ± 1.2</td></tr></tbody></table></table-wrap><fig id="Fig2"><label>Fig. 2</label><caption><p>Mer counts as a function of mer length. Hatched line, mouse; solid line, zebrafish. Panel <bold>a</bold>, Difference in average mer counts by group (OP vs. CP); Panel <bold>b</bold>, individual mer counts that were 5 time stdev of average of the CP; Panel <bold>c</bold>, is the same results as panel B except normalized to the number of possible mer combinations and shown as a percentage; Panel <bold>d</bold>, number of mer counts that were 5 times stdev of the average CP due to random chance; average ± stdev of 3 random selections (without replacement)</p></caption><graphic xlink:href="12864_2018_5042_Fig2_HTML" id="MO2"></graphic></fig></p><p id="Par43">With increasing mer length, the number of ‘unique’ mers (i.e., over−/under-represented mers in the OP) increased (Fig. <xref rid="Fig2" ref-type="fig">2b</xref>).</p><p id="Par44">To determine the number of FPs as a function of mer length and test the integrity of the experimental design, we randomly draw three additional sets of transcripts from the CP (without replacement) and retained only transcripts not used in the previous analyses. For the mouse, each set consisted of 333 transcripts, and for the zebrafish, each set consisted of 230 transcripts. In this experiment, ‘over−/under- represented’ mers are FPs because the transcripts originated from the control transcript pool (i.e., the CP). To help explain the results of this experiment, let us consider the mer ATACCGG in the mouse. This mer would be considered ‘unique’ if its count were more or less than 5 times the standard deviation of the average from the CP, which is based on of 30 sets of 333 transcripts (Additional file <xref rid="MOESM11" ref-type="media">11</xref>, columns A to I). The average and standard deviation in the CP was 8 ± 3.5, meaning one would expect to find it an average of 8 times in random draws of 333 mouse transcripts. Five times the standard deviation is 17.5, therefore the range of critical values for the mer count is: − 9.5 and 25.5. In the OP, the mer occurred 31 times and is therefore considered ‘unique’ based on the stated criterion (i.e., the count is greater than 25.5).</p><p id="Par45">To further test the experimental design and check for FPs, mers were counted in three additional random draws from the CP. The 7-mer ATACCGG, for example, occurred in 7 of the 333 transcripts in one set, 3 of the 333 transcripts in the second set, and 10 of the 333 transcripts in the third set (Additional file <xref rid="MOESM11" ref-type="media">11</xref>, columns K to V). Since none of these counts are outside the criterion (the average ± standard deviation for this mer was 8.1 ± 3.49), there is no FPs for this mer. Of note, this procedure was repeated for all unique mers in the transcript pools of the mouse and zebrafish, respectively.</p><p id="Par46">The results show that the number of FPs in the OP was close to zero for mer lengths of up to 8 bp (compare Fig. <xref rid="Fig2" ref-type="fig">2d</xref> to <xref rid="Fig2" ref-type="fig">b</xref>). Therefore, while there is a possibility that some mers in the OP are FPs, the number was small (e.g., 8-mers: 1.0% are FPs in the mouse and 8.9% are FPs in the zebrafish).</p><p id="Par47">When the length of mers was 9, however, the number of FPs significantly increased to an average (± std) abundance of 1240 ± 167.2 for the mouse (31.3% FP) and 571 ± 158.8 for the zebrafish (34.2% FP).</p><p id="Par48">We also calculated the false discovery rate (FDR) using the formula: number of false discoveries divided by number of discoveries for each mer size. For the mouse, the FDR of: 6 mers or smaller was zero, 7 mers was 0.0036, 8 mers was 0.0102, and 9 mers was 0.3124. For the zebrafish, the FDR of: 7-mers or smaller was zero, 8 mers was 0.0842, and 9 mers was 0.3415. These findings are aligned with the FP results shown above.</p><p id="Par49">The results are consistent with the notion that unique mers can be identified in the OP by comparing them to random draws of mers from the CP. However, FPs and FDR increased with mer length. Taken together, over- and under- represented mers were identified in the OP and many are 5 to 7 nt in length. Since some of the under-represented mers would have been calculated to have negative abundances, the remainder of the study focused on the over-represented ones.</p><p id="Par50">The important finding in this section is that it is possible to find mers that are over- and under-represented in the OP compared to CP with statistical confidence. These mers are the subject of further discussion.</p></sec><sec id="Sec14"><title>Survey of the unique mers</title><p id="Par51">The survey of the OP identified 5117 unique mers in the mouse and 2245 mers in the zebrafish (Table <xref rid="Tab2" ref-type="table">2</xref>). Normalized to the total number of combinations of 3- to 9-mers (<italic>n</italic> = 349,504), this represents ~ 1.5% of the total mers in the mouse and ~ 0.6% in the zebrafish. Of note, 47 of the unique mers were common to both organisms (Table <xref rid="Tab3" ref-type="table">3</xref>).<table-wrap id="Tab2"><label>Table 2</label><caption><p>Number of unique mers in transcripts of the OP by mer length and organism</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Mer length</th><th>Zebrafish</th><th>Mouse</th></tr></thead><tbody><tr><td>4</td><td>1</td><td>0</td></tr><tr><td>5</td><td>31</td><td>1</td></tr><tr><td>6</td><td>67</td><td>62</td></tr><tr><td>7</td><td>118</td><td>279</td></tr><tr><td>8</td><td>356</td><td>819</td></tr><tr><td>9</td><td>1672</td><td>3956</td></tr><tr><td>Sum</td><td>2245</td><td>5117</td></tr></tbody></table></table-wrap><table-wrap id="Tab3"><label>Table 3</label><caption><p>Unique mers common to transcripts of the OP for the zebrafish and mouse</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Mer length</th><th>Mer</th></tr></thead><tbody><tr><td>6</td><td>AAAUAC, AACGAA, ACAUAA</td></tr><tr><td>7</td><td>UGUGAAC, AUCUCCA, AAAUACA, UAGGUUA, CAUGAAA</td></tr><tr><td>8</td><td>CAGAAAGC, GUAAAGUC, GCACAAAG, ACGAAUAC, AGAAGAGU, CAUGUGAA, AAAUACAU, AUAGGUUA</td></tr><tr><td>9</td><td>CCAAUGUGG, CUAUGAAGG, AAGUCCCAG, CUGACAGUC, UUCUCUGUG, GUUUCUGUG, CUAUGUCUG, AUACAAGUG, GCAAGGUUC, CAUGUGAAC, UCUAUGAAG, AUAGGUUAC, UCUGGGGCA, CCUGCUGCU, UAUCAUCGA, AAAAGAUCA, AUUCAAUGU, AAGAAAUCA, ACAAAAUCA, CUUCUCCAU, CAGAACCAU, UUUAACCAA, CAUGCAGAA, CUGGAAGAA, AUACAUCAA, AAAGAUCAU, CAGUAUGAA, AGAAAUCAU, CCUACGAAU, GUCCUGAAA, AACAUGAAA</td></tr></tbody></table></table-wrap></p><p id="Par52">In fact, some of these mers are reverse complements to one another, which is of interest because they might form secondary structures and play roles in post-transcriptional regulation (Table <xref rid="Tab4" ref-type="table">4</xref>). In the mouse, 218 of the 5117 mers (4.3%) reverse complemented one another. In the zebrafish, 31 of the 2245 mers (1.4%) were reverse complements.<table-wrap id="Tab4"><label>Table 4</label><caption><p>Unique mers that were reverse complements by length and organism</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Organism</th><th>Mer length</th><th>Mers</th></tr></thead><tbody><tr><td rowspan="5">Zebrafish</td><td>5</td><td>UGUAU, AUACA</td></tr><tr><td>6</td><td>GUAUUU, UCAAAA, UUUUGA, AAAUAC, AUUUUU, AUAUAU, AAAAAU</td></tr><tr><td>7</td><td>UGUAUUU, CGUUGUU, CAUUUUG, CAAAAUG, ACAAAAU, AACAACG, AUUGUAU, AUACAAU, AUUUUGU, AAAUACA</td></tr><tr><td>8</td><td>CAUUUUGA, UCAAAAUG</td></tr><tr><td>9</td><td>GGCGGCAAG, CCAGGCUCA, CGUCUAGGU, GCUAGGGAC, GUCCCUAGC, CUUGCCGCC, ACCUAGACG, UGAGCCUGG, AGUAGGCUA, UAGCCUACU</td></tr><tr><td rowspan="4">Mouse</td><td>6</td><td>CUAUAG, AUGCAU</td></tr><tr><td>7</td><td>ACCUAUA, AUGACUG, CAGUCAU, GUCUCUA, UAGAGAC, UAUAGGU, UCUAGAA, UUCUAGA</td></tr><tr><td>8</td><td>CCAUGACU, GGUUACAU, CCUAUAGG, GUCUCUAC, GUAGAGAC, CUUCUAGA, CUAGAAGU, CUAUGACU, GUAUGAAU, CUAUAGGU, ACCUAUAG, UCUGCAGA, AGUCAUGG, AGUCAUAG, UCUAGAAG, ACUUCUAG, UUCUAGAA, AUGUAACC, UUUGCAAA, AAUGCAUU, AAAGCUUU, AUUCAUAC</td></tr><tr><td>9</td><td>CGGAGAGAA, GCGAAGACA, CCCUUCUUC, GCUGCUGCU, CCUGGAACU, CCAGUGUGA, CCUGAGUUC, CCUCUUCUG, GGUCUUCAA, CCUUGAACU, CCAACAUCA, GGUUUCUCU, GGUUACAGU, GGUUACAUU, GAGGGCAUC, GUGGCUCAC, CAGGGAAGA, CUGCUGCUG, CAGCAGCAG, CUCCAGCAU, CUGCUCUCU, CUGCAGAAG, CAGGAGAAA, CUCCUUCCU, CAGGAAGCA, CAGGAAGGA, GUGAGCCAC, GUCUCCUGU, GAGUGGUAG, GUCUCCAAA, CACAGAGAA, CUGAGUUCA, CAGAGAAAA, GUCUUCGCU, CAGAAGAGG, CUGAAGACA, GUCUUCAGA, CUGAAGAUG, CAGAAGAUG, CAGAAAGCA, CAGUAUGAA, GAUGCCCUC, CUUCCCAUC, GAUGGGAAG, CUACCACUC, CUUCCUCUU, GAACCUUUU, CUUCUGCAG, CUUGAGGAA, GAACUCAGG, GAACACACA, GAAGACACA, CUUCACUUG, CAAGUGAAG, GAAGAAGGG, CAUCUUCUG, CAUCUUCAG, CUUCUAGAA, GAAGAUGAU, GAAGAAGAA, CAAAGCCUU, CAAAGACUU, CAAACUUCU, GUUACAUUU, GUAAAGACU, CAAAUGUAA, CUUUUAAAA, UCCCAGCAA, UGGGAAGGA, AGCGAAGAC, UGCUGGGAA, AGCAGCAGC, UGCUGCUGU, UCCUGCAAA, UGCAGAAGA, UCCUUCCCA, UCCUUCCUG, UGCUUCCUG, AGGAAGGAG, UGCUUUCUG, ACAGCAGCA, ACAGGAGAC, ACACCAACA, AGAGAGCAG, UCACACUGG, UCUCUGUGU, ACACAGAGA, UGUGUCUUC, UGUGUGUUC, ACAGAGAAA, UGUCUUCGC, UCUGAAGAC, UGUCUUCAG, ACUGUAACC, AGAGAAACC, ACUGUUUCU, ACACAUACA, AGUCAUAGU, AGAGUUUCU, AGUCUUUAC, UCUUCCCUG, AGUUCCAGG, UGUUGGUGU, UCUUCUGCA, UGAACUCAG, UCUUCACAA, UGAUGUGUU, UGAUGUUGG, AGUUCAAGG, UCUAGAAGU, UCAUCUUCU, UCAUCUUCA, UGUUCUUCA, UGAAGUUCU, UGAAGAUGA, UGAAGAACA, ACUUCUAGA, AGUUGUUCU, AGAACUUCA, AGAACAACU, AGAAGAUGA, AGAAGUUUG, UGUAUGUGU, ACUAUGACU, AGAAACUCU, AGAAACAGU, AGUAUGAAU, UCAAAAACA, UGUUUUUGA, UGUUAUAAA, UUCCCAGCA, AAGGCUUUG, UUGCUGGGA, AUGCUGGAG, UUCCUCAAG, UUGCAGAAA, AAGCAGUUA, UUCCUUCUU, AAGAGGAAG, UUCUCUCCG, UUCUCUGUG, UUGUGAAGA, AACACAUCA, AAGUCUUUG, AAGAAGGAA, UUGAAGACC, UUCUUCUUC, UUCUAGAAG, AUCAUCUUC, AUCAUGAAA, UUCAUACUG, UUGUUUUUU, AAGUUAUUA, AAAGCCUUU, AAAGGCUUU, UUUGCAGGA, UUUGGAGAC, AAAGCUUUU, UUUCUCCUG, UUUCUGCAA, UAACUGCUU, UUUCUCUGU, AAAGUCUUU, AAAGACUUU, UUUCAUGAU, AAUGUAACC, AUUCAUACU, UUACAUUUG, UUUCAAAAA, UAUGUAUUU, AAAAGGUUC, AAAAGCUUU, UUUUCUCUG, AAAUGUAAC, UUUUCAAAA, UUUUGUUUU, UUUUGAAAA, AAAACAAAA, UUUUUGAAA, AAAUACAUA, UAAUAACUU, AAAAAACAA, UUUAUAACA, UUUUAAAAG, UUUUUUUAA, UUAAAAAAA</td></tr></tbody></table></table-wrap></p></sec><sec id="Sec15"><title>Number of unique mers per transcript</title><p id="Par53">The distribution of the unique mers was investigated to determine if they were found in all transcripts of the OP, or just a few. In other words, <italic>is the distribution of unique mers uniform across all transcripts?</italic> To address this question, we compared their distributions in both transcript pools (i.e., OP and CP). Here we assumed that the corresponding unique mers in transcripts of the CP should approximate a skewed (Poisson) distribution because they are relatively rare occurrences. The controls in this experiment were the three sets of random draws (with replacement) from the CP. We also examined the multiple occurrences of unique mers in the OP since a unique mer might occur multiple times in the same transcript.</p><p id="Par54">In the zebrafish, the frequencies of the unique mers per transcript varied between pools (Fig. <xref rid="Fig3" ref-type="fig">3</xref>). These findings indicate that not all transcripts in the OP have the same number of unique mers – i.e., the number of unique mers in a transcript was not uniform. In the OP, the maximum bin was 150 while the maximum bin in the CP was 100. Some transcripts of the OP have more than twice the number of unique mers in the 200, 250, and 300+ bins than those of the CP. Therefore, some zebrafish transcripts in the OP have many more unique mers than others.<fig id="Fig3"><label>Fig. 3</label><caption><p>Distribution of unique mers per gene transcript in the zebrafish. <bold>a</bold>, unique mers in OP; <bold>b</bold>, multiple unique mers in OP; <bold>c</bold>, unique mers in CP (3 independent random selections; each as a different shade of grey); <bold>d</bold>, multiple unique mers in CP (3 independent random selections)</p></caption><graphic xlink:href="12864_2018_5042_Fig3_HTML" id="MO3"></graphic></fig></p><p id="Par55">In terms of multiple occurrences of unique mers in the zebrafish, the distributions differed by pool also, with multiple unique mers occurring within the same transcript when compared to controls (Fig. <xref rid="Fig3" ref-type="fig">3b</xref>). For example, about 87 of the OP transcripts had more than 300 multiple unique mers compared to about 40 in the CP (Fig. <xref rid="Fig3" ref-type="fig">3d</xref>). Hence, not only are there many more unique mers in the OP but, in some cases, there are multiple occurrences of the same mer in the same transcript.</p><p id="Par56">In the mouse, the frequency distribution of unique mers per transcript was also different between the pools (Fig. <xref rid="Fig4" ref-type="fig">4a</xref>). Specifically, there was almost double the number of unique mers in the 200 bin of the CP than the OP, about the same number of unique mers in the 400 bin, and twice (or more) the number in the 600, 1000, and 1200+ bins of the OP than the CP (Fig. <xref rid="Fig4" ref-type="fig">4c</xref>). This finding is consistent with those of the zebrafish – i.e., there are many more unique mers in the OP than the CP.<fig id="Fig4"><label>Fig. 4</label><caption><p>Distribution of unique mers per gene transcript in the mouse. <bold>a</bold>, unique mers in OP of the mouse; <bold>b</bold>, multiple unique mers in OP of the mouse; <bold>c</bold>, unique mers in CP of mouse (3 independent random selections; each displayed as a different shade of grey); <bold>d</bold>, multiple unique mers in CP of the mouse (3 independent random selections)</p></caption><graphic xlink:href="12864_2018_5042_Fig4_HTML" id="MO4"></graphic></fig></p><p id="Par57">In terms of the multiple mer occurrences in the mouse, the results were different from the zebrafish; in general, there was little change between the histogram of the unique and multiple mers (compare Fig. <xref rid="Fig4" ref-type="fig">4a to b</xref>) – meaning that in contrast to the zebrafish, most of the unique mers did not occur multiple times in the same transcript sequence. Of note, this was not true for all cases as the 1200+ bin was somewhat bigger in the Fig. <xref rid="Fig4" ref-type="fig">4b</xref> than A. However, when compared to Fig. <xref rid="Fig3" ref-type="fig">3b</xref> to a, there is a substantial difference between unique and multiple mers in the zebrafish. The presumed reason for this disparity is that in the mouse, the unique mers tend to be longer in length than those in the zebrafish (Table <xref rid="Tab1" ref-type="table">1</xref>, Fig. <xref rid="Fig2" ref-type="fig">2c</xref>) and the longer the length, the less frequent its occurrence.</p><p id="Par58">Given the peculiar expression pattern of OPs (i.e., increase in abundance postmortem) it is perhaps quite indicative that the distribution of unique mers in the OP differs from those in the CP. Furthermore, there appears to be differences in multiple unique mers of these transcripts in the zebrafish but less so in the mouse.</p></sec><sec id="Sec16"><title>Groups of unique mers in the OP transcripts</title><p id="Par59">Based on the previous analyses, we rationalized that some transcripts in the OP might share the same unique mers. To investigate the relationships among the OP transcripts and the unique mers (in binary presence/absence format), we constructed matrices and then performed two-way hierarchical clustering. The matrix for the zebrafish consisted of 230 rows of transcripts by 2245 columns of unique mers (Additional file <xref rid="MOESM12" ref-type="media">12</xref>), and the matrix for the mouse consisted of 333 rows of transcripts by 5117 columns of unique mers (Additional file <xref rid="MOESM13" ref-type="media">13</xref>).</p><p id="Par60">The cluster analysis of the zebrafish identified 14 groups of transcripts and 20 groups of mers with high similarities, and the analysis of the mouse yielded 16 groups of transcripts and 20 groups of mers. The groups were collapsed by summation. For example, group A of the transcripts in the zebrafish consisted of 36 transcripts and Set 1 of the mers consisted of 25 unique mers. In total, 25 × 36 = 900 combinations, out of which 119 were actual occurrences of mers in the said transcripts (Additional file <xref rid="MOESM14" ref-type="media">14</xref>), meaning there were 119 occurrences in the collapsed group. We summed groups A to N and mer sets 1 to 20 to form a collapsed matrix of 14 columns of transcript groups by 20 rows of mer sets. The same procedure was repeated for the mouse. The collapsed groups were normalized by row (see Materials and Methods section) to produce the data for the heat maps. Note, the heat maps were turned 90 degrees to show transcripts as columns and mer sets as rows.</p><p id="Par61">The number of transcripts in a group and the number of mers in a set varied substantially for both organisms. Specifically, in the zebrafish, the number of transcripts in a group ranged from 1 to 59 (of the 230) (Fig. <xref rid="Fig5" ref-type="fig">5</xref>), and in the mouse, the number of transcripts by group ranged from 1 to 124 (of the total of 333) (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). Hence, some transcripts are very similar to one another in terms of unique mers, while others are distinctly different – there was no uniformity (i.e., equal number of mers distributed to equal number of transcripts).<fig id="Fig5"><label>Fig. 5</label><caption><p>Heatmap of transcript groups and mer sets for the zebrafish. M, count of mers in group; N, count of transcripts in group. White, high count; yellow-orange, median count; red, low count</p></caption><graphic xlink:href="12864_2018_5042_Fig5_HTML" id="MO5"></graphic></fig><fig id="Fig6"><label>Fig. 6</label><caption><p>Heatmap of Transcript groups and mer sets for the mouse. M, count of mers in group; N, count of transcripts in group. White, high count; yellow-orange, median count; red, low count</p></caption><graphic xlink:href="12864_2018_5042_Fig6_HTML" id="MO6"></graphic></fig></p><p id="Par62">The number of unique mers in a set ranged from 2 to 876 (of a total of 2245) in the zebrafish (Fig. <xref rid="Fig5" ref-type="fig">5</xref>) and from 40 to 1407 (of a total of 5117) in the mouse (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). Hence, some groups of mers are found in the same transcripts while others are found in different ones. Similar to the situation with the transcripts, the relationship among the mers was not straightforward– there appears to be a pattern.</p><p id="Par63">There are unifying features visible in the heatmaps. For example, all transcript groups in the zebrafish contained relatively similar counts of mers within the mer sets 5 as well as 19 (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). Similarly, in the mouse, all transcript groups had similar counts of mers within the mer set 2 (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). Hence, despite similarities and differences of the collapsed data, there are common sets of mers found within all transcripts.</p><sec id="Sec17"><title>Zebrafish heatmap</title><p id="Par64">In terms of differences, groups J and N are dissimilar from the other transcript groups (Fig. <xref rid="Fig5" ref-type="fig">5</xref>) and each group consists of a single transcript. Group J represents the transcript <italic>si_ch211-69b7.6</italic>, whose function is currently not known, and Group N represents the transcript <italic>Psd2</italic> (Pleckstrin and Sec7 domain containing 2), which is involved in regulating vesicle biogenesis in intracellular trafficking. The groups differed from the other groups in terms of the counts of mer sets 2, 3 and 8, which contain 876, 111, and 111 mers, respectively.</p><p id="Par65">There appears to be significant differences between transcript group D, G, C, L, K and M, which consist of 122 transcripts (of the 230 possible) and group A, B, H, I, E and F, which consist of 106 transcripts (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). These groups are distinct due to subtle differences in mer set 10, which consists of 9 mers and mer set 6, which consists of 60 mers.</p></sec><sec id="Sec18"><title>Mouse heatmap</title><p id="Par66">The heatmap of the mouse shows similar variation in the number of transcripts by group and mer set (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). Transcript group K, B, E, M, H, F, C and L, which represent 40 transcripts (of a possible 333) is different from group O, N, J, A, G, P, D, and I (293 of possible 333 transcripts). The mer set responsible for this difference is mer set 17, which contains 1407 mers. Group O, N, J, A, G, P, D, and I have a higher relative counts than group K, B, E, M, H, F, C and L. Interestingly, the 40 transcripts in group K, B, E, M, H, F, C and L are annotated as either zinc finger proteins or predicted coding genes – and not one of the transcripts encode a protein with known function.</p><p id="Par67">In the zebrafish, most (192) of the known functional gene transcripts are dispersed into many groups N, C, L, K, M, A, B, H, I and F, which represent 83% of the OP (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). In the mouse, most (245) of the known functional gene transcripts are found in groups A and G, which represent 74% of the OP (Fig. <xref rid="Fig6" ref-type="fig">6</xref>).</p><p id="Par68">Without stating any biological significance, there appears to be underlying patterns in the occurrence of unique mers in transcripts of the OP. At this stage of the analysis, it is difficult to quantify the patterns other than visually.</p></sec></sec><sec id="Sec19"><title>Density of multiple mers by transcript and organism</title><p id="Par69">We examined the number of ‘unique’ mers by transcript length since longer transcripts might have more mers (Additional file <xref rid="MOESM15" ref-type="media">15</xref>). Indeed, this was found true for the zebrafish -- there were more ‘unique’ mers with increasing transcript length (Pearson correlation coefficient, <italic>r</italic> = 0.55, <italic>p</italic> < 0.001). However, this relationship did not hold for the mouse (and we will show why below).</p><p id="Par70">The averaged (± stdev) density of multiple mers for the zebrafish was 0.14 ± 0.18 mers/nt (<italic>n</italic> = 230) and for the mouse was 0.40 ± 0.67 mers/nt (<italic>n</italic> = 333). That is, there are 14 unique mers for every 100 nucleotides in the transcripts of the zebrafish and 40 mers for every 100 nucleotides in the transcripts of the mouse. Note the high standard deviations indicating a wide variation in values.</p><p id="Par71">The highest and lowest densities of unique mers also differed between organisms. In the zebrafish, the highest density was ~ 1.0 mers/nt for <italic>Pimr</italic> gene transcripts, which corresponds to clusters B and H (Fig. <xref rid="Fig5" ref-type="fig">5</xref>), and the lowest density was ~ 0.04 mers/nt for transcripts found in cluster A. We plotted the relationship between multiple mers and transcript length to find that the <italic>Pimr</italic> gene transcripts are distinctly different (red dots) from those in the rest of the transcripts in the OP (black dots) (Fig. <xref rid="Fig7" ref-type="fig">7a</xref>). The <italic>Pimr</italic> genes encode proto-oncogene serine/threonine-protein kinases involved in regulating the cell cycle. The remaining transcripts have a linear relationship between multiple mers and transcript length (<italic>y</italic> = 0.1<italic>×</italic>; R<sup>2</sup> = 0.91; with <italic>x</italic> is transcript length and <italic>y</italic> is multiple mers).<fig id="Fig7"><label>Fig. 7</label><caption><p>Number of multiple unique mers in transcripts versus transcript length. <bold>a</bold>, zebrafish; <bold>b</bold>, mouse; Red, deviant transcripts. Red dots in the zebrafish correspond to Pimr transcripts; Red dots in the mouse represent 47 transcripts (see text)</p></caption><graphic xlink:href="12864_2018_5042_Fig7_HTML" id="MO7"></graphic></fig></p><p id="Par72">In the mouse, the highest density was ~ 2.6 mers/nt for annotated transcripts that do not have a canonical name (e.g., Gm14410, Gm14305. Gm14434, Gm2026, Gm11007, Gm2007, Gm4631) and were associated with Cluster B (Fig. <xref rid="Fig6" ref-type="fig">6</xref>) and the lowest density was ~ 0.04 mers/nt in transcripts associated with cluster A. A plot of the multiple mers by transcript length for the mouse revealed significant differences for a subset of the transcripts (red dots) when compared to the rest (black dots) (Fig. <xref rid="Fig7" ref-type="fig">7b</xref>). The red dots represent 47 annotated gene transcripts, many that do not have a canonical name and includes those with the highest mer densities per transcript (mentioned above). The red dots also include 25 transcripts annotated as zinc finger proteins, 3 Rik transcripts, 1 unprocessed pseudogene, 1 Fam containing transcript, and 10 functional gene transcripts. The remaining transcripts have a linear relationship between multiple mers and transcript length (<italic>y</italic> = 0.1<italic>×</italic>; R<sup>2</sup> = 0.95; with <italic>x</italic> is transcript length and <italic>y</italic> is multiple mers). Hence the reason for the poor correlation between multiple mers and transcript length in the mouse data (noted above) was due to 47 transcripts that deviated from the other 286 transcripts in terms of their mer density.</p><p id="Par73">We used <bold>RNAReg2</bold> to determine if there are any unique molecular features in the 10 functional gene transcripts: <italic>Bpifc, Ifitm7, Ms4a4c, Platr25, Rex2, Spag7, Styk1, Sva, Tmem239, Tnfrsf9</italic>. We specifically looked at the relationship between the unique mers in the transcripts and the tab-delimited output files from RegRNA2 (Additional file <xref rid="MOESM15" ref-type="media">15</xref>).</p><p id="Par74">While all of the transcripts have ‘ncRNA hybridization regions’ that matched the unique mers, no patterns could be found in the AU-rich elements, K-boxes, UNR boxes, untranslated region motifs, long stem loop structures or transcriptional regulatory motifs among the 10 functional genes. Therefore, we concluded that the gene transcripts contain putative ncRNA hybridization regions – but we have no supporting evidence that these regions are actually used by the transcriptional regulation.</p><p id="Par75">We rationalized that the transcripts with high mer densities might act as molecular sponges to RBPs and ncRNAs and thus alter their availability in the intracellular pools. If so, one would expect the profiles (i.e., transcript abundance by postmortem time) of transcripts with high densities and those transcripts affected by them to be highly correlated. Moreover, they should share similar unique mers that serve as putative binding sites. Principal component analysis was used to find patterns among transcripts with high mer densities using the correlations of their transcript abundance profiles to the rest of the profiles in the OP of the mouse brain. Network analysis was used to find shared mer binding sites.</p><p id="Par76">The two axes of the ordination plot accounted for 96% of the variability (Fig. <xref rid="Fig8" ref-type="fig">8a</xref>). There appears to be three distinct areas in the ordination plot. One location is occupied by Gm14399, the other location is populated by a group of 8 gene transcript and the third location is occupied by Gm14409. The correlations among the transcript profiles differed by high density transcripts suggesting that certain groups might regulate different sets of transcripts.<fig id="Fig8"><label>Fig. 8</label><caption><p>Ordination plot of transcripts with high mer densities (<bold>a</bold>) and network of transcripts with shared mers (<bold>b</bold>). The ordination was based on the correlations among mouse brain transcript profiles. The network was based on the number of shared mers in subset of the transcript profiles with high R<sup>2</sup> (> 0.95) to the transcripts with high mer densities. The network shows that the transcripts with high mer densities (i.e., molecular sponges) shared mers with many other transcripts</p></caption><graphic xlink:href="12864_2018_5042_Fig8_HTML" id="MO8"></graphic></fig></p><p id="Par77">To investigate the connections within the networks, we took a subset of the transcripts with high R<sup>2</sup>s (> 0.95), and counted the number of shared mers. A network plot revealed that transcripts with high mer densities are connected to many different transcripts and that some shared similar mers. For example, Gm14305 shared mers with Gm11007, Gm2007, Gm14308 and <italic>Hhmt1</italic> as wells as many other transcripts (Fig. <xref rid="Fig8" ref-type="fig">8b</xref>). This finding suggests that the number of possible transcripts (and pathways) that are affected by molecular sponges appears to be quite vast.</p><p id="Par78">The critical finding of the above analysis is that the transcripts with high mer densities have the potential to act as molecular sponges to other transcripts and thus regulate them post-transcriptionally.</p></sec><sec id="Sec20"><title>Multiple mer density by region (5’UTR, ORFs, 3’UTR)</title><p id="Par79">To investigate the density of unique mers by region, up to ten transcripts from each cluster (Figs. <xref rid="Fig5" ref-type="fig">5</xref> and <xref rid="Fig6" ref-type="fig">6</xref>) were compared to determine if there are significant differences in mer density by region (Additional file <xref rid="MOESM16" ref-type="media">16</xref>). Note that not all transcripts had 5’UTR and/or 3’UTR regions and some lacked ORFs (e.g., ncRNA).</p><p id="Par80">In the zebrafish, for the transcripts having all three regions, the 3’UTR region had significantly more mers/nt than the other two regions (Table <xref rid="Tab5" ref-type="table">5</xref>, Paired two-tailed T-tests, <italic>p</italic> < 0.0001). Transcripts lacking 5’UTR, 3’UTR, or ORFs have low densities (i.e., ~ 0.1 mers/nt), indicating regional effects.<table-wrap id="Tab5"><label>Table 5</label><caption><p>Number of unique mers by nucleotide (transcript length), region and organism. Two-way paired t-test across rows: a,b, <italic>p</italic> < 0.0001; c,d <italic>p</italic> < 0.01</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Organism</th><th>Regions</th><th>Number of transcripts</th><th>5’UTR</th><th>ORF</th><th>3’UTR</th><th>Non-coding</th></tr></thead><tbody><tr><td rowspan="5">Zebrafish</td><td>5’UTR, ORF, 3’UTR</td><td>70</td><td>0.2 ± 0.26<sup>a</sup></td><td>0.2 ± 0.28<sup>a</sup></td><td>0.3 ± 0.35<sup>b</sup></td><td>–</td></tr><tr><td>ORF, 3’UTR</td><td>4</td><td>–</td><td>0.1 ± 0.01</td><td>0.1 ± 0.03</td><td>–</td></tr><tr><td>5’UTR, ORF</td><td>1</td><td>0.1</td><td>0.1</td><td></td><td>–</td></tr><tr><td>ORF</td><td>3</td><td>–</td><td>0.1 ± 0.03</td><td>–</td><td>–</td></tr><tr><td>Non-coding</td><td>1</td><td>–</td><td>–</td><td>–</td><td>0.1</td></tr><tr><td rowspan="5">Mouse</td><td>5’UTR, ORF, 3’UTR</td><td>65</td><td>0.4 ± 0.60<sup>a</sup></td><td>0.6 ± 0.70<sup>b</sup></td><td>0.5 ± 0.80</td><td>–</td></tr><tr><td>ORF, 3’UTR</td><td>2</td><td>–</td><td>0.1 ± 0.00</td><td>0.1 ± 0.00</td><td>–</td></tr><tr><td>5’UTR, ORF</td><td>16</td><td>1.4 ± 0.80<sup>c</sup></td><td>2.0 ± 0.70<sup>d</sup></td><td>–</td><td>–</td></tr><tr><td>ORF</td><td>11</td><td>–</td><td>2.3 ± 0.50</td><td>–</td><td>–</td></tr><tr><td>Non-coding</td><td>10</td><td>–</td><td>–</td><td>–</td><td>0.4 ± 0.50</td></tr></tbody></table></table-wrap></p><p id="Par81">In contrast, the highest unique mer densities in the mouse were found in the ORFs of transcripts – not the 3’UTR region as in the zebrafish (Table <xref rid="Tab5" ref-type="table">5</xref>). In transcripts having all three regions, the ORFs had significantly higher densities than the 5’UTR (Paired two-tailed T-test, <italic>p</italic> < 0.0001). In gene transcripts that have both 5’UTR and ORFs (no 3’UTR), or those having neither 5’UTR nor 3’UTR regions (i.e., ORF only) had twice the mer densities than transcripts having all three regions. Moreover, higher mer densities were found in the ORFs than the 5’UTR (Table <xref rid="Tab5" ref-type="table">5</xref>, Paired two-tailed T-tests, <italic>p</italic> < 0.01). One possible reason for these differences is that the 16 samples having no 3’UTR and the 11 samples lacking untranslated regions (i.e., they were all ORFs) consist of genes annotated as ‘predicted coding gene’ or ‘zinc finger protein gene’. Hence, gene function might play a role in these differences.</p><p id="Par82">In summary, the highest mer densities in the zebrafish were found in the 3’UTR while the highest densities in the mouse were found in the ORFs. Perhaps this is not a surprise given the evolutionary distance between zebrafish and mice.</p></sec><sec id="Sec21"><title>Known motifs</title><p id="Par83">The following motifs are associated with increased mRNA stability or gene expression: the <italic>Hud</italic> binding site, YUNNYUY [<xref ref-type="bibr" rid="CR21">21</xref>]; the <italic>Rbfox</italic> binding site, UGCAUG [<xref ref-type="bibr" rid="CR10">10</xref>]; and UAUUUAU, GAGAAAA, AGAGAAA, UUUGCAC, AUGUGAA, UUGCACA, GGGAAGA [<xref ref-type="bibr" rid="CR22">22</xref>]. We screened these motifs against the unique mers to identify transcripts in the OP that might have increased stability or gene expression due to these motifs.</p><p id="Par84">Three hundred and fourteen of the 333 OP transcripts (94%) in the mouse and 189 of the 230 transcripts (82%) in the zebrafish contained one or more of the known binding motifs associated with increased mRNA stability or gene expression (Table <xref rid="Tab6" ref-type="table">6</xref>). Most of the transcripts in the OP of the mouse and zebrafish had at least two different motifs (Fig. <xref rid="Fig9" ref-type="fig">9</xref>). The number of previously reported motifs represents a small fraction of the total number of unique mers found in our study (180 of the 5117 unique mouse mers (3.5%) and 54 of the 2245 zebrafish mers (2.4%)). Hence, our study identified 4937 and 2191 putatively new motifs in transcripts of the OP of the mouse and zebrafish, respectively. It remains to be determined if these new motifs are functional or not.<table-wrap id="Tab6"><label>Table 6</label><caption><p>Number of transcripts by known protein binding site and organism. Hud binding site, YUNNYUY [<xref ref-type="bibr" rid="CR21">21</xref>]; Rbfox binding site, UGCAUG [<xref ref-type="bibr" rid="CR10">10</xref>]; and UAUUUAU, GAGAAAA, AGAGAAA, UUUGCAC, AUGUGAA, UUGCACA, GGGAAGA [<xref ref-type="bibr" rid="CR22">22</xref>]</p></caption><table frame="hsides" rules="groups"><thead><tr><th rowspan="2">Organism</th><th rowspan="2"><italic>n</italic> transcripts</th><th colspan="4">Protein binding sites</th></tr><tr><th><italic>Hud</italic></th><th><italic>Rbfox</italic></th><th>Jacobsen et al.</th><th>All three</th></tr></thead><tbody><tr><td>Mouse</td><td char="." align="char">333</td><td char="." align="char">287</td><td>126</td><td>258</td><td>314</td></tr><tr><td>Zebrafish</td><td char="." align="char">230</td><td char="." align="char">185</td><td>0</td><td>106</td><td>189</td></tr></tbody></table></table-wrap><fig id="Fig9"><label>Fig. 9</label><caption><p>Number of known binding sites per transcripts for the mouse (<bold>a</bold>) and zebrafish (<bold>b</bold>). Total number of transcripts for the mouse, <italic>n</italic> = 333 and for the zebrafish, <italic>n</italic> = 230. The following binding sites were examined: Hud binding site, YUNNYUY [<xref ref-type="bibr" rid="CR21">21</xref>]; Rbfox binding site, UGCAUG [<xref ref-type="bibr" rid="CR10">10</xref>]; and UAUUUAU, GAGAAAA, AGAGAAA, UUUGCAC, AUGUGAA, UUGCACA, GGGAAGA [<xref ref-type="bibr" rid="CR22">22</xref>]. Note: the zebrafish did not have Rbfox binding sites</p></caption><graphic xlink:href="12864_2018_5042_Fig9_HTML" id="MO9"></graphic></fig></p><p id="Par85">Finding the known motifs in the OP transcripts gives the credence to the notion that the OP transcripts are indeed undergoing some form of postmortem “regulation” or stabilization.</p></sec></sec><sec id="Sec22"><title>Discussion</title><p id="Par86">The motivation for our study was driven by curiosity into possible mechanisms responsible for the increase in transcript abundances with postmortem time, which have now been reported to occur in the zebrafish, mouse, and humans [<xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR4">4</xref>]. There is a need to understand regulatory features and how they influence transcriptional dynamics in order to comprehend the response of biological systems to stress. Yet, to our knowledge, no study has investigated possible reasons for increases in transcript abundance after organismal death. Such information is needed to provide baseline data for gene expression studies involving stressful conditions such as disease, starvation, and cancer.</p><sec id="Sec23"><title>Unique mers identified in the OP</title><p id="Par87">Our initial hypothesis was that among multiple reasons, there must be a signal, i.e., a nucleotide sequence that is responsible for postmortem activation of certain transcripts. Instead, we find sets of ‘unique’ mers in different groups of transcripts, with most sets consisting of ten to hundreds of different mers -- not just one or two.</p><p id="Par88">The total number of unique mers in the OP was relatively small compared to all possible mers, ~ 1.5% of the total combinations of 3- to 9- mers in the mouse and ~ 0.6% in the zebrafish. These small percentages are presumably due to the arbitrary criterion used to identify unique mers. The reason the criterion was set to 5 times the standard deviation of the average count of the mer in the CP was to ensure that the identified mers were not due to random chance (i.e., false positives, FPs). Our results indicate that chance of a random mer having a count exceeding the criterion was relatively rare -- but FPs did occur and their occurrence increased with mer length (Fig. <xref rid="Fig2" ref-type="fig">2d</xref>).</p><p id="Par89">The fact that several mers identified in our study have been previously reported to be involved with increased gene expression and/or mRNA stability (e.g., <italic>Hud</italic>, <italic>Rbfox, ARE</italic> binding sites; [<xref ref-type="bibr" rid="CR10">10</xref>, <xref ref-type="bibr" rid="CR21">21</xref>, <xref ref-type="bibr" rid="CR23">23</xref>]) is consistent with the idea that our experimental design was effective at identifying ‘unique’ mers in the postmortem transcriptome of two different organisms.</p></sec><sec id="Sec24"><title>Unique mers by transcript, region, and organism</title><p id="Par90">The number of unique mers in each transcript of the OP varied considerably. Some transcripts have a disproportionately high number of mers, while others have much lower numbers. Interestingly, in the mouse, several of the transcripts with high multiple mer densities have an ORF with no known function. Other transcripts have known functions, including: <italic>Bpifc</italic>, which is involved in innate immune response; <italic>Fam160b2,</italic> which is involved in phosphorylation of <italic>Hsp70</italic> [<xref ref-type="bibr" rid="CR24">24</xref>]; <italic>Ifitm7,</italic> which is involved in regulation of cell proliferation and immune response [<xref ref-type="bibr" rid="CR25">25</xref>]; <italic>Ms4a4c,</italic> which regulates receptor signaling and recycling [<xref ref-type="bibr" rid="CR26">26</xref>]; <italic>Spag7,</italic> which is involved in antiviral and inflammatory response [<xref ref-type="bibr" rid="CR27">27</xref>–<xref ref-type="bibr" rid="CR29">29</xref>]; <italic>Styk1,</italic> which is associated with cancer progression and promotes the Warburg effect through signaling of the PI3K/AKT pathway [<xref ref-type="bibr" rid="CR30">30</xref>–<xref ref-type="bibr" rid="CR32">32</xref>]; and <italic>Tnfrsf9,</italic> which is involved in positive regulation of immune system functions and leukocyte activation [<xref ref-type="bibr" rid="CR33">33</xref>]. In the zebrafish, a disproportionately high number of mers occurred in the <italic>Pimr</italic> gene transcripts, which are involved in cell cycling. These gene transcripts have common functions: cell survival, proliferation, cycling, stress compensation, and/or defense. It is enticing to speculate that the other transcripts (i.e., those with no known functions but with high mer densities) might also be involved in these functions.</p><p id="Par91">The density of multiple unique mers was higher in the ORFs than the 3’UTR in the mouse -- but quite the opposite was true in the zebrafish (Table <xref rid="Tab5" ref-type="table">5</xref>). That is, the zebrafish had a higher mer density in the 3’UTR than the other regions. In general, the 3′ UTR is involved in subcellular localization and mRNA stability, while the 5′ UTR play roles in translational control [<xref ref-type="bibr" rid="CR34">34</xref>]. Motifs within the UTR regions are thought to control functions by interacting with RBPs [<xref ref-type="bibr" rid="CR34">34</xref>]. Yet, the highest density of mers (2.3 ± 0.50 mers/nt) was in 11 transcripts that lacked UTRs (i.e., they were all ORFs). These findings are aligned with the notion that binding sites can exist all along the transcripts and not necessarily restricted to the UTRs [<xref ref-type="bibr" rid="CR35">35</xref>]. It is possible that these 11 transcripts act as large “molecular sponges” in stressful conditions, providing an additional layer of complexity to post-transcriptional regulation (which we discuss below).</p><p id="Par92">While the two organisms share 47 unique mers, there were significant differences in terms of their mer counts, the multiple mer densities by region, and the number of mers per transcript by organism. This finding suggests that post-transcriptional regulation varies significantly by organism – but this is not surprising since our original study [<xref ref-type="bibr" rid="CR4">4</xref>] sampled mRNAs in whole organisms in the case of the zebrafish and the organ/tissues of the brains and livers in the case of the mouse. Hence, the samples are not comparable and we would not expect post-transcriptional regulation to be the same in different organisms or organ/tissues.</p></sec><sec id="Sec25"><title>Unique Mers and known binding sites</title><p id="Par93">One set of unique mers with the sequence YUNNYUY apparently binds <italic>Hud</italic> proteins (Table <xref rid="Tab6" ref-type="table">6</xref>). <italic>Hud</italic> proteins stabilize mRNA by binding to AU-rich instability elements (AREs) in the 3’UTR and they target transcripts involved in neuronal differentiation, protein phosphatase regulation, ubiquitin ligation, and the transport, processing and translation of mRNAs [<xref ref-type="bibr" rid="CR21">21</xref>]. Interestingly, <italic>Hud</italic> proteins not only target their own mRNA but those of other RBPs, which suggests that it forms a network of post-transcriptional regulators [<xref ref-type="bibr" rid="CR21">21</xref>]. In the mouse, data from our previous study [<xref ref-type="bibr" rid="CR4">4</xref>] showed that <italic>Hud</italic> transcript abundance increased upon organismal death to reach maxima at 12 to 48 h postmortem (Fig. <xref rid="Fig10" ref-type="fig">10a</xref>). In the zebrafish, the <italic>Hud</italic> transcript abundance was about the same as the live controls for up to 4 h postmortem and then it declined and abruptly increased after 48 h (Fig. <xref rid="Fig10" ref-type="fig">10b</xref>). These findings are aligned with the notion that <italic>Hud</italic> proteins were involved in stabilizing mRNA in our study.<fig id="Fig10"><label>Fig. 10</label><caption><p>Gene transcript abundances measured by a calibrated microarray [<xref ref-type="bibr" rid="CR41">41</xref>, <xref ref-type="bibr" rid="CR42">42</xref>] (log transformed) by postmortem time. Abundances were normalized to flash frozen live controls (L). Black line, average. (<bold>a</bold>) Hud transcript in mouse; black dots, averaged abundance measured by probe A_55_P1990309 (n = 3 replicates for each dot except 48 h where n = 2 replicates); white dots, average abundance measured by probe A_55_P1990314; (<bold>b</bold>) Rbfox transcript in mouse; black dots, average abundance measured by probe A_55_P195339` (n = 3 replicates for each dot except last where n = 2 replicates); white dots, average abundance of probe A_55_P1953400; (<bold>c</bold>) Hud transcript in zebrafish; black dots, average abundance of probe A_15_P119510 (n = 2 replicates for each dot); white dots, average abundance of probe A_15_P120793. Data are from ref. [<xref ref-type="bibr" rid="CR4">4</xref>]</p></caption><graphic xlink:href="12864_2018_5042_Fig10_HTML" id="MO10"></graphic></fig></p><p id="Par94">Another unique mer with the sequence UGCAUG has previously been reported to serve as the binding site for <italic>Rbfox</italic> proteins that regulate splicing networks, mRNA stability and miRNA biogenesis [<xref ref-type="bibr" rid="CR10">10</xref>]. Apparently, the binding of <italic>Rbfox</italic> proteins to transcripts inhibits processing of the pri-microRNAs to pre-microRNAs, reduces expression of the mature miRNAs, and increases expression of targets normally downregulated by miRNAs [<xref ref-type="bibr" rid="CR10">10</xref>]. A previous study has shown that the abundance of transcripts with UGCAUG motifs in the 3’UTR positively correlates with <italic>Rbfox</italic> expression, and that knockdown of <italic>Rbfox</italic> decreases transcript abundances [<xref ref-type="bibr" rid="CR36">36</xref>]. These findings support the hypothesis that <italic>Rbfox</italic> enhances mRNA stability as well as gene expression. In our study, a little more than a third of the transcripts in the OP of the mouse have this binding site, but none were found in the OP of the zebrafish (Table <xref rid="Tab6" ref-type="table">6</xref>). In the mouse, data from our previous study [<xref ref-type="bibr" rid="CR4">4</xref>] showed that <italic>Rbfox</italic> transcript abundance increased after 30 min postmortem to reach a maximum at 48 h (Fig. <xref rid="Fig10" ref-type="fig">10c</xref>). These findings suggest that <italic>Rbfox</italic> proteins were interacting with some of the mouse transcripts in our previous study.</p><p id="Par95">The following 7 unique mers found in the OP have recently been reported as putative binding sites: UAUUUAU, GAGAAAA, AGAGAAA, UUUGCAC, AUGUGAA, UUGCACA, GGGAAGA [<xref ref-type="bibr" rid="CR34">34</xref>]. These sites have been correlated with increased gene expression in HeLa cells transfected with miRNAs. The UAUUUAU binding site is reported to be an ARE that signals rapid degradation or increased stability of mRNAs in response to stress [<xref ref-type="bibr" rid="CR36">36</xref>]. The Jacobsen et al. [<xref ref-type="bibr" rid="CR34">34</xref>] study showed that ARE binding sites and miRNA mediated regulation are interlinked, which is aligned with a similar study in <italic>Drosophila</italic> cells [<xref ref-type="bibr" rid="CR37">37</xref>]. While the significance and mechanistic insights of the 6 other putative binding sites were not discussed in the Jacobsen et al. study [<xref ref-type="bibr" rid="CR34">34</xref>], at least one of the seven binding sites was found in 258 of the 333 transcripts of the mouse and 106 of the 230 of the zebrafish, indicating that miRNAs might be involved in “regulating” the postmortem transcriptome (Table <xref rid="Tab6" ref-type="table">6</xref>).</p></sec><sec id="Sec26"><title>Post-transcriptional regulation of the postmortem transcriptome</title><p id="Par96">Several possible scenarios could be working in spatial and temporal combination to increase transcript stability and/or increase transcript abundance in the postmortem transcriptome. These scenarios are based, in part, on the “Competing endogenous RNA hypothesis”, which is provided at the end of the Discussion. However, without experimental evidence, we caution that these scenarios are speculative at best.</p><p id="Par97">One scenario is transcript stability is increased in the OP because they have more unique mers than the CP and RBPs bind to regulatory sites of transcripts of the OP blocking the binding of miRNAs, which are linked to degradation pathways. As a consequence, transcript stability is increased because the transcripts accumulate in the cells over time.</p><p id="Par98">A second scenario is postmortem genes are upregulated due to miRNA inhibition. Take, for example, transcripts regulated by p53 tumor suppressor that increase in abundance in response to miR-21 inhibition [<xref ref-type="bibr" rid="CR38">38</xref>].</p><p id="Par99">A third scenario is that some of the transcripts containing high multiple densities of mers act as molecular sponges that bind miRNAs and/or RBPs and therefore affect post-transcriptional regulation <italic>in trans</italic>. An example of this in our study was the 11 gene transcripts in the mouse with unknown functions and the <italic>Pimr</italic> transcripts in the zebrafish that had high densities in terms of mers per nucleotide (~ 2.4 mer/nt and ~ 1.0 mer/nt, respectively). Such high densities indicate that they contained many unique binding sites to sponge RBPs and/or ncRNAs. According to the data from our previous paper [<xref ref-type="bibr" rid="CR4">4</xref>], all the transcripts with high mer densities in the mouse increase in abundance right after death (0.5 h) and continued to increase, reaching a maximum abundance at 12 h, and then slowly decline (Fig. <xref rid="Fig11" ref-type="fig">11a</xref>). In the zebrafish, the <italic>Pimr</italic> gene transcripts increased slightly after death (relative to live controls) and abruptly increased after 12 h to maximize at 24 h (Fig. <xref rid="Fig11" ref-type="fig">11b</xref>). One-way to interpret these phenomena are that the transcripts are depleting the miRNA and/or RBP pools. In response to the decrease, a select group of genes involved in survival and stress compensation were passively transcribed, which accounts for the increases in transcript abundances in our original study.<fig id="Fig11"><label>Fig. 11</label><caption><p>Gene transcript abundances measured by a calibrated microarray [<xref ref-type="bibr" rid="CR41">41</xref>, <xref ref-type="bibr" rid="CR42">42</xref>] (log transformed) by postmortem time. Abundances were normalized to flash frozen live controls (L). Black line, average. (<bold>a</bold>) Mouse: Open circle, represents Gm11007, Gm2007, Gm4631, Gm14434, Gm2026, Gm14305, Gm14399, Gm14325, Zfp969, Gm4724, Gm14326 transcripts; closed circle, Zfp967, Zfp969, Zfp968; open square, Gm14410; closed square, Gm14305; open triangle, Gm14322; closed triangle, Gm14308; closed diamond, Gm14412. All points are the average of 3 replicates per sample time except the 48 h, which is the average of 2 replicates. (<bold>b</bold>) Zebrafish: Pimr transcript. Each point in the zebrafish represents the average of two individuals per sample time. Data are from ref. [<xref ref-type="bibr" rid="CR4">4</xref>]</p></caption><graphic xlink:href="12864_2018_5042_Fig11_HTML" id="MO11"></graphic></fig></p><p id="Par100">Further support for this scenario comes from the fact that most of the functional genes involved in survival and stress compensation were found in two clusters in the mouse: Groups A and G (59% of the OP) with low mer densities of 0.11 ± 0.12 mers/nt and 0.11 ± 0.05 mers/nt, respectively (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). In the zebrafish, most of the known functional gene transcripts are dispersed into groups A, C, F, K L, M, and N (93% of the OP) (Fig. <xref rid="Fig5" ref-type="fig">5</xref>), which have low mer densities (e.g., 0.10 ± 0.02 mers/nt). It is these genes that might have been passively upregulated due to lack of miRNA and RBPs to prevent them. This scenario makes sense for an evolutionary perspective because post-transcriptional regulation facilitates fast changes in response to stress so that cells can adapt to environmental change.</p></sec><sec id="Sec27"><title>Alternative splicing sites might differ under stress</title><p id="Par101">We assumed that the mRNA transcripts downloaded from NCBI represent dominant isoforms one would expect to find in nature. However, a recent study [<xref ref-type="bibr" rid="CR3">3</xref>] suggests that stress increases the production of different isoforms through alternative splicing. In other words, the composition of the transcripts might change in stressful conditions (i.e., different isoforms are produced). Our analysis did not account for this, however repeating our experiment using next-generation-sequencing methods might indeed provide additional insight into post-transcriptional regulation in postmortem gene expression, which is the focus of our future research.</p></sec><sec id="Sec28"><title>Competing endogenous RNA hypothesis</title><p id="Par102">According to the ‘competing endogenous RNA’ hypothesis, all types of RNA transcripts communicate through regulatory-binding sites and it is these interactions that regulate gene expression [<xref ref-type="bibr" rid="CR39">39</xref>]. The binding of miRNAs to sites represses translation and destabilizes the mRNA, thus having an overall negative regulatory role on gene expression. However, in the case when there is a limited pool of miRNAs to bind the sites or an overabundance of binding sites in transcripts, there is competition between targets to sequester miRNA. Thus, a surplus of binding sites dilutes the miRNA pool and gene expression resumes passively. Pseudogenes (i.e., those resembling known genes but are nonfunctional) as well as other transcripts can dilute the miRNA pool and thereby regulate their availability, and thus have an overall positive regulatory role on gene expression.</p><p id="Par103">Missing from the competing endogenous RNA hypothesis is the role of RBPs to compete with miRNA for regulatory binding sites. The presumed reason for this omission was at the time (i.e., 2011) there was a paucity of information supporting the idea that molecular sponges interact with proteins. However, proof exists today [<xref ref-type="bibr" rid="CR22">22</xref>]. A recent study reanalyzed highthroughput crosslinking and immunoprecipitation experiments in Human Embryonic Kidney Cells 293 to show that RBPs and miRNA often bind to the same or overlapping regulatory binding sites. The significance of this finding is twofold: (i) it suggests competition among the regulators (RBPs, miRNA, binding sites in different targets) and (ii) it suggests the relative concentrations of the RBPs and miRNAs to the regulatory binding sites might determine a transcript’s fate [<xref ref-type="bibr" rid="CR40">40</xref>].</p><p id="Par104">A third significant finding from the same study was the introduction of ‘hotspot’ binding sites that have high sequence conservation, accessibility, and enrichment in AU-rich elements (AREs) (i.e., devoid of guanines) and function by favoring competition among regulators [<xref ref-type="bibr" rid="CR40">40</xref>]. Apparently, target sites outside of hotspots have increased expression levels compared to targets sites within hotspots. Hence ‘hotspots’ are considered functional regulatory elements that provide an extra layer of regulation of post-transcriptional regulatory networks.</p></sec></sec><sec id="Sec29"><title>Conclusions</title><p id="Par105">This is the first study to investigate over-abundant mers in transcriptomic profiles after organismal death and raises interesting questions relative to post-transcriptional regulation and molecular biology.</p></sec><sec sec-type="supplementary-material"><title>Additional files</title><sec id="Sec30"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="12864_2018_5042_MOESM1_ESM.docx"><label>Additional file 1:</label><caption><p>Proof that using the ‘Chaos Genome Representation’ method to extract mers from transcript sequences is more practical (computational efficient) than string-based search algorithms. (DOCX 20 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM2"><media xlink:href="12864_2018_5042_MOESM2_ESM.xlsx"><label>Additional file 2:</label><caption><p>Two sheets in MS Excel file: (i) zebrafish_probes_PM, and (ii) mouse_probes_PM. PM, perfect match probes. Each sheet has two columns: first column is Agilent Probe ID and second column is DNA sequence. (XLSX 31 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM3"><media xlink:href="12864_2018_5042_MOESM3_ESM.fna"><label>Additional file 3:</label><caption><p>Two columns in text file of the over-abundant pool (OP) for the mouse. One column is Agilent Probe ID linked to Annotated Gene Name and second column is cDNA sequence. Total of 330 rows. (FNA 926 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM4"><media xlink:href="12864_2018_5042_MOESM4_ESM.fna"><label>Additional file 4:</label><caption><p>Two columns in text file of the over-abundant pool (OP) for the zebrafish. One column is Agilent Probe ID linked to Annotated Gene Name and second column is cDNA sequence. Total of 230 rows. (FNA 546 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM5"><media xlink:href="12864_2018_5042_MOESM5_ESM.fna"><label>Additional file 5:</label><caption><p>Two columns in text file of the control pool (CP) for the mouse. One column is Agilent Probe ID linked to Annotated Gene Name and second column is cDNA sequence. Total of 32,611 rows. (FNA 81851 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM6"><media xlink:href="12864_2018_5042_MOESM6_ESM.fna"><label>Additional file 6:</label><caption><p>Two columns in text file of the control pool (CP) for the zebrafish. One column is Agilent Probe ID linked to Annotated Gene Name and second column is cDNA sequence. Total of 27,433 rows. (FNA 129817 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM7"><media xlink:href="12864_2018_5042_MOESM7_ESM.xls"><label>Additional file 7:</label><caption><p>Two sheets in MS Excel file: (i) zebrafish, and (ii) mouse. Each sheet has 4 columns: first column is string length (strlen) of OP transcript; second column is blank; third column is string length of the corresponding CP transcript; fourth column is string length of corresponding CP2 transcript. Rows 1 to 231 in the zebrafish sheet contain the strlen of 230 transcripts in both OP and CP1 and CP2; rows 233 and 234 contains average and standard deviations of the columns; row 236 contains the two-tailed t-test results for OP vs CP1 and OP vs CP2. Rows 1 to 334 in the mouse sheet contain the strlen of 333 transcripts in both OP and CP1 and CP2; rows 336 and 337 contains average and standard deviations of the columns; row 339 contains the two-tailed t-test results for OP vs CP1 and OP vs CP2. (XLS 82 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM8"><media xlink:href="12864_2018_5042_MOESM8_ESM.xlsx"><label>Additional file 8:</label><caption><p>Nine sheets in MS Excel file of zebrafish data. The first sheet provides a detailed Readme that describes the sheets. Basically, first column is abundance of mer in OP, second column is average abundance in CP, third column is standard deviation in CP, and remaining 30 columns are abundances of 30 random draws from CP. Rows differ by mer length. (XLSX 61046 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM9"><media xlink:href="12864_2018_5042_MOESM9_ESM.xlsx"><label>Additional file 9:</label><caption><p>Nine sheets in MS Excel file of mouse data. The first sheet provides a detailed Readme that describes the sheets. Basically, first column is abundance of mer in OP, second column is average abundance in CP, third column is standard deviation in CP, and remaining 30 columns are abundances of 30 random draws from CP. Rows differ by mer length. (XLSX 62431 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM10"><media xlink:href="12864_2018_5042_MOESM10_ESM.xlsx"><label>Additional file 10:</label><caption><p>11 sheets in MS Excel file. It is similar to the Online Resource_8 zebrafish file except the raw data is missing to reduce matrix size. The purpose of the sheets is to calculate over- and under-abundant mers that are 5 X the standard deviation of the CP for each mer. The first sheet provides a detailed Readme that describes the sheets. Rows differ by mer length. (XLSX 65834 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM11"><media xlink:href="12864_2018_5042_MOESM11_ESM.xlsx"><label>Additional file 11:</label><caption><p>11 sheets in MS Excel file. It is similar to the Online Resource_9 mouse file except the raw data is missing to reduce matrix size. The purpose of the sheets is to calculate over- and under-abundant mers that are 5 X the standard deviation of the CP for each mer. The first sheet provides a detailed Readme that describes the sheets. Rows differ by mer length. (XLSX 69035 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM12"><media xlink:href="12864_2018_5042_MOESM12_ESM.xlsx"><label>Additional file 12:</label><caption><p>Multiple sheets in MS Excel file. The first sheet provides a detailed Readme that describes the sheets. Rows differ by mer length. The matrix file consists of 2245 columns and 230 rows. (XLSX 12475 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM13"><media xlink:href="12864_2018_5042_MOESM13_ESM.xlsx"><label>Additional file 13:</label><caption><p>Multiple sheets in MS Excel file. The first sheet provides a detailed Readme that describes the sheets. Rows differ by mer length. The matrix file consists of 5117 columns and 333 rows. (XLSX 40047 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM14"><media xlink:href="12864_2018_5042_MOESM14_ESM.xls"><label>Additional file 14:</label><caption><p>Two sheets in MS Excel file. The first sheet is the collapsed data of the Zebrafish and the second sheet is the collapsed data of the Mouse. Each sheet shows how the data was log normalized for making the heatmaps. The collapsed data was based on two way cluster groups using Wards linkage methods. (XLS 108 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM15"><media xlink:href="12864_2018_5042_MOESM15_ESM.xls"><label>Additional file 15:</label><caption><p>Four sheets in MS Excel file. The first sheet provides a detailed Readme that describes the sheets. The second and third sheets have the number of mer hits by transcript sequence length for the zebrafish and mouse. The fourth sheet has the summarize RegRNA2 output for 10 samples. (XLS 117 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM16"><media xlink:href="12864_2018_5042_MOESM16_ESM.xls"><label>Additional file 16:</label><caption><p>Two sheets in MS Excel file. The first sheet is the number of mers by region for the zebrafish and the second sheet is the same for the mouse. (XLS 86 kb)</p></caption></media></supplementary-material></p></sec></sec></body><back><glossary><title>Abbreviations</title><def-list><def-item><term>CP</term><def><p id="Par4">Control transcript pool</p></def></def-item><def-item><term>FP</term><def><p id="Par5">False positive</p></def></def-item><def-item><term>miRNA</term><def><p id="Par6">microRNA</p></def></def-item><def-item><term>ncRNA</term><def><p id="Par7">Non-coding RNA</p></def></def-item><def-item><term>OP</term><def><p id="Par8">Overabundant transcript pool</p></def></def-item><def-item><term>ORFs</term><def><p id="Par9">Open reading frames</p></def></def-item><def-item><term>RBP</term><def><p id="Par10">RNA binding proteins</p></def></def-item><def-item><term>UTR</term><def><p id="Par11">Untranslated regions</p></def></def-item></def-list></glossary><ack><title>Acknowledgements</title><p>We thank Shivani Soni for reviewing earlier versions of the manuscript. We also thank the anonymous reviewers for their insights in earlier versions of the manuscript.</p><sec id="FPar1"><title>Funding</title><p id="Par106">There was no funding source for this study. The authors personally paid the publication costs for this article.</p></sec><sec id="FPar2"><title>Availability of data and materials</title><p id="Par107">All data generated or analyzed during this study are included in this published article and its supplementary information files.</p></sec></ack><notes notes-type="author-contribution"><title>Authors’ contributions</title><p>PAN and AEP designed the study. PAN performed the bioinformatic and statistical analyses, and made the figures and tables. Both authors read and approved the final manuscript.</p></notes><notes notes-type="COI-statement"><sec id="FPar3"><title>Ethics approval and consent to participate</title><p>Not applicable, all data publicly available.</p></sec><sec id="FPar4"><title>Consent for publication</title><p>Not applicable.</p></sec><sec id="FPar5"><title>Competing interests</title><p>The authors declare that they have no competing interests.</p></sec><sec id="FPar6"><title>Publisher’s Note</title><p>Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p></sec></notes><ref-list id="Bib1"><title>References</title><ref id="CR1"><label>1.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Harvey</surname><given-names>R</given-names></name><name><surname>Dezi</surname><given-names>V</given-names></name><name><surname>Pizzinga</surname><given-names>M</given-names></name><name><surname>Willis</surname><given-names>AE</given-names></name></person-group><article-title>Post-transcriptional control of gene expression following stress: the role of RNA-binding proteins</article-title><source>Biochem Soc T</source><year>2017</year><volume>45</volume><fpage>1007</fpage><lpage>1014</lpage><pub-id pub-id-type="doi">10.1042/BST20160364</pub-id></element-citation></ref><ref id="CR2"><label>2.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hogg</surname><given-names>JR</given-names></name><name><surname>Collins</surname><given-names>K</given-names></name></person-group><article-title>Structured non-coding RNAs and the RNP renaissance</article-title><source>Curr Opin Chem Biol</source><year>2008</year><volume>12</volume><fpage>684</fpage><lpage>689</lpage><pub-id pub-id-type="doi">10.1016/j.cbpa.2008.09.027</pub-id><pub-id pub-id-type="pmid">18950732</pub-id></element-citation></ref><ref id="CR3"><label>3.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ferreira</surname><given-names>PG</given-names></name><name><surname>Muñoz-Aguirre</surname><given-names>M</given-names></name><name><surname>Reverter</surname><given-names>F</given-names></name><name><surname>Sá Godinho</surname><given-names>CP</given-names></name><name><surname>Sousa</surname><given-names>A</given-names></name><name><surname>Amadoz</surname><given-names>A</given-names></name><name><surname>Sodaei</surname><given-names>R</given-names></name><name><surname>Hidalgo</surname><given-names>MR</given-names></name><name><surname>Pervouchine</surname></name><name><surname>Dmitri</surname></name><name><surname>Carbonell-Caballero</surname><given-names>J</given-names></name><name><surname>Nurtdinov</surname><given-names>R</given-names></name><name><surname>Breschi</surname><given-names>A</given-names></name><name><surname>Amador</surname><given-names>R</given-names></name><name><surname>Oliveira</surname><given-names>P</given-names></name><name><surname>Çubuk</surname><given-names>C</given-names></name><name><surname>Curado</surname><given-names>J</given-names></name><name><surname>Aguet</surname><given-names>F</given-names></name><name><surname>Oliveira</surname><given-names>C</given-names></name><name><surname>Dopazo</surname><given-names>J</given-names></name><name><surname>Sammeth</surname><given-names>M</given-names></name><name><surname>Ardlie</surname><given-names>KG</given-names></name><name><surname>Guigó</surname><given-names>R</given-names></name></person-group><article-title>The effects of death and post-mortem cold ischemia on human tissue transcriptomes</article-title><source>Nat Commun</source><year>2018</year><volume>490</volume><fpage>2041</fpage><lpage>1723</lpage><pub-id pub-id-type="doi">10.1038/s41467-017-02772-x</pub-id></element-citation></ref><ref id="CR4"><label>4.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pozhitkov</surname><given-names>AE</given-names></name><name><surname>Neme</surname><given-names>R</given-names></name><name><surname>Domazet-Loo</surname><given-names>T</given-names></name><name><surname>Leroux</surname><given-names>BG</given-names></name><name><surname>Soni</surname><given-names>S</given-names></name><name><surname>Tautz</surname><given-names>D</given-names></name><name><surname>Noble</surname><given-names>PA</given-names></name></person-group><article-title>Tracing the dynamics of gene transcripts after organismal death</article-title><source>Open Biol</source><year>2017</year><volume>7</volume><fpage>160267</fpage><pub-id pub-id-type="doi">10.1098/rsob.160267</pub-id><pub-id pub-id-type="pmid">28123054</pub-id></element-citation></ref><ref id="CR5"><label>5.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dykes</surname><given-names>IM</given-names></name><name><surname>Emanueli</surname><given-names>C</given-names></name></person-group><article-title>Transcriptional and post-transcriptional gene regulation by long non-coding RNA</article-title><source>Genomics Proteomics Bioinformatics</source><year>2017</year><volume>15</volume><fpage>177</fpage><lpage>186</lpage><pub-id pub-id-type="doi">10.1016/j.gpb.2016.12.005</pub-id><pub-id pub-id-type="pmid">28529100</pub-id></element-citation></ref><ref id="CR6"><label>6.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>He</surname><given-names>JH</given-names></name><name><surname>Han</surname><given-names>ZP</given-names></name><name><surname>Liu</surname><given-names>JM</given-names></name><name><surname>Zhou</surname><given-names>JB</given-names></name><name><surname>Zou</surname><given-names>MX</given-names></name><name><surname>Lv</surname><given-names>YB</given-names></name><name><surname>Li</surname><given-names>YG</given-names></name><name><surname>Cao</surname><given-names>MR</given-names></name></person-group><article-title>Overexpression of long non-coding RNA Meg3 inhibits proliferation of hepatocellular carcinoma Huh7 cells via negative modulation of miRNA-664</article-title><source>J Cell Biochem</source><year>2017</year><volume>118</volume><fpage>3713</fpage><lpage>3721</lpage><pub-id pub-id-type="doi">10.1002/jcb.26018</pub-id><pub-id pub-id-type="pmid">28374914</pub-id></element-citation></ref><ref id="CR7"><label>7.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gerstberger</surname><given-names>S</given-names></name><name><surname>Hafner</surname><given-names>M</given-names></name><name><surname>Tuschl</surname><given-names>T</given-names></name></person-group><article-title>A census of human RNA-binding proteins</article-title><source>Nat Rev Gene</source><year>2014</year><volume>15</volume><fpage>829</fpage><lpage>845</lpage><pub-id pub-id-type="doi">10.1038/nrg3813</pub-id></element-citation></ref><ref id="CR8"><label>8.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maraia</surname><given-names>Richard J.</given-names></name><name><surname>Mattijssen</surname><given-names>Sandy</given-names></name><name><surname>Cruz-Gallardo</surname><given-names>Isabel</given-names></name><name><surname>Conte</surname><given-names>Maria R.</given-names></name></person-group><article-title>The La and related RNA-binding proteins (LARPs): structures, functions, and evolving perspectives</article-title><source>Wiley Interdisciplinary Reviews: RNA</source><year>2017</year><volume>8</volume><issue>6</issue><fpage>e1430</fpage><pub-id pub-id-type="doi">10.1002/wrna.1430</pub-id></element-citation></ref><ref id="CR9"><label>9.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Castella</surname><given-names>S</given-names></name><name><surname>Bernard</surname><given-names>R</given-names></name><name><surname>Corno</surname><given-names>M</given-names></name><name><surname>Fradin</surname><given-names>A</given-names></name><name><surname>Larcher</surname><given-names>JC</given-names></name></person-group><article-title>Ilf3 and NF90 functions in RNA biology</article-title><source>Wiley Interdiscip Rev RNA</source><year>2015</year><volume>6</volume><fpage>243</fpage><lpage>256</lpage><pub-id pub-id-type="doi">10.1002/wrna.1270</pub-id><pub-id pub-id-type="pmid">25327818</pub-id></element-citation></ref><ref id="CR10"><label>10.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Conboy</surname><given-names>John G.</given-names></name></person-group><article-title>Developmental regulation of RNA processing by Rbfox proteins</article-title><source>Wiley Interdisciplinary Reviews: RNA</source><year>2016</year><volume>8</volume><issue>2</issue><fpage>e1398</fpage><pub-id pub-id-type="doi">10.1002/wrna.1398</pub-id></element-citation></ref><ref id="CR11"><label>11.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Soengas</surname><given-names>MS</given-names></name><name><surname>Hernando</surname><given-names>E</given-names></name></person-group><article-title>Tyrp1 mRNA goes fishing for miRNAs in melanoma</article-title><source>Nat Cell Biol</source><year>2017</year><volume>19</volume><fpage>1311</fpage><lpage>1312</lpage><pub-id pub-id-type="doi">10.1038/ncb3637</pub-id><pub-id pub-id-type="pmid">29087386</pub-id></element-citation></ref><ref id="CR12"><label>12.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gilot</surname><given-names>D</given-names></name><name><surname>Migault</surname><given-names>M</given-names></name><name><surname>Bachelot</surname><given-names>L</given-names></name><name><surname>Journé</surname><given-names>F</given-names></name><name><surname>Rogiers</surname><given-names>A</given-names></name><name><surname>Donnou-Fournet</surname><given-names>E</given-names></name><name><surname>Mogha</surname><given-names>A</given-names></name><name><surname>Mouchet</surname><given-names>N</given-names></name><name><surname>Pinel-Marie</surname><given-names>ML</given-names></name><name><surname>Mari</surname><given-names>B</given-names></name><name><surname>Montier</surname><given-names>T</given-names></name><name><surname>Corre</surname><given-names>S</given-names></name><name><surname>Gautron</surname><given-names>A</given-names></name><name><surname>Rambow</surname><given-names>F</given-names></name><name><surname>El Hajj</surname><given-names>P</given-names></name><name><surname>Ben Jouira</surname><given-names>R</given-names></name><name><surname>Tartare-Deckert</surname><given-names>S</given-names></name><name><surname>Marine</surname><given-names>JC</given-names></name><name><surname>Felden</surname><given-names>B</given-names></name><name><surname>Ghanem</surname><given-names>G</given-names></name><name><surname>Galibert</surname><given-names>MD</given-names></name></person-group><article-title>A non-coding function of TYRP1 mRNA promotes melanoma growth</article-title><source>Nat Cell Biol</source><year>2017</year><volume>19</volume><fpage>1348</fpage><lpage>1357</lpage><pub-id pub-id-type="doi">10.1038/ncb3623</pub-id><pub-id pub-id-type="pmid">28991221</pub-id></element-citation></ref><ref id="CR13"><label>13.</label><mixed-citation publication-type="other">Liu J, Liu T, Wang X, He A. Circles reshaping the RNA world: from waste to treasure. Mol Cancer. 2017;16(58) 10.1186/s12943-017-0630-y. PMID: 28279183</mixed-citation></ref><ref id="CR14"><label>14.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Łabno</surname><given-names>A</given-names></name><name><surname>Tomecki</surname><given-names>R</given-names></name><name><surname>Dziembowski</surname><given-names>A</given-names></name></person-group><article-title>Cytoplasmic RNA decay pathways - enzymes and mechanisms</article-title><source>Biochim Biophys Acta</source><year>2016</year><volume>1863</volume><fpage>3125</fpage><lpage>3147</lpage><pub-id pub-id-type="doi">10.1016/j.bbamcr.2016.09.023</pub-id><pub-id pub-id-type="pmid">27713097</pub-id></element-citation></ref><ref id="CR15"><label>15.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wilusz</surname><given-names>JE</given-names></name><name><surname>Sunwoo</surname><given-names>H</given-names></name><name><surname>Spector</surname><given-names>DL</given-names></name></person-group><article-title>Long noncoding RNAs: functional surprises from the RNA world</article-title><source>Genes Dev</source><year>2009</year><volume>23</volume><fpage>1494</fpage><lpage>1504</lpage><pub-id pub-id-type="doi">10.1101/gad.1800909</pub-id><pub-id pub-id-type="pmid">19571179</pub-id></element-citation></ref><ref id="CR16"><label>16.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eddy</surname><given-names>SR</given-names></name></person-group><article-title>Non-coding RNA genes and the modern RNA world</article-title><source>Nat Rev Genet</source><year>2001</year><volume>2</volume><fpage>919</fpage><lpage>929</lpage><pub-id pub-id-type="doi">10.1038/35103511</pub-id><pub-id pub-id-type="pmid">11733745</pub-id></element-citation></ref><ref id="CR17"><label>17.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jeffrey</surname><given-names>HJ</given-names></name></person-group><article-title>Chaos game representation of gene structure</article-title><source>Nucleic Acids Res</source><year>1990</year><volume>18</volume><fpage>2163</fpage><lpage>2170</lpage><pub-id pub-id-type="doi">10.1093/nar/18.8.2163</pub-id><pub-id pub-id-type="pmid">2336393</pub-id></element-citation></ref><ref id="CR18"><label>18.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Noble</surname><given-names>PA</given-names></name><name><surname>Citek</surname><given-names>RW</given-names></name><name><surname>Ogunseitan</surname><given-names>OA</given-names></name></person-group><article-title>Tetranucleotide frequencies in microbial genomes</article-title><source>Electrophoresis</source><year>1998</year><volume>19</volume><fpage>528</fpage><lpage>535</lpage><pub-id pub-id-type="doi">10.1002/elps.1150190412</pub-id><pub-id pub-id-type="pmid">9588798</pub-id></element-citation></ref><ref id="CR19"><label>19.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Almeida</surname><given-names>JS</given-names></name><name><surname>Carriço</surname><given-names>JA</given-names></name><name><surname>Maretzek</surname><given-names>A</given-names></name><name><surname>Noble</surname><given-names>PA</given-names></name><name><surname>Fletcher</surname><given-names>M</given-names></name></person-group><article-title>Analysis of genomic sequences by Chaos game representation</article-title><source>Bioinformatics</source><year>2001</year><volume>17</volume><fpage>429</fpage><lpage>437</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/17.5.429</pub-id><pub-id pub-id-type="pmid">11331237</pub-id></element-citation></ref><ref id="CR20"><label>20.</label><mixed-citation publication-type="other">Chang TH, Huang HY, Hsu JB, Weng SL, Horng JT, Huang HD. An enhanced computational platform for investigating the roles of regulatory RNA and for identifying functional RNA motifs. BMC Bioinformatics Suppl 2013;2:S4. 10.1186/1471-2105-14-S2-S4 PMID: 23369107.</mixed-citation></ref><ref id="CR21"><label>21.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bolognani</surname><given-names>F</given-names></name><name><surname>Contente-Cuomo</surname><given-names>T</given-names></name><name><surname>Perrone-Bizzozero</surname><given-names>NI</given-names></name></person-group><article-title>Novel recognition motifs and biological functions of the RNA-binding protein HuD revealed by genome-wide identification of its targets</article-title><source>Nucleic Acids Res</source><year>2010</year><volume>38</volume><fpage>117</fpage><lpage>130</lpage><pub-id pub-id-type="doi">10.1093/nar/gkp863</pub-id><pub-id pub-id-type="pmid">19846595</pub-id></element-citation></ref><ref id="CR22"><label>22.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jacobsen</surname><given-names>A</given-names></name><name><surname>Wen</surname><given-names>J</given-names></name><name><surname>Marks</surname><given-names>DS</given-names></name><name><surname>Krogh</surname><given-names>A</given-names></name></person-group><article-title>Signatures of RNA binding proteins globally coupled to effective microRNA target sites</article-title><source>Genome Res</source><year>2010</year><volume>20</volume><fpage>1010</fpage><lpage>1019</lpage><pub-id pub-id-type="doi">10.1101/gr.103259.109</pub-id><pub-id pub-id-type="pmid">20508147</pub-id></element-citation></ref><ref id="CR23"><label>23.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barreau</surname><given-names>C</given-names></name><name><surname>Paillard</surname><given-names>L</given-names></name><name><surname>Osborne</surname><given-names>HB</given-names></name></person-group><article-title>AU-rich elements and associated factors: are there unifying principles?</article-title><source>Nucleic Acids Res</source><year>2006</year><volume>33</volume><fpage>7138</fpage><lpage>7150</lpage><pub-id pub-id-type="doi">10.1093/nar/gki1012</pub-id><pub-id pub-id-type="pmid">16391004</pub-id></element-citation></ref><ref id="CR24"><label>24.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ding</surname><given-names>CL</given-names></name><name><surname>Xu</surname><given-names>G</given-names></name><name><surname>Tang</surname><given-names>HL</given-names></name><name><surname>Zhu</surname><given-names>SY</given-names></name><name><surname>Zhao</surname><given-names>LJ</given-names></name><name><surname>Ren</surname><given-names>H</given-names></name><name><surname>Zhao</surname><given-names>P</given-names></name><name><surname>Qi</surname><given-names>ZT</given-names></name><name><surname>Wang</surname><given-names>W</given-names></name></person-group><article-title>Anchoring of both PKA-RIIα and 14-3-3θ regulates retinoic acid induced 16 mediated phosphorylation of heat shock protein 70</article-title><source>Oncotarget</source><year>2015</year><volume>6</volume><fpage>15540</fpage><lpage>15550</lpage><pub-id pub-id-type="doi">10.18632/oncotarget.3702</pub-id><pub-id pub-id-type="pmid">25900241</pub-id></element-citation></ref><ref id="CR25"><label>25.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname><given-names>LC</given-names></name><name><surname>Jiang</surname><given-names>Z</given-names></name><name><surname>Sang</surname><given-names>Y</given-names></name><name><surname>Harhay</surname><given-names>GP</given-names></name><name><surname>Lager</surname><given-names>KM</given-names></name></person-group><article-title>Evolutionary characterization of pig interferon-inducible transmembrane gene family and member expression dynamics in tracheobronchial lymph nodes of pigs infected with swine respiratory disease viruses</article-title><source>Vet Immunol Immunopathol</source><year>2014</year><volume>159</volume><fpage>180</fpage><lpage>191</lpage><pub-id pub-id-type="doi">10.1016/j.vetimm.2014.02.015</pub-id><pub-id pub-id-type="pmid">24656980</pub-id></element-citation></ref><ref id="CR26"><label>26.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cruse</surname><given-names>G</given-names></name><name><surname>Beaven</surname><given-names>MA</given-names></name><name><surname>Music</surname><given-names>SC</given-names></name><name><surname>Bradding</surname><given-names>P</given-names></name><name><surname>Gilfillan</surname><given-names>AM</given-names></name><name><surname>Metcalfe</surname><given-names>DD</given-names></name></person-group><article-title>The CD20 homologue MS4A4 directs trafficking of KIT toward clathrin-independent endocytosis pathways and thus regulates receptor signaling and recycling</article-title><source>Mol Biol Cell</source><year>2015</year><volume>26</volume><fpage>1711</fpage><lpage>1727</lpage><pub-id pub-id-type="doi">10.1091/mbc.E14-07-1221</pub-id><pub-id pub-id-type="pmid">25717186</pub-id></element-citation></ref><ref id="CR27"><label>27.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bens</surname><given-names>S</given-names></name><name><surname>Zichner</surname><given-names>T</given-names></name><name><surname>Stütz</surname><given-names>AM</given-names></name><name><surname>Caliebe</surname><given-names>A</given-names></name><name><surname>Wagener</surname><given-names>R</given-names></name><name><surname>Hoff</surname><given-names>K</given-names></name><name><surname>Korbel</surname><given-names>JO</given-names></name><name><surname>von Bismarck</surname><given-names>P</given-names></name><name><surname>Siebert</surname><given-names>R</given-names></name></person-group><article-title>SPAG7 is a candidate gene for the periodic fever, aphthous stomatitis, pharyngitis and adenopathy (PFAPA) syndrome</article-title><source>Genes Immun</source><year>2014</year><volume>15</volume><fpage>190</fpage><lpage>194</lpage><pub-id pub-id-type="doi">10.1038/gene.2013.73</pub-id><pub-id pub-id-type="pmid">24452265</pub-id></element-citation></ref><ref id="CR28"><label>28.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fernebro</surname><given-names>J</given-names></name><name><surname>Francis</surname><given-names>P</given-names></name><name><surname>Edén</surname><given-names>P</given-names></name><name><surname>Borg</surname><given-names>A</given-names></name><name><surname>Panagopoulos</surname><given-names>I</given-names></name><name><surname>Mertens</surname><given-names>F</given-names></name><name><surname>Vallon-Christersson</surname><given-names>J</given-names></name><name><surname>Akerman</surname><given-names>M</given-names></name><name><surname>Rydholm</surname><given-names>A</given-names></name><name><surname>Bauer</surname><given-names>HC</given-names></name><name><surname>Mandahl</surname><given-names>N</given-names></name><name><surname>Nilbert</surname><given-names>M</given-names></name></person-group><article-title>Gene expression profiles relate to SS18/SSX fusion type in synovial sarcoma</article-title><source>Int J Cancer</source><year>2006</year><volume>118</volume><fpage>1165</fpage><lpage>1172</lpage><pub-id pub-id-type="doi">10.1002/ijc.21475</pub-id><pub-id pub-id-type="pmid">16152617</pub-id></element-citation></ref><ref id="CR29"><label>29.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kerr</surname><given-names>JR</given-names></name><name><surname>Kaushik</surname><given-names>N</given-names></name><name><surname>Fear</surname><given-names>D</given-names></name><name><surname>Baldwin</surname><given-names>DA</given-names></name><name><surname>Nuwaysir</surname><given-names>EF</given-names></name><name><surname>Adcock</surname><given-names>IM</given-names></name></person-group><article-title>Single-nucleotide polymorphisms associated with symptomatic infection and differential human gene expression in healthy seropositive persons each implicate the cytoskeleton, integrin signaling, and oncosuppression in the pathogenesis of human parvovirus B19 infection</article-title><source>J Infect Dis</source><year>2005</year><volume>192</volume><fpage>276</fpage><lpage>286</lpage><pub-id pub-id-type="doi">10.1086/430950</pub-id><pub-id pub-id-type="pmid">15962222</pub-id></element-citation></ref><ref id="CR30"><label>30.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>L</given-names></name><name><surname>Ma</surname><given-names>C</given-names></name><name><surname>Bian</surname><given-names>Y</given-names></name><name><surname>Shao</surname><given-names>C</given-names></name><name><surname>Wang</surname><given-names>T</given-names></name><name><surname>Li</surname><given-names>J</given-names></name><name><surname>Chong</surname><given-names>X</given-names></name><name><surname>Su</surname><given-names>L</given-names></name><name><surname>Lu</surname><given-names>J</given-names></name></person-group><article-title>Aberrant expression of STYK1 and E-cadherin confer a poor prognosis for pancreatic cancer patients</article-title><source>Oncotarget</source><year>2017</year><volume>8</volume><fpage>111333</fpage><lpage>111345</lpage><pub-id pub-id-type="doi">10.18632/oncotarget.22794</pub-id><pub-id pub-id-type="pmid">29340057</pub-id></element-citation></ref><ref id="CR31"><label>31.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname><given-names>Yunfei</given-names></name><name><surname>Yang</surname><given-names>Ling</given-names></name><name><surname>He</surname><given-names>Jiao</given-names></name><name><surname>Yang</surname><given-names>Huai</given-names></name></person-group><article-title>STYK1 promotes Warburg effect through PI3K/AKT signaling and predicts a poor prognosis in nasopharyngeal carcinoma</article-title><source>Tumor Biology</source><year>2017</year><volume>39</volume><issue>7</issue><fpage>101042831771164</fpage><pub-id pub-id-type="doi">10.1177/1010428317711644</pub-id></element-citation></ref><ref id="CR32"><label>32.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>J</given-names></name><name><surname>Wang</surname><given-names>F</given-names></name><name><surname>Liu</surname><given-names>B</given-names></name><name><surname>Yang</surname><given-names>L</given-names></name><name><surname>Wang</surname><given-names>X</given-names></name><name><surname>Liu</surname><given-names>Y</given-names></name></person-group><article-title>Knockdown of serine threonine tyrosine kinase 1 (STYK1) inhibits the migration and tumorigenesis in glioma cells</article-title><source>Oncol Res</source><year>2017</year><volume>25</volume><fpage>931</fpage><lpage>937</lpage><pub-id pub-id-type="doi">10.3727/096504016X14772424117423</pub-id><pub-id pub-id-type="pmid">27983928</pub-id></element-citation></ref><ref id="CR33"><label>33.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arch</surname><given-names>RH</given-names></name><name><surname>Thompson</surname><given-names>CB</given-names></name></person-group><article-title>4-1BB and Ox40 are members of a tumor necrosis factor (TNF)-nerve growth factor receptor subfamily that bind TNF receptor-associated factors and activate nuclear factor kappaB</article-title><source>Mol Cell Biol</source><year>1998</year><volume>18</volume><fpage>558</fpage><lpage>565</lpage><pub-id pub-id-type="doi">10.1128/MCB.18.1.558</pub-id><pub-id pub-id-type="pmid">9418902</pub-id></element-citation></ref><ref id="CR34"><label>34.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mignone</surname><given-names>Flavio</given-names></name><name><surname>Gissi</surname><given-names>Carmela</given-names></name><name><surname>Liuni</surname><given-names>Sabino</given-names></name><name><surname>Pesole</surname><given-names>Graziano</given-names></name></person-group><source>Genome Biology</source><year>2002</year><volume>3</volume><issue>3</issue><fpage>reviews0004.1</fpage><pub-id pub-id-type="doi">10.1186/gb-2002-3-3-reviews0004</pub-id><pub-id pub-id-type="pmid">11897027</pub-id></element-citation></ref><ref id="CR35"><label>35.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tay</surname><given-names>Y</given-names></name><name><surname>Zhang</surname><given-names>J</given-names></name><name><surname>Thomson</surname><given-names>AM</given-names></name><name><surname>Lim</surname><given-names>B</given-names></name><name><surname>Rigoutsos</surname><given-names>I</given-names></name></person-group><article-title>MicroRNAs to Nanog, Oct4 and Sox2 coding regions modulate embryonic stem cell differentiation</article-title><source>Nature</source><year>2008</year><volume>455</volume><fpage>1124</fpage><lpage>1128</lpage><pub-id pub-id-type="doi">10.1038/nature07299</pub-id><pub-id pub-id-type="pmid">18806776</pub-id></element-citation></ref><ref id="CR36"><label>36.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ray</surname><given-names>D</given-names></name><name><surname>Kazan</surname><given-names>H</given-names></name><name><surname>Cook</surname><given-names>KB</given-names></name><name><surname>Weirauch</surname><given-names>MT</given-names></name><name><surname>Najafabadi</surname><given-names>HS</given-names></name><name><surname>Li</surname><given-names>X</given-names></name><name><surname>Gueroussov</surname><given-names>S</given-names></name><name><surname>Albu</surname><given-names>M</given-names></name><name><surname>Zheng</surname><given-names>H</given-names></name><name><surname>Yang</surname><given-names>A</given-names></name><name><surname>Na</surname><given-names>H</given-names></name><name><surname>Irimia</surname><given-names>M</given-names></name><name><surname>Matzat</surname><given-names>LH</given-names></name><name><surname>Dale</surname><given-names>RK</given-names></name><name><surname>Smith</surname><given-names>SA</given-names></name><name><surname>Yarosh</surname><given-names>CA</given-names></name><name><surname>Kelly</surname><given-names>SM</given-names></name><name><surname>Nabet</surname><given-names>B</given-names></name><name><surname>Mecenas</surname><given-names>D</given-names></name><name><surname>Li</surname><given-names>W</given-names></name><name><surname>Laishram</surname><given-names>RS</given-names></name><name><surname>Qiao</surname><given-names>M</given-names></name><name><surname>Lipshitz</surname><given-names>HD</given-names></name><name><surname>Piano</surname><given-names>F</given-names></name><name><surname>Corbett</surname><given-names>AH</given-names></name><name><surname>Carstens</surname><given-names>RP</given-names></name><name><surname>Frey</surname><given-names>BJ</given-names></name><name><surname>Anderson</surname><given-names>RA</given-names></name><name><surname>Lynch</surname><given-names>KW</given-names></name><name><surname>Penalva</surname><given-names>LO</given-names></name><name><surname>Lei</surname><given-names>EP</given-names></name><name><surname>Fraser</surname><given-names>AG</given-names></name><name><surname>Blencowe</surname><given-names>BJ</given-names></name><name><surname>Morris</surname><given-names>QD</given-names></name><name><surname>Hughes</surname><given-names>TR</given-names></name></person-group><article-title>A compendium of RNA-binding motifs for decoding gene regulation</article-title><source>Nature</source><year>2013</year><volume>499</volume><fpage>172</fpage><lpage>177</lpage><pub-id pub-id-type="doi">10.1038/nature12311</pub-id><pub-id pub-id-type="pmid">23846655</pub-id></element-citation></ref><ref id="CR37"><label>37.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Geisberg</surname><given-names>JV</given-names></name><name><surname>Moqtaderi</surname><given-names>Z</given-names></name><name><surname>Fan</surname><given-names>X</given-names></name><name><surname>Ozsolak</surname><given-names>F</given-names></name><name><surname>Struhl</surname><given-names>K</given-names></name></person-group><article-title>Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast</article-title><source>Cell</source><year>2014</year><volume>156</volume><fpage>812</fpage><lpage>824</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2013.12.026</pub-id><pub-id pub-id-type="pmid">24529382</pub-id></element-citation></ref><ref id="CR38"><label>38.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Frankel</surname><given-names>LB</given-names></name><name><surname>Christoffersen</surname><given-names>NR</given-names></name><name><surname>Jacobsen</surname><given-names>A</given-names></name><name><surname>Lindow</surname><given-names>M</given-names></name><name><surname>Krogh</surname><given-names>A</given-names></name><name><surname>Lund</surname><given-names>AH</given-names></name></person-group><article-title>Programmed cell death 4 (PDCD4) is an important functional target of the microRNA miR-21 in breast cancer cells</article-title><source>J Biol Chem</source><year>2008</year><volume>283</volume><fpage>1026</fpage><lpage>1033</lpage><pub-id pub-id-type="doi">10.1074/jbc.M707224200</pub-id><pub-id pub-id-type="pmid">17991735</pub-id></element-citation></ref><ref id="CR39"><label>39.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salmena</surname><given-names>L</given-names></name><name><surname>Poliseno</surname><given-names>L</given-names></name><name><surname>Tay</surname><given-names>Y</given-names></name><name><surname>Kats</surname><given-names>L</given-names></name><name><surname>Pandolfi</surname><given-names>PP</given-names></name></person-group><article-title>A ceRNA hypothesis: the Rosetta stone of a hidden RNA language?</article-title><source>Cell</source><year>2011</year><volume>146</volume><fpage>353</fpage><lpage>358</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2011.07.014</pub-id><pub-id pub-id-type="pmid">21802130</pub-id></element-citation></ref><ref id="CR40"><label>40.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plass</surname><given-names>M</given-names></name><name><surname>Rasmussen</surname><given-names>SH</given-names></name><name><surname>Krogh</surname><given-names>A</given-names></name></person-group><article-title>Highly accessible AU-rich regions in 3′ untranslated regions are hotspots for binding of regulatory factors</article-title><source>PLoS Comput Biol</source><year>2017</year><volume>13</volume><fpage>e1005460</fpage><pub-id pub-id-type="doi">10.1371/journal.pcbi.1005460</pub-id><pub-id pub-id-type="pmid">28410363</pub-id></element-citation></ref><ref id="CR41"><label>41.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pozhitkov</surname><given-names>AE</given-names></name><name><surname>Noble</surname><given-names>PA</given-names></name><name><surname>Bryk</surname><given-names>J</given-names></name><name><surname>Tautz</surname><given-names>D</given-names></name></person-group><article-title>A revised design for microarray experiments to account for experimental noise and uncertainty of probe response</article-title><source>PLoS One</source><year>2014</year><volume>9</volume><fpage>e91295</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0091295</pub-id><pub-id pub-id-type="pmid">24618910</pub-id></element-citation></ref><ref id="CR42"><label>42.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pozhitkov</surname><given-names>AE</given-names></name><name><surname>Noble</surname><given-names>PA</given-names></name></person-group><article-title>Gene meter: accurate abundance calculations of gene expression</article-title><source>Commun Integr Biol</source><year>2017</year><volume>10</volume><fpage>e1329785</fpage><pub-id pub-id-type="doi">10.1080/19420889.2017.1329785</pub-id><pub-id pub-id-type="pmid">28919937</pub-id></element-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0002970 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0002970 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= MersV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |