Linear Mixed-Effects Models in chemistry: A tutorial.

A common goal in chemistry is to study the relationship between a measured signal and the variability of certain factors. To this end, researchers often use Design of Experiment to decide which experiments to conduct and (Multiple) Linear Regression, and/or Analysis of Variance to analyze the collected data. Among the assumptions to the very foundation of this strategy, all the experiments are independent, conditional on the settings of the factors. Unfortunately, due to the presence of uncontrollable factors, real-life experiments often deviate from this assumption, making the data analysis results unreliable. In these cases, Mixed-Effects modeling, despite not being widely used in chemometrics, represents a solid data analysis framework to obtain reliable results. Here we provide a tutorial for Linear Mixed-Effects models. We gently introduce the reader to these models by showing some motivating examples. Then, we discuss the theory behind Linear Mixed-Effect models, and we show how to fit these models by making use of real-life data obtained from an exposome study. Throughout the paper we provide R code so that each researcher is able to implement these useful model themselves.

New insight from MALDI-TOF MS and multivariate data analysis on the botanical origin of polysaccharide-based paint binders in ancient Egypt.

Polysaccharide-based materials of plant origin are known to have been used as binding media in paint and ground layers of artifacts from ancient Egypt, including wall paintings, cartonnages and sarcophagi. The use of gums from Acacia, Astragalus and Prunus genera has been suggested in the literature on the basis of their qualitative or quantitative monosaccharide profile after complete chemical hydrolysis. The introduction of partial enzymatic digestion of the polysaccharide material, followed by analysis of the released oligosaccharides by matrix assisted laser desorption ionization-time-of-flight mass spectrometry, has proved effective in discriminating among gums from different genera, as well as among species within the Acacia genus. In this study, the previously built Acacia database was expanded, principal component analysis (PCA) was used to aid in grouping of the samples, and data interpretation was refined following a modified acacieae taxonomy. Application of the analytical strategy to investigate the paint binders in artworks from ancient Egypt allowed qualitative discrimination of gums at a species level, and provided new insights into the artists' material choices.

Systematic reduction of hyperspectral images for high-throughput plastic characterization.

Hyperspectral Imaging (HSI) combines microscopy and spectroscopy to assess the spatial distribution of spectroscopically active compounds in objects, and has diverse applications in food quality control, pharmaceutical processes, and waste sorting. However, due to the large size of HSI datasets, it can be challenging to analyze and store them within a reasonable digital infrastructure, especially in waste sorting where speed and data storage resources are limited. Additionally, as with most spectroscopic data, there is significant redundancy, making pixel and variable selection crucial for retaining chemical information. Recent high-tech developments in chemometrics enable automated and evidence-based data reduction, which can substantially enhance the speed and performance of Non-Negative Matrix Factorization (NMF), a widely used algorithm for chemical resolution of HSI data. By recovering the pure contribution maps and spectral profiles of distributed compounds, NMF can provide evidence-based sorting decisions for efficient waste management. To improve the quality and efficiency of data analysis on hyperspectral imaging (HSI) data, we apply a convex-hull method to select essential pixels and wavelengths and remove uninformative and redundant information. This process minimizes computational strain and effectively eliminates highly mixed pixels. By reducing data redundancy, data investigation and analysis become more straightforward, as demonstrated in both simulated and real HSI data for plastic sorting.

Water quality monitoring based on chemometric analysis of high-resolution phytoplankton data measured with flow cytometry.

River water is an important source of Dutch drinking water. For this reason, continuous monitoring of river water quality is needed. However, comprehensive chemical analyses with high-resolution gas chromatography [GC]-mass spectrometry [MS]/liquid chromatography [LC]-MS are quite tedious and time consuming; this makes them poorly fit for routine water quality monitoring and, therefore, many pollution events are missed. Phytoplankton are highly sensitive and responsive to toxicity, which makes them highly usable for effect-based water quality monitoring. Flow cytometry can measure the optical properties of phytoplankton every hour, generating a large amount of information-rich data in one year. However, this requires chemometrics, as the resulting fingerprints need to be processed into information about abnormal phytoplankton behaviour. We developed Discriminant Analysis of Multi-Aspect CYtometry (DAMACY) to model the "normal condition" of the phytoplankton community imposed by diurnal, meteorological, and other exogenous influences. DAMACY first describes the cellular variability and distribution of phytoplankton in each measurement using principal component analysis, and then aims to find subtle differences in these phytoplankton distributions that predict normal environmental conditions. Deviations from these normal environmental conditions indicated abnormal phytoplankton behaviour that happened alongside pollution events measured with the GC/MS and LC/MS systems. Thus, our results demonstrate that flow cytometry in combination with chemometrics may be used for an automated hourly assessment of river water quality and as a near real-time early warning for detecting harmful known or unknown contaminants. Finally, both the flow cytometer and the DAMACY algorithm run completely autonomous and only requires maintenance once or twice per year. The warning system results may be uploaded automatically, so that drinking water companies may temporary stop pumping water whenever abnormal phytoplankton behaviour is detected. In the case of prolonged abnormal phytoplankton behaviour, comprehensive analysis may still be used to identify the chemical compound, its origin, and toxicity.

