Introduction Batch effects in large untargeted metabolomics experiments are almost unavoidable, especially when sensitive detection techniques like mass spectrometry (MS) are employed. In order to obtain peak intensities that are comparable across all batches, corrections need to be performed. Since non-detects, i.e., signals with an intensity too low to be detected with certainty, are common in metabolomics studies, the batch correction methods need to take these into account. Objectives This paper aims to compare several batch correction methods, and investigates the effect of different strategies for handling non-detects. Methods Batch correction methods usually consist of regression models, possibly also accounting for trends within batches. To fit these models quality control samples (QCs), injected at regular intervals, can be used. Also study samples can be used, provided that the injection order is properly randomized. Normalization methods, not using information on batch labels or injection order, can correct for batch effects as well. Introducing two easy-to-use quality criteria, we assess the merits of these batch correction strategies using three large LC–MS and GC–MS data sets of samples from Arabidopsis thaliana. Results The three data sets have very different characteristics, leading to clearly distinct behaviour of the batch correction strategies studied. Explicit inclusion of information on batch and injection order in general leads to very good corrections; when enough QCs are available, also general normalization approaches perform well. Several approaches are shown to be able to handle non-detects—replacing them with very small numbers such as zero seems the worst of the approaches considered. Conclusion The use of quality control samples for batch correction leads to good results when enough QCs are available. If an experiment is properly set up, batch correction using the study samples usually leads to a similar high-quality correction, but has the advantage that more metabolites are corrected. The strategy for handling non-detects is important: choosing small values like zero can lead to suboptimal batch corrections.
Crop yield stability requires an attenuation of the reduction of yield losses caused by environmental stresses such as drought. Using a combination of metabolomics and high-throughput colorimetric assays, we analysed central metabolism and oxidative stress status in the flag leaf of 292 indica rice (Oryza sativa) accessions. Plants were grown in the field and were, at the reproductive stage, exposed to either well-watered or drought conditions to identify the metabolic processes associated with drought-induced grain yield loss. Photorespiration, protein degradation, and nitrogen recycling were the main processes involved in the drought-induced leaf metabolic reprogramming. Molecular markers of drought tolerance and sensitivity in terms of grain yield were identified using a multivariate model based on the values of the metabolites and enzyme activities across the population. The model highlights the central role of the ascorbate–glutathione cycle, particularly dehydroascorbate reductase, in minimizing drought-induced grain yield loss. In contrast, malondialdehyde was an accurate biomarker for grain yield loss, suggesting that drought-induced lipid peroxidation is the major constraint under these conditions. These findings highlight new breeding targets for improved rice grain yield stability under drought.
In order to better understand the milk proteome and its changes from colostrum to mature milk, samples taken at seven time points in the first 9 days from 4 individual cows were analyzed using proteomic techniques. Both the similarity in changes from day 0 to day 9 in the quantitative milk proteome, and the differences in specific protein abundance, were observed among four cows. One third of the quantified proteins showed a significant decrease in concentration over the first 9 days after calving, especially in the immune proteins (as much as 40 fold). Three relative high abundant enzymes (XDH, LPL, and RNASE1) and cell division and proliferation protein (CREG1) may be involved in the maturation of the gastro-intestinal tract. In addition, high correlations between proteins involved in complement and blood coagulation cascades illustrates the complex nature of biological interrelationships between milk proteins. The linear decrease of protease inhibitors and proteins involved in innate and adaptive immune system implies a protective role for protease inhibitor against degradation. In conclusion, the results found in this study not only improve our understanding of the role of colostrum in both host defense and development of the newborn calf but also provides guidance for the improvement of infant formula through better understanding of the complex interactions between milk proteins.
Fortification of food with iron is considered to be an effective approach to counter the global health problem caused by iron deficiency. However, reactivity of iron with the catechol moiety of food phenolics leads to discolouration and impairs bioavailability. In this study, we investigated the interplay between intrinsic and extrinsic factors on food discolouration caused by iron-catechol complexation. To this end, a three-level fractional factorial design was implemented. Absorbance spectra were analysed using statistical methods, including PCA, HCA, and ANOVA. Furthermore, a direct link between absorbance spectra and stoichiometry of the iron-catechol complexes was confirmed by ESI-Q-TOF-MS. All statistical methods confirm that the main effects affecting discolouration were type of iron salt, pH, and temperature. Additionally, several two-way interactions, such as type of iron salt × pH, pH × temperature, and type of iron salt × concentration significantly affected iron-catechol complexation. Our findings provide insight into iron-phenolic complexation-mediated discolouration, and facilitate the design of iron-fortified foods.