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.

Time variation of some selected topics in bioethical publications

Objective: To analyse the time variation of topics in bioethical publications as a proxy of the relative importance. Methods: We searched the Medline database for bioethics publications using the words ""ethics or bioethics'', and for 360 specific topics publications, associating Medical Subject Heading topic descriptors to those words. We calculated the ratio of bioethics publications to the total publications of Medline, and the ratio of each topic publications to the total bioethics publications, for five-year intervals, from 1970 to 2004. We calculated the time variation of ratios, dividing the difference between the highest and lowest ratio of each topic by its highest ratio. Four topics were described, selected to illustrate different patterns of variation: ""Induced Abortion'', ""Conflict of Interest'', ""Acquired Immunodeficiency Syndrome'', ""Medical Education.'' Results: The ratio of bioethics publications to total Medline publications increased from 0.003 to 0.012. The variation of the topic's ratios was higher than 0.7 for 68% of the topics. The Induced Abortion ratios decreased from 0.12 to 0.02. Conflict of Interest ratios increased from zero to 0.07. The Acquired Immunodeficiency Syndrome ratios were nearly zero in the first three intervals, had a peak of 0.06 during 1985-9, followed by a decrease to 0.01. Medical Education ratios varied few, from 0.04 to 0.03. Conclusions: There was an increase of bioethical publications in the Medline database. The topics in bioethics literature have an important time variation. Some factors were suggested to explain this variation: current legal cases, resolution of the issue, saturation of a discussion and epidemiologic importance.

Effect of a Single Application of TiF(4) and NaF Varnishes and Solutions Combined with Nd:YAG Laser Irradiation on Enamel Erosion in Vitro

Objective: This in vitro study aimed to analyze the influence of neodymium-doped yttrium aluminum garnet (Nd:YAG) laser irradiation on the efficacy of titanium tetrafluoride (TiF(4)) and sodium fluoride (NaF) varnishes and solutions to protect enamel against erosion. Background data: The effect of Nd:YAG laser irradiation on NaF and AmF was analyzed; however, there is no available data on the interaction between Nd:YAG laser irradiation and TiF(4). Methods: Bovine enamel specimens were pre-treated with NaF varnish, TiF(4) varnish, NaF solution, TiF(4) solution, placebo varnish, Nd:YAG (84.9 J/cm(2)), Nd:YAG prior to or through NaF varnish, Nd:YAG prior to or through TiF(4) varnish, Nd:YAG prior to or through NaF solution, Nd:YAG prior to or through TiF(4) solution, and Nd:YAG prior to or through placebo varnish. Controls remained untreated. Ten specimens in each group were then subjected to an erosive demineralization (Sprite Zero, 4x90 s/day) and remineralization (artificial saliva, between the erosive cycles) cycling for 5 days. Enamel loss was measured profilometrically (mu m). Additionally, treated but non-eroded specimens were additionally analyzed by scanning electron microscope (SEM) (each group n-2). The data were statistically analyzed by ANOVA and Tukey's post-hoc test (p < 0.05). Results: Only TiF(4) varnish (1.8 +/- 0.6 mu m), laser prior to TiF(4) varnish (1.7 +/- 0.3 mu m) and laser prior to TiF(4) solution (1.4 +/- 0.3 mu m) significantly reduced enamel erosion compared to the control (4.1 +/- 0.6 mu m). SEM pictures showed that specimens treated with TiF(4) varnish presented a surface coating. Conclusions: Nd:YAG laser irradiation was not effective against enamel erosion and it did not have any influence on the efficacy of F, except for TiF(4) solution. On the other hand, TiF(4) varnish protected against enamel erosion, without the influence of laser irradiation.

On the mass determination of super-Earths orbiting active stars: the CoRoT-7 system

