INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Genetic diversity among alfalfa genotypes (Medicago sativa L.) of non-dormant cultivars using SSR markers and agronomic traits

The aim of this study was to assess genetic diversity among 40 alfalfa (Medicago sativa L.) genotypes of different non-dormant (FD=8) cultivars. Biomass yield, regrowth speed and reaction to spring black stem, lepto leaf spot, and rust were evaluated. Analyses of variances were performed using a mixed model to examine the agronomic variation among individuals. A principal component analysis on standardized agronomic data was performed. Agronomic data were also used to calculate Gower's distance and UPGMA algorithm. For the molecular analysis, six SSR markers were evaluated and 84 alleles were identified. The genetic distance was estimated using standard Nei's distance. Average standard genetic diversity was 0.843, indicating a high degree of variability among genotypes. Finally, a generalized procrustes analysis was performed to calculate the correlation between molecular and agronomic distance, indicating a 65.4% of consensus. This value is likely related to the low number of individuals included in the study, which might have underestimated the real phenotypic variability among genotypes. Despite the low number of individuals and SSR markers analyzed, this study provides a baseline for future diversity studies to identify genetically distant alfalfa individuals or cultivars.

General descriptions, organic contents and ratios of geolipid components at DSDP Holes 75-530B and 75-532

Distributions of free and bound n-alkanes, n-alkanoic acids, and n-alkanols were determined in order to compare the character of organic matter contained in organic-carbon-rich sediments from two sites sampled by the hydraulic piston corer. Two diatomaceous debris-flow samples of Pleistocene age were obtained from Hole 530B in the Angola Basin. A sample of bioturbated Pleistocene diatomaceous clay and another of bioturbated late Miocene nannofossil clay were collected from Hole 532 on the Walvis Ridge. Geolipid distributions of all samples contain large terrigenous contributions and lesser amounts of marine components. Similarities in organic matter contents of Hole 530B and Hole 532 sediments suggest that a common depositional setting, probably on the Walvis Ridge, was the original source of these sediments through Quaternary, and possibly late Neogene, times and that downslope relocation of these biogenic deposits has frequently occurred.

Lithology, stage, organic carbon, d13C, hydrocarbons and fatty acids at DSDP Hole 93-603B

Organic-matter-rich Upper Cretaceous claystones from DSDP Hole 603B, lower continental rise, had organic carbon values ranging from 1.7 to 13.7%, C/N ratios from 32 to 72, and d13C values from -23.5 to -27.1 per mil. Lipid class maxima for the unbound alkanes (C29 and C31), unbound fatty acids (C28 and C30), and bound fatty acids (C24, C26 , and C28) and the strong odd-carbon and even-carbon preferences, respectively, suggested that the organic matter in these sediments was partially the result of input from continental plant waxes. Transport of the organic-matter-rich sediments to the deep sea from the near-shore environment probably resulted from turbiditic flow under oxygen-stressed conditions.

Novel application for exergy and thermoeconomic analysis of processes simulated with Aspen Plus

A visual basic application for Microsoft® Excel 2007 has been developed as a helpful tool to perform mass, energy, exergy and thermoeconomic (MHBT) calculations during the systematic analysis of energy processes simulated with Aspen Plus®. The application reads an Excel workbook containing three sheets with the matter, work and heat streams results of an Aspen Plus® simulation. The required information from the Aspen Plus® simulation and the algorithm/calculations of the application are described and applied to an Air Separation Unit (ASU). This application helps the designer when MHBT analyses are performed, as it increases the knowledge of the process simulated with Aspen Plus®. It’s a valuable tool not only because of the calculations performed, but also because it creates a new Excel workbook where the results and the formulae written on the cells are fully visible and editable. There is free access to the application and it has no protection allowing changes and improvements to be done.

Relating goal scheduling, precedence, and memory management in and-parallel execution of logic programs

