A Grid Supercomputing Enviroment for High Demand Computational Applications

Despite the huge increase in processor and interprocessor network performace, many computational problems remain unsolved due to lack of some critical resources such as floating point sustained performance, memory bandwidth, etc... Examples of these problems are found in areas of climate research, biology, astrophysics, high energy physics (montecarlo simulations) and artificial intelligence, among others. For some of these problems, computing resources of a single supercomputing facility can be 1 or 2 orders of magnitude apart from the resources needed to solve some them. Supercomputer centers have to face an increasing demand on processing performance, with the direct consequence of an increasing number of processors and systems, resulting in a more difficult administration of HPC resources and the need for more physical space, higher electrical power consumption and improved air conditioning, among other problems. Some of the previous problems can´t be easily solved, so grid computing, intended as a technology enabling the addition and consolidation of computing power, can help in solving large scale supercomputing problems. In this document, we describe how 2 supercomputing facilities in Spain joined their resources to solve a problem of this kind. The objectives of this experience were, among others, to demonstrate that such a cooperation can enable the solution of bigger dimension problems and to measure the efficiency that could be achieved. In this document we show some preliminary results of this experience and to what extend these objectives were achieved.

On the Relationship between Market Power and Bank Risk Taking

We analyse risk-taking behaviour of banks in the context of spatial competition. Banks mobilise unsecured deposits by offering deposit rates, which they invest either in a prudent or a gambling asset. Limited liability along with high return of a successful gamble induce moral hazard at the bank level. We show that when the market power is low, banks invest in the gambling asset. On the other hand, for sufficiently high levels of market power, all banks choose the prudent asset to invest in. We further show that a merger of two neighboring banks increases the likelihood of prudent behaviour. Finally, introduction of a deposit insurance scheme exacerbates banks’ moral hazard problem.

Peroxisomal polyhydroxyalkanoate biosynthesis is a promising strategy for bioplastic production in high biomass crops.

Polyhydroxyalkanoates (PHAs) are bacterial carbon storage polymers with diverse plastic-like properties. PHA biosynthesis in transgenic plants is being developed as a way to reduce the cost and increase the sustainability of industrial PHA production. The homopolymer polyhydroxybutyrate (PHB) is the simplest form of these biodegradable polyesters. Plant peroxisomes contain the substrate molecules and necessary reducing power for PHB biosynthesis, but peroxisomal PHB production has not been explored in whole soil-grown transgenic plants to date. We generated transgenic sugarcane (Saccharum sp.) with the three-enzyme Ralstonia eutropha PHA biosynthetic pathway targeted to peroxisomes. We also introduced the pathway into Arabidopsis thaliana, as a model system for studying and manipulating peroxisomal PHB production. PHB, at levels up to 1.6%-1.8% dry weight, accumulated in sugarcane leaves and A. thaliana seedlings, respectively. In sugarcane, PHB accumulated throughout most leaf cell types in both peroxisomes and vacuoles. A small percentage of total polymer was also identified as the copolymer poly (3-hydroxybutyrate-co-3-hydroxyvalerate) in both plant species. No obvious deleterious effect was observed on plant growth because of peroxisomal PHA biosynthesis at these levels. This study highlights how using peroxisomal metabolism for PHA biosynthesis could significantly contribute to reaching commercial production levels of PHAs in crop plants.

A robust genomic signature for the detection of colorectal cancer patients with microsatellite instability phenotype and high mutation frequency.

