Analyzing postprandial metabolomics data using multiway models: a simulation study.

BACKGROUND: Analysis of time-resolved postprandial metabolomics data can improve the understanding of metabolic mechanisms, potentially revealing biomarkers for early diagnosis of metabolic diseases and advancing precision nutrition and medicine. Postprandial metabolomics measurements at several time points from multiple subjects can be arranged as a subjects by metabolites by time points array. Traditional analysis methods are limited in terms of revealing subject groups, related metabolites, and temporal patterns simultaneously from such three-way data. RESULTS: We introduce an unsupervised multiway analysis approach based on the CANDECOMP/PARAFAC (CP) model for improved analysis of postprandial metabolomics data guided by a simulation study. Because of the lack of ground truth in real data, we generate simulated data using a comprehensive human metabolic model. This allows us to assess the performance of CP models in terms of revealing subject groups and underlying metabolic processes. We study three analysis approaches: analysis of fasting-state data using principal component analysis, T0-corrected data (i.e., data corrected by subtracting fasting-state data) using a CP model and full-dynamic (i.e., full postprandial) data using CP. Through extensive simulations, we demonstrate that CP models capture meaningful and stable patterns from simulated meal challenge data, revealing underlying mechanisms and differences between diseased versus healthy groups. CONCLUSIONS: Our experiments show that it is crucial to analyze both fasting-state and T0-corrected data for understanding metabolic differences among subject groups. Depending on the nature of the subject group structure, the best group separation may be achieved by CP models of T0-corrected or full-dynamic data. This study introduces an improved analysis approach for postprandial metabolomics data while also shedding light on the debate about correcting baseline values in longitudinal data analysis.

From affinity selection to kinetic selection in Germinal Centre modelling.

Affinity maturation is an evolutionary process by which the affinity of antibodies (Abs) against specific antigens (Ags) increases through rounds of B-cell proliferation, somatic hypermutation, and positive selection in germinal centres (GC). The positive selection of B cells depends on affinity, but the underlying mechanisms of affinity discrimination and affinity-based selection are not well understood. It has been suggested that selection in GC depends on both rapid binding of B-cell receptors (BcRs) to Ags which is kinetically favourable and tight binding of BcRs to Ags, which is thermodynamically favourable; however, it has not been shown whether a selection bias for kinetic properties is present in the GC. To investigate the GC selection bias towards rapid and tight binding, we developed an agent-based model of GC and compared the evolution of founder B cells with initially identical low affinities but with different association/dissociation rates for Ag presented by follicular dendritic cells in three Ag collection mechanisms. We compared an Ag collection mechanism based on association/dissociation rates of B-cell interaction with presented Ag, which includes a probabilistic rupture of bonds between the B-cell and Ag (Scenario-1) with a reference scenario based on an affinity-based Ag collection mechanism (Scenario-0). Simulations showed that the mechanism of Ag collection affects the GC dynamics and the GC outputs concerning fast/slow (un)binding of B cells to FDC-presented Ags. In particular, clones with lower dissociation rates outcompete clones with higher association rates in Scenario-1, while remaining B cells from clones with higher association rates reach higher affinities. Accordingly, plasma cell and memory B cell populations were biased towards B-cell clones with lower dissociation rates. Without such probabilistic ruptures during the Ag extraction process (Scenario-2), the selective advantage for clones with very low dissociation rates diminished, and the affinity maturation level of all clones decreased to the reference level.

Exploring dynamic metabolomics data with multiway data analysis: a simulation study.

BACKGROUND: Analysis of dynamic metabolomics data holds the promise to improve our understanding of underlying mechanisms in metabolism. For example, it may detect changes in metabolism due to the onset of a disease. Dynamic or time-resolved metabolomics data can be arranged as a three-way array with entries organized according to a subjects mode, a metabolites mode and a time mode. While such time-evolving multiway data sets are increasingly collected, revealing the underlying mechanisms and their dynamics from such data remains challenging. For such data, one of the complexities is the presence of a superposition of several sources of variation: induced variation (due to experimental conditions or inborn errors), individual variation, and measurement error. Multiway data analysis (also known as tensor factorizations) has been successfully used in data mining to find the underlying patterns in multiway data. To explore the performance of multiway data analysis methods in terms of revealing the underlying mechanisms in dynamic metabolomics data, simulated data with known ground truth can be studied. RESULTS: We focus on simulated data arising from different dynamic models of increasing complexity, i.e., a simple linear system, a yeast glycolysis model, and a human cholesterol model. We generate data with induced variation as well as individual variation. Systematic experiments are performed to demonstrate the advantages and limitations of multiway data analysis in analyzing such dynamic metabolomics data and their capacity to disentangle the different sources of variations. We choose to use simulations since we want to understand the capability of multiway data analysis methods which is facilitated by knowing the ground truth. CONCLUSION: Our numerical experiments demonstrate that despite the increasing complexity of the studied dynamic metabolic models, tensor factorization methods CANDECOMP/PARAFAC(CP) and Parallel Profiles with Linear Dependences (Paralind) can disentangle the sources of variations and thereby reveal the underlying mechanisms and their dynamics.

Human Blood Lipoprotein Predictions from 1H NMR Spectra: Protocol, Model Performances, and Cage of Covariance.

Lipoprotein subfractions are biomarkers for the early diagnosis of cardiovascular diseases. The reference method, ultracentrifugation, for measuring lipoproteins is time-consuming, and there is a need to develop a rapid method for cohort screenings. This study presents partial least-squares regression models developed using 1H nuclear magnetic resonance (NMR) spectra and concentrations of lipoproteins as measured by ultracentrifugation on 316 healthy Danes. This study explores, for the first time, different regions of the 1H NMR spectrum representing signals of molecules in lipoprotein particles and different lipid species to develop parsimonious, reliable, and optimal prediction models. A total of 65 lipoprotein main and subfractions were predictable with high accuracy, Q2 of >0.6, using an optimal spectral region (1.4-0.6 ppm) containing methylene and methyl signals from lipids. The models were subsequently tested on an independent cohort of 290 healthy Swedes with predicted and reference values matching by up to 85-95%. In addition, an open software tool was developed to predict lipoproteins concentrations in human blood from standardized 1H NMR spectral recordings.

Genetic Components of Root Architecture Remodeling in Response to Salt Stress.

Salinity of the soil is highly detrimental to plant growth. Plants respond by a redistribution of root mass between main and lateral roots, yet the genetic machinery underlying this process is still largely unknown. Here, we describe the natural variation among 347 Arabidopsis thaliana accessions in root system architecture (RSA) and identify the traits with highest natural variation in their response to salt. Salt-induced changes in RSA were associated with 100 genetic loci using genome-wide association studies. Two candidate loci associated with lateral root development were validated and further investigated. Changes in CYP79B2 expression in salt stress positively correlated with lateral root development in accessions, and cyp79b2 cyp79b3 double mutants developed fewer and shorter lateral roots under salt stress, but not in control conditions. By contrast, high HKT1 expression in the root repressed lateral root development, which could be partially rescued by addition of potassium. The collected data and multivariate analysis of multiple RSA traits, available through the Salt_NV_Root App, capture root responses to salinity. Together, our results provide a better understanding of effective RSA remodeling responses, and the genetic components involved, for plant performance in stress conditions.

Training for translation between disciplines: a philosophy for life and data sciences curricula.

Motivation: Our society has become data-rich to the extent that research in many areas has become impossible without computational approaches. Educational programmes seem to be lagging behind this development. At the same time, there is a growing need not only for strong data science skills, but foremost for the ability to both translate between tools and methods on the one hand, and application and problems on the other. Results: Here we present our experiences with shaping and running a masters' programme in bioinformatics and systems biology in Amsterdam. From this, we have developed a comprehensive philosophy on how translation in training may be achieved in a dynamic and multidisciplinary research area, which is described here. We furthermore describe two requirements that enable translation, which we have found to be crucial: sufficient depth and focus on multidisciplinary topic areas, coupled with a balanced breadth from adjacent disciplines. Finally, we present concrete suggestions on how this may be implemented in practice, which may be relevant for the effectiveness of life science and data science curricula in general, and of particular interest to those who are in the process of setting up such curricula. Supplementary information: Supplementary data are available at Bioinformatics online.

Repeatability and reproducibility of lipoprotein particle profile measurements in plasma samples by ultracentrifugation.

Background Characterization of lipoprotein particle profiles (LPPs) (including main classes and subclasses) by means of ultracentrifugation (UC) is highly requested given its clinical potential. However, rapid methods are required to replace the very labor-intensive UC method and one solution is to calibrate rapid nuclear magnetic resonance (NMR)-based prediction models, but the reliability of the UC-response method required for the NMR calibration has been largely overlooked. Methods This study provides a comprehensive repeatability and reproducibility study of various UC-based lipid measurements (cholesterol, triglycerides [TGs], free cholesterol, phospholipids, apolipoprotein [apo]A1 and apoB) in different main classes and subclasses of 25 duplicated fresh plasma samples and of 42 quality control (QC) frozen pooled plasma samples of healthy individuals. Results Cholesterol, apoA1 and apoB measurements were very repeatable in all classes (intraclass correlation coefficient [ICC]: 92.93%-99.54%). Free cholesterol and phospholipid concentrations in main classes and subclasses and TG concentrations in high-density lipoproteins (HDL), HDL subclasses and low-density lipoproteins (LDL) subclasses, showed worse repeatability (ICC: 19.21%-99.08%) attributable to low concentrations, variability introduced during UC and assay limitations. On frozen QC samples, the reproducibility of cholesterol, apoA1 and apoB concentrations was found to be better than for the free cholesterol, phospholipids and TGs concentrations. Conclusions This study shows that for LPPs measurements near or below the limit of detection (LOD) in some of the subclasses, as well as the use of frozen samples, results in worsened repeatability and reproducibility. Furthermore, we show that the analytical assay coupled to UC for free cholesterol and phospholipids have different repeatability and reproducibility. All of this needs to be taken into account when calibrating future NMR-based models.

Stoichiometry, Absolute Abundance, and Localization of Proteins in the Bacillus cereus Spore Coat Insoluble Fraction Determined Using a QconCAT Approach.

Spores of Bacillus cereus pose a threat to food safety due to their high resistance to the heat or acid treatments commonly used to make food microbiologically safe. Spores may survive these treatments and later resume growth either on foodstuffs or, after ingestion, upon entering the gut they are capable of producing toxins, which cause either vomiting or diarrhea. The outer layers of the spore, the spore coat and exosporium, consist primarily of proteins that may serve as potential biomarkers for detection. The major morphogenetic protein CotE is important for correct assembly and attachment of the outermost layer, the exosporium, and by extension retention of many proteins. However, characterization of the proteins affected by deletion of CotE has been limited to electrophoretic patterns. Here we report the effect of CotE deletion on the insoluble fraction of the spore proteome through liquid chromatography-Fourier transform tandem mass spectrometry (LC-FTMS/MS) analysis. A total of 560 proteins have been identified in both mutant and wild-type spore coat isolates. A further 163 proteins were identified exclusively in wild-type spore isolates indicating that they are dependent on CotE for their association with the spore. Several of these are newly confirmed as associated with the exosporium, namely BC_2569 (BclF), BC_3345, BC_2427, BC_2878, BC_0666, BC_2984, BC_3481, and BC_2570. A total of 153 proteins were only identified in ΔCotE spore isolates. This was observed for proteins that are known or likely to be interacting with or are encased by CotE. Crucial spore proteins were quantified using a QconCAT reference standard, the first time this was used in a biochemically heterogeneous system. This allowed us to determine the absolute abundance of 21 proteins, which spanned across three orders of magnitude and together covered 5.66% ± 0.51 of the total spore weight. Applying the QconCAT methodology to the ΔCotE mutant allowed us to quantify 4.13% ± 0.14 of the spore total weight and revealed a reduction in abundance for most known exosporium associated proteins upon CotE deletion. In contrast, several proteins, either known or likely to be interacting with or encased by CotE (i.e., GerQ), were more abundant. The results obtained provide deeper insight into the layered spore structure such as which proteins are exposed on the outside of the spore. This information is important for developing detection methods for targeting spores in a food safety setting. Furthermore, protein stoichiometry and determination of the abundance of germination mediating enzymes provides useful information for germination and outgrowth model development.

Global testing of shifts in metabolic phenotype.

INTRODUCTION: Current metabolomics approaches to unravel impact of diet- or lifestyle induced phenotype variation and shifts predominantly deploy univariate or multivariate approaches, with a posteriori interpretation at pathway level. This however often provides only a fragmented view on the involved metabolic pathways. OBJECTIVES: To demonstrate the feasibility of using Goeman's global test (GGT) for assessment of variation and shifts in metabolic phenotype at the level of a priori defined pathways. METHODS: Two intervention studies with identified phenotype variations and shifts were examined. In a weight loss (WL) intervention study obese subjects received a mixed meal challenge before and after WL. In a polyphenol (PP) intervention study obese subjects received a high fat mixed meal challenge (61E% fat) before and after a PP intervention. Plasma samples were obtained at fasting and during the postprandial response. Besides WL- and PP-induced phenotype shifts, also correlation of plasma metabolome with phenotype descriptors was assessed at pathway level. The plasma metabolome covered organic acids, amino acids, biogenic amines, acylcarnitines and oxylipins. RESULTS: For the population of the WL study, GGT revealed that HOMA correlated with the fasting levels of the TCA cycle, BCAA catabolism, the lactate, arginine-proline and phenylalanine-tyrosine pathways. For the population of the PP study, HOMA correlated with fasting metabolite levels of TCA cycle, fatty acid oxidation and phenylalanine-tyrosine pathways. These correlations were more pronounced for metabolic pathways in the fasting state, than during the postprandial response. The effect of the WL and PP intervention on a priori defined metabolic pathways, and correlation of pathways with insulin sensitivity as described by HOMA was in line with previous studies. CONCLUSION: GGT confirmed earlier biological findings in a hypothesis led approach. A main advantage of GGT is that it provides a direct view on involvement of a priori defined pathways in phenotype shifts.

Toward Reliable Lipoprotein Particle Predictions from NMR Spectra of Human Blood: An Interlaboratory Ring Test.

Lipoprotein profiling of human blood by 1H nuclear magnetic resonance (NMR) spectroscopy is a rapid and promising approach to monitor health and disease states in medicine and nutrition. However, lack of standardization of measurement protocols has prevented the use of NMR-based lipoprotein profiling in metastudies. In this study, a standardized NMR measurement protocol was applied in a ring test performed across three different laboratories in Europe on plasma and serum samples from 28 individuals. Data was evaluated in terms of (i) spectral differences, (ii) differences in LPD predictions obtained using an existing prediction model, and (iii) agreement of predictions with cholesterol concentrations in high- and low-density lipoproteins (HDL and LDL) particles measured by standardized clinical assays. ANOVA-simultaneous component analysis (ASCA) of the ring test spectral ensemble that contains methylene and methyl peaks (1.4-0.6 ppm) showed that 97.99% of the variance in the data is related to subject, 1.62% to sample type (serum or plasma), and 0.39% to laboratory. This interlaboratory variation is in fact smaller than the maximum acceptable intralaboratory variation on quality control samples. It is also shown that the reproducibility between laboratories is good enough for the LPD predictions to be exchangeable when the standardized NMR measurement protocol is followed. With the successful implementation of this protocol, which results in reproducible prediction of lipoprotein distributions across laboratories, a step is taken toward bringing NMR more into scope of prognostic and diagnostic biomarkers, reducing the need for less efficient methods such as ultracentrifugation or high-performance liquid chromatography (HPLC).

Time-series analysis of the transcriptome and proteome of Escherichia coli upon glucose repression.

Time-series transcript- and protein-profiles were measured upon initiation of carbon catabolite repression in Escherichia coli, in order to investigate the extent of post-transcriptional control in this prototypical response. A glucose-limited chemostat culture was used as the CCR-free reference condition. Stopping the pump and simultaneously adding a pulse of glucose, that saturated the cells for at least 1h, was used to initiate the glucose response. Samples were collected and subjected to quantitative time-series analysis of both the transcriptome (using microarray analysis) and the proteome (through a combination of 15N-metabolic labeling and mass spectrometry). Changes in the transcriptome and corresponding proteome were analyzed using statistical procedures designed specifically for time-series data. By comparison of the two sets of data, a total of 96 genes were identified that are post-transcriptionally regulated. This gene list provides candidates for future in-depth investigation of the molecular mechanisms involved in post-transcriptional regulation during carbon catabolite repression in E. coli, like the involvement of small RNAs.

Identification of Analytical Factors Affecting Complex Proteomics Profiles Acquired in a Factorial Design Study with Analysis of Variance: Simultaneous Component Analysis.

Complex shotgun proteomics peptide profiles obtained in quantitative differential protein expression studies, such as in biomarker discovery, may be affected by multiple experimental factors. These preanalytical factors may affect the measured protein abundances which in turn influence the outcome of the associated statistical analysis and validation. It is therefore important to determine which factors influence the abundance of peptides in a complex proteomics experiment and to identify those peptides that are most influenced by these factors. In the current study we analyzed depleted human serum samples to evaluate experimental factors that may influence the resulting peptide profile such as the residence time in the autosampler at 4 °C, stopping or not stopping the trypsin digestion with acid, the type of blood collection tube, different hemolysis levels, differences in clotting times, the number of freeze-thaw cycles, and different trypsin/protein ratios. To this end we used a two-level fractional factorial design of resolution IV (2(IV)(7-3)). The design required analysis of 16 samples in which the main effects were not confounded by two-factor interactions. Data preprocessing using the Threshold Avoiding Proteomics Pipeline (Suits, F.; Hoekman, B.; Rosenling, T.; Bischoff, R.; Horvatovich, P. Anal. Chem. 2011, 83, 7786-7794, ref 1) produced a data-matrix containing quantitative information on 2,559 peaks. The intensity of the peaks was log-transformed, and peaks having intensities of a low t-test significance (p-value > 0.05) and a low absolute fold ratio (<2) between the two levels of each factor were removed. The remaining peaks were subjected to analysis of variance (ANOVA)-simultaneous component analysis (ASCA). Permutation tests were used to identify which of the preanalytical factors influenced the abundance of the measured peptides most significantly. The most important preanalytical factors affecting peptide intensity were (1) the hemolysis level, (2) stopping trypsin digestion with acid, and (3) the trypsin/protein ratio. This provides guidelines for the experimentalist to keep the ratio of trypsin/protein constant and to control the trypsin reaction by stopping it with acid at an accurately set pH. The hemolysis level cannot be controlled tightly as it depends on the status of a patient's blood (e.g., red blood cells are more fragile in patients undergoing chemotherapy) and the care with which blood was sampled (e.g., by avoiding shear stress). However, its level can be determined with a simple UV spectrophotometric measurement and samples with extreme levels or the peaks affected by hemolysis can be discarded from further analysis. The loadings of the ASCA model led to peptide peaks that were most affected by a given factor, for example, to hemoglobin-derived peptides in the case of the hemolysis level. Peak intensity differences for these peptides were assessed by means of extracted ion chromatograms confirming the results of the ASCA model.

Culturing Synechocystis sp. Strain PCC 6803 with N2 and CO2 in a Diel Regime Reveals Multiphase Glycogen Dynamics with Low Maintenance Costs.

UNLABELLED: Investigating the physiology of cyanobacteria cultured under a diel light regime is relevant for a better understanding of the resulting growth characteristics and for specific biotechnological applications that are foreseen for these photosynthetic organisms. Here, we present the results of a multiomics study of the model cyanobacterium Synechocystis sp. strain PCC 6803, cultured in a lab-scale photobioreactor in physiological conditions relevant for large-scale culturing. The culture was sparged with N2 and CO2, leading to an anoxic environment during the dark period. Growth followed the availability of light. Metabolite analysis performed with (1)H nuclear magnetic resonance analysis showed that amino acids involved in nitrogen and sulfur assimilation showed elevated levels in the light. Most protein levels, analyzed through mass spectrometry, remained rather stable. However, several high-light-response proteins and stress-response proteins showed distinct changes at the onset of the light period. Microarray-based transcript analysis found common patterns of ∼56% of the transcriptome following the diel regime. These oscillating transcripts could be grouped coarsely into genes that were upregulated and downregulated in the dark period. The accumulated glycogen was degraded in the anaerobic environment in the dark. A small part was degraded gradually, reflecting basic maintenance requirements of the cells in darkness. Surprisingly, the largest part was degraded rapidly in a short time span at the end of the dark period. This degradation could allow rapid formation of metabolic intermediates at the end of the dark period, preparing the cells for the resumption of growth at the start of the light period. IMPORTANCE: Industrial-scale biotechnological applications are anticipated for cyanobacteria. We simulated large-scale high-cell-density culturing of Synechocystis sp. PCC 6803 under a diel light regime in a lab-scale photobioreactor. In BG-11 medium, Synechocystis grew only in the light. Metabolite analysis grouped the collected samples according to the light and dark conditions. Proteome analysis suggested that the majority of enzyme-activity regulation was not hierarchical but rather occurred through enzyme activity regulation. An abrupt light-on condition induced high-light-stress proteins. Transcript analysis showed distinct patterns for the light and dark periods. Glycogen gradually accumulated in the light and was rapidly consumed in the last quarter of the dark period. This suggests that the circadian clock primed the cellular machinery for immediate resumption of growth in the light.

Bioinformatics and systems biology: bridging the gap between heterogeneous student backgrounds.

Teaching students with very diverse backgrounds can be extremely challenging. This article uses the Bioinformatics and Systems Biology MSc in Amsterdam as a case study to describe how the knowledge gap for students with heterogeneous backgrounds can be bridged. We show that a mix in backgrounds can be turned into an advantage by creating a stimulating learning environment for the students. In the MSc Programme, conversion classes help to bridge differences between students, by mending initial knowledge and skill gaps. Mixing students from different backgrounds in a group to solve a complex task creates an opportunity for the students to reflect on their own abilities. We explain how a truly interdisciplinary approach to teaching helps students of all backgrounds to achieve the MSc end terms. Moreover, transferable skills obtained by the students in such a mixed study environment are invaluable for their later careers.

Capturing Arabidopsis root architecture dynamics with ROOT-FIT reveals diversity in responses to salinity.

The plant root is the first organ to encounter salinity stress, but the effect of salinity on root system architecture (RSA) remains elusive. Both the reduction in main root (MR) elongation and the redistribution of the root mass between MRs and lateral roots (LRs) are likely to play crucial roles in water extraction efficiency and ion exclusion. To establish which RSA parameters are responsive to salt stress, we performed a detailed time course experiment in which Arabidopsis (Arabidopsis thaliana) seedlings were grown on agar plates under different salt stress conditions. We captured RSA dynamics with quadratic growth functions (root-fit) and summarized the salt-induced differences in RSA dynamics in three growth parameters: MR elongation, average LR elongation, and increase in number of LRs. In the ecotype Columbia-0 accession of Arabidopsis, salt stress affected MR elongation more severely than LR elongation and an increase in LRs, leading to a significantly altered RSA. By quantifying RSA dynamics of 31 different Arabidopsis accessions in control and mild salt stress conditions, different strategies for regulation of MR and LR meristems and root branching were revealed. Different RSA strategies partially correlated with natural variation in abscisic acid sensitivity and different Na(+)/K(+) ratios in shoots of seedlings grown under mild salt stress. Applying root-fit to describe the dynamics of RSA allowed us to uncover the natural diversity in root morphology and cluster it into four response types that otherwise would have been overlooked.

A critical assessment of feature selection methods for biomarker discovery in clinical proteomics.

In this paper, we compare the performance of six different feature selection methods for LC-MS-based proteomics and metabolomics biomarker discovery-t test, the Mann-Whitney-Wilcoxon test (mww test), nearest shrunken centroid (NSC), linear support vector machine-recursive features elimination (SVM-RFE), principal component discriminant analysis (PCDA), and partial least squares discriminant analysis (PLSDA)-using human urine and porcine cerebrospinal fluid samples that were spiked with a range of peptides at different concentration levels. The ideal feature selection method should select the complete list of discriminating features that are related to the spiked peptides without selecting unrelated features. Whereas many studies have to rely on classification error to judge the reliability of the selected biomarker candidates, we assessed the accuracy of selection directly from the list of spiked peptides. The feature selection methods were applied to data sets with different sample sizes and extents of sample class separation determined by the concentration level of spiked compounds. For each feature selection method and data set, the performance for selecting a set of features related to spiked compounds was assessed using the harmonic mean of the recall and the precision (f-score) and the geometric mean of the recall and the true negative rate (g-score). We conclude that the univariate t test and the mww test with multiple testing corrections are not applicable to data sets with small sample sizes (n = 6), but their performance improves markedly with increasing sample size up to a point (n > 12) at which they outperform the other methods. PCDA and PLSDA select small feature sets with high precision but miss many true positive features related to the spiked peptides. NSC strikes a reasonable compromise between recall and precision for all data sets independent of spiking level and number of samples. Linear SVM-RFE performs poorly for selecting features related to the spiked compounds, even though the classification error is relatively low.

Simulation of the rate of transfer of antibiotic resistance between Escherichia coli strains cultured under well controlled environmental conditions.

It was demonstrated that the tetracycline resistance plasmid in Escherichia coli resembling K-12 23:06 containing the E. coli plasmid DM0133 could be transferred to tetracycline sensitive E. coli K-12 MG1655 YFP. The sensitive recipient strain has a slight metabolic advantage in continuous fermentation in absence of tetracycline pressure and as a result the numbers of the resistant recipient strain increase during fermentation. In presence of tetracycline pressure the sensitive strain is eliminated, but when it acquires tetracycline resistance the strain has still the same metabolic advantage as its sensitive parent strain in absence of tetracycline. Here a model will be shown that could explain the rate of transformation of a sensitive into a resistant recipient strain and its subsequent growth during continuous fermentation. According to the model the probability of formation of mutants would be much higher at the dilution rate of 0.09 compared to 0.28, whereas the growth of mutants would be much faster at high dilution rate. The growth model shows how the recipient mutants and the donor cells behave in relation to the dilution rate and the number of mutants. Apart from a deterministic model describing the growth rate of both the donor strain and the resistant recipient strain a stochastic model was developed that is particularly useful when low numbers of mutants are formed.

Generic framework for high-dimensional fixed-effects ANOVA.

In functional genomics it is more rule than exception that experimental designs are used to generate the data. The samples of the resulting data sets are thus organized according to this design and for each sample many biochemical compounds are measured, e.g. typically thousands of gene-expressions or hundreds of metabolites. This results in high-dimensional data sets with an underlying experimental design. Several methods have recently become available for analyzing such data while utilizing the underlying design. We review these methods by putting them in a unifying and general framework to facilitate understanding the (dis-)similarities between the methods. The biological question dictates which method to use and the framework allows for building new methods to accommodate a range of such biological questions. The framework is built on well known fixed-effect ANOVA models and subsequent dimension reduction. We present the framework both in matrix algebra as well as in more insightful geometrical terms. We show the workings of the different special cases of our framework with a real-life metabolomics example from nutritional research and a gene-expression example from the field of virology.

Understanding repertoire sequencing data through a multiscale computational model of the germinal center.

Sequencing of B-cell and T-cell immune receptor repertoires helps us to understand the adaptive immune response, although it only provides information about the clonotypes (lineages) and their frequencies and not about, for example, their affinity or antigen (Ag) specificity. To further characterize the identified clones, usually with special attention to the particularly abundant ones (dominant), additional time-consuming or expensive experiments are generally required. Here, we present an extension of a multiscale model of the germinal center (GC) that we previously developed to gain more insight in B-cell repertoires. We compare the extent that these simulated repertoires deviate from experimental repertoires established from single GCs, blood, or tissue. Our simulations show that there is a limited correlation between clonal abundance and affinity and that there is large affinity variability among same-ancestor (same-clone) subclones. Our simulations suggest that low-abundance clones and subclones, might also be of interest since they may have high affinity for the Ag. We show that the fraction of plasma cells (PCs) with high B-cell receptor (BcR) mRNA content in the GC does not significantly affect the number of dominant clones derived from single GCs by sequencing BcR mRNAs. Results from these simulations guide data interpretation and the design of follow-up experiments.