Context. Star activity makes the mass determination of CoRoT-7b and CoRoT 7c uncertain. Investigators of the CoRoT team proposed several solutions, but all but one of them are larger than the initial determinations of 4.8 +/- 0.8 M(Earth) for CoRoT-7b and 8.4 +/- 0.9 M(Earth) for CoRoT 7c. Aims. This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to redetermine the planet masses and to explore techniques for determining mass and orbital elements of planets discovered around active stars when the relative variation in the radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. Methods. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009, A&A, 506, 303) in the first mass determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) the approach proposed by Hatzes et al. (2010, A&A, 520, A93) using only those nights in which two or three observations were done; (2) a pure Fourier analysis. In all cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system. The periods are also kept fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2 454 847-873 are used because they include many nights with multiple observations; otherwise, it is not possible to separate the effects of the rotation fourth harmonic (5.91 d = P(rot)/4) from the alias of the orbital period of CoRoT-7b (0.853585 d). Results. The results of the various approaches are combined to give planet mass values 8.0 +/- 1.2 M(Earth) for CoRoT-7b and 13.6 +/- 1.4 M(Earth) for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given. Conclusions. The results obtained with three different approaches agree to give higher masses than those in previous determinations. From the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 11 +/- 3.5 g cm(-3), so CoRoT-7b may be rocky with a large iron core.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Kinematics of galaxies in compact groups Studying the B-band Tully-Fischer relation

We obtained new Fabry-Perot data cubes and derived velocity fields, monochromatic, and velocity dispersion maps for 28 galaxies in the Hickson compact groups 37, 40, 47, 49, 54, 56, 68, 79, and 93. We also derived rotation curves for 9 of the studied galaxies, 6 of which are strongly asymmetric. Combining these new data with previously published 2D kinematic maps of compact group galaxies, we investigated the differences between the kinematic and morphological position angles for a sample of 46 galaxies. We find that one third of the unbarred compact group galaxies have position angle misalignments between the stellar and gaseous components. This and the asymmetric rotation curves are clear signatures of kinematic perturbations, probably because of interactions among compact group galaxies. A comparison between the B-band Tully-Fisher relation for compact group galaxies and for the GHASP field-galaxy sample shows that, despite the high fraction of compact group galaxies with asymmetric rotation curves, these lay on the TF relation defined by galaxies in less dense environments, although with more scatter. This agrees with previous results, but now confirmed for a larger sample of 41 galaxies. We confirm the tendency for compact group galaxies at the low-mass end of the Tully-Fisher relation (HCG 49b, 89d, 96c, 96d, and 100c) to have either a magnitude that is too bright for its mass (suggesting brightening by star formation) and/or a low maximum rotational velocity for its luminosity (suggesting tidal stripping). These galaxies are outside the Tully Fisher relation at the 1 sigma level, even when the minimum acceptable values of inclinations are used to compute their maximum velocities. Including such galaxies with nu < 100 km s(-1) in the determination of the zero point and slope of the compact group B-band Tully-Fisher relation would strongly change the fit, making it different from the relation for field galaxies, which has to be kept in mind when studying scaling relations of interacting galaxies, especially at high redshifts.

Spiral-like structure at the centre of nearby clusters of galaxies

Context. X-ray data analysis have found that fairly complex structures at cluster centres are more common than expected. Many of these structures have similar morphologies, which exhibit spiral-like substructure. Aims. It is not yet well known how these structures are formed or maintained. Understanding the origin of these spiral-like features at the centre of some clusters is the major motivation behind this work. Methods. We analyse deep Chandra observations of 15 nearby galaxy clusters ( 0.01 < z < 0.06), and use X-ray temperature and substructure maps to detect small features at the cores of the clusters. Results. We detect spiral-like features at the centre of 7 clusters: A85, A426, A496, Hydra A cluster, Centaurus, Ophiuchus, and A4059. These patterns are similar to those found in numerical hydrodynamic simulations of cluster mergers with non-zero impact parameter. In some clusters of our sample, a strong radio source also occupies the inner region of the cluster, which indicates a possible connection between the two. Our investigation implies that these spiral-like structures may be caused by off-axis minor mergers. Since these features occur in regions of high density, they may confine radio emission from the central galaxy producing, in some cases, unusual radio morphology.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

A comparison of dry and wet season aerosol number fluxes over the Amazon rain forest

