An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Anti-Anopheles darlingi saliva antibodies as marker of Plasmodium vivax infection and clinical immunity in the Brazilian Amazon

Background: Despite governmental and private efforts on providing malaria control, this disease continues to be a major health threat. Thus, innovative strategies are needed to reduce disease burden. The malaria vectors, through the injection of saliva into the host skin, play important role on disease transmission and may influence malaria morbidity. This study describes the humoral immune response against Anopheles (An.) darlingi saliva in volunteers from the Brazilian Amazon and addresses the association between levels of specific antibodies and clinical presentation of Plasmodium (P.) vivax infection. Methods: Adult volunteers from communities in the Rondonia State, Brazil, were screened in order to assess the presence of P. vivax infection by light microscopy and nested PCR. Non-infected volunteers and individuals with symptomatic or symptomless infection were randomly selected and plasma collected. An. darlingi salivary gland sonicates (SGS) were prepared and used to measure anti-saliva antibody levels. Plasma interleukin (IL)-10 and interferon (IFN)-gamma levels were also estimated and correlated to anti-SGS levels. Results: Individuals infected with P. vivax presented higher levels of anti-SGS than non-infected individuals and antibody levels could discriminate infection. Furthermore, anti-saliva antibody measurement was also useful to distinguish asymptomatic infection from non-infection, with a high likelihood ratio. Interestingly, individuals with asymptomatic parasitaemia presented higher titers of anti-SGS and lower IFN-gamma/IL-10 ratio than symptomatic ones. In P. vivax-infected asymptomatic individuals, the IFN-gamma/IL-10 ratio was inversely correlated to anti-SGS titers, although not for while in symptomatic volunteers. Conclusion: The estimation of anti-An. darlingi antibody levels can indicate the probable P. vivax infection status and also could serve as a marker of disease severity in this region of Brazilian Amazon.

CHARACTERIZATION OF EXPONENTIAL DIVERGENCE OF THE KALMAN FILTER FOR TIME-VARYING SYSTEMS

This paper studies semistability of the recursive Kalman filter in the context of linear time-varying (LTV), possibly nondetectable systems with incorrect noise information. Semistability is a key property, as it ensures that the actual estimation error does not diverge exponentially. We explore structural properties of the filter to obtain a necessary and sufficient condition for the filter to be semistable. The condition does not involve limiting gains nor the solution of Riccati equations, as they can be difficult to obtain numerically and may not exist. We also compare semistability with the notions of stability and stability w.r.t. the initial error covariance, and we show that semistability in a sense makes no distinction between persistent and nonpersistent incorrect noise models, as opposed to stability. In the linear time invariant scenario we obtain algebraic, easy to test conditions for semistability and stability, which complement results available in the context of detectable systems. Illustrative examples are included.

Measurement of the underground atmospheric muon charge ratio using the MINOS Near Detector

The magnetized MINOS Near Detector, at a depth of 225 mwe, is used to measure the atmospheric muon charge ratio. The ratio of observed positive to negative atmospheric muon rates, using 301 days of data, is measured to be 1.266 +/- 0.001(stat)(-0.014)(+0.015)(syst). This measurement is consistent with previous results from other shallow underground detectors and is 0.108 +/- 0.019(stat + syst) lower than the measurement at the functionally identical MINOS Far Detector at a depth of 2070 mwe. This increase in charge ratio as a function of depth is consistent with an increase in the fraction of muons arising from kaon decay for increasing muon surface energies.

Measurement of neutrino oscillations with the MINOS detectors in the NuMI beam

This Letter reports new results from the MINOS experiment based on a two-year exposure to muon neutrinos from the Fermilab NuMI beam. Our data are consistent with quantum-mechanical oscillations of neutrino flavor with mass splitting vertical bar Delta m(2)vertical bar = (2.43 +/- 0.13) x 10(-3) eV(2) (68% C.L.) and mixing angle sin(2)(2 theta) > 0.90 (90% C.L.). Our data disfavor two alternative explanations for the disappearance of neutrinos in flight: namely, neutrino decays into lighter particles and quantum decoherence of neutrinos, at the 3.7 and 5.7 standard-deviation levels, respectively.

Measurement of the Neutrino Mass Splitting and Flavor Mixing by MINOS

