Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

Multi-site evaluation of an urban land-surface model: intra-urban heterogeneity, seasonality and parameter complexity requirements

An extensive off-line evaluation of the Noah/Single Layer Urban Canopy Model (Noah/SLUCM) urban land-surface model is presented using data from 15 sites to assess (1) the ability of the scheme to reproduce the surface energy balance observed in a range of urban environments, including seasonal changes, and (2) the impact of increasing complexity of input parameter information. Model performance is found to be most dependent on representation of vegetated surface area cover; refinement of other parameter values leads to smaller improvements. Model biases in net all-wave radiation and trade-offs between turbulent heat fluxes are highlighted using an optimization algorithm. Here we use the Urban Zones to characterize Energy partitioning (UZE) as the basis to assign default SLUCM parameter values. A methodology (FRAISE) to assign sites (or areas) to one of these categories based on surface characteristics is evaluated. Using three urban sites from the Basel Urban Boundary Layer Experiment (BUBBLE) dataset, an independent evaluation of the model performance with the parameter values representative of each class is performed. The scheme copes well with both seasonal changes in the surface characteristics and intra-urban heterogeneities in energy flux partitioning, with RMSE performance comparable to similar state-of-the-art models for all fluxes, sites and seasons. The potential of the methodology for high-resolution atmospheric modelling application using the Weather Research and Forecasting (WRF) model is highlighted. This analysis supports the recommendations that (1) three classes are appropriate to characterize the urban environment, and (2) that the parameter values identified should be adopted as default values in WRF.

Dysregulation of REST-regulated coding and non-coding RNAs in a cellular model of Huntington's disease

Huntingtin (Htt) protein interacts with many transcriptional regulators, with widespread disruption to the transcriptome in Huntington's disease (HD) brought about by altered interactions with the mutant Htt (muHtt) protein. Repressor Element-1 Silencing Transcription Factor (REST) is a repressor whose association with Htt in the cytoplasm is disrupted in HD, leading to increased nuclear REST and concomitant repression of several neuronal-specific genes, including brain-derived neurotrophic factor (Bdnf). Here, we explored a wide set of HD dysregulated genes to identify direct REST targets whose expression is altered in a cellular model of HD but that can be rescued by knock-down of REST activity. We found many direct REST target genes encoding proteins important for nervous system development, including a cohort involved in synaptic transmission, at least two of which can be rescued at the protein level by REST knock-down. We also identified several microRNAs (miRNAs) whose aberrant repression is directly mediated by REST, including miR-137, which has not previously been shown to be a direct REST target in mouse. These data provide evidence of the contribution of inappropriate REST-mediated transcriptional repression to the widespread changes in coding and non-coding gene expression in a cellular model of HD that may affect normal neuronal function and survival.

Rescue of gene expression by modified REST decoy oligonucleotides in a cellular model of Huntington's disease

Transcriptional dysfunction is a prominent hallmark of Huntington's disease (HD). Several transcription factors have been implicated in the aetiology of HD progression and one of the most prominent is repressor element 1 (RE1) silencing transcription factor (REST). REST is a global repressor of neuronal gene expression and in the presence of mutant Huntingtin increased nuclear REST levels lead to elevated RE1 occupancy and a concomitant increase in target gene repression, including brain-derived neurotrophic factor. It is of great interest to devise strategies to reverse transcriptional dysregulation caused by increased nuclear REST and determine the consequences in HD. Thus far, such strategies have involved RNAi or mutant REST constructs. Decoys are double-stranded oligodeoxynucleotides corresponding to the DNA-binding element of a transcription factor and act to sequester it, thereby abrogating its transcriptional activity. Here, we report the use of a novel decoy strategy to rescue REST target gene expression in a cellular model of HD. We show that delivery of the decoy in cells expressing mutant Huntingtin leads to its specific interaction with REST, a reduction in REST occupancy of RE1s and rescue of target gene expression, including Bdnf. These data point to an alternative strategy for rebalancing the transcriptional dysregulation in HD.

Transcriptional dysregulation of coding and non-coding genes in cellular models of Huntington's disease

