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.

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus ( Sigmodontinae), Dobrava virus ( Murinae), Puumala virus ( Arvicolinae), and Tula virus ( Arvicolinae). Our results reveal that hantaviruses exhibit shortterm substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

Lossless filter for multiple repeats with bounded edit distance

Background: Identifying local similarity between two or more sequences, or identifying repeats occurring at least twice in a sequence, is an essential part in the analysis of biological sequences and of their phylogenetic relationship. Finding such fragments while allowing for a certain number of insertions, deletions, and substitutions, is however known to be a computationally expensive task, and consequently exact methods can usually not be applied in practice. Results: The filter TUIUIU that we introduce in this paper provides a possible solution to this problem. It can be used as a preprocessing step to any multiple alignment or repeats inference method, eliminating a possibly large fraction of the input that is guaranteed not to contain any approximate repeat. It consists in the verification of several strong necessary conditions that can be checked in a fast way. We implemented three versions of the filter. The first is simply a straightforward extension to the case of multiple sequences of an application of conditions already existing in the literature. The second uses a stronger condition which, as our results show, enable to filter sensibly more with negligible (if any) additional time. The third version uses an additional condition and pushes the sensibility of the filter even further with a non negligible additional time in many circumstances; our experiments show that it is particularly useful with large error rates. The latter version was applied as a preprocessing of a multiple alignment tool, obtaining an overall time (filter plus alignment) on average 63 and at best 530 times smaller than before (direct alignment), with in most cases a better quality alignment. Conclusion: To the best of our knowledge, TUIUIU is the first filter designed for multiple repeats and for dealing with error rates greater than 10% of the repeats length.

Anemia screening in potential female blood donors: comparison of two different quantitative methods

Anemia screening before blood donation requires an accurate, quick, practical, and easy method with minimal discomfort for the donors. The aim of this study was to compare the accuracy of two quantitative methods of anemia screening: the HemoCue 201(+) (Aktiebolaget Leo Diagnostics) hemoglobin (Hb) and microhematocrit (micro-Hct) tests. Two blood samples of a single fingerstick were obtained from 969 unselected potential female donors to determine the Hb by HemoCue 201(+) and micro-Hct using HemataSTAT II (Separation Technology, Inc.), in alternating order. From each participant, a venous blood sample was drawn and run in an automatic hematology analyzer (ABX Pentra 60, ABX Diagnostics). Considering results of ABX Pentra 60 as true values, the sensitivity and specificity of HemoCue 201(+) and micro-Hct as screening methods were compared, using a venous Hb level of 12.0 g per dL as cutoff for anemia. The sensitivities of the HemoCue 201(+) and HemataSTAT II in detecting anemia were 56 percent (95% confidence interval [CI], 46.1%-65.5%) and 39.5 percent (95% CI, 30.2%-49.3%), respectively (p < 0.001). Analyzing only candidates with a venous Hb level lower than 11.0 g per dL, the deferral rate was 100 percent by HemoCue 201(+) and 77 percent by HemataSTAT II. The specificities of the methods were 93.5 and 93.2 percent, respectively. The HemoCue 201(+) showed greater discriminating power for detecting anemia in prospective blood donors than the micro-Hct method. Both presented equivalent deferral error rates of nonanemic potential donors. Compared to the micro-Hct, HemoCue 201(+) reduces the risk of anemic female donors giving blood, specially for those with lower Hb levels, without increasing the deferral of nonanemic potential donors.

Cutpoint selection for discretizing a continuous covariate for generalized estimating equations

We consider consider the problem of dichotomizing a continuous covariate when performing a regression analysis based on a generalized estimation approach. The problem involves estimation of the cutpoint for the covariate and testing the hypothesis that the binary covariate constructed from the continuous covariate has a significant impact on the outcome. Due to the multiple testing used to find the optimal cutpoint, we need to make an adjustment to the usual significance test to preserve the type-I error rates. We illustrate the techniques on one data set of patients given unrelated hematopoietic stem cell transplantation. Here the question is whether the CD34 cell dose given to patient affects the outcome of the transplant and what is the smallest cell dose which is needed for good outcomes. (C) 2010 Elsevier BM. All rights reserved.

Melatonin in maturation media fails to improve oocyte maturation, embryo development rates and DNA damage of bovine embryos

