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.

PAR evaluation of treated Class I extraction patients

Objective: To evaluate treatment changes and quality of finishing occlusion in Class I patients treated with four premolar extractions. Material and Methods: Dental casts of 94 subjects (50 males and 44 females) were evaluated. Mean pretreatment age was 13.46 years, and mean treatment time was 2.09 years. The peer assessment rating (PAR) index was obtained from pretreatment and posttreatment dental casts. Results: The mean pretreatment PAR index of 29.46 was reduced to 6.32 at posttreatment stage, achieving a reduction of 78.54% with treatment. There was correlation between the initial PAR and correction during treatment, that is, the more severe the malocclusion the greater the treatment changes. Conclusion: The cases evaluated showed a high-standard orthodontic finishing.

The metal content of bulge field stars from FLAMES-GIRAFFE spectra - I. Stellar parameters and iron abundances

Aims. We determine the iron distribution function (IDF) for bulge field stars, in three different fields along the Galactic minor axis and at latitudes b = -4 degrees, b = -6 degrees, and b = -12 degrees. A fourth field including NGC 6553 is also included in the discussion. Methods. About 800 bulge field K giants were observed with the GIRAFFE spectrograph of FLAMES@VLT at spectral resolution R similar to 20 000. Several of them were observed again with UVES at R similar to 45 000 to insure the accuracy of the measurements. The LTE abundance analysis yielded stellar parameters and iron abundances that allowed us to construct an IDF for the bulge that, for the first time, is based on high-resolution spectroscopy for each individual star. Results. The IDF derived here is centered on solar metallicity, and extends from [Fe/H] similar to -1.5 to [Fe/H] similar to + 0.5. The distribution is asymmetric, with a sharper cutoff on the high-metallicity side, and it is narrower than previously measured. A variation in the mean metallicity along the bulge minor axis is clearly between b = -4 degrees and b = -6 degrees ([Fe/H] decreasing similar to by 0.6 dex per kpc). The field at b = -12 degrees. is consistent with the presence of a gradient, but its quantification is complicated by the higher disk/bulge fraction in this field. Conclusions. Our findings support a scenario in which both infall and outflow were important during the bulge formation, and then suggest the presence of a radial gradient, which poses some challenges to the scenario in which the bulge would result solely from the vertical heating of the bar.

PARALLAXES OF SOUTHERN EXTREMELY COOL OBJECTS. I. TARGETS, PROPER MOTIONS, AND FIRST RESULTS

We present results from the PARallaxes of Southern Extremely Cool objects ( PARSEC) program, an observational program begun in 2007 April to determine parallaxes for 122 L and 28 T southern hemisphere dwarfs using the Wide Field Imager on the ESO 2.2 m telescope. The results presented here include parallaxes of 10 targets from observations over 18 months and a first version proper motion catalog. The proper motions were obtained by combining PARSEC observations astrometrically reduced with respect to the Second US Naval Observatory CCD Astrograph Catalog, and the Two Micron All Sky Survey Point Source Catalog. The resulting median proper motion precision is 5 mas yr(-1) for 195,700 sources. The 140 0.3 deg(2) fields sample the southern hemisphere in an unbiased fashion with the exception of the galactic plane due to the small number of targets in that region. The proper motion distributions are shown to be statistically well behaved. External comparisons are also fully consistent. We will continue to update this catalog until the end of the program, and we plan to improve it including also observations from the GSC2.3 database. We present preliminary parallaxes with a 4.2 mas median precision for 10 brown dwarfs, two of which are within 10 pc. These increase the present number of L dwarfs by 20% with published parallaxes. Of the 10 targets, seven have been previously discussed in the literature: two were thought to be binary, but the PARSEC observations show them to be single; one has been confirmed as a binary companion and another has been found to be part of a binary system, both of which will make good benchmark systems. These results confirm that the foreseen precision of PARSEC can be achieved and that the large field of view will allow us to identify wide binary systems. Observations for the PARSEC program will end in early 2011 providing three to four years of coverage for all targets. The main expected outputs are: more than a 100% increase in the number of L dwarfs with parallaxes, increment in the number of objects per spectral subclass up to L9-in conjunction with published results-to at least 10, and to put sensible limits on the general binary fraction of brown dwarfs. We aim to contribute significantly to the understanding of the faint end of the H-R diagram and of the L/T transition region.

Star formation history of Canis Major R1 I. Wide-Field X-ray study of the young stellar population

Aims. The CMa R1 star-forming region contains several compact clusters as well as many young early-B stars. It is associated with a well-known bright rimmed nebula, the nature of which is unclear (fossil HII region or supernova remnant). To help elucidate the nature of the nebula, our goal was to reconstruct the star-formation history of the CMa R1 region, including the previously unknown older, fainter low-mass stellar population, using X-rays. Methods. We analyzed images obtained with the ROSAT satellite, covering similar to 5 sq. deg. Complementary VRI photometry was performed with the Gemini South telescope. Colour-magnitude and colour-colour diagrams were used in conjunction with pre-main sequence evolutionary tracks to derive the masses and ages of the X-ray sources. Results. The ROSAT images show two distinct clusters. One is associated with the known optical clusters near Z CMa, to which similar to 40 members are added. The other, which we name the ""GU CMa"" cluster, is new, and contains similar to 60 members. The ROSAT sources are young stars with masses down to M(star) similar to 0.5 M(circle dot), and ages up to 10 Myr. The mass functions of the two clusters are similar, but the GU CMa cluster is older than the cluster around Z CMa by at least a few Myr. Also, the GU CMa cluster is away from any molecular cloud, implying that star formation must have ceased; on the contrary (as already known), star formation is very active in the Z CMa region.

