Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin
Identifieur interne : 000637 ( Pmc/Corpus ); précédent : 000636; suivant : 000638Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin
Auteurs : Maria W. Smith ; Lydie Herfort ; Kaitlin Tyrol ; Dominic Suciu ; Victoria Campbell ; Byron C. Crump ; Tawnya D. Peterson ; Peter Zuber ; Antonio M. Baptista ; Holly M. SimonSource :
- PLoS ONE [ 1932-6203 ] ; 2010.
Abstract
Through their metabolic activities, microbial populations mediate the impact of high gradient regions on ecological function and productivity of the highly dynamic Columbia River coastal margin (CRCM). A 2226-probe oligonucleotide DNA microarray was developed to investigate expression patterns for microbial genes involved in nitrogen and carbon metabolism in the CRCM. Initial experiments with the environmental microarrays were directed toward validation of the platform and yielded high reproducibility in multiple tests. Bioinformatic and experimental validation also indicated that >85% of the microarray probes were specific for their corresponding target genes and for a few homologs within the same microbial family. The validated probe set was used to query gene expression responses by microbial assemblages to environmental variability. Sixty-four samples from the river, estuary, plume, and adjacent ocean were collected in different seasons and analyzed to correlate the measured variability in chemical, physical and biological water parameters to differences in global gene expression profiles. The method produced robust seasonal profiles corresponding to pre-freshet spring (April) and late summer (August). Overall relative gene expression was high in both seasons and was consistent with high microbial abundance measured by total RNA, heterotrophic bacterial production, and chlorophyll
Url:
DOI: 10.1371/journal.pone.0013312
PubMed: 20967204
PubMed Central: 2954162
Links to Exploration step
PMC:2954162Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin</title><author><name sortKey="Smith, Maria W" sort="Smith, Maria W" uniqKey="Smith M" first="Maria W." last="Smith">Maria W. Smith</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Herfort, Lydie" sort="Herfort, Lydie" uniqKey="Herfort L" first="Lydie" last="Herfort">Lydie Herfort</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Tyrol, Kaitlin" sort="Tyrol, Kaitlin" uniqKey="Tyrol K" first="Kaitlin" last="Tyrol">Kaitlin Tyrol</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Suciu, Dominic" sort="Suciu, Dominic" uniqKey="Suciu D" first="Dominic" last="Suciu">Dominic Suciu</name><affiliation><nlm:aff id="aff2"><addr-line>CombiMatrix Corporation, Mukilteo, Washington, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Campbell, Victoria" sort="Campbell, Victoria" uniqKey="Campbell V" first="Victoria" last="Campbell">Victoria Campbell</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Crump, Byron C" sort="Crump, Byron C" uniqKey="Crump B" first="Byron C." last="Crump">Byron C. Crump</name><affiliation><nlm:aff id="aff3"><addr-line>Horn Point Laboratory, University of Maryland Center for Environmental Science, Cambridge, Maryland, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Peterson, Tawnya D" sort="Peterson, Tawnya D" uniqKey="Peterson T" first="Tawnya D." last="Peterson">Tawnya D. Peterson</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Zuber, Peter" sort="Zuber, Peter" uniqKey="Zuber P" first="Peter" last="Zuber">Peter Zuber</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Baptista, Antonio M" sort="Baptista, Antonio M" uniqKey="Baptista A" first="Antonio M." last="Baptista">Antonio M. Baptista</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Simon, Holly M" sort="Simon, Holly M" uniqKey="Simon H" first="Holly M." last="Simon">Holly M. Simon</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">20967204</idno><idno type="pmc">2954162</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2954162</idno><idno type="RBID">PMC:2954162</idno><idno type="doi">10.1371/journal.pone.0013312</idno><date when="2010">2010</date><idno type="wicri:Area/Pmc/Corpus">000637</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin</title><author><name sortKey="Smith, Maria W" sort="Smith, Maria W" uniqKey="Smith M" first="Maria W." last="Smith">Maria W. Smith</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Herfort, Lydie" sort="Herfort, Lydie" uniqKey="Herfort L" first="Lydie" last="Herfort">Lydie Herfort</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Tyrol, Kaitlin" sort="Tyrol, Kaitlin" uniqKey="Tyrol K" first="Kaitlin" last="Tyrol">Kaitlin Tyrol</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Suciu, Dominic" sort="Suciu, Dominic" uniqKey="Suciu D" first="Dominic" last="Suciu">Dominic Suciu</name><affiliation><nlm:aff id="aff2"><addr-line>CombiMatrix Corporation, Mukilteo, Washington, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Campbell, Victoria" sort="Campbell, Victoria" uniqKey="Campbell V" first="Victoria" last="Campbell">Victoria Campbell</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Crump, Byron C" sort="Crump, Byron C" uniqKey="Crump B" first="Byron C." last="Crump">Byron C. Crump</name><affiliation><nlm:aff id="aff3"><addr-line>Horn Point Laboratory, University of Maryland Center for Environmental Science, Cambridge, Maryland, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Peterson, Tawnya D" sort="Peterson, Tawnya D" uniqKey="Peterson T" first="Tawnya D." last="Peterson">Tawnya D. Peterson</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Zuber, Peter" sort="Zuber, Peter" uniqKey="Zuber P" first="Peter" last="Zuber">Peter Zuber</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Baptista, Antonio M" sort="Baptista, Antonio M" uniqKey="Baptista A" first="Antonio M." last="Baptista">Antonio M. Baptista</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Simon, Holly M" sort="Simon, Holly M" uniqKey="Simon H" first="Holly M." last="Simon">Holly M. Simon</name><affiliation><nlm:aff id="aff1"><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2010">2010</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Through their metabolic activities, microbial populations mediate the impact of high gradient regions on ecological function and productivity of the highly dynamic Columbia River coastal margin (CRCM). A 2226-probe oligonucleotide DNA microarray was developed to investigate expression patterns for microbial genes involved in nitrogen and carbon metabolism in the CRCM. Initial experiments with the environmental microarrays were directed toward validation of the platform and yielded high reproducibility in multiple tests. Bioinformatic and experimental validation also indicated that >85% of the microarray probes were specific for their corresponding target genes and for a few homologs within the same microbial family. The validated probe set was used to query gene expression responses by microbial assemblages to environmental variability. Sixty-four samples from the river, estuary, plume, and adjacent ocean were collected in different seasons and analyzed to correlate the measured variability in chemical, physical and biological water parameters to differences in global gene expression profiles. The method produced robust seasonal profiles corresponding to pre-freshet spring (April) and late summer (August). Overall relative gene expression was high in both seasons and was consistent with high microbial abundance measured by total RNA, heterotrophic bacterial production, and chlorophyll <italic>a</italic>. Both seasonal patterns involved large numbers of genes that were highly expressed relative to background, yet each produced very different gene expression profiles. April patterns revealed high differential gene expression in the coastal margin samples (estuary, plume and adjacent ocean) relative to freshwater, while little differential gene expression was observed along the river-to-ocean transition in August. Microbial gene expression profiles appeared to relate, in part, to seasonal differences in nutrient availability and potential resource competition. Furthermore, our results suggest that highly-active particle-attached microbiota in the Columbia River water column may perform dissimilatory nitrate reduction (both dentrification and DNRA) within anoxic particle microniches.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Simenstad, Ca" uniqKey="Simenstad C">CA Simenstad</name></author><author><name sortKey="Small, Lf" uniqKey="Small L">LF Small</name></author><author><name sortKey="David Mcintire, C" uniqKey="David Mcintire C">C David McIntire</name></author><author><name sortKey="Jay, Da" uniqKey="Jay D">DA Jay</name></author><author><name sortKey="Sherwood, C" uniqKey="Sherwood C">C Sherwood</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hickey, B" uniqKey="Hickey B">B Hickey</name></author><author><name sortKey="Banas, N" uniqKey="Banas N">N Banas</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Barnes, C" uniqKey="Barnes C">C Barnes</name></author><author><name sortKey="Duxbury, Ac" uniqKey="Duxbury A">AC Duxbury</name></author><author><name sortKey="Morse, B" uniqKey="Morse B">B Morse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bruland, Kw" uniqKey="Bruland K">KW Bruland</name></author><author><name sortKey="Lohan, Mc" uniqKey="Lohan M">MC Lohan</name></author><author><name sortKey="Aguilar Islas, Am" uniqKey="Aguilar Islas A">AM Aguilar-Islas</name></author><author><name sortKey="Smith, Gj" uniqKey="Smith G">GJ Smith</name></author><author><name sortKey="Sohst, B" uniqKey="Sohst B">B Sohst</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jay, Da" uniqKey="Jay D">DA Jay</name></author><author><name sortKey="Giese, Bs" uniqKey="Giese B">BS Giese</name></author><author><name sortKey="Sherwood, Cr" uniqKey="Sherwood C">CR Sherwood</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jay, D" uniqKey="Jay D">D Jay</name></author><author><name sortKey="Uncles, R" uniqKey="Uncles R">R Uncles</name></author><author><name sortKey="Largeir, J" uniqKey="Largeir J">J Largeir</name></author><author><name sortKey="Geyer, W" uniqKey="Geyer W">W Geyer</name></author><author><name sortKey="Vallino, J" uniqKey="Vallino J">J Vallino</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chawla, A" uniqKey="Chawla A">A Chawla</name></author><author><name sortKey="Jay, D" uniqKey="Jay D">D Jay</name></author><author><name sortKey="Baptista, A" uniqKey="Baptista A">A Baptista</name></author><author><name sortKey="Wilkin, M" uniqKey="Wilkin M">M Wilkin</name></author><author><name sortKey="Seaton, C" uniqKey="Seaton C">C Seaton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Small, Lf" uniqKey="Small L">LF Small</name></author><author><name sortKey="Prahl, F" uniqKey="Prahl F">F Prahl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Giovannoni, S" uniqKey="Giovannoni S">S Giovannoni</name></author><author><name sortKey="Rappe, Ms" uniqKey="Rappe M">MS Rappe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Small, Lf" uniqKey="Small L">LF Small</name></author><author><name sortKey="Mcintire, Cd" uniqKey="Mcintire C">CD McIntire</name></author><author><name sortKey="Macdonald, Kb" uniqKey="Macdonald K">KB MacDonald</name></author><author><name sortKey="Lara Lara, Jr" uniqKey="Lara Lara J">JR Lara-Lara</name></author><author><name sortKey="Frey, Be" uniqKey="Frey B">BE Frey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baross, Ja" uniqKey="Baross J">JA Baross</name></author><author><name sortKey="Crump, Bc" uniqKey="Crump B">BC Crump</name></author><author><name sortKey="Simenstad, Ca" uniqKey="Simenstad C">CA Simenstad</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Prahl, F" uniqKey="Prahl F">F Prahl</name></author><author><name sortKey="Small, Lf" uniqKey="Small L">LF Small</name></author><author><name sortKey="Eversmeyer, B" uniqKey="Eversmeyer B">B Eversmeyer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crump, Bc" uniqKey="Crump B">BC Crump</name></author><author><name sortKey="Baross, Ja" uniqKey="Baross J">JA Baross</name></author><author><name sortKey="Simenstad, Ca" uniqKey="Simenstad C">CA Simenstad</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crump, Bc" uniqKey="Crump B">BC Crump</name></author><author><name sortKey="Armbrust, Ev" uniqKey="Armbrust E">EV Armbrust</name></author><author><name sortKey="Baross, Ja" uniqKey="Baross J">JA Baross</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crump, Bc" uniqKey="Crump B">BC Crump</name></author><author><name sortKey="Baross, Ja" uniqKey="Baross J">JA Baross</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Griffiths, Ri" uniqKey="Griffiths R">RI Griffiths</name></author><author><name sortKey="Whiteley, As" uniqKey="Whiteley A">AS Whiteley</name></author><author><name sortKey="O Donnell, Ag" uniqKey="O Donnell A">AG O'Donnell</name></author><author><name sortKey="Bailey, Mj" uniqKey="Bailey M">MJ Bailey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Grant, S" uniqKey="Grant S">S Grant</name></author><author><name sortKey="Grant, Wd" uniqKey="Grant W">WD Grant</name></author><author><name sortKey="Cowan, Da" uniqKey="Cowan D">DA Cowan</name></author><author><name sortKey="Jones, Be" uniqKey="Jones B">BE Jones</name></author><author><name sortKey="Ma, Y" uniqKey="Ma Y">Y Ma</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Frias Lopez, J" uniqKey="Frias Lopez J">J Frias-Lopez</name></author><author><name sortKey="Shi, Y" uniqKey="Shi Y">Y Shi</name></author><author><name sortKey="Tyson, Gw" uniqKey="Tyson G">GW Tyson</name></author><author><name sortKey="Coleman, Ml" uniqKey="Coleman M">ML Coleman</name></author><author><name sortKey="Schuster, Sc" uniqKey="Schuster S">SC Schuster</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></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, J" uniqKey="Zhou J">J Zhou</name></author><author><name sortKey="Thompson, Dk" uniqKey="Thompson D">DK Thompson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Taroncher Oldenburg, G" uniqKey="Taroncher Oldenburg G">G Taroncher-Oldenburg</name></author><author><name sortKey="Griner, Em" uniqKey="Griner E">EM Griner</name></author><author><name sortKey="Francis, Ca" uniqKey="Francis C">CA Francis</name></author><author><name sortKey="Ward, Bb" uniqKey="Ward B">BB Ward</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Parro, V" uniqKey="Parro V">V Parro</name></author><author><name sortKey="Moreno Paz, M" uniqKey="Moreno Paz M">M Moreno-Paz</name></author><author><name sortKey="Gonzalez Toril, E" uniqKey="Gonzalez Toril E">E Gonzalez-Toril</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lodes, Mj" uniqKey="Lodes M">MJ Lodes</name></author><author><name sortKey="Suciu, D" uniqKey="Suciu D">D Suciu</name></author><author><name sortKey="Elliott, M" uniqKey="Elliott M">M Elliott</name></author><author><name sortKey="Stover, Ag" uniqKey="Stover A">AG Stover</name></author><author><name sortKey="Ross, M" uniqKey="Ross M">M Ross</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ghindilis, Al" uniqKey="Ghindilis A">AL Ghindilis</name></author><author><name sortKey="Smith, Mw" uniqKey="Smith M">MW Smith</name></author><author><name sortKey="Schwarzkopf, Kr" uniqKey="Schwarzkopf K">KR Schwarzkopf</name></author><author><name sortKey="Roth, Km" uniqKey="Roth K">KM Roth</name></author><author><name sortKey="Peyvan, K" uniqKey="Peyvan K">K Peyvan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bolotin, S" uniqKey="Bolotin S">S Bolotin</name></author><author><name sortKey="Lombos, E" uniqKey="Lombos E">E Lombos</name></author><author><name sortKey="Yeung, R" uniqKey="Yeung R">R Yeung</name></author><author><name sortKey="Eshaghi, A" uniqKey="Eshaghi A">A Eshaghi</name></author><author><name sortKey="Blair, J" uniqKey="Blair J">J Blair</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="He, S" uniqKey="He S">S He</name></author><author><name sortKey="Kunin, V" uniqKey="Kunin V">V Kunin</name></author><author><name sortKey="Haynes, M" uniqKey="Haynes M">M Haynes</name></author><author><name sortKey="Martin, Hg" uniqKey="Martin H">HG Martin</name></author><author><name sortKey="Ivanova, N" uniqKey="Ivanova N">N Ivanova</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Rh" uniqKey="Liu R">RH Liu</name></author><author><name sortKey="Munro, Sb" uniqKey="Munro S">SB Munro</name></author><author><name sortKey="Nguyen, T" uniqKey="Nguyen T">T Nguyen</name></author><author><name sortKey="Siuda, T" uniqKey="Siuda T">T Siuda</name></author><author><name sortKey="Suciu, D" uniqKey="Suciu D">D Suciu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Markowitz, Vm" uniqKey="Markowitz V">VM Markowitz</name></author><author><name sortKey="Szeto, E" uniqKey="Szeto E">E Szeto</name></author><author><name sortKey="Palaniappan, K" uniqKey="Palaniappan K">K Palaniappan</name></author><author><name sortKey="Grechkin, Y" uniqKey="Grechkin Y">Y Grechkin</name></author><author><name sortKey="Chu, K" uniqKey="Chu K">K Chu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Aguilar Islas, Am" uniqKey="Aguilar Islas A">AM Aguilar-Islas</name></author><author><name sortKey="Bruland, Kw" uniqKey="Bruland K">KW Bruland</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kudela, Rm" uniqKey="Kudela R">RM Kudela</name></author><author><name sortKey="Peterson, Td" uniqKey="Peterson T">TD Peterson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turner, Re" uniqKey="Turner R">RE Turner</name></author><author><name sortKey="Rabalais, Nn" uniqKey="Rabalais N">NN Rabalais</name></author><author><name sortKey="Justic, D" uniqKey="Justic D">D Justic</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Colbert, D" uniqKey="Colbert D">D Colbert</name></author><author><name sortKey="Mcmanus, J" uniqKey="Mcmanus J">J McManus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brown, C" uniqKey="Brown C">C Brown</name></author><author><name sortKey="Ozretich, R" uniqKey="Ozretich R">R Ozretich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Arrigo, Kr" uniqKey="Arrigo K">KR Arrigo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Conomos, T" uniqKey="Conomos T">T Conomos</name></author><author><name sortKey="Gross, Mg" uniqKey="Gross M">MG Gross</name></author><author><name sortKey="Barnes, Ca" uniqKey="Barnes C">CA Barnes</name></author><author><name sortKey="Richards, Fa" uniqKey="Richards F">FA Richards</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Arhonditsis, G" uniqKey="Arhonditsis G">G Arhonditsis</name></author><author><name sortKey="Brett, Mt" uniqKey="Brett M">MT Brett</name></author><author><name sortKey="Frodge, J" uniqKey="Frodge J">J Frodge</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sala, Mm" uniqKey="Sala M">MM Sala</name></author><author><name sortKey="Terrado, R" uniqKey="Terrado R">R Terrado</name></author><author><name sortKey="Lovejoy, C" uniqKey="Lovejoy C">C Lovejoy</name></author><author><name sortKey="Unrein, F" uniqKey="Unrein F">F Unrein</name></author><author><name sortKey="Pedros Alio, C" uniqKey="Pedros Alio C">C Pedros-Alio</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Teira, E" uniqKey="Teira E">E Teira</name></author><author><name sortKey="Gasol, Jm" uniqKey="Gasol J">JM Gasol</name></author><author><name sortKey="Aranguren Gassis, M" uniqKey="Aranguren Gassis M">M Aranguren-Gassis</name></author><author><name sortKey="Fernandez, A" uniqKey="Fernandez A">A Fernandez</name></author><author><name sortKey="Gonzalez, J" uniqKey="Gonzalez J">J Gonzalez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crump, Bc" uniqKey="Crump B">BC Crump</name></author><author><name sortKey="Hopkinson, Cs" uniqKey="Hopkinson C">CS Hopkinson</name></author><author><name sortKey="Sogin, Ml" uniqKey="Sogin M">ML Sogin</name></author><author><name sortKey="Hobbie, Je" uniqKey="Hobbie J">JE Hobbie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lodes, Mj" uniqKey="Lodes M">MJ Lodes</name></author><author><name sortKey="Caraballo, M" uniqKey="Caraballo M">M Caraballo</name></author><author><name sortKey="Suciu, D" uniqKey="Suciu D">D Suciu</name></author><author><name sortKey="Munro, S" uniqKey="Munro S">S Munro</name></author><author><name sortKey="Kumar, A" uniqKey="Kumar A">A Kumar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lodes, Mj" uniqKey="Lodes M">MJ Lodes</name></author><author><name sortKey="Suciu, D" uniqKey="Suciu D">D Suciu</name></author><author><name sortKey="Wilmoth, Jl" uniqKey="Wilmoth J">JL Wilmoth</name></author><author><name sortKey="Ross, M" uniqKey="Ross M">M Ross</name></author><author><name sortKey="Munro, S" uniqKey="Munro S">S Munro</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maurer, K" uniqKey="Maurer K">K Maurer</name></author><author><name sortKey="Cooper, J" uniqKey="Cooper J">J Cooper</name></author><author><name sortKey="Caraballo, M" uniqKey="Caraballo M">M Caraballo</name></author><author><name sortKey="Crye, J" uniqKey="Crye J">J Crye</name></author><author><name sortKey="Suciu, D" uniqKey="Suciu D">D Suciu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Allawi, Ht" uniqKey="Allawi H">HT Allawi</name></author><author><name sortKey="Santalucia, J" uniqKey="Santalucia J">J SantaLucia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Allawi, Ht" uniqKey="Allawi H">HT Allawi</name></author><author><name sortKey="Santalucia, J" uniqKey="Santalucia J">J SantaLucia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Allawi, Ht" uniqKey="Allawi H">HT Allawi</name></author><author><name sortKey="Santalucia, J" uniqKey="Santalucia J">J SantaLucia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Seshadri, R" uniqKey="Seshadri R">R Seshadri</name></author><author><name sortKey="Kravitz, Sa" uniqKey="Kravitz S">SA Kravitz</name></author><author><name sortKey="Smarr, L" uniqKey="Smarr L">L Smarr</name></author><author><name sortKey="Gilna, P" uniqKey="Gilna P">P Gilna</name></author><author><name sortKey="Frazier, M" uniqKey="Frazier M">M Frazier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Delong, Ef" uniqKey="Delong E">EF DeLong</name></author><author><name sortKey="Preston, Cm" uniqKey="Preston C">CM Preston</name></author><author><name sortKey="Mincer, T" uniqKey="Mincer T">T Mincer</name></author><author><name sortKey="Rich, V" uniqKey="Rich V">V Rich</name></author><author><name sortKey="Hallam, Sj" uniqKey="Hallam S">SJ Hallam</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rusch, Db" uniqKey="Rusch D">DB Rusch</name></author><author><name sortKey="Halpern, Al" uniqKey="Halpern A">AL Halpern</name></author><author><name sortKey="Sutton, G" uniqKey="Sutton G">G Sutton</name></author><author><name sortKey="Heidelberg, Kb" uniqKey="Heidelberg K">KB Heidelberg</name></author><author><name sortKey="Williamson, S" uniqKey="Williamson S">S Williamson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kan, J" uniqKey="Kan J">J Kan</name></author><author><name sortKey="Suzuki, Mt" uniqKey="Suzuki M">MT Suzuki</name></author><author><name sortKey="Wang, K" uniqKey="Wang K">K Wang</name></author><author><name sortKey="Evans, Se" uniqKey="Evans S">SE Evans</name></author><author><name sortKey="Chen, F" uniqKey="Chen F">F Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Glibert, Pm" uniqKey="Glibert P">PM Glibert</name></author><author><name sortKey="Garside, C" uniqKey="Garside C">C Garside</name></author><author><name sortKey="Fuhrman, Ja" uniqKey="Fuhrman J">JA Fuhrman</name></author><author><name sortKey="Roman, Mr" uniqKey="Roman M">MR Roman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tiedje, J" uniqKey="Tiedje J">J Tiedje</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robertson, La" uniqKey="Robertson L">LA Robertson</name></author><author><name sortKey="Dalsgaar, T" uniqKey="Dalsgaar T">T Dalsgaar</name></author><author><name sortKey="Revsbech, Np" uniqKey="Revsbech N">NP Revsbech</name></author><author><name sortKey="Kuenen, Jg" uniqKey="Kuenen J">JG Kuenen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robertson, La" uniqKey="Robertson L">LA Robertson</name></author><author><name sortKey="Van Niel, Ewj" uniqKey="Van Niel E">EWJ Van Niel</name></author><author><name sortKey="Torremans, Ram" uniqKey="Torremans R">RAM Torremans</name></author><author><name sortKey="Kuenen, Jg" uniqKey="Kuenen J">JG Kuenen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Takaya, N" uniqKey="Takaya N">N Takaya</name></author><author><name sortKey="Catalan Sakairi, Mab" uniqKey="Catalan Sakairi M">MAB Catalan-Sakairi</name></author><author><name sortKey="Sakaguchi, Y" uniqKey="Sakaguchi Y">Y Sakaguchi</name></author><author><name sortKey="Kato, I" uniqKey="Kato I">I Kato</name></author><author><name sortKey="Zhou, Z" uniqKey="Zhou Z">Z Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Michotey, V" uniqKey="Michotey V">V Michotey</name></author><author><name sortKey="Bonin, P" uniqKey="Bonin P">P Bonin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Omnes, P" uniqKey="Omnes P">P Omnes</name></author><author><name sortKey="Slawyk, G" uniqKey="Slawyk G">G Slawyk</name></author><author><name sortKey="Garcia, N" uniqKey="Garcia N">N Garcia</name></author><author><name sortKey="Bonin, P" uniqKey="Bonin P">P Bonin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="John, De" uniqKey="John D">DE John</name></author><author><name sortKey="Wang, Za" uniqKey="Wang Z">ZA Wang</name></author><author><name sortKey="Liu, X" uniqKey="Liu X">X Liu</name></author><author><name sortKey="Byrne, Rh" uniqKey="Byrne R">RH Byrne</name></author><author><name sortKey="Corredor, Je" uniqKey="Corredor J">JE Corredor</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Treusch, Ah" uniqKey="Treusch A">AH Treusch</name></author><author><name sortKey="Vergin, Kl" uniqKey="Vergin K">KL Vergin</name></author><author><name sortKey="Finlay, La" uniqKey="Finlay L">LA Finlay</name></author><author><name sortKey="Donatz, Mg" uniqKey="Donatz M">MG Donatz</name></author><author><name sortKey="Burton, Rm" uniqKey="Burton R">RM Burton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dubey, Ak" uniqKey="Dubey A">AK Dubey</name></author><author><name sortKey="Baker, Cs" uniqKey="Baker C">CS Baker</name></author><author><name sortKey="Suzuki, K" uniqKey="Suzuki K">K Suzuki</name></author><author><name sortKey="Jones, Ad" uniqKey="Jones A">AD Jones</name></author><author><name sortKey="Pandit, P" uniqKey="Pandit P">P Pandit</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schultz, Je" uniqKey="Schultz J">JE Schultz</name></author><author><name sortKey="Matin, A" uniqKey="Matin A">A Matin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bronk, Da" uniqKey="Bronk D">DA Bronk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ghindilis, Al" uniqKey="Ghindilis A">AL Ghindilis</name></author><author><name sortKey="Smith, Mw" uniqKey="Smith M">MW Smith</name></author><author><name sortKey="Schwarzkopf, Kr" uniqKey="Schwarzkopf K">KR Schwarzkopf</name></author><author><name sortKey="Zhan, C" uniqKey="Zhan C">C Zhan</name></author><author><name sortKey="Evans, Dr" uniqKey="Evans D">DR Evans</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Preston, Cm" uniqKey="Preston C">CM Preston</name></author><author><name sortKey="Marin, R" uniqKey="Marin R">R Marin</name></author><author><name sortKey="Jensen, Sd" uniqKey="Jensen S">SD Jensen</name></author><author><name sortKey="Feldman, J" uniqKey="Feldman J">J Feldman</name></author><author><name sortKey="Birch, Jm" uniqKey="Birch J">JM Birch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Parsons, T" uniqKey="Parsons T">T Parsons</name></author><author><name sortKey="Maita, Y" uniqKey="Maita Y">Y Maita</name></author><author><name sortKey="Lalli, Cm" uniqKey="Lalli C">CM Lalli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sugimura, Y" uniqKey="Sugimura Y">Y Sugimura</name></author><author><name sortKey="Suzuki, Y" uniqKey="Suzuki Y">Y Suzuki</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Verardo, Dj" uniqKey="Verardo D">DJ Verardo</name></author><author><name sortKey="Froelich, Pn" uniqKey="Froelich P">PN Froelich</name></author><author><name sortKey="Mcintyre, A" uniqKey="Mcintyre A">A McIntyre</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Quackenbush, J" uniqKey="Quackenbush J">J Quackenbush</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schena, M" uniqKey="Schena M">M Schena</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Eisen, Mb" uniqKey="Eisen M">MB Eisen</name></author><author><name sortKey="Spellman, Pt" uniqKey="Spellman P">PT Spellman</name></author><author><name sortKey="Brown, Po" uniqKey="Brown P">PO Brown</name></author><author><name sortKey="Botstein, D" uniqKey="Botstein D">D Botstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edgar, R" uniqKey="Edgar R">R Edgar</name></author><author><name sortKey="Domrachev, M" uniqKey="Domrachev M">M Domrachev</name></author><author><name sortKey="Lash, Ae" uniqKey="Lash A">AE Lash</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">20967204</article-id><article-id pub-id-type="pmc">2954162</article-id><article-id pub-id-type="publisher-id">10-PONE-RA-20286R1</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0013312</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline"><subject>Genetics and Genomics/Gene Expression</subject><subject>Microbiology/Environmental Microbiology</subject><subject>Marine and Aquatic Sciences/Biological Oceanography</subject></subj-group></article-categories><title-group><article-title>Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin</article-title><alt-title alt-title-type="running-head">Columbia River Gene Expression</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Smith</surname><given-names>Maria W.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Herfort</surname><given-names>Lydie</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Tyrol</surname><given-names>Kaitlin</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Suciu</surname><given-names>Dominic</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Campbell</surname><given-names>Victoria</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Crump</surname><given-names>Byron C.</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Peterson</surname><given-names>Tawnya D.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Zuber</surname><given-names>Peter</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Baptista</surname><given-names>Antonio M.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Simon</surname><given-names>Holly M.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Division of Environmental & Biomolecular Systems, Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Beaverton, Oregon, United States of America</addr-line></aff><aff id="aff2"><label>2</label><addr-line>CombiMatrix Corporation, Mukilteo, Washington, United States of America</addr-line></aff><aff id="aff3"><label>3</label><addr-line>Horn Point Laboratory, University of Maryland Center for Environmental Science, Cambridge, Maryland, United States of America</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Gilbert</surname><given-names>Jack Anthony</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1">Argonne National Laboratory, United States of America</aff><author-notes><corresp id="cor1">* E-mail: <email>simonh@ebs.ogi.edu</email></corresp><fn fn-type="con"><p>Conceived and designed the experiments: MWS LH DS BCC PZ AMB HMS. Performed the experiments: MWS KT VC. Analyzed the data: MWS LH KT DS VC BCC TDP PZ AMB HMS. Contributed reagents/materials/analysis tools: MWS LH KT DS VC BCC TDP PZ AMB HMS. Wrote the paper: MWS LH DS BCC TDP PZ AMB HMS. Cruise management, sample collection: LH BCC TDP. Bioinformatics support for microarray probe design: DS.</p></fn></author-notes><pub-date pub-type="collection"><year>2010</year></pub-date><pub-date pub-type="epub"><day>13</day><month>10</month><year>2010</year></pub-date><volume>5</volume><issue>10</issue><elocation-id>e13312</elocation-id><history><date date-type="received"><day>21</day><month>6</month><year>2010</year></date><date date-type="accepted"><day>10</day><month>9</month><year>2010</year></date></history><permissions><copyright-statement>Smith et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</copyright-statement></permissions><abstract><p>Through their metabolic activities, microbial populations mediate the impact of high gradient regions on ecological function and productivity of the highly dynamic Columbia River coastal margin (CRCM). A 2226-probe oligonucleotide DNA microarray was developed to investigate expression patterns for microbial genes involved in nitrogen and carbon metabolism in the CRCM. Initial experiments with the environmental microarrays were directed toward validation of the platform and yielded high reproducibility in multiple tests. Bioinformatic and experimental validation also indicated that >85% of the microarray probes were specific for their corresponding target genes and for a few homologs within the same microbial family. The validated probe set was used to query gene expression responses by microbial assemblages to environmental variability. Sixty-four samples from the river, estuary, plume, and adjacent ocean were collected in different seasons and analyzed to correlate the measured variability in chemical, physical and biological water parameters to differences in global gene expression profiles. The method produced robust seasonal profiles corresponding to pre-freshet spring (April) and late summer (August). Overall relative gene expression was high in both seasons and was consistent with high microbial abundance measured by total RNA, heterotrophic bacterial production, and chlorophyll <italic>a</italic>. Both seasonal patterns involved large numbers of genes that were highly expressed relative to background, yet each produced very different gene expression profiles. April patterns revealed high differential gene expression in the coastal margin samples (estuary, plume and adjacent ocean) relative to freshwater, while little differential gene expression was observed along the river-to-ocean transition in August. Microbial gene expression profiles appeared to relate, in part, to seasonal differences in nutrient availability and potential resource competition. Furthermore, our results suggest that highly-active particle-attached microbiota in the Columbia River water column may perform dissimilatory nitrate reduction (both dentrification and DNRA) within anoxic particle microniches.</p></abstract><counts><page-count count="19"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>The Columbia River and its tributaries, with a drainage basin of 660,480 km<sup>2</sup>, represent the second-largest freshwater discharge in the United States <xref ref-type="bibr" rid="pone.0013312-Simenstad1">[1]</xref>, <xref ref-type="bibr" rid="pone.0013312-Hickey1">[2]</xref>. Due largely to snow melt in the spring, discharge fluctuates seasonally from the highest volumes in April through June to the lowest in September and October. The Columbia River has a profound influence on biogeochemical processes in the coastal ocean through the delivery of nutrients in a massive plume that during times of high discharge reaches dozens to hundreds of kilometers from the river mouth <xref ref-type="bibr" rid="pone.0013312-Hickey1">[2]</xref>–<xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>. We refer in this paper to the Columbia River coastal margin (CRCM) as the continuum between the river, the estuary, the plume and the host continental shelf of the Eastern North Pacific Ocean.</p><p>The Columbia River estuary is characterized by a strong tidal cycle, high turbidity, and vertical stratification varying in strength with the tides and river discharge <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>–<xref ref-type="bibr" rid="pone.0013312-Chawla1">[7]</xref>. Estuarine waters contain biological and chemical gradients established by the mixing of freshwater and seawater <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>, <xref ref-type="bibr" rid="pone.0013312-Small1">[8]</xref> that are thought to deeply influence the composition of natural bacterioplankton communities <xref ref-type="bibr" rid="pone.0013312-Giovannoni1">[9]</xref>. A productive detrital food web is driven by heterotrophic bacterioplankton, and is based on allochthonous organic material and freshwater phytoplankton that develop seasonally in river impoundments <xref ref-type="bibr" rid="pone.0013312-Small2">[10]</xref>. Tidally-driven estuarine turbidity maxima (ETM) events trap and re-suspend both mineral and organic particles transported through the estuary and extend their residence time several-fold <xref ref-type="bibr" rid="pone.0013312-Small1">[8]</xref>, <xref ref-type="bibr" rid="pone.0013312-Baross1">[11]</xref>–<xref ref-type="bibr" rid="pone.0013312-Crump1">[13]</xref>. As a result, particle-attached bacteria trapped in the estuary by ETM account for approximately 90% of the heterotrophic bacterial activity in the water column <xref ref-type="bibr" rid="pone.0013312-Crump1">[13]</xref>. Thus, microbial assemblages are influenced not only by water chemistry, which varies temporally as a function of season, oceanographic conditions, tidal phase and river discharge <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>, but also by transport from both the river and adjacent coastal ocean.</p><p>Phylogenetic diversity of the estuarine microbial community has been analyzed for both bacteria <xref ref-type="bibr" rid="pone.0013312-Crump2">[14]</xref> and archaea <xref ref-type="bibr" rid="pone.0013312-Crump3">[15]</xref> using 16S rRNA gene sequence libraries. The majority of free-living estuarine bacterial sequences were found to be similar to tidal freshwater and adjacent coastal ocean environmental clones, while the particle–attached fraction contained a higher proportion (up to 75%) of uniquely estuarine clones <xref ref-type="bibr" rid="pone.0013312-Crump2">[14]</xref>. While this particle-attached fraction is known to be more active than the free-living bacterioplankton <xref ref-type="bibr" rid="pone.0013312-Crump1">[13]</xref>, nothing is known about their specific metabolic properties.</p><p>Recent developments in molecular biological tools permit investigations of metabolic responses to fluctuating conditions without a requirement for microbial cultivation; for example, through analysis of gene expression activity. Improvements in protocols for RNA isolation and in techniques for the removal of ribosomal RNA yielded high-quality enrichments of mRNA from environmental samples in recent studies <xref ref-type="bibr" rid="pone.0013312-Griffiths1">[16]</xref>–<xref ref-type="bibr" rid="pone.0013312-Poretsky1">[19]</xref>, and subsequent large-scale sequence analysis revealed highly expressed genes <xref ref-type="bibr" rid="pone.0013312-FriasLopez1">[18]</xref>. However, the associated high cost of sequencing still limits the number of samples that can be analyzed in this manner, and, therefore, also the spatial and temporal resolution that can be achieved. Alternatively, microarray-based technologies allow for the analysis of a relatively large number of samples at a fraction of the cost of sequencing. DNA microarrays are, therefore, compelling functional tools for investigating microbial transcriptional activity in the environment <xref ref-type="bibr" rid="pone.0013312-Zhou1">[20]</xref>–<xref ref-type="bibr" rid="pone.0013312-Parro1">[22]</xref>. However, for successful microarray applications, several important issues need to be addressed, including (i) selection of representative sets of annotated functional genes, and (ii) specificity of microarray probe hybridization within the context of a large pool of unknown environmental mRNA <xref ref-type="bibr" rid="pone.0013312-Parro1">[22]</xref>.</p><p>We used DNA oligonucleotide microarrays to assess differential gene expression in bacterial and archaeal populations in response to physical and chemical gradients in the CRCM. Freshwater influx from the Columbia River and coastal upwelling from the continental shelf influence biogeochemical cycles by creating large physical and chemical gradients in the water column <xref ref-type="bibr" rid="pone.0013312-Prahl1">[12]</xref>. This study was undertaken as a first step toward understanding how microbial populations respond, at the transcriptional level, to such gradients.</p><p>For this work we employed the CombiMatrix CustomArray™ format, which, in recent studies, produced high quality gene expression data for a variety of complex clinical and field samples <xref ref-type="bibr" rid="pone.0013312-Parro1">[22]</xref>–<xref ref-type="bibr" rid="pone.0013312-Liu1">[27]</xref>. Using this system, we designed probes for a set of several thousand genes from bacteria and archaea that were selected from functional annotations in the Integrated Microbial Genomes (IMG) data management system of the U.S. Department of Energy Joint Genome Institute (DOE JGI; <xref ref-type="bibr" rid="pone.0013312-Markowitz1">[28]</xref>). Our results indicate that bacterial and archaeal gene expression in the CRCM varied primarily with seasonality in environmental characteristics. Trends involving light, phytoplankton biomass, proximity to the river mouth, and availability of nitrate, phosphate, and dissolved organic carbon (DOC) were observed in the gene expression profiles. An association of gene expression patterns with habitat (tidal freshwater, estuary, plume, or coastal) was also observed, but could not be accounted for by salinity differences alone. The existence of stable, repeating gene expression patterns over temporal and spatial gradients suggests that expression of microbial genes may vary in predictable ways, and, therefore, may be a useful indicator of environmental change. Overall, our results indicate that the highly dynamic CRCM, which is the focus of an NSF Science and Technology Center for Coastal Margin Observation and Prediction (CMOP), is an ideal testing ground for new tools aimed at elucidating microbial metabolic responses to changing environmental conditions.</p></sec><sec id="s2"><title>Results and Discussion</title><sec id="s2a"><title>Geochemical characterization of CRCM water samples</title><p>Our data set consisted of 110 water samples that were collected along transects of the tidal freshwater, estuary, plume and adjacent coastal ocean (<xref ref-type="fig" rid="pone-0013312-g001">Fig. 1A</xref>) on ship-based field campaigns. Campaigns were conducted from 7 to 14 days during summer with low river discharge (August 2007), late fall (November 2007), pre-freshet early spring (April 2008), and spring freshet at record high river discharge (June 2008) (<xref ref-type="fig" rid="pone-0013312-g001">Fig 1B</xref>). A detailed analysis of the properties of CRCM waters was undertaken to provide environmental context for the gene expression data. Physical and chemical data are plotted in <xref ref-type="fig" rid="pone-0013312-g002">Fig. 2</xref>. The X-axis indicates samples sorted (i) from left to right according to year and month starting from August 2007; (ii) then within each month by location, starting with the adjacent ocean through the plume and estuary to the tidal freshwater; and (iii) within each location, according to salinity, from high (left) to low (right). Temperature measurements showed the typical seasonal pattern, with a summer high of 20°C in the river in August, and an average of 8–10°C in spring and fall (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2A</xref>). Oxygen concentrations varied between 7–15 mg/L, and in the river and estuary they were inversely correlated with salinity (<italic>R<sup>2</sup></italic> of −0.86 in August, and −0.92 in both April and June). Thus, oxygen concentrations were higher in freshwater than in the estuary in all months sampled, except November. (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2A, 2B</xref>). This observation may be explained, in part, by release of oxygen during photosynthesis and growth by freshwater phytoplankton in the river.</p><fig id="pone-0013312-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g001</object-id><label>Figure 1</label><caption><title>Characterization of the Columbia River coastal margin.</title><p>(<bold>A</bold>) Sampling locations (indicated as stars) along the CRCM transects. (<bold>B</bold>) Columbia River daily water discharge (dashed and dotted lines indicate minimum and maximum values, respectively), and estuary residence times (solid line), calculated based on U.S. Geologic Survey's National Stream Water Quality Network data collected at Bonneville Dam. Residence times are defined as R = V/(86400*Q), where R is the residence time in days, V is the volume of the estuary in m<sup>3</sup> relative to mean sea level (defined from Beaver Army to the mouth), and Q is the river discharge in m<sup>3</sup>s<sup>−1</sup>. The data correspond to the 2-week long sample collection times for each month.</p></caption><graphic xlink:href="pone.0013312.g001"></graphic></fig><fig id="pone-0013312-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g002</object-id><label>Figure 2</label><caption><title>Selected physical and chemical characteristics corresponding to the seasonal sample sets collected for microarray analysis.</title><p>In all graphs, the X-axis indicates samples sorted first from left to right, according to year and month starting from August and November, 2007, then April and June, 2008; second, within each month by location, starting with the coastal ocean through the plume and estuary to the river; and third, according to salinity within each location, from high (left) to low (right). (<bold>A</bold>) Temperature (black line) and oxygen concentration (gray line); (<bold>B</bold>) Silicic acid concentration (black line) and salinity (gray line); (<bold>C</bold>) and (<bold>D</bold>), nitrate and phosphate concentrations, respectively. Black vertical lines divide seasons, dashed vertical lines within each season divide the ocean (O) and plume (P) from estuary (E) and river (R) samples. The months in which sampling occurred are indicated below.</p></caption><graphic xlink:href="pone.0013312.g002"></graphic></fig></sec><sec id="s2b"><title>Nutrients</title><p>Seasonally-fluctuating discharge volumes in the Columbia River are accompanied by changes in concentrations of macronutrients and organic detritus. Among major rivers, the Columbia River is unusual in being silicic acid rich and nitrate poor during the summer months, with nitrate concentrations of only 2–10 µM in June and July <xref ref-type="bibr" rid="pone.0013312-AguilarIslas1">[29]</xref>, <xref ref-type="bibr" rid="pone.0013312-Kudela1">[30]</xref>. This is in contrast, for example, to the Mississippi River, in which summertime nitrate concentrations have been measured at >100 µM <xref ref-type="bibr" rid="pone.0013312-Turner1">[31]</xref>. Our measurements were consistent with previous monthly analyses of nitrate and silicic acid at the Beaver Army Terminal station (River Mile 53) by the US Geological Survey over the last decade <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>. The silicic acid concentration in the river was high during all seasons and was inversely correlated (<italic>R<sup>2</sup></italic> = −0.93) with water salinity (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2B</xref>), similar to previous observations <xref ref-type="bibr" rid="pone.0013312-AguilarIslas1">[29]</xref>. Typical summer nitrate concentrations (∼7 µM) were observed for river water samples collected in August 2007 (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2C</xref>). In winter, nitrate concentrations can be as high as 50 µM, coincident with high winter rainfall and high flow from coastal tributaries <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>, <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>. Nitrate concentrations in November 2007 and April 2008 were typical of the ‘shoulder seasons’ before and after the winter nitrate peak, with the nitrate concentrations elevated approximately 3X above the summer levels (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2C</xref>). In April, November and August, the plume and coastal ocean samples had nitrate concentrations that were much lower than those in the tidal freshwater and estuary, consistent with the idea that nitrate flowing into the lower Columbia River from the Columbia basin watershed is typically used within the estuary <xref ref-type="bibr" rid="pone.0013312-Hickey1">[2]</xref>. Several of the August samples were apparently influenced by upwelling; the two most prominent of these are indicated with arrowheads in <xref ref-type="fig" rid="pone-0013312-g002">Fig. 2A</xref>. They displayed clear upwelling characteristics such as reduced water temperature, low oxygen concentration, and highly elevated levels of nitrate and phosphate (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2C and 2D</xref>, respectively). In contrast to the other sampling times, June nitrate levels were similar in the river, plume and a number of ocean samples, consistent with the observation that elevated nitrate concentrations from the watershed may be delivered to the ocean by high riverflow <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>. Finally, ammonium concentrations varied from 1 to 3 µM in the majority of samples, and did not show clear seasonal patterns (<xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>).</p><p>Phosphate concentrations in the CRCM also varied seasonally (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2D</xref>), with the highest levels in the adjacent coastal ocean in August, and in the plume and estuary in June. The higher August levels appeared to be at least partly explained by upwelling, while the higher June measurements were likely due to elevated river discharge. A similar dynamic was described for the Tillamook Bay estuary (Oregon, USA), for which elevated estuarine flushing rates produced by high river discharge limited nutrient uptake by phyto- and bacterioplankton, resulting in enhanced nutrient delivery to the coastal ocean <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>. On the other hand, reduced phosphate concentrations were observed in the tidal freshwater and estuary in April with the average N∶P ratio of 131∶1. In late April to early May of 2002 a similar observation was made in the upper portions of the Yaquina Bay estuary, with the N∶P ratio reaching as high as 176∶1 <xref ref-type="bibr" rid="pone.0013312-Brown1">[33]</xref>. We speculate that the apparent phosphate depletion might have been caused by an abundant phytoplankton community in the Columbia River in April (see below) under conditions of relatively high nitrate <xref ref-type="bibr" rid="pone.0013312-Arrigo1">[34]</xref>, low turbidity and high light. In total, the observed macronutrient dynamics were consistent with previous observations indicating that most nutrients from local watersheds are consumed within the estuary during summer months, rather than exported to the continental shelf <xref ref-type="bibr" rid="pone.0013312-Conomos1">[35]</xref>.</p></sec><sec id="s2c"><title>Microbial abundance and production rates in the CRCM</title><p>We measured autotrophic standing stocks (chlorophyll [chl] <italic>a</italic> concentrations), growth rates of heterotrophic plankton (production rates), and total RNA concentrations (total living microbial biomass) in the CRCM samples. The highest correlations among chl <italic>a</italic> concentration, production rate and total RNA concentration were observed in August, with the highest values for all three biological characteristics observed in the estuary and plume (<italic>R<sup>2</sup></italic> from 0.8 to 0.94), whereas ocean and freshwater end-members had relatively low values (<italic>R<sup>2</sup></italic> from 0.7 to 0.8) (<xref ref-type="fig" rid="pone-0013312-g003">Fig. 3A–C</xref>). Similar, but less pronounced trends for microbial abundance were also observed in June and November in the estuary and plume. In contrast to other sampling times, pre-freshet April samples also showed very high chl <italic>a</italic> and RNA concentrations in the tidal freshwater and at low salinities (0–5 PSU) in the estuary. High RNA and chl <italic>a</italic> concentrations may be at least partly explained by phytoplankton blooms developing in the river freshwater, which is a common occurrence in spring <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>, <xref ref-type="bibr" rid="pone.0013312-Arhonditsis1">[36]</xref>–<xref ref-type="bibr" rid="pone.0013312-Teira1">[38]</xref>. This freshwater phytoplankton is believed to perish in the estuary at high salinities, providing detritus for bacterial community development <xref ref-type="bibr" rid="pone.0013312-Small2">[10]</xref>. However, the pheophytin <italic>a</italic> concentration was low (<xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>), indicating that the majority of chl <italic>a</italic> corresponded to living, rather than detrital, phytoplankton biomass. This observation is consistent with observations from Tillamook Bay indicating phytoplankton biomass accumulation in the middle and lower estuary was especially high in spring and summer <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>. Heterotrophic plankton production rates in the estuary were similar in April and June (0.7 and 0.64 µg CL<sup>−1</sup>h<sup>−1</sup>, respectively), and were consistent with previous measurements indicating that heterotrophic activity in the estuary was higher than that in the adjacent coastal ocean or in the tidal freshwater at that time <xref ref-type="bibr" rid="pone.0013312-Crump2">[14]</xref>, <xref ref-type="bibr" rid="pone.0013312-Crump4">[39]</xref>.</p><fig id="pone-0013312-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g003</object-id><label>Figure 3</label><caption><title>Biological characteristics of the seasonal sample sets collected for microarray analysis.</title><p>In all graphs, the X-axis shows samples sorted as described in <xref ref-type="fig" rid="pone-0013312-g002">Fig. 2</xref>. Salinity plots are superimposed with the bar graphs. (<bold>A</bold>), chlorophyll <italic>a</italic> content; (<bold>B</bold>) bacterial carbon production measured by leucine incorporation; (<bold>C</bold>) total RNA concentrations normalized per liter of sampled water; and (<bold>D</bold>), dissolved organic carbon content. Black vertical lines divide seasons, dashed vertical lines within each season divide the ocean (O) and plume (P) from estuary (E) and river (R) samples. Seasons are shown as months below the graphs.</p></caption><graphic xlink:href="pone.0013312.g003"></graphic></fig><p>DOC concentrations were higher in June (4–4.5 mg/L) than in the other seasons (1.5–2 mg/L, <xref ref-type="fig" rid="pone-0013312-g003">Fig. 3D</xref>). Measurements of particulate organic carbon (POC) performed on a subset of samples from April, June, and November cruises showed that April POC values in the tidal freshwater and estuary (700 to >1000 µg/L) were almost 5 times higher than concentrations observed in the ocean samples collected in April (235 µg/L) and in November river (tidal freshwater and estuary) water samples (281 µg/L). POC concentrations in the tidal freshwater and estuary in June were also relatively high (670 µg/L).</p><p>Given these data, the relatively low microbial abundance (estimated from total RNA concentrations) observed in June compared to April (7 versus 20 µg/L of total RNA; <xref ref-type="fig" rid="pone-0013312-g003">Fig. 3</xref>) was unexpected, since macronutrient concentrations in both the estuary and plume were high (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2C and 2D</xref>). However, the June sample collection occurred during the peak of spring freshet at record-high river discharge (almost 3X higher than other sampling times) and under a considerably reduced water residence time of 2 days (<xref ref-type="fig" rid="pone-0013312-g001">Fig. 1B</xref>). This suggests that microbial populations require longer water residence times to fully develop in the estuary, even when high nutrient concentrations are present. Chl <italic>a</italic> production was likely also depressed in June due to high turbidity and subsequent low light levels that accompanied the high riverflow. Consistent with this hypothesis, the short water residence times during freshet were believed to prevent the development of the abundant estuarine microbiota typically observed during summer and fall within the Tillamook Bay <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>, and in East Coast estuaries (Parker River estuary and Plum Island Sound) <xref ref-type="bibr" rid="pone.0013312-Crump4">[39]</xref>.</p><p>Prior to the onset of seasonal upwelling, the plume tends to have a more freshwater character compared to upwelling-influenced plumes that entrain high-nutrient, high-salinity waters <xref ref-type="bibr" rid="pone.0013312-Bruland1">[4]</xref>. Consistent with this character and with the patterns observed in the river, chl <italic>a</italic> and total RNA were higher in the plume in April compared to June, despite favorable temperatures, summer light conditions, and abundant macronutrient concentrations in the latter month (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2C and D</xref>, <xref ref-type="fig" rid="pone-0013312-g003">Fig. 3D</xref>). We hypothesize that an abundant estuarine microbial community is needed for the development of active plankton in the plume, particularly for the large spring plume, which is characterized by low salinities (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2B</xref>). Plume salinities in June were as low as 10 PSU, in contrast to other seasons, when the plume had salinities of at least 18 to 20 PSU. Alternatively, higher salinities in the plume may also be important for microbial populations to develop, e.g., from an allochthonous marine inoculum.</p></sec><sec id="s2d"><title>Design and validation of functional probes for analysis of gene expression</title><p>Oligonucleotide probes for microarray analysis of gene expression were designed with CombiMatrix probe design software (Probe Weaver) based on: (i) uniqueness within a defined gene and sequence set; (ii) Tm and length within a specified range; 70–75°C, and 35–40-mer, respectively; and (iii) absence of stable secondary structures and repeat sequences. Several published studies describe evaluation of the sensitivity and specificity of various CombiMatrix microarray probe sets <xref ref-type="bibr" rid="pone.0013312-Lodes1">[23]</xref>, <xref ref-type="bibr" rid="pone.0013312-Bolotin1">[25]</xref>, <xref ref-type="bibr" rid="pone.0013312-He1">[26]</xref>, <xref ref-type="bibr" rid="pone.0013312-Lodes2">[40]</xref>–<xref ref-type="bibr" rid="pone.0013312-Maurer1">[42]</xref>. These studies indicated that CombiMatrix probe performance was correlated with theoretical calculations for DNA-DNA hybridizations predicting 1°C reduction in T<sub>m</sub> for every 1% sequence mismatch <xref ref-type="bibr" rid="pone.0013312-Allawi1">[43]</xref>–<xref ref-type="bibr" rid="pone.0013312-Allawi3">[45]</xref>. Theoretically, 3 mismatches in a 35-nt sequence (8.5% sequence mismatch), will result in an 8.5°C drop in T<sub>m</sub>. Empirical validation using the recommended hybridization conditions showed that 3 mismatches reduced T<sub>m</sub> by at least 15°C, resulting in hybridization signals that were at or below background levels (data not shown). Thus, selection of probes for uniqueness within a specified sequence set by the Probe Weaver software requires at least 15°C difference in T<sub>m</sub> to the closest hit/s for nontarget sequences in the set.</p><p>Little sequence information was available for bacterial and archaeal populations in the CRCM from which to design functional probes; the only relevant data corresponded to 16S rDNA clone libraries constructed by Byron Crump and colleagues <xref ref-type="bibr" rid="pone.0013312-Crump1">[13]</xref>, <xref ref-type="bibr" rid="pone.0013312-Crump2">[14]</xref>. Many of these 16S rRNA gene sequences corresponded to uncultivated organisms about which essentially nothing else was known. Thus, our initial strategy was to design microarray probes with somewhat broad specificity (‘common’ probes) that targeted well-annotated functional genes from sequenced genomes of multiple genera within a phylum. We reasoned that such probes would also recognize homologous genes from uncharacterized microorganisms of the same phylum that were likely to be present in environmental water samples. To evaluate this approach, the CombiMatrix probe design algorithm was applied for probe selection from all predicted open reading frames of a subset of 11 fully sequenced genomes from cultured isolates, including three Alphaproteobacteria (<italic>Pelagibacter ubique</italic>, <italic>Erythrobacter litoralis</italic>, <italic>Roseobacter denitrificans</italic>), three Betaproteobacteria (<italic>Polaromonas naphthalenivorans</italic>, <italic>Polynucleobacter</italic> sp, <italic>Rhodoferax ferrireducens</italic>), one Bacteroidetes (<italic>Flavobacterium johnsoniae</italic>), three Gammaproteobacteria (<italic>Marinomonas</italic> sp., <italic>Marinobacter aquaeolei</italic>, <italic>Nitrosococcus oceani</italic>), and one Actinobacteria (<italic>Salinispora arenicola</italic>).</p><p>In total 27,143 non-overlapping 35-mer probes with T<sub>m</sub> between 65°C and 75°C were designed for the 11 genomes, and then each probe was used in a BLAST search through all genomes using CombiMatrix Genotyper software <xref ref-type="bibr" rid="pone.0013312-Lodes1">[23]</xref>, <xref ref-type="bibr" rid="pone.0013312-Lodes3">[41]</xref>. In order to identify common probes targeting homologs from several genomes, Genotyper calculated the number of cross-hybridizing sequence hits with less than 15°C difference in T<sub>m</sub> for each designed probe. Unexpectedly, 80% (21,747) of all probes evaluated produced only one hit: a match of perfect sequence identity from the corresponding genome. Fifteen percent of the probes (4,047) produced 2 hits, the perfect sequence match and one additional partially complementary match from a different genome. Phylogenetic information is limited to the phylum level in Genotyper, thus for the 15% of probes producing more than one hit, two-thirds (2,825) produced a second hit within the same microbial phylum. Finally, 1,349 (5%) of the probes produced 3 or more hits. These probes corresponded to either highly conserved housekeeping genes (i.e. encoding ribosomal proteins), or to genes potentially involved in horizontal gene transfer events (many of these were prophage and transposon-related). Thus, greater than 90% of the microarray probes designed with standard criteria for T<sub>m</sub> and sequence complexity appeared to be specific for their corresponding target gene and a few homologs within the same phylum, even though they were not originally selected for uniqueness within a large multi-genome context. Although the Genotyper analysis was somewhat inconclusive due to the limited annotations produced by the software, it suggested that sequence diversity was too high to allow design of functional probes that would detect homologs across multiple genera within individual families. Therefore, we instead selected all annotated microbial genes for each functional category, and then designed and tested unique probes for each gene.</p><p>Functional genes of interest were selected using the IMG 2.5 system <xref ref-type="bibr" rid="pone.0013312-Markowitz1">[28]</xref>. Approximately 300 sequenced genomes from the Bacteria and Archaea (251 and 52, respectively), representing 246 species, 161 genera and 96 families were selected from taxonomic groups known to be present in the CRCM and/or in water and soil samples from temperate habitats of the Northern hemisphere. Approximately 5000 genes were selected from these genomes, and the corresponding microarray probes were designed for specificity within the set of genes of interest using the standard CombiMatrix Probe Weaver software as described above. In addition, more than 1000 ribosomal RNA sequences from both prokaryotic and eukaryotic microorganisms present in environmental water and soil samples were included into the design process to avoid cross-reactivity with rRNA. The probe sets were synthesized <italic>in situ</italic> on the surface of CombiMatrix CustomArray™ microarrays, and were tested by hybridization with three cDNA pools prepared from 8–10 RNA samples that were collected at different times from disparate locations across the CRCM. The data were normalized using the standard quantile algorithm, and the background signal value was calculated as the mean intensity of the lowest 5% of all signals. Finally, probes that produced significantly high microarray signals – exceeding the background value at least 3-fold in one or more of the pooled target samples – were selected.</p><p>The selected probe set was used to create a custom microarray design of 2226 different oligonucleotide probes based on the CombiMatrix 4X2K CustomArray™ platform.</p><p>Because we are interested in microbial transformations of nitrogen in the CRCM, the majority (1638, 73.6%) of probes corresponded to genes involved in nitrogen metabolism. Probes for genes involved in carbon metabolism (451, 20.3%), housekeeping functions (100, 4.5%) and light perception and utilization (37, 1.6%) were also included. Phylogenetic composition of the probe set is shown in <xref ref-type="table" rid="pone-0013312-t001">Table 1</xref>. Approximately 21% of the probes represented the Archaea (both Crenarchaeota and Euryarchaeota), 60% represented various Proteobacteria, and the remaining probes (19%) represented Bacteroidetes/Chlorobi, Firmicutes, Acidobacteria and “other” (miscellaneous), divisions. Probe coverage was limited to taxonomic groups containing at least one fully sequenced genome.</p><table-wrap id="pone-0013312-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.t001</object-id><label>Table 1</label><caption><title>Phylogenetic composition of the 2226 microarray probe set.</title></caption><alternatives><graphic id="pone-0013312-t001-1" xlink:href="pone.0013312.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Phylum</td><td align="left" rowspan="1" colspan="1">All probes</td><td align="left" rowspan="1" colspan="1">Expressed probes</td><td align="left" rowspan="1" colspan="1">% of expressed</td><td align="left" rowspan="1" colspan="1">April</td><td align="left" rowspan="1" colspan="1">June</td><td align="left" rowspan="1" colspan="1">August</td><td align="left" rowspan="1" colspan="1">November</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><bold>Euryarchaeota</bold></td><td align="left" rowspan="1" colspan="1"><bold>356</bold></td><td align="left" rowspan="1" colspan="1"><bold>296</bold></td><td align="left" rowspan="1" colspan="1"><bold>83.1</bold></td><td align="left" rowspan="1" colspan="1"><bold>0.99</bold></td><td align="left" rowspan="1" colspan="1"><bold>0.97</bold></td><td align="left" rowspan="1" colspan="1"><bold>0.97</bold></td><td align="left" rowspan="1" colspan="1"><bold>0.99</bold></td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>Crenarchaeota</bold></td><td align="left" rowspan="1" colspan="1"><bold>113</bold></td><td align="left" rowspan="1" colspan="1"><bold>100</bold></td><td align="left" rowspan="1" colspan="1"><bold>88.5</bold></td><td align="left" rowspan="1" colspan="1"><bold>1.09</bold></td><td align="left" rowspan="1" colspan="1"><bold>1.15</bold></td><td align="left" rowspan="1" colspan="1"><bold>1.15</bold></td><td align="left" rowspan="1" colspan="1"><bold>1.08</bold></td></tr><tr><td align="left" rowspan="1" colspan="1">Gammaproteobacteria</td><td align="left" rowspan="1" colspan="1">595</td><td align="left" rowspan="1" colspan="1">564</td><td align="left" rowspan="1" colspan="1">94.8</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">1.68</td><td align="left" rowspan="1" colspan="1">1.68</td><td align="left" rowspan="1" colspan="1">1.67</td></tr><tr><td align="left" rowspan="1" colspan="1">Alphaproteobacteria</td><td align="left" rowspan="1" colspan="1">426</td><td align="left" rowspan="1" colspan="1">418</td><td align="left" rowspan="1" colspan="1">98.1</td><td align="left" rowspan="1" colspan="1">2.13</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">2.12</td><td align="left" rowspan="1" colspan="1">2.09</td></tr><tr><td align="left" rowspan="1" colspan="1">Betaproteobacteria</td><td align="left" rowspan="1" colspan="1">187</td><td align="left" rowspan="1" colspan="1">184</td><td align="left" rowspan="1" colspan="1">98.4</td><td align="left" rowspan="1" colspan="1">1.6</td><td align="left" rowspan="1" colspan="1">1.59</td><td align="left" rowspan="1" colspan="1">1.62</td><td align="left" rowspan="1" colspan="1">1.61</td></tr><tr><td align="left" rowspan="1" colspan="1">Deltaproteobacteria</td><td align="left" rowspan="1" colspan="1">147</td><td align="left" rowspan="1" colspan="1">145</td><td align="left" rowspan="1" colspan="1">98.6</td><td align="left" rowspan="1" colspan="1">2.44</td><td align="left" rowspan="1" colspan="1">2.45</td><td align="left" rowspan="1" colspan="1">2.29</td><td align="left" rowspan="1" colspan="1">2.32</td></tr><tr><td align="left" rowspan="1" colspan="1">Chlorobi</td><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">1.99</td><td align="left" rowspan="1" colspan="1">2.05</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">1.94</td></tr><tr><td align="left" rowspan="1" colspan="1">Firmicutes</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">94.7</td><td align="left" rowspan="1" colspan="1">1.59</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">1.67</td><td align="left" rowspan="1" colspan="1">1.75</td></tr><tr><td align="left" rowspan="1" colspan="1">Actinobacteria</td><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">92.2</td><td align="left" rowspan="1" colspan="1">2.41</td><td align="left" rowspan="1" colspan="1">1.95</td><td align="left" rowspan="1" colspan="1">2.09</td><td align="left" rowspan="1" colspan="1">2.12</td></tr><tr><td align="left" rowspan="1" colspan="1">Bacteroidetes</td><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">39</td><td align="left" rowspan="1" colspan="1">83.0</td><td align="left" rowspan="1" colspan="1">1.27</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">1.14</td><td align="left" rowspan="1" colspan="1">1.29</td></tr><tr><td align="left" rowspan="1" colspan="1">Planctomycetes</td><td align="left" rowspan="1" colspan="1">37</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">97.3</td><td align="left" rowspan="1" colspan="1">1.73</td><td align="left" rowspan="1" colspan="1">1.93</td><td align="left" rowspan="1" colspan="1">1.67</td><td align="left" rowspan="1" colspan="1">1.67</td></tr><tr><td align="left" rowspan="1" colspan="1">Chloroflexi</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">96.4</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">2.57</td><td align="left" rowspan="1" colspan="1">2.5</td><td align="left" rowspan="1" colspan="1">2.57</td></tr><tr><td align="left" rowspan="1" colspan="1">Epsilonproteobacteria</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">88.9</td><td align="left" rowspan="1" colspan="1">1.42</td><td align="left" rowspan="1" colspan="1">1.33</td><td align="left" rowspan="1" colspan="1">1.31</td><td align="left" rowspan="1" colspan="1">1.47</td></tr><tr><td align="left" rowspan="1" colspan="1">Magnetococci</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">3.91</td><td align="left" rowspan="1" colspan="1">3.01</td><td align="left" rowspan="1" colspan="1">3.52</td></tr><tr><td align="left" rowspan="1" colspan="1">Thermotogae</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">1.44</td><td align="left" rowspan="1" colspan="1">1.86</td><td align="left" rowspan="1" colspan="1">1.87</td><td align="left" rowspan="1" colspan="1">1.57</td></tr><tr><td align="left" rowspan="1" colspan="1">Verrucomicrobia</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">2.92</td><td align="left" rowspan="1" colspan="1">2.97</td><td align="left" rowspan="1" colspan="1">4.64</td><td align="left" rowspan="1" colspan="1">2.68</td></tr><tr><td align="left" rowspan="1" colspan="1">Zetaproteobacteria</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">1.95</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">1.81</td><td align="left" rowspan="1" colspan="1">2.09</td></tr><tr><td align="left" rowspan="1" colspan="1">Aquificae</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">1.84</td><td align="left" rowspan="1" colspan="1">1.58</td><td align="left" rowspan="1" colspan="1">2.17</td></tr><tr><td align="left" rowspan="1" colspan="1">Acidobacteria</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">2.32</td><td align="left" rowspan="1" colspan="1">3.07</td><td align="left" rowspan="1" colspan="1">2.62</td><td align="left" rowspan="1" colspan="1">2.67</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>The probes were selected as expressed if corresponding signal intensities exceeded the baseline value (background plus 3X standard deviations) in at least 2 samples. The median expression levels were calculated for probes corresponding to expressed genes for each phylogenetic group, and then represented as fold changes over the median signal intensity of the whole microarray data set. Median expression levels for archaeal divisions are shown in bold.</p></fn></table-wrap-foot></table-wrap></sec><sec id="s2e"><title>Bioinformatic analysis of microarray probe specificity using large microbial and environmental sequence databases</title><p>Further assessment of probe specificity was done using large microbial metagenome sequence databases that recently became available for batch BLAST searches through CAMERA, the Community Cyberinfrastructure for Advanced Marine Microbial Ecology Research and Analysis (<ext-link ext-link-type="uri" xlink:href="http://camera.calit2.net/index.php">http://camera.calit2.net/index.php</ext-link> <xref ref-type="bibr" rid="pone.0013312-Seshadri1">[46]</xref>). CAMERA currently serves as a repository for approximately 80 microbial sequencing databases, including (i) metagenome data for oceans, lakes, rivers, hot springs and soils (i.e. <xref ref-type="bibr" rid="pone.0013312-DeLong1">[47]</xref>, <xref ref-type="bibr" rid="pone.0013312-Rusch1">[48]</xref>); (ii) human and animal microbiome sequencing data; and (iii) fully and partially sequenced microbial genomes deposited at the National Center for Biotechnology Information (NCBI). CAMERA's BLAST tool was used for specificity analyses of all 2226 probes against 9 environmental databases of assembled sequences, as well as “CAMERA's Non-Identical Nucleotide Sequences” database, the largest non-redundant database with 65 million sequences totaling approximately 207 gigabases. We applied the CAMERA BLAST search tool with subsequent selection of resulting hits for potential cross-hybridization with microarray probes, using the following criteria: (i) an alignment length of at least 30 nucleotides with <italic>E</italic> value of less than 0.05%, and (ii) at least 90% identity score for the query-match alignment. These criteria resulted in ≥27 perfectly matched positions for each hit, allowing up to 8 mismatches (≥77% overall sequence identity). Because we wished to evaluate the maximum number of potential cross-hybridizing hits, we intentionally chose a less stringent approach than that used by CombiMatrix <xref ref-type="bibr" rid="pone.0013312-Lodes3">[41]</xref> or in experimental work that determined the threshold value for probe specificity as 87% sequence identity <xref ref-type="bibr" rid="pone.0013312-TaroncherOldenburg1">[21]</xref>.</p><p>The CAMERA databases used for this bioinformatic evaluation and the numbers of corresponding hits are shown in <xref ref-type="table" rid="pone-0013312-t002">Table 2</xref>. Only four databases (shown with asterisks in <xref ref-type="table" rid="pone-0013312-t002">Table 2</xref>) provided annotation for the results, which was taxonomic, rather than functional, in nature. As expected, analysis with the environmental virus and eukaryotic microbial sequence databases produced either zero or very small numbers of cross-hybridizing hits (fewer than 10). Only 3 of our probes produced cross-hybridizing hits from eukaryotic genomes, despite the large amount (over 6 gigabases) of non-redundant eukaryotic coding sequences evaluated. Very small numbers of hits were also observed for the databases from Minnesota farm soil, or from open-ocean and deep-water samples from the geographically distant Hawaii ocean time-series. Only 118 probes from our set of 2226 produced hits in the 7 Gb NCBI environmental samples database, reflecting the fact that the Pacific Northwest coastal margin is poorly represented in metagenomic sequencing projects. In contrast, a relatively small 0.8 Gb Moore Foundation Marine Microbial Genomes database produced hits for 526 probes (24% of the total). Thus, the numbers of hits were apparently independent of the database size, and instead were determined by taxonomic, geographic and habitat relevance. Furthermore, almost all probes, 93 to 98% (2,077 and 2,186, respectively), produced hits in the large integrated databases of assembled prokaryotic sequences (bolded in <xref ref-type="table" rid="pone-0013312-t002">Table 2</xref>).</p><table-wrap id="pone-0013312-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.t002</object-id><label>Table 2</label><caption><title>Cross-hybridizing sequence hits resulting from bioinformatic evaluation of 2226 microarray probes against the CAMERA repository.</title></caption><alternatives><graphic id="pone-0013312-t002-2" xlink:href="pone.0013312.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Database name as defined in CAMERA</td><td align="left" rowspan="1" colspan="1">Total length (bp)</td><td align="left" rowspan="1" colspan="1"># sequences</td><td align="left" rowspan="1" colspan="1"># Cross-hyb hits</td><td align="left" rowspan="1" colspan="1"># probes with hits</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GOS: Site-specific 16S Sequences (N)</td><td align="left" rowspan="1" colspan="1">3,118,182</td><td align="left" rowspan="1" colspan="1">4,125</td><td align="left" rowspan="1" colspan="1">n/a</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">GOS: move858 Assembled 0.002-0.22 Chesapeake Bay (N)</td><td align="left" rowspan="1" colspan="1">8,669,804</td><td align="left" rowspan="1" colspan="1">5,357</td><td align="left" rowspan="1" colspan="1">n/a</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">FarmSoil: Assembled Sequences (N) (Minnesota farm soil)</td><td align="left" rowspan="1" colspan="1">144,897,582</td><td align="left" rowspan="1" colspan="1">139,340</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">HOT: All ORFs (N) (Hawaii Ocean Time-series ALOHA)</td><td align="left" rowspan="1" colspan="1">169,784,453</td><td align="left" rowspan="1" colspan="1">449,086</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">Moore Foundation Marine Microbial Genomes (N)*</td><td align="left" rowspan="1" colspan="1">856,811,427</td><td align="left" rowspan="1" colspan="1">12,886</td><td align="left" rowspan="1" colspan="1">630</td><td align="left" rowspan="1" colspan="1">526</td></tr><tr><td align="left" rowspan="1" colspan="1">GOS: Combined Assembly Coding Sequences (N)</td><td align="left" rowspan="1" colspan="1">3,668,987,939</td><td align="left" rowspan="1" colspan="1">6,115,750</td><td align="left" rowspan="1" colspan="1">137</td><td align="left" rowspan="1" colspan="1">62</td></tr><tr><td align="left" rowspan="1" colspan="1">Eukaryotic Microbial Genomes (N)*</td><td align="left" rowspan="1" colspan="1">6,342,658,807</td><td align="left" rowspan="1" colspan="1">1,453,409</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">3</td></tr><tr><td align="left" rowspan="1" colspan="1">All NCBI Environmental Samples (ENV_NT)</td><td align="left" rowspan="1" colspan="1">7,194,061,284</td><td align="left" rowspan="1" colspan="1">17,695,887</td><td align="left" rowspan="1" colspan="1">218</td><td align="left" rowspan="1" colspan="1">118</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>All Prokaryotic Genomes (N)*</bold></td><td align="left" rowspan="1" colspan="1">9,577,197,991</td><td align="left" rowspan="1" colspan="1">655,666</td><td align="left" rowspan="1" colspan="1">3770</td><td align="left" rowspan="1" colspan="1">2,186</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold>CAMERA's Non-Identical Nucleotide Sequences (N)*</bold></td><td align="left" rowspan="1" colspan="1">179,511,589,666</td><td align="left" rowspan="1" colspan="1">38,512,986</td><td align="left" rowspan="1" colspan="1">4,367</td><td align="left" rowspan="1" colspan="1">2,077</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Columns 2 and 3 show the total amount of sequence information and the number of individual sequences, respectively, within each database in Column 1. Column 4 shows the numbers of sequences selected as potentially cross-hybridizing with the microarray probes. Column 5 shows the numbers of microarray probes that produced at least one hit in the corresponding CAMERA databases. Asterisks show databases with taxonomic annotations. The two databases used for taxonomic analysis of probe hits are shown in bold.</p></fn></table-wrap-foot></table-wrap><p>The distribution of microarray probe hits within CAMERA's ‘All Prokaryotic Genomes” and ‘Non-Identical Nucleotide Sequence’ databases, the two largest databases from <xref ref-type="table" rid="pone-0013312-t002">Table 2</xref>, were analyzed to determine taxonomic rank, and the results are shown in the <xref ref-type="table" rid="pone-0013312-t003">Table 3</xref>. Only a small number of probes (1.8 to 7%, <xref ref-type="table" rid="pone-0013312-t003">Table 3</xref>) did not produce any significant hits. The majority of probes (86 to 91%) produced BLAST hits from the same genus that was used to design the probe (and typically, also from the same species, although this was not always determinable due to incomplete annotations). When matches to multiple genera within a single microbial family were allowed, 91% to 97% of the probes placed into this category (<xref ref-type="table" rid="pone-0013312-t003">Table 3</xref>). Approximately 1.5–2% (34 to 46, <xref ref-type="table" rid="pone-0013312-t003">Table 3</xref>) of probes hit genes from organisms outside of the microbial family containing the target organism. Among them fewer than 1% hit genes outside of the target organism's phylum (8 to 17 probes, for ‘All Prokaryotic’ and ‘CAMERA's Non-Identical’ databases, respectively). All but one of the probes in this latter group produced multiple hits, which typically corresponded to genes from the same genus/family as the target organism, and 1–2 hits from outside of the family. Taken together, our analyses indicated that the microarray probe set was largely specific for the target gene and organism, but approximately 10 to 15% of the probes also detected homologs from different genera within the corresponding target family, or, in very rare cases, outside of the targeted family. These results, therefore, suggested that the probe set might be useful not only for detecting gene expression from targeted organisms, but also from related, but currently uncharacterized microorganisms in the CRCM.</p><table-wrap id="pone-0013312-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.t003</object-id><label>Table 3</label><caption><title>Distribution of microarray probe hits within CAMERA's Non-Identical Nucleotide Sequence and All Prokaryotic Genomes databases.</title></caption><alternatives><graphic id="pone-0013312-t003-3" xlink:href="pone.0013312.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Database</td><td align="left" rowspan="1" colspan="1">CAMERA's Non-Identical</td><td align="left" rowspan="1" colspan="1">All Prokaryotic</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Total number of probes with hits</td><td align="left" rowspan="1" colspan="1">2,077 (93.3%)</td><td align="left" rowspan="1" colspan="1">2,186 (98.2%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with single hits</td><td align="left" rowspan="1" colspan="1">1,459 (65.5%)</td><td align="left" rowspan="1" colspan="1">1,595 (71.6%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with double hits</td><td align="left" rowspan="1" colspan="1">330 (14.8%)</td><td align="left" rowspan="1" colspan="1">339 (15.2%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with multiple hits</td><td align="left" rowspan="1" colspan="1">288 (12.9%)</td><td align="left" rowspan="1" colspan="1">252 (11.3%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with all hits from the same genus</td><td align="left" rowspan="1" colspan="1">1, 921 (86.3%)</td><td align="left" rowspan="1" colspan="1">2,019 (90.7%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with all hits from the same family</td><td align="left" rowspan="1" colspan="1">2,031 (91.4%)</td><td align="left" rowspan="1" colspan="1">2,152 (96.7%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Probes with some hits outside of the family</td><td align="left" rowspan="1" colspan="1">46 (2%)</td><td align="left" rowspan="1" colspan="1">34 (1.5%)</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt103"><label></label><p>The numbers in brackets are percentages of the total number of probes (2,226).</p></fn></table-wrap-foot></table-wrap></sec><sec id="s2f"><title>Experimental validation of microarray probe specificity</title><p>Two approaches were used for experimental validation of microarray probe specificity under the microarray hybridization conditions used in this study. The first was aimed at determining whether the probes cross-hybridized with rRNA-derived targets, which is considered to be a serious problem in microarray applications to environmental samples <xref ref-type="bibr" rid="pone.0013312-He1">[26]</xref>. For this reason, the probes were designed against cross-hybridization with over 1000 rRNA sequences (both prokaryotic and eukaryotic). In addition, application of the CAMERA Blast tool to evaluate probe matches to the ‘GOS Site-specific 16S Sequences’ database containing over 4,000 assembled sequences (over 3 Mb in total, <xref ref-type="table" rid="pone-0013312-t002">Table 2</xref>), did not produce any significant hits. To evaluate cross-hybridization of rRNA in the target preparations, we performed subtractive hybridization assays. For the majority of probes, the rRNA-depleted samples produced higher signals in comparison with untreated samples (data not shown). Thus, an increased proportion of mRNA in the target preparations resulted in an increase in signal intensity, indicating that our probes hybridized specifically with corresponding mRNA-derived targets (and did not cross-hybridize with rRNA).</p><p>Additional validation was performed using laboratory-grown cultures of <italic>Bacillus subtilis</italic> and <italic>Pseudomonas putida</italic>. Total RNA from each culture was isolated, converted to fluorescently labeled cDNA, and hybridized in duplicate to microarrays. A typical example of the raw signal intensity data sorted by organism annotation is shown for microarray hybridization with <italic>B. subtilis</italic> targets (<xref ref-type="supplementary-material" rid="pone.0013312.s001">Fig. S1</xref>). High signal intensities exceeding background at least 10X were observed for 9 out of 10 of the probes designed from <italic>B. subtilis</italic> functional genes. In contrast, only 3 probes from closely related <italic>Bacillus</italic> species (∼15% of probes in this category) produced similar signals with <italic>B. subtilis</italic> cDNA. For probes corresponding to unrelated species from other families, 0.2% produced signals exceeding 3X background, and 9% produced signals exceeding 1X background. Thus, 99.8% of unrelated probes had signal values that did not exceed 3X background (calculated as the lowest 5% of all signals), which was our cutoff for positive microarray signals. Similar data were obtained for <italic>P. putida</italic> laboratory cultures (not shown). Experimental validation of probe specificity closely mirrored results from the bioinformatics analysis, indicating that the majority of probes were specific for the intended target, while approximately 15% of probes to genes from closely related taxa, and only a very minor proportion (0.2%) of the probes to genes from organisms outside of the targeted taxonomic family, hybridized to <italic>B. subtilis</italic> cDNA.</p></sec><sec id="s2g"><title>Microarray data acquisition from environmental samples and analysis of reproducibility</title><p>Total RNA from 64 samples (14 from August, 2007; 17 from November, 2007; 17 from April, 2008; and 16 from June, 2008) was used to prepare microarray hybridization targets by reverse transcription into fluorescently labeled cDNA. Due to relatively high microbial abundance in the CRCM, the RNA yields were sufficient to generate labeled targets without additional RNA amplification steps, thus preventing potential bias introduced by this common technique <xref ref-type="bibr" rid="pone.0013312-Parro1">[22]</xref>. The targets were hybridized in duplicate, and replicate hybridizations were compared using scatter plots to estimate experimental variation (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2</xref>). The replicate data were highly consistent when the same target was hybridized to two different sectors of the same microarray chip, to different microarray chips, or to the same microarray chip upon re-use (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2A</xref>). The data points grouped around the 45° line and were well within the standard 2-fold cut-off lines (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2A</xref>). The corresponding correlation coefficients (<italic>R<sup>2</sup></italic>) were >0.98, demonstrating that inter- and intra-chip variation of the microarray platform was quite low. A similar approach was used to evaluate variability introduced by sampling and target preparation protocols (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2B</xref>). We observed high reproducibility (<italic>R<sup>2</sup></italic>≥0.96), for target preparations from the same RNA in two independent reverse transcription reactions, or from two different freshwater samples collected at the same location at approximately the same time (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2B</xref>). However, increase in the variation (<italic>R<sup>2</sup></italic> = 0.94) was observed if two water samples were collected at different river locations 50 km apart on the same day (<xref ref-type="supplementary-material" rid="pone.0013312.s002">Fig. S2B</xref>, right). Taken together, these data indicate that experimental variability was determined largely by sampling location, and not by sample handling, RNA isolation, or microarray hybridizations. Thus, duplicate hybridizations performed for the same sample were averaged to create mean expression values. The data were deposited according to the MIAME reporting guidelines and are accessible from the Gene Expression Omnibus (GEO) public database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/projects/geo/">http://www.ncbi.nlm.nih.gov/projects/geo/</ext-link>) with the accession number GSE18303.</p></sec><sec id="s2h"><title>Selection and analysis of significantly expressed genes</title><p>Genes were selected as significantly expressed if the corresponding probes produced normalized signal intensities exceeding the background value plus 3X standard deviation in at least 2 different samples. Among 2226 genes represented on the microarrays, 2076 (93%) fit this criterion, and were selected to compare gene expression levels of the major prokaryotic phyla represented on the microarrays (<xref ref-type="table" rid="pone-0013312-t001">Table 1</xref>). The median expression value was calculated for the whole dataset and for each phylum separately in four sample groups, representing different sampling times. Phylum-specific gene expression levels are shown in <xref ref-type="table" rid="pone-0013312-t001">Table 1</xref> as fold changes relative to the overall median expression value for all normalized microarrays. For bacteria, these levels varied from 1.3 to 4.6 fold (<xref ref-type="table" rid="pone-0013312-t001">Table 1</xref>), with the highest values (approximately 2.5 to 4.5-fold across all seasons) observed for genes from the Verrucomicrobia and unclassified Proteobacteria (in particular, Magnetococcus). Relatively high levels of expression (>2-fold across all seasons) were also observed for genes of Chloroflexi, Acidobacteria, Alphaproteobacteria and Deltaproteobacteria. In contrast, expression levels for both archaeal divisions were close to the median and much lower than for many of the bacterial groups (<xref ref-type="table" rid="pone-0013312-t001">Table 1</xref>). As is true in general for microarray approaches, observed changes in the signal for a particular transcript may be influenced by abundance of the organism expressing the transcript, by differential regulation of the corresponding gene, or both. In any case, among 246 prokaryotic species represented by 2226 microarray probes, 244 had at least one significantly expressed gene in multiple samples, validating our initial selection of genomes.</p><p>For analysis of gene expression patterns in the 64-sample set, the signal value for each probe in each sample was calculated as the log-transformed ratio of normalized intensity versus the background (common for all microarray data after normalization). Two-dimensional (2D) clustering of genes (X-axis) and samples (Y-axis) is shown in <xref ref-type="fig" rid="pone-0013312-g004">Fig. 4</xref>. The heat map clearly indicates that the majority of genes were significantly expressed in multiple samples, and changes were observed in the corresponding gene expression profiles with respect to different seasons (described below). These data suggest that similar microbial taxa were metabolically active in different habitats across the CRCM, and were consistent with the fact that the probes most likely detected similar gene expression patterns from different species within a microbial family across sampling sites.</p><fig id="pone-0013312-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g004</object-id><label>Figure 4</label><caption><title>Expression profiles in the seasonal sample sets, displayed as fold changes over the background.</title><p>Two-dimensional clustering of 2076 probes corresponding to expressed genes (columns) and 56 samples (rows) was done using BRB ArrayTools software based on similarity matrix calculated with an agglomerative algorithm, and complete link correlation. Expression values are colored according to the ratio of signal to the background (displayed using a logarithmic scale). White to black scale indicates magnitude (from low to high). Sample names are composed of season (Apr, Jun, Aug, and Nov), habitat (R, river and estuary; O, plume and coastal ocean), depth (S, surface 1–3 m; M, mid-depth 4–50 m), location, and unique sample number. Location codes: number shows distance from the coast in km; CR, Columbia River transect in the plume and coastal ocean; NH, Newport Hydroline transect in the coastal ocean at Newport; AST and HAM, estuary locations near Astoria (river mile 7–9) and Hammond (river mile 5), respectively; TID, estuary locations in the tidal basin (river mile 22–23); BA, river location at Beaver Army Dock (river mile 53) near Quincy, Oregon; UP, river at mile 74.</p></caption><graphic xlink:href="pone.0013312.g004"></graphic></fig><p>Clustering of samples by similarity of gene expression patterns is shown as a dendrogram on the right side of the heat map, which divides the samples into two major clusters (<xref ref-type="fig" rid="pone-0013312-g004">Fig. 4</xref>, upper and lower clusters are divided by a white line). June and November samples were present in both clusters. April and August samples, however, placed into lower and upper clusters, respectively, clearly separating gene expression profiles in April and August from one another. This pronounced seasonal difference in gene expression patterns evidently exceeded the differences in expression across habitats with different salinities (tidal freshwater, estuary, plume and adjacent coastal ocean). Consistent with this observation, calculated pair-wise correlation coefficients (<italic>R</italic><sup>2</sup>) for expression patterns in samples with variable salinities collected from the estuary and plume and analyzed for each season were 0.8 to 0.9 (<xref ref-type="supplementary-material" rid="pone.0013312.s003">Fig. S3</xref>). Correlation to freshwater end-member expression patterns decreased with increasing salinity, and expression patterns differed most in ocean samples at high salinities (<italic>R<sup>2</sup></italic> from 0.6 to 0.7). This indicates that estuarine microorganisms may carry out similar ecological functions over a range of salinities, and supports the idea put forth by Crump <italic>et al.</italic> <xref ref-type="bibr" rid="pone.0013312-Crump2">[14]</xref> that they form communities that are distinct from those in tidal freshwater and the adjacent coastal ocean. Also consistent with the results presented here, a recent report indicated that seasonal differences in the composition of Chesapeake Bay bacterioplankton assemblages were more significant than the variation accounted for by gradients in salinity <xref ref-type="bibr" rid="pone.0013312-Kan1">[49]</xref>.</p></sec><sec id="s2i"><title>Variation in differential gene expression patterns across salinity gradients in the CRCM with season and habitat</title><p>We evaluated season- and habitat-specific differential gene expression in archaeal and bacterial populations in the estuary, plume and adjacent coastal ocean compared to the tidal freshwater baseline. The baselines were calculated for each sampling season by averaging of the two freshwater samples collected at the same river location (0 PSU). For samples collected within a season, gene expression ratios were calculated relative to the corresponding freshwater baseline. Overall, 1496 (72% of the total) genes exhibited greater than 2-fold ratio changes in at least two samples, and these were selected for 2D clustering analysis (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>). The clustering diagram was dominated by five distinct patterns, designated “clusters A–E”; each consisting of hundreds of expressed genes (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>, divided by the white lines). The differential gene expression patterns could not be accounted for by one unique environmental factor, or by linear combinations of factors. However, use of Primer 6 software to generate multidimensional scaling plots of maximum variability among samples revealed a seasonal trend in the data set (<xref ref-type="supplementary-material" rid="pone.0013312.s004">Fig. S4</xref>), similar to results shown in <xref ref-type="fig" rid="pone-0013312-g004">Fig. 4</xref>, with April and August samples clustering separately. Also consistent with the previous analysis, MDS plots did not reveal distinct clustering of samples based on habitat (tidal freshwater, estuary, plume or adjacent ocean). Close examination of the 2D clustering diagram in <xref ref-type="fig" rid="pone-0013312-g005">Figure 5</xref>, however, suggested additional trends in the dataset. Cluster A was composed of samples from low-light habitats in the estuary, mainly from November. This cluster also included estuarine samples that were collected from below the surface in June (4) and August (1), and one sample collected at the surface in June when turbidity was extremely high (data not shown). The relatively low differential gene expression observed in this cluster compared to the freshwater baseline samples can be explained by the fact that the river also experienced low light levels in November, as well as high turbidity in June. Photosynthetically active radiation (PAR) data were collected only for June and April, but the values corresponding to samples collected at 2-m depth in the estuary were almost 3X higher in April than in June (2.7×10<sup>13</sup> versus 9.5×10<sup>12</sup> micromole/m<sup>2</sup>/s, respectively). Consistent with the idea that light levels may have influenced sample clustering, analysis of the small set of probes on our microarrays corresponding to genes involved in light perception and utilization indicated that many of them produced relatively low signals for the samples within cluster A, and for June and November baseline samples, while the corresponding microarray probe signals were 2–10X higher for April samples in clusters D and E.</p><fig id="pone-0013312-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g005</object-id><label>Figure 5</label><caption><title>Differential expression across salinity gradients of the CRCM compared to seasonal tidal freshwater baselines.</title><p>Gene expression ratios were calculated for each probe in each sample as log (base 2) transformed ratios over the corresponding seasonal average freshwater values. Genes were selected as differentially expressed if they had ratios over 2-fold in at least two different samples. Two-dimensional clustering of 1496 differentially expressed genes (columns) and 56 samples (rows) was performed using BRB ArrayTools software based on similarity matrix calculated with an agglomerative algorithm, and complete link correlation. Expression data are colored according to fold changes in the mRNA levels (displayed using a logarithmic scale). Green and red colors indicate decreased and increased values, respectively, within a sample relative to the seasonal freshwater baseline. The color scale indicates the magnitude of change. Sample names are colored according to the season: blue, April 2008; green, June 2008; red, August 2007; and black, November 2007. Samples are designated as described in <xref ref-type="fig" rid="pone-0013312-g004">Fig. 4</xref>.</p></caption><graphic xlink:href="pone.0013312.g005"></graphic></fig><p>Cluster B was composed of samples collected entirely in August across habitats characterized by variable salinities. Microbial biomass was dominated by phytoplankton throughout the CRCM in August (<xref ref-type="fig" rid="pone-0013312-g003">Fig. 3A and C</xref>), which could account for the similarities in gene expression patterns compared to the freshwater baseline. In fact, the two August samples that did not fall into this cluster contained the lowest concentrations of chl <italic>a</italic> (<xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>). Also consistent with this idea, a large bloom of the ciliated protist <italic>Myrionecta rubra</italic>, containing cryptophyte chloroplasts, was observed in the estuary during the August cruise in 2007. Samples in cluster B were also characterized by higher temperatures compared to those in other groups, but because this was also the case for the freshwater baseline samples, increased temperature alone would not have resulted in the differential gene expression patterns that were observed.</p><p>Clusters C–E contained samples exhibiting the highest levels of differential gene expression. Cluster C was composed of samples collected at high salinities in the plume and adjacent coastal ocean in June, at locations where chl <italic>a</italic> and total microbial abundance (total RNA) were very low (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2</xref>) and DOC was high (<xref ref-type="fig" rid="pone-0013312-g003">Fig. 3</xref>). This was in contrast to the June estuary samples, which did not fall into this cluster, and the freshwater baseline samples, both of which had higher concentrations of chl <italic>a</italic> and lower DOC. One August and two November ocean samples included in cluster C also had low total RNA and chl <italic>a</italic> concentrations (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2</xref>), and intermediate DOC levels (<xref ref-type="fig" rid="pone-0013312-g003">Fig. 3</xref>). Thus, both salinity and availability of DOC may have influenced the grouping of samples into this cluster.</p><p>Clusters D and E contained samples from seasons when nitrate levels were high in the tidal freshwater baseline. The samples had nitrate concentrations that were, on average, 25% (cluster D) and 47% (cluster E) of the concentrations in corresponding seasonal freshwater baseline samples (<xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>), suggesting nitrate uptake by active microorganisms. Nitrate concentrations in cluster C samples were also about 70%, on average, of the levels in the corresponding freshwater baseline samples. This interpretation is consistent with the relatively high levels of differential gene expression observed for samples in these three clusters. For cluster D (April) samples, the nitrate:phosphate ratios were additionally 4–20X greater in the river baseline samples compared to the estuary and plume (240 versus 59 and 11, respectively, <xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref> and <xref ref-type="fig" rid="pone-0013312-g003">Fig. 3</xref>), which may correspond to the grouping of samples into this cluster.</p><p>Cluster E was composed of samples collected near the mouth of the Columbia River (with one exception) in November, June, and April. Close proximity of cluster E samples to the mouth of the river, and the failure to determine trends from other measured variables suggests that gene expression patterns may be correlated with factor(s) that we did not measure. All together, the trends described above suggest that bacterial and archaeal gene expression in the CRCM may be influenced (either directly or indirectly) by light (clusters A, D, E), phytoplankton biomass (cluster B), salinity (cluster C), proximity to the river mouth (cluster E), and availability of DOC (cluster C), nitrate, and phosphate (clusters C, D, E).</p></sec><sec id="s2j"><title>Stability in gene expression patterns over time and across salinity gradients</title><p>Within each cluster, the highest similarity in gene expression patterns was observed for samples collected close together in space (in the estuary and plume) and time (1–2 day intervals). This was the case even though water salinity in these samples varied widely, e.g. from 10 to 30 PSU in April (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>, cluster D) and 6 to 25 PSU in August (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>, cluster B). Some of these samples had nearly identical gene expression patterns, despite being collected at different stages of the tidal cycle. With the exception of salinity, the environmental data for samples with such highly similar gene expression patterns were also very similar, with calculated pair-wise correlation coefficients in the range of 0.8–0.85 (for all numeric environmental data in <xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>). In November and August, some samples collected 5 to 10 days apart were also highly similar in terms of gene expression patterns, and, again, the temporal stability of these patterns corresponded to relatively similar environmental conditions (<italic>R<sup>2</sup></italic> from 0.7 to 0.75). This suggests that gene expression, and, potentially, metabolic activity of microbial assemblages in the estuary and plume are stable over time and across a wide range of salinities, as long as other physical and chemical factors remain relatively unchanged.</p></sec><sec id="s2k"><title>Functional gene expression signatures associated with salinity gradients of the CRCM</title><p>The largest and most uniform common gene expression pattern encompassed all of the April samples and many November and June samples, and it is represented by clusters C, D, and E (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>). Self-organizing maps (SOM) and Analysis of Variation (ANOVA at <italic>P</italic> values below 0.001) were used to select probes pertaining to this prominent expression pattern in the corresponding samples. The results of the two analyses were nearly identical, and included approximately 500 probes with higher relative signals compared to the freshwater baseline.</p><p>Several functional gene groups were represented by the probe subset, including genes involved in nitrate assimilation (<xref ref-type="fig" rid="pone-0013312-g006">Fig. 6</xref>), dissimilatory nitrate/nitrite reduction (e.g., in denitrification and dissimilatory nitrate reduction to ammonia pathways, <xref ref-type="fig" rid="pone-0013312-g007">Fig. 7</xref>), and carbon utilization pathways (<xref ref-type="fig" rid="pone-0013312-g008">Fig. 8</xref>). The corresponding 143 probes used to generate the data for <xref ref-type="fig" rid="pone-0013312-g006">Fig. 6</xref>, <xref ref-type="fig" rid="pone-0013312-g007">7</xref>, <xref ref-type="fig" rid="pone-0013312-g008">8</xref> were additionally evaluated by BLAST searches against the entire NCBI nucleotide sequence database using the criteria described above. The majority of probes (138 out of 143) produced hits only to the intended functional gene target within the same microbial genus. Five probes out of 143 produced an additional hit from a different genus within the same family. These results were consistent with the bioinformatic and experimental probe validations described above, and indicated that the selected probe subset hybridized with specific functional targets >96% of the time.</p><fig id="pone-0013312-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g006</object-id><label>Figure 6</label><caption><title>High relative expression of nitrate assimilation genes in comparison to the corresponding seasonal freshwater baseline.</title><p>Gene sets were selected from clusters C–E of <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>. Expression ratio calculations, 2D clustering, and the heat map display were prepared as described in <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>, except that samples are shown as columns, and genes as rows. Gene names include the corresponding organism and phylum information. The sample cluster that has the upregulation pattern is separated from the rest of heat map by the while line, and the corresponding sample names are shown in bold.</p></caption><graphic xlink:href="pone.0013312.g006"></graphic></fig><fig id="pone-0013312-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g007</object-id><label>Figure 7</label><caption><title>High relative expression of denitrification pathway genes in comparison to the corresponding seasonal freshwater baseline.</title><p>Gene sets were selected from clusters C–E of <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>. Expression ratio calculations, 2D clustering, and heat map display were done as described in <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref> and <xref ref-type="fig" rid="pone-0013312-g006">6</xref>.</p></caption><graphic xlink:href="pone.0013312.g007"></graphic></fig><fig id="pone-0013312-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g008</object-id><label>Figure 8</label><caption><title>High relative expression of genes involved in carbon fixation, scavenging, and storage.</title><p>Gene sets were selected from clusters C–E of <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>. Expression ratio calculations were performed relative to the corresponding seasonal freshwater baseline; 2D clustering, and heat map display were done as described in <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref> and <xref ref-type="fig" rid="pone-0013312-g006">6</xref>.</p></caption><graphic xlink:href="pone.0013312.g008"></graphic></fig><p>The complete probe set used for microarray analysis contained all annotated genes for nitrogen metabolism for each included genome, and in many cases, genes from both nitrate assimilation and dissimilatory nitrate/nitrite reduction pathways were present in the same organism. However, the phylogenetic composition of the two associated probe subsets showed very little overlap (less than 15%), indicating that the enhanced signals for genes in these pathways likely corresponded to their expression in different organisms. The relatively high expression observed for both of these gene sets appeared to correspond to lower nitrate concentrations in the estuary, plume and adjacent coastal ocean compared to the freshwater baseline in April, November and June (discussed above) and to the (inferred) increase in nitrate uptake. Similar to our results, Glibert <italic>et al.</italic> <xref ref-type="bibr" rid="pone.0013312-Glibert1">[50]</xref> also found that nitrate uptake was relatively higher in April in the plume of the Chesapeake Bay estuary compared to other months during late spring, summer and winter.</p><p>Although both denitrification and DNRA are prominent processes in marine and estuarine sediments, expression of the genes involved in these pathways in our water samples was surprising because low levels of oxygen are required for cells to respire nitrate, and oxygen concentrations in the CRCM were relatively high during all analyzed seasons (<xref ref-type="fig" rid="pone-0013312-g002">Fig. 2A</xref>). We therefore hypothesize that denitrification and DNRA may be carried out by particle-attached bacteria that reside within suspended sediment particles, where oxygen may be limiting. Assuming microarray signals relate directly to enzymatic activity, our working hypothesis is that particle-attached organisms, which encounter low-oxygen microenvironments, expressed the genes for dissimilatory nitrate reduction in the water column, while free-living, or both particle-attached and free-living microorganisms, expressed nitrate assimilation genes. Alternatively, high expression of DNRA and denitrification genes may have resulted from residual activity of organisms transported from sub-oxic habitats, including estuarine sediments or the adjacent hyporheic zone. However, denitrifiers, at least, are often facultative anaerobes that do not require strictly anoxic conditions for respiration of nitrate <xref ref-type="bibr" rid="pone.0013312-Tiedje1">[51]</xref>, and the isolation of bacterial strains capable of denitrification under aerobic conditions has been reported <xref ref-type="bibr" rid="pone.0013312-Robertson1">[52]</xref>–<xref ref-type="bibr" rid="pone.0013312-Takaya1">[54]</xref>. Future biogeochemical rate measurements are needed to determine if these processes are occurring on suspended particulate matter in the Columbia River water column, a phenomenon that has been observed in the River Rhone plume and coastal waters of the northwestern Mediterranean Sea <xref ref-type="bibr" rid="pone.0013312-Michotey1">[55]</xref>, <xref ref-type="bibr" rid="pone.0013312-Omnes1">[56]</xref>.</p><p>Another probe subset with elevated signals corresponded to differential expression of genes involved in carbon utilization pathways (<xref ref-type="fig" rid="pone-0013312-g008">Fig. 8</xref>). These genes were relatively highly expressed in 13 out of 16 April samples, and also in 5 ocean samples from June. The gene group included carbon fixation genes encoding ribulose-1,5-bisphosphate carboxylase oxygenase (RuBisCO), the enzyme used to catalyze the first major step of carbon fixation in the Calvin cycle, and carbonic anhydrase (CA), the key enzyme of the CO<sub>2</sub> concentrating mechanism <xref ref-type="bibr" rid="pone.0013312-John1">[57]</xref>. Based on the relatively high expression observed in April samples of such genes involved in prokaryotic primary production, it appears that blooms of autotrophic prokaryotes may occur in the CRCM in spring, similar to the eukaryotic phytoplankton blooms that occur there at that time <xref ref-type="bibr" rid="pone.0013312-Small2">[10]</xref>. Consistent with this idea, a recent study showed that a spring phytoplankton bloom in the Sargasso Sea was accompanied by a bloom of Alphaproteobacteria <xref ref-type="bibr" rid="pone.0013312-Treusch1">[58]</xref>.</p><p>In addition, the highly expressed probe subset contained genes for the carbon starvation protein, CstA, and for the carbon storage regulator protein, CsrA. Both are key members of the global carbon storage regulatory (Csr) system involved in nutrient scavenging, utilization of alternative carbon sources, and quorum sensing <xref ref-type="bibr" rid="pone.0013312-Dubey1">[59]</xref>. Expression of the <italic>cstA</italic> gene, shown to encode a peptide transporter in <italic>Escherichia coli</italic> <xref ref-type="bibr" rid="pone.0013312-Schultz1">[60]</xref>, is consistent with reports that dissolved combined amino acids (DCAA) are components of estuarine and marine dissolved organic matter and are utilized as a source of carbon and nitrogen by heterotrophic bacteria <xref ref-type="bibr" rid="pone.0013312-Bronk1">[61]</xref>. Because coverage of genes involved in carbon metabolism was limited to a few groups in the current probe set, future iterations will include genes for more pathways, with the goal of examining the relationship between microbial gene expression and dissolved organic matter utilization in the CRCM.</p></sec><sec id="s2l"><title>Phylogenetic analysis of differential gene expression patterns</title><p>Consistent with our probe validation results, we analyzed the phylogenetic composition of the differentially expressed gene sets corresponding to clusters A through E (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>) at the family level (<xref ref-type="fig" rid="pone-0013312-g009">Fig. 9</xref>). For samples within each cluster we selected the probes that exhibited high signals relative to the freshwater baseline (the numbers of probes ranged from 135 to 578, depending on the cluster). Finally, we calculated the relative percentages of corresponding differentially expressed genes representing each microbial family. The relative representation of each family is also shown for the entire probe set, as well as for all probes producing signals above background (the 2076 probe set in <xref ref-type="fig" rid="pone-0013312-g004">Fig. 4</xref>). The composition of families represented in each cluster differed from the initial probe set and from one another. Although variability in the numbers of families represented by differentially-expressed genes in the different clusters was not large, there was a general trend toward higher diversity with increasing differential gene expression. Thus, cluster A contained 35 families and cluster B, 34 families with differentially expressed genes, while clusters C–E, showing more elevated levels of differential gene expression, contained 37, 39 and 40 families, respectively. Some families in the initial probe set were not detected in particular clusters. A striking example of this was the absence of methanogen families in cluster A. Of the five families of methanogens represented by the probe set, three were not detected in cluster A, and a fourth was reduced relative to its representation in the other clusters. In contrast, differential expression was detected from all five methanogen families in cluster E, and from four of the five families in clusters B, C and D (<xref ref-type="fig" rid="pone-0013312-g009">Fig. 9A</xref>). Interestingly, two other families from the top 40 representative families were not detected in cluster A, and one of these was Haloarcula, another member of the Euryarchaeota.</p><fig id="pone-0013312-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0013312.g009</object-id><label>Figure 9</label><caption><title>Phylogenetic analysis of differentially expressed microarray gene sets.</title><p>For each cluster (corresponding to clusters A–E in <xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>), all genes upregulated over 2-fold relative to the seasonal freshwater baseline in at least 2 samples were selected (for cluster A, 160; B, 135; C, 301; D, 372; and E, 575 genes). The percentages of different prokaryotic families within each gene set are shown by different colors within the bar graphs. To the left is shown family representation of the entire probe set and of the probe set corresponding to expressed genes. Only families represented by a minimum of 15 probes are shown, these comprise 40 of the 96 total families. (A) The top 20 microbial families represented by >61% of the 2226 probe set; and (B) the next 20 microbial families represented by >16% of all probes. Euryarchaeota families with relatively high differential gene expression in cluster E are indicated in blue.</p></caption><graphic xlink:href="pone.0013312.g009"></graphic></fig><p>The highest similarity in phylogenetic composition was observed between cluster D and E gene sets (<italic>R<sup>2</sup></italic> = 0.96), consistent with the sample clustering analysis described above (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>). In contrast to the other three clusters, D and E gene sets contained samples with relatively low differential expression of genes from Chloroflexaceae and Moraxellaceae. Distinguishing the two, however, was the relatively high differential expression of Euryarchaeota genes corresponding to several different families in cluster E (<xref ref-type="fig" rid="pone-0013312-g009">Fig. 9</xref>, shown in blue), and relatively low expression of Opitutaceae genes (data not shown). Finally, cluster B (August cluster) contained samples with relatively high differential expression of Chlorobiaceae, Oxalobacteraceae, and Aurantimonadaceae genes (<xref ref-type="fig" rid="pone-0013312-g009">Fig. 9</xref>). Taken together, ourdata indicated that gene expression changes relative to freshwater and seasons occurred in phylogenetically diverse microorganisms (<xref ref-type="fig" rid="pone-0013312-g005">Fig. 5</xref>, <xref ref-type="fig" rid="pone-0013312-g009">9</xref>).</p></sec><sec id="s2m"><title>Conclusions</title><p>In this study we employed a functional gene microarray probe set for detection of target genes and their closely related homologs in bacteria and archaea from environmental samples. This initial probe set was not intended to provide comprehensive coverage, since so little genome sequence information was available for CRCM microorganisms. Our aim, instead, was to test the efficacy of gene expression microarrays for detecting microbial responses to environmental change. Since many of the analyzed genes showed relatively high expression, coverage by the current probe set is likely to be representative of several important functional processes. Expansion of the probe set and additional improvements are planned for future work.</p><p>Application of multiple bioinformatic and experimental probe validation analyses consistently yielded 85–90% specificity for detection of designated targets, with 10–15% recognition of homologs from organisms within the same microbial family. Because of the likelihood that our probes hybridized to homologs of their target gene in at least some environmental samples, we limited our interpretations of trends in the data to those involving relatively large gene sets (e.g., <xref ref-type="fig" rid="pone-0013312-g006">Fig. 6</xref>, <xref ref-type="fig" rid="pone-0013312-g007">7</xref>, <xref ref-type="fig" rid="pone-0013312-g008">8</xref>). We have furthermore avoided conclusions relating functional processes to specific genera or species, and have instead analyzed gene expression profiles at higher taxonomic levels. If gene expression in the CRCM by members of a particular genus or species is of particular interest, additional analyses will need to be performed to provide exact identification. A more general interpretation of gene expression profiles, however, provided us with information about large-scale responses of prokaryotic assemblages during different seasons in the CRCM.</p><p>The strong seasonal shift observed in gene expression profiles between April and August is likely to be related to seasonal differences in nutrient inputs into the CRCM. In the Columbia River, as in other estuaries in this region <xref ref-type="bibr" rid="pone.0013312-Hickey1">[2]</xref>, <xref ref-type="bibr" rid="pone.0013312-Brown1">[33]</xref>, river flows dominate nitrogen inputs into the estuary in winter and spring, while upwelling has a greater influence in summer and fall. Our gene expression data imply that the microbial community response to these nutrient inputs may also be different. Interestingly, the upwelling index for August 2007 was relatively low (data not shown), potentially corresponding to the relatively low levels of differential gene expression observed in the estuary, plume and adjacent coastal ocean relative to freshwater during that time. Given that differential gene expression in April was comparatively elevated, these observations may suggest that higher nutrient inputs into the estuary contribute to elevated levels of differential gene expression by bacteria and archaea relative to expression upriver, while lower nutrient inputs result in more similarity in gene expression throughout the system. Multi-year analyses to assess variation in gene expression in response to variability in composition, quantity and source of nutrient inputs will be necessary to test this hypothesis.</p><p>Taken together, our results suggest that river-derived nutrients, which build up over winter, have large effects on gene expression of microbial assemblages in the estuary, plume and adjacent coastal waters in the spring when temperature and light become favorable for growth. Both April, prior to the spring freshet, and August in late summer were productive by a number of criteria, including measurements of total RNA, chl <italic>a</italic>, and overall gene expression. The presence of abundant microbial populations in the estuary in April coincided with elevated differential gene expression relative to the end-member freshwater samples and involved large numbers of probes (370–500). Microbial biomass and gene expression was also relatively high in the river and plume in August, but for reasons discussed above related to lower nutrient inputs and perhaps also because phytoplankton were abundant during this time, little differential gene expression was observed. Our data also suggest that generally unfavorable conditions in the estuary (winter conditions in November, and record high river discharge in June) resulted in relatively lower microarray signals and limited elevation of gene expression (observed only for 160 genes) for most of the corresponding samples. Thus, our data support the hypothesis that during periods of high riverflow both estuarine nutrient gradients and estuarine microbial communities are strongly influenced by river discharge <xref ref-type="bibr" rid="pone.0013312-Colbert1">[32]</xref>.</p><p>Our analysis of geochemical gradients indicated that no single factor was limiting for the development of productive microbial assemblages. Instead, microbial activity may have been determined by multiple combinations of environmental factors, such as described in a large body of literature discussing resource co-limitation in oceans (for review see <xref ref-type="bibr" rid="pone.0013312-Arrigo1">[34]</xref>). Nevertheless, we were able to discern certain trends in differential gene expression related to combinations of light, location, salinity, nutrient concentrations, and phytoplankton biomass. Somewhat surprising was the fact that, within the estuary, salinity by itself did not appear to more strongly influence gene expression profiles. This may be a consequence of adaptation by microbial populations to constantly fluctuating salinities resulting from the dynamic mixing of seawater and freshwater.</p><p>The seasonal and inter-annual variability of external forcing (e.g., <xref ref-type="supplementary-material" rid="pone.0013312.s005">Fig. S5</xref>) and circulation conditions in the Columbia River are too high to allow an exhaustive characterization of microbial gene expression based on two years of field campaigns. Continuing campaigns (e.g., 2009–2010) will help, but our ultimate goal is the development of in-situ continuous microbial observations. Sensors to accomplish these types of observations are being developed at CMOP <xref ref-type="bibr" rid="pone.0013312-Ghindilis2">[62]</xref> and may additionally be supplemented with advanced instrumentation developed elsewhere <xref ref-type="bibr" rid="pone.0013312-Preston1">[63]</xref>. Although the microarray profiles generated from this work are complex and their interpretation with respect to growth and activity of microbial populations in the CRCM will require ongoing analyses, our results indicate that environmental data and gene expression data may correspond in ways that facilitate our understanding of the drivers – physical, chemical, and microbial – of estuarine processes.</p></sec></sec><sec sec-type="materials|methods" id="s3"><title>Materials and Methods</title><sec id="s3a"><title>Collection of water samples</title><p>Water samples from the tidal freshwater, estuary, plume and adjacent coastal ocean were collected during four CMOP research cruises in August 2007, November 2007, April 2008, and June 2008, as described in detail in <ext-link ext-link-type="uri" xlink:href="http://www.stccmop.org/research/cruise">http://www.stccmop.org/research/cruise</ext-link>. Locations of the sampling stations are shown in the map in <xref ref-type="fig" rid="pone-0013312-g001">Fig. 1A</xref>. For each location, water was collected at the surface (1–3 m), and at a mid-depth ranging from 4 to 50 m depending on salinity, and, whenever possible, in association with a peak in chl <italic>a</italic> fluorescence. Samples were collected using 10-L Niskin bottles mounted to a CTD (conductivity-temperature-depth meter) rosette. All microorganisms present in the water column were collected by filtering water through Sterivex 0.22 µm filter units (Fisher Scientific, Waltham, MA). The volume of filtered water varied from 1 to 6 liters depending on rapidity of filter clogging (typically, 4–6 liters in the plume and coastal ocean, and 1–2 liters in the river), and was recorded for each sample. The filters were preserved in 2 mL of RNAlater reagent (Ambion, Austin, TX), frozen, and stored at -80°C.</p></sec><sec id="s3b"><title>Environmental data acquisition</title><p>During sample collection, water temperature, salinity, and depth were recorded with a Sea-Bird 911+ CTD (conductivity-temperature-depth) profiler (Sea-Bird Electronics Inc., Bellevue, WA). Dissolved oxygen SBE43 (Sea-Bird) and chl <italic>a</italic> fluorescence (WetStar, WET Labs, Philomath, OR) sensors provided additional data. Aliquots of each water sample were used to obtain chemical and biological data, including nutrient concentrations, bacterial production, and chl <italic>a</italic> and pheophytin <italic>a</italic> fluorescence. Selected chemical data for the analyzed samples sets are shown in <xref ref-type="fig" rid="pone-0013312-g001">Fig. 1</xref>–<xref ref-type="fig" rid="pone-0013312-g002">2</xref>, while the entire suite of measurements is shown in <xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>. Daily values of river discharge at the Bonneville Dam, near Stevenson, Oregon (<xref ref-type="fig" rid="pone-0013312-g001">Fig. 1B</xref>) were acquired from the U.S. Geologic Survey's National Stream Water Quality Network (<ext-link ext-link-type="uri" xlink:href="http://water.usgs.gov/nasqan/">http://water.usgs.gov/nasqan/</ext-link>). Macronutrients (nitrate + nitrite – referred to hereafter as nitrate, ammonium, silicic acid and phosphate) were measured on a Lachat QuikChem 8000 Flow Injection Analysis system using standard colorimetric methods <xref ref-type="bibr" rid="pone.0013312-Parsons1">[64]</xref>. Bacterial production was measured as the rate of incorporation of L-[<sup>3</sup>H]leucine (50 nmol/L final concentration) using methods described elsewhere <xref ref-type="bibr" rid="pone.0013312-Crump4">[39]</xref>. Dissolved organic carbon was measured by filtering 20 mL of water through a GFF filter (ø 25 mm, Whatman, Piscataway, NJ). The filters were analyzed by HPL analytical services as described <xref ref-type="bibr" rid="pone.0013312-Sugimura1">[65]</xref> using a TOC-5000 total organic carbon analyzer (Shimadzu, Columbia, MD). Particulate organic carbon (POC) was measured by filtering 100–300 mL of water onto a pre-combusted (12 hours at 500°C) GFF filter (ø 25 mm, Whatman) to collect suspended particulate matter for elemental analysis. The POC content on the acid-fumed filters was determined using a Carlo Erba NA-1500 Elemental Analyzer (Thermo Fisher Scientific, Waltham, MA), configured and operated as described <xref ref-type="bibr" rid="pone.0013312-Verardo1">[66]</xref>.</p></sec><sec id="s3c"><title>Total RNA isolation from filter samples</title><p>For each filtered water sample, total RNA was isolated from all microorganisms, both prokaryotic and eukaryotic, that were retained on 0.2 µm filters. The RNA isolation protocol described in Griffiths <italic>et al</italic> <xref ref-type="bibr" rid="pone.0013312-Griffiths1">[16]</xref> was modified in the following manner. Each Sterivex filter was removed from the holder, cut longitudinally and placed into 2 ml screw-top eppendorf tubes containing CTAB extraction buffer (5% CTAB, 0.8 M NaCl in 0.1 M phosphate buffer, pH 8.0) and 0.5 mL zirconia/silica beads (1∶1 mixture of 0.1 and 1 mm, BioSpec Products, Bartlesville, OK). Bead beating for microbial cell disruption was performed for 1 min using a FastPrep FP120 machine (Thermo Fisher Scientific, Waltham, MA). The filter strips were re-extracted 3 times to maximize RNA yield, after which the extracts were pooled together and supplemented with an equal volume of a 1∶1 mixture of phenol:chloroform, and 0.1 volume each of 10% SDS, and 10% sodium lauryl sarcosine. Sample extracts were incubated for 20 min at room temperature with intense shaking (80 rpm), followed by standard phenol:chloroform extraction and isopropanol precipitation. The extracted RNA was treated with RNase-free DNase (Qiagen, Germantown, MD), and then purified using the RNeasy Kit (Qiagen) as recommended by the manufacturer. The RNA content was measured spectrophotometrically with a Nanodrop 3300 instrument (Thermo Fisher Scientific), and RNA quality was evaluated by capillary electrophoresis using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). Prominent rRNA peaks (both prokaryotic and eukaryotic) were observed in all samples. RNA concentrations were normalized per liter of filtered water to generate the RNA values shown in <xref ref-type="fig" rid="pone-0013312-g003">Fig. 3C</xref>. We used several sets of replicate samples to analyze the reproducibility of our RNA isolation procedure and the variability in total RNA content introduced by sample handling during water collection, filtering, and on-board storage. In all cases, RNA concentrations calculated for independently processed replicate samples had low coefficients of variation (4 to 14%, data not shown). Thus, the normalized RNA concentrations were used as a measurement of total living microbial biomass in the samples.</p></sec><sec id="s3d"><title>Preparation of fluorescently labeled cDNA targets for microarray hybridization</title><p>Total RNA samples were converted into cDNA using the Superscript III Reverse Transcriptase Kit (Invitrogen, Carlsbad, CA) with random hexamer primers. RNA was removed with RNase H according to the manufacturer's protocol. The cDNA samples were purified using QIAQuick Nucleotide Removal Kit (Qiagen), and fluorescently labeled using LabelIT uArray Cy5 Kit (Mirus Bio, Madison WI) as described by the manufacturer. The Cy5-labeled cDNA samples were purified using Illustra ProbeQuant G-50 Micro Columns (GE Healthcare, Piscataway, NJ), dried in a SpeedVac centrifuge and applied for microarray hybridization.</p></sec><sec id="s3e"><title>Subtractive hybridization assays</title><p>Partial depletion of rRNA from total RNA was performed using MICROBExpress™ Bacterial mRNA Enrichment Kit (Ambion, Austin, TX) as recommended by the manufacturer. A pooled sample of total RNA was divided into two aliquots, one of which was subjected to rRNA subtractive hybridization and removal using universal rRNA-specific primers and magnetic beads. Target cDNA samples were prepared independently for both untreated and partially depleted aliquots and equal amounts of each target preparation were hybridized to replicate microarrays. Results indicated that the relative rRNA content of this aliquot decreased to approximately 80% (from >95%), as estimated by capillary electrophoresis (data not shown).</p></sec><sec id="s3f"><title>Selection of genes of interest and microarray design</title><p>Genes of interest were selected using the Integrated Microbial Genomes (IMG) 2.5 system from DOE Joint Genome Institute <xref ref-type="bibr" rid="pone.0013312-Markowitz1">[28]</xref>. Approximately 300 sequenced genomes from environmental bacteria and archaea (251 and 48, respectively) were used for keyword-based selection of genes from those genomes with well-defined functional annotation related to carbon and nitrogen metabolism (which made up 48% and 27% of the selected genes, respectively). The remaining 25% of selected genes encoded well-characterized enzymes of central metabolism (housekeeping functions) and a small number (37) of light perception genes.</p><p>Oligonucleotide probes were designed from the list of genes of interest using Probe Weaver software (CombiMatrix Corporation, Mukilteo, WA). The recommended standard settings were used for selection of 35- to 45-nt oligomers in the sense orientation, with the melting temperature of target:probe hybrids varying from 70°C to 75°C. The software was used to select probes based on their specificity within the submitted gene set. Thus, to prevent cross-hybridization of the selected probes with rRNA, 146 and 882 rRNA genes from eukaryotic and prokaryotic microorganisms, respectively, were added to the submitted gene list. These rRNA sequences were not used to design probes, but rather to filter the selected probes against cross-hybridization. After completion of the probe design, the probe sets were synthesized in situ on the surface of CombiMatrix CustomArray™ 12K and 4X2K oligonucleotide microarrays (described in <ext-link ext-link-type="uri" xlink:href="http://www.combimatrix.com">www.combimatrix.com</ext-link>).</p></sec><sec id="s3g"><title>Microarray hybridization and re-use</title><p>Hybridization of CustomArray™ microarrays was performed as recommended by the manufacturer (<ext-link ext-link-type="uri" xlink:href="http://www.combimatrix.com/docs/PTL005_00_4x2K_Hyb_Imaging.pdf">http://www.combimatrix.com/docs/PTL005_00_4x2K_Hyb_Imaging.pdf</ext-link>), using the single-color experimental scheme. The microarrays were hybridized at 46°C for 16 hours, using Cy5-labeled cDNA targets, 3–4 µg per microarray for the 12K format, or 1.5–2 µg per sector for the 4X2K format. Hybridization images were obtained using a ScanArray 4000 fluorescent scanner (PerkinElmer, Waltham, MA) at 5 micron resolution, and quantified with Microarray Imager software (CombiMatrix). For microarray re-use, the labeled targets were chemically denatured and removed from microarray probes as described in the CombiMatrix stripping protocol (<ext-link ext-link-type="uri" xlink:href="http://www.combimatrix.com/docs/PTL002_01_4x2K_StrippingReHyb.pdf">http://www.combimatrix.com/docs/PTL002_01_4x2K_StrippingReHyb.pdf</ext-link>). The quality of target removal was evaluated by scanning at 5 micron resolution. In total, each microarray was used four times, as recommended by the manufacturer.</p></sec><sec id="s3h"><title>Microarray data analysis and deposition</title><p>Data normalization based on the quantile algorithm was done using the Probe Weaver software from CombiMatrix. The background signal value was calculated as the mean intensity of the lowest 5% of all signals. Genes were selected as significantly expressed if the corresponding probes produced normalized signal intensities exceeding the background value plus 3X standard deviation in at least 2 different samples. Gene expression ratios for each seasonal sample set were calculated versus corresponding freshwater river baselines. For each gene represented on the microarrays, the baseline was calculated as the average of the normalized signal intensities in two independent samples collected at the Beaver Army Terminal (River Mile 53, <xref ref-type="fig" rid="pone-0013312-g001">Fig. 1A</xref>, the corresponding samples are shown in bold in <xref ref-type="supplementary-material" rid="pone.0013312.s006">Table S1</xref>). Ratios of differential gene expression were then calculated for each gene in the remaining 56 samples as log<sub>2</sub> (signal intensity of the gene in the sample) minus log<sub>2</sub> (intensity of the gene in the corresponding seasonal river baseline). Finally, to select for statistically significant changes, genes were designated as ‘differentially expressed’ if the corresponding ratios varied by at least 2-fold in at least 2 out of 56 samples (consistent with recommended criteria <xref ref-type="bibr" rid="pone.0013312-Quackenbush1">[67]</xref>, <xref ref-type="bibr" rid="pone.0013312-Schena1">[68]</xref>). Scatter plot analysis of duplicate hybridizations, calculation of pair-wise correlation coefficients among samples, and ANOVA were done using NCI-supported software BRB-ArrayTools (<ext-link ext-link-type="uri" xlink:href="http://linus.nci.nih.gov/pilot/index.html">http://linus.nci.nih.gov/pilot/index.html</ext-link>). Clustering and self-organizing map (SOM) analyses were performed using Cluster and TreeView software (<ext-link ext-link-type="uri" xlink:href="http://rana.lbl.gov/eisen/">http://rana.lbl.gov/eisen/</ext-link><xref ref-type="bibr" rid="pone.0013312-Eisen1">[69]</xref>). The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus <xref ref-type="bibr" rid="pone.0013312-Edgar1">[70]</xref>, and are accessible through GEO Series accession number GSE18303 (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE18303">http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE18303</ext-link>).</p></sec></sec><sec sec-type="supplementary-material" id="s4"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0013312.s001"><label>Figure S1</label><caption><p>Validation of microarray probe specificity with Bacillus subtilis RNA. Total RNA was isolated from laboratory cultures of Bacillus subtilis, converted to fluorescently labeled cDNA, and hybridized in duplicate to microarrays. (A) Representative plot of the raw signal intensity data sorted by organism annotations; (B) raw signal intensity data for a subset of probes corresponding to Bacilli/Bacillales probes only.</p><p>(0.18 MB TIF)</p></caption><media xlink:href="pone.0013312.s001.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0013312.s002"><label>Figure S2</label><caption><p>Experimental variation and reproducibility of microarray data. Microarray data were analyzed using scatter plots of replicate experiments. Each data point on a scatter plot represents an individual probe. X and Y values are the normalized signal intensities (log scale) in the first and second hybridizations, respectively. The angled lines show the two-fold cut-offs for the ratios between X and Y values. Data points located outside of the lines show significant (over 2-fold) difference between values obtained in the first and second hybridizations, and they are indicated with closed circles. The data points located within the cut-off lines are indicated with open circles.</p><p>(0.82 MB TIF)</p></caption><media xlink:href="pone.0013312.s002.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0013312.s003"><label>Figure S3</label><caption><p>Pair-wise correlation (Spearman rank) coefficients calculated for four sampling seasons. Within each season,microarray profiles for each sample were compared to those for the corresponding seasonal freshwater end-member (collected at Beaver Army Dock, <xref ref-type="fig" rid="pone-0013312-g001">Fig. 1A</xref>). The resulting correlation coefficients were plotted according to water salinities. Black and gray squares represent ocean and river samples, respectively. Trend lines were generated in Excel.</p><p>(0.12 MB TIF)</p></caption><media xlink:href="pone.0013312.s003.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0013312.s004"><label>Figure S4</label><caption><p>Multidimensional scaling (MDS) plots of maximum variability among samples. MDS plots were constructed using Primer6 software. Samples are displayed according to sampling time (A, April; J, June; AG, August; N, November).</p><p>(0.31 MB TIF)</p></caption><media xlink:href="pone.0013312.s004.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0013312.s005"><label>Figure S5</label><caption><p>Seasonal and inter-annual variability of Columbia River discharges. Recent years (1997 and 2001) with extreme discharge levels are shown for context relative to the years of CMOP campaigns to date (2007–2009; campaign periods marked at the bottom of the graph). In light gray, discharges for all other years since 1997 are shown. Inter-annual variations are most marked in winter and during spring freshets (May-June), but are present throughout the year. Q, river discharge in m3s-1.</p><p>(0.20 MB TIF)</p></caption><media xlink:href="pone.0013312.s005.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0013312.s006"><label>Table S1</label><caption><p>Corresponding physical, chemical and biological characteristics of the seasonal sample sets analyzed by microarray analysis. Location codes: number shows distance from the coast in km; CR, Columbia River transect in the plume and coastal ocean; NH, Newport Hydroline transect in the coastal ocean at Newport, Oregon; AST and HAM, CRE locations near Astoria (river mile 7–9) and Hammond (river mile 5), respectively; TID, CRE locations in the tidal basin (river mile 22–23); BA, river location at Beaver Army Dock (river mile 53) near Quincy, Oregon; UP, river at mile 74. River baseline samples are shown in bold. Leucine incorp., leucine incorporation by heterotrophic plankton.</p><p>(0.64 MB TIF)</p></caption><media xlink:href="pone.0013312.s006.tif" mimetype="image" mime-subtype="tiff"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>We gratefully acknowledge the chief scientists of CMOP cruises, and the captains and crews of the R/V <italic>Wecoma</italic> and the R/V <italic>Barnes</italic> for their assistance on this research. We thank P. Kahn and R. Letelier for chl <italic>a</italic> analysis, J. Jennings for nutrient analysis, F. Prahl for pigment and POC data, and University of Maryland Center for Environmental Science for DOC analyses. We are very grateful to K. Geszvain and M. Nakano for RNA samples from bacterial cultures. We are also grateful to the CMOP Cyberinfrastructure team for providing access to environmental data.</p></ack><fn-group><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The research reported in this article was funded by a grant from the National Science Foundation (NSF grant OCE 0424602); <ext-link ext-link-type="uri" xlink:href="http://www.nsf.gov/">http://www.nsf.gov/</ext-link>. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group><ref-list><title>References</title><ref id="pone.0013312-Simenstad1"><label>1</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simenstad</surname><given-names>CA</given-names></name><name><surname>Small</surname><given-names>LF</given-names></name><name><surname>David McIntire</surname><given-names>C</given-names></name><name><surname>Jay</surname><given-names>DA</given-names></name><name><surname>Sherwood</surname><given-names>C</given-names></name></person-group><year>1990</year><article-title>Columbia river estuary studies: An introduction to the estuary, a brief history, and prior studies.</article-title><source>Progress In Oceanography</source><volume>25</volume><fpage>1</fpage><lpage>13</lpage></mixed-citation></ref><ref id="pone.0013312-Hickey1"><label>2</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hickey</surname><given-names>B</given-names></name><name><surname>Banas</surname><given-names>N</given-names></name></person-group><year>2003</year><article-title>Oceanography of the U.S. Pacific Northwest Coastal Ocean and estuaries with application to coastal ecology.</article-title><source>Estuaries and Coasts</source><volume>26</volume><fpage>1010</fpage><lpage>1031</lpage></mixed-citation></ref><ref id="pone.0013312-Barnes1"><label>3</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Barnes</surname><given-names>C</given-names></name><name><surname>Duxbury</surname><given-names>AC</given-names></name><name><surname>Morse</surname><given-names>B</given-names></name></person-group><year>1972</year><article-title>Circulation and selected properties of the Columbia River effluent at sea.</article-title><person-group person-group-type="editor"><name><surname>Pruter</surname><given-names>T</given-names></name><name><surname>Alverson</surname><given-names>DL</given-names></name></person-group><source>The Columbia River estuary and adjacent ocean waters: bioenvironmental studies</source><publisher-loc>Seattle</publisher-loc><publisher-name>University of Washington Press</publisher-name><fpage>71</fpage><lpage>80</lpage></mixed-citation></ref><ref id="pone.0013312-Bruland1"><label>4</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bruland</surname><given-names>KW</given-names></name><name><surname>Lohan</surname><given-names>MC</given-names></name><name><surname>Aguilar-Islas</surname><given-names>AM</given-names></name><name><surname>Smith</surname><given-names>GJ</given-names></name><name><surname>Sohst</surname><given-names>B</given-names></name><etal></etal></person-group><year>2008</year><article-title>Factors influencing the chemistry of the near-field Columbia River plume: nitrate, silicic acid, dissolved Fe, and dissolved Mn.</article-title><source>J Geophys Res</source><volume>113</volume><fpage>C00B02</fpage></mixed-citation></ref><ref id="pone.0013312-Jay1"><label>5</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jay</surname><given-names>DA</given-names></name><name><surname>Giese</surname><given-names>BS</given-names></name><name><surname>Sherwood</surname><given-names>CR</given-names></name></person-group><year>1990</year><article-title>Energetics and sedimentary processes in the Columbia River Estuary.</article-title><source>Progress In Oceanography</source><volume>25</volume><fpage>157</fpage><lpage>174</lpage></mixed-citation></ref><ref id="pone.0013312-Jay2"><label>6</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jay</surname><given-names>D</given-names></name><name><surname>Uncles</surname><given-names>R</given-names></name><name><surname>Largeir</surname><given-names>J</given-names></name><name><surname>Geyer</surname><given-names>W</given-names></name><name><surname>Vallino</surname><given-names>J</given-names></name><etal></etal></person-group><year>1997</year><article-title>A review of recent developments in estuarine scalar flux estimation.</article-title><source>Estuaries and Coasts</source><volume>20</volume><fpage>262</fpage><lpage>280</lpage></mixed-citation></ref><ref id="pone.0013312-Chawla1"><label>7</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chawla</surname><given-names>A</given-names></name><name><surname>Jay</surname><given-names>D</given-names></name><name><surname>Baptista</surname><given-names>A</given-names></name><name><surname>Wilkin</surname><given-names>M</given-names></name><name><surname>Seaton</surname><given-names>C</given-names></name></person-group><year>2008</year><article-title>Seasonal variability and estuary–shelf interactions in circulation dynamics of a river-dominated estuary.</article-title><source>Estuaries and Coasts</source><volume>31</volume><fpage>269</fpage><lpage>288</lpage></mixed-citation></ref><ref id="pone.0013312-Small1"><label>8</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Small</surname><given-names>LF</given-names></name><name><surname>Prahl</surname><given-names>F</given-names></name></person-group><year>2004</year><article-title>A particle conveyor belt process in the Columbia River estuary: Evidence from chlorophyll a and particulate organic carbon.</article-title><source>Estuaries and Coasts</source><volume>27</volume><fpage>999</fpage><lpage>1013</lpage></mixed-citation></ref><ref id="pone.0013312-Giovannoni1"><label>9</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Giovannoni</surname><given-names>S</given-names></name><name><surname>Rappe</surname><given-names>MS</given-names></name></person-group><year>2000</year><article-title>Evolution diversity and molecular ecology of marine prokaryotes.</article-title><person-group person-group-type="editor"><name><surname>Kirchman</surname><given-names>D</given-names></name></person-group><source>Microbial ecology of the oceans</source><publisher-loc>New York</publisher-loc><publisher-name>Wiley-Liss</publisher-name><fpage>47</fpage><lpage>84</lpage></mixed-citation></ref><ref id="pone.0013312-Small2"><label>10</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Small</surname><given-names>LF</given-names></name><name><surname>McIntire</surname><given-names>CD</given-names></name><name><surname>MacDonald</surname><given-names>KB</given-names></name><name><surname>Lara-Lara</surname><given-names>JR</given-names></name><name><surname>Frey</surname><given-names>BE</given-names></name><etal></etal></person-group><year>1990</year><article-title>Primary production, plant and detrital biomass, and particle transport in the Columbia River Estuary.</article-title><source>Progress In Oceanography</source><volume>25</volume><fpage>175</fpage><lpage>210</lpage></mixed-citation></ref><ref id="pone.0013312-Baross1"><label>11</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Baross</surname><given-names>JA</given-names></name><name><surname>Crump</surname><given-names>BC</given-names></name><name><surname>Simenstad</surname><given-names>CA</given-names></name></person-group><year>1994</year><article-title>Elevated ‘microbial loop’ activities in the Columbia River estuary turbidity maximum.</article-title><person-group person-group-type="editor"><name><surname>Dyer</surname><given-names>K</given-names></name><name><surname>Orth</surname><given-names>R</given-names></name></person-group><source>Changing particle fluxes in estuaries: implications from science to management</source><publisher-loc>Friedensborg</publisher-loc><publisher-name>Olsen and Olsen Press</publisher-name><fpage>459</fpage><lpage>464</lpage></mixed-citation></ref><ref id="pone.0013312-Prahl1"><label>12</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Prahl</surname><given-names>F</given-names></name><name><surname>Small</surname><given-names>LF</given-names></name><name><surname>Eversmeyer</surname><given-names>B</given-names></name></person-group><year>1997</year><article-title>Biogeochemical characterization of suspended particulate matter in the Columbia River estuary.</article-title><source>Marine Ecology Progress Series</source><volume>160</volume><fpage>173</fpage><lpage>184</lpage></mixed-citation></ref><ref id="pone.0013312-Crump1"><label>13</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crump</surname><given-names>BC</given-names></name><name><surname>Baross</surname><given-names>JA</given-names></name><name><surname>Simenstad</surname><given-names>CA</given-names></name></person-group><year>1998</year><article-title>Dominance of particle-attached bacteria in the Columbia River estuary, USA.</article-title><source>Aquatic Microbial Ecology</source><volume>14</volume><fpage>7</fpage><lpage>18</lpage></mixed-citation></ref><ref id="pone.0013312-Crump2"><label>14</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crump</surname><given-names>BC</given-names></name><name><surname>Armbrust</surname><given-names>EV</given-names></name><name><surname>Baross</surname><given-names>JA</given-names></name></person-group><year>1999</year><article-title>Phylogenetic analysis of particle-attached and free-living bacterial communities in the Columbia river, its estuary, and the adjacent coastal ocean.</article-title><source>Appl Environ Microbiol</source><volume>65</volume><fpage>3192</fpage><lpage>3204</lpage><pub-id pub-id-type="pmid">10388721</pub-id></mixed-citation></ref><ref id="pone.0013312-Crump3"><label>15</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crump</surname><given-names>BC</given-names></name><name><surname>Baross</surname><given-names>JA</given-names></name></person-group><year>2000</year><article-title>Archaeaplankton in the Columbia River, its estuary and the adjacent coastal ocean, USA.</article-title><source>FEMS Microbiol Ecol</source><volume>31</volume><fpage>231</fpage><lpage>239</lpage><pub-id pub-id-type="pmid">10719204</pub-id></mixed-citation></ref><ref id="pone.0013312-Griffiths1"><label>16</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Griffiths</surname><given-names>RI</given-names></name><name><surname>Whiteley</surname><given-names>AS</given-names></name><name><surname>O'Donnell</surname><given-names>AG</given-names></name><name><surname>Bailey</surname><given-names>MJ</given-names></name></person-group><year>2000</year><article-title>Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA- and rRNA-based microbial community composition.</article-title><source>Appl Environ Microbiol</source><volume>66</volume><fpage>5488</fpage><lpage>5491</lpage><pub-id pub-id-type="pmid">11097934</pub-id></mixed-citation></ref><ref id="pone.0013312-Grant1"><label>17</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grant</surname><given-names>S</given-names></name><name><surname>Grant</surname><given-names>WD</given-names></name><name><surname>Cowan</surname><given-names>DA</given-names></name><name><surname>Jones</surname><given-names>BE</given-names></name><name><surname>Ma</surname><given-names>Y</given-names></name><etal></etal></person-group><year>2006</year><article-title>Identification of eukaryotic open reading frames in metagenomic cDNA libraries made from environmental samples.</article-title><source>Appl Environ Microbiol</source><volume>72</volume><fpage>135</fpage><lpage>143</lpage><pub-id pub-id-type="pmid">16391035</pub-id></mixed-citation></ref><ref id="pone.0013312-FriasLopez1"><label>18</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Frias-Lopez</surname><given-names>J</given-names></name><name><surname>Shi</surname><given-names>Y</given-names></name><name><surname>Tyson</surname><given-names>GW</given-names></name><name><surname>Coleman</surname><given-names>ML</given-names></name><name><surname>Schuster</surname><given-names>SC</given-names></name><etal></etal></person-group><year>2008</year><article-title>Microbial community gene expression in ocean surface waters.</article-title><source>Proc Natl Acad Sci U S A</source><volume>105</volume><fpage>3805</fpage><lpage>3810</lpage><pub-id pub-id-type="pmid">18316740</pub-id></mixed-citation></ref><ref id="pone.0013312-Poretsky1"><label>19</label><mixed-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><etal></etal></person-group><year>2009</year><article-title>Comparative day/night metatranscriptomic analysis of microbial communities in the North Pacific subtropical gyre.</article-title><source>Environ Microbiol</source><volume>11</volume><fpage>1358</fpage><lpage>1375</lpage><pub-id pub-id-type="pmid">19207571</pub-id></mixed-citation></ref><ref id="pone.0013312-Zhou1"><label>20</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>J</given-names></name><name><surname>Thompson</surname><given-names>DK</given-names></name></person-group><year>2002</year><article-title>Challenges in applying microarrays to environmental studies.</article-title><source>Current Opinion in Biotechnology</source><volume>13</volume><fpage>204</fpage><lpage>207</lpage><pub-id pub-id-type="pmid">12180093</pub-id></mixed-citation></ref><ref id="pone.0013312-TaroncherOldenburg1"><label>21</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Taroncher-Oldenburg</surname><given-names>G</given-names></name><name><surname>Griner</surname><given-names>EM</given-names></name><name><surname>Francis</surname><given-names>CA</given-names></name><name><surname>Ward</surname><given-names>BB</given-names></name></person-group><year>2003</year><article-title>Oligonucleotide microarray for the study of functional gene diversity in the nitrogen cycle in the environment.</article-title><source>Applied and Environmental Microbiology</source><volume>69</volume><fpage>1159</fpage><lpage>1171</lpage><pub-id pub-id-type="pmid">12571043</pub-id></mixed-citation></ref><ref id="pone.0013312-Parro1"><label>22</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parro</surname><given-names>V</given-names></name><name><surname>Moreno-Paz</surname><given-names>M</given-names></name><name><surname>Gonzalez-Toril</surname><given-names>E</given-names></name></person-group><year>2007</year><article-title>Analysis of environmental transcriptomes by DNA microarrays.</article-title><source>Environ Microbiol</source><volume>9</volume><fpage>453</fpage><lpage>464</lpage><pub-id pub-id-type="pmid">17222143</pub-id></mixed-citation></ref><ref id="pone.0013312-Lodes1"><label>23</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lodes</surname><given-names>MJ</given-names></name><name><surname>Suciu</surname><given-names>D</given-names></name><name><surname>Elliott</surname><given-names>M</given-names></name><name><surname>Stover</surname><given-names>AG</given-names></name><name><surname>Ross</surname><given-names>M</given-names></name><etal></etal></person-group><year>2006</year><article-title>Use of Semiconductor-Based Oligonucleotide Microarrays for Influenza A Virus Subtype Identification and Sequencing.</article-title><source>J Clinical Microbiol</source><volume>44</volume><fpage>1209</fpage><lpage>1218</lpage><pub-id pub-id-type="pmid">16597840</pub-id></mixed-citation></ref><ref id="pone.0013312-Ghindilis1"><label>24</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghindilis</surname><given-names>AL</given-names></name><name><surname>Smith</surname><given-names>MW</given-names></name><name><surname>Schwarzkopf</surname><given-names>KR</given-names></name><name><surname>Roth</surname><given-names>KM</given-names></name><name><surname>Peyvan</surname><given-names>K</given-names></name><etal></etal></person-group><year>2007</year><article-title>CombiMatrix oligonucleotide arrays: genotyping and gene expression assays employing electrochemical detection.</article-title><source>Biosens Bioelectron</source><volume>22</volume><fpage>1853</fpage><lpage>1860</lpage><pub-id pub-id-type="pmid">16891109</pub-id></mixed-citation></ref><ref id="pone.0013312-Bolotin1"><label>25</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bolotin</surname><given-names>S</given-names></name><name><surname>Lombos</surname><given-names>E</given-names></name><name><surname>Yeung</surname><given-names>R</given-names></name><name><surname>Eshaghi</surname><given-names>A</given-names></name><name><surname>Blair</surname><given-names>J</given-names></name><etal></etal></person-group><year>2009</year><article-title>Verification of the Combimatrix influenza detection assay for the detection of influenza A subtype during the 2007-2008 influenza season in Toronto, Canada.</article-title><source>Virol J</source><volume>6</volume><fpage>37</fpage><lpage>45</lpage><pub-id pub-id-type="pmid">19321005</pub-id></mixed-citation></ref><ref id="pone.0013312-He1"><label>26</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>He</surname><given-names>S</given-names></name><name><surname>Kunin</surname><given-names>V</given-names></name><name><surname>Haynes</surname><given-names>M</given-names></name><name><surname>Martin</surname><given-names>HG</given-names></name><name><surname>Ivanova</surname><given-names>N</given-names></name><etal></etal></person-group><year>2010</year><article-title>Metatranscriptomic array analysis of ‘Candidatus Accumulibacter phosphatis’-enriched enhanced biological phosphorus removal sludge.</article-title><source>Environ Microbiol</source><volume>12</volume><fpage>1205</fpage><lpage>1217</lpage><pub-id pub-id-type="pmid">20148930</pub-id></mixed-citation></ref><ref id="pone.0013312-Liu1"><label>27</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>RH</given-names></name><name><surname>Munro</surname><given-names>SB</given-names></name><name><surname>Nguyen</surname><given-names>T</given-names></name><name><surname>Siuda</surname><given-names>T</given-names></name><name><surname>Suciu</surname><given-names>D</given-names></name><etal></etal></person-group><year>2006</year><article-title>Integrated microfluidic CustomArray device for bacterial genotyping and identification.</article-title><source>Journal of the Association for Laboratory Automation</source><volume>11</volume><fpage>360</fpage><lpage>367</lpage></mixed-citation></ref><ref id="pone.0013312-Markowitz1"><label>28</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Markowitz</surname><given-names>VM</given-names></name><name><surname>Szeto</surname><given-names>E</given-names></name><name><surname>Palaniappan</surname><given-names>K</given-names></name><name><surname>Grechkin</surname><given-names>Y</given-names></name><name><surname>Chu</surname><given-names>K</given-names></name><etal></etal></person-group><year>2008</year><article-title>The integrated microbial genomes (IMG) system in 2007: data content and analysis tool extensions.</article-title><source>Nucleic Acids Research</source><volume>36</volume><fpage>D528</fpage><lpage>533</lpage><pub-id pub-id-type="pmid">17933782</pub-id></mixed-citation></ref><ref id="pone.0013312-AguilarIslas1"><label>29</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aguilar-Islas</surname><given-names>AM</given-names></name><name><surname>Bruland</surname><given-names>KW</given-names></name></person-group><year>2006</year><article-title>Dissolved manganese and silicic acid in the Columbia River plume: A major source to the California current and coastal waters off Washington and Oregon.</article-title><source>Marine Chemistry</source><volume>101</volume><fpage>233</fpage><lpage>247</lpage></mixed-citation></ref><ref id="pone.0013312-Kudela1"><label>30</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kudela</surname><given-names>RM</given-names></name><name><surname>Peterson</surname><given-names>TD</given-names></name></person-group><year>2009</year><article-title>Influence of a buoyant river plume on phytoplankton nutrient dynamics: What controls standing stocks and productivity?</article-title><source>J Geophys Res</source><volume>114</volume><fpage>C00B11</fpage></mixed-citation></ref><ref id="pone.0013312-Turner1"><label>31</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turner</surname><given-names>RE</given-names></name><name><surname>Rabalais</surname><given-names>NN</given-names></name><name><surname>Justic</surname><given-names>D</given-names></name></person-group><year>2008</year><article-title>Gulf of Mexico hypoxia: alternate states and a legacy.</article-title><source>Environmental Science & Technology</source><volume>42</volume><fpage>2323</fpage><lpage>2327</lpage><pub-id pub-id-type="pmid">18504960</pub-id></mixed-citation></ref><ref id="pone.0013312-Colbert1"><label>32</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colbert</surname><given-names>D</given-names></name><name><surname>McManus</surname><given-names>J</given-names></name></person-group><year>2003</year><article-title>Nutrient biogeochemistry in an upwelling-influenced estuary of the Pacific Northwest (Tillamook Bay, Oregon, USA).</article-title><source>Estuaries and Coasts</source><volume>26</volume><fpage>1205</fpage><lpage>1219</lpage></mixed-citation></ref><ref id="pone.0013312-Brown1"><label>33</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>C</given-names></name><name><surname>Ozretich</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Coupling between the coastal ocean and Yaquina Bay, Oregon: importance of oceanic inputs relative to other nitrogen sources.</article-title><source>Estuaries and Coasts</source><volume>32</volume><fpage>219</fpage><lpage>237</lpage></mixed-citation></ref><ref id="pone.0013312-Arrigo1"><label>34</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arrigo</surname><given-names>KR</given-names></name></person-group><year>2005</year><article-title>Marine microorganisms and global nutrient cycles.</article-title><source>Nature</source><volume>437</volume><fpage>349</fpage><lpage>355</lpage><pub-id pub-id-type="pmid">16163345</pub-id></mixed-citation></ref><ref id="pone.0013312-Conomos1"><label>35</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Conomos</surname><given-names>T</given-names></name><name><surname>Gross</surname><given-names>MG</given-names></name><name><surname>Barnes</surname><given-names>CA</given-names></name><name><surname>Richards</surname><given-names>FA</given-names></name></person-group><year>1972</year><article-title>River-ocean nutrient relations in summer.</article-title><person-group person-group-type="editor"><name><surname>Pruter</surname><given-names>A</given-names></name><name><surname>Alverson</surname><given-names>DL</given-names></name></person-group><source>The Columbia River and adjacent ocean waters</source><publisher-loc>Seattle</publisher-loc><publisher-name>University of Washington Press</publisher-name><fpage>151</fpage><lpage>175</lpage></mixed-citation></ref><ref id="pone.0013312-Arhonditsis1"><label>36</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arhonditsis</surname><given-names>G</given-names></name><name><surname>Brett</surname><given-names>MT</given-names></name><name><surname>Frodge</surname><given-names>J</given-names></name></person-group><year>2003</year><article-title>Environmental control and limnological impacts of a large recurrent spring bloom in Lake Washington, USA.</article-title><source>Environ Manage</source><volume>31</volume><fpage>603</fpage><lpage>618</lpage><pub-id pub-id-type="pmid">12719891</pub-id></mixed-citation></ref><ref id="pone.0013312-Sala1"><label>37</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sala</surname><given-names>MM</given-names></name><name><surname>Terrado</surname><given-names>R</given-names></name><name><surname>Lovejoy</surname><given-names>C</given-names></name><name><surname>Unrein</surname><given-names>F</given-names></name><name><surname>Pedros-Alio</surname><given-names>C</given-names></name></person-group><year>2008</year><article-title>Metabolic diversity of heterotrophic bacterioplankton over winter and spring in the coastal Arctic Ocean.</article-title><source>Environ Microbiol</source><volume>10</volume><fpage>942</fpage><lpage>949</lpage><pub-id pub-id-type="pmid">18218033</pub-id></mixed-citation></ref><ref id="pone.0013312-Teira1"><label>38</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Teira</surname><given-names>E</given-names></name><name><surname>Gasol</surname><given-names>JM</given-names></name><name><surname>Aranguren-Gassis</surname><given-names>M</given-names></name><name><surname>Fernandez</surname><given-names>A</given-names></name><name><surname>Gonzalez</surname><given-names>J</given-names></name><etal></etal></person-group><year>2008</year><article-title>Linkages between bacterioplankton community composition, heterotrophic carbon cycling and environmental conditions in a highly dynamic coastal ecosystem.</article-title><source>Environ Microbiol</source><volume>10</volume><fpage>906</fpage><lpage>917</lpage><pub-id pub-id-type="pmid">18215158</pub-id></mixed-citation></ref><ref id="pone.0013312-Crump4"><label>39</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crump</surname><given-names>BC</given-names></name><name><surname>Hopkinson</surname><given-names>CS</given-names></name><name><surname>Sogin</surname><given-names>ML</given-names></name><name><surname>Hobbie</surname><given-names>JE</given-names></name></person-group><year>2004</year><article-title>Microbial biogeography along an estuarine salinity gradient: combined influences of bacterial growth and residence time.</article-title><source>Appl Environ Microbiol</source><volume>70</volume><fpage>1494</fpage><lpage>1505</lpage><pub-id pub-id-type="pmid">15006771</pub-id></mixed-citation></ref><ref id="pone.0013312-Lodes2"><label>40</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lodes</surname><given-names>MJ</given-names></name><name><surname>Caraballo</surname><given-names>M</given-names></name><name><surname>Suciu</surname><given-names>D</given-names></name><name><surname>Munro</surname><given-names>S</given-names></name><name><surname>Kumar</surname><given-names>A</given-names></name><etal></etal></person-group><year>2009</year><article-title>Detection of cancer with serum miRNAs on an oligonucleotide microarray.</article-title><source>PLoS One</source><volume>4</volume><fpage>e6229</fpage><pub-id pub-id-type="pmid">19597549</pub-id></mixed-citation></ref><ref id="pone.0013312-Lodes3"><label>41</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lodes</surname><given-names>MJ</given-names></name><name><surname>Suciu</surname><given-names>D</given-names></name><name><surname>Wilmoth</surname><given-names>JL</given-names></name><name><surname>Ross</surname><given-names>M</given-names></name><name><surname>Munro</surname><given-names>S</given-names></name><etal></etal></person-group><year>2007</year><article-title>Identification of upper respiratory tract pathogens using electrochemical detection on an oligonucleotide microarray.</article-title><source>PLoS One</source><volume>2</volume><fpage>e924</fpage><pub-id pub-id-type="pmid">17895966</pub-id></mixed-citation></ref><ref id="pone.0013312-Maurer1"><label>42</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maurer</surname><given-names>K</given-names></name><name><surname>Cooper</surname><given-names>J</given-names></name><name><surname>Caraballo</surname><given-names>M</given-names></name><name><surname>Crye</surname><given-names>J</given-names></name><name><surname>Suciu</surname><given-names>D</given-names></name><etal></etal></person-group><year>2006</year><article-title>Electrochemically generated acid and its containment to 100 micron reaction areas for the production of DNA microarrays.</article-title><source>PLoS One</source><volume>1</volume><fpage>e34</fpage><pub-id pub-id-type="pmid">17183662</pub-id></mixed-citation></ref><ref id="pone.0013312-Allawi1"><label>43</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Allawi</surname><given-names>HT</given-names></name><name><surname>SantaLucia</surname><given-names>J</given-names></name></person-group><year>1998</year><article-title>Nearest neighbor thermodynamic parameters for internal G.A mismatches in DNA.</article-title><source>Biochemistry</source><volume>37</volume><fpage>2170</fpage><lpage>2179</lpage><pub-id pub-id-type="pmid">9485363</pub-id></mixed-citation></ref><ref id="pone.0013312-Allawi2"><label>44</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Allawi</surname><given-names>HT</given-names></name><name><surname>SantaLucia</surname><given-names>J</given-names></name></person-group><year>1998</year><article-title>Nearest-neighbor thermodynamics of internal A.C mismatches in DNA: sequence dependence and pH Effects.</article-title><source>Biochemistry</source><volume>37</volume><fpage>9435</fpage><lpage>9444</lpage><pub-id pub-id-type="pmid">9649326</pub-id></mixed-citation></ref><ref id="pone.0013312-Allawi3"><label>45</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Allawi</surname><given-names>HT</given-names></name><name><surname>SantaLucia</surname><given-names>J</given-names></name></person-group><year>1998</year><article-title>Thermodynamics of internal C.T mismatches in DNA.</article-title><source>Nucleic Acids Res</source><volume>26</volume><fpage>2694</fpage><lpage>2701</lpage><pub-id pub-id-type="pmid">9592156</pub-id></mixed-citation></ref><ref id="pone.0013312-Seshadri1"><label>46</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Seshadri</surname><given-names>R</given-names></name><name><surname>Kravitz</surname><given-names>SA</given-names></name><name><surname>Smarr</surname><given-names>L</given-names></name><name><surname>Gilna</surname><given-names>P</given-names></name><name><surname>Frazier</surname><given-names>M</given-names></name></person-group><year>2007</year><article-title>CAMERA: A Community Resource for Metagenomics.</article-title><source>PLoS Biol</source><volume>5</volume><fpage>e75</fpage><pub-id pub-id-type="pmid">17355175</pub-id></mixed-citation></ref><ref id="pone.0013312-DeLong1"><label>47</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>DeLong</surname><given-names>EF</given-names></name><name><surname>Preston</surname><given-names>CM</given-names></name><name><surname>Mincer</surname><given-names>T</given-names></name><name><surname>Rich</surname><given-names>V</given-names></name><name><surname>Hallam</surname><given-names>SJ</given-names></name><etal></etal></person-group><year>2006</year><article-title>Community genomics among stratified microbial assemblages in the ocean's interior.</article-title><source>Science</source><volume>311</volume><fpage>496</fpage><lpage>503</lpage><pub-id pub-id-type="pmid">16439655</pub-id></mixed-citation></ref><ref id="pone.0013312-Rusch1"><label>48</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rusch</surname><given-names>DB</given-names></name><name><surname>Halpern</surname><given-names>AL</given-names></name><name><surname>Sutton</surname><given-names>G</given-names></name><name><surname>Heidelberg</surname><given-names>KB</given-names></name><name><surname>Williamson</surname><given-names>S</given-names></name><etal></etal></person-group><year>2007</year><article-title>The Sorcerer II Global Ocean Sampling Expedition: Northwest Atlantic through Eastern Tropical Pacific.</article-title><source>PLoS Biol</source><volume>5</volume><fpage>e77</fpage><pub-id pub-id-type="pmid">17355176</pub-id></mixed-citation></ref><ref id="pone.0013312-Kan1"><label>49</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kan</surname><given-names>J</given-names></name><name><surname>Suzuki</surname><given-names>MT</given-names></name><name><surname>Wang</surname><given-names>K</given-names></name><name><surname>Evans</surname><given-names>SE</given-names></name><name><surname>Chen</surname><given-names>F</given-names></name></person-group><year>2007</year><article-title>High temporal but low spatial heterogeneity of bacterioplankton in the Chesapeake Bay.</article-title><source>Appl Environ Microbiol</source><volume>73</volume><fpage>6776</fpage><lpage>6789</lpage><pub-id pub-id-type="pmid">17827310</pub-id></mixed-citation></ref><ref id="pone.0013312-Glibert1"><label>50</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Glibert</surname><given-names>PM</given-names></name><name><surname>Garside</surname><given-names>C</given-names></name><name><surname>Fuhrman</surname><given-names>JA</given-names></name><name><surname>Roman</surname><given-names>MR</given-names></name></person-group><year>1991</year><article-title>Time-dependent coupling of inorganic and organic nitrogen uptake and regeneration in the plume of the Chesapeake Bay estuary and its regulation by large heterotrophs.</article-title><source>Limnol Oceanogr</source><volume>36</volume><fpage>895</fpage><lpage>909</lpage></mixed-citation></ref><ref id="pone.0013312-Tiedje1"><label>51</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Tiedje</surname><given-names>J</given-names></name></person-group><year>1988</year><article-title>Ecology of denitrification and dissimilatory nitrate reduction to ammonium.</article-title><person-group person-group-type="editor"><name><surname>Zehnder</surname><given-names>A</given-names></name></person-group><source>Biology of anaerobic microorganisms</source><publisher-loc>New York</publisher-loc><publisher-name>John Wiley & Sons</publisher-name><fpage>179</fpage><lpage>244</lpage></mixed-citation></ref><ref id="pone.0013312-Robertson1"><label>52</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>LA</given-names></name><name><surname>Dalsgaar</surname><given-names>T</given-names></name><name><surname>Revsbech</surname><given-names>NP</given-names></name><name><surname>Kuenen</surname><given-names>JG</given-names></name></person-group><year>1995</year><article-title>Confirmation of ‘aerobic denitrification’ in batch cultures using gas chromatography and 15N mass spectrometry.</article-title><source>FEMS Microbiol Ecol</source><volume>18</volume><fpage>113</fpage><lpage>120</lpage></mixed-citation></ref><ref id="pone.0013312-Robertson2"><label>53</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>LA</given-names></name><name><surname>Van Niel</surname><given-names>EWJ</given-names></name><name><surname>Torremans</surname><given-names>RAM</given-names></name><name><surname>Kuenen</surname><given-names>JG</given-names></name></person-group><year>1988</year><article-title>Simultaneous nitrification and denitrification in aerobic chemostat cultures of Thiosphaera pantotropha.</article-title><source>Appl Environ Microbiol</source><volume>54</volume><fpage>2812</fpage><lpage>2818</lpage><pub-id pub-id-type="pmid">16347780</pub-id></mixed-citation></ref><ref id="pone.0013312-Takaya1"><label>54</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Takaya</surname><given-names>N</given-names></name><name><surname>Catalan-Sakairi</surname><given-names>MAB</given-names></name><name><surname>Sakaguchi</surname><given-names>Y</given-names></name><name><surname>Kato</surname><given-names>I</given-names></name><name><surname>Zhou</surname><given-names>Z</given-names></name><etal></etal></person-group><year>2003</year><article-title>Aerobic denitrifying bacteria that produce low levels of nitrous oxide.</article-title><source>Appl Environ Microbiol</source><volume>69</volume><fpage>3152</fpage><lpage>3157</lpage><pub-id pub-id-type="pmid">12788710</pub-id></mixed-citation></ref><ref id="pone.0013312-Michotey1"><label>55</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Michotey</surname><given-names>V</given-names></name><name><surname>Bonin</surname><given-names>P</given-names></name></person-group><year>1997</year><article-title>Evidence for anaerobic bacterial processes in the water column: denitrification and dissimilatory nitrate ammonification in the northwestern Mediterranean Sea.</article-title><source>Marine Ecology Progress Series</source><volume>160</volume><fpage>47</fpage><lpage>56</lpage></mixed-citation></ref><ref id="pone.0013312-Omnes1"><label>56</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Omnes</surname><given-names>P</given-names></name><name><surname>Slawyk</surname><given-names>G</given-names></name><name><surname>Garcia</surname><given-names>N</given-names></name><name><surname>Bonin</surname><given-names>P</given-names></name></person-group><year>1996</year><article-title>Evidence of denitrification and nitrate ammonification in the River Rhone plume (northwestern Mediterranean Sea).</article-title><source>Marine Ecology Progress Series</source><volume>141</volume><fpage>275</fpage><lpage>281</lpage></mixed-citation></ref><ref id="pone.0013312-John1"><label>57</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>John</surname><given-names>DE</given-names></name><name><surname>Wang</surname><given-names>ZA</given-names></name><name><surname>Liu</surname><given-names>X</given-names></name><name><surname>Byrne</surname><given-names>RH</given-names></name><name><surname>Corredor</surname><given-names>JE</given-names></name><etal></etal></person-group><year>2007</year><article-title>Phytoplankton carbon fixation gene (RuBisCO) transcripts and air-sea CO2 flux in the Mississippi River plume.</article-title><source>ISME J</source><volume>1</volume><fpage>517</fpage><lpage>531</lpage><pub-id pub-id-type="pmid">18043653</pub-id></mixed-citation></ref><ref id="pone.0013312-Treusch1"><label>58</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Treusch</surname><given-names>AH</given-names></name><name><surname>Vergin</surname><given-names>KL</given-names></name><name><surname>Finlay</surname><given-names>LA</given-names></name><name><surname>Donatz</surname><given-names>MG</given-names></name><name><surname>Burton</surname><given-names>RM</given-names></name><etal></etal></person-group><year>2009</year><article-title>Seasonality and vertical structure of microbial communities in an ocean gyre.</article-title><source>ISME J</source><volume>3</volume><fpage>1148</fpage><lpage>1163</lpage><pub-id pub-id-type="pmid">19494846</pub-id></mixed-citation></ref><ref id="pone.0013312-Dubey1"><label>59</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dubey</surname><given-names>AK</given-names></name><name><surname>Baker</surname><given-names>CS</given-names></name><name><surname>Suzuki</surname><given-names>K</given-names></name><name><surname>Jones</surname><given-names>AD</given-names></name><name><surname>Pandit</surname><given-names>P</given-names></name><etal></etal></person-group><year>2003</year><article-title>CsrA regulates translation of the Escherichia coli carbon starvation gene, cstA, by blocking ribosome access to the cstA transcript.</article-title><source>J Bacteriol</source><volume>185</volume><fpage>4450</fpage><lpage>4460</lpage><pub-id pub-id-type="pmid">12867454</pub-id></mixed-citation></ref><ref id="pone.0013312-Schultz1"><label>60</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schultz</surname><given-names>JE</given-names></name><name><surname>Matin</surname><given-names>A</given-names></name></person-group><year>1991</year><article-title>Molecular and functional characterization of a carbon starvation gene of Escherichia coli.</article-title><source>J Mol Biol</source><volume>218</volume><fpage>129</fpage><lpage>140</lpage><pub-id pub-id-type="pmid">1848300</pub-id></mixed-citation></ref><ref id="pone.0013312-Bronk1"><label>61</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Bronk</surname><given-names>DA</given-names></name></person-group><year>2002</year><article-title>Dynamics of DON.</article-title><person-group person-group-type="editor"><name><surname>Hansell</surname><given-names>DA</given-names></name><name><surname>Carlson</surname><given-names>CA</given-names></name></person-group><source>Biogeochemistry of marine dissolved organic matter</source><publisher-loc>San Diego</publisher-loc><publisher-name>Academic Press</publisher-name><fpage>153</fpage><lpage>231</lpage></mixed-citation></ref><ref id="pone.0013312-Ghindilis2"><label>62</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghindilis</surname><given-names>AL</given-names></name><name><surname>Smith</surname><given-names>MW</given-names></name><name><surname>Schwarzkopf</surname><given-names>KR</given-names></name><name><surname>Zhan</surname><given-names>C</given-names></name><name><surname>Evans</surname><given-names>DR</given-names></name><etal></etal></person-group><year>2009</year><article-title>Sensor array: impedimetric label-free sensing of DNA hybridization in real time for rapid, PCR-based detection of microorganisms.</article-title><source>Electroanalysis</source><volume>21</volume><fpage>1459</fpage><lpage>1468</lpage></mixed-citation></ref><ref id="pone.0013312-Preston1"><label>63</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Preston</surname><given-names>CM</given-names></name><name><surname>Marin</surname><given-names>R</given-names><suffix>3rd</suffix></name><name><surname>Jensen</surname><given-names>SD</given-names></name><name><surname>Feldman</surname><given-names>J</given-names></name><name><surname>Birch</surname><given-names>JM</given-names></name><etal></etal></person-group><year>2009</year><article-title>Near real-time, autonomous detection of marine bacterioplankton on a coastal mooring in Monterey Bay, California, using rRNA-targeted DNA probes.</article-title><source>Environ Microbiol</source><volume>11</volume><fpage>1168</fpage><lpage>1180</lpage><pub-id pub-id-type="pmid">19210704</pub-id></mixed-citation></ref><ref id="pone.0013312-Parsons1"><label>64</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Parsons</surname><given-names>T</given-names></name><name><surname>Maita</surname><given-names>Y</given-names></name><name><surname>Lalli</surname><given-names>CM</given-names></name></person-group><year>1984</year><article-title>A Manual of chemical and biological methods for seawater analysis.</article-title><publisher-loc>New York</publisher-loc><publisher-name>Pergamon Press</publisher-name><size units="page"></size></mixed-citation></ref><ref id="pone.0013312-Sugimura1"><label>65</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sugimura</surname><given-names>Y</given-names></name><name><surname>Suzuki</surname><given-names>Y</given-names></name></person-group><year>1988</year><article-title>A high-temperature catalytic oxidation method for the determination of non-volatile dissolved organic carbon in seawater by direct injection of a liquid sample.</article-title><source>Marine Chemistry</source><volume>24</volume><fpage>105</fpage><lpage>131</lpage></mixed-citation></ref><ref id="pone.0013312-Verardo1"><label>66</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Verardo</surname><given-names>DJ</given-names></name><name><surname>Froelich</surname><given-names>PN</given-names></name><name><surname>McIntyre</surname><given-names>A</given-names></name></person-group><year>1990</year><article-title>Determination of organic carbon and nitrogen in marine sediments using the Carlo Erba NA-1500 analyzer. Deep Sea Research Part A.</article-title><source>Oceanographic Research Papers</source><volume>37</volume><fpage>157</fpage><lpage>165</lpage></mixed-citation></ref><ref id="pone.0013312-Quackenbush1"><label>67</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Quackenbush</surname><given-names>J</given-names></name></person-group><year>2001</year><article-title>Computational analysis of microarray data.</article-title><source>Nat Rev Genet</source><volume>2</volume><fpage>418</fpage><lpage>427</lpage><pub-id pub-id-type="pmid">11389458</pub-id></mixed-citation></ref><ref id="pone.0013312-Schena1"><label>68</label><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Schena</surname><given-names>M</given-names></name></person-group><year>2003</year><article-title>Microarray analysis.</article-title><publisher-loc>Hoboken, NJ</publisher-loc><publisher-name>John Wiley & Sons</publisher-name></mixed-citation></ref><ref id="pone.0013312-Eisen1"><label>69</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eisen</surname><given-names>MB</given-names></name><name><surname>Spellman</surname><given-names>PT</given-names></name><name><surname>Brown</surname><given-names>PO</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name></person-group><year>1998</year><article-title>Cluster analysis and display of genome-wide expression patterns.</article-title><source>Proc Natl Acad Sci U S A</source><volume>95</volume><fpage>14863</fpage><lpage>14868</lpage><pub-id pub-id-type="pmid">9843981</pub-id></mixed-citation></ref><ref id="pone.0013312-Edgar1"><label>70</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname><given-names>R</given-names></name><name><surname>Domrachev</surname><given-names>M</given-names></name><name><surname>Lash</surname><given-names>AE</given-names></name></person-group><year>2002</year><article-title>Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.</article-title><source>Nucleic Acids Res</source><volume>30</volume><fpage>207</fpage><lpage>210</lpage><pub-id pub-id-type="pmid">11752295</pub-id></mixed-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 000637 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000637 | 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:2954162
|texte= Seasonal Changes in Bacterial and Archaeal Gene Expression Patterns across Salinity Gradients in the Columbia River Coastal Margin
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:20967204" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a CyberinfraV1
| This area was generated with Dilib version V0.6.25. |  |