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.

Impact of outdoor biomass air pollution on hypertension hospital admissions

Background This study aimed to evaluate the association between the total suspended particles (TSP) generated from burning sugar cane plantations and the incidence of hospital admissions from hypertension in the city of Araraquara. Methods The study was an ecological time-series study. Total daily records of hypertension (ICD 10th I10-15) were obtained from admitted patients of all ages in a hospital in Araraquara, Sao Paulo State, Brazil, from 23 March 2003 to 27 July 2004. The daily concentration of TSP (mu g/m(3)) was obtained using a Handi-Vol sampler placed in downtown Araraquara. The local airport provided daily measures of temperature and humidity. In generalised linear Poisson regression models, the daily number of hospital admissions for hypertension was considered to be the dependent variable and the daily TSP concentration the independent variable. Results TSP presented a lagged effect on hypertension admissions, which was first observed 1 day after a TSP increase and remained almost unchanged for the following 2 days. A 10 mu g/m(3) increase in the TSP 3 day moving average lagged in 1 day led to an increase in hypertension-related hospital admissions during the harvest period (12.5%, 95% CI 5.6% to 19.9%) that was almost 30% higher than during non-harvest periods (9.0%, 95% CI 4.0% to 14.3%). Conclusions Increases in TSP concentrations were associated with hypertension-related hospital admissions. Despite the benefits of reduced air pollution in urban cities achieved by using ethanol produced from sugar cane to power automobiles, areas where the sugar cane is produced and harvested were found to have increased public health risk.

A Novel Hepatitis C Virus Genotyping Method Based on Liquid Microarray

The strategy used to treat HCV infection depends on the genotype involved. An accurate and reliable genotyping method is therefore of paramount importance. We describe here, for the first time, the use of a liquid microarray for HCV genotyping. This liquid microarray is based on the 5'UTR - the most highly conserved region of HCV - and the variable region NS5B sequence. The simultaneous genotyping of two regions can be used to confirm findings and should detect inter-genotypic recombination. Plasma samples from 78 patients infected with viruses with genotypes and subtypes determined in the Versant (TM) HCV Genotype Assay LiPA (version I; Siemens Medical Solutions, Diagnostics Division, Fernwald, Germany) were tested with our new liquid microarray method. This method successfully determined the genotypes of 74 of the 78 samples previously genotyped in the Versant (TM) HCV Genotype Assay LiPA (74/78, 95%). The concordance between the two methods was 100% for genotype determination (74/74). At the subtype level, all 3a and 2b samples gave identical results with both methods (17/17 and 7/7, respectively). Two 2c samples were correctly identified by microarray, but could only be determined to the genotype level with the Versant (TM) HCV assay. Genotype ""1'' subtypes (1a and 1b) were correctly identified by the Versant (TM) HCV assay and the microarray in 68% and 40% of cases, respectively. No genotype discordance was found for any sample. HCV was successfully genotyped with both methods, and this is of prime importance for treatment planning. Liquid microarray assays may therefore be added to the list of methods suitable for HCV genotyping. It provides comparable results and may readily be adapted for the detection of other viruses frequently co-infecting HCV patients. Liquid array technology is thus a reliable and promising platform for HCV genotyping.

Ultra-Deep Sequencing Reveals the microRNA Expression Pattern of the Human Stomach

Background: While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia. Methodology/Principal Findings: A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue. Conclusions/Significance: This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

Psychometric properties of the Child Health Assessment Questionnaire (CHAQ) applied to children and adolescents with cerebral palsy

