Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities
Identifieur interne : 000168 ( Pmc/Corpus ); précédent : 000167; suivant : 000169Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities
Auteurs : Leyden Fernandez ; Josep M. Mercader ; Mercè Planas-Fèlix ; David TorrentsSource :
- BMC Genomics [ 1471-2164 ] ; 2014.
Abstract
It has been shown in a number of metagenomic studies that the addition and removal of specific genes have allowed microbiomes to adapt to specific environmental conditions by losing and gaining specific functions. But it is not known whether and how the regulation of gene expression also contributes to adaptation.
We have here characterized and analyzed the metaregulome of three different environments, as well as their impact in the adaptation to particular variable physico-chemical conditions. For this, we have developed a computational protocol to extract regulatory regions and their corresponding transcription factors binding sites directly from metagenomic reads and applied it to three well known environments: Acid Mine, Whale Fall, and Waseca Farm. Taking the density of regulatory sites in promoters as a measure of the potential and complexity of gene regulation, we found it to be quantitatively the same in all three environments, despite their different physico-chemical conditions and species composition. However, we found that each environment distributes their regulatory potential differently across their functional space. Among the functions with highest regulatory potential in each niche, we found significant enrichment of processes related to sensing and buffering external variable factors specific to each environment, like for example, the availability of co-factors in deep sea, of oligosaccharides in soil and the regulation of pH in the acid mine.
These results highlight the potential impact of gene regulation in the adaptation of bacteria to the different habitats through the distribution of their regulatory potential among specific functions, and point to critical environmental factors that challenge the growth of any microbial community.
The online version of this article (doi:10.1186/1471-2164-15-877) contains supplementary material, which is available to authorized users.
Url:
DOI: 10.1186/1471-2164-15-877
PubMed: 25294412
PubMed Central: 4287501
Links to Exploration step
PMC:4287501Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities</title>
<author><name sortKey="Fernandez, Leyden" sort="Fernandez, Leyden" uniqKey="Fernandez L" first="Leyden" last="Fernandez">Leyden Fernandez</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mercader, Josep M" sort="Mercader, Josep M" uniqKey="Mercader J" first="Josep M" last="Mercader">Josep M. Mercader</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Planas Felix, Merce" sort="Planas Felix, Merce" uniqKey="Planas Felix M" first="Mercè" last="Planas-Fèlix">Mercè Planas-Fèlix</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Torrents, David" sort="Torrents, David" uniqKey="Torrents D" first="David" last="Torrents">David Torrents</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff42">Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg Luís Companys 23, Barcelona, 08010 USA</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">25294412</idno>
<idno type="pmc">4287501</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4287501</idno>
<idno type="RBID">PMC:4287501</idno>
<idno type="doi">10.1186/1471-2164-15-877</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">000168</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities</title>
<author><name sortKey="Fernandez, Leyden" sort="Fernandez, Leyden" uniqKey="Fernandez L" first="Leyden" last="Fernandez">Leyden Fernandez</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mercader, Josep M" sort="Mercader, Josep M" uniqKey="Mercader J" first="Josep M" last="Mercader">Josep M. Mercader</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Planas Felix, Merce" sort="Planas Felix, Merce" uniqKey="Planas Felix M" first="Mercè" last="Planas-Fèlix">Mercè Planas-Fèlix</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Torrents, David" sort="Torrents, David" uniqKey="Torrents D" first="David" last="Torrents">David Torrents</name>
<affiliation><nlm:aff id="Aff41">Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff42">Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg Luís Companys 23, Barcelona, 08010 USA</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">BMC Genomics</title>
<idno type="eISSN">1471-2164</idno>
<imprint><date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>It has been shown in a number of metagenomic studies that the addition and removal of specific genes have allowed microbiomes to adapt to specific environmental conditions by losing and gaining specific functions. But it is not known whether and how the regulation of gene expression also contributes to adaptation.</p>
</sec>
<sec><title>Results</title>
<p>We have here characterized and analyzed the metaregulome of three different environments, as well as their impact in the adaptation to particular variable physico-chemical conditions. For this, we have developed a computational protocol to extract regulatory regions and their corresponding transcription factors binding sites directly from metagenomic reads and applied it to three well known environments: Acid Mine, Whale Fall, and Waseca Farm. Taking the density of regulatory sites in promoters as a measure of the potential and complexity of gene regulation, we found it to be quantitatively the same in all three environments, despite their different physico-chemical conditions and species composition. However, we found that each environment distributes their regulatory potential differently across their functional space. Among the functions with highest regulatory potential in each niche, we found significant enrichment of processes related to sensing and buffering external variable factors specific to each environment, like for example, the availability of co-factors in deep sea, of oligosaccharides in soil and the regulation of pH in the acid mine.</p>
</sec>
<sec><title>Conclusions</title>
<p>These results highlight the potential impact of gene regulation in the adaptation of bacteria to the different habitats through the distribution of their regulatory potential among specific functions, and point to critical environmental factors that challenge the growth of any microbial community.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/1471-2164-15-877) contains supplementary material, which is available to authorized users.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Tringe, Sg" uniqKey="Tringe S">SG Tringe</name>
</author>
<author><name sortKey="Von Mering, C" uniqKey="Von Mering C">C von Mering</name>
</author>
<author><name sortKey="Kobayashi, A" uniqKey="Kobayashi A">A Kobayashi</name>
</author>
<author><name sortKey="Salamov, Aa" uniqKey="Salamov A">AA Salamov</name>
</author>
<author><name sortKey="Chen, K" uniqKey="Chen K">K Chen</name>
</author>
<author><name sortKey="Chang, Hw" uniqKey="Chang H">HW Chang</name>
</author>
<author><name sortKey="Podar, M" uniqKey="Podar M">M Podar</name>
</author>
<author><name sortKey="Short, Jm" uniqKey="Short J">JM Short</name>
</author>
<author><name sortKey="Mathur, Ej" uniqKey="Mathur E">EJ Mathur</name>
</author>
<author><name sortKey="Detter, Jc" uniqKey="Detter J">JC Detter</name>
</author>
<author><name sortKey="Bork, P" uniqKey="Bork P">P Bork</name>
</author>
<author><name sortKey="Hugenholtz, P" uniqKey="Hugenholtz P">P Hugenholtz</name>
</author>
<author><name sortKey="Rubin, Em" uniqKey="Rubin E">EM Rubin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gianoulis, Ta" uniqKey="Gianoulis T">TA Gianoulis</name>
</author>
<author><name sortKey="Raes, J" uniqKey="Raes J">J Raes</name>
</author>
<author><name sortKey="Patel, Pv" uniqKey="Patel P">PV Patel</name>
</author>
<author><name sortKey="Bjornson, R" uniqKey="Bjornson R">R Bjornson</name>
</author>
<author><name sortKey="Korbel, Jo" uniqKey="Korbel J">JO Korbel</name>
</author>
<author><name sortKey="Letunic, I" uniqKey="Letunic I">I Letunic</name>
</author>
<author><name sortKey="Yamada, T" uniqKey="Yamada T">T Yamada</name>
</author>
<author><name sortKey="Paccanaro, A" uniqKey="Paccanaro A">A Paccanaro</name>
</author>
<author><name sortKey="Jensen, Lj" uniqKey="Jensen L">LJ Jensen</name>
</author>
<author><name sortKey="Snyder, M" uniqKey="Snyder M">M Snyder</name>
</author>
<author><name sortKey="Bork, P" uniqKey="Bork P">P Bork</name>
</author>
<author><name sortKey="Gerstein, Mb" uniqKey="Gerstein M">MB Gerstein</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Poretsky, Rs" uniqKey="Poretsky R">RS Poretsky</name>
</author>
<author><name sortKey="Bano, N" uniqKey="Bano N">N Bano</name>
</author>
<author><name sortKey="Buchan, A" uniqKey="Buchan A">A Buchan</name>
</author>
<author><name sortKey="Lecleir, G" uniqKey="Lecleir G">G LeCleir</name>
</author>
<author><name sortKey="Kleikemper, J" uniqKey="Kleikemper J">J Kleikemper</name>
</author>
<author><name sortKey="Pickering, M" uniqKey="Pickering M">M Pickering</name>
</author>
<author><name sortKey="Pate, Wm" uniqKey="Pate W">WM Pate</name>
</author>
<author><name sortKey="Moran, Ma" uniqKey="Moran M">MA Moran</name>
</author>
<author><name sortKey="Hollibaugh, Jt" uniqKey="Hollibaugh J">JT Hollibaugh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Poretsky, Rs" uniqKey="Poretsky R">RS Poretsky</name>
</author>
<author><name sortKey="Hewson, I" uniqKey="Hewson I">I Hewson</name>
</author>
<author><name sortKey="Sun, S" uniqKey="Sun S">S Sun</name>
</author>
<author><name sortKey="Allen, Ae" uniqKey="Allen A">AE Allen</name>
</author>
<author><name sortKey="Zehr, Jp" uniqKey="Zehr J">JP Zehr</name>
</author>
<author><name sortKey="Moran, Ma" uniqKey="Moran M">MA Moran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gifford, Sm" uniqKey="Gifford S">SM Gifford</name>
</author>
<author><name sortKey="Sharma, S" uniqKey="Sharma S">S Sharma</name>
</author>
<author><name sortKey="Rinta Kanto, Jm" uniqKey="Rinta Kanto J">JM Rinta-Kanto</name>
</author>
<author><name sortKey="Moran, Ma" uniqKey="Moran M">MA Moran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Farre, D" uniqKey="Farre D">D Farre</name>
</author>
<author><name sortKey="Bellora, N" uniqKey="Bellora N">N Bellora</name>
</author>
<author><name sortKey="Mularoni, L" uniqKey="Mularoni L">L Mularoni</name>
</author>
<author><name sortKey="Messeguer, X" uniqKey="Messeguer X">X Messeguer</name>
</author>
<author><name sortKey="Alba, Mm" uniqKey="Alba M">MM Alba</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lin, Z" uniqKey="Lin Z">Z Lin</name>
</author>
<author><name sortKey="Wu, Ws" uniqKey="Wu W">WS Wu</name>
</author>
<author><name sortKey="Liang, H" uniqKey="Liang H">H Liang</name>
</author>
<author><name sortKey="Woo, Y" uniqKey="Woo Y">Y Woo</name>
</author>
<author><name sortKey="Li, Wh" uniqKey="Li W">WH Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Merino, E" uniqKey="Merino E">E Merino</name>
</author>
<author><name sortKey="Jensen, Ra" uniqKey="Jensen R">RA Jensen</name>
</author>
<author><name sortKey="Yanofsky, C" uniqKey="Yanofsky C">C Yanofsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
<author><name sortKey="Rhodius, V" uniqKey="Rhodius V">V Rhodius</name>
</author>
<author><name sortKey="Gross, C" uniqKey="Gross C">C Gross</name>
</author>
<author><name sortKey="Siggia, Ed" uniqKey="Siggia E">ED Siggia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mccue, L" uniqKey="Mccue L">L McCue</name>
</author>
<author><name sortKey="Thompson, W" uniqKey="Thompson W">W Thompson</name>
</author>
<author><name sortKey="Carmack, C" uniqKey="Carmack C">C Carmack</name>
</author>
<author><name sortKey="Ryan, Mp" uniqKey="Ryan M">MP Ryan</name>
</author>
<author><name sortKey="Liu, Js" uniqKey="Liu J">JS Liu</name>
</author>
<author><name sortKey="Derbyshire, V" uniqKey="Derbyshire V">V Derbyshire</name>
</author>
<author><name sortKey="Lawrence, Ce" uniqKey="Lawrence C">CE Lawrence</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huffman, Jl" uniqKey="Huffman J">JL Huffman</name>
</author>
<author><name sortKey="Brennan, Rg" uniqKey="Brennan R">RG Brennan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rodionov, Da" uniqKey="Rodionov D">DA Rodionov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gelfand, Ms" uniqKey="Gelfand M">MS Gelfand</name>
</author>
<author><name sortKey="Koonin, Ev" uniqKey="Koonin E">EV Koonin</name>
</author>
<author><name sortKey="Mironov, Aa" uniqKey="Mironov A">AA Mironov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Laing, E" uniqKey="Laing E">E Laing</name>
</author>
<author><name sortKey="Sidhu, K" uniqKey="Sidhu K">K Sidhu</name>
</author>
<author><name sortKey="Hubbard, Sj" uniqKey="Hubbard S">SJ Hubbard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iqbal, M" uniqKey="Iqbal M">M Iqbal</name>
</author>
<author><name sortKey="Mast, Y" uniqKey="Mast Y">Y Mast</name>
</author>
<author><name sortKey="Amin, R" uniqKey="Amin R">R Amin</name>
</author>
<author><name sortKey="Hodgson, Da" uniqKey="Hodgson D">DA Hodgson</name>
</author>
<author><name sortKey="Consortium, S" uniqKey="Consortium S">S Consortium</name>
</author>
<author><name sortKey="Wohlleben, W" uniqKey="Wohlleben W">W Wohlleben</name>
</author>
<author><name sortKey="Burroughs, Nj" uniqKey="Burroughs N">NJ Burroughs</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, L" uniqKey="Li L">L Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Liu, J" uniqKey="Liu J">J Liu</name>
</author>
<author><name sortKey="Xu, X" uniqKey="Xu X">X Xu</name>
</author>
<author><name sortKey="Stormo, Gd" uniqKey="Stormo G">GD Stormo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sun, J" uniqKey="Sun J">J Sun</name>
</author>
<author><name sortKey="Tuncay, K" uniqKey="Tuncay K">K Tuncay</name>
</author>
<author><name sortKey="Haidar, Aa" uniqKey="Haidar A">AA Haidar</name>
</author>
<author><name sortKey="Ensman, L" uniqKey="Ensman L">L Ensman</name>
</author>
<author><name sortKey="Stanley, F" uniqKey="Stanley F">F Stanley</name>
</author>
<author><name sortKey="Trelinski, M" uniqKey="Trelinski M">M Trelinski</name>
</author>
<author><name sortKey="Ortoleva, P" uniqKey="Ortoleva P">P Ortoleva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Leuze, Mr" uniqKey="Leuze M">MR Leuze</name>
</author>
<author><name sortKey="Karpinets, Tv" uniqKey="Karpinets T">TV Karpinets</name>
</author>
<author><name sortKey="Syed, Mh" uniqKey="Syed M">MH Syed</name>
</author>
<author><name sortKey="Beliaev, As" uniqKey="Beliaev A">AS Beliaev</name>
</author>
<author><name sortKey="Uberbacher, Ec" uniqKey="Uberbacher E">EC Uberbacher</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, S" uniqKey="Zhang S">S Zhang</name>
</author>
<author><name sortKey="Li, S" uniqKey="Li S">S Li</name>
</author>
<author><name sortKey="Niu, M" uniqKey="Niu M">M Niu</name>
</author>
<author><name sortKey="Pham, Pt" uniqKey="Pham P">PT Pham</name>
</author>
<author><name sortKey="Su, Z" uniqKey="Su Z">Z Su</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Novichkov, Ps" uniqKey="Novichkov P">PS Novichkov</name>
</author>
<author><name sortKey="Laikova, On" uniqKey="Laikova O">ON Laikova</name>
</author>
<author><name sortKey="Novichkova, Es" uniqKey="Novichkova E">ES Novichkova</name>
</author>
<author><name sortKey="Gelfand, Ms" uniqKey="Gelfand M">MS Gelfand</name>
</author>
<author><name sortKey="Arkin, Ap" uniqKey="Arkin A">AP Arkin</name>
</author>
<author><name sortKey="Dubchak, I" uniqKey="Dubchak I">I Dubchak</name>
</author>
<author><name sortKey="Rodionov, Da" uniqKey="Rodionov D">DA Rodionov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grissa, I" uniqKey="Grissa I">I Grissa</name>
</author>
<author><name sortKey="Vergnaud, G" uniqKey="Vergnaud G">G Vergnaud</name>
</author>
<author><name sortKey="Pourcel, C" uniqKey="Pourcel C">C Pourcel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grissa, I" uniqKey="Grissa I">I Grissa</name>
</author>
<author><name sortKey="Vergnaud, G" uniqKey="Vergnaud G">G Vergnaud</name>
</author>
<author><name sortKey="Pourcel, C" uniqKey="Pourcel C">C Pourcel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Overbeek, R" uniqKey="Overbeek R">R Overbeek</name>
</author>
<author><name sortKey="Begley, T" uniqKey="Begley T">T Begley</name>
</author>
<author><name sortKey="Butler, Rm" uniqKey="Butler R">RM Butler</name>
</author>
<author><name sortKey="Choudhuri, Jv" uniqKey="Choudhuri J">JV Choudhuri</name>
</author>
<author><name sortKey="Chuang, Hy" uniqKey="Chuang H">HY Chuang</name>
</author>
<author><name sortKey="Cohoon, M" uniqKey="Cohoon M">M Cohoon</name>
</author>
<author><name sortKey="De Crecy Lagard, V" uniqKey="De Crecy Lagard V">V de Crecy-Lagard</name>
</author>
<author><name sortKey="Diaz, N" uniqKey="Diaz N">N Diaz</name>
</author>
<author><name sortKey="Disz, T" uniqKey="Disz T">T Disz</name>
</author>
<author><name sortKey="Edwards, R" uniqKey="Edwards R">R Edwards</name>
</author>
<author><name sortKey="Fonstein, M" uniqKey="Fonstein M">M Fonstein</name>
</author>
<author><name sortKey="Frank, Ed" uniqKey="Frank E">ED Frank</name>
</author>
<author><name sortKey="Gerdes, S" uniqKey="Gerdes S">S Gerdes</name>
</author>
<author><name sortKey="Glass, Em" uniqKey="Glass E">EM Glass</name>
</author>
<author><name sortKey="Goesmann, A" uniqKey="Goesmann A">A Goesmann</name>
</author>
<author><name sortKey="Hanson, A" uniqKey="Hanson A">A Hanson</name>
</author>
<author><name sortKey="Iwata Reuyl, D" uniqKey="Iwata Reuyl D">D Iwata-Reuyl</name>
</author>
<author><name sortKey="Jensen, R" uniqKey="Jensen R">R Jensen</name>
</author>
<author><name sortKey="Jamshidi, N" uniqKey="Jamshidi N">N Jamshidi</name>
</author>
<author><name sortKey="Krause, L" uniqKey="Krause L">L Krause</name>
</author>
<author><name sortKey="Kubal, M" uniqKey="Kubal M">M Kubal</name>
</author>
<author><name sortKey="Larsen, N" uniqKey="Larsen N">N Larsen</name>
</author>
<author><name sortKey="Linke, B" uniqKey="Linke B">B Linke</name>
</author>
<author><name sortKey="Mchardy, Ac" uniqKey="Mchardy A">AC McHardy</name>
</author>
<author><name sortKey="Meyer, F" uniqKey="Meyer F">F Meyer</name>
</author>
<author><name sortKey="Neuweger, H" uniqKey="Neuweger H">H Neuweger</name>
</author>
<author><name sortKey="Olsen, G" uniqKey="Olsen G">G Olsen</name>
</author>
<author><name sortKey="Olson, R" uniqKey="Olson R">R Olson</name>
</author>
<author><name sortKey="Osterman, A" uniqKey="Osterman A">A Osterman</name>
</author>
<author><name sortKey="Portnoy, V" uniqKey="Portnoy V">V Portnoy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huson, Dh" uniqKey="Huson D">DH Huson</name>
</author>
<author><name sortKey="Auch, Af" uniqKey="Auch A">AF Auch</name>
</author>
<author><name sortKey="Qi, J" uniqKey="Qi J">J Qi</name>
</author>
<author><name sortKey="Schuster, Sc" uniqKey="Schuster S">SC Schuster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chin, Cs" uniqKey="Chin C">CS Chin</name>
</author>
<author><name sortKey="Chuang, Jh" uniqKey="Chuang J">JH Chuang</name>
</author>
<author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Abele, D" uniqKey="Abele D">D Abele</name>
</author>
<author><name sortKey="Puntarulo, S" uniqKey="Puntarulo S">S Puntarulo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Storz, G" uniqKey="Storz G">G Storz</name>
</author>
<author><name sortKey="Imlay, Ja" uniqKey="Imlay J">JA Imlay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edwards, Kj" uniqKey="Edwards K">KJ Edwards</name>
</author>
<author><name sortKey="Gihring, Tm" uniqKey="Gihring T">TM Gihring</name>
</author>
<author><name sortKey="Banfield, Jf" uniqKey="Banfield J">JF Banfield</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Albers, S V" uniqKey="Albers S">S-V Albers</name>
</author>
<author><name sortKey="Koning, Sm" uniqKey="Koning S">SM Koning</name>
</author>
<author><name sortKey="Konings, Wn" uniqKey="Konings W">WN Konings</name>
</author>
<author><name sortKey="Driessen, Ajm" uniqKey="Driessen A">AJM Driessen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kempf, B" uniqKey="Kempf B">B Kempf</name>
</author>
<author><name sortKey="Bremer, E" uniqKey="Bremer E">E Bremer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Osorio, H" uniqKey="Osorio H">H Osorio</name>
</author>
<author><name sortKey="Martinez, V" uniqKey="Martinez V">V Martinez</name>
</author>
<author><name sortKey="Nieto, Pa" uniqKey="Nieto P">PA Nieto</name>
</author>
<author><name sortKey="Holmes, Ds" uniqKey="Holmes D">DS Holmes</name>
</author>
<author><name sortKey="Quatrini, R" uniqKey="Quatrini R">R Quatrini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schrenk, Mo" uniqKey="Schrenk M">MO Schrenk</name>
</author>
<author><name sortKey="Edwards, Kj" uniqKey="Edwards K">KJ Edwards</name>
</author>
<author><name sortKey="Goodman, Rm" uniqKey="Goodman R">RM Goodman</name>
</author>
<author><name sortKey="Hamers, Rj" uniqKey="Hamers R">RJ Hamers</name>
</author>
<author><name sortKey="Banfield, Jf" uniqKey="Banfield J">JF Banfield</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Guazzaroni, Me" uniqKey="Guazzaroni M">ME Guazzaroni</name>
</author>
<author><name sortKey="Morgante, V" uniqKey="Morgante V">V Morgante</name>
</author>
<author><name sortKey="Mirete, S" uniqKey="Mirete S">S Mirete</name>
</author>
<author><name sortKey="Gonzalez Pastor, Je" uniqKey="Gonzalez Pastor J">JE Gonzalez-Pastor</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sun, S" uniqKey="Sun S">S Sun</name>
</author>
<author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name>
</author>
<author><name sortKey="Li, W" uniqKey="Li W">W Li</name>
</author>
<author><name sortKey="Altintas, I" uniqKey="Altintas I">I Altintas</name>
</author>
<author><name sortKey="Lin, A" uniqKey="Lin A">A Lin</name>
</author>
<author><name sortKey="Peltier, S" uniqKey="Peltier S">S Peltier</name>
</author>
<author><name sortKey="Stocks, K" uniqKey="Stocks K">K Stocks</name>
</author>
<author><name sortKey="Allen, Ee" uniqKey="Allen E">EE Allen</name>
</author>
<author><name sortKey="Ellisman, M" uniqKey="Ellisman M">M Ellisman</name>
</author>
<author><name sortKey="Grethe, J" uniqKey="Grethe J">J Grethe</name>
</author>
<author><name sortKey="Wooley, J" uniqKey="Wooley J">J Wooley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name>
</author>
<author><name sortKey="Gish, W" uniqKey="Gish W">W Gish</name>
</author>
<author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name>
</author>
<author><name sortKey="Myers, Ew" uniqKey="Myers E">EW Myers</name>
</author>
<author><name sortKey="Lipman, Dj" uniqKey="Lipman D">DJ Lipman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Giovannoni, Sj" uniqKey="Giovannoni S">SJ Giovannoni</name>
</author>
<author><name sortKey="Tripp, Hj" uniqKey="Tripp H">HJ Tripp</name>
</author>
<author><name sortKey="Givan, S" uniqKey="Givan S">S Givan</name>
</author>
<author><name sortKey="Podar, M" uniqKey="Podar M">M Podar</name>
</author>
<author><name sortKey="Vergin, Kl" uniqKey="Vergin K">KL Vergin</name>
</author>
<author><name sortKey="Baptista, D" uniqKey="Baptista D">D Baptista</name>
</author>
<author><name sortKey="Bibbs, L" uniqKey="Bibbs L">L Bibbs</name>
</author>
<author><name sortKey="Eads, J" uniqKey="Eads J">J Eads</name>
</author>
<author><name sortKey="Richardson, Th" uniqKey="Richardson T">TH Richardson</name>
</author>
<author><name sortKey="Noordewier, M" uniqKey="Noordewier M">M Noordewier</name>
</author>
<author><name sortKey="Rappe, Ms" uniqKey="Rappe M">MS Rappé</name>
</author>
<author><name sortKey="Short, Jm" uniqKey="Short J">JM Short</name>
</author>
<author><name sortKey="Carrington, Jc" uniqKey="Carrington J">JC Carrington</name>
</author>
<author><name sortKey="Mathur, Ej" uniqKey="Mathur E">EJ Mathur</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burge, Sw" uniqKey="Burge S">SW Burge</name>
</author>
<author><name sortKey="Daub, J" uniqKey="Daub J">J Daub</name>
</author>
<author><name sortKey="Eberhardt, R" uniqKey="Eberhardt R">R Eberhardt</name>
</author>
<author><name sortKey="Tate, J" uniqKey="Tate J">J Tate</name>
</author>
<author><name sortKey="Barquist, L" uniqKey="Barquist L">L Barquist</name>
</author>
<author><name sortKey="Nawrocki, Ep" uniqKey="Nawrocki E">EP Nawrocki</name>
</author>
<author><name sortKey="Eddy, Sr" uniqKey="Eddy S">SR Eddy</name>
</author>
<author><name sortKey="Gardner, Pp" uniqKey="Gardner P">PP Gardner</name>
</author>
<author><name sortKey="Bateman, A" uniqKey="Bateman A">A Bateman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hyatt, D" uniqKey="Hyatt D">D Hyatt</name>
</author>
<author><name sortKey="Chen, Gl" uniqKey="Chen G">GL Chen</name>
</author>
<author><name sortKey="Locascio, Pf" uniqKey="Locascio P">PF LoCascio</name>
</author>
<author><name sortKey="Land, Ml" uniqKey="Land M">ML Land</name>
</author>
<author><name sortKey="Larimer, Fw" uniqKey="Larimer F">FW Larimer</name>
</author>
<author><name sortKey="Hauser, Lj" uniqKey="Hauser L">LJ Hauser</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">25294412</article-id>
<article-id pub-id-type="pmc">4287501</article-id>
<article-id pub-id-type="publisher-id">6783</article-id>
<article-id pub-id-type="doi">10.1186/1471-2164-15-877</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Fernandez</surname>
<given-names>Leyden</given-names>
</name>
<address><email>leyden.fernandez@bsc.es</email>
</address>
<xref ref-type="aff" rid="Aff41"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Mercader</surname>
<given-names>Josep M</given-names>
</name>
<address><email>josep.mercader@bsc.es</email>
</address>
<xref ref-type="aff" rid="Aff41"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Planas-Fèlix</surname>
<given-names>Mercè</given-names>
</name>
<address><email>merce.planas@bsc.es</email>
</address>
<xref ref-type="aff" rid="Aff41"></xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><name><surname>Torrents</surname>
<given-names>David</given-names>
</name>
<address><email>david.torrents@bsc.es</email>
</address>
<xref ref-type="aff" rid="Aff41"></xref>
<xref ref-type="aff" rid="Aff42"></xref>
</contrib>
<aff id="Aff41"><label></label>
Joint IRB-BSC program on Computational Biology. BSC, Jordi Girona, 29, 08034 Barcelona, Spain</aff>
<aff id="Aff42"><label></label>
Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg Luís Companys 23, Barcelona, 08010 USA</aff>
</contrib-group>
<pub-date pub-type="epub"><day>8</day>
<month>10</month>
<year>2014</year>
</pub-date>
<pub-date pub-type="collection"><year>2014</year>
</pub-date>
<volume>15</volume>
<issue>1</issue>
<elocation-id>877</elocation-id>
<history><date date-type="received"><day>4</day>
<month>2</month>
<year>2014</year>
</date>
<date date-type="accepted"><day>24</day>
<month>9</month>
<year>2014</year>
</date>
</history>
<permissions><copyright-statement>© Fernandez et al.; licensee BioMed Central Ltd. 2014</copyright-statement>
<license license-type="open-access"><license-p>This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution 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 the original work is properly credited. 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>It has been shown in a number of metagenomic studies that the addition and removal of specific genes have allowed microbiomes to adapt to specific environmental conditions by losing and gaining specific functions. But it is not known whether and how the regulation of gene expression also contributes to adaptation.</p>
</sec>
<sec><title>Results</title>
<p>We have here characterized and analyzed the metaregulome of three different environments, as well as their impact in the adaptation to particular variable physico-chemical conditions. For this, we have developed a computational protocol to extract regulatory regions and their corresponding transcription factors binding sites directly from metagenomic reads and applied it to three well known environments: Acid Mine, Whale Fall, and Waseca Farm. Taking the density of regulatory sites in promoters as a measure of the potential and complexity of gene regulation, we found it to be quantitatively the same in all three environments, despite their different physico-chemical conditions and species composition. However, we found that each environment distributes their regulatory potential differently across their functional space. Among the functions with highest regulatory potential in each niche, we found significant enrichment of processes related to sensing and buffering external variable factors specific to each environment, like for example, the availability of co-factors in deep sea, of oligosaccharides in soil and the regulation of pH in the acid mine.</p>
</sec>
<sec><title>Conclusions</title>
<p>These results highlight the potential impact of gene regulation in the adaptation of bacteria to the different habitats through the distribution of their regulatory potential among specific functions, and point to critical environmental factors that challenge the growth of any microbial community.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/1471-2164-15-877) contains supplementary material, which is available to authorized users.</p>
</sec>
</abstract>
<kwd-group xml:lang="en"><title>Keywords</title>
<kwd>Adaptation</kwd>
<kwd>Environment</kwd>
<kwd>Gene regulation</kwd>
<kwd>Metagenomes</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2014</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1"><title>Background</title>
<p>Metagenomic studies generate a massive amount of sequence information of communities of organisms living in different physicochemical conditions. This allows, for the first time, to search for the molecular and genetic basis of adaptation through the comparison and the study of genomes of different species sharing the same environment, and of similar species living in different conditions. The comparative studies of the potential protein content in many of these datasets have already provided interesting examples of specific functions that correlate with specific characteristics of the environment. For example, in the search of functional fingerprints related to specific habitats, a comparative analysis between soil, and deep and superficial aquatic environments found abundant orthologous groups specific of these particular habitats
[<xref ref-type="bibr" rid="CR1">1</xref>
]. In this case, the examination of higher order processes reveals differences in energy production between these three niches, such as starch and sucrose metabolism in soil or photosynthesis in oligotrophic surface waters
[<xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR2">2</xref>
].</p>
<p>More recently, metagenomic studies have gone beyond the sequencing of DNA and the counting of genes, and have incorporated techniques and protocols to detect, measure and analyze their transcriptome. While the sequencing of metagenomes provides an overview of the genes present in specific environments that can potentially play a role in adaptation, the analysis of expression provides a more precise picture of what functions are expressed and active in a particular moment of the environment. Even though the techniques for mRNA isolation and sequencing from metagenomic samples are still not able to provide comprehensive pictures of expression profiles, there have been important progresses in this direction and some interesting findings. For example, one of the first studies of metatranscriptome, despite it covered a small fraction of the expressed genes, identified specific biological processes active in bacterioplankton communities that could be correlated with either marine or freshwater conditions
[<xref ref-type="bibr" rid="CR3">3</xref>
]. As the coverage and accuracy of these analysis increased (mostly by including next generation sequencing techniques), more active processes have been linked to variable environmental conditions. For instance, an expression time-series performed on microbial communities living in surface oceanic showed that processes of energy production were active in hours with light, while anabolic housekeeping processes were predominant during the night
[<xref ref-type="bibr" rid="CR4">4</xref>
]. Despite the underlying methodology behind, metatranscriptomics still needs to overcome several challenges
[<xref ref-type="bibr" rid="CR5">5</xref>
]. But the rapid progress in this field is promising and we will soon have the opportunity of building accurate expression profiles and compare them across environments, as well as exploring the interaction of processes of different organisms within specific environments.</p>
<p>In the present study we have conducted a novel approach that complements and bridges metagenomic and metranscriptomic concepts. The rationale behind this study relies on the hypothesis that the regulation of the expression of those biological functions that confer adaptation to variable environmental conditions will show higher complexity, i.e. they will have complex regulatory regions.</p>
<p>Previous studies
[<xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR7">7</xref>
] have shown that genes with complex regulation requirements show higher number of transcription factor binding sites (TFBSs) in their upstream cis-regulatory regions compared to housekeeping genes. For example, stress-response genes in yeast need a precise regulation of their expressions patterns to adapt to drastic changes of environmental conditions and also show a significantly higher number of different TFBSs in their upstream regulatory regions. Beyond the extensive analysis of the regulatory characteristics of particular functions
[<xref ref-type="bibr" rid="CR8">8</xref>
], up to now, there are not global approaches and studies on how the regulatory potential of entire microbial communities is influenced and organized in natural habitats.</p>
<p>In particular, and using the same rationale, we have measured and compared the complexity of gene regulation in bacteria and archaea living in environments with distinct underlying physico-chemical conditions. For that purpose, we searched within each of the environments for specific functional signatures predicted to have high regulatory potential. These are correlated with specific and also dynamic physico-chemical stress factors of each of the niches. The functional significance of the differences detected highlights the existence of adaptation strategies that rely on the regulatory potential of regions that control the expression of specific fitness genes.</p>
</sec>
<sec id="Sec2"><title>Results and discussion</title>
<p>With the ultimate goal of identifying and characterizing the extend, to which environmental factors influence the organization of the regulatory potential of particular microbial communities, we have studied and compared the regulome of three fundamentally different ecological niches using whole metagenomic data. We next provide details on the major results and findings of this study: (1) The development of a new pipeline for the identification and prediction of proximal regulatory regions and their TFBS from metagenomic data; (2) and the generation of a collection of regulatory regions from three well studied and reference metagenomic samples (Whale Fall, Waseca Farm and Acid Mine). The comparative analysis of this data has shown that, while (3) the overall distribution of TFBS on promoters is the same across environments, their distribution across their functional space is significantly different, as (4) promoters with higher number of TFBS tend to regulate environment specific functions, and (5) a fraction of these are environment specific and can be linked to characteristic external physicochemical factors (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S1).</p>
<sec id="Sec3"><title>Identification and classification of proximal regulatory regions from metagenomic data</title>
<p>We first characterized and analyzed the gene regulatory space from metagenomic data obtained from three well-known sequenced environments with clearly different physico-chemical properties: Whale Fall Community, Acid Mine and Waseca County Farm Soil
[<xref ref-type="bibr" rid="CR1">1</xref>
]. For that, we started by identifying and defining gene regulatory regions to later characterize them, as to their levels of TF binding, i.e. their regulatory potential. For the design of a search strategy, we followed two major considerations: first (1) avoiding biases in favor of most abundant and well-known bacteria (and closely related species), as well as, (2) ensuring an equal coverage through all the phyla detectable in those samples. As a result, we developed a pipeline that consists of two major steps: (1) first the identification of proximal regulatory regions and then, within each of them, (2) the prediction of potential regulatory transcription factor binding sites. The complete pipeline is detailed in the Methods and summarized graphically in Figure
<xref rid="Fig1" ref-type="fig">1</xref>
.<fig id="Fig1"><label>Figure 1</label>
<caption><p><bold>Computational pipeline generated for the identification and characterization of proximal regulatory regions from metagenomic data.</bold>
The figure summarizes the methods used and results obtained in each step of the study.</p>
</caption>
<graphic xlink:href="12864_2014_6783_Fig1_HTML" id="d30e312"></graphic>
</fig>
</p>
<p>Through extensive homology searches, our procedure identified putative proximal regulatory regions in Waseca Farm Soil, in Acid Mine Drainage, and in the Whale Falls Sample (a complete catalogue of these regions can be found in Additional files
<xref rid="MOESM2" ref-type="media">2</xref>
,
<xref rid="MOESM3" ref-type="media">3</xref>
and
<xref rid="MOESM4" ref-type="media">4</xref>
). A first and basic taxonomical analysis of these sequences shows that these promoters cover all phyla (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S2) that were previously described in these environments
[<xref ref-type="bibr" rid="CR1">1</xref>
].</p>
<p>Next, we estimated the level of regulatory potential for each of these promoter regions through the prediction of their transcription factor binding sites. In order to minimize possible biases favoring promoters from well-studied bacteria (or close species), we did not consider TFBSs prediction strategies that rely on the homology mapping of described TFBSs. Instead, we used a <italic>de novo</italic>
prediction protocol that relies on the identification of palindromic repeats
[<xref ref-type="bibr" rid="CR9">9</xref>
], which have been previously determined as preferred binding sites for transcription regulators in bacteria and archaea
[<xref ref-type="bibr" rid="CR9">9</xref>
–<xref ref-type="bibr" rid="CR14">14</xref>
]. Because this method was originally developed for the analysis of single genomes
[<xref ref-type="bibr" rid="CR9">9</xref>
] and, although it has been applied to a wide variety of bacterial sequences and studies
[<xref ref-type="bibr" rid="CR14">14</xref>
–<xref ref-type="bibr" rid="CR16">16</xref>
], we needed to adapt it to cope with the heterogeneity and redundancy of metagenomic data by including some modification in the scoring system.</p>
</sec>
<sec id="Sec4"><title>Evaluation of predicted promoters and TFBSs</title>
<p>Like any other <italic>de novo</italic>
prediction method in sequence analysis, we have to initially assume the presence of false positive TFBS models among correct predictions. To assess for the reliability of all of our predictions and to put our strategy and results into the context of our goals and of the current knowledge about regulatory regions in prokaryotes, we performed different quantitative and qualitative comparisons with available independent data and methodologies.</p>
<p>From a quantitative point of view, we (1) first observed that the global average of 10 TFBS per promoter (with 0 as minimum and 25 as maximum values) that we identify from all three environments is in agreement with previous estimates obtained with different bacterial species and methodologies. For example, using genome comparative analysis, an average of 11–13 TFBS motifs per promoter was found for <italic>Shewanella</italic>
[<xref ref-type="bibr" rid="CR17">17</xref>
]. In addition, a study of the transcription regulatory network of <italic>E. Coli K12</italic>
predicted up to 16 sites per promoter
[<xref ref-type="bibr" rid="CR18">18</xref>
], and up to 20 through the identification of half-sites motifs
[<xref ref-type="bibr" rid="CR19">19</xref>
]. (2) We also evaluated the performance of our methodology by comparing our results with those obtained with an independent method, MotifClick, that predicts cis-regulatory regions using a graph-based polynomial-time algorithm
[<xref ref-type="bibr" rid="CR20">20</xref>
]. After running both predictors over intergenic <italic>E. Coli</italic>
regions, we observed that the densities of TFBS resulting from one or the other strategy showed high correlation values (rho = 0.52, p-value < 2.2 × 10<sup>−6</sup>
; (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S3).</p>
<p>From a qualitative point of view, we first (1) assessed the biological significance of our predictions by carrying out a randomization test consisting in applying the same prediction pipeline to our collection of promoters with their nucleotide sequence completely shuffled, i.e. with no biological information. We observed that the distributions of the number of motifs per promoter were significantly different between the real and the randomized sample (Figure
<xref rid="Fig2" ref-type="fig">2</xref>
). (2) Furthermore, we screened for coincidences between our predicted TFBSs and those reported in the RegPrecise database
[<xref ref-type="bibr" rid="CR21">21</xref>
], which consist on manually curated site reconstructions in various bacteria genomes. This comparison showed that 28% of our predicted binding sites include, at least, one possible binding sequence of the matrices for each of the 38 TFs included in RegPrecise (Additional file
<xref rid="MOESM5" ref-type="media">5</xref>
). (3) Finally, we also searched for a particular type of false predictions, which consist on regulatory palindromic repeats with no binding potential, named Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)
[<xref ref-type="bibr" rid="CR22">22</xref>
]. The results that we obtained using the CRISPRFinder web tool
[<xref ref-type="bibr" rid="CR23">23</xref>
] showed a negligible amount of these regions (less than 1% of our set of promoters), which were subsequently removed from the analysis.<fig id="Fig2"><label>Figure 2</label>
<caption><p><bold>Distribution of regulatory regions predicted per each environments analyzed in this study according to their content of TFBSs.</bold>
Solid lines represent the distributions of TFBS obtained using real data, whereas dashed lines show predictions using randomized DNA sequences, i.e. with no biological information. Three major groups of regulatory potential are also shown: low, medium and high, which corresponds to the three bins analyzed in Additional file <xref rid="MOESM4" ref-type="media">1</xref>
: Figure S4.</p>
</caption>
<graphic xlink:href="12864_2014_6783_Fig2_HTML" id="d30e424"></graphic>
</fig>
</p>
<p>In summary, all these evaluation tests suggest that our set of promoters is both, quantitatively and qualitatively reliable, as they show a significant fraction of reported TFBS, and a small portion of false positives. But, most importantly, the presence of this small fraction of false positives is not expected to affect our final conclusions, as these come from comparisons within and between environments and do not rely on absolute TFBSs counts.</p>
</sec>
<sec id="Sec5"><title>Functional organization of regulatory potentials within each environment</title>
<p>We then studied how microbial communities living in these environments organize and distribute their regulatory potential through the different biological functions and to which extend this could obey to specific adaptation needs. It is interesting to observe that, whereas the range of density of predicted sites per promoter is wide within each of the environments, the overall distribution and the averages are similar: 9.98 (±3.29), 9.58 (±3.49) and 10.28 (±3.35) for Acid Mine Drainage, Waseca Farm Soil and Whale Falls samples, respectively (Figure
<xref rid="Fig2" ref-type="fig">2</xref>
). This indicates that, although these three environments present (1) different sequence coverage, (2) different physicochemical characteristics and (3) different species composition, the overall regulatory potential, as to the total number of different TFBS, and their distribution across the promoters follow a similar pattern.</p>
<p>To go beyond simple counts and to explore whether or not this regulatory potential is distributed equally through all the functions of each of the metagenomes, we first identified the functions under the control of our collection of proximal regulatory regions. For this, we assigned to each promoter the functional category (from SEED database)
[<xref ref-type="bibr" rid="CR24">24</xref>
] of the corresponding downstream coding region using MEGAN
[<xref ref-type="bibr" rid="CR25">25</xref>
] see (Additional files
<xref rid="MOESM6" ref-type="media">6</xref>
,
<xref rid="MOESM7" ref-type="media">7</xref>
and
<xref rid="MOESM8" ref-type="media">8</xref>
, for a complete list of functions and TFBS densities). We first investigated whether the regulatory potential is organized differently over the functional space of each of the environments. For this, we ranked all promoters of each sample according to their TFBS density and count, for each density group, how many associated functions are specific of that particular environment, or co-occur in one or in the other two samples. This analysis showed significant differences between promoters. Interestingly, the functions under the control of promoters with high number of TFBSs show significantly less co-occurrences between environments, than those regulated by promoters with lower regulatory potential. The fact that promoters with high density of TFBS are enriched in environment specific functions provides the first hint that processes that require complex regulation might provide adaptation to environment specific variable external factors. (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S4). We expect that a large fraction of functions that showed a higher co-occurrence among environments likely correspond to housekeeping roles.</p>
<p>To study this further, we next investigated which functions are specifically enriched among the highly regulated ones in each of the environments. For this, we zoomed into the fraction of the 33% highest regulated promoters (i.e. with more than 12 TFBSs/promoter) and subdivided it further into subgroups covering the 1, 5, 10, 20, 30 and 40 top percentages of TFBS density, to finally analyze the functional enrichments within each of them. This analysis highlighted different enriched functions in each of the environments (see Additional file
<xref rid="MOESM4" ref-type="media">1</xref>
: Figure S5 (Acid Mine), S6 (Waseca Soil), S7 (Whale Falls)). These enriched functions cover different types of processes, the majority of them involved in sensing and buffering external factors, such as, receptors and transporters in Acid Mine and stress response systems in Whale Falls.</p>
</sec>
<sec id="Sec6"><title>Potential environment-gene regulation relationships</title>
<p>In order to finally highlight potential points of interaction between highly regulated functions that could provide adaptation to variable conditions specific to each of the environments, we first selected for each habitat, those functions that show stronger enrichment, i.e. with pvalue < 0.05, among the top 1, 5, 10 and 20% groups and with clear orthologous functions in the other two samples. This subgroup of functions include (virulence, cell cycle, carbohydrates metabolism, stress response and cofactors metabolism), which we then compared among environments and evaluated their relationship with the niche specific variable factors. For this, we carried out extensive literature searches on different biochemical mechanisms of adaptation guided by these functions and the characteristics of the environment. Despite the limited information about the environment physico-chemical factors characteristic of available metagenomic studies, we propose in the following sections potential adaptive scenarios by correlating highly regulated functions with known variable external factors in each of the environments.</p>
</sec>
<sec id="Sec7"><title>Waseca Farm Soil</title>
<p>In Waseca Farm Soil, carbohydrates metabolism related functions appear as highly regulated, more precisely di and oligosaccharides metabolism (pvalue = 1×10<sup>−16</sup>
, within environment and adjusted pvalue (Bonferroni) = 9 × 10<sup>−13</sup>
for Fisher’s exact test between groups). This fact could be in concordance with the fluctuations in organic matter concentrations in the soil, such as, plant debris, which has also been previously proposed as an explanation for the presence of other carbohydrate metabolism functions specific of this environment
[<xref ref-type="bibr" rid="CR1">1</xref>
]. This further agrees with the behavior observed in lower eukaryotes abundant in soil, like yeast, where high complexity in their transcriptional regulation were found upstream of genes that play a role in carbohydrates metabolism
[<xref ref-type="bibr" rid="CR26">26</xref>
]; and with the fact that, in this niche, the upstream region of the FruR gene, a known TF that regulates carbohydrate metabolism, appears as highly regulated, with the highest number of predicted TFBS (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S8).</p>
</sec>
<sec id="Sec8"><title>Whale Falls samples</title>
<p>A different scenario is observed in Whale Fall where, even though each of the subsamples were collected in a specific moment of decomposition from two different whales and at different depths, they all share similar general physico-chemical patterns, predominating the drastic fluctuations of nutrient availability
[<xref ref-type="bibr" rid="CR1">1</xref>
]. In agreement of what would be expected for microorganisms living in these kind of environments, most of the highly regulated functions that are enriched in whale falls samples are related to adaptation capabilities to starving periods (Figure
<xref rid="Fig3" ref-type="fig">3</xref>
). Particularly, we found TFBS rich promoters upstream of genes that are involved in cell cycle and growth, i.e. the control of basic macromolecular synthesis operon. This is in contrast to what happens in Waseca and Acid Mine, where the same functions present lower density of TFBSs (Figure
<xref rid="Fig3" ref-type="fig">3</xref>
).<fig id="Fig3"><label>Figure 3</label>
<caption><p><bold>Box-plots showing the TFBSs density per promoter related to highly regulated process in Acid Mine, Waseca Soil and Whale Falls samples.</bold>
The * indicates pvalue ≤ 0.05, and the ** indicates pvalue ≤ 0.01 (Fisher’s exact test) between the groups where the line is drawn. The numbers and functions listed at the right side represent the genes with highest regulatory potential. Their respective positions in the box-plot obey to their average level of regulation in a particular environment. In the X axis, are written the parental (above) and specific enriched group of functions (below) according to SEED classification.</p>
</caption>
<graphic xlink:href="12864_2014_6783_Fig3_HTML" id="d30e511"></graphic>
</fig>
</p>
<p>Moreover, bacterial communities living in cold water are also exposed to high concentrations of oxidant reagents
[<xref ref-type="bibr" rid="CR27">27</xref>
] causing an increase in the metabolic costs associated with the activation of antioxidant defenses. In fact, functions related with the response to oxidative stress appear specifically enriched in this environment compared to others. These functions comprise, for example, the hydrogen peroxide-inducible gene activator and a hem and copper containing membrane protein (NnrS), that needs to respond to external NO concentrations. Additionally, parts involved in the machinery that protects genomic DNA during prolonged non-growing phases
[<xref ref-type="bibr" rid="CR28">28</xref>
], like the non-specific DNA binding protein (Dps), also appear as highly regulated in this niche.</p>
<p>It has been also pointed before, that the uptake and metabolisms of cofactors and amino acids are particularly variable in marine environments, essential to adapt to typical oceanic oligotrophic conditions
[<xref ref-type="bibr" rid="CR2">2</xref>
]. In agreement with this, cofactor metabolism related functions are also enriched (adjusted pvalue (between groups) ≤ 0.05). In particular, we found enrichment for enzymes involved in the metabolism of molybdenum cofactors, pterin and folate (Figure
<xref rid="Fig3" ref-type="fig">3</xref>
). These findings were further confirmed by the overrepresentation of TyrR and ArgR binding sites in this niche, both known to be TFs involved in the control of amino acid transport for the synthesis of proteins (according to the RegPrecise database; see Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S8).</p>
</sec>
<sec id="Sec9"><title>Acid Mine</title>
<p>The acid mine is characterized by extreme physico-chemical conditions, showing low pH records and fluctuating temperature, conductivity and rainfall (see Figure
<xref rid="Fig4" ref-type="fig">4</xref>
A)
[<xref ref-type="bibr" rid="CR29">29</xref>
]. Among the functions with high regulatory potential that appear enriched in this niche are those known to play a role in the adaptation to changes in external osmolarity, typical of environments with variable distribution of rainfall across the year
[<xref ref-type="bibr" rid="CR30">30</xref>
, <xref ref-type="bibr" rid="CR31">31</xref>
] (Figure
<xref rid="Fig4" ref-type="fig">4</xref>
A). It is worth mentioning the high regulatory potential of some genes related to the TonB transport system (Figures
<xref rid="Fig3" ref-type="fig">3</xref>
,
<xref rid="Fig4" ref-type="fig">4</xref>
B), which are also involved in avoiding toxicity by keeping metal homeostasis inside the cell
[<xref ref-type="bibr" rid="CR32">32</xref>
], in particular of iron. The high regulatory potential of the TonB-dependent receptor and the iron chelator utilization protein (Figure
<xref rid="Fig3" ref-type="fig">3</xref>
) might provide homeostasis (i.e. plasticity) to acid mine bacteria living under variable ferric concentrations, which is further confirmed by the fact that a significant fraction of homeostasis-related promoters could be assigned to <italic>Leptospirillum</italic>
(genus known to be adapted to low pH
[<xref ref-type="bibr" rid="CR33">33</xref>
]) (Figure
<xref rid="Fig4" ref-type="fig">4</xref>
). In addition, we found overrepresentation of binding sites for LexA transcription factor in this niche (see Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S8), and, specifically in Ton and Tol transport systems related promoters (the sequence for LexA binding site is in Figure
<xref rid="Fig4" ref-type="fig">4</xref>
B, colored in red). LexA transcription factor is known to be involved in the response to DNA damage and external pH fluctuations
[<xref ref-type="bibr" rid="CR34">34</xref>
]. In fact, when we evaluated the fraction of binding sites shared between two members of the Ton and Tol system (iron chelator utilization protein and TonB dependent receptor), we found a high number of coincidences for other sites besides LexA (i.e. sites for the transcription factors ModR and ModE involved in metal metabolism) (Figure
<xref rid="Fig4" ref-type="fig">4</xref>
B).<fig id="Fig4"><label>Figure 4</label>
<caption><p><bold>Example of the relationship between environmental factors and organization of regulatory potential in Acid Mine. A)</bold>
Seasonal variation in pH, conductivity and iron concentration across the year. The left panel represents samples from “UBA” and right panel represents the samples collected from “5WAY” (data taken from [<xref ref-type="bibr" rid="CR29">29</xref>
]). <bold>B)</bold>
TFBSs found in the promoters which regulate Ton and Tol transport system genes. The left panel shows a graphic illustrating the regulatory potential related to Ton and Tol transport system in Acid Mine (“UBA”). The promoters are grouped by function and species (X axis). Each bar in the graphic represents a promoter and, the red color in the bars indicates that the number of sites per promoter is equal or more than fifteen. Below are represented the most abundant palindromic sites found in those promoters, some of them corresponding to known TF, such as, LexA, ModR, AgaR, etc. The right panel represents a network of TFBSs shared between promoters of TonB dependent receptor and iron chelator utilization protein; both proteins are represented by white circles. The TFBSs are represented by colored circles with size proportional to the number of sites found.</p>
</caption>
<graphic xlink:href="12864_2014_6783_Fig4_HTML" id="d30e603"></graphic>
</fig>
</p>
<p>Taken together, the fact that highly regulated functions are not the same between the different environments agrees with previous metatranscriptome studies
[<xref ref-type="bibr" rid="CR3">3</xref>
, <xref ref-type="bibr" rid="CR4">4</xref>
] and indicates that the organization of the regulatory potential between the functional space of each niche is different and influenced by the environmental physico-chemical conditions. This could reflect organism-environment interaction points where gene regulation should be able to provide enough plasticity to the functional network for the adaptation to variable external parameters.</p>
</sec>
</sec>
<sec id="Sec10" sec-type="conclusions"><title>Conclusions</title>
<p>We have here studied how variable physico-chemical conditions of the environment can shape the regulome of microbial communities living there to provide adaptation. We have combined existing and novel methodologies and applied it to three environments (Acid Mine, Whale Fall and Waseca Farm Soil) to identify and characterize, for the first time, their regulatory space, i.e. proximal promoters and their corresponding TFBSs. Taking the density of TFBSs as a measure of the level of regulatory potential, we first observed that, despite the differences of the living conditions of each of the environments studied here, their distribution of the regulatory potential, at quantitative level, appears to be nearly identical. However, when we went beyond simple counts we observed that the associated cellular functions in different groups related to the regulatory potential tend to be environment specific. This supports our hypothesis and expectation that point to a role of gene regulation in the adaptation of organisms to particular and variable external factors. Also in this direction, we have found specific functional enrichments among highly regulated functions in each of the metagenomas, suggesting potentials interaction points between gene regulation and dynamic environmental conditions. In particular, we have identified points of interaction between signatures of significant functional enrichment and specific characteristics of the marine and terrestrial environments. These results highlight the impact of gene regulation in the adaptation of microbes to their habitat. Beyond contributing to the general understanding of how wild bacterial communities interact with the environment, our methodology can also be used to identify potential external factors to which bacteria are particularly sensitive in order to design efficient communities for therapeutic, or ecological needs.</p>
</sec>
<sec id="Sec11" sec-type="methods"><title>Methods</title>
<sec id="Sec12"><title>Datasets</title>
<p>Metagenomic samples (i.e. Sanger sequencing reads) were downloaded from the Camera Database
[<xref ref-type="bibr" rid="CR35">35</xref>
]. In particular, (1) samples of whale falls were obtained from three independent libraries named Whale falls: CAM_SMPL_WHALEFALLBONE (Whale fall carcass bone, W. Antarctic Peninsula Shelf), CAM_SMPL_WHALEFALLMAT (Whale fall carcass microbial mat, Santa Cruz Basin), and CAM_SMPL_WHALEFALLRIB (Whale fall carcass rib bone, Santa Cruz Basin) in the database. These three sets differ in the depth of the sampling and come from two different whale samples. (2) The Acid Mine dataset is formed by 5-Way (CG) Acid Mine Drainage Biofilm Metagenome and UBA Acid Mine Drainage Biofilm Metagenome reads. The first corresponds to a low-complexity microbial biofilm growing hundreds of feet underground within a pyrite (FeS<sub>2</sub>
) ore body. The UBA biofilm was subaerial, collected from the base of a ~2 m high pile of pyrite sediment. (3) The third environment corresponds to a surface soil (0–10 cm) collected from a Waseca County farm in Minnesota.</p>
</sec>
<sec id="Sec13"><title>Promoter identification</title>
<p>The prediction and classification of regulatory regions from metagenomic data relies on the extraction of DNA regions upstream of coding genes detected through homology searches directly from the sequencing reads. For this protocol we selected conservative filters to ensure the reliability of the putative promoters found. Simplified in Figure
<xref rid="Fig1" ref-type="fig">1</xref>
, our protocol consisted in: (1) filtering out reads shorter than 800 base pairs. This filter keeps up to 90% of all reads and ensures both, the detection of the coding region and the extraction of the putative promoter from the reads; (2) detection of reads with coding potential through the comparison of all the sequences of each metagenome with all bacterial and archaeal annotated proteins (NCBI; <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/Ftp/">http://www.ncbi.nlm.nih.gov/Ftp/</ext-link>
), using BLASTx (default parameters
[<xref ref-type="bibr" rid="CR36">36</xref>
], and selecting those reads with a match to a known protein over, at least, 150 amino acids and with more than 50% of sequence identity; (3) filtering out those positive reads that did not contain at least 300 nucleotide of non-coding sequence upstream of the region matching in BLASTx. This filter enriches our sampling in regions with regulatory potential by avoiding internal genes of operons, which are expected to have short upstream regions with no regulatory potential. Finally, from the remaining accepted reads (13572, 3017, and 3215, for Acid Mine Drainage, Waseca Farm Soil and Whale Falls Samples, respectively) we extracted 300 nucleotides upstream of the coding region as putative promoter sequence. We expect that the 300 base pairs criteria will affect equally all bacteria and environments and will not favor bacteria with largest genomes, as this length has been also described for <italic>Pelagibacter ubique</italic>
, the free living bacteria with the smallest genome known
[<xref ref-type="bibr" rid="CR37">37</xref>
]. Moreover, fixing this length also avoids short intergenic regions within operons, as their regulatory role is not yet well understood.</p>
<p>To avoid other possible biases favoring common species in these environments and to make possible comparative and qualitative analyses between them, we also removed the redundancy within these collections of putative promoter sequences using a cutoff of 98% of sequence identity. We also removed those reads that correspond to eukaryotic DNA, mostly from plant species in the Waseca sample, identified using MEGAN
[<xref ref-type="bibr" rid="CR25">25</xref>
]. To discard the inclusion of (parts of) ncRNA genes into the collection of promoters, we applied a second filter to remove ncRNAs that target untranslated 5’ portions of mRNAs by using Rfam
[<xref ref-type="bibr" rid="CR38">38</xref>
] and also we did a second prediction of coding region in our set of putative promoters using the software Prodigal
[<xref ref-type="bibr" rid="CR39">39</xref>
] that allows the identification of genes even if the specie is unknown.</p>
</sec>
<sec id="Sec14"><title>Prediction of transcription factor binding sites</title>
<p>We next searched for sequence motifs with binding potential within the putative promoters identified before. For this, we used a <italic>de novo</italic>
prediction method that is based on the identification of palindromic repeats separated by a spacer DNA region. In particular, we used the most recent adaptations of the method
[<xref ref-type="bibr" rid="CR14">14</xref>
] originally described by Li and coworkers
[<xref ref-type="bibr" rid="CR9">9</xref>
].</p>
<p>In order to identify putative cis regulatory elements, we screened each promoter sequence for <italic>W</italic>
<sub>1</sub>
<italic>NW</italic>
<sub>2</sub>
, DNA motifs, where <italic>W</italic>
<sub><italic>1</italic>
</sub>
and <italic>W</italic>
<sub>2</sub>
are 3–5 nucleotide long palindromic sequences separated by <italic>N</italic>
(0–30) arbitrary bases. This method relies on the fact that prokaryotic TFBSs are usually palindromes between 12 and 30 base pairs, which may facilitate the dimerization and binding of TFs
[<xref ref-type="bibr" rid="CR12">12</xref>
].</p>
<p>To assign a probabilistic values to all motifs found, we first calculated the probability of observing <italic>n</italic>
(<italic>D</italic>
) copies of a dimer <italic>D</italic>
by chance, by pooling all the promoters and calculating its expected frequency from the formula,
<disp-formula id="Equ1"><label>1</label>
<graphic xlink:href="12864_2014_6783_Equ1_HTML.gif" position="anchor"></graphic>
</disp-formula>
</p>
<p>where <italic>n</italic>
(<italic>W</italic>
1) and <italic>n</italic>
(<italic>W</italic>
2) are the total number of occurrences of <italic>W</italic>
<sub>1</sub>
and <italic>W</italic>
<sub>2</sub>
in the whole data set (all three environments together) and <italic>L</italic>
eff(D) = ∑<italic>r</italic>
(<italic>L</italic>
(<italic>r</italic>
) − <italic>L</italic>
(<italic>D</italic>
) + 1) is the number of independent positions in the data where a motif <italic>D</italic>
of length <italic>L</italic>
(<italic>D</italic>
) can be found. The summation is over all the occurrences among 11,614 promoters identified, each with a length <italic>L</italic>
(<italic>r</italic>
) (i.e. the estimated distance between coding regions). Finally, a <italic>P</italic>
-value is assigned to each of the motifs assuming that the background follows a Poisson distribution:
<disp-formula id="Equ2"><label>2</label>
<graphic xlink:href="12864_2014_6783_Equ2_HTML.gif" position="anchor"></graphic>
</disp-formula>
</p>
<p>and is considered significant if <italic>P</italic>
< 1/<italic>N</italic>
<sub><italic>motif</italic>
</sub>
, where <italic>N</italic>
<sub><italic>motif</italic>
</sub>
is the total number of positive motifs found. As <italic>W</italic>
<sub>1</sub>
is the reverse complement of W<sub>2</sub>
(palindrome), the cutoff on <italic>P</italic>
is corrected by the total number of palindromic dimers found
[<xref ref-type="bibr" rid="CR9">9</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
].</p>
<p>In order to identify environment specific enrichment of our know TFBSs (i.e. those present in the RegPrecise database), we run a Kruskal Wallis test to compare the density of each particular known TFBS among all three environments. The density of known TFBS per metagenomes is calculated as follows:
<disp-formula id="Equ3"><label>3</label>
<graphic xlink:href="12864_2014_6783_Equ3_HTML.gif" position="anchor"></graphic>
</disp-formula>
</p>
<p>where D (x) is the density of TFBSs per metagenome, N represents the number of promoters found in the x metagenome and Tbp is the number of base pairs per promoter (300 base pairs). The complete list of overrepresented TFBSs found in our selected promoter set are shown in Additional files
<xref rid="MOESM6" ref-type="media">6</xref>
,
<xref rid="MOESM8" ref-type="media">7</xref>
,
<xref rid="MOESM8" ref-type="media">8</xref>
,
<xref rid="MOESM9" ref-type="media">9</xref>
,
<xref rid="MOESM10" ref-type="media">10</xref>
and
<xref rid="MOESM11" ref-type="media">11</xref>
, for Acid Mine, Waseca Farm Soil and Whale Falls samples, respectively.</p>
</sec>
<sec id="Sec15"><title>Method validation</title>
<p>For the randomization test on TFBS prediction, we run the corresponding searching methodology on predicted promoter regions after shuffling their sequence using a 20 nt window to ensure the minimum variance of local nucleotide composition.</p>
<p>For the comparison with the MotifClick method
[<xref ref-type="bibr" rid="CR20">20</xref>
] we first downloaded intergenic regions from the <italic>Escherichia coli K12-W3110</italic>
genome from IMG database (<ext-link ext-link-type="uri" xlink:href="https://img.jgi.doe.gov">https://img.jgi.doe.gov</ext-link>
). We ran MotifClick (motif length = 14 nt) over these regions, specifically 300 nucleotide upstream annotated TSS and recorded the number of positive predictions per promoter. These values were then compared with the results provided by our method applied on the same set of <italic>E. Coli</italic>
regulatory regions (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
).</p>
</sec>
<sec id="Sec16"><title>Statistical procedure for the functional analysis</title>
<p>Functional assignment for all the data was performed by MEGAN software
[<xref ref-type="bibr" rid="CR25">25</xref>
] using the output of BLAST searches of our reads against databases of known bacterial proteins. Through this comparison we could identify up to 1646 (Whale Falls Samples), 4646 (Acid Mine Drainage) and 1514 (Minnesota Farm Soil) gene upstream segments with functional assignment. In order to roughly study up to which level low, medium and high regulated functions are shared among environment we have run a Spearman test for independence using R, for the rectangular plot and correspondence analysis we use the plot function included in R graphics (<ext-link ext-link-type="uri" xlink:href="http://www.r-project.org/">http://www.r-project.org/</ext-link>
) (see Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Figure S4).</p>
<p>In addition, functional enrichment analysis was done by first ranking all promoters as to their number of predicted TFBSs. Then, for each of the groups of interest, we ran a Fisher’s exact test for count data to see whether particular functions within each group (top 1%, 5%, 10%, 20%) were specifically enriched versus the total distribution of functions. For this, we have used “all intermediate” functional levels according to MEGAN classification. Heat maps for all function within environment were obtained using package ggplot2 for R (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
). Then, we retained significant cases based on two criteria 1) functions whose p < < 0.05 within environment and 2) functions with orthologous in the other three environments. Those selected groups were compared again, this time among environments, for this analysis we ran a Fisher’s exact test to see whether functional enrichment within environment were maintained among them.</p>
</sec>
</sec>
<sec sec-type="supplementary-material"><title>Electronic supplementary material</title>
<sec id="Sec17"><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="12864_2014_6783_MOESM1_ESM.pdf"><caption><p>Additional file 1: Figure S1: Shows the overview of the general results of this study. <bold>Figure S2.</bold>
shows the comparative analysis of the taxa obtained with MEGAN on our promoter regions compared with that obtained previously using 16S rRNA information from the same samples in Waseca soil (a), Whale falls (b), and Acid mine (c). <bold>Figure S3.</bold>
represents the correlation analysis between the TFBSs predictions per promoter using the method explained in this paper versus MotifClick predictions. <bold>Figure S4.</bold>
illustrates a global view of the relationship between regulatory potential and the level of co-occurring functions within each of the environments. <bold>Figure S5.</bold>
Results of the functional enrichment analysis for Acid Mine using the predefined bins. <bold>Figure S6.</bold>
Results of the functional enrichment analysis for Waseca Farm using predefined bins. <bold>Figure S7.</bold>
Results of the functional enrichment analysis for Whale Falls using predefined bins. <bold>Figure S8.</bold>
shows the relative abundances of our TFBS prediction that matched known TFBS. (PDF 4 MB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM2"><media xlink:href="12864_2014_6783_MOESM2_ESM.zip"><caption><p>Additional file 2: <bold>List of promoters selected after applying the methodology described in Figure</bold>
<xref rid="Fig1" ref-type="fig">1</xref>
<bold>on Waseca Farm Soil data.</bold>
(ZIP 261 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM3"><media xlink:href="12864_2014_6783_MOESM3_ESM.zip"><caption><p>Additional file 3: <bold>List of the promoters selected after applying the methodology described in Figure</bold>
<xref rid="Fig1" ref-type="fig">1</xref>
<bold>on Acid Mine data.</bold>
(ZIP 817 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM4"><media xlink:href="12864_2014_6783_MOESM4_ESM.zip"><caption><p>Additional file 4: <bold>List of the promoters selected after applying the methodology described in Figure</bold>
<xref rid="Fig1" ref-type="fig">1</xref>
<bold>on Whale Fall Samples data.</bold>
(ZIP 264 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM5"><media xlink:href="12864_2014_6783_MOESM5_ESM.pdf"><caption><p>Additional file 5:
<bold>Table listing the number of TFBSs per genomes found after applying our method versus the number of sites described in Regprecise database.</bold>
(PDF 88 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM6"><media xlink:href="12864_2014_6783_MOESM6_ESM.csv"><caption><p>Additional file 6:
<bold>Table in CSV format listing the number of TFBSs identified for each promoter and the function assigned to the corresponding downstream coding region in Acid Mine.</bold>
(CSV 542 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM7"><media xlink:href="12864_2014_6783_MOESM7_ESM.csv"><caption><p>Additional file 7:
<bold>Table in CSV format listing the number of TFBSs identified for each promoter and the function assigned to the corresponding downstream coding region in Waseca Soil.</bold>
(CSV 172 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM8"><media xlink:href="12864_2014_6783_MOESM8_ESM.csv"><caption><p>Additional file 8:
<bold>Table in CSV format listing the number of TFBSs identified for each promoter and the function assigned to the corresponding downstream coding region in Whale Falls.</bold>
(CSV 151 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM9"><media xlink:href="12864_2014_6783_MOESM9_ESM.csv"><caption><p>Additional file 9: <bold>A list (CSV MS-DOS format) of overrepresented TFBSs per promoter found in Acid Mine, Waseca Soils and Whale falls, respectively.</bold>
The abbreviated nomenclature used for the binding sites is the following: N, W, Sequence, where N is the number of variable nucleotides. W is the number of nucleotides defining the inverted repeat. Sequence is the actual sequence of the site. Example: 10 3 ATC, corresponds to: ATCNNNNNNNNNNGAT. (CSV 3 MB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM10"><media xlink:href="12864_2014_6783_MOESM10_ESM.csv"><caption><p>Additional file 10: <bold>A list (CSV MS-DOS format) of overrepresented TFBSs per promoter found in Acid Mine, Waseca Soils and Whale falls, respectively.</bold>
The abbreviated nomenclature used for the binding sites is the following: N, W, Sequence, where N is the number of variable nucleotides. W is the number of nucleotides defining the inverted repeat. Sequence is the actual sequence of the site. Example: 10 3 ATC, corresponds to: ATCNNNNNNNNNNGAT. (CSV 873 KB)</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM11"><media xlink:href="12864_2014_6783_MOESM11_ESM.csv"><caption><p>Additional file 11: <bold>A list (CSV MS-DOS format) of overrepresented TFBSs per promoter found in Acid Mine, Waseca Soils and Whale falls, respectively.</bold>
The abbreviated nomenclature used for the binding sites is the following: N, W, Sequence, where N is the number of variable nucleotides. W is the number of nucleotides defining the inverted repeat. Sequence is the actual sequence of the site. Example: 10 3 ATC, corresponds to: ATCNNNNNNNNNNGAT. (CSV 1 MB)</p>
</caption>
</media>
</supplementary-material>
</sec>
</sec>
</body>
<back><fn-group><fn><p><bold>Competing interests</bold>
</p>
<p>The authors declare that they have no competing interests.</p>
</fn>
<fn><p><bold>Authors’ contributions</bold>
</p>
<p>DT conceived this study. LF and DT designed the methodology. LF analyzed the data and performs the literature searches. JMM and LF participated in the statistical analysis. LF and MP performed the regulatory network analysis of TFBSs. JMM, DT, MP and LF have been involved in drafting the manuscript. All authors read and approved the final version of the manuscript.</p>
</fn>
</fn-group>
<ack><p>This work was supported by a grant from Ministerio de Economía y Competitividad through the project BES-2008-005973. Josep M. Mercader was supported by Sara Borrell Fellowship from the Instituto Carlos III. This work has been supported by the grant SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government. Mercè Planas is funded by the Obra Social Fundación la Caixa under the Severo Ochoa 2013 program.</p>
</ack>
<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>Tringe</surname>
<given-names>SG</given-names>
</name>
<name><surname>von Mering</surname>
<given-names>C</given-names>
</name>
<name><surname>Kobayashi</surname>
<given-names>A</given-names>
</name>
<name><surname>Salamov</surname>
<given-names>AA</given-names>
</name>
<name><surname>Chen</surname>
<given-names>K</given-names>
</name>
<name><surname>Chang</surname>
<given-names>HW</given-names>
</name>
<name><surname>Podar</surname>
<given-names>M</given-names>
</name>
<name><surname>Short</surname>
<given-names>JM</given-names>
</name>
<name><surname>Mathur</surname>
<given-names>EJ</given-names>
</name>
<name><surname>Detter</surname>
<given-names>JC</given-names>
</name>
<name><surname>Bork</surname>
<given-names>P</given-names>
</name>
<name><surname>Hugenholtz</surname>
<given-names>P</given-names>
</name>
<name><surname>Rubin</surname>
<given-names>EM</given-names>
</name>
</person-group>
<article-title>Comparative metagenomics of microbial communities</article-title>
<source>Science</source>
<year>2005</year>
<volume>308</volume>
<fpage>554</fpage>
<lpage>557</lpage>
<pub-id pub-id-type="doi">10.1126/science.1107851</pub-id>
<pub-id pub-id-type="pmid">15845853</pub-id>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gianoulis</surname>
<given-names>TA</given-names>
</name>
<name><surname>Raes</surname>
<given-names>J</given-names>
</name>
<name><surname>Patel</surname>
<given-names>PV</given-names>
</name>
<name><surname>Bjornson</surname>
<given-names>R</given-names>
</name>
<name><surname>Korbel</surname>
<given-names>JO</given-names>
</name>
<name><surname>Letunic</surname>
<given-names>I</given-names>
</name>
<name><surname>Yamada</surname>
<given-names>T</given-names>
</name>
<name><surname>Paccanaro</surname>
<given-names>A</given-names>
</name>
<name><surname>Jensen</surname>
<given-names>LJ</given-names>
</name>
<name><surname>Snyder</surname>
<given-names>M</given-names>
</name>
<name><surname>Bork</surname>
<given-names>P</given-names>
</name>
<name><surname>Gerstein</surname>
<given-names>MB</given-names>
</name>
</person-group>
<article-title>Quantifying environmental adaptation of metabolic pathways in metagenomics</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2009</year>
<volume>106</volume>
<fpage>1374</fpage>
<lpage>1379</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0808022106</pub-id>
<pub-id pub-id-type="pmid">19164758</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poretsky</surname>
<given-names>RS</given-names>
</name>
<name><surname>Bano</surname>
<given-names>N</given-names>
</name>
<name><surname>Buchan</surname>
<given-names>A</given-names>
</name>
<name><surname>LeCleir</surname>
<given-names>G</given-names>
</name>
<name><surname>Kleikemper</surname>
<given-names>J</given-names>
</name>
<name><surname>Pickering</surname>
<given-names>M</given-names>
</name>
<name><surname>Pate</surname>
<given-names>WM</given-names>
</name>
<name><surname>Moran</surname>
<given-names>MA</given-names>
</name>
<name><surname>Hollibaugh</surname>
<given-names>JT</given-names>
</name>
</person-group>
<article-title>Analysis of microbial gene transcripts in environmental samples</article-title>
<source>Appl Environ Microbiol</source>
<year>2005</year>
<volume>71</volume>
<fpage>4121</fpage>
<lpage>4126</lpage>
<pub-id pub-id-type="doi">10.1128/AEM.71.7.4121-4126.2005</pub-id>
<pub-id pub-id-type="pmid">16000831</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poretsky</surname>
<given-names>RS</given-names>
</name>
<name><surname>Hewson</surname>
<given-names>I</given-names>
</name>
<name><surname>Sun</surname>
<given-names>S</given-names>
</name>
<name><surname>Allen</surname>
<given-names>AE</given-names>
</name>
<name><surname>Zehr</surname>
<given-names>JP</given-names>
</name>
<name><surname>Moran</surname>
<given-names>MA</given-names>
</name>
</person-group>
<article-title>Comparative day/night metatranscriptomic analysis of microbial communities in the North Pacific subtropical gyre</article-title>
<source>Environ Microbiol</source>
<year>2009</year>
<volume>11</volume>
<fpage>1358</fpage>
<lpage>1375</lpage>
<pub-id pub-id-type="doi">10.1111/j.1462-2920.2008.01863.x</pub-id>
<pub-id pub-id-type="pmid">19207571</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gifford</surname>
<given-names>SM</given-names>
</name>
<name><surname>Sharma</surname>
<given-names>S</given-names>
</name>
<name><surname>Rinta-Kanto</surname>
<given-names>JM</given-names>
</name>
<name><surname>Moran</surname>
<given-names>MA</given-names>
</name>
</person-group>
<article-title>Quantitative analysis of a deeply sequenced marine microbial metatranscriptome</article-title>
<source>ISME J</source>
<year>2011</year>
<volume>5</volume>
<fpage>461</fpage>
<lpage>472</lpage>
<pub-id pub-id-type="doi">10.1038/ismej.2010.141</pub-id>
<pub-id pub-id-type="pmid">20844569</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Farre</surname>
<given-names>D</given-names>
</name>
<name><surname>Bellora</surname>
<given-names>N</given-names>
</name>
<name><surname>Mularoni</surname>
<given-names>L</given-names>
</name>
<name><surname>Messeguer</surname>
<given-names>X</given-names>
</name>
<name><surname>Alba</surname>
<given-names>MM</given-names>
</name>
</person-group>
<article-title>Housekeeping genes tend to show reduced upstream sequence conservation</article-title>
<source>Genome Biol</source>
<year>2007</year>
<volume>8</volume>
<fpage>R140</fpage>
<pub-id pub-id-type="doi">10.1186/gb-2007-8-7-r140</pub-id>
<pub-id pub-id-type="pmid">17626644</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname>
<given-names>Z</given-names>
</name>
<name><surname>Wu</surname>
<given-names>WS</given-names>
</name>
<name><surname>Liang</surname>
<given-names>H</given-names>
</name>
<name><surname>Woo</surname>
<given-names>Y</given-names>
</name>
<name><surname>Li</surname>
<given-names>WH</given-names>
</name>
</person-group>
<article-title>The spatial distribution of cis regulatory elements in yeast promoters and its implications for transcriptional regulation</article-title>
<source>BMC Genomics</source>
<year>2010</year>
<volume>11</volume>
<fpage>581</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-11-581</pub-id>
<pub-id pub-id-type="pmid">20958978</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Merino</surname>
<given-names>E</given-names>
</name>
<name><surname>Jensen</surname>
<given-names>RA</given-names>
</name>
<name><surname>Yanofsky</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>Evolution of bacterial trp operons and their regulation</article-title>
<source>Curr Opin Microbiol</source>
<year>2008</year>
<volume>11</volume>
<fpage>78</fpage>
<lpage>86</lpage>
<pub-id pub-id-type="doi">10.1016/j.mib.2008.02.005</pub-id>
<pub-id pub-id-type="pmid">18374625</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
<name><surname>Rhodius</surname>
<given-names>V</given-names>
</name>
<name><surname>Gross</surname>
<given-names>C</given-names>
</name>
<name><surname>Siggia</surname>
<given-names>ED</given-names>
</name>
</person-group>
<article-title>Identification of the binding sites of regulatory proteins in bacterial genomes</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2002</year>
<volume>99</volume>
<fpage>11772</fpage>
<lpage>11777</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.112341999</pub-id>
<pub-id pub-id-type="pmid">12181488</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McCue</surname>
<given-names>L</given-names>
</name>
<name><surname>Thompson</surname>
<given-names>W</given-names>
</name>
<name><surname>Carmack</surname>
<given-names>C</given-names>
</name>
<name><surname>Ryan</surname>
<given-names>MP</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JS</given-names>
</name>
<name><surname>Derbyshire</surname>
<given-names>V</given-names>
</name>
<name><surname>Lawrence</surname>
<given-names>CE</given-names>
</name>
</person-group>
<article-title>Phylogenetic footprinting of transcription factor binding sites in proteobacterial genomes</article-title>
<source>Nucleic Acids Res</source>
<year>2001</year>
<volume>29</volume>
<fpage>774</fpage>
<lpage>782</lpage>
<pub-id pub-id-type="doi">10.1093/nar/29.3.774</pub-id>
<pub-id pub-id-type="pmid">11160901</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huffman</surname>
<given-names>JL</given-names>
</name>
<name><surname>Brennan</surname>
<given-names>RG</given-names>
</name>
</person-group>
<article-title>Prokaryotic transcription regulators: more than just the helix-turn-helix motif</article-title>
<source>Curr Opin Struct Biol</source>
<year>2002</year>
<volume>12</volume>
<fpage>98</fpage>
<lpage>106</lpage>
<pub-id pub-id-type="doi">10.1016/S0959-440X(02)00295-6</pub-id>
<pub-id pub-id-type="pmid">11839496</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rodionov</surname>
<given-names>DA</given-names>
</name>
</person-group>
<article-title>Comparative genomic reconstruction of transcriptional regulatory networks in bacteria</article-title>
<source>Chem Rev</source>
<year>2007</year>
<volume>107</volume>
<fpage>3467</fpage>
<lpage>3497</lpage>
<pub-id pub-id-type="doi">10.1021/cr068309+</pub-id>
<pub-id pub-id-type="pmid">17636889</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gelfand</surname>
<given-names>MS</given-names>
</name>
<name><surname>Koonin</surname>
<given-names>EV</given-names>
</name>
<name><surname>Mironov</surname>
<given-names>AA</given-names>
</name>
</person-group>
<article-title>Prediction of transcription regulatory sites in Archaea by a comparative genomic approach</article-title>
<source>Nucleic Acids Res</source>
<year>2000</year>
<volume>28</volume>
<fpage>695</fpage>
<lpage>705</lpage>
<pub-id pub-id-type="doi">10.1093/nar/28.3.695</pub-id>
<pub-id pub-id-type="pmid">10637320</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Laing</surname>
<given-names>E</given-names>
</name>
<name><surname>Sidhu</surname>
<given-names>K</given-names>
</name>
<name><surname>Hubbard</surname>
<given-names>SJ</given-names>
</name>
</person-group>
<article-title>Predicted transcription factor binding sites as predictors of operons in Escherichia coli and Streptomyces coelicolor</article-title>
<source>BMC Genomics</source>
<year>2008</year>
<volume>9</volume>
<fpage>79</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-9-79</pub-id>
<pub-id pub-id-type="pmid">18269733</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iqbal</surname>
<given-names>M</given-names>
</name>
<name><surname>Mast</surname>
<given-names>Y</given-names>
</name>
<name><surname>Amin</surname>
<given-names>R</given-names>
</name>
<name><surname>Hodgson</surname>
<given-names>DA</given-names>
</name>
<name><surname>Consortium</surname>
<given-names>S</given-names>
</name>
<name><surname>Wohlleben</surname>
<given-names>W</given-names>
</name>
<name><surname>Burroughs</surname>
<given-names>NJ</given-names>
</name>
</person-group>
<article-title>Extracting regulator activity profiles by integration of de novo motifs and expression data: characterizing key regulators of nutrient depletion responses in Streptomyces coelicolor</article-title>
<source>Nucleic Acids Res</source>
<year>2012</year>
<volume>40</volume>
<fpage>5227</fpage>
<lpage>5239</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gks205</pub-id>
<pub-id pub-id-type="pmid">22406834</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>GADEM: a genetic algorithm guided formation of spaced dyads coupled with an EM algorithm for motif discovery</article-title>
<source>J Comput Biol</source>
<year>2009</year>
<volume>16</volume>
<fpage>317</fpage>
<lpage>329</lpage>
<pub-id pub-id-type="doi">10.1089/cmb.2008.16TT</pub-id>
<pub-id pub-id-type="pmid">19193149</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname>
<given-names>J</given-names>
</name>
<name><surname>Xu</surname>
<given-names>X</given-names>
</name>
<name><surname>Stormo</surname>
<given-names>GD</given-names>
</name>
</person-group>
<article-title>The cis-regulatory map of Shewanella genomes</article-title>
<source>Nucleic Acids Res</source>
<year>2008</year>
<volume>36</volume>
<fpage>5376</fpage>
<lpage>5390</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkn515</pub-id>
<pub-id pub-id-type="pmid">18701645</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname>
<given-names>J</given-names>
</name>
<name><surname>Tuncay</surname>
<given-names>K</given-names>
</name>
<name><surname>Haidar</surname>
<given-names>AA</given-names>
</name>
<name><surname>Ensman</surname>
<given-names>L</given-names>
</name>
<name><surname>Stanley</surname>
<given-names>F</given-names>
</name>
<name><surname>Trelinski</surname>
<given-names>M</given-names>
</name>
<name><surname>Ortoleva</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Transcriptional regulatory network discovery via multiple method integration: application to e. coli K12</article-title>
<source>Algorithms Mol Biol</source>
<year>2007</year>
<volume>2</volume>
<fpage>2</fpage>
<pub-id pub-id-type="doi">10.1186/1748-7188-2-2</pub-id>
<pub-id pub-id-type="pmid">17397539</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Leuze</surname>
<given-names>MR</given-names>
</name>
<name><surname>Karpinets</surname>
<given-names>TV</given-names>
</name>
<name><surname>Syed</surname>
<given-names>MH</given-names>
</name>
<name><surname>Beliaev</surname>
<given-names>AS</given-names>
</name>
<name><surname>Uberbacher</surname>
<given-names>EC</given-names>
</name>
</person-group>
<article-title>Binding Motifs in Bacterial Gene Promoters Modulate Transcriptional Effects of Global Regulators CRP and ArcA</article-title>
<source>Gene Regul Syst Bio</source>
<year>2012</year>
<volume>6</volume>
<fpage>93</fpage>
<lpage>107</lpage>
<pub-id pub-id-type="doi">10.4137/GRSB.S9357</pub-id>
<pub-id pub-id-type="pmid">22701314</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>S</given-names>
</name>
<name><surname>Li</surname>
<given-names>S</given-names>
</name>
<name><surname>Niu</surname>
<given-names>M</given-names>
</name>
<name><surname>Pham</surname>
<given-names>PT</given-names>
</name>
<name><surname>Su</surname>
<given-names>Z</given-names>
</name>
</person-group>
<article-title>MotifClick: prediction of cis-regulatory binding sites via merging cliques</article-title>
<source>BMC Bioinformatics</source>
<year>2011</year>
<volume>12</volume>
<fpage>238</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2105-12-238</pub-id>
<pub-id pub-id-type="pmid">21679436</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Novichkov</surname>
<given-names>PS</given-names>
</name>
<name><surname>Laikova</surname>
<given-names>ON</given-names>
</name>
<name><surname>Novichkova</surname>
<given-names>ES</given-names>
</name>
<name><surname>Gelfand</surname>
<given-names>MS</given-names>
</name>
<name><surname>Arkin</surname>
<given-names>AP</given-names>
</name>
<name><surname>Dubchak</surname>
<given-names>I</given-names>
</name>
<name><surname>Rodionov</surname>
<given-names>DA</given-names>
</name>
</person-group>
<article-title>RegPrecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes</article-title>
<source>Nucleic Acids Res</source>
<year>2010</year>
<volume>38</volume>
<fpage>D111</fpage>
<lpage>118</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkp894</pub-id>
<pub-id pub-id-type="pmid">19884135</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grissa</surname>
<given-names>I</given-names>
</name>
<name><surname>Vergnaud</surname>
<given-names>G</given-names>
</name>
<name><surname>Pourcel</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats</article-title>
<source>BMC Bioinformatics</source>
<year>2007</year>
<volume>8</volume>
<fpage>172</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2105-8-172</pub-id>
<pub-id pub-id-type="pmid">17521438</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grissa</surname>
<given-names>I</given-names>
</name>
<name><surname>Vergnaud</surname>
<given-names>G</given-names>
</name>
<name><surname>Pourcel</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats</article-title>
<source>Nucleic Acids Res</source>
<year>2007</year>
<volume>35</volume>
<fpage>W52</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkm360</pub-id>
<pub-id pub-id-type="pmid">17537822</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Overbeek</surname>
<given-names>R</given-names>
</name>
<name><surname>Begley</surname>
<given-names>T</given-names>
</name>
<name><surname>Butler</surname>
<given-names>RM</given-names>
</name>
<name><surname>Choudhuri</surname>
<given-names>JV</given-names>
</name>
<name><surname>Chuang</surname>
<given-names>HY</given-names>
</name>
<name><surname>Cohoon</surname>
<given-names>M</given-names>
</name>
<name><surname>de Crecy-Lagard</surname>
<given-names>V</given-names>
</name>
<name><surname>Diaz</surname>
<given-names>N</given-names>
</name>
<name><surname>Disz</surname>
<given-names>T</given-names>
</name>
<name><surname>Edwards</surname>
<given-names>R</given-names>
</name>
<name><surname>Fonstein</surname>
<given-names>M</given-names>
</name>
<name><surname>Frank</surname>
<given-names>ED</given-names>
</name>
<name><surname>Gerdes</surname>
<given-names>S</given-names>
</name>
<name><surname>Glass</surname>
<given-names>EM</given-names>
</name>
<name><surname>Goesmann</surname>
<given-names>A</given-names>
</name>
<name><surname>Hanson</surname>
<given-names>A</given-names>
</name>
<name><surname>Iwata-Reuyl</surname>
<given-names>D</given-names>
</name>
<name><surname>Jensen</surname>
<given-names>R</given-names>
</name>
<name><surname>Jamshidi</surname>
<given-names>N</given-names>
</name>
<name><surname>Krause</surname>
<given-names>L</given-names>
</name>
<name><surname>Kubal</surname>
<given-names>M</given-names>
</name>
<name><surname>Larsen</surname>
<given-names>N</given-names>
</name>
<name><surname>Linke</surname>
<given-names>B</given-names>
</name>
<name><surname>McHardy</surname>
<given-names>AC</given-names>
</name>
<name><surname>Meyer</surname>
<given-names>F</given-names>
</name>
<name><surname>Neuweger</surname>
<given-names>H</given-names>
</name>
<name><surname>Olsen</surname>
<given-names>G</given-names>
</name>
<name><surname>Olson</surname>
<given-names>R</given-names>
</name>
<name><surname>Osterman</surname>
<given-names>A</given-names>
</name>
<name><surname>Portnoy</surname>
<given-names>V</given-names>
</name>
</person-group>
<article-title>The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes</article-title>
<source>Nucleic Acids Res</source>
<year>2005</year>
<volume>33</volume>
<fpage>5691</fpage>
<lpage>5702</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gki866</pub-id>
<pub-id pub-id-type="pmid">16214803</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huson</surname>
<given-names>DH</given-names>
</name>
<name><surname>Auch</surname>
<given-names>AF</given-names>
</name>
<name><surname>Qi</surname>
<given-names>J</given-names>
</name>
<name><surname>Schuster</surname>
<given-names>SC</given-names>
</name>
</person-group>
<article-title>MEGAN analysis of metagenomic data</article-title>
<source>Genome Res</source>
<year>2007</year>
<volume>17</volume>
<fpage>377</fpage>
<lpage>386</lpage>
<pub-id pub-id-type="doi">10.1101/gr.5969107</pub-id>
<pub-id pub-id-type="pmid">17255551</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chin</surname>
<given-names>CS</given-names>
</name>
<name><surname>Chuang</surname>
<given-names>JH</given-names>
</name>
<name><surname>Li</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Genome-wide regulatory complexity in yeast promoters: separation of functionally conserved and neutral sequence</article-title>
<source>Genome Res</source>
<year>2005</year>
<volume>15</volume>
<fpage>205</fpage>
<lpage>213</lpage>
<pub-id pub-id-type="doi">10.1101/gr.3243305</pub-id>
<pub-id pub-id-type="pmid">15653830</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abele</surname>
<given-names>D</given-names>
</name>
<name><surname>Puntarulo</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Formation of reactive species and induction of antioxidant defence systems in polar and temperate marine invertebrates and fish</article-title>
<source>Comp Biochem Physiol A Mol Integr Physiol</source>
<year>2004</year>
<volume>138</volume>
<fpage>405</fpage>
<lpage>415</lpage>
<pub-id pub-id-type="doi">10.1016/j.cbpb.2004.05.013</pub-id>
<pub-id pub-id-type="pmid">15369829</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Storz</surname>
<given-names>G</given-names>
</name>
<name><surname>Imlay</surname>
<given-names>JA</given-names>
</name>
</person-group>
<article-title>Oxidative stress</article-title>
<source>Curr Opin Microbiol</source>
<year>1999</year>
<volume>2</volume>
<fpage>188</fpage>
<lpage>194</lpage>
<pub-id pub-id-type="doi">10.1016/S1369-5274(99)80033-2</pub-id>
<pub-id pub-id-type="pmid">10322176</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edwards</surname>
<given-names>KJ</given-names>
</name>
<name><surname>Gihring</surname>
<given-names>TM</given-names>
</name>
<name><surname>Banfield</surname>
<given-names>JF</given-names>
</name>
</person-group>
<article-title>Seasonal variations in microbial populations and environmental conditions in an extreme acid mine drainage environment</article-title>
<source>Appl Environ Microbiol</source>
<year>1999</year>
<volume>65</volume>
<fpage>3627</fpage>
<lpage>3632</lpage>
<pub-id pub-id-type="pmid">10427059</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Albers</surname>
<given-names>S-V</given-names>
</name>
<name><surname>Koning</surname>
<given-names>SM</given-names>
</name>
<name><surname>Konings</surname>
<given-names>WN</given-names>
</name>
<name><surname>Driessen</surname>
<given-names>AJM</given-names>
</name>
</person-group>
<article-title>Insights into ABC Transport in Archaea</article-title>
<source>J Bioenerg Biomembr</source>
<year>2004</year>
<volume>36</volume>
<fpage>5</fpage>
<lpage>15</lpage>
<pub-id pub-id-type="doi">10.1023/B:JOBB.0000019593.84933.e6</pub-id>
<pub-id pub-id-type="pmid">15168605</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kempf</surname>
<given-names>B</given-names>
</name>
<name><surname>Bremer</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>Uptake and synthesis of compatible solutes as microbial stress responses to high-osmolality environments</article-title>
<source>Arch Microbiol</source>
<year>1998</year>
<volume>170</volume>
<fpage>319</fpage>
<lpage>330</lpage>
<pub-id pub-id-type="doi">10.1007/s002030050649</pub-id>
<pub-id pub-id-type="pmid">9818351</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Osorio</surname>
<given-names>H</given-names>
</name>
<name><surname>Martinez</surname>
<given-names>V</given-names>
</name>
<name><surname>Nieto</surname>
<given-names>PA</given-names>
</name>
<name><surname>Holmes</surname>
<given-names>DS</given-names>
</name>
<name><surname>Quatrini</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Microbial iron management mechanisms in extremely acidic environments: comparative genomics evidence for diversity and versatility</article-title>
<source>BMC Microbiol</source>
<year>2008</year>
<volume>8</volume>
<fpage>203</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2180-8-203</pub-id>
<pub-id pub-id-type="pmid">19025650</pub-id>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schrenk</surname>
<given-names>MO</given-names>
</name>
<name><surname>Edwards</surname>
<given-names>KJ</given-names>
</name>
<name><surname>Goodman</surname>
<given-names>RM</given-names>
</name>
<name><surname>Hamers</surname>
<given-names>RJ</given-names>
</name>
<name><surname>Banfield</surname>
<given-names>JF</given-names>
</name>
</person-group>
<article-title>Distribution of Thiobacillus ferrooxidans and Leptospirillum ferrooxidans: Implications for generation of acid mine drainage</article-title>
<source>Science</source>
<year>1998</year>
<volume>279</volume>
<fpage>1519</fpage>
<lpage>1522</lpage>
<pub-id pub-id-type="doi">10.1126/science.279.5356.1519</pub-id>
<pub-id pub-id-type="pmid">9488647</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guazzaroni</surname>
<given-names>ME</given-names>
</name>
<name><surname>Morgante</surname>
<given-names>V</given-names>
</name>
<name><surname>Mirete</surname>
<given-names>S</given-names>
</name>
<name><surname>Gonzalez-Pastor</surname>
<given-names>JE</given-names>
</name>
</person-group>
<article-title>Novel acid resistance genes from the metagenome of the Tinto River, an extremely acidic environment</article-title>
<source>Environ Microbiol</source>
<year>2013</year>
<volume>15</volume>
<fpage>1088</fpage>
<lpage>1102</lpage>
<pub-id pub-id-type="doi">10.1111/1462-2920.12021</pub-id>
<pub-id pub-id-type="pmid">23145860</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname>
<given-names>S</given-names>
</name>
<name><surname>Chen</surname>
<given-names>J</given-names>
</name>
<name><surname>Li</surname>
<given-names>W</given-names>
</name>
<name><surname>Altintas</surname>
<given-names>I</given-names>
</name>
<name><surname>Lin</surname>
<given-names>A</given-names>
</name>
<name><surname>Peltier</surname>
<given-names>S</given-names>
</name>
<name><surname>Stocks</surname>
<given-names>K</given-names>
</name>
<name><surname>Allen</surname>
<given-names>EE</given-names>
</name>
<name><surname>Ellisman</surname>
<given-names>M</given-names>
</name>
<name><surname>Grethe</surname>
<given-names>J</given-names>
</name>
<name><surname>Wooley</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Community cyberinfrastructure for Advanced Microbial Ecology Research and Analysis: the CAMERA resource</article-title>
<source>Nucleic Acids Res</source>
<year>2011</year>
<volume>39</volume>
<fpage>D546</fpage>
<lpage>551</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkq1102</pub-id>
<pub-id pub-id-type="pmid">21045053</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname>
<given-names>SF</given-names>
</name>
<name><surname>Gish</surname>
<given-names>W</given-names>
</name>
<name><surname>Miller</surname>
<given-names>W</given-names>
</name>
<name><surname>Myers</surname>
<given-names>EW</given-names>
</name>
<name><surname>Lipman</surname>
<given-names>DJ</given-names>
</name>
</person-group>
<article-title>Basic local alignment search tool</article-title>
<source>J Mol Biol</source>
<year>1990</year>
<volume>215</volume>
<fpage>403</fpage>
<lpage>410</lpage>
<pub-id pub-id-type="doi">10.1016/S0022-2836(05)80360-2</pub-id>
<pub-id pub-id-type="pmid">2231712</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Giovannoni</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Tripp</surname>
<given-names>HJ</given-names>
</name>
<name><surname>Givan</surname>
<given-names>S</given-names>
</name>
<name><surname>Podar</surname>
<given-names>M</given-names>
</name>
<name><surname>Vergin</surname>
<given-names>KL</given-names>
</name>
<name><surname>Baptista</surname>
<given-names>D</given-names>
</name>
<name><surname>Bibbs</surname>
<given-names>L</given-names>
</name>
<name><surname>Eads</surname>
<given-names>J</given-names>
</name>
<name><surname>Richardson</surname>
<given-names>TH</given-names>
</name>
<name><surname>Noordewier</surname>
<given-names>M</given-names>
</name>
<name><surname>Rappé</surname>
<given-names>MS</given-names>
</name>
<name><surname>Short</surname>
<given-names>JM</given-names>
</name>
<name><surname>Carrington</surname>
<given-names>JC</given-names>
</name>
<name><surname>Mathur</surname>
<given-names>EJ</given-names>
</name>
</person-group>
<article-title>Genome streamlining in a cosmopolitan oceanic bacterium</article-title>
<source>Science</source>
<year>2005</year>
<volume>309</volume>
<fpage>1242</fpage>
<lpage>1245</lpage>
<pub-id pub-id-type="doi">10.1126/science.1114057</pub-id>
<pub-id pub-id-type="pmid">16109880</pub-id>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Burge</surname>
<given-names>SW</given-names>
</name>
<name><surname>Daub</surname>
<given-names>J</given-names>
</name>
<name><surname>Eberhardt</surname>
<given-names>R</given-names>
</name>
<name><surname>Tate</surname>
<given-names>J</given-names>
</name>
<name><surname>Barquist</surname>
<given-names>L</given-names>
</name>
<name><surname>Nawrocki</surname>
<given-names>EP</given-names>
</name>
<name><surname>Eddy</surname>
<given-names>SR</given-names>
</name>
<name><surname>Gardner</surname>
<given-names>PP</given-names>
</name>
<name><surname>Bateman</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Rfam 11.0: 10 years of RNA families</article-title>
<source>Nucleic Acids Res</source>
<year>2013</year>
<volume>41</volume>
<fpage>D226</fpage>
<lpage>D232</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gks1005</pub-id>
<pub-id pub-id-type="pmid">23125362</pub-id>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hyatt</surname>
<given-names>D</given-names>
</name>
<name><surname>Chen</surname>
<given-names>GL</given-names>
</name>
<name><surname>LoCascio</surname>
<given-names>PF</given-names>
</name>
<name><surname>Land</surname>
<given-names>ML</given-names>
</name>
<name><surname>Larimer</surname>
<given-names>FW</given-names>
</name>
<name><surname>Hauser</surname>
<given-names>LJ</given-names>
</name>
</person-group>
<article-title>Prodigal: prokaryotic gene recognition and translation initiation site identification</article-title>
<source>BMC Bioinformatics</source>
<year>2010</year>
<volume>11</volume>
<fpage>119</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2105-11-119</pub-id>
<pub-id pub-id-type="pmid">20211023</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/CyberinfraV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000168 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000168 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= CyberinfraV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:4287501 |texte= Adaptation to environmental factors shapes the organization of regulatory regions in microbial communities }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:25294412" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a CyberinfraV1
![]() | This area was generated with Dilib version V0.6.25. | ![]() |