Microsatellite instability (MSI) occurs in 10-20% of colorectal tumours and is associated with good prognosis. Here we describe the development and validation of a genomic signature that identifies colorectal cancer patients with MSI caused by DNA mismatch repair deficiency with high accuracy. Microsatellite status for 276 stage II and III colorectal tumours has been determined. Full-genome expression data was used to identify genes that correlate with MSI status. A subset of these samples (n = 73) had sequencing data for 615 genes available. An MSI gene signature of 64 genes was developed and validated in two independent validation sets: the first consisting of frozen samples from 132 stage II patients; and the second consisting of FFPE samples from the PETACC-3 trial (n = 625). The 64-gene MSI signature identified MSI patients in the first validation set with a sensitivity of 90.3% and an overall accuracy of 84.8%, with an AUC of 0.942 (95% CI, 0.888-0.975). In the second validation, the signature also showed excellent performance, with a sensitivity 94.3% and an overall accuracy of 90.6%, with an AUC of 0.965 (95% CI, 0.943-0.988). Besides correct identification of MSI patients, the gene signature identified a group of MSI-like patients that were MSS by standard assessment but MSI by signature assessment. The MSI-signature could be linked to a deficient MMR phenotype, as both MSI and MSI-like patients showed a high mutation frequency (8.2% and 6.4% of 615 genes assayed, respectively) as compared to patients classified as MSS (1.6% mutation frequency). The MSI signature showed prognostic power in stage II patients (n = 215) with a hazard ratio of 0.252 (p = 0.0145). Patients with an MSI-like phenotype had also an improved survival when compared to MSS patients. The MSI signature was translated to a diagnostic microarray and technically and clinically validated in FFPE and frozen samples.

Sentinel lymph node involvement and a high Breslow index are independent factors of risk for early relapse of melanoma.

AIM: The clinical relevance of sentinel lymph node (SLN) analysis was evaluated prospectively and compared with other known risk factors of relapse in early stage melanoma. METHODS: Surgery was guided by lymphoscintigraphy, blue dye and gamma probe detection. SLN were analysed by haematoxylin eosin (HE) histochemistry and multimarker immunohistochemistry (IHC). Disease free survival (DFS) was evaluated with Kaplan-Meier plots according to different parameters and Cox analyses of variance. RESULTS: From 210 patients a total of 381 SLN were excised. Lymphoscintigraphy identified all excised SLN with only 2 false positive lymphatic lakes. Fifty patients (24%) had tumour positive SLN. With a mean follow-up of 31.3 months, 29 tumour recurrences were observed, 19 (38%) in 50 SLN positive and 10 (6%) in 160 SLN negative patients. Strong predictive factors for early relapse (p < 0.0005) were SLN positivity and a high Breslow index. CONCLUSION: SLN tumour positivity is an independent factor of high risk for early relapse with a higher power of discrimination than the Breslow index.

Power laws and scaling of rain events and dry spells in the Catalonia region

We analyze the statistics of rain-event sizes, rain-event durations, and dry-spell durations in a network of 20 rain gauges scattered in an area situated close to the NW Mediterranean coast. Power-law distributions emerge clearly for the dryspell durations, with an exponent around 1.50 ± 0.05, although for event sizes and durations the power-law ranges are rather limited, in some cases. Deviations from power-law behavior are attributed to finite-size effects. A scaling analysis helps to elucidate the situation, providing support for the existence of scale invariance in these distributions. It is remarkable that rain data of not very high resolution yield findings in agreement with self-organized critical phenomena.

Simultaneous B(0)- and B(1)+-map acquisition for fast localized shim, frequency, and RF power determination in the heart at 3 T.

High-field (>or=3 T) cardiac MRI is challenged by inhomogeneities of both the static magnetic field (B(0)) and the transmit radiofrequency field (B(1)+). The inhomogeneous B fields not only demand improved shimming methods but also impede the correct determination of the zero-order terms, i.e., the local resonance frequency f(0) and the radiofrequency power to generate the intended local B(1)+ field. In this work, dual echo time B(0)-map and dual flip angle B(1)+-map acquisition methods are combined to acquire multislice B(0)- and B(1)+-maps simultaneously covering the entire heart in a single breath hold of 18 heartbeats. A previously proposed excitation pulse shape dependent slice profile correction is tested and applied to reduce systematic errors of the multislice B(1)+-map. Localized higher-order shim correction values including the zero-order terms for frequency f(0) and radiofrequency power can be determined based on the acquired B(0)- and B(1)+-maps. This method has been tested in 7 healthy adult human subjects at 3 T and improved the B(0) field homogeneity (standard deviation) from 60 Hz to 35 Hz and the average B(1)+ field from 77% to 100% of the desired B(1)+ field when compared to more commonly used preparation methods.

Predictors of HIV-protection behaviour in HIV-positive men who have sex with casual male partners: a test of the explanatory power of an extended Information-Motivation-Behavioural Skills model.