HD (Huntington's disease) is a late onset heritable neurodegenerative disorder that is characterized by neuronal dysfunction and death, particularly in the cerebral cortex and medium spiny neurons of the striatum. This is followed by progressive chorea, dementia and emotional dysfunction, eventually resulting in death. HD is caused by an expanded CAG repeat in the first exon of the HD gene that results in an abnormally elongated polyQ (polyglutamine) tract in its protein product, Htt (Huntingtin). Wild-type Htt is largely cytoplasmic; however, in HD, proteolytic N-terminal fragments of Htt form insoluble deposits in both the cytoplasm and nucleus, provoking the idea that mutHtt (mutant Htt) causes transcriptional dysfunction. While a number of specific transcription factors and co-factors have been proposed as mediators of mutHtt toxicity, the causal relationship between these Htt/transcription factor interactions and HD pathology remains unknown. Previous work has highlighted REST [RE1 (repressor element 1)-silencing transcription factor] as one such transcription factor. REST is a master regulator of neuronal genes, repressing their expression. Many of its direct target genes are known or suspected to have a role in HD pathogenesis, including BDNF (brain-derived neurotrophic factor). Recent evidence has also shown that REST regulates transcription of regulatory miRNAs (microRNAs), many of which are known to regulate neuronal gene expression and are dysregulated in HD. Thus repression of miRNAs constitutes a second, indirect mechanism by which REST can alter the neuronal transcriptome in HD. We will describe the evidence that disruption to the REST regulon brought about by a loss of interaction between REST and mutHtt may be a key contributory factor in the widespread dysregulation of gene expression in HD.

Intra-crystalline protein diagenesis (IcPD) in Patella vulgata. Part II: Breakdown and temperature sensitivity

Artificial diagenesis of the intra-crystalline proteins isolated from Patella vulgata was induced by isothermal heating at 140 °C, 110 °C and 80 °C. Protein breakdown was quantified for multiple amino acids, measuring the extent of peptide bond hydrolysis, amino acid racemisation and decomposition. The patterns of diagenesis are complex; therefore the kinetic parameters of the main reactions were estimated by two different methods: 1) a well-established approach based on fitting mathematical expressions to the experimental data, e.g. first-order rate equations for hydrolysis and power-transformed first-order rate equations for racemisation; and 2) an alternative model-free approach, which was developed by estimating a “scaling” factor for the independent variable (time) which produces the best alignment of the experimental data. This method allows the calculation of the relative reaction rates for the different temperatures of isothermal heating. High-temperature data were compared with the extent of degradation detected in sub-fossil Patella specimens of known age, and we evaluated the ability of kinetic experiments to mimic diagenesis at burial temperature. The results highlighted a difference between patterns of degradation at low and high temperature and therefore we recommend caution for the extrapolation of protein breakdown rates to low burial temperatures for geochronological purposes when relying solely on kinetic data.

MODIS-based vegetation index has sufficient sensitivity to indicate stand-level intra-seasonal climatic stress in oak and beech forests

Context: Variation in photosynthetic activity of trees induced by climatic stress can be effectively evaluated using remote sensing data. Although adverse effects of climate on temperate forests have been subjected to increased scrutiny, the suitability of remote sensing imagery for identification of drought stress in such forests has not been explored fully. Aim: To evaluate the sensitivity of MODIS-based vegetation index to heat and drought stress in temperate forests, and explore the differences in stress response of oaks and beech. Methods: We identified 8 oak and 13 beech pure and mature stands, each covering between 4 and 13 MODIS pixels. For each pixel, we extracted a time series of MODIS NDVI from 2000 to 2010. We identified all sequences of continuous unseasonal NDVI decline to be used as the response variable indicative of environmental stress. Neural Networks-based regression modelling was then applied to identify the climatic variables that best explain observed NDVI declines. Results: Tested variables explained 84–97% of the variation in NDVI, whilst air temperature-related climate extremes were found to be the most influential. Beech showed a linear response to the most influential climatic predictors, while oak responded in a unimodal pattern suggesting a better coping mechanism. Conclusions: MODIS NDVI has proved sufficiently sensitive as a stand-level indicator of climatic stress acting upon temperate broadleaf forests, leading to its potential use in predicting drought stress from meteorological observations and improving parameterisation of forest stress indices.

Using standby generators in support of an electricity network: logistics and experience

