Accelerating Bayesian hierarchical clustering of time series data with a randomised algorithm

We live in an era of abundant data. This has necessitated the development of new and innovative statistical algorithms to get the most from experimental data. For example, faster algorithms make practical the analysis of larger genomic data sets, allowing us to extend the utility of cutting-edge statistical methods. We present a randomised algorithm that accelerates the clustering of time series data using the Bayesian Hierarchical Clustering (BHC) statistical method. BHC is a general method for clustering any discretely sampled time series data. In this paper we focus on a particular application to microarray gene expression data. We define and analyse the randomised algorithm, before presenting results on both synthetic and real biological data sets. We show that the randomised algorithm leads to substantial gains in speed with minimal loss in clustering quality. The randomised time series BHC algorithm is available as part of the R package BHC, which is available for download from Bioconductor (version 2.10 and above) via http://bioconductor.org/packages/2.10/bioc/html/BHC.html. We have also made available a set of R scripts which can be used to reproduce the analyses carried out in this paper. These are available from the following URL. https://sites.google.com/site/randomisedbhc/

Atmospheric emissions from the deepwater Horizon spill constrain air-water partitioning, hydrocarbon fate, and leak rate

The fate of deepwater releases of gas and oil mixtures is initially determined by solubility and volatility of individual hydrocarbon species; these attributes determine partitioning between air and water. Quantifying this partitioning is necessary to constrain simulations of gas and oil transport, to predict marine bioavailability of different fractions of the gas-oil mixture, and to develop a comprehensive picture of the fate of leaked hydrocarbons in the marine environment. Analysis of airborne atmospheric data shows massive amounts (∼258,000 kg/day) of hydrocarbons evaporating promptly from the Deepwater Horizon spill; these data collected during two research flights constrain air-water partitioning, thus bioavailability and fate, of the leaked fluid. This analysis quantifies the fraction of surfacing hydrocarbons that dissolves in the water column (∼33% by mass), the fraction that does not dissolve, and the fraction that evaporates promptly after surfacing (∼14% by mass). We do not quantify the leaked fraction lacking a surface expression; therefore, calculation of atmospheric mass fluxes provides a lower limit to the total hydrocarbon leak rate of 32,600 to 47,700 barrels of fluid per day, depending on reservoir fluid composition information. This study demonstrates a new approach for rapid-response airborne assessment of future oil spills. Copyright 2011 by the American Geophysical Union

Adaptive Evolution of the Myo6 Gene in Old World Fruit Bats (Family: Pteropodidae)

PMCID: PMC3631194This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Elevated Ratio of Urinary Metabolites of Thromboxane and Prostacyclin Is Associated with Adverse Cardiovascular Events in ADAPT

Results from prevention trials, including the Alzheimer's Disease Anti-inflammatory Prevention Trial (ADAPT), have fueled discussion about the cardiovascular (CV) risks associated with non-steroidal anti-inflammatory drugs (NSAIDs). We tested the hypotheses that (i) adverse CV events reported among ADAPT participants (aged 70 years and older) are associated with increased ratio of urine 11-dehydrothromboxane B2 (Tx-M) to 2′3-donor–6-keto-PGF1 (PGI-M) attributable to NSAID treatments; (ii) coincident use of aspirin (ASA) would attenuate NSAID-induced changes in Tx-M/PGI-M ratio; and (iii) use of NSAIDs and/or ASA would not alter urine or plasma concentrations of F2-isoprostanes (IsoPs), in vivo biomarkers of free radical damage. We quantified urine Tx-M and PGI-M, and urine and plasma F2-IsoPs from 315 ADAPT participants using stable isotope dilution assays with gas chromatography/mass spectrometry, and analyzed these data by randomized drug assignment and self-report compliance as well as ASA use. Adverse CV events were significantly associated with higher urine Tx-M/PGI-M ratio, which seemed to derive mainly from lowered PGI-M. Participants taking ASA alone had reduced urine Tx-M/PGI-M compared to no ASA or NSAID; however, participants taking NSAIDs plus ASA did not have reduced urine Tx-M/PGI-M ratio compared to NSAIDs alone. Neither NSAID nor ASA use altered plasma or urine F2-IsoPs. These data suggest a possible mechanism for the increased risk of CV events reported in ADAPT participants assigned to NSAIDs, and suggest that the changes in the Tx-M/PGI-M ratio was not substantively mitigated by coincident use of ASA in individuals 70 years or older