This prospective study applies an extended Information-Motivation-Behavioural Skills (IMB) model to establish predictors of HIV-protection behaviour among HIV-positive men who have sex with men (MSM) during sex with casual partners. Data have been collected from anonymous, self-administered questionnaires and analysed by using descriptive and backward elimination regression analyses. In a sample of 165 HIV-positive MSM, 82 participants between the ages of 23 and 78 (M=46.4, SD=9.0) had sex with casual partners during the three-month period under investigation. About 62% (n=51) have always used a condom when having sex with casual partners. From the original IMB model, only subjective norm predicted condom use. More important predictors that increased condom use were low consumption of psychotropics, high satisfaction with sexuality, numerous changes in sexual behaviour after diagnosis, low social support from friends, alcohol use before sex and habitualised condom use with casual partner(s). The explanatory power of the calculated regression model was 49% (p<0.001). The study reveals the importance of personal and social resources and of routines for condom use, and provides information for the research-based conceptualisation of prevention offers addressing especially people living with HIV ("positive prevention").

High bi-atrial organization in patients with long-standing persistent atrial fibrillation terminated within the left atrium

Sustained atrial fibrillation (AF) is maintained by sites displaying high dominant frequency (DF). In patients (pts) with long-standing persistent AF (LS-pAF), their spatial distribution and the presence of a left-to-right atrial DF gradient remain poorly known. We hypothesized that the pre-ablation bi-atrial frequency characteristics of LS-pAF pts terminated within the left atrium (LT) are different from that of non terminated (NT) ones. Methods: 23 consecutive pts (59±7y, LS-pAF duration 19±12m) underwent stepwise catheter ablation (step-CA) consisting in pulmonary veins isolation, left atrial (LA) defragmentation, and right atrial (RA) ablations for non terminated AF. A quadripolar catheter (CAT) was placed into the RA appendage (RAA), a decapolar CAT into the coronary sinus (CS) and a duodecapolar CAT into the LA divided into 8 segments. For each segment, 20-sec of bipolar recording was acquired. The DF was defined as the largest peak in the power spectrum (3-15 Hz). The inter-atrial DF gradient was defined as the DF difference between LA and RA appendages. Results: LS-pAF was terminated in 83% (19/23) of the pts: 17 LT, 2 during RA ablation and 4 NT. The figure shows that before ablation bi-atrial DF values of LT pts are significantly lower than that of NT pts for each LA segment as well as for the RAA (p < 0.05). No significant LA-to-RA DF gradient was observed both for LT (0.3±0.5 Hz, p=ns) and NT (0.5±0.03 Hz, p=ns) pts. No significant difference in DF values was observed between LA segments. Conclusions: The lower DF of LT pts is suggestive of a higher organization within both atria compared to NT pts. Our findings suggest that low bi-atrial DF values, but not inter-atrial DF gradient, might be of interest for selecting LS-pAF candidates for sinus rhythm restoration by step-CA.

Application for fault location in electrical power distribution systems

Fault location has been studied deeply for transmission lines due to its importance in power systems. Nowadays the problem of fault location on distribution systems is receiving special attention mainly because of the power quality regulations. In this context, this paper presents an application software developed in Matlabtrade that automatically calculates the location of a fault in a distribution power system, starting from voltages and currents measured at the line terminal and the model of the distribution power system data. The application is based on a N-ary tree structure, which is suitable to be used in this application due to the highly branched and the non- homogeneity nature of the distribution systems, and has been developed for single-phase, two-phase, two-phase-to-ground, and three-phase faults. The implemented application is tested by using fault data in a real electrical distribution power system

Ultra high throughput sequencing in human DNA variation detection: a comparative study on the NDUFA3-PRPF31 region.