The interactions among three important issues involved in the implementation of logic programs in parallel (goal scheduling, precedence, and memory management) are discussed. A simplified, parallel memory management model and an efficient, load-balancing goal scheduling strategy are presented. It is shown how, for systems which support "don't know" non-determinism, special care has to be taken during goal scheduling if the space recovery characteristics of sequential systems are to be preserved. A solution based on selecting only "newer" goals for execution is described, and an algorithm is proposed for efficiently maintaining and determining precedence relationships and variable ages across parallel goals. It is argued that the proposed schemes and algorithms make it possible to extend the storage performance of sequential systems to parallel execution without the considerable overhead previously associated with it. The results are applicable to a wide class of parallel and coroutining systems, and they represent an efficient alternative to "all heap" or "spaghetti stack" allocation models.

Present and Future of Photovoltaic Solar Electricity

Ponencia Invitada presentada en el IEEE Region 8 Student Branch and GOLD Congress

Personalized rule-based closed-loop control algorithm for type 1 diabetes

Type 1 diabetes-mellitus implies a life-threatening absolute insulin deficiency. Artificial pancreas (CGM sensor, insulin pump and control algorithm) is promising to outperform current open-loop therapies.

Deciduous tree reconstruction algorithm based on cylinder fitting from mobile terrestrial laser scanned point clouds

Vector reconstruction of objects from an unstructured point cloud obtained with a LiDAR-based system (light detection and ranging) is one of the most promising methods to build three dimensional models of orchards. The cylinder fitting method for woody structure reconstruction of leafless trees from point clouds obtained with a mobile terrestrial laser scanner (MTLS) has been analysed. The advantage of this method is that it performs reconstruction in a single step. The most time consuming part of the algorithm is generation of the cylinder direction, which must be recalculated at the inclusion of each point in the cylinder. The tree skeleton is obtained at the same time as the cluster of cylinders is formed. The method does not guarantee a unique convergence and the reconstruction parameter values must be carefully chosen. A balanced processing of clusters has also been defined which has proven to be very efficient in terms of processing time by following the hierarchy of branches, predecessors and successors. The algorithm was applied to simulated MTLS of virtual orchard models and to MTLS data of real orchards. The constraints applied in the method have been reviewed to ensure better convergence and simpler use of parameters. The results obtained show a correct reconstruction of the woody structure of the trees and the algorithm runs in linear logarithmic time

Uncoupling of transfer of the presequence and unfolding of the mature domain in precursor translocation across the mitochondrial outer membrane

Translocation of mitochondrial precursor proteins across the mitochondrial outer membrane is facilitated by the translocase of the outer membrane (TOM) complex. By using site-specific photocrosslinking, we have mapped interactions between TOM proteins and a mitochondrial precursor protein arrested at two distinct stages, stage A (accumulated at 0°C) and stage B (accumulated at 30°C), in the translocation across the outer membrane at high resolution not achieved previously. Although the stage A and stage B intermediates were assigned previously to the forms bound to the cis site and the trans site of the TOM complex, respectively, the results of crosslinking indicate that the presequence of the intermediates at both stage A and stage B is already on the trans side of the outer membrane. The mature domain is unfolded and bound to Tom40 at stage B whereas it remains folded at stage A. After dissociation from the TOM complex, translocation of the stage B intermediate, but not of the stage A intermediate, across the inner membrane was promoted by the intermembrane-space domain of Tom22. We propose a new model for protein translocation across the outer membrane, where translocation of the presequence and unfolding of the mature domain are not necessarily coupled.

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.

Reciprocal regulation of Gsα by palmitate and the βγ subunit

Hormonal activation of Gs, the stimulatory regulator of adenylyl cyclase, promotes dissociation of αs from Gβγ, accelerates removal of covalently attached palmitate from the Gα subunit, and triggers release of a fraction of αs from the plasma membrane into the cytosol. To elucidate relations among these three events, we assessed biochemical effects in vitro of attached palmitate on recombinant αs prepared from Sf9 cells. In comparison to the unpalmitoylated protein (obtained from cytosol of Sf9 cells, treated with a palmitoyl esterase, or expressed as a mutant protein lacking the site for palmitoylation), palmitoylated αs (from Sf9 membranes, 50% palmitoylated) was more hydrophobic, as indicated by partitioning into TX-114, and bound βγ with 5-fold higher affinity. βγ protected GDP-bound αs, but not αs· GTP[γS], from depalmitoylation by a recombinant esterase. We conclude that βγ binding and palmitoylation reciprocally potentiate each other in promoting membrane attachment of αs and that dissociation of αs·GTP from βγ is likely to mediate receptor-induced αs depalmitoylation and translocation of the protein to cytosol in intact cells.