An intranasal selective antisense oligonucleotide impairs lung cyclooxygenase-2 production and improves inflammation, but worsens airway function, in house dust mite sensitive mice

<p>Abstract</p> <p>Background</p> <p>Despite its reported pro-inflammatory activity, cyclooxygenase (COX)-2 has been proposed to play a protective role in asthma. Accordingly, COX-2 might be down-regulated in the airway cells of asthmatics. This, together with results of experiments to assess the impact of COX-2 blockade in ovalbumin (OVA)-sensitized mice in vivo, led us to propose a novel experimental approach using house dust mite (HDM)-sensitized mice in which we mimicked altered regulation of COX-2.</p> <p>Methods</p> <p>Allergic inflammation was induced in BALBc mice by intranasal exposure to HDM for 10 consecutive days. This model reproduces spontaneous exposure to aeroallergens by asthmatic patients. In order to impair, but not fully block, COX-2 production in the airways, some of the animals received an intranasal antisense oligonucleotide. Lung COX-2 expression and activity were measured along with bronchovascular inflammation, airway reactivity, and prostaglandin production.</p> <p>Results</p> <p>We observed impaired COX-2 mRNA and protein expression in the lung tissue of selective oligonucleotide-treated sensitized mice. This was accompanied by diminished production of mPGE synthase and PGE₂in the airways. In sensitized mice, the oligonucleotide induced increased airway hyperreactivity (AHR) to methacholine, but a substantially reduced bronchovascular inflammation. Finally, mRNA levels of hPGD synthase remained unchanged.</p> <p>Conclusion</p> <p>Intranasal antisense therapy against COX-2 in vivo mimicked the reported impairment of COX-2 regulation in the airway cells of asthmatic patients. This strategy revealed an unexpected novel dual effect: inflammation was improved but AHR worsened. This approach will provide insights into the differential regulation of inflammation and lung function in asthma, and will help identify pharmacological targets within the COX-2/PG system.</p

Candida albicans Scavenges Host Zinc via Pra1 during Endothelial Invasion

The ability of pathogenic microorganisms to assimilate essential nutrients from their hosts is critical for pathogenesis. Here we report endothelial zinc sequestration by the major human fungal pathogen, Candida albicans. We hypothesised that, analogous to siderophore-mediated iron acquisition, C. albicans utilises an extracellular zinc scavenger for acquiring this essential metal. We postulated that such a “zincophore” system would consist of a secreted factor with zinc-binding properties, which can specifically reassociate with the fungal cell surface. In silico analysis of the C. albicans secretome for proteins with zinc binding motifs identified the pH-regulated antigen 1 (Pra1). Three-dimensional modelling of Pra1 indicated the presence of at least two zinc coordination sites. Indeed, recombinantly expressed Pra1 exhibited zinc binding properties in vitro. Deletion of PRA1 in C. albicans prevented fungal sequestration and utilisation of host zinc, and specifically blocked host cell damage in the absence of exogenous zinc. Phylogenetic analysis revealed that PRA1 arose in an ancient fungal lineage and developed synteny with ZRT1 (encoding a zinc transporter) before divergence of the Ascomycota and Basidiomycota. Structural modelling indicated physical interaction between Pra1 and Zrt1 and we confirmed this experimentally by demonstrating that Zrt1 was essential for binding of soluble Pra1 to the cell surface of C. albicans. Therefore, we have identified a novel metal acquisition system consisting of a secreted zinc scavenger (“zincophore”), which reassociates with the fungal cell. Furthermore, functional similarities with phylogenetically unrelated prokaryotic systems indicate that syntenic zinc acquisition loci have been independently selected during evolution

Axons Amplify Somatic Incomplete Spikes into Uniform Amplitudes in Mouse Cortical Pyramidal Neurons

