Discovering subgroups of patients from DNA copy number data using NMF on compacted matrices

In the study of complex genetic diseases, the identification of subgroups of patients sharing similar genetic characteristics represents a challenging task, for example, to improve treatment decision. One type of genetic lesion, frequently investigated in such disorders, is the change of the DNA copy number (CN) at specific genomic traits. Non-negative Matrix Factorization (NMF) is a standard technique to reduce the dimensionality of a data set and to cluster data samples, while keeping its most relevant information in meaningful components. Thus, it can be used to discover subgroups of patients from CN profiles. It is however computationally impractical for very high dimensional data, such as CN microarray data. Deciding the most suitable number of subgroups is also a challenging problem. The aim of this work is to derive a procedure to compact high dimensional data, in order to improve NMF applicability without compromising the quality of the clustering. This is particularly important for analyzing high-resolution microarray data. Many commonly used quality measures, as well as our own measures, are employed to decide the number of subgroups and to assess the quality of the results. Our measures are based on the idea of identifying robust subgroups, inspired by biologically/clinically relevance instead of simply aiming at well-separated clusters. We evaluate our procedure using four real independent data sets. In these data sets, our method was able to find accurate subgroups with individual molecular and clinical features and outperformed the standard NMF in terms of accuracy in the factorization fitness function. Hence, it can be useful for the discovery of subgroups of patients with similar CN profiles in the study of heterogeneous diseases.

Power and Energy Implications of the Number of Threads Used on the Intel Xeon Phi

Energy consumption has become an important area of research of late. With the advent of new manycore processors, situations have arisen where not all the processors need to be active to reach an optimal relation between performance and energy usage. In this paper a study of the power and energy usage of a series of benchmarks, the PARSEC and the SPLASH- 2X Benchmark Suites, on the Intel Xeon Phi for different threads configurations, is presented. To carry out this study, a tool was designed to monitor and record the power usage in real time during execution time and afterwards to compare the r

Mapping neutral, ion and electron number densities within laser ablated plasma plumes

Spatially and temporally varying neutral, ion and electron number densities have been mapped out within laser ablated plasma plumes expanding into vacuum. Ablation of a magnesium target was performed using a KrF laser, 30 ns pulse duration and 248 nm wavelength. During the initial stage of plasma expansion (t <EQ 100 ns) interferometry has been used to obtain line averaged electron number densities, for laser power densities on target in the range 1.3 - 3.0 X 108 W/cm2. Later in the plasma expansion (t equals 1 microsecond(s) ) simultaneous absorption and laser induced fluorescence spectroscopy has been used to determine 3D neutral and ion number densities, for a power density equal to 6.7 X 107 W/cm2. Two distinct regions within the plume were identified. One is a fast component (approximately 106 cm-1) consisting of ions and neutrals with maximum number densities observed to be approximately 30 and 4 X 1012 cm-3 respectively, and the second consists of slow moving neutral material at a number density of up to 1015 cm-3. Additionally a Langmuir probe has been used to obtain ion and electron number densities at very late times in the plasma expansion (1 microsecond(s) <EQ t <EQ 15 microsecond(s) ). A copper target was ablated using a Nd:YAG laser, 7.5 ns duration and 532 nm (2 (omega) ) wavelength, with a power density on target equal to 6 X 108 W/cm2. Two regions within the plume with different velocities were observed. Within a fast component (approximately 3 X 106 cms-1) electron and ion number densities of the order 5 X 1012 cm-3 were observed and within the second slower component (approximately 106 cms-1) electron and ion number densities of the order 1 - 2 X 1013 cm-3 were determined.

Distributed Massive MIMO in Cellular Networks: Impact of Imperfect Hardware and Number of Oscillators

Distributed massive multiple-input multiple-output (MIMO) combines the array gain of coherent MIMO processing with the proximity gains of distributed antenna setups. In this paper, we analyze how transceiver hardware impairments affect the downlink with maximum ratio transmission. We derive closed-form spectral efficiencies expressions and study their asymptotic behavior as the number of the antennas increases. We prove a scaling law on the hardware quality, which reveals that massive MIMO is resilient to additive distortions, while multiplicative phase noise is a limiting factor. It is also better to have separate oscillators at each antenna than one per BS.

The paradox of invasion: Reeves’ muntjac deer invade the British Isles from a limited number of founding females

High levels of genetic diversity and high propagule pressure are favoured by conservation biologists as the basis for successful reintroductions and ensuring the persistence of populations. However, invasion ecologists recognize the ‘paradox of invasion’, as successful species introductions may often be characterized by limited numbers of individuals and associated genetic bottlenecks. In the present study, we used a combination of high-resolution nuclear and mitochondrial genetic markers to investigate the invasion history of Reeves' muntjac deer in the British Isles. This invasion has caused severe economic and ecological damage, with secondary spread currently a concern throughout Europe and potentially globally. Microsatellite analysis based on eight loci grouped all 176 introduced individuals studied from across the species' range in the UK into one genetic cluster, and seven mitochondrial D-loop haplotypes were recovered, two of which were present at very low frequency and were related to more common haplotypes. Our results indicate that the entire invasion can be traced to a single founding event involving a low number of females. These findings highlight the fact that even small releases of species may, if ignored, result in irreversible and costly invasion, regardless of initial genetic diversity or continual genetic influx.

Magnetic Characteristics of a High-layer-number NiFe/FeMn Multilayer