Measurements of neutrino oscillations using the disappearance of muon neutrinos from the Fermilab NuMI neutrino beam as observed by the two MINOS detectors are reported. New analysis methods have been applied to an enlarged data sample from an exposure of 7.25 x 10(20) protons on target. A fit to neutrino oscillations yields values of vertical bar Delta m(2)vertical bar = (2.32(-0.08)(+0.12) x 10(-3) eV(2) for the atmospheric mass splitting and sin(2)(2 theta) > 0.90 (90% C.L.) for the mixing angle. Pure neutrino decay and quantum decoherence hypotheses are excluded at 7 and 9 standard deviations, respectively.

Measurement of neutral mesons in p plus p collisions at root s=200 GeV and scaling properties of hadron production

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K(S)(0), omega, eta', and phi mesons in p + p collisions at root s 200 GeV. Measurements of omega and phi production in different decay channels give consistent results. New results for the omega are in agreement with previously published data and extend the measured p(T) coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-p(T) and characterizing the low-p(T) regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

Measurement of transverse single-spin asymmetries for J/psi production in polarized p plus p collisions at root s=200 GeV

We report the first measurement of transverse single-spin asymmetries in J/psi production from transversely polarized p + p collisions at root s = 200 GeV with data taken by the PHENIX experiment in 2006 and 2008. The measurement was performed over the rapidity ranges 1.2 < vertical bar y vertical bar < 2.2 and vertical bar y vertical bar < 0.35 for transverse momenta up to 6 GeV/c. J/psi production at the Relativistic Heavy Ion Collider is dominated by processes involving initial-state gluons, and transverse single-spin asymmetries of the J/psi can provide access to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing question in QCD of the J/psi production mechanism. Asymmetries were obtained as a function of J/psi transverse momentum and Feynman-x, with a value of -0.086 +/- 0.026(stat) +/- 0.003(syst) in the forward region. This result suggests possible nonzero trigluon correlation functions in transversely polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function.

High p(T) direct photon and pi(0) triggered azimuthal jet correlations and measurement of k(T) for isolated direct photons in p plus p collisions at root s=200 GeV

Correlations of charged hadrons of 1< p(T) < 10 Gev/c with high pT direct photons and pi(0) mesons in the range 5< p(T) < 15 Gev/c are used to study jet fragmentation in the gamma + jet and dijet channels, respectively. The magnitude of the partonic transverse momentum, k(T), is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

Detailed measurement of the e(+)e(-) pair continuum in p plus p and Au plus Au collisions at root s(NN)=200 GeV and implications for direct photon production

PHENIX has measured the e(+)e(-) pair continuum in root s(NN) = 200 GeV Au+Au and p+p collisions over a wide range of mass and transverse momenta. The e(+)e(-) yield is compared to the expectations from hadronic sources, based on PHENIX measurements. In the intermediate-mass region, between the masses of the phi and the J/psi meson, the yield is consistent with expectations from correlated c (c) over bar production, although other mechanisms are not ruled out. In the low-mass region, below the phi, the p+p inclusive mass spectrum is well described by known contributions from light meson decays. In contrast, the Au+Au minimum bias inclusive mass spectrum in this region shows an enhancement by a factor of 4.7 +/- 0.4(stat) +/- 1.5(syst) +/- 0.9(model). At low mass (m(ee) < 0.3 GeV/c(2)) and high p(T) (1 < p(T) < 5 GeV/c) an enhanced e(+)e(-) pair yield is observed that is consistent with production of virtual direct photons. This excess is used to infer the yield of real direct photons. In central Au+Au collisions, the excess of the direct photon yield over the p+p is exponential in p(T), with inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar or equal to 300-600 MeV at times of 0.6-0.15 fm/c after the collision are in qualitative agreement with the direct photon data in Au+Au. For low p(T) < 1 GeV/c the low-mass region shows a further significant enhancement that increases with centrality and has an inverse slope of T similar or equal to 100 MeV. Theoretical models underpredict the low-mass, low-p(T) enhancement.

Measurement of Bottom Versus Charm as a Function of Transverse Momentum with Electron-Hadron Correlations in p plus p Collisions at s=200 GeV

The momentum distribution of electrons from semileptonic decays of charm and bottom quarks for midrapidity |y|< 0.35 in p+p collisions at s=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range 2 < p(T)< 7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D -> e(+/-)K(-/+)X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in p(T). A fixed-order-plus-next-to-leading-log perturbative quantum chromodynamics calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is sigma(bb)=3.2(-1.1)(+1.2)(stat)(-1.3)(+1.4)(syst)mu b.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Direct Measurement of the Bose-Einstein Condensation Universality Class in NiCl(2)-4SC(NH(2))(2) at Ultralow Temperatures

In this work, we demonstrate field-induced Bose-Einstein condensation (BEC) in the organic compound NiCl(2)-4SC(NH(2))(2) using ac susceptibility measurements down to 1 mK. The Ni S=1 spins exhibit 3D XY antiferromagnetism between a lower critical field H(c1)similar to 2 T and a upper critical field H(c2)similar to 12 T. The results show a power-law temperature dependence of the phase transition line H(c1)(T)-H(c1)(0)=aT(alpha) with alpha=1.47 +/- 0.10 and H(c1)(0)=2.053 T, consistent with the 3D BEC universality class. Near H(c2), a kink was found in the phase boundary at approximately 150 mK.

Direct measurement of spin correlations using magnetostriction

We demonstrate that the short-range spin correlator < S(i)center dot S(j)>, a fundamental measure of the interaction between adjacent spins, can be directly measured in certain insulating magnets. We present magnetostriction data for the insulating organic compound NiCl(2)-4SC(NH(2))(2), and show that the magnetostriction as a function of field is proportional to the dominant short-range spin correlator. Furthermore, the constant of proportionality between the magnetostriction and the spin correlator gives information about the spin-lattice interaction. Combining these results with the measured Young's modulus, we are able to extract dJ/dz, the dependence of the superexchange constant J on the Ni interionic distance z.

Measurement of the Parity-Violating Longitudinal Single-Spin Asymmetry for W(+/-) Boson Production in Polarized Proton-Proton Collisions at root s=500 GeV

We report the first measurement of the parity-violating single-spin asymmetries for midrapidity decay positrons and electrons from W(+) and W(-) boson production in longitudinally polarized proton-proton collisions at root s = 500 GeV by the STAR experiment at RHIC. The measured asymmetries, A(L)(W+) = -0.27 +/- 0.10(stat.) +/- 0.02(syst.) +/- 0.03(norm.) and A(L)(W-) = 0.14 +/- 0.19(stat.) +/- 0.02(syst.) +/- 0.01(norm.), are consistent with theory predictions, which are large and of opposite sign. These predictions are based on polarized quark and antiquark distribution functions constrained by polarized deep-inelastic scattering measurements.