BACKGROUND: Ultra high throughput sequencing (UHTS) technologies find an important application in targeted resequencing of candidate genes or of genomic intervals from genetic association studies. Despite the extraordinary power of these new methods, they are still rarely used in routine analysis of human genomic variants, in part because of the absence of specific standard procedures. The aim of this work is to provide human molecular geneticists with a tool to evaluate the best UHTS methodology for efficiently detecting DNA changes, from common SNPs to rare mutations. METHODOLOGY/PRINCIPAL FINDINGS: We tested the three most widespread UHTS platforms (Roche/454 GS FLX Titanium, Illumina/Solexa Genome Analyzer II and Applied Biosystems/SOLiD System 3) on a well-studied region of the human genome containing many polymorphisms and a very rare heterozygous mutation located within an intronic repetitive DNA element. We identify the qualities and the limitations of each platform and describe some peculiarities of UHTS in resequencing projects. CONCLUSIONS/SIGNIFICANCE: When appropriate filtering and mapping procedures are applied UHTS technology can be safely and efficiently used as a tool for targeted human DNA variations detection. Unless particular and platform-dependent characteristics are needed for specific projects, the most relevant parameter to consider in mainstream human genome resequencing procedures is the cost per sequenced base-pair associated to each machine.

Simultaneous EEG-fMRI at ultra-high field: Artifact prevention and safety assessment.

The simultaneous recording of scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can provide unique insights into the dynamics of human brain function, and the increased functional sensitivity offered by ultra-high field fMRI opens exciting perspectives for the future of this multimodal approach. However, simultaneous recordings are susceptible to various types of artifacts, many of which scale with magnetic field strength and can seriously compromise both EEG and fMRI data quality in recordings above 3T. The aim of the present study was to implement and characterize an optimized setup for simultaneous EEG-fMRI in humans at 7T. The effects of EEG cable length and geometry for signal transmission between the cap and amplifiers were assessed in a phantom model, with specific attention to noise contributions from the MR scanner coldheads. Cable shortening (down to 12cm from cap to amplifiers) and bundling effectively reduced environment noise by up to 84% in average power and 91% in inter-channel power variability. Subject safety was assessed and confirmed via numerical simulations of RF power distribution and temperature measurements on a phantom model, building on the limited existing literature at ultra-high field. MRI data degradation effects due to the EEG system were characterized via B0 and B1(+) field mapping on a human volunteer, demonstrating important, although not prohibitive, B1 disruption effects. With the optimized setup, simultaneous EEG-fMRI acquisitions were performed on 5 healthy volunteers undergoing two visual paradigms: an eyes-open/eyes-closed task, and a visual evoked potential (VEP) paradigm using reversing-checkerboard stimulation. EEG data exhibited clear occipital alpha modulation and average VEPs, respectively, with concomitant BOLD signal changes. On a single-trial level, alpha power variations could be observed with relative confidence on all trials; VEP detection was more limited, although statistically significant responses could be detected in more than 50% of trials for every subject. Overall, we conclude that the proposed setup is well suited for simultaneous EEG-fMRI at 7T.

Exploring Cross-layer power management for PGAS applications on the SCC platform

Technological limitations and power constraints are resulting in high-performance parallel computing architectures that are based on large numbers of high-core-count processors. Commercially available processors are now at 8 and 16 cores and experimental platforms, such as the many-core Intel Single-chip Cloud Computer (SCC) platform, provide much higher core counts. These trends are presenting new sets of challenges to HPC applications including programming complexity and the need for extreme energy efficiency.In this work, we first investigate the power behavior of scientific PGAS application kernels on the SCC platform, and explore opportunities and challenges for power management within the PGAS framework. Results obtained via empirical evaluation of Unified Parallel C (UPC) applications on the SCC platform under different constraints, show that, for specific operations, the potential for energy savings in PGAS is large; and power/performance trade-offs can be effectively managed using a cross-layerapproach. We investigate cross-layer power management using PGAS language extensions and runtime mechanisms that manipulate power/performance tradeoffs. Specifically, we present the design, implementation and evaluation of such a middleware for application-aware cross-layer power management of UPC applications on the SCC platform. Finally, based on our observations, we provide a set of recommendations and insights that can be used to support similar power management for PGAS applications on other many-core platforms.

A survey on voltage sag events in power systems

This paper aims to survey the techniques and methods described in literature to analyse and characterise voltage sags and the corresponding objectives of these works. The study has been performed from a data mining point of view

Unusual voltage sag event detection in power systems

Three multivariate statistical tools (principal component analysis, factor analysis, analysis discriminant) have been tested to characterize and model the sags registered in distribution substations. Those models use several features to represent the magnitude, duration and unbalanced grade of sags. They have been obtained from voltage and current waveforms. The techniques are tested and compared using 69 registers of sags. The advantages and drawbacks of each technique are listed