Preanalytical Pitfalls in Untargeted Plasma Nuclear Magnetic Resonance Metabolomics of Endocrine Hypertension.

Despite considerable morbidity and mortality, numerous cases of endocrine hypertension (EHT) forms, including primary aldosteronism (PA), pheochromocytoma and functional paraganglioma (PPGL), and Cushing's syndrome (CS), remain undetected. We aimed to establish signatures for the different forms of EHT, investigate potentially confounding effects and establish unbiased disease biomarkers. Plasma samples were obtained from 13 biobanks across seven countries and analyzed using untargeted NMR metabolomics. We compared unstratified samples of 106 PHT patients to 231 EHT patients, including 104 PA, 94 PPGL and 33 CS patients. Spectra were subjected to a multivariate statistical comparison of PHT to EHT forms and the associated signatures were obtained. Three approaches were applied to investigate and correct confounding effects. Though we found signatures that could separate PHT from EHT forms, there were also key similarities with the signatures of sample center of origin and sample age. The study design restricted the applicability of the corrections employed. With the samples that were available, no biomarkers for PHT vs. EHT could be identified. The complexity of the confounding effects, evidenced by their robustness to correction approaches, highlighted the need for a consensus on how to deal with variabilities probably attributed to preanalytical factors in retrospective, multicenter metabolomics studies.

Pre- versus post-operative untargeted plasma nuclear magnetic resonance spectroscopy metabolomics of pheochromocytoma and paraganglioma.

PURPOSE: Pheochromocytomas and Paragangliomas (PPGL) result in chronic catecholamine excess and serious health complications. A recent study obtained a metabolic signature in plasma from PPGL patients; however, its targeted nature may have generated an incomplete picture and a broader approach could provide additional insights. We aimed to characterize the plasma metabolome of PPGL patients before and after surgery, using an untargeted approach, and to broaden the scope of the investigated metabolic impact of these tumors. DESIGN: A cohort of 36 PPGL patients was investigated. Blood plasma samples were collected before and after surgical tumor removal, in association with clinical and tumor characteristics. METHODS: Plasma samples were analyzed using untargeted nuclear magnetic resonance (NMR) spectroscopy metabolomics. The data were evaluated using a combination of uni- and multi-variate statistical methods. RESULTS: Before surgery, patients with a nonadrenergic tumor could be distinguished from those with an adrenergic tumor based on their metabolic profiles. Tyrosine levels were significantly higher in patients with high compared to those with low BMI. Comparing subgroups of pre-operative samples with their post-operative counterparts, we found a metabolic signature that included ketone bodies, glucose, organic acids, methanol, dimethyl sulfone and amino acids. Three signals with unclear identities were found to be affected. CONCLUSIONS: Our study suggests that the pathways of glucose and ketone body homeostasis are affected in PPGL patients. BMI-related metabolite levels were also found to be altered, potentially linking muscle atrophy to PPGL. At baseline, patient metabolomes could be discriminated based on their catecholamine phenotype.

Improved modelling for low-correlated multiple responses by common-subset-of-independent-variables partial-least-squares.