Melatonin (MEL) acts as a powerful scavenger of free radicals and direct gonadal responses to melatonin have been reported in the literature. Few studies, however, have evaluated the effect of MEL during in vitro maturation (IVM) on bovine embryos. This study tested the addition of MEL to maturation medium (MM) with no gonadotropins on nuclear maturation and embryo development rates and the incidence of DNA damage in resulting embryos. Cumulus-oocyte complexes were aspirated from abattoir ovaries and cultured in MM (TCM-199 medium supplemented with 10% fetal calf serum - FCS) at 39ºC and 5% CO2 in air. After 24 hours of culture in MM with 0.5 µg mL-1 FSH and 5.0 µg mL-1 LH; 10-9 M MEL) or 10-9 M MEL, 0.5 µg mL-1 FSH and 5.0 µg mL-1 LH, the oocytes were stained with Hoechst 33342 to evaluate nuclear maturation rate. After in vitro fertilization and embryo culture, development rates were evaluated and the blastocysts were assessed for DNA damage by Comet assay. There was no effect of melatonin added to the MM, alone or in combination with gonadotropins, on nuclear maturation, cleavage and blastocyst rates. These rates ranged between 88% to 90%, 85% to 88% and 42% to 46%, respectively. The extent of DNA damage in embryos was also not affected by MEL supplementation during IVM. The addition of 10-9 M MEL to the MM failed to improve nuclear maturation and embryo development rates and the incidence of DNA damage in resulting embryos, but was able to properly substitute for gonadotropins during IVM.

Deleterious effect of TRIS buffer on growth rates and pigment content of Gracilaria birdiae Plastino & E.C. Oliveira (Gracilariales, Rhodophyta)

This work evaluated the effects of Tris (hydroxymethyl)-aminomethane (TRIS) buffer and its interaction with nutrient concentration on the development of Gracilaria birdie, a common species on the Brazilian coast that has been exploited for agar production. Responses to different conditions were assessed through growth rates and pigment content (chlorophyll a, phycoerythrin, phycocyanin and allophycocyanin). Provasoli's nutrient solution with and without TRIS addition was tested at concentrations of 12.5, 25 and 50%. The pH was also monitored. G. birdiae grew better in the absence of TRIS and at low nutrient concentrations, 12.5 and 25% (growth rates of 10.8-11.3%.day-1). Higher contents of phycoerythrin and chlorophyll a were observed without TRIS at 12.5 and 25% (Phycoerythrin, 649.6-698.0 μg g-1 fresh biomass; Chlorophyll a, 156.0-168.6 μg g-1 fresh biomass). These findings highlight the deleterious effect of TRIS on growth and phycoerythrin and chlorophyll a content. They also demonstrate the importance of appropriate nutrient concentration for laboratory cultures, depending on the intrinsic characteristics of each species.

Effects of temperature and salinity on the survival rates of coxicerberus ramosae (Albuquerque, 1978), an interstitial isopod of a Sandy Beach on the coast of Brazil

The tolerance to the combined effects of temperature and salinity was investigated in the interstitial isopod Coxicerberus ramosae (Albuquerque, 1978), a species of intertidal zone of sandy beaches in Rio de Janeiro, Brazil. The animals were collected on Praia Vermelha Beach. The experiments lasted 24 h and nine salinities and seven temperatures were used for a total of 63 combinations. Thirty animals were tested in each combination. The species showed high survival in most of the combinations. The temperature of 35 ºC was lethal and at 5 ºC, the animals tolerated only a narrow range of salinities. The statistical analyses showed that the effects of temperature and salinity were significant on the survival, which confirmed the euryhalinity and eurythermy of this species.

Rates of antimicrobial resistance in latin america (2004-2007) and in vitro activity of the glycylcycline tigecycline and of other antibiotics

As a part of the Tigecycline Evaluation and Surveillance Trial (T.E.S.T.), Gram-positive and Gram-negative bacterial isolates were collected from 33 centers in Latin America (centers in Argentina, Brazil, Chile, Colombia, Guatemala, Honduras, Jamaica, Mexico, Panama, Puerto Rico, and Venezuela) from January 2004 to September 2007. Argentina and Mexico were the greatest contributors of isolates to this study. Susceptibilities were determined according to Clinical Laboratory Standards Institute guidelines. Resistance levels were high for most key organisms across Latin America: 48.3% of Staphylococcus aureus isolates were methicillin-resistant while 21.4% of Acinetobacter spp. isolates were imipenem-resistant. Extended-spectrum β-lactamase were reported in 36.7% of Klebsiella pneumoniae and 20.8% of E. coli isolates. Tigecycline was the most active agent against Gram-positive isolates. Tigecycline was also highly active against all Gram-negative organisms, with the exception of Pseuodomonas aeruginosa, against which piperacillin-tazobactam was the most active agent tested (79.3% of isolates susceptible). The in vitro activity of tigecycline against both Gram-positive and Gram-negative isolates indicates that it may be an useful tool for the treatment of nosocomial infections, even those caused by organisms that are resistant to other antibacterial agents.

