Serveur d'exploration sur la visibilité du Havre

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.
***** Acces problem to record *****\

Identifieur interne : 000142 ( Pmc/Corpus ); précédent : 0001419; suivant : 0001430 ***** probable Xml problem with record *****

Links to Exploration step


Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Data processing and classification analysis of proteomic changes: a case study of oil pollution in the mussel,
<italic>Mytilus edulis</italic>
</title>
<author>
<name sortKey="Monsinjon, Tiphaine" sort="Monsinjon, Tiphaine" uniqKey="Monsinjon T" first="Tiphaine" last="Monsinjon">Tiphaine Monsinjon</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Andersen, Odd Ketil" sort="Andersen, Odd Ketil" uniqKey="Andersen O" first="Odd Ketil" last="Andersen">Odd Ketil Andersen</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leboulenger, Francois" sort="Leboulenger, Francois" uniqKey="Leboulenger F" first="François" last="Leboulenger">François Leboulenger</name>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Knigge, Thomas" sort="Knigge, Thomas" uniqKey="Knigge T" first="Thomas" last="Knigge">Thomas Knigge</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">16970821</idno>
<idno type="pmc">1592071</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1592071</idno>
<idno type="RBID">PMC:1592071</idno>
<idno type="doi">10.1186/1477-5956-4-17</idno>
<date when="2006">2006</date>
<idno type="wicri:Area/Pmc/Corpus">000142</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Data processing and classification analysis of proteomic changes: a case study of oil pollution in the mussel,
<italic>Mytilus edulis</italic>
</title>
<author>
<name sortKey="Monsinjon, Tiphaine" sort="Monsinjon, Tiphaine" uniqKey="Monsinjon T" first="Tiphaine" last="Monsinjon">Tiphaine Monsinjon</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Andersen, Odd Ketil" sort="Andersen, Odd Ketil" uniqKey="Andersen O" first="Odd Ketil" last="Andersen">Odd Ketil Andersen</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leboulenger, Francois" sort="Leboulenger, Francois" uniqKey="Leboulenger F" first="François" last="Leboulenger">François Leboulenger</name>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Knigge, Thomas" sort="Knigge, Thomas" uniqKey="Knigge T" first="Thomas" last="Knigge">Thomas Knigge</name>
<affiliation>
<nlm:aff id="I1">IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="I2">Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">Proteome Science</title>
<idno type="eISSN">1477-5956</idno>
<imprint>
<date when="2006">2006</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<sec>
<title>Background</title>
<p>Proteomics may help to detect subtle pollution-related changes, such as responses to mixture pollution at low concentrations, where clear signs of toxicity are absent. The challenges associated with the analysis of large-scale multivariate proteomic datasets have been widely discussed in medical research and biomarker discovery. This concept has been introduced to ecotoxicology only recently, so data processing and classification analysis need to be refined before they can be readily applied in biomarker discovery and monitoring studies.</p>
</sec>
<sec>
<title>Results</title>
<p>Data sets obtained from a case study of oil pollution in the Blue mussel were investigated for differential protein expression by retentate chromatography-mass spectrometry and decision tree classification. Different tissues and different settings were used to evaluate classifiers towards their discriminatory power. It was found that, due the intrinsic variability of the data sets, reliable classification of unknown samples could only be achieved on a broad statistical basis (n > 60) with the observed expression changes comprising high statistical significance and sufficient amplitude. The application of stringent criteria to guard against overfitting of the models eventually allowed satisfactory classification for only one of the investigated data sets and settings.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>Machine learning techniques provide a promising approach to process and extract informative expression signatures from high-dimensional mass-spectrometry data. Even though characterisation of the proteins forming the expression signatures would be ideal, knowledge of the specific proteins is not mandatory for effective class discrimination. This may constitute a new biomarker approach in ecotoxicology, where working with organisms, which do not have sequenced genomes render protein identification by database searching problematic. However, data processing has to be critically evaluated and statistical constraints have to be considered before supervised classification algorithms are employed.</p>
</sec>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></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">Proteome Sci</journal-id>
<journal-title>Proteome Science</journal-title>
<issn pub-type="epub">1477-5956</issn>
<publisher>
<publisher-name>BioMed Central</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">16970821</article-id>
<article-id pub-id-type="pmc">1592071</article-id>
<article-id pub-id-type="publisher-id">1477-5956-4-17</article-id>
<article-id pub-id-type="doi">10.1186/1477-5956-4-17</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Methodology</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Data processing and classification analysis of proteomic changes: a case study of oil pollution in the mussel,
<italic>Mytilus edulis</italic>
</article-title>
</title-group>
<contrib-group>
<contrib id="A1" corresp="yes" contrib-type="author">
<name>
<surname>Monsinjon</surname>
<given-names>Tiphaine</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<xref ref-type="aff" rid="I2">2</xref>
<email>Tiphaine.Monsinjon@univ-lehavre.fr</email>
</contrib>
<contrib id="A2" contrib-type="author">
<name>
<surname>Andersen</surname>
<given-names>Odd Ketil</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>oka@irisresearch.no</email>
</contrib>
<contrib id="A3" contrib-type="author">
<name>
<surname>Leboulenger</surname>
<given-names>François</given-names>
</name>
<xref ref-type="aff" rid="I2">2</xref>
<email>Francois.Leboulenger@univ-lehavre.fr</email>
</contrib>
<contrib id="A4" contrib-type="author">
<name>
<surname>Knigge</surname>
<given-names>Thomas</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<xref ref-type="aff" rid="I2">2</xref>
<email>Thomas.Knigge@univ-lehavre.fr</email>
</contrib>
</contrib-group>
<aff id="I1">
<label>1</label>
IRIS – International Research Institute of Stavanger AS, Randaberg, Norway</aff>
<aff id="I2">
<label>2</label>
Laboratoire d'Ecotoxicologie – Milieux Aquatiques, Université du Havre, Le Havre, France</aff>
<pub-date pub-type="collection">
<year>2006</year>
</pub-date>
<pub-date pub-type="epub">
<day>13</day>
<month>9</month>
<year>2006</year>
</pub-date>
<volume>4</volume>
<fpage>17</fpage>
<lpage>17</lpage>
<ext-link ext-link-type="uri" xlink:href="http://www.proteomesci.com/content/4/1/17"></ext-link>
<history>
<date date-type="received">
<day>11</day>
<month>4</month>
<year>2006</year>
</date>
<date date-type="accepted">
<day>13</day>
<month>9</month>
<year>2006</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright © 2006 Monsinjon et al; licensee BioMed Central Ltd.</copyright-statement>
<copyright-year>2006</copyright-year>
<copyright-holder>Monsinjon et al; licensee BioMed Central Ltd.</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0">
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"></ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p>
<pmc-comment> Monsinjon Tiphaine Tiphaine.Monsinjon@univ-lehavre.fr Data processing and classification analysis of proteomic changes: a case study of oil pollution in the mussel, Mytilus edulis 2006Proteome Science 4(1): 17-. (2006)1477-5956(2006)4:1<17>urn:ISSN:1477-5956</pmc-comment>
</license>
</permissions>
<abstract>
<sec>
<title>Background</title>
<p>Proteomics may help to detect subtle pollution-related changes, such as responses to mixture pollution at low concentrations, where clear signs of toxicity are absent. The challenges associated with the analysis of large-scale multivariate proteomic datasets have been widely discussed in medical research and biomarker discovery. This concept has been introduced to ecotoxicology only recently, so data processing and classification analysis need to be refined before they can be readily applied in biomarker discovery and monitoring studies.</p>
</sec>
<sec>
<title>Results</title>
<p>Data sets obtained from a case study of oil pollution in the Blue mussel were investigated for differential protein expression by retentate chromatography-mass spectrometry and decision tree classification. Different tissues and different settings were used to evaluate classifiers towards their discriminatory power. It was found that, due the intrinsic variability of the data sets, reliable classification of unknown samples could only be achieved on a broad statistical basis (n > 60) with the observed expression changes comprising high statistical significance and sufficient amplitude. The application of stringent criteria to guard against overfitting of the models eventually allowed satisfactory classification for only one of the investigated data sets and settings.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>Machine learning techniques provide a promising approach to process and extract informative expression signatures from high-dimensional mass-spectrometry data. Even though characterisation of the proteins forming the expression signatures would be ideal, knowledge of the specific proteins is not mandatory for effective class discrimination. This may constitute a new biomarker approach in ecotoxicology, where working with organisms, which do not have sequenced genomes render protein identification by database searching problematic. However, data processing has to be critically evaluated and statistical constraints have to be considered before supervised classification algorithms are employed.</p>
</sec>
</abstract>
</article-meta>
</front>
<body>
<sec>
<title>Background</title>
<p>In ecotoxicology, proteomics has been increasingly applied for the screening of protein expression changes caused by pollutants [
<xref ref-type="bibr" rid="B1">1</xref>
-
<xref ref-type="bibr" rid="B10">10</xref>
]. Proteomic profiles are generally altered upon stress and hence should also be distinctive of certain types of toxic exposure [
<xref ref-type="bibr" rid="B2">2</xref>
]. 2-DE is still the prevailing tool for protein separation and analysis, followed by MS to identify the proteins of interest. MS-based protein profiling is rapidly emerging within the current applications of proteomic technologies. Especially in medical research, it has received increasing attention for the development of new diagnostic criteria and disease monitoring. Consequently, there exists a large body of literature on how to best exploit datasets from proteomic MS, comprising large sample sizes and high dimensionality. Owing to the amount and complexity of the generated data, machine-learning techniques have been considered the methods of choice for analysing such type of multivariate data [
<xref ref-type="bibr" rid="B11">11</xref>
,
<xref ref-type="bibr" rid="B12">12</xref>
]. A plethora of supervised learning methods, such as partial least squares, discriminant and logistic regression analysis, genetic algorithms, artificial neural networks,
<italic>k</italic>
-nearest-neighbour, support vector machines and decision trees have been evaluated for this purpose [reviewed in [
<xref ref-type="bibr" rid="B11">11</xref>
,
<xref ref-type="bibr" rid="B13">13</xref>
,
<xref ref-type="bibr" rid="B14">14</xref>
]].</p>
<p>In conjunction with SELDI, this concept has also been introduced into ecotoxicology [
<xref ref-type="bibr" rid="B5">5</xref>
]. After Hogstrand
<italic>et al</italic>
. employed SELDI for protein profiling of rainbow trout gills during exposure to waterborne zinc [
<xref ref-type="bibr" rid="B2">2</xref>
], Knigge
<italic>et al</italic>
. have suggested the use of classifiers based on ProteinChip mass-spectra to identify field-exposure to copper or PAH in the bivalve
<italic>Mytilus edulis </italic>
(L.) [
<xref ref-type="bibr" rid="B5">5</xref>
], widely used as an environmental sentinel. Such classifiers may constitute a so-called "biomarker pattern" which can be of high discriminatory power, regardless of the identity of the specific proteins that form it.</p>
<p>The global analysis of cellular constituents potentially provides a more comprehensive view of toxicity, since toxicity generally involves not only single interactions but also triggers a cascade of alterations [
<xref ref-type="bibr" rid="B15">15</xref>
]. The ability to display a multitude of alterations renders it particularly suitable for the evaluation of combined exposure to toxicants [
<xref ref-type="bibr" rid="B16">16</xref>
]. Moreover, proteomics allows molecular fingerprinting of protein expression changes following the exposures to low levels of contaminants, where conventional methods of toxicology may not be sufficiently sensitive [
<xref ref-type="bibr" rid="B17">17</xref>
]. More importantly, perhaps, a set of proteins can potentially achieve higher accuracy and specificity than any single biomarker alone [
<xref ref-type="bibr" rid="B18">18</xref>
].</p>
<p>Oil is a complex mixture of hydrocarbons including PAHs which are strong carcinogens to humans and wildlife [
<xref ref-type="bibr" rid="B19">19</xref>
] and APs which were demonstrated to elicit oestrogen agonistic and/or antagonistic properties, thereby exerting hormone-modulating effects [
<xref ref-type="bibr" rid="B20">20</xref>
]. Even if these chemicals are rapidly diluted in the vicinity of offshore installations just after their discharge, the effluents may still exert low concentration effects [
<xref ref-type="bibr" rid="B21">21</xref>
], giving reason to evaluate whether these discharges may harm the biological resources in the sea.</p>
<p>As datasets from proteomic MS are typically characterised by large numbers of variables but relatively small sample sizes, untypical machine learning problems are encountered and supervised training of classifiers may be considered problematic [
<xref ref-type="bibr" rid="B22">22</xref>
]. In particular, overfitting of the multivariable models may seriously undermine the biological relevance of the sets of proteins constituting the classifier [
<xref ref-type="bibr" rid="B23">23</xref>
,
<xref ref-type="bibr" rid="B24">24</xref>
]. Owing to these concerns, the present study focuses on pre-processing and analysis of proteomic MS data generated with SELDI, thereby taking a refined and more conservative approach then reported previously [
<xref ref-type="bibr" rid="B5">5</xref>
,
<xref ref-type="bibr" rid="B8">8</xref>
,
<xref ref-type="bibr" rid="B9">9</xref>
]. To evaluate the developed protocol for data processing and classifier generation, different datasets derived from a 21-day exposure of Blue mussels in a laboratory flow-through system to crude oil with and without a spike of PAHs and APs were used. We included feature selection based on univariate statistical methods prior to decision tree classification. For final evaluation, a minimum variation threshold for the definite peaks was considered in order to advance biological relevance and robustness of the protein expression signatures [
<xref ref-type="bibr" rid="B25">25</xref>
,
<xref ref-type="bibr" rid="B26">26</xref>
].</p>
</sec>
<sec>
<title>Results</title>
<sec>
<title>SELDI profiling of gill proteins</title>
<p>A panel of gill-spectra (n = 51 for C, n = 66 for oil and n = 55 for sO; note: the discrepancy between the final number of spectra analysed and the number of animals sampled in total results from the removal of low quality spectra prior to analysis; see material and methods section) yielded 55 peaks with differential expression. Statistical analysis revealed 20 of these to be significantly altered, with levels of differential expression below 50%. Just one peak (
<italic>m/z </italic>
4755) displayed approximately two-fold alteration, albeit with rather high variability and hence low significance (0.05 ≥
<italic>p </italic>
> 0.01). Of the 20 peaks showing significantly altered peak intensities, only six were classified highly significant (
<italic>p </italic>
< 0.001). The most prominent peaks, combining high statistical significance with at least 50 % change in expression were
<italic>m/z </italic>
6696 found in oil and sO as well as
<italic>m/z </italic>
9661 and
<italic>m/z </italic>
14 847 in oil only, all of which were down-regulated (Table
<xref ref-type="table" rid="T1">1A</xref>
). Remarkably, the effect of oil exposure on gills tended to be more pronounced than that of sO, as generally more peaks were found at a highly significant level and the expression changes were slightly stronger (data not shown). In addition, eight peaks were unique to oil exposure (data not shown), amongst them the one at
<italic>m/z </italic>
9661.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption>
<p>Differentially expressed peptides/proteins from gills (A) and digestive gland (DG, B) of Blue mussels obtained from controls (C) and 21d exposures to 0.5 mg/L crude oil (oil) and 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO); mean intensities (arbitrary values) and SD of controls (Gills, n = 51; DG, n = 74), oil (Gills, n = 66; DG, n = 71) and sO (Gills, n = 55; DG, n = 69). Expression changes relative to the controls are significantly different at a level of
<italic>p </italic>
< 0.001 (one-way ANOVA for C
<italic>vs</italic>
. oil
<italic>vs</italic>
. sO) at least for one of the exposures. Peaks are listed in decreasing order of their respective
<italic>p </italic>
values from top to bottom.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td align="center">Peak [m/z]</td>
<td align="center" colspan="3">Mean intensity</td>
</tr>
<tr>
<td></td>
<td align="center">
<bold>C</bold>
</td>
<td align="center">
<bold>oil</bold>
</td>
<td align="center">
<bold>sO</bold>
</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" colspan="4">
<bold>A –
<italic>Gills</italic>
</bold>
</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="center" colspan="4">Down-regulation</td>
</tr>
<tr>
<td align="center">14847</td>
<td align="center">
<bold>1.01 </bold>
±
<italic>0.32</italic>
</td>
<td align="center">
<bold>0.61 </bold>
±
<italic>0.26</italic>
</td>
<td align="center">
<bold>0.81 </bold>
±
<italic>0.31</italic>
</td>
</tr>
<tr>
<td align="center">6696</td>
<td align="center">
<bold>2.76 </bold>
±
<italic>1.19</italic>
</td>
<td align="center">
<bold>1.86 </bold>
±
<italic>1.34</italic>
</td>
<td align="center">
<bold>1.64 </bold>
±
<italic>1.33</italic>
</td>
</tr>
<tr>
<td align="center">9661</td>
<td align="center">
<bold>8.11 </bold>
±
<italic>3.82</italic>
</td>
<td align="center">
<bold>5.35 </bold>
±
<italic>3.63</italic>
</td>
<td align="center">
<bold>8.35 </bold>
±
<italic>5.06</italic>
</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="left" colspan="4">
<bold>B –
<italic>DG</italic>
</bold>
</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="center" colspan="4">Up-regulation</td>
</tr>
<tr>
<td align="center">7811</td>
<td align="center">
<bold>0.83 </bold>
±
<italic>0.74</italic>
</td>
<td align="center">
<bold>1.46 </bold>
±
<italic>1.22</italic>
</td>
<td align="center">
<bold>2.76 </bold>
±
<italic>1.84</italic>
</td>
</tr>
<tr>
<td align="center">16196</td>
<td align="center">
<bold>0.75 </bold>
±
<italic>1.20</italic>
</td>
<td align="center">
<bold>1.23 </bold>
±
<italic>1.46</italic>
</td>
<td align="center">
<bold>2.18 </bold>
±
<italic>2.06</italic>
</td>
</tr>
<tr>
<td align="center">4744</td>
<td align="center">
<bold>0.43 </bold>
±
<italic>1.05</italic>
</td>
<td align="center">
<bold>0.65 </bold>
±
<italic>1.60</italic>
</td>
<td align="center">
<bold>1.37 </bold>
±
<italic>1.80</italic>
</td>
</tr>
<tr>
<td align="center">3795</td>
<td align="center">
<bold>0.84 </bold>
±
<italic>1.42</italic>
</td>
<td align="center">
<bold>1.33 </bold>
±
<italic>2.19</italic>
</td>
<td align="center">
<bold>2.64 </bold>
±
<italic>3.39</italic>
</td>
</tr>
<tr>
<td align="center">3984</td>
<td align="center">
<bold>1.04 </bold>
±
<italic>2.03</italic>
</td>
<td align="center">
<bold>1.04 </bold>
±
<italic>1.97</italic>
</td>
<td align="center">
<bold>2.67 </bold>
±
<italic>3.05</italic>
</td>
</tr>
<tr>
<td align="center">4108</td>
<td align="center">
<bold>0.49 </bold>
±
<italic>1.36</italic>
</td>
<td align="center">
<bold>0.44 </bold>
±
<italic>1.32</italic>
</td>
<td align="center">
<bold>1.34 </bold>
±
<italic>1.95</italic>
</td>
</tr>
<tr>
<td align="center">3552</td>
<td align="center">
<bold>0.69 </bold>
±
<italic>2.00</italic>
</td>
<td align="center">
<bold>0.54 </bold>
±
<italic>1.34</italic>
</td>
<td align="center">
<bold>1.67 </bold>
±
<italic>2.66</italic>
</td>
</tr>
<tr>
<td align="center">28764</td>
<td align="center">
<bold>0.76 </bold>
±
<italic>0.7</italic>
7</td>
<td align="center">
<bold>1.78 </bold>
±
<italic>1.45</italic>
</td>
<td align="center">
<bold>1.76 </bold>
±
<italic>1.23</italic>
</td>
</tr>
<tr>
<td align="center">52306</td>
<td align="center">
<bold>0.34 </bold>
±
<italic>0.25</italic>
</td>
<td align="center">
<bold>0.72 </bold>
±
<italic>0.73</italic>
</td>
<td align="center">
<bold>0.46 </bold>
±
<italic>0.44</italic>
</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="center" colspan="4">Down-regulation</td>
</tr>
<tr>
<td align="center">9172</td>
<td align="center">
<bold>3.78 </bold>
±
<italic>4.26</italic>
</td>
<td align="center">
<bold>2.46 </bold>
±
<italic>2.21</italic>
</td>
<td align="center">
<bold>1.01 </bold>
±
<italic>0.86</italic>
</td>
</tr>
<tr>
<td align="center">18250</td>
<td align="center">
<bold>0.47 </bold>
±
<italic>0.40</italic>
</td>
<td align="center">
<bold>0.23 </bold>
±
<italic>0.26</italic>
</td>
<td align="center">
<bold>0.15 </bold>
±
<italic>0.17</italic>
</td>
</tr>
<tr>
<td align="center">7410</td>
<td align="center">
<bold>1.09 </bold>
±
<italic>0.59</italic>
</td>
<td align="center">
<bold>0.74 </bold>
±
<italic>0.57</italic>
</td>
<td align="center">
<bold>0.55 </bold>
±
<italic>0.32</italic>
</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec>
<title>SELDI profiling of digestive gland proteins</title>
<p>From the panel of DG-spectra (n = 74 for C, n = 71 for oil and n = 69 for sO), 88 peaks were found to be differentially expressed. Of the 49 peaks showing significantly altered peak intensities, more than 60 % were highly significant at a level of
<italic>p </italic>
< 0.001 and the expression levels of almost 30 % were altered more than two-fold, the majority of which were up-regulated (Table
<xref ref-type="table" rid="T1">1B</xref>
). In terms of significance and level of expression change, the peaks at
<italic>m/z </italic>
7811 and
<italic>m/z </italic>
9172 should be highlighted, especially for their response to additional PAHs and APs.</p>
<p>In general, increases or decreases of expression detected in both types of exposure behaved consistently, commonly being less pronounced in the oil-exposure. Evidently, spiking crude oil with APs and PAHs augmented responses compared to oil alone in a somewhat 'dose dependent' manner (Figure
<xref ref-type="fig" rid="F1">1</xref>
). Moreover, some peaks appeared to be altered solely in sO-exposure (dotted bars in Figure
<xref ref-type="fig" rid="F1">1</xref>
). Spike-related intensification of differential protein expression was also characterised by higher levels of significance for certain sO-peaks (data not shown).</p>
<fig position="float" id="F1">
<label>Figure 1</label>
<caption>
<p>Relative expression changes of differentially expressed peptides and proteins in the digestive gland of Blue mussels following 21d exposure to 0.5 mg/L crude oil (oil) and 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO) in % of control intensity; *, significantly different at 0.05 ≥
<italic>p </italic>
≥ 0.01, **, at 0.01 ≥
<italic>p </italic>
≥ 0.001, ***, at
<italic>p </italic>
< 0.001. Dotted bars indicate differential expression occurring in one of two exposures only.</p>
</caption>
<graphic xlink:href="1477-5956-4-17-1"></graphic>
</fig>
</sec>
<sec>
<title>SELDI profiling of sex-dependent responses</title>
<p>To obtain information on differential expression with respect to the gender of the mussels, the datasets were split into six groups, a male and a female one for the controls and the two exposures respectively. Consequently, in the face of high inter-individual variability, the statistics were weakened. This was reflected in fewer percentages of highly significant peaks (
<italic>p </italic>
< 0.001; gills: 14 % for females, 21 % for males; DG: 27 % for females, 42 % for males). The majority of these peaks were already displayed in the overall analysis in which gender information was not included. For example, the most prominent peaks for DG,
<italic>m/z </italic>
7811 and
<italic>m/z </italic>
9172 are found to be highly significant in males and in females with a strong degree of up- or down-regulation. However, by splitting the data, additional information on these peaks could be acquired, as for instance, the peak at
<italic>m/z </italic>
7811 obviously responded more strongly in females than in males: in oil-female, it shows a two-fold increase (
<italic>p </italic>
< 0.001) compared to a non-significant 1.5 fold in the oil-male as well as a 4.5 fold increase in sO-female compared to 2.5 fold in sO-male, both at
<italic>p </italic>
< 0.001. Similarly, the peak at
<italic>m/z </italic>
18 250 shows a significant 3 fold down-regulation for the females of the oil exposure plus a highly significant 5 fold down-regulation for the sO-females but this is less pronounced in males with a 1.5 fold decrease for oil (
<italic>p </italic>
≤ 0.05) and 2.5 fold decrease for sO (
<italic>p </italic>
< 0.001). Accordingly, these incidences characterise the sex-differentiated variations of a general response.</p>
<p>To extract truly male- or female-specific responses, the univariate statistics carried out for males and females separately have been 'subtracted' from each other, in order to keep only peaks appearing in one of the genders. In DG, where most of the sex-specific alterations were obtained, 22 peaks showing differential expression were unique to females and 20 to males. Regarding the levels of differential expression, 55 % of these male-specific peaks were altered more than two-fold compared to only 36 % in females. Eventually, the sex-specific peaks were crosschecked for similar but non-significant differential expression, and if such occurred, those responses were regarded less likely to be genuinely sex-specific and thus the peaks were excluded. Figure
<xref ref-type="fig" rid="F2">2</xref>
shows the most definite peaks for DG, which were attributed to a clear sex-specific response; notably they were all found to be in the peptide range. Another five peptides represented some kind of sex-differentiated response in which up-regulation was found in one gender and down-regulation in the other, four of which were significantly altered only in males. The majority of these alterations were triggered by spiking the crude oil with additional APs and PAHs.</p>
<fig position="float" id="F2">
<label>Figure 2</label>
<caption>
<p>Histogram of sex dependent expression changes (intensities as arbitrary units) obtained for digestive gland of
<italic>M. edulis </italic>
upon 21d exposure to 0.5 mg/L crude oil (oil) and 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO); CF = control females, CM = control males, sOF = spiked oil females, sOM = spiked oil males, oilF = oil females, oilM = oil females; males are in blue and females in orange; the upper five represent female-specific and lower three male-specific peaks. Significantly different at 0.05 ≥
<italic>p </italic>
≥ 0.01, **, at 0.01 ≥
<italic>p </italic>
≥ 0.001, ***, at
<italic>p </italic>
< 0.001.</p>
</caption>
<graphic xlink:href="1477-5956-4-17-2"></graphic>
</fig>
<p>In addition, the datasets for gills were examined for sex-specific responses. Although this tissue was expected to be less prone to respond in a sex-specific manner, some proteins were revealed; the most important ones are depicted in Figure
<xref ref-type="fig" rid="F3">3</xref>
. These proteins were predominantly found in males, the majority of which exceeded 10 kDa. There was but one remarkable exception: more detailed analysis of the data revealed a peptide of
<italic>m/z </italic>
4185 that was practically absent in the controls and induced in oil (~17-fold) and sO (~10-fold), although only significant in the latter at a level of
<italic>p </italic>
≤ 0.05 due to an enormous variability as indicated by a high SD.</p>
<fig position="float" id="F3">
<label>Figure 3</label>
<caption>
<p>Histogram of sex dependent expression changes (intensities as arbitrary units) obtained for gills of
<italic>M. edulis </italic>
upon 21d exposure to 0.5 mg/L crude oil (oil) and 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO); CF = control females, CM = control males, sOF = spiked oil females, sOM = spiked oil males, oilF = oil females, oilM = oil females; males are depicted as blue and females as orange; the upper two represent female-specific and lower four male-specific peaks. Significantly different at 0.05 ≥
<italic>p </italic>
≥ 0.01, **, at 0.01 ≥
<italic>p </italic>
≥ 0.001, ***, at
<italic>p </italic>
< 0.001.</p>
</caption>
<graphic xlink:href="1477-5956-4-17-3"></graphic>
</fig>
</sec>
<sec>
<title>Classifiers</title>
<p>The above data analysis was complemented by generating classification models using pattern recognition software based on CART to validate the predictive value of the differentially expressed proteins. In principal, the models on the gill data required more nodes in the decision tree associated with higher error costs to carry out the classification to a satisfactory degree than the ones for the DG data. This is likely to diminish their predictive reliability and accounts for the absence of robustness as was confirmed by comparing the estimated classification success on the basis of random cross-validation of 10% of the LS with the actual classification success as validated with the TS. The two best models for gills predicted a relatively good classification with 73 % of the oil-samples correctly attributed, 81 % of the sO-samples and 88 % for the controls each. Only the sO-model for DG assumed to perform a better classification of 91 %. Still, the overall classification success for gills, when tested with unknown samples, turned out not to be as estimated and less good than for DG (Table
<xref ref-type="table" rid="T2">2</xref>
). For example, with the gill-TS, merely 50 % of the sO-samples could be assigned correctly and for oil-exposure, only 60 % of the controls were recognised (Table
<xref ref-type="table" rid="T2">2A</xref>
). In other words, both models generated for gills showed sufficient classification success only for one of two groups. On the contrary, the models for DG yielded relatively good classification success for the controls and sO-samples; just the oil samples showed less than 70 % correct classification (Table
<xref ref-type="table" rid="T2">2B</xref>
). However, these models also displayed significant discrepancy between estimated and actual classification successes (maximum 11 percentage points), although not in the same extent as for the gills (maximum 31 percentage points).</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption>
<p>Prediction success of the classifiers for gills (A) and digestive gland (B) from the Blue mussel for pairwise comparison of controls (C) with either of the two different exposures to i) 0.5 mg/L crude oil (oil) and ii) 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO) for 21 days as validated with an independent sample set not used in model construction. Sample size for learning and testing set, amount of samples attributed to each class, individual and overall classification success in % of testing set is presented.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td align="center" colspan="8">
<bold>Prediction success</bold>
</td>
</tr>
<tr>
<td align="left" colspan="8">
<bold>A </bold>
<italic>Gills</italic>
</td>
</tr>
</thead>
<tbody>
<tr>
<td align="center">
<bold>Class</bold>
</td>
<td align="center">
<bold>Learning</bold>
</td>
<td align="center">
<bold>Testing</bold>
</td>
<td align="center">
<bold>%</bold>
</td>
<td align="center">
<bold>C</bold>
</td>
<td align="center">
<bold>oil</bold>
</td>
<td align="center">
<bold>sO</bold>
</td>
<td align="center">
<bold>Tot %</bold>
</td>
</tr>
<tr>
<td colspan="8">
<hr></hr>
</td>
</tr>
<tr>
<td align="center">
<bold>C</bold>
</td>
<td align="center">34</td>
<td align="center">17</td>
<td align="center">65</td>
<td align="center">11</td>
<td align="center">6</td>
<td></td>
<td align="center">72</td>
</tr>
<tr>
<td align="center">
<bold>oil</bold>
</td>
<td align="center">44</td>
<td align="center">22</td>
<td align="center">77</td>
<td align="center">5</td>
<td align="center">17</td>
<td></td>
<td></td>
</tr>
<tr>
<td align="center">
<bold>C</bold>
</td>
<td align="center">34</td>
<td align="center">17</td>
<td align="center">82</td>
<td align="center">14</td>
<td></td>
<td align="center">3</td>
<td align="center">66</td>
</tr>
<tr>
<td align="center">
<bold>sO</bold>
</td>
<td align="center">37</td>
<td align="center">18</td>
<td align="center">50</td>
<td align="center">9</td>
<td></td>
<td align="center">9</td>
<td></td>
</tr>
<tr>
<td colspan="8">
<hr></hr>
</td>
</tr>
<tr>
<td align="left" colspan="8">
<bold>B </bold>
<italic>DG</italic>
</td>
</tr>
<tr>
<td colspan="8">
<hr></hr>
</td>
</tr>
<tr>
<td align="center">
<bold>Class</bold>
</td>
<td align="center">
<bold>Learning</bold>
</td>
<td align="center">
<bold>Testing</bold>
</td>
<td align="center">
<bold>%</bold>
</td>
<td align="center">
<bold>C</bold>
</td>
<td align="center">
<bold>oil</bold>
</td>
<td align="center">
<bold>sO</bold>
</td>
<td align="center">
<bold>Tot %</bold>
</td>
</tr>
<tr>
<td colspan="8">
<hr></hr>
</td>
</tr>
<tr>
<td align="center">
<bold>C</bold>
</td>
<td align="center">48</td>
<td align="center">26</td>
<td align="center">85</td>
<td align="center">22</td>
<td align="center">4</td>
<td></td>
<td align="center">76</td>
</tr>
<tr>
<td align="center">
<bold>oil</bold>
</td>
<td align="center">47</td>
<td align="center">24</td>
<td align="center">67</td>
<td align="center">7</td>
<td align="center">17</td>
<td></td>
<td></td>
</tr>
<tr>
<td align="center">
<bold>C</bold>
</td>
<td align="center">48</td>
<td align="center">26</td>
<td align="center">77</td>
<td align="center">20</td>
<td></td>
<td align="center">6</td>
<td align="center">80</td>
</tr>
<tr>
<td align="center">
<bold>sO</bold>
</td>
<td align="center">45</td>
<td align="center">24</td>
<td align="center">83</td>
<td align="center">4</td>
<td></td>
<td align="center">20</td>
<td></td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The models for gills and for DG did well reflect the significance of certain peaks as ranked by univariate statistics and by their levels of up- or down-regulation (Table
<xref ref-type="table" rid="T1">1</xref>
and Table
<xref ref-type="table" rid="T3">3</xref>
). Since the models for gills and for DG were generated for two groups only (C
<italic>vs</italic>
. oil plus C
<italic>vs</italic>
. sO), the most important peptides and proteins as listed in Table
<xref ref-type="table" rid="T1">1</xref>
were not necessarily represented in both models. Whenever this was the case, they did not exactly display the same discriminatory power for classification of oil and sO samples respectively (Table
<xref ref-type="table" rid="T3">3A</xref>
and
<xref ref-type="table" rid="T3">3B</xref>
). Attempts to generate classifiers with three groups including controls and both of the exposures failed. Those classifiers were unable to effectively distinguish controls from exposures and thus were omitted. If taken together, though, most of the prominent peaks determined by Kruskal-Wallis One-way ANOVA can be retrieved in either of the two classifiers for oil and sO.</p>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption>
<p>Peak constituents of the optimal classification models for gills (A) and digestive gland (B) of Blue mussels exposed to 0.5 mg/L crude oil (oil) and 0.5 mg/L crude oil spiked with a mixture of APs and PAHs (0.1 mg/L; sO) for 21 days. The score equals the discriminatory weight of the variable within the classifier. A pre-selection of highly significant peaks (
<italic>p </italic>
< 0.001) has been carried out prior to model construction.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td align="left" colspan="4">
<bold>A </bold>
<italic>Gills</italic>
</td>
</tr>
</thead>
<tbody>
<tr>
<td align="center" colspan="2">oil
<italic>vs</italic>
. C</td>
<td align="center" colspan="2">sO
<italic>vs</italic>
. C</td>
</tr>
<tr>
<td align="center">Peak [m/z]</td>
<td align="center">Score [%]</td>
<td align="center">Peak [m/z]</td>
<td align="center">Score [%]</td>
</tr>
<tr>
<td colspan="2">
<hr></hr>
</td>
<td colspan="2">
<hr></hr>
</td>
</tr>
<tr>
<td align="center">14847</td>
<td align="center">100.00</td>
<td align="center">6696</td>
<td align="center">100.00</td>
</tr>
<tr>
<td align="center">9661</td>
<td align="center">86.57</td>
<td align="center">28017</td>
<td align="center">45.98</td>
</tr>
<tr>
<td align="center">6696</td>
<td align="center">65.34</td>
<td align="center">40456</td>
<td align="center">35.26</td>
</tr>
<tr>
<td align="center">44667</td>
<td align="center">54.41</td>
<td align="center">14847</td>
<td align="center">27.95</td>
</tr>
<tr>
<td align="center">28017</td>
<td align="center">52.31</td>
<td align="center">110694</td>
<td align="center">19.99</td>
</tr>
<tr>
<td align="center">40456</td>
<td align="center">31.73</td>
<td align="center">4755</td>
<td align="center">18.53</td>
</tr>
<tr>
<td align="center">12848</td>
<td align="center">12.86</td>
<td align="center">31830</td>
<td align="center">17.12</td>
</tr>
<tr>
<td align="center">20316</td>
<td align="center">10.60</td>
<td align="center">6103</td>
<td align="center">16.51</td>
</tr>
<tr>
<td align="center">29340</td>
<td align="center">9.09</td>
<td align="center">88989</td>
<td align="center">14.53</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">44667</td>
<td align="center">12.89</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="left" colspan="4">
<bold>B </bold>
<italic>DG</italic>
</td>
</tr>
<tr>
<td colspan="4">
<hr></hr>
</td>
</tr>
<tr>
<td align="center" colspan="2">oil
<italic>vs</italic>
. C</td>
<td align="center" colspan="2">sO
<italic>vs</italic>
. C</td>
</tr>
<tr>
<td align="center">Peak [m/z]</td>
<td align="center">Score [%]</td>
<td align="center">Peak [m/z]</td>
<td align="center">Score [%]</td>
</tr>
<tr>
<td colspan="2">
<hr></hr>
</td>
<td colspan="2">
<hr></hr>
</td>
</tr>
<tr>
<td align="center">16413</td>
<td align="center">100.00</td>
<td align="center">7811</td>
<td align="center">100.00</td>
</tr>
<tr>
<td align="center">28764</td>
<td align="center">89.59</td>
<td align="center">9172</td>
<td align="center">55.46</td>
</tr>
<tr>
<td align="center">18250</td>
<td align="center">77.33</td>
<td align="center">18250</td>
<td align="center">54.69</td>
</tr>
<tr>
<td align="center">12481</td>
<td align="center">54.11</td>
<td align="center">7410</td>
<td align="center">47.04</td>
</tr>
<tr>
<td align="center">7410</td>
<td align="center">48.11</td>
<td align="center">4744</td>
<td align="center">40.56</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">8489</td>
<td align="center">21.94</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">3984</td>
<td align="center">20.99</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">16196</td>
<td align="center">13.41</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">12481</td>
<td align="center">11.73</td>
</tr>
<tr>
<td></td>
<td></td>
<td align="center">28764</td>
<td align="center">8.17</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec>
<title>Discussion</title>
<p>By using one surface-chemistry for gills and DG, which were assumed to be most affected targets in mussels, it was possible to compare differences in the response of two dissimilar tissues to oil exposure. At first, it was found that both tissues differ remarkably in number, type and amplitude of differentially expressed peptide and protein peaks. Such differences may have an important influence on the performance of classifier generation and could render certain datasets unsuitable for tree-structured analysis (see below). In consequence, the choice of tissue investigated by MS proteomic profiling for the construction of decision tree models is not trivial. It cannot be excluded that some of this observed difference may be attributed to minor alterations in the protein extraction procedure, which were necessary due to the high amounts of lipids in the DG. However, we do not believe these were the most important factors in causing overall differences in the profiles. The DG represents a major site of synthesis and detoxification [
<xref ref-type="bibr" rid="B27">27</xref>
,
<xref ref-type="bibr" rid="B28">28</xref>
]. With respect to the complexity of its functions, it is not surprising that a higher amount of differentially expressed peptides and proteins in response to toxicant exposure could be obtained.</p>
<p>Nonetheless, classifiers have been generated for both gill and DG data sets by subjecting all proteins and peptides with highly significant expression changes to pattern recognition analysis, an artificial learning algorithm that generates supervised classifiers. Generally spoken, the dual aim of carrying out class prediction is i) to evaluate the discriminatory power of differentially expressed proteins and ii) to identify a so-called protein expression signature indicative of the type of exposure. This is based on the following suppositions: first, it is expected that the higher the discriminatory power of a variable, the more relevant it should be as a potential biomarker; second, as suggested earlier [
<xref ref-type="bibr" rid="B1">1</xref>
], the set of proteins itself may serve as a biomarker (thus also-called 'biomarker pattern') as it is able to distinguish between control and exposed animals. The concept, that a set of multiple marker proteins acting in concert produces better classifiers containing a higher level of discriminatory power than could be obtained by a single biomarker alone has been widely acknowledged [
<xref ref-type="bibr" rid="B18">18</xref>
,
<xref ref-type="bibr" rid="B23">23</xref>
].</p>
<p>Unlike the hypothesis driven approach, which investigates known molecules, proteomics may reveal associations between proteins and exposure to contaminants that have not been described previously. However, the identification of several proteins in question is a laborious and time-consuming operation. Even worse, it is derogated by the lack of database information on non-model organisms [
<xref ref-type="bibr" rid="B29">29</xref>
,
<xref ref-type="bibr" rid="B30">30</xref>
], in particular for many invertebrate species such as mussels, which to date are poorly characterised at the genome and proteome level [
<xref ref-type="bibr" rid="B31">31</xref>
]. Previous attempts to identify the constituents of expression signatures in bivalves obtained from the exposure to pollutants documented difficulties in identifying key proteins [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B7">7</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]. Some evidence suggests that a bias towards cytoskeletal proteins emerged [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
,
<xref ref-type="bibr" rid="B32">32</xref>
], presumably due to their relative abundance and prevalence in databases [
<xref ref-type="bibr" rid="B3">3</xref>
]. Therefore, even if informative protein markers can be extracted from proteomic datasets, their identification may remain elusive. While it will definitely aid understanding of the mechanistic role and hence improve the diagnostic reliability [
<xref ref-type="bibr" rid="B1">1</xref>
], it is not imperative to know the identity of an effective biomarker of exposure, given that it would constantly appear under certain conditions [
<xref ref-type="bibr" rid="B2">2</xref>
]. Taking this idea further, a set of proteins and peptides, specific to a particular stressor would constitute the diagnostic marker, independent of their identity. Biomarker pattern recognition based on tree-structured analysis of large-scale proteomic MS data certainly represents a powerful approach to achieve this goal.</p>
<p>It should also be noted that even though the correlation of proteomic alterations with biochemical pathways would render proteomic assays more powerful and should eventually be aspired where possible, in many cases the significance of expression changes is not validated prior to protein identification. Indeed protein identification often represents a means of evaluating such significance by attempting to reveal a biologically plausible function, which in turn would demonstrate their involvement in responses related to pollution exposure [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B7">7</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]. Oberemm
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B17">17</xref>
], for instance, could confirm the differentially expressed proteins of the thymus tissue of marmoset exposed to TCDD to be related to immune function, which is particularly affected by this substance. In ecotoxicology, however, where non-model organisms are frequently employed for monitoring purposes, such confirmation often fails as database matches cannot be obtained for the key proteins of the exposure derived expression signatures [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B6">6</xref>
]. It has been suggested, though discussed critically, that the evaluation of proteomic patterns from MS proteomic profiling can proceed independently from the identities of their proteins [
<xref ref-type="bibr" rid="B18">18</xref>
,
<xref ref-type="bibr" rid="B33">33</xref>
-
<xref ref-type="bibr" rid="B35">35</xref>
] and that classification algorithms would represent a means to extract informative proteins from such data [
<xref ref-type="bibr" rid="B11">11</xref>
,
<xref ref-type="bibr" rid="B13">13</xref>
].</p>
<p>In comparison to earlier studies [
<xref ref-type="bibr" rid="B5">5</xref>
,
<xref ref-type="bibr" rid="B8">8</xref>
,
<xref ref-type="bibr" rid="B9">9</xref>
], however, we experienced limitations to the generation of classifiers resulting in lower sensitivity and specificity of the classification success. Bjørnstad
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B8">8</xref>
] presented discrimination models using mussel haemolymph, which were able to classify oil-exposed mussels with more than 90 % accuracy. Yet, statistical constraints resulting from too many variables in combination with too few samples may have impaired classifier construction and resulted in overfitting of the models [
<xref ref-type="bibr" rid="B22">22</xref>
,
<xref ref-type="bibr" rid="B24">24</xref>
]. In general, few sample numbers per class make it easy to produce seemingly robust classifiers that give excellent results for both LS and TS [
<xref ref-type="bibr" rid="B13">13</xref>
,
<xref ref-type="bibr" rid="B23">23</xref>
]. Accordingly, the reliability of such 'optimal' models may be illusory, notwithstanding their very good classification success. In contrast, when reducing the dimensionality of the dataset prior to classifier construction, the highest classification success was 80% for the DG dataset and sO exposure. To improve the extraction of biologically relevant predictors and to guard against overfitting, we included feature selection as proposed by Levner [
<xref ref-type="bibr" rid="B22">22</xref>
]. Using univariate statistical tests, which rank the variables according to their significance, we decided to only integrate highly significant peaks (
<italic>p </italic>
< 0.001) into the models, thereby also minimising the risk of including falsely differentiating peaks. In addition, we ranked the robustness of the models before maximum classification success by considering only those decision trees with the lowest misclassification costs (
<italic>i.e</italic>
. low error rates). Consequently, the poorer classification success as compared to the previous studies is owed to improved and more prudent data processing.</p>
<p>From the comparison of gill and DG datasets, as well as from the investigation of proteins and peptides with sex-specific differential expression, two major conclusions could be inferred: (1) As the overall expression changes observed with gills resulted in fewer peaks with a significance level of
<italic>p </italic>
< 0.001 (which, in addition, showed notably lower expression changes than those of DG) the gill dataset performed relatively poor in the actual classification of unknown samples. Besides, none of the
<italic>m/z </italic>
values contained in the classifier reaches a twofold expression change. Conversely, the DG dataset was able to construct better decision rules, which were more robust, resulting in a higher actual classification success. It can thus be concluded, that a sufficient number of peaks statistically corroborated, possibly combined with minimum amplitudes of differential expression, is required to enhance the discriminatory power and reliability of a biomarker pattern. Not all datasets will be equally suitable. (2) Splitting the datasets into male and female for each of the tissue, plus the requirement to separate them into LS and independent TS rendered them too small to result in appropriate class prediction. This clearly emphasises the importance of sufficient biological replicates and hence large enough sample sets. Incorporation of gender information into the entire dataset, did not give any satisfactory results as to specific sex-dependent proteins and peptides. Eventually, sex-specifically altered proteins and peptides could only be deduced from statistical significances, which as well have been substantially weakened by reduced sample numbers. Statistical significance by itself, however, does not provide any indication of the diagnostic value of those peaks. Concerning the number of biological replicates, it should be noted that the tree-growing methodology is data intensive. Thus, even though reducing the number of variables for the input space (
<italic>i.e. m/z </italic>
peaks) by introducing a filter, eventually much larger datasets would be needed for classifier validation.</p>
<p>Similar ecotoxicological case-control studies have been analysed recently, either by SELDI or 2-DE [
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B8">8</xref>
-
<xref ref-type="bibr" rid="B10">10</xref>
]. Some of the complex data resulting from those analyses have then been subjected to various data mining techniques such as principal component analysis, hierarchical clustering, non-metric multidimensional scaling and CART, resulting in the conclusion that a set of proteins specific or indicative of the treatment has been obtained, some of which were identified as related to the metabolism of xenobiotics [
<xref ref-type="bibr" rid="B7">7</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]. The strength of those conclusions largely depends on the study design as well as on data-acquisition and data-mining methods. With the ebullient expectations of applying proteomics to ecotoxicology, the consideration of factors compromising these findings may have been insufficient. For instance, constitutive proteomic expression changes are not well investigated and expected to vary widely with differences in age, diet or developmental stage. Accordingly, the changes due to toxicant exposure may not be greater than the noise of protein expression variability [
<xref ref-type="bibr" rid="B36">36</xref>
]. Secondary experimental effects, such as reduced food intake and energy deprivation, may also account for the observed proteomic changes. Many of the alterations are likely to be of a more general nature and would become manifested with quite different types of pollutant-derived stress as well as with otherwise adverse environmental conditions. This might also be one of the reasons for repeatedly identifying ubiquitous genes and proteins in varying experiments, such as cytoskeletal proteins (actin, mysion, tropomysion, tubulin [
<xref ref-type="bibr" rid="B3">3</xref>
,
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]), malate dehydrogenase [
<xref ref-type="bibr" rid="B6">6</xref>
,
<xref ref-type="bibr" rid="B37">37</xref>
], glutathione S-transferase [
<xref ref-type="bibr" rid="B38">38</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
] or proteins associated with physiological pathways of respiration (
<italic>e.g</italic>
. carbonic anhydrase [
<xref ref-type="bibr" rid="B39">39</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
]) and oxidative stress (variants of superoxide dismutase [
<xref ref-type="bibr" rid="B7">7</xref>
,
<xref ref-type="bibr" rid="B10">10</xref>
,
<xref ref-type="bibr" rid="B38">38</xref>
]). Experimental repeatability and reproducibility also represent important information that needs to be investigated and monitored more thoroughly [
<xref ref-type="bibr" rid="B40">40</xref>
]. Eventually, data analysis protocols including data processing steps (
<italic>i.e</italic>
. filtering of noise, data normalisation, peak/spot matching, peak/spot detection and quantification) as well as the data-mining methodology can be largely varied. This in turn will have a significant impact on the protein patterns generated.</p>
<p>Once established in the laboratory where natural variation is reduced, as for any other biomarker, an expression signature has to be confirmed in field situations, most likely with varying sources and compositions of chemicals as it is typically the case for mixture pollution as well as involving different populations. Natural biotic and abiotic conditions will influence the differential expression of proteins at the various sites. These confounding factors already complicate the application of many biomarkers in the field. As an example, Arts
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B41">41</xref>
] have reported very high variability of heat shock protein levels and discussed their suitability for field studies. With mussels, tidal, diurnal and seasonal effects are likely to change baseline proteomic expressions similar to the expression changes observed with specific proteins [
<xref ref-type="bibr" rid="B42">42</xref>
-
<xref ref-type="bibr" rid="B44">44</xref>
]. Even though a pattern of multiple marker proteins is expected to have a higher predictive value than a single biomarker [
<xref ref-type="bibr" rid="B18">18</xref>
], in ecotoxicology the problematic of multifactorial action of confounding factors along with possible non-monotonous concentration response relationships [
<xref ref-type="bibr" rid="B41">41</xref>
] may be actually magnified in the multivariate proteomic profiles. On the other hand, supervised learning methods could possibly ignore non-informative biological variance in the datasets and segregate unspecific from specific responses, thereby selecting those variables only that are able to indicate the particular character of an exposure. Additionally, randomisation and matching of potential confounding factors (age, sex, etc.) prior to data analysis would be likely to prevent biases in the obtained results [
<xref ref-type="bibr" rid="B36">36</xref>
]. As such, the proteomics based biomarker pattern, could integrate specific as well as secondary effects, all of which are part of the organisms' combat against toxic action.</p>
<p>Despite the potential of proteomic patterns as sensitive markers that comprise multiple molecular endpoints with high discriminatory and ideally explanatory power, a careful approach has to be taken in each step of proteomic profiling, biomarker discovery and validation. Errors in study design and execution can lead to misleading results, especially when exploiting the vast datasets using complex multivariate analyses [
<xref ref-type="bibr" rid="B13">13</xref>
]. It then has to be demonstrated whether these patterns are robust enough to be retrieved in the face of biological variability and if they can be related causally or at least linked statistically to higher levels of biological organisation [
<xref ref-type="bibr" rid="B45">45</xref>
], before proteomics may be involved into risk assessment.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>Machine learning and classification algorithms may represent a powerful means to extract relevant information from the large data sets obtained by MS proteomics. However, supervised training of classifiers is prone to overfitting, resulting in excellent classification success and thus has to be conducted with caution. Moreover, CART performs better with larger sample sizes, which could be processed by SELDI due to its high-throughput capacities and the possibilities for standardisation and automatisation of procedures, but may not be available from sample collection. Consequently, the optimal strategy to screen proteomics MS data from similar ecotoxicological studies has yet to be elaborated.</p>
</sec>
<sec sec-type="methods">
<title>Methods</title>
<sec>
<title>Organisms</title>
<p>Adult
<italic>M. edulis </italic>
were collected in November 2002 along the Førlandsfjord nearby Stavanger, Norway. Norwegian authorities had previously classified the Fjord as non-contaminated. The animals were transferred to the laboratory and placed in 123 L tanks per group with a flow rate of 3 L seawater per min. For the continuous flow-through system, natural fjord water of 4°C was taken from a depth of 80 m in the water column [
<xref ref-type="bibr" rid="B46">46</xref>
]. Each tank, one per exposure and control, contained 250 individuals with a size range of 6.2 to 9.0 cm. The mussels were fed every second day with an algae mixture consisting of
<italic>Rhodomonas sp</italic>
. and
<italic>Isochrisis sp</italic>
. Acclimation to laboratory conditions lasted for 13 days prior to starting the experiment.</p>
</sec>
<sec>
<title>Exposure</title>
<p>The Blue mussels were allocated into three groups: i) control mussels (C), ii) mussels exposed to 0.5 mg/L of crude oil (oil) and iii) mussels exposed to 0.5 mg/L crude oil spiked with a mixture of 0.1 mg/L combined APs and PAHs (sO) as nominal concentrations for 21 days. Prior to the insertion of the animals, the surfaces of the exposure tanks have been pre-absorbed with the exposure media for three days. A minimum of 30 males and 30 females were sampled for each group. The determination of the gender was carried out by smeared probes of the gonads observed by light microscopy.</p>
<p>North Sea crude oil (Stratfjord B) was dispersed into seawater mechanically by a mixing valve and a Dispax
<sup>® </sup>
rotator, running at a speed of 10 000 rpm. The droplet size was monitored frequently by a Coulter
<sup>® </sup>
II particle size analyser. Oil concentrations in the water were calculated from the estimated particle size and number [
<xref ref-type="bibr" rid="B47">47</xref>
]. The spike of APs + PAHs represented an additional boost of low molecular size PAHs (Naphthalene plus C1–C3 alkyl homologues, Fluorene, Phenanthrene plus C1–C2 alkyl homologues, Dibenzothiophene plus C1–C2 alkyl homologues; sumPAH: 0.018 mg/L]) and the most common APs found in produced water (p-cresol, m-ethylphenol, 3,5-dimethylphenol, 2,4,6-trimethylphenol, 2-tert-butylphenol, 3-tert-butylphenol, 4-n-butylphenol, 4-pentylphenol; sumAPs: 0.082 mg/L). Both the APs and the PAHs were based upon previous studies of their composition in PW discharges and monitored by HPLC analysis [
<xref ref-type="bibr" rid="B47">47</xref>
,
<xref ref-type="bibr" rid="B48">48</xref>
]. More detailed analysis information can be obtained from Bjørnstadt
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B8">8</xref>
]. A general description of the experimental design is given in Sundt
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B49">49</xref>
].</p>
</sec>
<sec>
<title>Protein extraction and chip preparation</title>
<p>Gills and DG were dissected from live mussels, snap frozen and stored at -80°C until further processing. Previous method optimisation revealed strong anionic exchange chip surface to yield in a maximum number of peaks with both tissue types und thus represented the surface chemistry of choice for proteomic expression profiling [
<xref ref-type="bibr" rid="B5">5</xref>
]. Although not assessed in this study, reproducibility has been evaluated prior to method implementation and was found to be satisfactory. Coefficients of variance for SELDI applications have been reported to be reasonably below 20% indicating that SELDI experiments are consistent at the level of common proteins expressed [
<xref ref-type="bibr" rid="B8">8</xref>
,
<xref ref-type="bibr" rid="B9">9</xref>
,
<xref ref-type="bibr" rid="B50">50</xref>
]. However, according to Listgarten and Emili [
<xref ref-type="bibr" rid="B14">14</xref>
] these are no absolute indicators of quality as coefficients of variance are less informative when feature detection is part of the analytical process. In a quality assessment study conducted by Hong
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B50">50</xref>
] systematic variability across spot position as well as across chips or plates has not been detected, demonstrating good reproducibility of SELDI experiments.</p>
<p>Protein extraction and chip preparation were carried out as described previously [
<xref ref-type="bibr" rid="B5">5</xref>
]. Briefly, tissues were homogenised with 50 mM Tris, pH 7.6, 150 mM NaCl, 0.1% Triton X-100, 10 mM DTT (1:4 w/v for gills and 1:8 w/v for DG) on ice and centrifuged three times for 20 min at 20 000
<italic>g</italic>
, 4°C. To avoid protein breakdown during preparation, a protease inhibitor cocktail (Sigma-Aldrich) was added to the samples. Because lipids generally disturb the binding of proteins to the chip surface, they were absorbed in the lipid rich DG homogenates by Liposorp™ (Calbiochem; 1.5:1 v/v). Total protein content of the supernatant was determined according to the Bradford method [
<xref ref-type="bibr" rid="B51">51</xref>
]. Protein concentration was adjusted to 1 mg/mL prior to sample dilution with binding buffer (50 mM sodium acetate, 0.1% Triton X-100, pH 5.5) resulting in a final protein concentration of 1 μg/mL. Per spot, 20 μg of protein was incubated overnight at 4°C in a 96-well bioprocessor (Ciphergen Biosystems Inc.) on a platform shaker. Three washes with binding buffer devoid of Triton X-100 and two quick washes with ultra pure H
<sub>2</sub>
O subsequently removed unspecific bound proteins. Sinapic acid resolved in 50% ACN/0.1% TFA was applied as a matrix. The arrays were processed according to an automated data collection protocol in a PBS-II protein chip MS reader. Mass accuracy was calibrated externally using the 'All-in-1' molecular mass standard (Ciphergen Biosystems Inc.).</p>
</sec>
<sec>
<title>Spectral processing</title>
<p>The raw intensity data were pre-processed prior to subsequent protein expression profiling using ProteinChip
<sup>® </sup>
software version 3.1.1 (Ciphergen Biosystems Inc.). All spectra of the three different groups were assembled and examined for two separate regions on account of different noise levels [
<xref ref-type="bibr" rid="B52">52</xref>
,
<xref ref-type="bibr" rid="B53">53</xref>
]. Normalisation of peak intensities to total ion current was from
<italic>m/z </italic>
3000 to 20 000 for the low molecular weight area and from
<italic>m/z </italic>
20 000 to 100% of the spectrum size for the high molecular weight range. Despite generally good overall consistency among spectra for SELDI experiments, the reproducibility does not define the quality of individual spectra [
<xref ref-type="bibr" rid="B50">50</xref>
]. Hence, low quality spectra need to be identified [
<xref ref-type="bibr" rid="B40">40</xref>
,
<xref ref-type="bibr" rid="B50">50</xref>
] by defining outliers and should be removed prior to data analyses. This was carried out by using the quartile method on the basis of calculated normalisation factors. For mass-normalisation, identical peaks were utilised in the low and high molecular weight area to assure the compatibility of peak detection. Gill spectra were internally calibrated with three peaks found in all of the spectra (
<italic>m/z </italic>
9242, 16 284 and 40 454). For the calibration of DG spectra four prominent peaks could be utilised (
<italic>m/z </italic>
7206, 9538, 21 643 and 35 687). Peak detection was similar for low and high molecular weight ranges, except for the cluster mass window which was 0.5 % up to
<italic>m/z </italic>
20 000 and 1 % for the rest of the spectral region. The S/N was set at three for the first and two for the second pass with baseline subtraction on. Peaks were required to be present in a minimum of all spectra, equivalent to approximately two-thirds of the spectra of one group (
<italic>i.e</italic>
. 20 % when peaks of C, oil and sO where clustered or 12 % if the spectra where additionally divided into male and female groups). Estimated peaks were added by the software.</p>
</sec>
<sec>
<title>Data analysis</title>
<p>For an initial statistical evaluation of the two datasets, nonparametric Kruskal-Wallis One-way ANOVA was performed on all profiles (C
<italic>vs</italic>
. oil
<italic>vs</italic>
. sO). For the identification of sex-specific peaks, the datasets were split into the male and female fraction of each group. A protein or peptide was considered to be differentially expressed if a statistically significant alteration in its intensity was observed when compared to the control group. The overall significance level was set at 5 % and the variables were ranked according to their statistical significance with three levels of significance (0.05 ≥
<italic>p </italic>
> 0.01; 0.01 ≥
<italic>p </italic>
≥ 0.001 and
<italic>p </italic>
< 0.001).</p>
<p>In view of the variance of peak intensities between individual spectra, and the prevalent low expression changes, we focussed on the highest level of significance (
<italic>p </italic>
< 0.001) for the most informative and robust variables [
<xref ref-type="bibr" rid="B54">54</xref>
]. According to the range of expression changes found in this and earlier studies with SELDI and Blue mussels [
<xref ref-type="bibr" rid="B5">5</xref>
], the ones with more than two-fold up- or-down-regulation were considered more likely to be steadily discriminable from random variation [
<xref ref-type="bibr" rid="B55">55</xref>
]. This empirical cut-off value, corresponding to a log ratio with an absolute value bigger than 0.3, is widely found in literature [
<italic>e.g</italic>
. [
<xref ref-type="bibr" rid="B25">25</xref>
,
<xref ref-type="bibr" rid="B55">55</xref>
]] and has been statistically confirmed by Sabatti
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B26">26</xref>
].</p>
<p>To determine the diagnostic value of marker candidates, the LS was subjected to supervised classification analysis using BiomarkerPattern™ software, version 4.01 (Ciphergen Biosystems Inc.). BiomarkerPattern™ software is an implementation of CART, a nonparametric statistical procedure based on the binary recursive partitioning algorithm introduced by Breiman
<italic>et al</italic>
. [
<xref ref-type="bibr" rid="B56">56</xref>
], with cost-complexity pruning by 10-fold cross-validation. Details of CART analysis have been described elsewhere [
<xref ref-type="bibr" rid="B13">13</xref>
], [
<xref ref-type="bibr" rid="B57">57</xref>
]. Briefly, CART begins with a root node and, through a series of yes/no questions, generates descendent nodes until final classification is reached or further splitting is terminated (
<italic>e.g</italic>
. too few cases). Each split separates a parent node into exactly two child nodes using one rule at a time. Here, the splitting rules depend on the normalised intensity levels of the
<italic>m/z </italic>
values obtained from the SELDI protein expression profile. Once maximal trees are grown, smaller sub-trees are generated by pruning away the branches of the maximal tree and the best tree is determined by testing for error rates (
<italic>i.e</italic>
. costs of misclassification). The intuitive presentation of decision rules in a tree-structured form facilitates the interpretation and applicability of the obtained classifiers.</p>
<p>For the construction of decision tree models, the datasets were divided randomly into a LS which comprised about two thirds of the spectra of the respective groups and a TS consisting of the remaining third (Table
<xref ref-type="table" rid="T2">2</xref>
). Sample statistics, defining the input variables for classifier generation were performed on the LS for C
<italic>vs</italic>
. oil and C
<italic>vs</italic>
. sO by Mann-Whitney U-test. In order to reduce the dimensionality for model construction and to obviate false positive discoveries of potentially discriminating variables [
<xref ref-type="bibr" rid="B13">13</xref>
] the input matrix was restricted to normalised intensity levels with
<italic>p </italic>
< 0.001.</p>
<p>Different models were constructed by varying user-defined criteria (
<italic>e.g</italic>
. Gini power), whereby CART selects the variables in an independent manner for each model. Optimal model selection was carried out by the analyst, not merely for maximum predictive success but obligatory involves low 'error costs' and few decision nodes within the tree-building algorithm in order to ascertain the robustness of the model. Moreover, estimated classification success was supposed not to differ too much between the respective groups. Subsequently, the chosen models were independently tested with the TS, a set of spectra not involved in the generation of the classifier. This was done to verify whether the actual and estimated classification success would be in good agreement, thus signifying the concluding criterion for applicable decision models.</p>
<p>In this study we have examined tree-structured classification with respect to two-class problems only. Theoretically, there is no limitation to the number of categories for classifier generation; also additional information, such as gender can be included. However, in our trials the outcome resulted in poor prediction success or did not provide any valid additional information; consequently those approaches were not included.</p>
</sec>
</sec>
<sec>
<title>Abbreviations</title>
<p>
<bold>ACN</bold>
, acetonitrile;
<bold>APs</bold>
, alkylated phenols;
<bold>C</bold>
, controls;
<bold>CART</bold>
, classification and regression trees,
<bold>2-DE</bold>
, two-dimensional gel electrophoresis;
<bold>DG</bold>
, digestive gland;
<bold>DTT</bold>
, dithiotreitol;
<bold>HPLC</bold>
, high-performance liquid chromatography;
<bold>LS</bold>
, learning set;
<bold>MS</bold>
, mass-spectrometry,
<bold>m/z</bold>
, mass-to-charge ratio;
<bold>oil</bold>
, oil exposure;
<bold>PAHs</bold>
, polycyclic aromatic hydrocarbons;
<bold>PW</bold>
, produced water;
<bold>SD</bold>
, standard deviation;
<bold>SELDI</bold>
, surface-enhanced laser desorption/ionisation time-of-flight mass-spectrometry;
<bold>S/N</bold>
, signal-to-noise ratio;
<bold>sO</bold>
, spiked oil exposure;
<bold>TFA</bold>
, trifluoroacetic acid;
<bold>TCDD</bold>
, 2,3,7,8-Tetrachlorodibenzo-
<italic>p</italic>
-dioxin;
<bold>TS</bold>
, testing set</p>
</sec>
<sec>
<title>Competing interests</title>
<p>The author(s) declare that they have no competing interests.</p>
</sec>
<sec>
<title>Authors' contributions</title>
<p>OKA conceived the study, rose funding and participated in design and coordination of the experiments. TM and TK were involved in the experimental work; carried out sampling, sample preparation and SELDI-TOF mass spectrometry as well as data analyses. FL provided helpful suggestions and assisted in finalising the manuscript. TM and TK established the procedures of data processing and classifier generation. TM drafted the manuscript and TK wrote the final version.</p>
</sec>
</body>
<back>
<ack>
<sec>
<title>Acknowledgements</title>
<p>Financial support was given by the Norwegian Research Council (Grant no. 133724/420) and the EU commission (Grant no. EVK3-CT-2000-00025). We would like to thank Ciphergen Biosystems Inc. for their helpful comments on ProteinChip
<sup>® </sup>
and BiomarkerPattern™ analysis. Thanks are due to Dawn Hallidy for kindly revising the English. We are also indebted to Dr. H.-R. Köhler and Dr. R. Triebskorn, and their group at the Animal Physiological Ecology at the University of Tübingen, Germany for critical discussion of the manuscript.</p>
</sec>
</ack>
<ref-list>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bradley</surname>
<given-names>BP</given-names>
</name>
<name>
<surname>Shrader</surname>
<given-names>EA</given-names>
</name>
<name>
<surname>Kimmel</surname>
<given-names>DG</given-names>
</name>
<name>
<surname>Meiller</surname>
<given-names>JC</given-names>
</name>
</person-group>
<article-title>Protein expression signatures: an application of proteomics</article-title>
<source>Mar Environ Res</source>
<year>2002</year>
<volume>54</volume>
<fpage>373</fpage>
<lpage>377</lpage>
<pub-id pub-id-type="pmid">12408590</pub-id>
<pub-id pub-id-type="doi">10.1016/S0141-1136(02)00115-0</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hogstrand</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Balesaria</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Glover</surname>
<given-names>CN</given-names>
</name>
</person-group>
<article-title>Application of genomics and proteomics for study of the integrated response to zinc exposure in a non-model fish species, the rainbow trout</article-title>
<source>Comp Biochem Physiol</source>
<year>2002</year>
<volume>133B</volume>
<fpage>523</fpage>
<lpage>35</lpage>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodríguez-Ortega</surname>
<given-names>MJ</given-names>
</name>
<name>
<surname>Grøsvik</surname>
<given-names>BE</given-names>
</name>
<name>
<surname>Rodríguez-Ariza</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Goksøyr</surname>
<given-names>A</given-names>
</name>
<name>
<surname>López-Barea</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Changes in protein expression profiles in bivalve molluscs (
<italic>Chamaelea gallina</italic>
) exposed to four model environmental pollutants</article-title>
<source>Proteomics</source>
<year>2003</year>
<volume>3</volume>
<fpage>1535</fpage>
<lpage>1543</lpage>
<pub-id pub-id-type="pmid">12923779</pub-id>
<pub-id pub-id-type="doi">10.1002/pmic.200300491</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shrader</surname>
<given-names>EA</given-names>
</name>
<name>
<surname>Henry</surname>
<given-names>TR</given-names>
</name>
<name>
<surname>Greeley</surname>
<given-names>MS</given-names>
</name>
<name>
<surname>Bradley</surname>
<given-names>BP</given-names>
</name>
</person-group>
<article-title>Proteomics in zebrafish exposed to endocrine disrupting chemicals</article-title>
<source>Ecotoxicology</source>
<year>2003</year>
<volume>12</volume>
<fpage>485</fpage>
<lpage>488</lpage>
<pub-id pub-id-type="pmid">14680328</pub-id>
<pub-id pub-id-type="doi">10.1023/B:ECTX.0000003034.69538.eb</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Knigge</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Monsinjon</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
</person-group>
<article-title>Surface-enhanced laser desorption/ionization-time of flight-mass spectrometry approach to biomarker discovery in Blue mussels (
<italic>Mytilus edulis</italic>
) exposed to polyaromatic hydrocarbons and heavy metals under field conditions</article-title>
<source>Proteomics</source>
<year>2004</year>
<volume>4</volume>
<fpage>2722</fpage>
<lpage>2727</lpage>
<pub-id pub-id-type="pmid">15352246</pub-id>
<pub-id pub-id-type="doi">10.1002/pmic.200300828</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manduzio</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Cosette</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Gricourt</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Jouenne</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Lenz</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
<name>
<surname>Leboulenger</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Rocher</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Proteome modifications of Blue mussel (
<italic>Mytilus edulis </italic>
L.) gills as an effect of water pollution</article-title>
<source>Proteomics</source>
<year>2005</year>
<volume>5</volume>
<fpage>4958</fpage>
<lpage>4963</lpage>
<pub-id pub-id-type="pmid">16267813</pub-id>
<pub-id pub-id-type="doi">10.1002/pmic.200401328</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mi</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Orbea</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Syme</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Ahmed</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Cajaraville</surname>
<given-names>MP</given-names>
</name>
<name>
<surname>Cristóbal</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Peroxisomal proteomics, a new tool for risk assessment of peroxisome proliferating pollutants in the marine environment</article-title>
<source>Proteomics</source>
<year>2005</year>
<volume>5</volume>
<fpage>3954</fpage>
<lpage>3965</lpage>
<pub-id pub-id-type="pmid">16130170</pub-id>
<pub-id pub-id-type="doi">10.1002/pmic.200401243</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bjørnstad</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Larsen</surname>
<given-names>BK</given-names>
</name>
<name>
<surname>Skadsheim</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>MB</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
</person-group>
<article-title>The potential of ecotoxicoproteomics in environmental monitoring: biomarker profiling in mussel plasma using proteinchip array technology</article-title>
<source>J Toxicol Environ Health</source>
<year>2006</year>
<volume>69A</volume>
<fpage>77</fpage>
<lpage>96</lpage>
<pub-id pub-id-type="doi">10.1080/15287390500259277</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gomiero</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Pampanin</surname>
<given-names>DM</given-names>
</name>
<name>
<surname>Bjørnstad</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Larsen</surname>
<given-names>BK</given-names>
</name>
<name>
<surname>Provan</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Lyng</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
</person-group>
<article-title>An ecotoxicoproteomic approach (SELDI-TOF mass spectrometry) to biomarker discovery in crab exposed to pollutants under laboratory conditions</article-title>
<source>Aquat Toxicol</source>
<year>2006</year>
<volume>78S</volume>
<fpage>S34</fpage>
<lpage>S41</lpage>
<pub-id pub-id-type="pmid">16631935</pub-id>
<pub-id pub-id-type="doi">10.1016/j.aquatox.2006.02.013</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Apraiz</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Mi</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Cristobal</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Identification of proteomic signatures of exposure to marine pollutants in mussels (
<italic>Mytilus edulis</italic>
)</article-title>
<source>Mol Cell Proteomics</source>
<year>2006</year>
<volume>5</volume>
<fpage>1274</fpage>
<lpage>1285</lpage>
<pub-id pub-id-type="pmid">16603574</pub-id>
<pub-id pub-id-type="doi">10.1074/mcp.M500333-MCP200</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kell</surname>
<given-names>DB</given-names>
</name>
<name>
<surname>Darby</surname>
<given-names>RM</given-names>
</name>
<name>
<surname>Draper</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Genomic computing. Explanatory analysis of plant expression profiling data using machine learning</article-title>
<source>Plant Physiol</source>
<year>2001</year>
<volume>126</volume>
<fpage>943</fpage>
<lpage>951</lpage>
<pub-id pub-id-type="pmid">11457944</pub-id>
<pub-id pub-id-type="doi">10.1104/pp.126.3.943</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Ackermann</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Jolly</surname>
<given-names>RA</given-names>
</name>
<name>
<surname>Eckstein</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Kulkarni</surname>
<given-names>NH</given-names>
</name>
<name>
<surname>Helvering</surname>
<given-names>LM</given-names>
</name>
<name>
<surname>Goldstein</surname>
<given-names>KM</given-names>
</name>
<name>
<surname>Shou</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Estrem</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Ryan</surname>
<given-names>TP</given-names>
</name>
<name>
<surname>Colet</surname>
<given-names>JM</given-names>
</name>
<name>
<surname>Thomas</surname>
<given-names>CE</given-names>
</name>
<name>
<surname>Stevens</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Onyia</surname>
<given-names>JE</given-names>
</name>
</person-group>
<article-title>Data-driven analysis approach for biomarker discovery using molecular-profiling technologies</article-title>
<source>Biomarkers</source>
<year>2005</year>
<volume>10</volume>
<fpage>153</fpage>
<lpage>172</lpage>
<pub-id pub-id-type="pmid">16076730</pub-id>
<pub-id pub-id-type="doi">10.1080/13547500500107430</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fung</surname>
<given-names>ET</given-names>
</name>
<name>
<surname>Weinberger</surname>
<given-names>SR</given-names>
</name>
<name>
<surname>Gavin</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>F</given-names>
</name>
</person-group>
<article-title>Bioinformatics approaches in clinical proteomics</article-title>
<source>Expert Rev Proteom</source>
<year>2005</year>
<volume>2</volume>
<fpage>847</fpage>
<lpage>862</lpage>
<pub-id pub-id-type="doi">10.1586/14789450.2.6.847</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Listgarten</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Emili</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Statistical and computational methods for comparative proteomic profiling using liquid chromatography-tandem mass spectrometry</article-title>
<source>Mol Cell Proteom</source>
<year>2005</year>
<volume>4</volume>
<fpage>419</fpage>
<lpage>434</lpage>
<pub-id pub-id-type="doi">10.1074/mcp.R500005-MCP200</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aardema</surname>
<given-names>MJ</given-names>
</name>
<name>
<surname>MacGregor</surname>
<given-names>JT</given-names>
</name>
</person-group>
<article-title>Toxicology and genetic toxicology in the new era of 'toxicogenomics': impact of '-omics' technologies</article-title>
<source>Mutat Res</source>
<year>2002</year>
<volume>499</volume>
<fpage>13</fpage>
<lpage>25</lpage>
<pub-id pub-id-type="pmid">11804602</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feron</surname>
<given-names>VJ</given-names>
</name>
<name>
<surname>Groten</surname>
<given-names>JP</given-names>
</name>
</person-group>
<article-title>Toxicological evaluation of chemical mixtures</article-title>
<source>Food Chem Toxicol</source>
<year>2002</year>
<volume>40</volume>
<fpage>825</fpage>
<lpage>839</lpage>
<pub-id pub-id-type="pmid">11983277</pub-id>
<pub-id pub-id-type="doi">10.1016/S0278-6915(02)00021-2</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oberemm</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Meckert</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Brandenburger</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Herzig</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Lindner</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Kalenberg</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Krause</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Ittrich</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Kopp-Schneider</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Stahlmann</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Richter-Reichhelm</surname>
<given-names>HB</given-names>
</name>
<name>
<surname>Gundert-Remy</surname>
<given-names>U</given-names>
</name>
</person-group>
<article-title>Differential signatures of protein expression in marmoset liver and thymus induced by single-dose TCDD treatment</article-title>
<source>Toxicology</source>
<year>2005</year>
<volume>206</volume>
<fpage>33</fpage>
<lpage>48</lpage>
<pub-id pub-id-type="pmid">15590107</pub-id>
<pub-id pub-id-type="doi">10.1016/j.tox.2004.06.061</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Petricoin</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Liotta</surname>
<given-names>LA</given-names>
</name>
</person-group>
<article-title>The vision for a new diagnostic paradigm</article-title>
<source>Clinic Chem</source>
<year>2003</year>
<volume>49</volume>
<fpage>1276</fpage>
<lpage>1278</lpage>
<pub-id pub-id-type="doi">10.1373/49.8.1276</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Varanasi</surname>
<given-names>U</given-names>
</name>
</person-group>
<source>Metabolism of polycyclic aromatic hydrocarbons in the aquatic environment</source>
<year>1989</year>
<publisher-name>Boca Raton: CRC Press</publisher-name>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arukwe</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Celius</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Walther</surname>
<given-names>BT</given-names>
</name>
<name>
<surname>Goksøyr</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Plasma levels of vitellogenin and eggshell
<italic>zona radiata </italic>
proteins in 4-nonylphenol and
<italic>o, p'</italic>
-DDT treated juvenile Atlantic salmon (
<italic>Salmo salar</italic>
)</article-title>
<source>Mar Environ Res</source>
<year>1998</year>
<volume>46</volume>
<fpage>133</fpage>
<lpage>136</lpage>
<pub-id pub-id-type="doi">10.1016/S0141-1136(98)00002-6</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Washburn</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Stone</surname>
<given-names>S</given-names>
</name>
<name>
<surname>MacIntyre</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Dispersion of produced water in a coastal environment and its biological implications</article-title>
<source>Cont Shelf Res</source>
<year>1999</year>
<volume>19</volume>
<fpage>57</fpage>
<lpage>78</lpage>
<pub-id pub-id-type="doi">10.1016/S0278-4343(98)00068-5</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levner</surname>
<given-names>I</given-names>
</name>
</person-group>
<article-title>Feature selection and nearest centroid classification for protein mass spectrometry</article-title>
<source>BMC Bioinformatics</source>
<year>2005</year>
<volume>6</volume>
<fpage>68</fpage>
<pub-id pub-id-type="pmid">15788095</pub-id>
<pub-id pub-id-type="doi">10.1186/1471-2105-6-68</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Somorjai</surname>
<given-names>RL</given-names>
</name>
<name>
<surname>Dolenko</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Baumgartner</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Class prediction and discovery using gene microarray and proteomics mass spectroscopy data: curses, caveats, cautions</article-title>
<source>Bioinformatics</source>
<year>2003</year>
<volume>19</volume>
<fpage>1484</fpage>
<lpage>1491</lpage>
<pub-id pub-id-type="pmid">12912828</pub-id>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btg182</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ransohoff</surname>
<given-names>DF</given-names>
</name>
</person-group>
<article-title>Lessons from controversy: ovarian cancer screening and serum proteomics</article-title>
<source>J Nat Canc Inst</source>
<year>2005</year>
<volume>97</volume>
<fpage>315</fpage>
<lpage>319</lpage>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Klenø</surname>
<given-names>TG</given-names>
</name>
<name>
<surname>Kiehr</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Baunsgaard</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Siedelmann</surname>
<given-names>UG</given-names>
</name>
</person-group>
<article-title>Combination of 'omics' data to investigate the mechanism(s) of hydrazine-induced hepatotoxicity in rats and to identify potential biomarkers</article-title>
<source>Biomarkers</source>
<year>2004</year>
<volume>9</volume>
<fpage>116</fpage>
<lpage>138</lpage>
<pub-id pub-id-type="pmid">15370871</pub-id>
<pub-id pub-id-type="doi">10.1080/13547500410001728408</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sabatti</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Karsten</surname>
<given-names>SL</given-names>
</name>
<name>
<surname>Geschwind</surname>
<given-names>DH</given-names>
</name>
</person-group>
<article-title>Thresholding rules for recovering a sparse signal from microarray experiments</article-title>
<source>Math Biosci</source>
<year>2002</year>
<volume>176</volume>
<fpage>17</fpage>
<lpage>34</lpage>
<pub-id pub-id-type="pmid">11867081</pub-id>
<pub-id pub-id-type="doi">10.1016/S0025-5564(01)00102-X</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kress</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Schmekel</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Nott</surname>
<given-names>JA</given-names>
</name>
</person-group>
<article-title>Ultrastructure of the digestive gland in the opisthobranch mollusk, Runcina. Veliger</article-title>
<source>The Veliger</source>
<year>1994</year>
<volume>37</volume>
<fpage>358</fpage>
<lpage>373</lpage>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lobo-da-Cunha</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>The digestive cells of the hepatopancreas in
<italic>Aplysia depilans </italic>
(Mollusca, Opisthobranchia): ultrastructural and cytochemical study</article-title>
<source>Tissue Cell</source>
<year>2000</year>
<volume>32</volume>
<fpage>49</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="pmid">10798317</pub-id>
<pub-id pub-id-type="doi">10.1054/tice.1999.0082</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Graham</surname>
<given-names>DRM</given-names>
</name>
<name>
<surname>Elliott</surname>
<given-names>ST</given-names>
</name>
<name>
<surname>Eyk</surname>
<given-names>JEV</given-names>
</name>
</person-group>
<article-title>Broad-based proteomic strategies: a practical guide to proteomics and functional screening</article-title>
<source>J Physiol</source>
<year>2005</year>
<volume>563</volume>
<fpage>1</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="pmid">15611015</pub-id>
<pub-id pub-id-type="doi">10.1113/jphysiol.2004.080341</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barret</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Brophy</surname>
<given-names>PM</given-names>
</name>
<name>
<surname>Hamilton</surname>
<given-names>JV</given-names>
</name>
</person-group>
<article-title>Analysing proteomic data</article-title>
<source>Int J Parasitol</source>
<year>2005</year>
<volume>35</volume>
<fpage>543</fpage>
<lpage>553</lpage>
<pub-id pub-id-type="pmid">15826646</pub-id>
<pub-id pub-id-type="doi">10.1016/j.ijpara.2005.01.013</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>López</surname>
<given-names>JL</given-names>
</name>
</person-group>
<article-title>Role of proteomics in taxonomy: the Mytilus complex as a model of study</article-title>
<source>J Chromatogr B Analyt Technol Biomed Life Sci</source>
<year>2005</year>
<volume>815</volume>
<fpage>261</fpage>
<lpage>74</lpage>
<pub-id pub-id-type="pmid">15652815</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>López</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Marina</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Álvarez</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Vázquez</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Application of proteomics for fast identification of species-specific peptides from marine species</article-title>
<source>Proteomics</source>
<year>2002</year>
<volume>2</volume>
<fpage>1658</fpage>
<lpage>65</lpage>
<pub-id pub-id-type="pmid">12469335</pub-id>
<pub-id pub-id-type="doi">10.1002/1615-9861(200212)2:12<1658::AID-PROT1658>3.0.CO;2-4</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Diamandis</surname>
<given-names>EP</given-names>
</name>
</person-group>
<article-title>Proteomic patterns in biological fluids: do they represent the future of cancer diagnosis?</article-title>
<source>Clinic Chem</source>
<year>2003</year>
<volume>49</volume>
<fpage>1272</fpage>
<lpage>1275</lpage>
<pub-id pub-id-type="doi">10.1373/49.8.1272</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Baker</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>In biomarkers we trust?</article-title>
<source>Nat Biotech</source>
<year>2005</year>
<volume>23</volume>
<fpage>297</fpage>
<lpage>304</lpage>
<pub-id pub-id-type="doi">10.1038/nbt0305-297</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robbins</surname>
<given-names>RJ</given-names>
</name>
<name>
<surname>Villanueva</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Tempst</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Distilling cancer biomarkers from the serum peptidome: high technology reading of tea leves or insight to clinical systems biology?</article-title>
<source>J Clinic Oncology</source>
<volume>23</volume>
<fpage>4835</fpage>
<lpage>4837</lpage>
<pub-id pub-id-type="doi">10.1200/JCO.2005.02.912</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>K-M</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J-H</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Design issues in toxiconomics using DNA microarry experiment</article-title>
<source>Toxicol Appl Pharmacol</source>
<year>2005</year>
<volume>207</volume>
<fpage>S200</fpage>
<lpage>208</lpage>
<pub-id pub-id-type="doi">10.1016/j.taap.2005.01.045</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>YK</given-names>
</name>
<name>
<surname>WI</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>SH</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>MY</given-names>
</name>
</person-group>
<article-title>Proteomic analysis of cadmium-induced protein profile alterations from marine alga
<italic>Nannochloropsis oculata</italic>
</article-title>
<source>Ecotoxicology</source>
<year>2005</year>
<volume>14</volume>
<fpage>589</fpage>
<lpage>596</lpage>
<pub-id pub-id-type="pmid">16215694</pub-id>
<pub-id pub-id-type="doi">10.1007/s10646-005-0009-5</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boutet</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Tanguy</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Moraga</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Response of the Pacific oyster
<italic>Crassostrea gigas </italic>
to hydrocarbon contamination under experimental conditions</article-title>
<source>Gene</source>
<year>2004</year>
<volume>329</volume>
<fpage>147</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="pmid">15033537</pub-id>
<pub-id pub-id-type="doi">10.1016/j.gene.2003.12.027</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>David</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Tanguy</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Pichavant</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Moraga</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Response of the Pacific oyster
<italic>Crassostrea gigas </italic>
to hypoxia exposure under experimental conditions</article-title>
<source>FEBS J</source>
<year>2005</year>
<volume>272</volume>
<fpage>5635</fpage>
<lpage>52</lpage>
<pub-id pub-id-type="pmid">16262701</pub-id>
<pub-id pub-id-type="doi">10.1111/j.1742-4658.2005.04960.x</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ligett</surname>
<given-names>WS</given-names>
</name>
<name>
<surname>Barker</surname>
<given-names>PE</given-names>
</name>
<name>
<surname>Semmers</surname>
<given-names>OJ</given-names>
</name>
<name>
<surname>Cazares</surname>
<given-names>LH</given-names>
</name>
</person-group>
<article-title>Measurement reproducibility in the early stages of biomareker development</article-title>
<source>Dis Markers</source>
<year>2004</year>
<volume>20</volume>
<fpage>295</fpage>
<lpage>307</lpage>
<pub-id pub-id-type="pmid">15665389</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arts</surname>
<given-names>M-JSJ</given-names>
</name>
<name>
<surname>Schill</surname>
<given-names>RO</given-names>
</name>
<name>
<surname>Knigge</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Eckwert</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Kammenga</surname>
<given-names>JE</given-names>
</name>
<name>
<surname>Köhler</surname>
<given-names>H-R</given-names>
</name>
</person-group>
<article-title>Stress proteins (hsp70, hsp60) induced in isopods and nematodes by field exposure to metals in a gradient near Avonmouth, UK</article-title>
<source>Ecotoxicology</source>
<year>2004</year>
<volume>13</volume>
<fpage>739</fpage>
<lpage>755</lpage>
<pub-id pub-id-type="pmid">15736846</pub-id>
<pub-id pub-id-type="doi">10.1007/s10646-003-4473-5</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schill</surname>
<given-names>RO</given-names>
</name>
<name>
<surname>Gayle</surname>
<given-names>PM</given-names>
</name>
<name>
<surname>Köhler</surname>
<given-names>HR</given-names>
</name>
</person-group>
<article-title>Daily stress protein (hsp70) cycle in chitons (
<italic>Acanthopleura granulata </italic>
Gmelin, 1791) which inhabit the rocky intertidal shoreline in a tropical ecosystem</article-title>
<source>Comp Biochem Physiol</source>
<year>2002</year>
<volume>131C</volume>
<fpage>253</fpage>
<lpage>258</lpage>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luedeking</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Koehler</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Regulation of expression of multixenobiotic resistance (MXR) genes by environmental factors in the Blue mussel
<italic>Mytilus edulis</italic>
</article-title>
<source>Aquat Toxicol</source>
<year>2004</year>
<volume>69</volume>
<fpage>1</fpage>
<lpage>10</lpage>
<pub-id pub-id-type="pmid">15210293</pub-id>
<pub-id pub-id-type="doi">10.1016/j.aquatox.2004.03.003</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bodin</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Burgeot</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Stanisiere</surname>
<given-names>JY</given-names>
</name>
<name>
<surname>Bocquené</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Menard</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Minier</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Boutet</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Amat</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Cherel</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Budzinski</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Seasonal variations of a battery of biomarkers and physiological indices for the mussel
<italic>Mytilus galloprovincialis </italic>
transplanted into the northwest Mediterranean Sea</article-title>
<source>Comp Biochem Physiol</source>
<year>2004</year>
<volume>138C</volume>
<fpage>411</fpage>
<lpage>427</lpage>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adams</surname>
<given-names>SM</given-names>
</name>
</person-group>
<article-title>Assessing cause and effect of multiple stressors on marine systems</article-title>
<source>Mar Pollut Bull</source>
<year>2005</year>
<volume>51</volume>
<fpage>649</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="pmid">16291182</pub-id>
<pub-id pub-id-type="doi">10.1016/j.marpolbul.2004.11.040</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sanni</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Øysaed</surname>
<given-names>KB</given-names>
</name>
<name>
<surname>Høivangli</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Gaudebert</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>A continuous flow system (CFS) for chronic exposure of aquatic organisms</article-title>
<source>Mar Environ Res</source>
<year>1998</year>
<volume>46</volume>
<fpage>97</fpage>
<lpage>101</lpage>
<pub-id pub-id-type="doi">10.1016/S0141-1136(97)00086-X</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aas</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Baussant</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Balk</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Liewenborg</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
</person-group>
<article-title>PAH metabolites in bile, cytochrome P4501A and DNA adducts as environmental risk parameters for chronic oil exposure: a laboratory experiment with Atlantic cod</article-title>
<source>Aquat Toxicol</source>
<year>2000</year>
<volume>51</volume>
<fpage>241</fpage>
<lpage>58</lpage>
<pub-id pub-id-type="pmid">11064127</pub-id>
<pub-id pub-id-type="doi">10.1016/S0166-445X(00)00108-9</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bechmann</surname>
<given-names>RK</given-names>
</name>
</person-group>
<article-title>Effect of the endocrine disrupter nonylphenol on the marine copepode
<italic>Tisbe battagliai</italic>
</article-title>
<source>Sci Tot Environ</source>
<year>1999</year>
<volume>233</volume>
<fpage>33</fpage>
<lpage>46</lpage>
<pub-id pub-id-type="doi">10.1016/S0048-9697(99)00177-1</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sundt</surname>
<given-names>RC</given-names>
</name>
<name>
<surname>Pampanin</surname>
<given-names>DM</given-names>
</name>
<name>
<surname>Larsen</surname>
<given-names>BK</given-names>
</name>
<name>
<surname>Brede</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Herzke</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Bjørnstad</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>OK</given-names>
</name>
</person-group>
<article-title>The BEEP Stavanger workshop: Mesocosm exposures</article-title>
<source>Aquat Toxicol</source>
<year>2006</year>
<volume>S78</volume>
<fpage>S5</fpage>
<lpage>S12</lpage>
<pub-id pub-id-type="pmid">16635530</pub-id>
<pub-id pub-id-type="doi">10.1016/j.aquatox.2006.02.012</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hong</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Dragan</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Epstein</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Teitel</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Q</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Perkins</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>W</given-names>
</name>
</person-group>
<article-title>Quality control and quality assessment of data from surface-enhanced laser desorption/ionization (SELDI) time-of-flight (TOF) mass spectrometry (MS)</article-title>
<source>BMC Bioinformatics</source>
<year>2005</year>
<volume>S6</volume>
<fpage>S5</fpage>
<pub-id pub-id-type="pmid">16026602</pub-id>
<pub-id pub-id-type="doi">10.1186/1471-2105-6-S2-S5</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bradford</surname>
<given-names>MM</given-names>
</name>
</person-group>
<article-title>A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding</article-title>
<source>Anal Biochem</source>
<year>1976</year>
<volume>72</volume>
<fpage>248</fpage>
<lpage>254</lpage>
<pub-id pub-id-type="pmid">942051</pub-id>
<pub-id pub-id-type="doi">10.1016/0003-2697(76)90527-3</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kozak</surname>
<given-names>KR</given-names>
</name>
<name>
<surname>Amneus</surname>
<given-names>MW</given-names>
</name>
<name>
<surname>Pusey</surname>
<given-names>SM</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Luong</surname>
<given-names>MN</given-names>
</name>
<name>
<surname>Luong</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Reddy</surname>
<given-names>ST</given-names>
</name>
<name>
<surname>Farias-Eisner</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Identification of biomarkers for ovarian cancer using strong anion-exchange ProteinChips: potential use in diagnosis and prognosis</article-title>
<source>PNAS</source>
<year>2003</year>
<volume>100</volume>
<fpage>12343</fpage>
<lpage>12348</lpage>
<pub-id pub-id-type="pmid">14523236</pub-id>
<pub-id pub-id-type="doi">10.1073/pnas.2033602100</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rogers</surname>
<given-names>MA</given-names>
</name>
<name>
<surname>Clarke</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Noble</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Munro</surname>
<given-names>NP</given-names>
</name>
<name>
<surname>Paul</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Selby</surname>
<given-names>PJ</given-names>
</name>
<name>
<surname>Banks</surname>
<given-names>RE</given-names>
</name>
</person-group>
<article-title>Proteomic profiling of urinary proteins in renal cancer by surface-enhanced laser desorption ionization and neural-network analysis: identification of key issues affecting potential clinical utility</article-title>
<source>Cancer Res</source>
<year>2003</year>
<volume>63</volume>
<fpage>6971</fpage>
<lpage>6983</lpage>
<pub-id pub-id-type="pmid">14583499</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Clarke</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Chan</surname>
<given-names>DW</given-names>
</name>
</person-group>
<article-title>The application of clinical proteomics to cancer and other diseases</article-title>
<source>Clin Chem Lab Med</source>
<year>2003</year>
<volume>41</volume>
<fpage>1562</fpage>
<lpage>1570</lpage>
<pub-id pub-id-type="pmid">14708880</pub-id>
<pub-id pub-id-type="doi">10.1515/CCLM.2003.239</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Mohan</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Miyakoshi</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Baylink</surname>
<given-names>DJ</given-names>
</name>
</person-group>
<article-title>Differential protein profile in the ear-punched tissue of regeneration and non-regeneration strains of mice: a novel approach to explore the candidate genes for soft-tissue regeneration</article-title>
<source>Biochim Biophys Acta</source>
<year>2000</year>
<volume>1524</volume>
<fpage>102</fpage>
<lpage>109</lpage>
<pub-id pub-id-type="pmid">11113556</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Breiman</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Friedman</surname>
<given-names>JH</given-names>
</name>
<name>
<surname>Olshen</surname>
<given-names>RA</given-names>
</name>
<name>
<surname>Stone</surname>
<given-names>CJ</given-names>
</name>
</person-group>
<person-group person-group-type="editor">
<name>
<surname>Bickel P, Cleveland W, Dudley R</surname>
</name>
</person-group>
<article-title>Classification and regression trees</article-title>
<source>Wadsworth Statistics/Probability Series</source>
<year>1984</year>
<publisher-name>Wadsworth International Group, TN, USA</publisher-name>
</citation>
</ref>
<ref id="B57">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Steinberg</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Colla</surname>
<given-names>P</given-names>
</name>
</person-group>
<source>CART – Classification and regression trees</source>
<year>1997</year>
<publisher-name>Salford Systems, San Diego, USA</publisher-name>
</citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/France/explor/LeHavreV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000142  | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000142  | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/France
   |area=    LeHavreV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     
   |texte=   
}}

Wicri

This area was generated with Dilib version V0.6.25.
Data generation: Sat Dec 3 14:37:02 2016. Site generation: Tue Mar 5 08:25:07 2024