Picky with peakpicking: assessing chromatographic peak quality with simple metrics in metabolomics.

BACKGROUND: Chromatographic peakpicking continues to represent a significant bottleneck in automated LC-MS workflows. Uncontrolled false discovery rates and the lack of manually-calibrated quality metrics require researchers to visually evaluate individual peaks, requiring large amounts of time and breaking replicability. This problem is exacerbated in noisy environmental datasets and for novel separation methods such as hydrophilic interaction columns in metabolomics, creating a demand for a simple, intuitive, and robust metric of peak quality. RESULTS: Here, we manually labeled four HILIC oceanographic particulate metabolite datasets to assess the performance of individual peak quality metrics. We used these datasets to construct a predictive model calibrated to the likelihood that visual inspection by an MS expert would include a given mass feature in the downstream analysis. We implemented two novel peak quality metrics, a custom signal-to-noise metric and a test of similarity to a bell curve, both calculated from the raw data in the extracted ion chromatogram, and found that these outperformed existing measurements of peak quality. A simple logistic regression model built on two metrics reduced the fraction of false positives in the analysis from 70-80% down to 1-5% and showed minimal overfitting when applied to novel datasets. We then explored the implications of this quality thresholding on the conclusions obtained by the downstream analysis and found that while only 10% of the variance in the dataset could be explained by depth in the default output from the peakpicker, approximately 40% of the variance was explained when restricted to high-quality peaks alone. CONCLUSIONS: We conclude that the poor performance of peakpicking algorithms significantly reduces the power of both univariate and multivariate statistical analyses to detect environmental differences. We demonstrate that simple models built on intuitive metrics and derived from the raw data are more robust and can outperform more complex models when applied to new data. Finally, we show that in properly curated datasets, depth is a major driver of variability in the marine microbial metabolome and identify several interesting metabolite trends for future investigation.

Glycine betaine uptake and metabolism in marine microbial communities.

Glycine betaine (GBT) is a compatible solute in high concentrations in marine microorganisms. As a component of labile organic matter, GBT has complex biochemical potential as a substrate for microbial use that is unconstrained in the environment. Here we determine the uptake kinetics and metabolic fate of GBT in two natural microbial communities in the North Pacific characterized by different nitrate concentrations. Dissolved GBT had maximum uptake rates of 0.36 and 0.56 nM h-1 with half-saturation constants of 79 and 11 nM in the high nitrate and low nitrate stations respectively. During multiday incubations, most GBT taken into cells was retained as a compatible solute. Stable isotopes derived from the added GBT were also observed in other metabolites, including choline, carnitine and sarcosine, suggesting that GBT was used for biosynthesis and for catabolism to pyruvate and ammonium. Where nitrate was scarce, GBT was primarily metabolized via demethylation to glycine. Gene transcript data were consistent with SAR11 using GBT as a source of methyl groups to fuel the methionine cycle. Where nitrate concentrations were higher, more GBT was partitioned for lipid biosynthesis by both bacteria and eukaryotic phytoplankton. Our data highlight unexpected metabolic pathways and potential routes of microbial metabolite exchange.

Does Formation of Multicellular Colonies by Choanoflagellates Affect Their Susceptibility to Capture by Passive Protozoan Predators?

Microbial eukaryotes, critical links in aquatic food webs, are unicellular, but some, such as choanoflagellates, form multicellular colonies. Are there consequences to predator avoidance of being unicellular vs. forming larger colonies? Choanoflagellates share a common ancestor with animals and are used as model organisms to study the evolution of multicellularity. Escape in size from protozoan predators is suggested as a selective factor favoring evolution of multicellularity. Heterotrophic protozoans are categorized as suspension feeders, motile raptors, or passive predators that eat swimming prey which bump into them. We focused on passive predation and measured the mechanisms responsible for the susceptibility of unicellular vs. multicellular choanoflagellates, Salpingoeca helianthica, to capture by passive heliozoan predators, Actinosphaerium nucleofilum, which trap prey on axopodia radiating from the cell body. Microvideography showed that unicellular and colonial choanoflagellates entered the predator's capture zone at similar frequencies, but a greater proportion of colonies contacted axopodia. However, more colonies than single cells were lost during transport by axopodia to the cell body. Thus, feeding efficiency (proportion of prey entering the capture zone that were engulfed in phagosomes) was the same for unicellular and multicellular prey, suggesting that colony formation is not an effective defense against such passive predators.