Background: Cerebral palsy (CP) patients have motor limitations that can affect functionality and abilities for activities of daily living (ADL). Health related quality of life and health status instruments validated to be applied to these patients do not directly approach the concepts of functionality or ADL. The Child Health Assessment Questionnaire (CHAQ) seems to be a good instrument to approach this dimension, but it was never used for CP patients. The purpose of the study was to verify the psychometric properties of CHAQ applied to children and adolescents with CP. Methods: Parents or guardians of children and adolescents with CP, aged 5 to 18 years, answered the CHAQ. A healthy group of 314 children and adolescents was recruited during the validation of the CHAQ Brazilian-version. Data quality, reliability and validity were studied. The motor function was evaluated by the Gross Motor Function Measure (GMFM). Results: Ninety-six parents/guardians answered the questionnaire. The age of the patients ranged from 5 to 17.9 years (average: 9.3). The rate of missing data was low(< 9.3%). The floor effect was observed in two domains, being higher only in the visual analogue scales (<= 35.5%). The ceiling effect was significant in all domains and particularly high in patients with quadriplegia (81.8 to 90.9%) and extrapyramidal (45.4 to 91.0%). The Cronbach alpha coefficient ranged from 0.85 to 0.95. The validity was appropriate: for the discriminant validity the correlation of the disability index with the visual analogue scales was not significant; for the convergent validity CHAQ disability index had a strong correlation with the GMFM (0.77); for the divergent validity there was no correlation between GMFM and the pain and overall evaluation scales; for the criterion validity GMFM as well as CHAQ detected differences in the scores among the clinical type of CP (p < 0.01); for the construct validity, the patients' disability index score (mean: 2.16; SD: 0.72) was higher than the healthy group ( mean: 0.12; SD: 0.23)(p < 0.01). Conclusion: CHAQ reliability and validity were adequate to this population. However, further studies are necessary to verify the influence of the ceiling effect on the responsiveness of the instrument.

Ablation Rate and Morphology of Superficial and Deep Dentin Irradiated with Different Er:YAG Laser Energy Levels

Objective: In this study we evaluated the ablation rate of superficial and deep dentin irradiated with different Er:YAG laser energy levels, and observed the micromorphological aspects of the lased substrates with a scanning electron microscope (SEM). Background Data: Little is known about the effect of Er: YAG laser irradiation on different dentin depths. Materials and Methods: Sixty molar crowns were bisected, providing 120 specimens, which were randomly assigned into two groups ( superficial or deep dentin), and later into five subgroups (160, 200, 260, 300, or 360 mJ). Initial masses of the specimens were obtained. After laser irradiation, the final masses were obtained and mass losses were calculated followed by the preparation of specimens for SEM examination. Mass-loss values were subjected to two-way ANOVA and Fisher's least significant difference multiple-comparison tests (p < 0.05). Results: There was no difference between superficial and deep dentin. A significant and gradual increase in the mass-loss values was reached when energies were raised, regardless of the dentin depth. The energy level of 360 mJ showed the highest values and was statistically significantly different from the other energy levels. The SEM images showed that deep dentin was more selectively ablated, especially intertubular dentin, promoting tubule protrusion. At 360 mJ the micromorphological features were similar for both dentin depths. Conclusion: The ablation rate did not depend on the depth of the dentin, and an energy level lower than 360 mJ is recommended to ablate both superficial and deep dentin effectively without causing tissue damage.

Influence of pulse repetition rate on temperature rise and working time during composite filling removal with the Er : YAG laser

Objective: The purpose of this study was to assess the efficacy of Er:YAG laser energy for composite resin removal and the influence of pulse repetition rate on the thermal alterations occurring during laser ablation. Materials and Methods: Composite resin filling was placed in cavities (1.0 mm deep) prepared in bovine teeth and the specimens were randomly assigned to five groups according to the technique used for composite filling removal. In group I (controls), the restorations were removed using a high-speed diamond bur. In the other groups, the composite fillings were removed using an Er: YAG laser with different pulse repetition rates: group 2-2 Hz; group 3-4 Hz; group 4-6 Hz; and group 5-10 Hz. The time required for complete removal of the restorative material and the temperature changes were recorded. Results: Temperature rise during composite resin removal with the Er: YAG laser occurred in the substrate underneath the restoration and was directly proportional to the increase in pulse repetition rate. None of the groups had a temperature increase during composite filling removal of more than 5.6 degrees C, which is considered the critical point above which irreversible thermal damage to the pulp may result. Regarding the time for composite filling removal, all the laser-ablated groups (except for group 5 [10 Hz]) required more time than the control group for complete elimination of the material from the cavity walls. Conclusion: Under the tested conditions, Er: YAG laser irradiation was efficient for composite resin ablation and did not cause a temperature increase above the limit considered safe for the pulp. Among the tested pulse repetition rates, 6 Hz produced minimal temperature change compared to the control group (high-speed bur), and allowed composite filling removal within a time period that is acceptable for clinical conditions.

The qWR star HD 45166 - II. Fundamental stellar parameters and evidence of a latitude-dependent wind

