Primer design and marker clustering for multiplex SNP-IT primer extension genotyping assay using statistical modeling

MOTIVATION: The optimization of the primer design is critical for the development of high-throughput SNP genotyping methods. Recently developed statistical models of the SNP-IT primer extension genotyping reaction allow further improvement of primer quality for the assay. RESULTS: Here we describe how the statistical models can be used to improve primer design for the assay. We also show how to optimize clustering of the SNP markers into multiplex panels using statistical model for multiplex SNP-IT. The primer set failure probability calculated by a model is used as a minimization function for both primer selection and primers clustering. Three clustering algorithms for the multiplex genotyping SNP-IT assay are described and their relative performance is evaluated. We also describe the approaches to improve the speed of primer design and clustering calculations when using the statistical models. Our clustering decreases the average failure probability of the marker set by 7-25%. The experimental marker failure rate in the multiplex reaction was reduced dramatically and success rate can be achieved as high as 96%. AVAILABILITY: The primer design using statistical models is freely available from www.autoprimer.com.     

Modeling the relationship between LVAD support time and gene expression changes in the human heart by penalized partial least squares

MOTIVATION: Heart failure affects more than 20 million people in the world. Heart transplantation is the most effective therapy, but the number of eligible patients far outweighs the number of available donor hearts. The left mechanical ventricular assist device (LVAD) has been developed as a successful substitution therapy that aids the failing ventricle while a patient is waiting for the donor heart. We obtained genomics data from paired human heart samples harvested at the time of LVAD implant and explant. The heart failure patients in our study were supported by the LVAD for various periods of time. The goal of this study is to model the relationship between the time of LVAD support and gene expression changes. RESULTS: To serve the purpose, we propose a novel penalized partial least squares (PPLS) method to build a regression model. Compared with partial least squares and Breiman's random forest method, PPLS gives the best prediction results for the LVAD data.     

Cellware--a multi-algorithmic software for computational systems biology

The intracellular environment of a cell hosts a wide variety of enzymatic reactions, diffusion events, molecular binding, polymerization and metabolic channeling. To transform these biological events into a computational framework, distinct modeling strategies are required. While currently no tool is capable of capturing all these events, progress is being made to create an integrated environment for the modeling community. To address this niche requirement, Cellware has been developed to offer a multi-algorithmic environment for modeling and simulating both deterministic and stochastic events in the cell. AVAILABILITY: The software is available for free and can be downloaded from http://www.bii.a-star.edu.sg/sbg/cellware     

Proteolytic degradation of the nuclear isoform of uracil-DNA glycosylase occurs during the S phase of the cell cycle

Uracil-DNA glycosylases are enzymes that remove uracil from DNA and initiate base-excision repair. These enzymes play a key role in maintaining genomic integrity by reducing the mutagenic events caused by G:C to A:T transition mutations. The recent finding that a family of RNA editing enzymes (APOBECs) can deaminate cytosine in DNA has raised the interest in these base-excision repair enzymes. This research focuses on the regulation of the nuclear isoform of uracil-DNA glycosylase, a 36000 Da protein that contains a unique 44 amino acid N-terminus. In synchronized HeLa cells, UDG1A protein levels decrease to barely detectable levels during the S phase of the cell cycle. Immunoblot analysis of immunoprecipitated or affinity-isolated UDG1A reveals ubiquitin-conjugated UDG1A when proteolysis is inhibited using N-acetyl-leu-leu-norleu-al or MG132, inhibitors of proteosomal dependent protein degradation. Transient transfection experiments, with histidine-tagged ubiquitin, were used to confirm that endogenous UDG1A is ubiquitinated in vivo. Addition of the nuclear export inhibitor, leptomycin B, prevents ubiquitination and degradation of UDG1A. This indicates that translocation from the nucleus may be a step in UDG1A turnover. Finally, UDG1A protein degradation is prevented when cells are incubated with the cyclin-dependent kinase inhibitor, roscovitine. These results suggest that protein phosphorylation and/or nuclear export participate in the post-translational regulation of UDG1A protein levels.     

Determination of NG,NG-dimethyl-L-arginine, an endogenous NO synthase inhibitor, by gas chromatography-mass spectrometry

