Accelerating epistasis analysis in human genetics with consumer graphics hardware

BACKGROUND: Human geneticists are now capable of measuring more than one million DNA sequence variations from across the human genome. The new challenge is to develop computationally feasible methods capable of analyzing these data for associations with common human disease, particularly in the context of epistasis. Epistasis describes the situation where multiple genes interact in a complex non-linear manner to determine an individual's disease risk and is thought to be ubiquitous for common diseases. Multifactor Dimensionality Reduction (MDR) is an algorithm capable of detecting epistasis. An exhaustive analysis with MDR is often computationally expensive, particularly for high order interactions. This challenge has previously been met with parallel computation and expensive hardware. The option we examine here exploits commodity hardware designed for computer graphics. In modern computers Graphics Processing Units (GPUs) have more memory bandwidth and computational capability than Central Processing Units (CPUs) and are well suited to this problem. Advances in the video game industry have led to an economy of scale creating a situation where these powerful components are readily available at very low cost. Here we implement and evaluate the performance of the MDR algorithm on GPUs. Of primary interest are the time required for an epistasis analysis and the price to performance ratio of available solutions. FINDINGS: We found that using MDR on GPUs consistently increased performance per machine over both a feature rich Java software package and a C++ cluster implementation. The performance of a GPU workstation running a GPU implementation reduces computation time by a factor of 160 compared to an 8-core workstation running the Java implementation on CPUs. This GPU workstation performs similarly to 150 cores running an optimized C++ implementation on a Beowulf cluster. Furthermore this GPU system provides extremely cost effective performance while leaving the CPU available for other tasks. The GPU workstation containing three GPUs costs $2000 while obtaining similar performance on a Beowulf cluster requires 150 CPU cores which, including the added infrastructure and support cost of the cluster system, cost approximately$82,500. CONCLUSION: Graphics hardware based computing provides a cost effective means to perform genetic analysis of epistasis using MDR on large datasets without the infrastructure of a computing cluster

Absence of Street Lighting May Prevent Vehicle Crime, but Spatial and Temporal Displacement Remains a Concern

OBJECTIVES: This paper estimates the effect of changes in street lighting at night on levels of crime at street-level. Analyses investigate spatial and temporal displacement of crime into adjacent streets. METHODS: Offense data (burglaries, robberies, theft of and theft from vehicles, and violent crime) were obtained from Thames Valley Police, UK. Street lighting data (switching lights off at midnight, dimming, and white light) were obtained from local authorities. Monthly counts of crime at street-level were analyzed using a conditional fixed-effects Poisson regression model, adjusting for seasonal and temporal variation. Two sets of models analyzed: (1) changes in night-time crimes adjusting for changes in day-time crimes and (2) changes in crimes at all times of the day. RESULTS: Switching lights off at midnight was strongly associated with a reduction in night-time theft from vehicles relative to daytime (rate ratio RR 0.56; 0.41–0.78). Adjusted for changes in daytime, night-time theft from vehicles increased (RR 1.55; 1.14–2.11) in adjacent roads where street lighting remained unchanged. CONCLUSION: Theft from vehicle offenses reduced in streets where street lighting was switched off at midnight but may have been displaced to better-lit adjacent streets. Relative to daytime, night-time theft from vehicle offenses reduced in streets with dimming while theft from vehicles at all times of the day increased, thus suggesting temporal displacement. These findings suggest that the absence of street lighting may prevent theft from vehicles, but there is a danger of offenses being temporally or spatially displaced

Growth and dislocation studies of β-HMX

Background: The defect structure of organic materials is important as it plays a major role in their crystal growth properties. It also can play a subcritical role in “hot-spot” detonation processes of energetics and one such energetic is cyclotetramethylene-tetranitramine, in the commonly used beta form (β-HMX). Results: The as-grown crystals grown by evaporation from acetone show prismatic, tabular and columnar habits, all with {011}, {110}, (010) and (101) faces. Etching on (010) surfaces revealed three different types of etch pits, two of which could be identified with either pure screw or pure edge dislocations, the third is shown to be an artifact of the twinning process that this material undergoes. Examination of the {011} and {110} surfaces show only one type of etch pit on each surface; however their natural asymmetry precludes the easy identification of their Burgers vector or dislocation type. Etching of cleaved {011} surfaces demonstrates that the etch pits can be associated with line dislocations. All dislocations appear randomly on the crystal surfaces and do not form alignments characteristic of mechanical deformation by dislocation slip. Conclusions: Crystals of β-HMX grown from acetone show good morphological agreement with that predicted by modelling, with three distinct crystal habits observed depending upon the supersaturation of the growth solution. Prismatic habit was favoured at low supersaturation, while tabular and columnar crystals were predominant at higher super saturations. The twin plane in β-HMX was identified as a (101) reflection plane. The low plasticity of β-HMX is shown by the lack of etch pit alignments corresponding to mechanically induced dislocation arrays. On untwinned {010} faces, two types of dislocations exist, pure edge dislocations with b = [010] and pure screw dislocations with b = [010]. On twinned (010) faces, a third dislocation type exists and it is proposed that these pits are associated with pure screw dislocations with b = [010]

Comparing videofluoroscopy and endoscopy to assess swallowing in bottle-fed young infants in the neonatal intensive care unit

