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

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]

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

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

Mechanical properties of β-HMX

Background: For a full understanding of the mechanical properties of a material, it is essential to understand the defect structures and associated properties and microhardness indentation is a technique that can aid this understanding. Results: The Vickers hardness on (010), {011} and {110} faces lay in the range of 304-363 MPa. The Knoop Hardnesses on the same faces lay in the range 314-482 MPa. From etching of three indented surfaces, the preferred slip planes have been identified as (001) and (101). For a dislocation glide, the most likely configuration for dislocation movement on the (001) planes is (001) [100] (|b| = 0.65 nm) and for the (101) plane as (101) 101~(|b| = 1.084 nm) although (101) [010] (|b| = 1.105 nm) is possible. Tensile testing showed that at a stress value of 2.3 MPa primary twinning occurred and grew with increasing stress. When the stress was relaxed, the twins decreased in size, but did not disappear. The twinning shear strain was calculated to be 0.353 for the (101) twin plane. Conclusions: HMX is considered to be brittle, compared to other secondary explosives. Comparing HMX with a range of organic solids, the values for hardness numbers are similar to those of other brittle systems. Under the conditions developed beneath a pyramidal indenter, dislocation slip plays a major part in accommodating the local deformation stresses. © 2015 Gallagher et al.; licensee Springer

Histone deacetylases as new therapy targets for platinum-resistant epithelial ovarian cancer

Introduction: In developed countries, ovarian cancer is the fourth most common cancer in women. Due to the nonspecific symptomatology associated with the disease many patients with ovarian cancer are diagnosed late, which leads to significantly poorer prognosis. Apart from surgery and radiotherapy, a substantial number of ovarian cancer patients will undergo chemotherapy and platinum based agents are the mainstream first-line therapy for this disease. Despite the initial efficacy of these therapies, many women relapse; therefore, strategies for second-line therapies are required. Regulation of DNA transcription is crucial for tumour progression, metastasis and chemoresistance which offers potential for novel drug targets. Methods: We have reviewed the existing literature on the role of histone deacetylases, nuclear enzymes regulating gene transcription. Results and conclusion: Analysis of available data suggests that a signifant proportion of drug resistance stems from abberant gene expression, therefore HDAC inhibitors are amongst the most promising therapeutic targets for cancer treatment. Together with genetic testing, they may have a potential to serve as base for patient-adapted therapies

Problematic social media use: results from a large-scale nationally representative adolescent sample

Despite social media use being one of the most popular activities among adolescents, prevalence estimates among teenage samples of social media (problematic) use are lacking in the field. The present study surveyed a nationally representative Hungarian sample comprising 5,961 adolescents as part of the European School Survey Project on Alcohol and Other Drugs (ESPAD). Using the Bergen Social Media Addiction Scale (BSMAS) and based on latent profile analysis, 4.5% of the adolescents belonged to the at-risk group, and reported low self-esteem, high level of depression symptoms, and elevated social media use. Results also demonstrated that BSMAS has appropriate psychometric properties. It is concluded that adolescents at-risk of problematic social media use should be targeted by school-based prevention and intervention programs

Periodic eclipses of the young star PDS 110 discovered with WASP and KELT photometry

We report the discovery of eclipses by circumstellar disc material associated with the young star PDS 110 in the Ori OB1a association using the SuperWASP and Kilodegree Extremely Little Telescope surveys. PDS 110 (HD 290380, IRAS 05209-0107) is a rare Fe/Ge-type star, an similar to 10 Myr-old accreting intermediate-mass star showing strong infrared excess (L-IR/L-bol similar or equal to 0.25). Two extremely similar eclipses with a depth of 30 per cent and duration similar to 25 d were observed in 2008 November and 2011 January. We interpret the eclipses as caused by the same structure with an orbital period of 808 +/- 2 d. Shearing over a single orbit rules out diffuse dust clumps as the cause, favouring the hypothesis of a companion at similar to 2 au. The characteristics of the eclipses are consistent with transits by an unseen low-mass (1.8-70M(Jup)) planet or brown dwarf with a circumsecondary disc of diameter similar to 0.3 au. The next eclipse event is predicted to take place in 2017 September and could be monitored by amateur and professional observatories across the world

First Measurement of Transferred Polarization in the Exclusive e p --> e' K+ Lambda Reaction

The first measurements of the transferred polarization for the exclusive ep --> e'K+ Lambda reaction have been performed in Hall B at the Thomas Jefferson National Accelerator Facility using the CLAS spectrometer. A 2.567 GeV electron beam was used to measure the hyperon polarization over a range of Q2 from 0.3 to 1.5 (GeV/c)2, W from 1.6 to 2.15 GeV, and over the full center-of-mass angular range of the K+ meson. Comparison with predictions of hadrodynamic models indicates strong sensitivity to the underlying resonance contributions. A non-relativistic quark model interpretation of our data suggests that the s-sbar quark pair is produced with spins predominantly anti-aligned. Implications for the validity of the widely used 3P0 quark-pair creation operator are discussed.Comment: 6 pages, 4 figure