BACKGROUND: Action potentials are the essential unit of neuronal encoding. Somatic sequential spikes in the central nervous system appear various in amplitudes. To be effective neuronal codes, these spikes should be propagated to axonal terminals where they activate the synapses and drive postsynaptic neurons. It remains unclear whether these effective neuronal codes are based on spike timing orders and/or amplitudes. METHODOLOGY/PRINCIPAL FINDINGS: We investigated this fundamental issue by simultaneously recording the axon versus soma of identical neurons and presynaptic vs. postsynaptic neurons in the cortical slices. The axons enable somatic spikes in low amplitude be enlarged, which activate synaptic transmission in consistent patterns. This facilitation in the propagation of sequential spikes through the axons is mechanistically founded by the short refractory periods, large currents and high opening probability of axonal voltage-gated sodium channels. CONCLUSION/SIGNIFICANCE: An amplification of somatic incomplete spikes into axonal complete ones makes sequential spikes to activate consistent synaptic transmission. Therefore, neuronal encoding is likely based on spike timing order, instead of graded analogues

Associations of ATR and CHEK1 Single Nucleotide Polymorphisms with Breast Cancer

DNA damage and replication checkpoints mediated by the ATR-CHEK1 pathway are key to the maintenance of genome stability, and both ATR and CHEK1 have been proposed as potential breast cancer susceptibility genes. Many novel variants recently identified by the large resequencing projects have not yet been thoroughly tested in genome-wide association studies for breast cancer susceptibility. We therefore used a tagging SNP (tagSNP) approach based on recent SNP data available from the 1000 genomes projects, to investigate the roles of ATR and CHEK1 in breast cancer risk and survival. ATR and CHEK1 tagSNPs were genotyped in the Sheffield Breast Cancer Study (SBCS; 1011 cases and 1024 controls) using Illumina GoldenGate assays. Untyped SNPs were imputed using IMPUTE2, and associations between genotype and breast cancer risk and survival were evaluated using logistic and Cox proportional hazard regression models respectively on a per allele basis. Significant associations were further examined in a meta-analysis of published data or confirmed in the Utah Breast Cancer Study (UBCS). The most significant associations for breast cancer risk in SBCS came from rs6805118 in ATR (p=7.6x10-5) and rs2155388 in CHEK1 (p=3.1x10-6), but neither remained significant after meta-analysis with other studies. However, meta-analysis of published data revealed a weak association between the ATR SNP rs1802904 (minor allele frequency is 12%) and breast cancer risk, with a summary odds ratio (confidence interval) of 0.90 (0.83-0.98) [p=0.0185] for the minor allele. Further replication of this SNP in larger studies is warranted since it is located in the target region of 2 microRNAs. No evidence of any survival effects of ATR or CHEK1 SNPs were identified. We conclude that common alleles of ATR and CHEK1 are not implicated in breast cancer risk or survival, but we cannot exclude effects of rare alleles and of common alleles with very small effect sizes

Novel M tuberculosis Antigen-Specific T-Cells Are Early Markers of Infection and Disease Progression

Mycobacterium tuberculosis Region-of-Difference-1 gene products present opportunities for specific diagnosis of M. tuberculosis infection, yet immune responses to only two gene-products, Early Secretory Antigenic Target-6 (ESAT-6) and Culture Filtrate Protein-10 (CFP-10), have been comprehensively investigated.T-cell responses to Rv3873, Rv3878 and Rv3879c were quantified by IFN-γ-enzyme-linked-immunospot (ELISpot) in 846 children with recent household tuberculosis exposure and correlated with kinetics of tuberculin skin test (TST) and ESAT-6/CFP-10-ELISpot conversion over six months and clinical outcome over two years.Responses to Rv3873, Rv3878, and Rv3879c were present in 20-25% of contacts at enrolment. Rv3873 and Rv3879c responses were associated with and preceded TST conversion (P=0.02 and P=0.04 respectively), identifying these antigens as early targets of cell-mediated immunity following M. tuberculosis exposure. Responses to Rv3873 were additionally associated with subsequent ESAT-6/CFP-10-ELISpot conversion (P=0.04). Responses to Rv3873 and Rv3878 predicted progression to active disease (adjusted incidence rate ratio [95% CI] 3.06 [1.05,8.95; P=0.04], and 3.32 [1.14,9.71; P=0.03], respectively). Presence of a BCG-vaccination scar was associated with a 67% (P=0.03) relative risk reduction for progression to active tuberculosis.These RD1-derived antigens are early targets of cellular immunity following tuberculosis exposure and T-cells specific for these antigens predict progression to active tuberculosis suggesting diagnostic and prognostic utility