Contains fulltext : 235605.pdf (Publisher’s version ) (Closed access

Mangarara Formation: exhumed remnants of a middle Miocene, temperate carbonate, submarine channel-fan system on the eastern margin of Taranaki Basin, New Zealand

The middle Miocene Mangarara Formation is a thin (1–60 m), laterally discontinuous unit of moderately to highly calcareous (40–90%) facies of sandy to pure limestone, bioclastic sandstone, and conglomerate that crops out in a few valleys in North Taranaki across the transition from King Country Basin into offshore Taranaki Basin. The unit occurs within hemipelagic (slope) mudstone of Manganui Formation, is stratigraphically associated with redeposited sandstone of Moki Formation, and is overlain by redeposited volcaniclastic sandstone of Mohakatino Formation. The calcareous facies of the Mangarara Formation are interpreted to be mainly mass-emplaced deposits having channelised and sheet-like geometries, sedimentary structures supportive of redeposition, mixed environment fossil associations, and stratigraphic enclosure within bathyal mudrocks and flysch. The carbonate component of the deposits consists mainly of bivalves, larger benthic foraminifers (especially Amphistegina), coralline red algae including rhodoliths (Lithothamnion and Mesophyllum), and bryozoans, a warm-temperate, shallow marine skeletal association. While sediment derivation was partly from an eastern contemporary shelf, the bulk of the skeletal carbonate is inferred to have been sourced from shoal carbonate factories around and upon isolated basement highs (Patea-Tongaporutu High) to the south. The Mangarara sediments were redeposited within slope gullies and broad open submarine channels and lobes in the vicinity of the channel-lobe transition zone of a submarine fan system. Different phases of sediment transport and deposition (lateral-accretion and aggradation stages) are identified in the channel infilling. Dual fan systems likely co-existed, one dominating and predominantly siliciclastic in nature (Moki Formation), and the other infrequent and involving the temperate calcareous deposits of Mangarara Formation. The Mangarara Formation is an outcrop analogue for middle Miocene-age carbonate slope-fan deposits elsewhere in subsurface Taranaki Basin, New Zealand

Fluctuation induces evolutionary branching in a modeled microbial ecosystem

The impact of environmental fluctuation on species diversity is studied with a model of the evolutionary ecology of microorganisms. We show that environmental fluctuation induces evolutionary branching and assures the consequential coexistence of multiple species. Pairwise invasibility analysis is applied to illustrate the speciation process. We also discuss how fluctuation affects species diversity.Comment: 4 pages, 4 figures. Submitted to Physical Review Letter

Uncovering the most kinetically influential reaction pathway driving the generation of HCN from oxyma/DIC adduct: a theoretical study

The combination of ethyl (hydroxyimino)cyanoacetate (Oxyma) and diisopropylcarbodiimide (DIC) has demonstrated superior performance in amino acid activation for peptide synthesis. However, it was recently reported that Oxyma and DIC could react to generate undesired hydrogen cyanide (HCN) at 20 °C, raising safety concerns for the practical use of this activation strategy. To help minimize the risks, there is a need for a comprehensive investigation of the mechanism and kinetics of the generation of HCN. Here we show the results of the first systematic computational study of the underpinning mechanism, including comparisons with experimental data. Two pathways for the decomposition of the Oxyma/DIC adduct are derived to account for the generation of HCN and its accompanying cyclic product. These two mechanisms differ in the electrophilic carbon atom attacked by the nucleophilic sp2-nitrogen in the cyclization step and in the cyclic product generated. On the basis of computed “observed” activation energies, ΔGobs⧧, the mechanism that proceeds via the attack of the sp2-nitrogen at the oxime carbon is identified as the most kinetically favorable one, a conclusion that is supported by closer agreement between predicted and experimental 13C NMR data. These results can provide a theoretical basis to develop a design strategy for suppressing HCN generation when using Oxyma/DIC for amino acid activation

Environmental controls on the elemental composition of a Southern Hemisphere strain of the coccolithophore Emiliania huxleyi

A series of semi-continuous incubation experiments were conducted with the coccolithophore Emiliania huxleyi strain NIWA1108 (Southern Ocean isolate) to examine the effects of five environmental drivers (nitrate and phosphate concentrations, irradiance, temperature, and partial pressure of CO2 (pCO2)) on both the physiological rates and elemental composition of the coccolithophore. Here, we report the alteration of the elemental composition of E. huxleyi in response to the changes in these environmental drivers. A series of dose–response curves for the cellular elemental composition of E. huxleyi were fitted for each of the five drivers across an environmentally representative gradient. The importance of each driver in regulating the elemental composition of E. huxleyi was ranked using a semi-quantitative approach. The percentage variations in elemental composition arising from the change in each driver between present-day and model-projected conditions for the year 2100 were calculated. Temperature was the most important driver controlling both cellular particulate organic and inorganic carbon content, whereas nutrient concentrations were the most important regulator of cellular particulate nitrogen and phosphorus of E. huxleyi. In contrast, elevated pCO2 had the greatest influence on cellular particulate inorganic carbon to organic carbon ratio, resulting in a decrease in the ratio. Our results indicate that the different environmental drivers play specific roles in regulating the elemental composition of E. huxleyi with wide-reaching implications for coccolithophore-related marine biogeochemical cycles, as a consequence of the regulation of E. huxleyi physiological processes

Genetic variation in the HLA region is associated with susceptibility to herpes zoster.

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22,981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10(-8)). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine