Statistical methods for comparative phenomics using high-throughput phenotype microarrays

Permutation-Based Inference for Function-on-Scalar Regression With an Application in PET Brain Imaging

The density of various proteins throughout the human brain can be studied through the use of positron emission tomography (PET) imaging. We report here on data from a study of serotonin transporter (5-HTT) binding. While PET imaging data analysis is most commonly performed on data that are aggregated into several discrete a priori regions of interest, in this study, primary interest is on measures of 5-HTT binding potential that are made at many locations along a continuous anatomically defined tract, one that was chosen to follow serotonergic axons. Our goal is to characterize the binding patterns along this tract and also to determine how such patterns differ between control subjects and depressed patients. Due to the nature of our data, we utilize function-on-scalar regression modeling to make optimal use of our data. Inference on both main effects (position along the tract; diagnostic group) and their interactions is made using permutation testing strategies that do not require distributional assumptions. Also, to investigate the question of homogeneity we implement a permutation testing strategy, which adapts a "block bootstrapping" approach from time series analysis to the functional data setting.

Bonferroni-Holm and permutation tests to compare health data: methodological and applicative issues

BACKGROUND

Statistical methodology is a powerful tool in the health research; however, there is wide accord that statistical methodologies are not usually used properly. In particular when multiple comparisons are needed, it is necessary to check the rate of false positive results and the potential inflation of type I errors. In this case, permutation testing methods are useful to check the simultaneous significance level and identify the most significant factors.

METHODS

In this paper an application of permutation tests, in the medical context of Inflammatory Bowel Diseases, is performed. The main goal is to assess the existence of significant differences between Crohn's Disease (CD) and Ulcerative Colitis (UC). The Sequentially Rejective Multiple Test (Bonferroni-Holm procedure) is used to find which of the partial tests are effectively significant and solve the problem of the multiplicity control.

RESULTS

Applying Non-Parametric Combination (NPC) Test for partial and combined tests we conclude that Crohn's Disease patients and Ulcerative Colitis patients differ between them for most examined variables. UC patients compared with the CD patients, have a higher diagnosis age, not show smoking status, proportion of patients treated with immunosuppressants or with biological drugs is lower than the CD patients, even if the duration of such therapies is longer. CD patients have a higher rate of re-hospitalization. Diabetes is more present in the sub-population of UC patients. Analyzing the Charlson score we can highlight that UC patients have a more severe clinical situation than CD patients. Finally, CD patients are more frequently subject to surgery compared to UC. Appling of the Bonferroni Holm procedure, which provided adjusted p-values, we note that only nine of the examined variables are statistically significant: Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Re-hospitalization, Gender and Duration of Immunosoppressive Therapy. Therefore, we can conclude that these are the specific variables that can discriminate effectively the Crohn's Disease and Ulcerative Colitis groups.

CONCLUSIONS

We identified significant variables that discriminate the two groups, satisfying the multiplicity problem, in fact we can affirm that Smoking habit, Immunosuppressive therapy, Surgery, Biological Drug, Diabetes, Adverse Events, Hospitalization, Gender and Duration of Immunosoppressive Therapy are the effectively significant variables.

Comparison of Chemical Sensitivity of Fresh and Long-Stored Heat Resistant Neosartorya fischeri Environmental Isolates Using BIOLOG Phenotype MicroArray System

Spoilage of heat processed food and beverage by heat resistant fungi (HRF) is a major problem for food industry in many countries. Neosartorya fischeri is the leading source of spoilage in thermally processed products. Its resistance to heat processing and toxigenicity makes studies about Neosartorya fischeri metabolism and chemical sensitivity essential. In this study chemical sensitivity of two environmental Neosartorya fischeri isolates were compared. One was isolated from canned apples in 1923 (DSM3700), the other from thermal processed strawberry product in 2012 (KC179765), used as long-stored and fresh isolate, respectively. The study was conducted using Biolog Phenotype MicroArray platforms of chemical sensitivity panel and traditional hole-plate method. The study allowed for obtaining data about Neosartorya fischeri growth inhibitors. The fresh isolate appeared to be much more resistant to chemical agents than the long-stored isolate. Based on phenotype microarray assay nitrogen compounds, toxic cations and membrane function compounds were the most effective in growth inhibition of N. fischeri isolates. According to the study zaragozic acid A, thallium(I) acetate and sodium selenate were potent and promising N. fischeri oriented fungicides which was confirmed by both chemical sensitivity microplates panel and traditional hole-plate methods.