We report the static & dynamic magnetic characteristics of a high-layer-number NiFe/FeMn multilayer test structure with potential applications in broadband absorber and filter devices. To allow fine control over the absorption linewidths and to understand the mechanisms governing the resonances in a tailored structure similar to that expected to be used in real world applications, the multilayer was intentionally designed to have layer thickness and interface roughness variations. Magnetometry measurements show the sample has complex hysteresis loops with features consistent with single ferromagnetic film reversals. Structural characterisation by transmission electron microscopy allows us to correlate the magnetic properties with structural features. Analysis of resonance frequencies from broadband ferromagnetic resonance measurements as a function of field magnitude and orientation provide values of the local exchange bias, rotatable anisotropy, and uniaxial anisotropy fields for specific layers in the stack and explain the observed mode softening. The linewidths of the multilayer are adjustable around the bias field, approaching twice that seen at larger fields, allowing control over the bandwidth of devices formed from the structure.

Low alanine aminotransferase levels and higher number of cardiovascular events in people with Type 2 diabetes: analysis of the Fenofibrate Intervention and Event Lowering in Diabetes (FIELD) study.

AIMS: To determine whether alanine aminotransferase or gamma-glutamyltransferase levels, as markers of liver health and non-alcoholic fatty liver disease, might predict cardiovascular events in people with Type 2 diabetes.

METHODS: Data from the Fenofibrate Intervention and Event Lowering in Diabetes study were analysed to examine the relationship between liver enzymes and incident cardiovascular events (non-fatal myocardial infarction, stroke, coronary and other cardiovascular death, coronary or carotid revascularization) over 5 years.

RESULTS: Alanine aminotransferase level had a linear inverse relationship with the first cardiovascular event occurring in participants during the study period. After adjustment, for every 1 sd higher baseline alanine aminotransferase value (13.2 U/l), the risk of a cardiovascular event was 7% lower (95% CI 4-13; P=0.02). Participants with alanine aminotransferase levels below and above the reference range 8-41 U/l for women and 9-59 U/l for men, had hazard ratios for a cardiovascular event of 1.86 (95% CI 1.12-3.09) and 0.65 (95% CI 0.49-0.87), respectively (P=0.001). No relationship was found for gamma-glutamyltransferase.

CONCLUSIONS: The data may indicate that in people with Type 2 diabetes, which is associated with higher alanine aminotransferase levels because of prevalent non-alcoholic fatty liver disease, a low alanine aminotransferase level is a marker of hepatic or systemic frailty rather than health. This article is protected by copyright. All rights reserved.

Integration of copy number and transcriptomics provides risk stratification in prostate cancer: A discovery and validation cohort study

BACKGROUND: Understanding the heterogeneous genotypes and phenotypes of prostate cancer is fundamental to improving the way we treat this disease. As yet, there are no validated descriptions of prostate cancer subgroups derived from integrated genomics linked with clinical outcome.

METHODS: In a study of 482 tumour, benign and germline samples from 259 men with primary prostate cancer, we used integrative analysis of copy number alterations (CNA) and array transcriptomics to identify genomic loci that affect expression levels of mRNA in an expression quantitative trait loci (eQTL) approach, to stratify patients into subgroups that we then associated with future clinical behaviour, and compared with either CNA or transcriptomics alone.

FINDINGS: We identified five separate patient subgroups with distinct genomic alterations and expression profiles based on 100 discriminating genes in our separate discovery and validation sets of 125 and 103 men. These subgroups were able to consistently predict biochemical relapse (p = 0.0017 and p = 0.016 respectively) and were further validated in a third cohort with long-term follow-up (p = 0.027). We show the relative contributions of gene expression and copy number data on phenotype, and demonstrate the improved power gained from integrative analyses. We confirm alterations in six genes previously associated with prostate cancer (MAP3K7, MELK, RCBTB2, ELAC2, TPD52, ZBTB4), and also identify 94 genes not previously linked to prostate cancer progression that would not have been detected using either transcript or copy number data alone. We confirm a number of previously published molecular changes associated with high risk disease, including MYC amplification, and NKX3-1, RB1 and PTEN deletions, as well as over-expression of PCA3 and AMACR, and loss of MSMB in tumour tissue. A subset of the 100 genes outperforms established clinical predictors of poor prognosis (PSA, Gleason score), as well as previously published gene signatures (p = 0.0001). We further show how our molecular profiles can be used for the early detection of aggressive cases in a clinical setting, and inform treatment decisions.

INTERPRETATION: For the first time in prostate cancer this study demonstrates the importance of integrated genomic analyses incorporating both benign and tumour tissue data in identifying molecular alterations leading to the generation of robust gene sets that are predictive of clinical outcome in independent patient cohorts.

The Application of Geophysics to a Number of Threats to Irish Soils

In its 2002 Communication "Towards a Thematic Strategy on Soil Protection" (COM 2002, 179), the Commission identified the main threats to which soils in the EU are confronted. More recently, the EU published a Thematic Strategy for Soil Protection (COM 2006, 231) as well as a proposal for a Soil Framework Directive (COM 2006, 232). The proposed directive, if implemented, will require states to identify a number of specific soil degradation processes that have occurred or are likely to occur in the near future. This project applies a number of recent geophysical developments, in addition to a range of traditional approaches, to a number of areas of Irish concern to Thematic Strategy for Soil Protection and the proposed Soil Framework Directive.

A generalization of the Hoffman-Lovász upper bound on the independence number of a regular graph

A family of quadratic programming problems whose optimal values are upper bounds on the independence number of a graph is introduced. Among this family, the quadratic programming problem which gives the best upper bound is identified. Also the proof that the upper bound introduced by Hoffman and Lovász for regular graphs is a particular case of this family is given. In addition, some new results characterizing the class of graphs for which the independence number attains the optimal value of the above best upper bound are given. Finally a polynomial-time algorithm for approximating the size of the maximum independent set of an arbitrary graph is described and the computational experiments carried out on 36 DIMACS clique benchmark instances are reported.