The B0.5IVe CoRoT target HD 49330 I. Photometric analysis from CoRoT data

Context. Be stars undergo outbursts producing a circumstellar disk from the ejected material. The beating of non-radial pulsations has been put forward as a possible mechanism of ejection. Aims. We analyze the pulsational behavior of the early B0.5IVe star HD 49330 observed during the first CoRoT long run towards the Galactical anticenter (LRA1). This Be star is located close to the lower edge of the beta Cephei instability strip in the HR diagram and showed a 0.03 mag outburst during the CoRoT observations. It is thus an ideal case for testing the aforementioned hypothesis. Methods. We analyze the CoRoT light curve of HD 49330 using Fourier methods and non-linear least square fitting. Results. In this star, we find pulsation modes typical of beta Cep stars (p modes) and SPB stars (g modes) with amplitude variations along the run directly correlated with the outburst. These results provide new clues about the origin of the Be phenomenon as well as strong constraints on the seismic modelling of Be stars.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri I. Long-term monitoring observations

Context. Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red (V/R) ratio of their primary components. The underlying structural and dynamical variations of the disks are only partly understood. From observations of the bright Be-shell star. Tau, the possibly broadest and longest data set illustrating the prototype of this behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, and interferometric observations. Aims. The dense, long-time monitoring permits a better separation of repetitive and ephemeral variations. The broad wavelength coverage includes lines formed under different physical conditions, i.e. different locations in the disk, so that the dynamics can be probed throughout much of the disk. Polarimetry and interferometry constrain the spatial structure. All together, the objective is a better understand the dynamics and life cycle of decretion disks. Methods. Standard methods of data acquisition, reduction, and analysis were applied. Results. From 3 V/R cycles between 1997 and 2008, a mean cycle length in Ha of 1400-1430 days was derived. After each minimum in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase may make the strongest contribution to the variability in cycle length. There is no obvious connection between the V/R cycle and the 133-day orbital period of the not otherwise detected companion. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and variability of the equivalent width. They seem to form in separate regions of the disk, probably crossing the line of sight at different times. The interferometry has resolved the continuum and the line emission in Br gamma and HeI 2.06. The phasing of the Br gamma emission shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The plane of the disk is constant throughout the observed V/R cycles. The observations lay the foundation for the fully self-consistent, one-armed, disk-oscillation model developed in Paper II.

Mn, Cu, and Zn abundances in barium stars and their correlations with neutron capture elements

Barium stars are optimal sites for studying the correlations between the neutron-capture elements and other species that may be depleted or enhanced, because they act as neutron seeds or poisons during the operation of the s-process. These data are necessary to help constrain the modeling of the neutron-capture paths and explain the s-process abundance curve of the solar system. Chemical abundances for a large number of barium stars with different degrees of s-process excesses, masses, metallicities, and evolutionary states are a crucial step towards this goal. We present abundances of Mn, Cu, Zn, and various light and heavy elements for a sample of barium and normal giant stars, and present correlations between abundances contributed to different degrees by the weak-s, mains, and r-processes of neutron capture, between Fe-peak elements and heavy elements. Data from the literature are also considered in order to better study the abundance pattern of peculiar stars. The stellar spectra were observed with FEROS/ESO. The stellar atmospheric parameters of the eight barium giant stars and six normal giants that we analyzed lie in the range 4300 < T(eff)/K < 5300, -0.7 < [Fe/H] <= 0.12 and 1.5 <= log g < 2.9. Carbon and nitrogen abundances were derived by spectral synthesis of the molecular bands of C(2), CH, and CN. For all other elements we used the atomic lines to perform the spectral synthesis. A very large scatter was found mainly for the Mn abundances when data from the literature were considered. We found that [Zn/Fe] correlates well with the heavy element excesses, its abundance clearly increasing as the heavy element excesses increase, a trend not shown by the [Cu/Fe] and [Mn/Fe] ratios. Also, the ratios involving Mn, Cu, and Zn and heavy elements usually show an increasing trend toward higher metallicities. Our results suggest that a larger fraction of the Zn synthesis than of Cu is owed to massive stars, and that the contribution of the main-s process to the synthesis of both elements is small. We also conclude that Mn is mostly synthesized by SN Ia, and that a non-negligible fraction of the synthesis of Mn, Cu, and Zn is owed to the weak s-process.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

Nuclear incoherent photoproduction of pi(0) and eta from 4 to 12 GeV