Yersinia pestis in the Age of Big Data

As omics-driven technologies developed rapidly, genomics, transcriptomics, proteomics, metabolomics and other omics-based data have been accumulated in unprecedented speed. Omics-driven big data in biology have changed our way of research. "Big science" has promoted our understanding of biology in a holistic overview that is impossibly achieved by traditional hypothesis-driven research. In this chapter, we gave an overview of omics-driven research on Y. pestis, provided a way of thinking on Yersinia pestis research in the age of big data, and made some suggestions to integrate omics-based data for systems understanding of Y. pestis.

Novel R pipeline for analyzing Biolog Phenotypic MicroArray data

Data produced by Biolog Phenotype MicroArrays are longitudinal measurements of cells’ respiration on distinct substrates. We introduce a three-step pipeline to analyze phenotypic microarray data with novel procedures for grouping, normalization and effect identification. Grouping and normalization are standard problems in the analysis of phenotype microarrays defined as categorizing bacterial responses into active and non-active, and removing systematic errors from the experimental data, respectively. We expand existing solutions by introducing an important assumption that active and non-active bacteria manifest completely different metabolism and thus should be treated separately. Effect identification, in turn, provides new insights into detecting differing respiration patterns between experimental conditions, e.g. between different combinations of strains and temperatures, as not only the main effects but also their interactions can be evaluated. In the effect identification, the multilevel data are effectively processed by a hierarchical model in the Bayesian framework. The pipeline is tested on a data set of 12 phenotypic plates with bacterium Yersinia enterocolitica. Our pipeline is implemented in R language on the top of opm R package and is freely available for research purposes.

Characterization of amygdalin-degrading Lactobacillus species

AIMS

Cyanogenic glycosides are phytotoxic secondary metabolites produced by some crop plants. The aim of this study was to identify lactic acid bacteria (LAB) capable of catabolizing amygdalin, a model cyanogenic glycoside, for use in the biodetoxification of amygdalin-containing foods and feeds.

METHODS AND RESULTS

Amygdalin-catabolizing lactobacilli were characterized using a combination of cultivation-dependent and molecular assays. Lactobacillus paraplantarum and Lactobacillus plantarum grew robustly on amygdalin (Amg(+)), while other LAB species typically failed to catabolize amygdalin (Amg(-)). Interestingly, high concentrations of amygdalin and two of its metabolic derivatives (mandelonitrile and benzaldehyde) inhibited the growth of Lact. plantarum RENO 0093. The differential regulation of genes tentatively involved in cyanohydrin metabolism illustrated that the metabolism of amygdalin- and glucose-grown cultures also differed significantly.

CONCLUSIONS

Amygdalin fermentation was a relatively uncommon phenotype among the LAB and generally limited to strains from the Lact. plantarum group. Phenotype microarrays (PM) enabled strain-level discrimination between closely related strains within a species and suggested that phenotypic differences might affect niche specialization.

SIGNIFICANCE AND IMPACT OF THE STUDY

Amygdalin-degrading lactobacilli with practical application in the biodetoxification of amygdalin were characterized. These strains show potential for use as starter cultures to improve the safety of foods and feeds.

Diet complexity and estrogen receptor β status affect the composition of the murine intestinal microbiota

Intestinal microbial dysbiosis contributes to the dysmetabolism of luminal factors, including steroid hormones (sterones) that affect the development of chronic gastrointestinal inflammation and the incidence of sterone-responsive cancers of the breast, prostate, and colon. Little is known, however, about the role of specific host sterone nucleoreceptors, including estrogen receptor β (ERβ), in microbiota maintenance. Herein, we test the hypothesis that ERβ status affects microbiota composition and determine if such compositionally distinct microbiota respond differently to changes in diet complexity that favor Proteobacteria enrichment. To this end, conventionally raised female ERβ(+/+) and ERβ(-/-) C57BL/6J mice (mean age of 27 weeks) were initially reared on 8604, a complex diet containing estrogenic isoflavones, and then fed AIN-76, an isoflavone-free semisynthetic diet, for 2 weeks. 16S rRNA gene surveys revealed that the fecal microbiota of 8604-fed mice and AIN-76-fed mice differed, as expected. The relative diversity of Proteobacteria, especially the Alphaproteobacteria and Gammaproteobacteria, increased significantly following the transition to AIN-76. Distinct patterns for beneficial Lactobacillales were exclusive to and highly abundant among 8604-fed mice, whereas several Proteobacteria were exclusive to AIN-76-fed mice. Interestingly, representative orders of the phyla Proteobacteria, Bacteroidetes, and Firmicutes, including the Lactobacillales, also differed as a function of murine ERβ status. Overall, these interactions suggest that sterone nucleoreceptor status and diet complexity may play important roles in microbiota maintenance. Furthermore, we envision that this model for gastrointestinal dysbiosis may be used to identify novel probiotics, prebiotics, nutritional strategies, and pharmaceuticals for the prevention and resolution of Proteobacteria-rich dysbiosis.