Simplified Mathematical Model for an Anaerobic Sequencing Batch Biofilm Reactor Treating Lipid-Rich Wastewater Subject to Rising Organic Loading Rates

This study proposes a simplified mathematical model to describe the processes occurring in an anaerobic sequencing batch biofilm reactor (ASBBR) treating lipid-rich wastewater. The reactor, subjected to rising organic loading rates, contained biomass immobilized cubic polyurethane foam matrices, and was operated at 32 degrees C +/- 2 degrees C, using 24-h batch cycles. In the adaptation period, the reactor was fed with synthetic substrate for 46 days and was operated without agitation. Whereas agitation was raised to 500 rpm, the organic loading rate (OLR) rose from 0.3 g chemical oxygen demand (COD) . L(-1) . day(-1) to 1.2 g COD . L(-1) . day(-1). The ASBBR was fed fat-rich wastewater (dairy wastewater), in an operation period lasting for 116 days, during which four operational conditions (OCs) were tested: 1.1 +/- 0.2 g COD . L(-1) . day(-1) (OC1), 4.5 +/- 0.4 g COD . L(-1) . day(-1) (OC2), 8.0 +/- 0.8 g COD . L(-1) . day(-1) (OC3), and 12.1 +/- 2.4 g COD . L(-1) . day(-1) (OC4). The bicarbonate alkalinity (BA)/COD supplementation ratio was 1:1 at OC1, 1:2 at OC2, and 1:3 at OC3 and OC4. Total COD removal efficiencies were higher than 90%, with a constant production of bicarbonate alkalinity, in all OCs tested. After the process reached stability, temporal profiles of substrate consumption were obtained. Based on these experimental data a simplified first-order model was fit, making possible the inference of kinetic parameters. A simplified mathematical model correlating soluble COD with volatile fatty acids (VFA) was also proposed, and through it the consumption rates of intermediate products as propionic and acetic acid were inferred. Results showed that the microbial consortium worked properly and high efficiencies were obtained, even with high initial substrate concentrations, which led to the accumulation of intermediate metabolites and caused low specific consumption rates.

Continental passive microwave-based rainfall estimation algorithm: Application to the Amazon Basin

This paper presents a new statistical algorithm to estimate rainfall over the Amazon Basin region using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm relies on empirical relationships derived for different raining-type systems between coincident measurements of surface rainfall rate and 85-GHz polarization-corrected brightness temperature as observed by the precipitation radar (PR) and TMI on board the TRMM satellite. The scheme includes rain/no-rain area delineation (screening) and system-type classification routines for rain retrieval. The algorithm is validated against independent measurements of the TRMM-PR and S-band dual-polarization Doppler radar (S-Pol) surface rainfall data for two different periods. Moreover, the performance of this rainfall estimation technique is evaluated against well-known methods, namely, the TRMM-2A12 [ the Goddard profiling algorithm (GPROF)], the Goddard scattering algorithm (GSCAT), and the National Environmental Satellite, Data, and Information Service (NESDIS) algorithms. The proposed algorithm shows a normalized bias of approximately 23% for both PR and S-Pol ground truth datasets and a mean error of 0.244 mm h(-1) ( PR) and -0.157 mm h(-1)(S-Pol). For rain volume estimates using PR as reference, a correlation coefficient of 0.939 and a normalized bias of 0.039 were found. With respect to rainfall distributions and rain area comparisons, the results showed that the formulation proposed is efficient and compatible with the physics and dynamics of the observed systems over the area of interest. The performance of the other algorithms showed that GSCAT presented low normalized bias for rain areas and rain volume [0.346 ( PR) and 0.361 (S-Pol)], and GPROF showed rainfall distribution similar to that of the PR and S-Pol but with a bimodal distribution. Last, the five algorithms were evaluated during the TRMM-Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) 1999 field campaign to verify the precipitation characteristics observed during the easterly and westerly Amazon wind flow regimes. The proposed algorithm presented a cumulative rainfall distribution similar to the observations during the easterly regime, but it underestimated for the westerly period for rainfall rates above 5 mm h(-1). NESDIS(1) overestimated for both wind regimes but presented the best westerly representation. NESDIS(2), GSCAT, and GPROF underestimated in both regimes, but GPROF was closer to the observations during the easterly flow.

Unsuccessful attacks dominate a drone-preying wasp's hunting performance near stingless bee nests

Bee males (drones) of stingless bees tend to congregate near entrances of conspecific nests, where they wait for virgin queens that initiate their nuptial flight. We observed that the Neotropical solitary wasp Trachypus boharti (Hymenoptera, Cabronidae) specifically preys on males of the stingless bee Scaptotrigona postica (Hymenoptera, Apidae); these wasps captured up to 50 males per day near the entrance of a single hive. Over 90% of the wasp attacks were unsuccessful; such erroneous attacks often involved conspecific wasps and worker bees. After the capture of non-male prey, wasps almost immediately released these individuals unharmed and continued hunting. A simple behavioral experiment showed that at short distances wasps were not specifically attracted to S. postica males nor were they repelled by workers of the same species. Likely, short-range prey detection near the bees' nest is achieved mainly by vision whereas close-range prey recognition is based principally on chemical and/or mechanical cues. We argue that the dependence on the wasp's visual perception during attack and the crowded and dynamic hunting conditions caused wasps to make many preying attempts that failed. Two wasp-density-related factors, wasp-prey distance and wasp-wasp encounters, may account for the fact that the highest male capture and unsuccessful wasp bee encounter rates occurred at intermediate wasp numbers.

Late Entry into HIV Care: Estimated Impact on AIDS Mortality Rates in Brazil, 2003-2006

Background: Worldwide, a high proportion of HIV-infected individuals enter into HIV care late. Here, our objective was to estimate the impact that late entry into HIV care has had on AIDS mortality rates in Brazil. Methodology/Principal Findings: We analyzed data from information systems regarding HIV-infected adults who sought treatment at public health care facilities in Brazil from 2003 to 2006. We initially estimated the prevalence of late entry into HIV care, as well as the probability of death in the first 12 months, the percentage of the risk of death attributable to late entry, and the number of avoidable deaths. We subsequently adjusted the annual AIDS mortality rate by excluding such deaths. Of the 115,369 patients evaluated, 50,358 (43.6%) had entered HIV care late, and 18,002 died in the first 12 months, representing a 16.5% probability of death in the first 12 months (95% CI: 16.3-16.7). By comparing patients who entered HIV care late with those who gained timely access, we found that the risk ratio for death was 49.5 (95% CI: 45.1-54.2). The percentage of the risk of death attributable to late entry was 95.5%, translating to 17,189 potentially avoidable deaths. Averting those deaths would have lowered the 2003-2006 AIDS mortality rate by 39.5%. Including asymptomatic patients with CD4(+) T cell counts >200 and <= 350 cells/mm(3) in the group who entered HIV care late increased this proportion by 1.8%. Conclusions/Significance: In Brazil, antiretroviral drugs reduced AIDS mortality by 43%. Timely entry would reduce that rate by a similar proportion, as well as resulting in a 45.2% increase in the effectiveness of the program for HIV care. The World Health Organization recommendation that asymptomatic patients with CD4(+) T cell counts <= 350 cells/mm(3) be treated would not have a significant impact on this scenario.

The impact of measurement errors in the identification of regulatory networks

Background: There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. Results: This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. Conclusions: Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.

Functional specialization of Eucalyptus fine roots: contrasting potential uptake rates for nitrogen, potassium and calcium tracers at varying soil depths

1. Little is known about the role of deep roots in the nutrition of forest trees and their ability to provide a safety-net service taking up nutrients leached from the topsoil. 2. To address this issue, we studied the potential uptake of N, K and Ca by Eucalyptus grandis trees (6 years of age - 25 m mean height), in Brazil, as a function of soil depth, texture and water content. We injected NO(3)(-)- (15)N, Rb(+) (analogue of K(+)) and Sr(2+) (analogue of Ca(2+)) tracers simultaneously in a solution through plastic tubes at 10, 50, 150 and 300 cm in depth in a sandy and a clayey Ferralsol soil. A complete randomized design was set up with three replicates of paired trees per injection depth and soil type. Recently expanded leaves were sampled at various times after tracer injection in the summer, and the experiment was repeated in the winter. Soil water contents were continuously monitored at the different depths in the two soils. 3. Determination of foliar Rb and Sr concentrations and (15)N atom % made it possible to estimate the relative uptake potential (RUP) of tracer injections from the four soil depths and the specific RUP (SRUP), defined as RUP, per unit of fine root length density in the corresponding soil layer. 4. The highest tracer uptake rates were found in the topsoil, but contrasting RUP distributions were observed for the three tracers. Whilst the RUP was higher for NO(3)(-)- (15)N than for Rb(+) and Sr(2+) in the upper 50 cm of soil, the highest SRUP values for Sr(2+) and Rb(+) were found at a depth of 300 cm in the sandy soil, as well as in the clayey soil when gravitational solutions reached that depth. 5. Our results suggest that the fine roots of E. grandis trees exhibit contrasting potential uptake rates with depth depending on the nutrient. This functional specialization of roots might contribute to the high growth rates of E. grandis trees, efficiently providing the large amounts of nutrients required throughout the development of these fast-growing plantations.