Context. The enigmatic object HD 45166 is a qWR star in a binary system with an orbital period of 1.596 day, and presents a rich emission-line spectrum in addition to absorption lines from the companion star (B7 V). As the system inclination is very small (i = 0.77 degrees +/- 0.09 degrees), HD 45166 is an ideal laboratory for wind-structure studies. Aims. The goal of the present paper is to determine the fundamental stellar and wind parameters of the qWR star. Methods. A radiative transfer model for the wind and photosphere of the qWR star was calculated using the non-LTE code CMFGEN. The wind asymmetry was also analyzed using a recently-developed version of CMFGEN to compute the emerging spectrum in two-dimensional geometry. The temporal-variance spectrum (TVS) was calculated to study the line-profile variations. Results. Abundances and stellar and wind parameters of the qWR star were obtained. The qWR star has an effective temperature of T(eff) = 50 000 +/- 2000 K, a luminosity of log(L/L(circle dot)) = 3.75 +/- 0.08, and a corresponding photospheric radius of R(phot) = 1.00 R(circle dot). The star is helium-rich (N(H)/N(He) = 2.0), while the CNO abundances are anomalous when compared either to solar values, to planetary nebulae, or to WR stars. The mass-loss rate is. M = 2.2 x 10(-7) M(circle dot) yr(-1), and the wind terminal velocity is v(infinity) = 425 km s(-1). The comparison between the observed line profiles and models computed under different latitude-dependent wind densities strongly suggests the presence of an oblate wind density enhancement, with a density contrast of at least 8: 1 from equator to pole. If a high velocity polar wind is present (similar to 1200 km s(-1)), the minimum density contrast is reduced to 4:1. Conclusions. The wind parameters determined are unusual when compared to O-type stars or to typical WR stars. While for WR stars v(infinity)/v(esc) > 1.5, in the case of HD 45166 it is much smaller (v(infinity)/v(esc) = 0.32). In addition, the efficiency of momentum transfer is eta = 0.74, which is at least 4 times smaller than in a typical WR. We find evidence for the presence of a wind compression zone, since the equatorial wind density is significantly higher than the polar wind. The TVS supports the presence of such a latitude-dependent wind and a variable absorption/scattering gas near the equator.

Chemical evolution models for spiral disks: the Milky Way, M31, and M33

Context. The distribution of chemical abundances and their variation with time are important tools for understanding the chemical evolution of galaxies. In particular, the study of chemical evolution models can improve our understanding of the basic assumptions made when modelling our Galaxy and other spirals. Aims. We test a standard chemical evolution model for spiral disks in the Local Universe and study the influence of a threshold gas density and different efficiencies in the star formation rate (SFR) law on radial gradients of abundance, gas, and SFR. The model is then applied to specific galaxies. Methods. We adopt a one-infall chemical evolution model where the Galactic disk forms inside-out by means of infall of gas, and we test different thresholds and efficiencies in the SFR. The model is scaled to the disk properties of three Local Group galaxies (the Milky Way, M31 and M33) by varying its dependence on the star formation efficiency and the timescale for the infall of gas onto the disk. Results. Using this simple model, we are able to reproduce most of the observed constraints available in the literature for the studied galaxies. The radial oxygen abundance gradients and their time evolution are studied in detail. The present day abundance gradients are more sensitive to the threshold than to other parameters, while their temporal evolutions are more dependent on the chosen SFR efficiency. A variable efficiency along the galaxy radius can reproduce the present day gas distribution in the disk of spirals with prominent arms. The steepness in the distribution of stellar surface density differs from massive to lower mass disks, owing to the different star formation histories. Conclusions. The most massive disks seem to have evolved faster (i.e., with more efficient star formation) than the less massive ones, thus suggesting a downsizing in star formation for spirals. The threshold and the efficiency of star formation play a very important role in the chemical evolution of spiral disks. For instance, an efficiency varying with radius can be used to regulate the star formation. The oxygen abundance gradient can steepen or flatten in time depending on the choice of this parameter.

BONA FIDE, STRONG-VARIABLE GALACTIC LUMINOUS BLUE VARIABLE STARS ARE FAST ROTATORS: DETECTION OF A HIGH ROTATIONAL VELOCITY IN HR CARINAE