In this study, a new approach for PLS modelling for low-correlated multiple responses, called Common-Subset-of-Independent-Variables Partial-Least-Squares, denoted as CSIV-PLS1, is proposed and evaluated. In CSIV-PLS1, for each response vector, individual PLS1 models with individual model complexities are developed, based on one common set of independent variables, obtained after variable selection by the Final Complexity Adapted Models method, using the absolute values of the PLS regression coefficients, denoted as FCAM-REG. CSIV-PLS1 combines a common variable set for all response vectors, which is a characteristic of PLS2, with the individual model complexity for each response, which is a characteristic of PLS1. These characteristics make CSIV-PLS1 more flexible than PLS2. The selective and predictive abilities of the proposed CSIV-PLS1 method are investigated using one simulated and four real data sets with low-correlated multiple responses from different sources. The simulated data set is used to test the general applicability of the CSIV-PLS1 method. The predictive abilities, measured by the RMSEP values, resulting from CSIV-PLS1 models, are statistically compared with those of the corresponding PLS1 and PLS2 models, using one-tailed paired t-tests. The selective ability of the CSIV-PLS1 method is good, because mostly variables with an informative meaning to the responses are selected. The RMSEP values resulting from the CSIV-PLS1 method are (i) significantly lower at the 95% confidence level than those of the corresponding PLS2 method, and (ii) borderline significantly lower at the 90-95% confidence level than those of the corresponding PLS1 methods. In case of low-correlated multiple responses, the predictive ability of the CSIV-PLS1 method is significantly better than that of the PLS2 method, and borderline significantly better than those of the corresponding PLS1 methods. Therefore, CSIV-PLS1 modelling may be an alternative for PLS1 or PLS2.

High-throughput single cell data analysis - A tutorial.

White blood cells protect the body against disease but may also cause chronic inflammation, auto-immune diseases or leukemia. There are many different white blood cell types whose identity and function can be studied by measuring their protein expression. Therefore, high-throughput analytical instruments were developed to measure multiple proteins on millions of single cells. The information-rich biochemistry information may only be fully extracted using multivariate statistics. Here we show an overview of the most essential steps for multivariate data analysis of single cell data. We used white blood cells (immunology) as a case study, but a similar approach may be used in environment or biotech research. The first step is analyzing the study design and subsequently formulating a research question. The three main designs are immunophenotyping (finding different cell types), cell activation and rare cell discovery. When preparing the data it is essential to consider the design and focus on the cell type of interest by removing all unwanted events. After pre-processing, the ten-thousands to millions of single cells per sample need to be converted into a cellular distribution. For immunophenotyping a clustering method such as Self-Organizing Maps is useful and for cell activation a model that describes the covariance such as Principal Component Analysis is useful. In rare cell discovery it is useful to first model all common cells and remove them to find the rare cells. Finally discriminant analysis based on the cellular distribution may highlight which cell (sub)types are different between groups.

Nicotinamide Riboside Improves Ataxia Scores and Immunoglobulin Levels in Ataxia Telangiectasia.

BACKGROUND: Treatment of animal models with ataxia telangiectasia (A-T) with nicotinamide riboside (NR) improved their neurological outcome and survival. OBJECTIVE: The aim of this study is to investigate the effects of NR in patients with A-T. METHODS: In this open-label, proof-of-concept study, 24 patients with A-T were treated with NR during four consecutive months. The effects of NR on ataxia, dysarthria, quality of life, and laboratory parameters were analyzed. RESULTS: During treatment, ataxia scores improved; mean total Scale for the Assessment and Rating of Ataxia and International Cooperative Ataxia Rating Scale scores decreased to 2.4 and 10.1 points, respectively. After NR withdrawal, ataxia scores worsened. In immunodeficient patients, the mean serum IgG concentration increased substantially until the end of the study period with 0.52 g/L. Untargeted metabolomics analysis revealed increased plasma levels of NR metabolites and purine nucleosides during treatment. Adverse effects did not occur. CONCLUSIONS: Treatment with NR is tolerated well and associated with improvement in ataxia and serum immunoglobulin concentrations in patients with A-T. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Flow cytometric evaluation of the neutrophil compartment in COVID-19 at hospital presentation: A normal response to an abnormal situation.

Coronavirus disease 2019 (COVID-19) is a rapidly emerging pandemic disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Critical COVID-19 is thought to be associated with a hyper-inflammatory process that can develop into acute respiratory distress syndrome, a critical disease normally mediated by dysfunctional neutrophils. This study tested the hypothesis whether the neutrophil compartment displays characteristics of hyperinflammation in COVID-19 patients. Therefore, a prospective study was performed on all patients with suspected COVID-19 presenting at the emergency room of a large academic hospital. Blood drawn within 2 d after hospital presentation was analyzed by point-of-care automated flow cytometry and compared with blood samples collected at later time points. COVID-19 patients did not exhibit neutrophilia or eosinopenia. Unexpectedly neutrophil activation markers (CD11b, CD16, CD10, and CD62L) did not differ between COVID-19-positive patients and COVID-19-negative patients diagnosed with other bacterial/viral infections, or between COVID-19 severity groups. In all patients, a decrease was found in the neutrophil maturation markers indicating an inflammation-induced left shift of the neutrophil compartment. In COVID-19 this was associated with disease severity.