Electricity load shifting is becoming a big topic in the world of ‘green’ retail. Marks & Spencer (M&S) aim to become the world’s most sustainable retailer (1) and part of that commitment means contributing to the future electricity network. While intelligent operation of fridges and Heating, Ventilation and Air Conditioning (HVAC) systems are a wide area of research, standby generators should be considered too, as they are the most widely adopted form of distributed generation. In this paper, the experience of using standby generators in Northern Ireland to support the grid is shared and the logistics of future projects are discussed. Interactions with maintenance schedules, electricity costs, grid code, staffing and store opening times are discussed as well as the financial implications associated with running generators for grid support.

Intra-tumor heterogeneity: lessons from microbial evolution and clinical implications

Multiple subclonal populations of tumor cells can coexist within the same tumor. This intra-tumor heterogeneity will have clinical implications and it is therefore important to identify factors that drive or suppress such heterogeneous tumor progression. Evolutionary biology can provide important insights into this process. In particular, experimental evolution studies of microbial populations, which exist as clonal populations that can diversify into multiple subclones, have revealed important evolutionary processes driving heterogeneity within a population. There are transferrable lessons that can be learnt from these studies that will help us to understand the process of intra-tumor heterogeneity in the clinical setting. In this review, we summarize drivers of microbial diversity that have been identified, such as mutation rate and environmental influences, and discuss how knowledge gained from microbial experimental evolution studies may guide us to identify and understand important selective factors that promote intra-tumor heterogeneity. Furthermore, we discuss how these factors could be used to direct and optimize research efforts to improve patient care, focusing on therapeutic resistance. Finally, we emphasize the need for longitudinal studies to address the impact of these potential tumor heterogeneity-promoting factors on drug resistance, metastatic potential and clinical outcome.

Simple method for sub-diffraction resolution imaging of cellular structures on standard confocal microscopes by three-photon absorption of quantum dots

This study describes a simple technique that improves a recently developed 3D sub-diffraction imaging method based on three-photon absorption of commercially available quantum dots. The method combines imaging of biological samples via tri-exciton generation in quantum dots with deconvolution and spectral multiplexing, resulting in a novel approach for multi-color imaging of even thick biological samples at a 1.4 to 1.9-fold better spatial resolution. This approach is realized on a conventional confocal microscope equipped with standard continuous-wave lasers. We demonstrate the potential of multi-color tri-exciton imaging of quantum dots combined with deconvolution on viral vesicles in lentivirally transduced cells as well as intermediate filaments in three-dimensional clusters of mouse-derived neural stem cells (neurospheres) and dense microtubuli arrays in myotubes formed by stacks of differentiated C2C12 myoblasts.

Comparison of bioactivities, binding properties and intra-follicular levels of bovine follistatins

Five isoforms of follistatin (FST) (Mr 31, 33, 35, 37, 41kDa) were purified from bovine follicular fluid (bFF). Comparison of their activin- and heparan sulphate proteoglycan (HSP)-binding properties and bio-potencies in neutralization of activin-A action in vitro revealed that all five isoforms bound activin-A, but with different affinities. Only the 31kDa isoform (FST-288) bound to HSP. FST-288 also showed the greatest biopotency with 35 and 41kDa isoforms being least potent. To determine whether bovine follicle development is associated with changing intrafollicular FST and activin profiles, we analyzed bFF from dominant (DF) and subordinate (SF) follicles collected at strategic times during a synchronized estrous cycle. Total FST, activin-A and activin-AB were measured by immunoassay while individual FST isoforms were quantified by immunoblotting. Follicle diameter was positively correlated with estrogen:progesterone ratio (r=0.56) in bFF but negatively correlated with activin-A (r=-0.34), activin-AB (r=-0.80) and ‘total’ FST (r=-0.70) levels. Follicle diameter was positively correlated with abundance of the 41 kDa isoform (r=0.59) but negatively correlated with abundance of 33 and 31 kDa isoforms (r=-0.56, -0.41). Both follicle status (DF vs SF) and cycle stage affected total FST, activin-A, activin-B levels while follicle status, but not cycle stage, affected abundance of 41, 37, 33 and 31kDa FST isoforms. Collectively, these findings indicate that intrafollicular FST isoforms that differ in their ability to bind and neutralise activins and associate with cell-surface proteoglycans, show divergent changes during follicle development. Enhanced FST production may have an important negative role, either directly or via inhibition of the positive effects of activins, on follicle growth and function during follicular waves.