A fully validated gas chromatographic-mass spectrometric (GC-MS) method for the accurate and precise quantification of NG,NG-dimethyl-L-arginine (asymmetric dimethylarginine, ADMA), an endogenous inhibitor of the NO synthase, in cell culture supernatants and in small volumes of plasma is described. ADMA was concentrated by solid phase extraction and converted to its methyl ester pentafluoropropionic amide derivative. The derivatives were analyzed without any further purification. Using gas chromatography-chemical ionization mass spectrometry, fragment ions at m/z 634 and m/z 640 were obtained for ADMA and for NG,NG-[2H6]-dimethyl-L-arginine ([2H6]-ADMA) as internal standard, respectively. [2H6]-ADMA was synthesized by reaction of L-ornithine fastened at bromcyan-agarose with dimethylamine. The limit of detection of the method was 2 fmol, while the limit of quantitation for cell culture supernatants was 0.05 microM. The method was validated in a concentration range of 0-1.2 microM in cell culture medium and 0-2 microM in 50 microl aliquots of human plasma. The precision was > or =97% and the accuracy was determined to be > or =94%. This method is fast, rugged and an alternative to high performance liquid chromatography (HPLC) analysis of ADMA in cell culture supernatants and small volumes of human plasma.     

Proteomic analysis of an orthotopic neuroblastoma xenograft animal model

Neuroblastoma is the most common extracranial solid tumour of childhood and comprises up to 50% of malignancies among infants. There is a great need of designing novel therapeutic strategies and proteome analysis is one approach for defining markers useful for tumour diagnosis, as well as molecular targets for novel experimental therapies. We started by comparing healthy adrenal glands (which are the election organs developing primary neuroblastoma, NB, tumours) and adrenal glands carrying primary NB tumours, taken from nude mice. Standard maps of healthy and tumour samples were generated by analysis with the PDQuest software. The comparison between such maps showed up- and down-regulation of 84 polypeptide chains, out of a total of 700 spots detected by a fluorescent stain, Sypro Ruby. Spots that were differentially expressed between the two groups, were analysed by MALDI-TOF mass spectrometry and 14 of these spots were identified so far. Among these proteins, of particular interest are the down-regulated proteins adrenodoxin (21-folds), carbonic anhydrase III (eight-folds) and aldose reductase related protein I (eight-folds), as well as the up-regulated protein peptidyl-propyl cis-trans isomerase A (five-folds). Moreover new proteins, which were absent in control samples, were expressed in tumour samples, such as nucleophosmin (NPM) and stathmin (oncoprotein 18).     

Two-dimensional gel electrophoresis database of murine R1 embryonic stem cells

Embryonic stem (ES) cell lines represent a population of undifferentiated pluripotent cells capable of multilineage differentiation in vitro. Although very useful for studying developmental processes, human ES cell lines have also been suggested as a potential and unlimited source for cellular transplantation and the treatment of human disease. The proteomic basis of embryonic stemness (pluripotentiality and multilineage differentiation) and the transitions that lead to specific cell lineages however, remain to be defined. As an important first step in defining these processes, we have performed a proteomic analysis of undifferentiated mouse R1 ES cell lines using pH 3-10, 4-7 and 6-11 two-dimensional electrophoresis gels, matrix-assisted laser desorption/ionization and tandem mass spectrometry. Of the 700 gel spots analyzed, 241 distinct protein species were identified corresponding to 218 unique proteins, with a significant proportion functionally related to protein expression.     

Remotely Sensed Interannual Variations and Trends in Terrestrial Net Primary Productivity 1981–2000

Spatial and temporal variations in net primary production (NPP) are of great importance to ecological studies, natural resource management, and terrestrial carbon sink estimates. However, most of the existing estimates of interannual variation in NPP at regional and global scales were made at coarse resolutions with climate-driven process models. In this study, we quantified global NPP variation at an 8 km and 10-day resolution from 1981 to 2000 based on satellite observations. The high resolution was achieved using the GLObal Production Efficiency Model (GLO-PEM), which was driven with variables derived almost entirely from satellite remote sensing. The results show that there was an increasing trend toward enhanced terrestrial NPP that was superimposed on high seasonal and interannual variations associated with climate variability and that the increase was occurring in both northern and tropical latitudes. NPP generally decreased in El Niño season and increased in La Niña seasons, but the magnitude and spatial pattern of the response varied widely between individual events. Our estimates also indicate that the increases in NPP during the period were caused mainly by increases in atmospheric carbon dioxide and precipitation. The enhancement of NPP by warming was limited to northern high latitudes (above 50°N); in other regions, the interannual variations in NPP were correlated negatively with temperature and positively with precipitation.