Variable Selection in Untargeted Metabolomics and the Danger of Sparsity.

The goal of metabolomics is to measure as many metabolites as possible in order to capture biomarkers that may indicate disease mechanisms. Variable selection in chemometric methods can be divided into the following two groups: (1) sparse methods that find the minimal set of variables to discriminate between groups and (2) methods that find all variables important for discrimination. Such important variables can be summarized into metabolic pathways using pathway analysis tools like Mummichog. As a test case, we studied the metabolic effects of treatment with nicotinamide riboside, a form of vitamin B3, in a cohort of patients with ataxia-telangiectasia. Vitamin B3 is an important co-factor for many enzymatic reactions in the human body. Thus, the variable selection method was expected to find vitamin B3 metabolites and also other secondary metabolic changes during treatment. However, sparse methods did not select any vitamin B3 metabolites despite the fact that these metabolites showed a large difference when comparing intensity before and during treatment. Univariate analysis or significance multivariate correlation (sMC) in combination with pathway analysis using Mummichog were able to select vitamin B3 metabolites. Moreover, sMC analysis found additional metabolites. Therefore, in our comparative study, sMC displayed the best performance for selection of relevant variables.

Multi-set Pre-processing of Multicolor Flow Cytometry Data.

Flow Cytometry is an analytical technology to simultaneously measure multiple markers per single cell. Ten thousands to millions of single cells can be measured per sample and each sample may contain a different number of cells. All samples may be bundled together, leading to a 'multi-set' structure. Many multivariate methods have been developed for Flow Cytometry data but none of them considers this structure in their quantitative handling of the data. The standard pre-processing used by existing multivariate methods provides models mainly influenced by the samples with more cells, while such a model should provide a balanced view of the biomedical information within all measurements. We propose an alternative 'multi-set' preprocessing that corrects for the difference in number of cells measured, balancing the relative importance of each multi-cell sample in the data while using all data collected from these expensive analyses. Moreover, one case example shows how multi-set pre-processing may benefit removal of undesired measurement-to-measurement variability and another where class-based multi-set pre-processing enhances the studied response upon comparison to the control reference samples. Our results show that adjusting data analysis algorithms to consider this multi-set structure may greatly benefit immunological insight and classification performance of Flow Cytometry data.

A novel data fusion method for the effective analysis of multiple panels of flow cytometry data.

Multicolour flow cytometry (MFC) is used to measure multiple cellular markers at the single-cell level. Cellular markers may be coloured with different panels of fluorescently-labelled antibodies to enable cell identification or the detection of activated cells in pre-defined, 'gated' specific cell subsets. The number of markers that can be used per measurement is technologically limited however, requiring every panel to be analysed in a separate aliquot measurement. The combined analyses of these dedicated panels may enhance the predictive ability of these measurements and could enrich the interpretation of the immunological information. Here we introduce a fusion method for MFC data, based on DAMACY (Discriminant Analysis of Multi-Aspect Cytometry data), which can combine information from complementary panels. This approach leads to both enhanced predictions and clearer interpretations in comparison with the analysis of separate measurements. We illustrate this method using two datasets: the response of neutrophils evoked by a systemic endotoxin challenge and the activated immune status of the innate cells, T cells and B cells in obese versus lean individuals. The data fusion approach was able to detect cells that do not individually show a difference between clinical phenotypes but do play a role in combination with other cells.

Multi-dimensional flow cytometry analysis reveals increasing changes in the systemic neutrophil compartment during seven consecutive days of endurance exercise.