Vertical number fluxes of aerosol particles and vertical fluxes of CO(2) were measured with the eddy covariance method at the top of a 53 m high tower in the Amazon rain forest as part of the LBA (The Large Scale Biosphere Atmosphere Experiment in Amazonia) experiment. The observed aerosol number fluxes included particles with sizes down to 10 nm in diameter. The measurements were carried out during the wet and dry season in 2008. In this study focus is on the dry season aerosol fluxes, with significant influence from biomass burning, and these are compared with aerosol fluxes measured during the wet season. Net particle deposition fluxes dominated in daytime in both seasons and the deposition flux was considerably larger in the dry season due to the much higher dry season particle concentration. The particle transfer velocity increased linearly with increasing friction velocity in both seasons. The difference in transfer velocity between the two seasons was small, indicating that the seasonal change in aerosol number size distribution is not enough for causing any significant change in deposition velocity. In general, particle transfer velocities in this study are low compared to studies over boreal forests. The reasons are probably the high percentage of accumulation mode particles and the low percentage of nucleation mode particles in the Amazon boundary layer, both in the dry and wet season, and low wind speeds in the tropics compared to the midlatitudes. In the dry season, nocturnal particle fluxes behaved very similar to the nocturnal CO(2) fluxes. Throughout the night, the measured particle flux at the top of the tower was close to zero, but early in the morning there was an upward particle flux peak that is not likely a result of entrainment or local pollution. It is possible that these morning upward particle fluxes are associated with emission of primary biogenic particles from the rain forest. Emitted particles may be stored within the canopy during stable conditions at nighttime, similarly to CO(2), and being released from the canopy when conditions become more turbulent in the morning.

Transverse momentum dependence of eta meson suppression in Au plus Au collisions at root s(NN)=200 GeV

New measurements by the PHENIX experiment at the Relativistic Heavy Ion Collider for. production at midrapidity as a function of transverse momentum ((PT)) and collision centrality in root s(NN) = 200 GeV Au + Au and p + p collisions are presented. They indicate nuclear modification factors (R(AA)) which are similar in both magnitude and trend to those found in earlier pi(0) measurements. Linear fits to R(AA) as a function of (PT) in 5-20 GeV/c show that the slope is consistent with zero within two standard deviations at all centralities, although a slow rise cannot be excluded. Having different statistical and systematic uncertainties, the pi(0) and eta measurements are complementary at high (PT); thus, along with the extended (PT) range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

Transverse momentum dependence of J/psi polarization at midrapidity in p plus p collisions at root s=200 GeV

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p + p collisions at root s = 200 GeV performed by the PHENIX Experiment at the Relativistic Heavy Ion Collider. The J/psi polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p(T) < 5 GeV/c and vertical bar y vertical bar < 0.35. The polarization in the helicity and Gottfried-Jackson frames is consistent with zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance in the Collins-Soper frame and the uncertainties of the current data. The results are compared to observations for other collision systems and center of mass energies and to different quarkonia production models.

Charged Kaon Interferometric Probes of Space-Time Evolution in Au plus Au Collisions at s(NN)=200 GeV

Bose-Einstein correlations of charged kaons are used to probe Au+Au collisions at s(NN)=200 GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. Three-dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission source function. The centrality dependences of the three Gaussian radii are well described by a single linear function of N(part)(1/3) with a zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r greater than or similar to 10 fm, although the bulk emission at lower radius is well described by a Gaussian. The presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a heavy-ion collision is extended, and that similar measurements with pions are not solely due to the decay of long-lived resonances.

Lepton sector of a fourth generation

In extensions of the standard model with a heavy fourth generation, one important question is what makes the fourth-generation lepton sector, particularly the neutrinos, so different from the lighter three generations. We study this question in the context of models of electroweak symmetry breaking in warped extra dimensions, where the flavor hierarchy is generated by choosing the localization of the zero-mode fermions in the extra dimension. In this setup the Higgs sector is localized near the infrared brane, whereas the Majorana mass term is localized at the ultraviolet brane. As a result, light neutrinos are almost entirely Majorana particles, whereas the fourth-generation neutrino is mostly a Dirac fermion. We show that it is possible to obtain heavy fourth-generation leptons in regions of parameter space where the light neutrino masses and mixings are compatible with observation. We study the impact of these bounds, as well as the ones from lepton flavor violation, on the phenomenology of these models.

Temperature-driven collapse of a coherent binary mixture of condensates

We present a temperature- dependent Hartree- Fock- Bogoliubov- Popov theory to analyze the properties of the equilibrium states of an homogeneous mixture of bosonic atoms in two different hyperfine states and in the presence of an internal Josephson coupling. In our calculation we show that the bistable structure of the equilibrium states at zero temperature changes when we increase the temperature of the system. We investigate two mechanisms of the disappearance of bistability. In one, near the collapse of one of the equilibrium states, the acoustical branch becomes unstable and the gap of the optical branch goes to zero. In the other, there is no divergent behavior of the system and bistability disappears at a temperature in which the two equilibrium states merge at a zero- population fraction imbalance. When we further increase the temperature, this state remains as a unique equilibrium configuration.