Application of two machine learning algorithms to genetic association studies in the presence of covariates

BACKGROUND: Population-based investigations aimed at uncovering genotype-trait associations often involve high-dimensional genetic polymorphism data as well as information on multiple environmental and clinical parameters. Machine learning (ML) algorithms offer a straightforward analytic approach for selecting subsets of these inputs that are most predictive of a pre-defined trait. The performance of these algorithms, however, in the presence of covariates is not well characterized. METHODS AND RESULTS: In this manuscript, we investigate two approaches: Random Forests (RFs) and Multivariate Adaptive Regression Splines (MARS). Through multiple simulation studies, the performance under several underlying models is evaluated. An application to a cohort of HIV-1 infected individuals receiving anti-retroviral therapies is also provided. CONCLUSION: Consistent with more traditional regression modeling theory, our findings highlight the importance of considering the nature of underlying gene-covariate-trait relationships before applying ML algorithms, particularly when there is potential confounding or effect mediation

Production of phi mesons at mid-rapidity in sqrt(s_NN) = 200 GeV Au+Au collisions at RHIC

We present the first results of meson production in the K^+K^- decay channel from Au+Au collisions at sqrt(s_NN) = 200 GeV as measured at mid-rapidity by the PHENIX detector at RHIC. Precision resonance centroid and width values are extracted as a function of collision centrality. No significant variation from the PDG accepted values is observed. The transverse mass spectra are fitted with a linear exponential function for which the derived inverse slope parameter is seen to be constant as a function of centrality. These data are also fitted by a hydrodynamic model with the result that the freeze-out temperature and the expansion velocity values are consistent with the values previously derived from fitting single hadron inclusive data. As a function of transverse momentum the collisions scaled peripheral.to.central yield ratio RCP for the is comparable to that of pions rather than that of protons. This result lends support to theoretical models which distinguish between baryons and mesons instead of particle mass for explaining the anomalous proton yield.Comment: 326 authors, 24 pages text, 23 figures, 6 tables, RevTeX 4. To be submitted to Physical Review C as a regular article. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Diffractive Dijet Production at sqrt(s)=630 and 1800 GeV at the Fermilab Tevatron

We report a measurement of the diffractive structure function $F_{jj}^D$ of the antiproton obtained from a study of dijet events produced in association with a leading antiproton in $\bar pp$ collisions at $\sqrt s=630$ GeV at the Fermilab Tevatron. The ratio of $F_{jj}^D$ at $\sqrt s=630$ GeV to $F_{jj}^D$ obtained from a similar measurement at $\sqrt s=1800$ GeV is compared with expectations from QCD factorization and with theoretical predictions. We also report a measurement of the $\xi$ ($x$-Pomeron) and $\beta$ ($x$ of parton in Pomeron) dependence of $F_{jj}^D$ at $\sqrt s=1800$ GeV. In the region $0.035<\xi<0.095$, $|t|<1$ GeV$^2$ and $\beta<0.5$, $F_{jj}^D(\beta,\xi)$ is found to be of the form $\beta^{-1.0\pm 0.1} \xi^{-0.9\pm 0.1}$, which obeys $\beta$-$\xi$ factorization.Comment: LaTeX, 9 pages, Submitted to Phys. Rev. Letter

A Study of B0 -> J/psi K(*)0 pi+ pi- Decays with the Collider Detector at Fermilab

We report a study of the decays B0 -> J/psi K(*)0 pi+ pi-, which involve the creation of a u u-bar or d d-bar quark pair in addition to a b-bar -> c-bar(c s-bar) decay. The data sample consists of 110 1/pb of p p-bar collisions at sqrt{s} = 1.8 TeV collected by the CDF detector at the Fermilab Tevatron collider during 1992-1995. We measure the branching ratios to be BR(B0 -> J/psi K*0 pi+ pi-) = (8.0 +- 2.2 +- 1.5) * 10^{-4} and BR(B0 -> J/psi K0 pi+ pi-) = (1.1 +- 0.4 +- 0.2) * 10^{-3}. Contributions to these decays are seen from psi(2S) K(*)0, J/psi K0 rho0, J/psi K*+ pi-, and J/psi K1(1270)

Search for Kaluza-Klein Graviton Emission in ppˉp\bar{p} Collisions at s=1.8\sqrt{s}=1.8 TeV using the Missing Energy Signature

We report on a search for direct Kaluza-Klein graviton production in a data sample of 84 ${pb}^{-1}$ of \ppb collisions at $\sqrt{s}$ = 1.8 TeV, recorded by the Collider Detector at Fermilab. We investigate the final state of large missing transverse energy and one or two high energy jets. We compare the data with the predictions from a $3+1+n$-dimensional Kaluza-Klein scenario in which gravity becomes strong at the TeV scale. At 95% confidence level (C.L.) for $n$=2, 4, and 6 we exclude an effective Planck scale below 1.0, 0.77, and 0.71 TeV, respectively.Comment: Submitted to PRL, 7 pages 4 figures/Revision includes 5 figure

Measurement of the average time-integrated mixing probability of b-flavored hadrons produced at the Tevatron

We have measured the number of like-sign (LS) and opposite-sign (OS) lepton pairs arising from double semileptonic decays of $b$ and $\bar{b}$-hadrons, pair-produced at the Fermilab Tevatron collider. The data samples were collected with the Collider Detector at Fermilab (CDF) during the 1992-1995 collider run by triggering on the existence of $\mu \mu$ and $e \mu$ candidates in an event. The observed ratio of LS to OS dileptons leads to a measurement of the average time-integrated mixing probability of all produced $b$-flavored hadrons which decay weakly, $\bar{\chi} = 0.152 \pm 0.007$ (stat.) $\pm 0.011$ (syst.), that is significantly larger than the world average $\bar{\chi} = 0.118 \pm 0.005$.Comment: 47 pages, 10 figures, 15 tables Submitted to Phys. Rev.

