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 Intersection of Religious Autonomy and Religious Symbols: Setting the Stage

Two sets of issues in the area of law and religion have generated a large share of attention and controversy across a wide number of countries and jurisdictions in recent years. The first set of issues relates to the autonomy of churches and other religiously affiliated entities such as schools and social service organisations in their hiring and personnel decisions, involving the question of how far, if at all, such entities should be free from the influence and oversight of the state. The second set of issues involves the presence of religious symbols in the public sphere, such as in state schools or on public lands, involving the question of how far the state should be free from the influence of religion. Although these issues – freedom of religion from the state, and freedom of the state from religion – could be viewed as opposite sides of the same coin, they are almost always treated as separate lines of inquiry, and the implications of each for the other have not been the subject of much scrutiny. In this Introduction, we consider whether insights might be drawn from thinking about these issues both from a comparative law perspective and also from considering these two lines of cases together.

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.

TAP Reactor Development and Application in the Kinetic Characterization of Catalysts for Heterogeneously Catalyzed Gas Phase Reactions

This manuscript describes the application and further development of the TAP technique in kinetic characterization of heterogeneous catalysis. The major application of TAP systems is to study mechanisms, kinetics and transport phenomena in heterogeneous catalysis, all of which is made possible by the sub-millisecond time resolution. Furthermore, the kinetic information obtained can be used to gain an insight into the mechanism occurring over the catalyst system. This is advantageous as heterogeneous catalysts with an improved efficiency can be developed as a result. TAP kinetic studies are carried out at low pressure (~1x10-7 mbar) and TAP pulses are sufficiently small (1013-1015 molecules) so as to maintain this low pressure. The use of a small number of molecules in comparison to the total number of active sites means the state of the catalyst remains relatively unchanged. The use of the low intensity pulses also makes the pressure gradient negligible and so allows the TAP reactor system to operate in the Knudsen Diffusion regime, where gas-gas reactions are eliminated. Hence only gas-catalyst reactions are investigated and, by the use of moment analysis of observed exit flow, rate constants of elementary steps of the reaction can be obtained.

In this manuscript, two attempts to further the TAP technique are reported. Firstly, the work undertaken at QUB to attempt to control the number of molecules of condensable reagents that can be pulsed during a TAP pulse experiment is disclosed. Secondly, a collaborative project with SAI Ltd Manchester is discussed in a separate chapter, where technical details and validation of a customised time of flight mass spectrometer (ToF MS) for the QUB TAP-1 system are reported. A collaborative project with Cardiff Catalysis Institute focusing on the study of CO oxidation over hopcalite catalysts is also reported. The analysis of the experimental results has provided an insight into the possible mechanism of the oxidation of CO over these catalysts. A correction function has also been derived which accounts for the adsorption of reactant molecules over inert materials that are used for the reactor packing in TAP experiments. This function was then applied to the selective reduction of O2 in a H2 rich ethene feed, so that more accurate TAP moment based analysis could be conducted.

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.

High Performance Real-time Pre-Processing for Fixed-Complexity Sphere Decoder

Pre-processing (PP) of received symbol vector and channel matrices is an essential pre-requisite operation for Sphere Decoder (SD)-based detection of Multiple-Input Multiple-Output (MIMO) wireless systems. PP is a highly complex operation, but relative to the total SD workload it represents a relatively small fraction of the overall computational cost of detecting an OFDM MIMO frame in standards such as 802.11n. Despite this, real-time PP architectures are highly inefficient, dominating the resource cost of real-time SD architectures. This paper resolves this issue. By reorganising the ordering and QR decomposition sub operations of PP, we describe a Field Programmable Gate Array (FPGA)-based PP architecture for the Fixed Complexity Sphere Decoder (FSD) applied to 4 × 4 802.11n MIMO which reduces resource cost by 50% as compared to state-of-the-art solutions whilst maintaining real-time performance.

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.