Metabolic model refinement using phenotypic microarray data

Phenotypic microarray (PM) is a standardized, high-throughput technology for profiling phenotypes of microorganisms, which allows for characterization on around 2,000 different media conditions. The data generated using PM can be incorporated into genome-scale metabolic models to improve their predictive capability. In addition, a comparison of phenotypic profiles of wild-type and gene knockout mutants can give essential information about gene functions of unknown genes. In this chapter, we present a protocol to refine preconstructed metabolic models using the PM data. Both manual refinement and algorithmic approaches for integrating the PM data into metabolic models have been discussed.

Biolog phenotype microarrays

Phenotype microarrays nicely complement traditional genomic, transcriptomic, and proteomic analysis by offering opportunities for researchers to ground microbial systems analysis and modeling in a broad yet quantitative assessment of the organism's physiological response to different metabolites and environments. Biolog phenotype assays achieve this by coupling tetrazolium dyes with minimally defined nutrients to measure the impact of hundreds of carbon, nitrogen, phosphorous, and sulfur sources on redox reactions that result from compound-induced effects on the electron transport chain. Over the years, we have used Biolog's reproducible and highly sensitive assays to distinguish closely related bacterial isolates, to understand their metabolic differences, and to model their metabolic behavior using flux balance analysis. This chapter describes Biolog phenotype microarray system components, reagents, and methods, particularly as they apply to bacterial identification, characterization, and metabolic analysis.

Visualization and curve-parameter estimation strategies for efficient exploration of phenotype microarray kinetics

Background

The Phenotype MicroArray (OmniLog® PM) system is able to simultaneously capture a large number of phenotypes by recording an organism's respiration over time on distinct substrates. This technique targets the object of natural selection itself, the phenotype, whereas previously addressed ‘-omics’ techniques merely study components that finally contribute to it. The recording of respiration over time, however, adds a longitudinal dimension to the data. To optimally exploit this information, it must be extracted from the shapes of the recorded curves and displayed in analogy to conventional growth curves.

Methodology

The free software environment R was explored for both visualizing and fitting of PM respiration curves. Approaches using either a model fit (and commonly applied growth models) or a smoothing spline were evaluated. Their reliability in inferring curve parameters and confidence intervals was compared to the native OmniLog® PM analysis software. We consider the post-processing of the estimated parameters, the optimal classification of curve shapes and the detection of significant differences between them, as well as practically relevant questions such as detecting the impact of cultivation times and the minimum required number of experimental repeats.

Conclusions

We provide a comprehensive framework for data visualization and parameter estimation according to user choices. A flexible graphical representation strategy for displaying the results is proposed, including 95% confidence intervals for the estimated parameters. The spline approach is less prone to irregular curve shapes than fitting any of the considered models or using the native PM software for calculating both point estimates and confidence intervals. These can serve as a starting point for the automated post-processing of PM data, providing much more information than the strict dichotomization into positive and negative reactions. Our results form the basis for a freely available R package for the analysis of PM data.

Application of phenotypic microarrays to environmental microbiology

Environmental organisms are extremely diverse and only a small fraction has been successfully cultured in the laboratory. Culture in micro wells provides a method for rapid screening of a wide variety of growth conditions and commercially available plates contain a large number of substrates, nutrient sources, and inhibitors, which can provide an assessment of the phenotype of an organism. This review describes applications of phenotype arrays to anaerobic and thermophilic microorganisms, use of the plates in stress response studies, in development of culture media for newly discovered strains, and for assessment of phenotype of environmental communities. Also discussed are considerations and challenges in data interpretation and visualization, including data normalization, statistics, and curve fitting.