The quorum-sensing transcriptional regulator TraR requires its cognate signaling ligand for protein folding, protease resistance, and dimerization

Complexes between the quorum-sensing regulator TraR and its inducing ligand autoinducer (AAI) are soluble in Escherichia coli, whereas apo-TraR is almost completely insoluble. Here we show that the lack of soluble TraR is due in large part to rapid proteolysis, inasmuch as apo-TraR accumulated to high levels in an E. coli strain deficient in Clp and Lon proteases. In pulse labeling experiments, AAI protected TraR against proteolysis only when it was added before the radiolabel. This observation indicates that TraR proteins can productively bind AAI only during their own synthesis on polysomes, whereas fully synthesized apo-TraR proteins are not functional AAI receptors. Purified apo-TraR was rapidly degraded by trypsin to oligopeptides, whereas TraR–AAI complexes were more resistant to trypsin and were cleaved at discrete interdomain linkers, indicating that TraR requires AAI to attain its mature tertiary structure. TraR–AAI complexes eluted from a gel filtration column as dimers and bound DNA as dimers. In contrast, apo-TraR was monomeric, and incubation with AAI under a variety of conditions did not cause dimerization. We conclude that AAI is critical for the folding of nascent TraR protein into its mature tertiary structure and that full-length apo-TraR cannot productively bind AAI and is consequently targeted for rapid proteolysis.

Switch from Myc/Max to Mad1/Max binding and decrease in histone acetylation at the telomerase reverse transcriptase promoter during differentiation of HL60 cells

Recent evidence suggests that the Myc and Mad1 proteins are implicated in the regulation of the gene encoding the human telomerase reverse transcriptase (hTERT), the catalytic subunit of telomerase. We have analyzed the in vivo interaction between endogenous c-Myc and Mad1 proteins and the hTERT promoter in HL60 cells with the use of the chromatin immunoprecipitation assay. The E-boxes at the hTERT proximal promoter were occupied in vivo by c-Myc in exponentially proliferating HL60 cells but not in cells induced to differentiate by DMSO. In contrast, Mad1 protein was induced and bound to the hTERT promoter in differentiated HL60 cells. Concomitantly, the acetylation of the histones at the promoter was significantly reduced. These data suggest that the reciprocal E-box occupancy by c-Myc and Mad1 is responsible for activation and repression of the hTERT gene in proliferating and differentiated HL60 cells, respectively. Furthermore, the histone deacetylase inhibitor trichostatin A inhibited deacetylation of histones at the hTERT promoter and attenuated the repression of hTERT transcription during HL60 cell differentiation. In addition, trichostatin A treatment activated hTERT transcription in resting human lymphocytes and fibroblasts. Taken together, these results indicate that acetylation/deacetylation of histones is operative in the regulation of hTERT expression.

Comparative study of RPSALG algorithm for convex semi-infinite programming

The Remez penalty and smoothing algorithm (RPSALG) is a unified framework for penalty and smoothing methods for solving min-max convex semi-infinite programing problems, whose convergence was analyzed in a previous paper of three of the authors. In this paper we consider a partial implementation of RPSALG for solving ordinary convex semi-infinite programming problems. Each iteration of RPSALG involves two types of auxiliary optimization problems: the first one consists of obtaining an approximate solution of some discretized convex problem, while the second one requires to solve a non-convex optimization problem involving the parametric constraints as objective function with the parameter as variable. In this paper we tackle the latter problem with a variant of the cutting angle method called ECAM, a global optimization procedure for solving Lipschitz programming problems. We implement different variants of RPSALG which are compared with the unique publicly available SIP solver, NSIPS, on a battery of test problems.