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.

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.

Relationship Between Surface Topography and Energy Density Distribution of Er,Cr:YSGG Beam on Irradiated Dentin: An Atomic Force Microscopy Study

Objective: The aim of this study was to assess by atomic force microscopy (AFM) the effect of Er,Cr:YSGG laser application on the surface microtopography of radicular dentin. Background: Lasers have been used for various purposes in dentistry, where they are clinically effective when used in an appropriate manner. The Er, Cr: YSGG laser can be used for caries prevention when settings are below the ablation threshold. Materials and Methods: Four specimens of bovine dentin were irradiated using an Er, Cr:YSGG laser (lambda = 2.78 mu m), at a repetition rate of 20 Hz, with a 750-mu m-diameter sapphire tip and energy density of 2.8 J/cm(2) (12.5 mJ/pulse). After irradiation, surface topography was analyzed by AFM using a Si probe in tapping mode. Quantitative and qualitative information concerning the arithmetic average roughness (Ra) and power spectral density analyses were obtained from central, intermediate, and peripheral areas of laser pulses and compared with data from nonirradiated samples. Results: Dentin Ra for different areas were as follows: central, 261.26 (+/- 21.65) nm; intermediate, 83.48 (+/- 6.34) nm; peripheral, 45.8 (+/- 13.47) nm; and nonirradiated, 35.18 (+/- 2.9) nm. The central region of laser pulses presented higher ablation of intertubular dentin, with about 340-760 nm height, while intermediate, peripheral, and nonirradiated regions presented no difference in height of peritubular and interperitubular dentin. Conclusion: According to these results, we can assume that even when used at a low-energy density parameter, Er, Cr: YSGG laser can significantly alter the microtopography of radicular dentin, which is an important characteristic to be considered when laser is used for clinical applications.

Gold-implanted shallow conducting layers in polymethylmethacrylate

PMMA (polymethylmethacrylate) was ion implanted with gold at very low energy and over a range of different doses using a filtered cathodic arc metal plasma system. A nanometer scale conducting layer was formed, fully buried below the polymer surface at low implantation dose, and evolving to include a gold surface layer as the dose was increased. Depth profiles of the implanted material were calculated using the Dynamic TRIM computer simulation program. The electrical conductivity of the gold-implanted PMMA was measured in situ as a function of dose. Samples formed at a number of different doses were subsequently characterized by Rutherford backscattering spectrometry, and test patterns were formed on the polymer by electron beam lithography. Lithographic patterns were imaged by atomic force microscopy and demonstrated that the contrast properties of the lithography were well maintained in the surface-modified PMMA.

Conducting polymer formed by low energy gold ion implantation

A buried conducting layer of metal/polymer nanocomposite was formed by very low energy gold ion implantation into polymethylmethacrylate. The conducting layer is similar to 3 nm deep and of width similar to 1 nm. In situ resistivity measurements were performed as the implantation proceeded, and the conductivity thus obtained as a function of buried gold concentration. The measured conductivity obeys the behavior well established for composites in the percolation regime. The critical concentration, below which the polymer remains an insulator, is attained at a dose similar to 1.0 x 10(16) atoms/cm(2) of implanted gold ions. (C) 2008 American Institute of Physics.

Neutrino and antineutrino inclusive charged-current cross section measurements with the MINOS near detector

The energy dependence of the neutrino-iron and antineutrino-iron inclusive charged-current cross sections and their ratio have been measured using a high-statistics sample with the MINOS near detector exposed to the NuMI beam from the main injector at Fermilab. Neutrino and antineutrino fluxes were determined using a low hadronic energy subsample of charged-current events. We report measurements of nu-Fe ((nu) over bar - Fe) cross section in the energy range 3-50 GeV (5-50 GeV) with precision of 2%-8% (3%-9%) and their ratio which is measured with precision 2%-8%. The data set spans the region from low energy, where accurate measurements are sparse, up to the high-energy scaling region where the cross section is well understood.

Search for sterile neutrino mixing in the MINOS long-baseline experiment

A search for depletion of the combined flux of active neutrino species over a 735 km baseline is reported using neutral-current interaction data recorded by the MINOS detectors in the NuMI neutrino beam. Such a depletion is not expected according to conventional interpretations of neutrino oscillation data involving the three known neutrino flavors. A depletion would be a signature of oscillations or decay to postulated noninteracting sterile neutrinos, scenarios not ruled out by existing data. From an exposure of 3.18 x 10(20) protons on target in which neutrinos of energies between similar to 500 MeV and 120 GeV are produced predominantly as nu(mu), the visible energy spectrum of candidate neutral-current reactions in the MINOS far detector is reconstructed. Comparison of this spectrum to that inferred from a similarly selected near-detector sample shows that of the portion of the nu(mu) flux observed to disappear in charged-current interaction data, the fraction that could be converting to a sterile state is less than 52% at 90% confidence level (C. L.). The hypothesis that active neutrinos mix with a single sterile neutrino via oscillations is tested by fitting the data to various models. In the particular four-neutrino models considered, the mixing angles theta(24) and theta(34) are constrained to be less than 11 degrees and 56 degrees at 90% C. L., respectively. The possibility that active neutrinos may decay to sterile neutrinos is also investigated. Pure neutrino decay without oscillations is ruled out at 5.4 standard deviations. For the scenario in which active neutrinos decay into sterile states concurrently with neutrino oscillations, a lower limit is established for the neutrino decay lifetime tau(3)/m(3) > 2.1 x 10(-12) s/eV at 90% C.L.

Search for Muon-Neutrino to Electron-Neutrino Transitions in MINOS

This Letter reports on a search for nu(mu)->nu(e) transitions by the MINOS experiment based on a 3.14x10(20) protons-on-target exposure in the Fermilab NuMI beam. We observe 35 events in the Far Detector with a background of 27 +/- 5(stat)+/- 2(syst) events predicted by the measurements in the Near Detector. If interpreted in terms of nu(mu)->nu(e) oscillations, this 1.5 sigma excess of events is consistent with sin(2)(2 theta(13)) comparable to the CHOOZ limit when |Delta m(2)|=2.43x10(-3) eV(2) and sin(2)(2 theta(23))=1.0 are assumed.