We report optical observations of the luminous blue variable (LBV) HR Carinae which show that the star has reached a visual minimum phase in 2009. More importantly, we detected absorptions due to Si lambda lambda 4088-4116. To match their observed line profiles from 2009 May, a high rotational velocity of nu(rot) similar or equal to 150 +/- 20 km s(-1) is needed (assuming an inclination angle of 30 degrees), implying that HR Car rotates at similar or equal to 0.88 +/- 0.2 of its critical velocity for breakup (nu(crit)). Our results suggest that fast rotation is typical in all strong-variable, bona fide galactic LBVs, which present S-Dor-type variability. Strong-variable LBVs are located in a well-defined region of the HR diagram during visual minimum (the ""LBV minimum instability strip""). We suggest this region corresponds to where nu(crit) is reached. To the left of this strip, a forbidden zone with nu(rot)/nu(crit) > 1 is present, explaining why no LBVs are detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low nu(rot), we propose that LBVs can be separated into two groups: fast-rotating, strong-variable stars showing S-Dor cycles (such as AG Car and HR Car) and slow-rotating stars with much less variability (such as P Cygni and HD 168625). We speculate that supernova (SN) progenitors which had S-Dor cycles before exploding (such as in SN 2001ig, SN 2003bg, and SN 2005gj) could have been fast rotators. We suggest that the potential difficulty of fast-rotating Galactic LBVs to lose angular momentum is additional evidence that such stars could explode during the LBV phase.

Automated supervised classification of variable stars in the CoRoT programme Method and application to the first four exoplanet fields

Aims. In this work, we describe the pipeline for the fast supervised classification of light curves observed by the CoRoT exoplanet CCDs. We present the classification results obtained for the first four measured fields, which represent a one-year in-orbit operation. Methods. The basis of the adopted supervised classification methodology has been described in detail in a previous paper, as is its application to the OGLE database. Here, we present the modifications of the algorithms and of the training set to optimize the performance when applied to the CoRoT data. Results. Classification results are presented for the observed fields IRa01, SRc01, LRc01, and LRa01 of the CoRoT mission. Statistics on the number of variables and the number of objects per class are given and typical light curves of high-probability candidates are shown. We also report on new stellar variability types discovered in the CoRoT data. The full classification results are publicly available.

Populations III.1 and III.2 gamma-ray bursts: constraints on the event rate for future radio and X-ray surveys

Aims. We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to E(iso) greater than or similar to 10(55-57) erg, providing a unique probe of the high-redshift Universe. Methods. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. Results. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect similar to 10-10(4) radio afterglows above similar to 0.3 mJy on the sky with similar to 1 year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of N < 20 GRBs per year integrated over at z > 6 for Pop III.2 and N < 0.08 per year integrated over at z > 10 for Pop III.1 with EXIST, and N < 0.2 for Pop III.2 GRBs per year integrated over at z > 6 with Swift.

Exogenous Cx43 expression decrease cell proliferation rate in rat hepatocarcinoma cells independently of functional gap junction

Background: Gap junction intercellular communication (GJIC) is considered to play a role in the regulation of homeostasis because it regulates important processes, such as cell proliferation and cell differentiation. A reduced or lost GJIC capacity has been observed in solid tumors and studies have demonstrated that GJIC restoration in tumor cells contribute to reversion of the transformed phenotype. This observation supports the idea that restoration of the functional channel is essential in this process. However, in the last years, reports have proposed that just the increase in the expression of specific connexins can contribute to reversion of the malign phenotype in some tumor cells. In the present work, we studied the effects of exogenous Connexin 43 (Cx43) expression on the proliferative behavior and phenotype of rat hepatocarcinoma cells. Results: The exogenous Cx43 did not increase GJIC capacity of transfected cells, but it was critical to decrease the cell proliferation rate as well as reorganization of the actin filaments and cell flattening. We also observed more adhesion capacity to substrate after Cx43 transfection. Conclusion: Cx43 expression leads to a decrease of the growth of the rat hepatocellular carcinoma cells and it contributes to the reversion of the transformed phenotype. These effects were independent of the GJIC and were probably associated with the phosphorylation pattern changes and redistribution of the Cx43 protein.

Electron dephasing scattering rate in two-dimensional GaAs/InGaAs heterostructures with embedded InAs quantum dots

We report a comprehensive study of weak-localization and electron-electron interaction effects in a GaAs/InGaAs two-dimensional electron system with nearby InAs quantum dots, using measurements of the electrical conductivity with and without magnetic field. Although both the effects introduce temperature dependent corrections to the zero magnetic field conductivity at low temperatures, the magnetic field dependence of conductivity is dominated by the weak-localization correction. We observed that the electron dephasing scattering rate tau(-1)(phi), obtained from the magnetoconductivity data, is enhanced by introducing quantum dots in the structure, as expected, and obeys a linear dependence on the temperature and elastic mean free path, which is against the Fermi-liquid model. (c) 2008 American Institute of Physics. [DOI: 10.1063/1.2996034]

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.