Photoperiod Regulates Lean Mass Accretion, but Not Adiposity, in Growing F344 Rats Fed a High Fat Diet

yesIn this study the effects of photoperiod and diet, and their interaction, were examined for their effects on growth and body composition in juvenile F344 rats over a 4-week period. On long (16L:8D), relative to short (8L:16D), photoperiod food intake and growth rate were increased, but percentage adiposity remained constant (ca 3-4%). On a high fat diet (HFD), containing 22.8% fat (45% energy as fat), food intake was reduced, but energy intake increased on both photoperiods. This led to a small increase in adiposity (up to 10%) without overt change in body weight. These changes were also reflected in plasma leptin and lipid levels. Importantly while both lean and adipose tissue were strongly regulated by photoperiod on a chow diet, this regulation was lost for adipose, but not lean tissue, on HFD. This implies that a primary effect of photoperiod is the regulation of growth and lean mass accretion. Consistent with this both hypothalamic GHRH gene expression and serum IGF-1 levels were photoperiod dependent. As for other animals and humans, there was evidence of central hyposomatotropism in response to obesity, as GHRH gene expression was suppressed by the HFD. Gene expression of hypothalamic AgRP and CRH, but not NPY nor POMC, accorded with the energy balance status on long and short photoperiod. However, there was a general dissociation between plasma leptin levels and expression of these hypothalamic energy balance genes. Similarly there was no interaction between the HFD and photoperiod at the level of the genes involved in thyroid hormone metabolism (Dio2, Dio3, TSHβ or NMU), which are important mediators of the photoperiodic response. These data suggest that photoperiod and HFD influence body weight and body composition through independent mechanisms but in each case the role of the hypothalamic energy balance genes is not predictable based on their known function.Scottish Government (Rural and Environment Science and Analytical Services Division, http://www.scotland.gov.uk/), AWR LR LMT PJM and the BBSRC, (http://www.bbsrc.ac.uk/home/home.aspx, grant BB/K001043/1), AWR GH PJ

Vitamins A, C and E and the risk of breast cancer: results from a case-control study in Greece

Although several dietary compounds are hypothesized to have anticarcinogenic properties, the role ofpecific micronutrients in the development of breast cancer remains unclear. To address this issue, we assessed intake of retinol, β-carotene, vitamin C and vitamin E in relation to breast cancer risk in a case–control study in Greece. Eight hundrednd twenty women with histologically confirmed breast cancer were compared with 1548 control women. Dietary data were collectedhrough a 115-item semiquantitative food frequency questionnaire. Data were modelled by logistic regression, with adjustment forotal energy intake and established breast cancer risk factors, as well as mutual adjustment among the micronutrients. Amongost-menopausal women, there was no association between any of the micronutrients evaluated and risk of breast cancer. Amongremenopausal women, β-carotene, vitamin C and vitamin E were each inversely associated with breast cancer risk, but afterutual adjustment among the three nutrients only β-carotene remained significant; the odds ratio (OR) for a one-quintilencrease in β-carotene intake was 0.84 (95% confidence interval 0.73–0.97). The inverse association observed with β-carotene intake, however, is slightly weaker than the association previously observed with vegetable intake in these data,aising the possibility that the observed β-carotene effect is accounted for by another component of vegetables. ©1999 Cancer Research Campaig

Cytotoxic and Pathogenic Properties of Klebsiella oxytoca Isolated from Laboratory Animals

Klebsiella oxytoca is an opportunistic pathogen implicated in various clinical diseases in animals and humans. Studies suggest that in humans K. oxytoca exerts its pathogenicity in part through a cytotoxin. However, cytotoxin production in animal isolates of K. oxytoca and its pathogenic properties have not been characterized. Furthermore, neither the identity of the toxin nor a complete repertoire of genes involved in K. oxytoca pathogenesis have been fully elucidated. Here, we showed that several animal isolates of K. oxytoca, including the clinical isolates, produced secreted products in bacterial culture supernatant that display cytotoxicity on HEp-2 and HeLa cells, indicating the ability to produce cytotoxin. Cytotoxin production appears to be regulated by the environment, and soy based product was found to have a strong toxin induction property. The toxin was identified, by liquid chromatography-mass spectrometry and NMR spectroscopy, as low molecular weight heat labile benzodiazepine, tilivalline, previously shown to cause cytotoxicity in several cell lines, including mouse L1210 leukemic cells. Genome sequencing and analyses of a cytotoxin positive K. oxytoca strain isolated from an abscess of a mouse, identified genes previously shown to promote pathogenesis in other enteric bacterial pathogens including ecotin, several genes encoding for type IV and type VI secretion systems, and proteins that show sequence similarity to known bacterial toxins including cholera toxin. To our knowledge, these results demonstrate for the first time, that animal isolates of K. oxytoca, produces a cytotoxin, and that cytotoxin production is under strict environmental regulation. We also confirmed tilivalline as the cytotoxin present in animal K. oxytoca strains. These findings, along with the discovery of a repertoire of genes with virulence potential, provide important insights into the pathogenesis of K. oxytoca. As a novel diagnostic tool, tilivalline may serve as a biomarker for K oxytoca-induced cytotoxicity in humans and animals through detection in various samples from food to diseased samples using LC-MS/MS. Induction of K. oxytoca cytotoxin by consumption of soy may be in part involved in the pathogenesis of gastrointestinal disease