The mechanism of incoherent pi(0) and eta photoproduction from complex nuclei is investigated from 4 to 12 GeV with an extended version of the multicollisional Monte Carlo (MCMC) intranuclear cascade model. The calculations take into account the elementary photoproduction amplitudes via a Regge model and the nuclear effects of photon shadowing, Pauli blocking, and meson-nucleus final-state interactions. The results for pi(0) photoproduction reproduced for the first time the magnitude and energy dependence of the measured rations sigma(gamma A)/sigma(gamma N) for several nuclei (Be, C, Al, Cu, Ag, and Pb) from a Cornell experiment. The results for eta photoproduction fitted the inelastic background in Cornell's yields remarkably well, which is clearly not isotropic as previously considered in Cornell's analysis. With this constraint for the background, the eta -> gamma gamma. decay width was extracted using the Primakoff method, combining Be and Cu data [Gamma(eta ->gamma gamma) = 0.476(62) keV] and using Be data only [Gamma(eta ->gamma gamma) = 0.512(90) keV]; where the errors are only statistical. These results are in sharp contrast (similar to 50-60%) with the value reported by the Cornell group [Gamma(eta ->gamma gamma). = 0.324(46) keV] and in line with the Particle Data Group average of 0.510(26) keV.

eta ->gamma gamma decay width via the Primakoff cross section

Incoherent eta photoproduction in nuclei is evaluated at forward angles within 4 to 9 GeV using a multiple scattering Monte Carlo cascade calculation with full eta-nucleus final-state interactions. The Primakoff, nuclear coherent and nuclear incoherent components of the cross sections fit remarkably well previous measurements for Be and Cu from Cornell, suggesting a destructive interference between the Coulomb and nuclear coherent amplitudes for Cu. The inelastic background of the data is consistently attributed to the nuclear incoherent part, which is clearly not isotropic as previously considered in Cornell's analysis. The respective Primakoff cross sections from Be and Cu give Gamma(eta ->gamma gamma)=0.476(62) keV, where the quoted error is only statistical. This result is consistent with the Particle Data Group average of 0.510(26) keV and in sharp contrast (similar to 50%) with the value of 0.324(46) keV obtained at Cornell.

Dihadron azimuthal correlations in Au+Au collisions at root s(NN)=200 GeV

Azimuthal angle (Delta phi) correlations are presented for a broad range of transverse momentum (0.4 < p(T) < 10 GeV/c) and centrality (0-92%) selections for charged hadrons from dijets in Au+Au collisions at root s(NN) = 200 GeV. With increasing p(T), the away-side Delta phi distribution evolves from a broad and relatively flat shape to a concave shape, then to a convex shape. Comparisons with p + p data suggest that the away-side distribution can be divided into a partially suppressed ""head"" region centered at Delta phi similar to pi, and an enhanced ""shoulder"" region centered at Delta phi similar to pi +/- 1.1. The p(T) spectrum for the associated hadrons in the head region softens toward central collisions. The spectral slope for the shoulder region is independent of centrality and trigger p(T). The properties of the near-side distributions are also modified relative to those in p + p collisions, reflected by the broadening of the jet shape in Delta phi and Delta eta, and an enhancement of the per-trigger yield. However, these modifications seem to be limited to p(T)less than or similar to 4 GeV/c, above which both the hadron pair shape and per-trigger yield become similar to p + p collisions. These observations suggest that both the away- and near-side distributions contain a jet fragmentation component which dominates for p(T) greater than or similar to 5 GeV/c and a medium-induced component which is important for p(T) less than or similar to 4 GeV/c. We also quantify the role of jets at intermediate and low p(T) through the yield of jet-induced pairs in comparison with binary scaled p + p pair yield. The yield of jet-induced pairs is suppressed at high pair proxy energy (sum of the p(T) magnitudes of the two hadrons) and is enhanced at low pair proxy energy. The former is consistent with jet quenching; the latter is consistent with the enhancement of soft hadron pairs due to transport of lost energy to lower p(T).

Scaling properties at freeze-out in relativistic heavy-ion collisions

Identified charged pion, kaon, and proton spectra are used to explore the system size dependence of bulk freeze-out properties in Cu + Cu collisions at root s(NN) = 200 and 62.4 GeV. The data are studied with hydrodynamically motivated blast-wave and statistical model frameworks in order to characterize the freeze-out properties of the system. The dependence of freeze-out parameters on beam energy and collision centrality is discussed. Using the existing results from Au + Au and pp collisions, the dependence of freeze-out parameters on the system size is also explored. This multidimensional systematic study furthers our understanding of the QCD phase diagram revealing the importance of the initial geometrical overlap of the colliding ions. The analysis of Cu + Cu collisions expands the system size dependence studies from Au + Au data with detailed measurements in the smaller system. The systematic trends of the bulk freeze-out properties of charged particles is studied with respect to the total charged particle multiplicity at midrapidity, exploring the influence of initial state effects.

Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions

Parity (P)-odd domains, corresponding to nontrivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in noncentral collisions. To study this effect, we investigate a three-particle mixed-harmonics azimuthal correlator which is a P-even observable, but directly sensitive to the charge-separation effect. We report measurements of this observable using the STAR detector in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62 GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators and discuss in detail possible contributions from other effects that are not related to P violation.