Stabilization of the i-motif structure by the intra-strand cross-link formation

The i-motif structures are formed by oligonucleotides containing cytosine tracts under acidic conditions. The folding of the i-motif under physiological conditions is of great interest because of its biological role. In this study, we investigated the effect of the intra-strand cross-link on the stability of the i-motif structure. The 4-vinyl-substituted analog of thymidine (T-vinyl) was incorporated into the 5′-end of the human telomere complementary strand, which formed the intra-strand cross-link with the internal adenine. The intra-strand cross-linked i-motif displayed CD spectra similar to that of the natural i-motif at acidic pH, which was transformed into a random coil with the increasing pH. The pH midpoint for the transition from the i-motif to random coil increased from pH 6.1 for the natural one to pH 6.8 for the cross-linked one. The thermodynamic parameters were obtained by measuring the thermal melting behaviors by CD and UV, and it was determined that the intra-strand cross-linked i-motif is stabilized due to a favorable entropy effect. Thus, this study has clearly indicated the validity of the intra-strand cross-linking for stabilization of the i-motif structure.

Differential epigenetic reprogramming in response to specific endocrine therapies promotes cholesterol biosynthesis and cellular invasion

Endocrine therapies target the activation of the oestrogen receptor alpha (ERα) via distinct mechanisms, but it is not clear whether breast cancer cells can adapt to treatment using drug-specific mechanisms. Here we demonstrate that resistance emerges via drug-specific epigenetic reprogramming. Resistant cells display a spectrum of phenotypical changes with invasive phenotypes evolving in lines resistant to the aromatase inhibitor (AI). Orthogonal genomics analysis of reprogrammed regulatory regions identifies individual drug-induced epigenetic states involving large topologically associating domains (TADs) and the activation of super-enhancers. AI-resistant cells activate endogenous cholesterol biosynthesis (CB) through stable epigenetic activation in vitro and in vivo. Mechanistically, CB sparks the constitutive activation of oestrogen receptors alpha (ERα) in AI-resistant cells, partly via the biosynthesis of 27-hydroxycholesterol. By targeting CB using statins, ERα binding is reduced and cell invasion is prevented. Epigenomic-led stratification can predict resistance to AI in a subset of ERα-positive patients

Cellular mechanisms of cardiac hypertrophy.

Hypertrophy of myocytes in the heart ventricles is an important adaptation that in vivo occurs in response to a requirement for increased contractile power. It involves changes at the level of gene transcription, stimulation of the rate of protein synthesis (translation), and increased assembly of myofibrils. There is mounting evidence of the involvement of reversible protein phosphorylation and dephosphorylation in most of these processes. Protein kinase C, mitogen-activated protein kinases, and transcription factors have been implicated in the modulation of the transcriptional changes. Activation of translation may also be mediated through protein phosphorylation/dephosphorylation, although this has not been clearly established in the heart. Here we provide a critical overview of the signalling pathways involved in the hypertrophic response and provide a scheme to account for many of its features.

Regulation of gene and protein expression in cardiac myocyte hypertrophy and apoptosis

Considerable efforts have been expended in elucidating the inter-cellular and intra-cellular signaling pathways which elicit cardiac myocyte hypertrophy or apoptosis, and in identifying the changes which are associated with the end-stage of the response. The challenge now is to link the two. Although some of the signaling effects will be the acute modulation of existing protein function, long-term effects which bring about and maintain the hypertrophic state or which culminate in cell death are mediated at the level of gene and protein expression. With the advances in micro-array technology and genome sequencing, it is now possible to obtain a picture of the global gene expression profile in myocytes or in whole heart which dictates the proteins which could be made. This is not the final picture since additional regulation at the level of translation modulates the relative proportions of each protein that can be made from the transcriptome. Even here, further regulation of protein stability and turnover means that ultimately it is still necessary to examine the proteome to determine what may cause the functional changes in a cell. Thus, in order to gain a full picture of events which regulate the response and gain some insight into possible points of intervention for therapy, it is necessary to examine gene expression, mRNA translation and protein expression in concert.