Endurance exercise is associated with a transient increase in neutrophil counts in the peripheral blood. Here we investigate the impact of intensified endurance exercise on the neutrophil compartment. We hypothesized that intensified endurance exercise leads to mobilization of neutrophil subsets, which are normally absent in the blood. Furthermore, we followed the potential build-up of neutrophil activation and the impact on overnight recovery of the neutrophil compartment during a seven-day cycling tour. The neutrophil compartment was studied in 28 healthy amateur cyclists participating in an eight-day strenuous cycling tour. Blood samples were taken at baseline, after 4 days and after 7 days of cycling. The neutrophil compartment was analyzed in terms of numbers and its phenotype by deep phenotyping of flow cytometry data with the multi-dimensional analysis method FLOOD. Repeated endurance exercise led to a gradual increase in total neutrophil counts over the days leading to a 1.26 fold-increase (95%CI 1.01-1.51 p = 0.0431) in the morning of day 8. Flow cytometric measurements revealed the appearance of 2 additional neutrophil subsets: CD16brightCD62Ldim and CD16dimCD62Lbright. A complex change in neutrophil phenotypes was present characterized by decreased expression of both CD11b and CD62L and marked increased expression of LAIR-1, VLA-4 and CBRM1/5. The changes in expression were found on all neutrophils present in the blood. Strikingly, in strong contrast to our findings during acute inflammation evoked by LPS challenge, these neutrophils did not upregulate classical degranulation markers. In fact, our FLOOD analysis revealed that the exercise induced neutrophil phenotype did not overlap with the neutrophil subsets arising upon acute inflammation. In conclusion, during multiple days of endurance exercise the neutrophil compartment does not regain homeostasis overnight. Thereby our study supports the concept of a build-up of inflammatory cues during repeated endurance exercise training, causing a prolonged change of the systemic neutrophil compartment.

Automated flow cytometric identification of disease-specific cells by the ECLIPSE algorithm.

Multicolor Flow Cytometry (MFC)-based gating allows the selection of cellular (pheno)types based on their unique marker expression. Current manual gating practice is highly subjective and may remove relevant information to preclude discovery of cell populations with specific co-expression of multiple markers. Only multivariate approaches can extract such aspects of cell variability from multi-dimensional MFC data. We describe the novel method ECLIPSE (Elimination of Cells Lying in Patterns Similar to Endogeneity) to identify and characterize aberrant cells present in individuals out of homeostasis. ECLIPSE combines dimensionality reduction by Simultaneous Component Analysis with Kernel Density Estimates. A Difference between Densities (DbD) is used to eliminate cells in responder samples that overlap in marker expression with cells of controls. Thereby, subsequent data analyses focus on the immune response-specific cells, leading to more informative and focused models. To prove the power of ECLIPSE, we applied the method to study two distinct datasets: the in vivo neutrophil response induced by systemic endotoxin challenge and in studying the heterogeneous immune-response of asthmatics. ECLIPSE described the well-characterized common response in the LPS challenge insightfully, while identifying slight differences between responders. Also, ECLIPSE enabled characterization of the immune response associated to asthma, where the co-expressions between all markers were used to stratify patients according to disease-specific cell profiles.

Novel data analysis method for multicolour flow cytometry links variability of multiple markers on single cells to a clinical phenotype.

Multicolour Flow Cytometry (MFC) produces multidimensional analytical data on the quantitative expression of multiple markers on single cells. This data contains invaluable biomedical information on (1) the marker expressions per cell, (2) the variation in such expression across cells, (3) the variability of cell marker expression across samples that (4) may vary systematically between cells collected from donors and patients. Current conventional and even advanced data analysis methods for MFC data explore only a subset of these levels. The Discriminant Analysis of MultiAspect CYtometry (DAMACY) we present here provides a comprehensive view on health and disease responses by integrating all four levels. We validate DAMACY by using three distinct datasets: in vivo response of neutrophils evoked by systemic endotoxin challenge, the clonal response of leukocytes in bone marrow of acute myeloid leukaemia (AML) patients, and the complex immune response in blood of asthmatics. DAMACY provided good accuracy 91-100% in the discrimination between health and disease, on par with literature values. Additionally, the method provides figures that give insight into the marker expression and cell variability for more in-depth interpretation, that can benefit both physicians and biomedical researchers to better diagnose and monitor diseases that are reflected by changes in blood leukocytes.

Increasing conclusiveness of clinical breath analysis by improved baseline correction of multi capillary column - ion mobility spectrometry (MCC-IMS) data.

Current challenges of clinical breath analysis include large data size and non-clinically relevant variations observed in exhaled breath measurements, which should be urgently addressed with competent scientific data tools. In this study, three different baseline correction methods are evaluated within a previously developed data size reduction strategy for multi capillary column - ion mobility spectrometry (MCC-IMS) datasets. Introduced for the first time in breath data analysis, the Top-hat method is presented as the optimum baseline correction method. A refined data size reduction strategy is employed in the analysis of a large breathomic dataset on a healthy and respiratory disease population. New insights into MCC-IMS spectra differences associated with respiratory diseases are provided, demonstrating the additional value of the refined data analysis strategy in clinical breath analysis.