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.

Differential expression of mRNAs encoding co-receptor (betaglycan) and activin type-I preovulatory and prehierarchical follicles of the putative inhibin and type-II receptors in the laying hen ovary

Ovarian follicle development is primarily regulated by an interplay between the pituitary gonadotrophins, LH and FSH, and ovary-derived steroids. Increasing evidence implicates regulatory roles of transforming growth factor-beta (TGF beta) superfamily members, including inhibins and activins. The aim of this study was to identify the expression of mRNAs encoding key receptors of the inhibin/activin system in ovarian follicles ranging from 4 mm in diameter to the dominant F1 follicle (similar to 40 turn). Ovaries were collected (n=16) from inid-sequence hens maintained on a long-day photoschedule (16h of light:8 h of darkness). All follicles removed were dissected into individual granulosa and thecal layers. RNA was extracted and cDNA synthesized. Real-time quantitative PCR was used to quantify the expression of niRNA encoding betaglycan, activin receptor (ActR) subtypes (type-I, -IIA and -IIB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH); receptor expression data were normalized to GAPDH expression. Detectable levels of ActRI, -IIA and -IIB and the inhibin co-receptor (betaglycan) expression were found in all granulosa and thecal layers analysed. Granulosa ActRI mRNA peaked (P < 0(.)05) in 8-9(.)9 mm follicles, whereas ActRIIA rose significantly from 6-7(.)9 mm to 8-9(.)9 nun, before filling to F3/2; levels then rose sharply (3-fold) to F1 levels. Granulosa betaglycan niRNA expression rose 3-fold from 4-5(.)9 min to 8-9(.)9 mm, before falling 4-fold to F3/2; levels then rose sharply (4-fold) to F1 levels. ActRIIB levels did not vary significantly during follicular development. Thecal ActRI mRNA expression was similar from 4-7(.)9 mm then decreased significantly to a nadir at the F4 position, before increasing 2-fold to the F1 (P < 0(.)05). Although thecal ActRIIB and -IIA expression did not vary significantly from 4 nim to F3, ActRIIB expression increased significantly (2-fold) from F3 to F1 and ActIIA, increased 22-fold from F2 to F1 (P < 0(.)05). Thecal betaglycan fell to a nadir at F6 after follicle selection; levels then increased significantly to F2, before filling similar to 50% in the F I. In all follicles studied expression of betaglycan and ActRI (granulosa: 1-0(.)65, P < 0-001, n=144/group; theca: r=0(.)49, P < 0-001, n=144/group) was well correlated. No significant correlations were identified between betaglycan and ActRIIA or -IIB. Considering all follicles analysed, granulosa mRNA expression of betaglycan, ActRI ActRIIA and ActRIIB were all significantly lower than in corresponding thecal tissue (betaglycan, 11(.)4-fold; ActRIIB, 5(.)1-fold; ActR(.) 3-8-fold: ActRIIA, 2(.)8-fold). The co-localization of type-I and -II activin receptors and betaglycan on granulosa and thecal cells are consistent with a local auto/paracrine role of inhibins and activins in modulating ovarian follicle development, selection and progression in the domestic fowl.

Rapid parallel expression in E-coli and insect cells: Analysis of five lef gene products of the Autographa californica multiple nuclear polyhedrosis virus (AcMNPV)

A number of strategies are emerging for the high throughput (HTP) expression of recombinant proteins to enable structural and functional study. Here we describe a workable HTP strategy based on parallel protein expression in E. coli and insect cells. Using this system we provide comparative expression data for five proteins derived from the Autographa californica polyhedrosis virus genome that vary in amino acid composition and in molecular weight. Although the proteins are part of a set of factors known to be required for viral late gene expression, the precise function of three of the five, late expression factors (lefs) 6, 7 and 10, is unknown. Rapid expression and characterisation has allowed the determination of their ability to bind DNA and shown a cellular location consistent with their properties. Our data point to the utility of a parallel expression strategy to rapidly obtain workable protein expression levels from many open reading frames (ORFs).

Differential expression of mRNAs encoding co-receptor (betaglycan) and activin type-I preovulatory and prehierarchical follicles of the putative inhibin and type-II receptors in the laying hen ovary

Ovarian follicle development is primarily regulated by an interplay between the pituitary gonadotrophins, LH and FSH, and ovary-derived steroids. Increasing evidence implicates regulatory roles of transforming growth factor-beta (TGF beta) superfamily members, including inhibins and activins. The aim of this study was to identify the expression of mRNAs encoding key receptors of the inhibin/activin system in ovarian follicles ranging from 4 mm in diameter to the dominant F1 follicle (similar to 40 turn). Ovaries were collected (n=16) from inid-sequence hens maintained on a long-day photoschedule (16h of light:8 h of darkness). All follicles removed were dissected into individual granulosa and thecal layers. RNA was extracted and cDNA synthesized. Real-time quantitative PCR was used to quantify the expression of niRNA encoding betaglycan, activin receptor (ActR) subtypes (type-I, -IIA and -IIB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH); receptor expression data were normalized to GAPDH expression. Detectable levels of ActRI, -IIA and -IIB and the inhibin co-receptor (betaglycan) expression were found in all granulosa and thecal layers analysed. Granulosa ActRI mRNA peaked (P < 0(.)05) in 8-9(.)9 mm follicles, whereas ActRIIA rose significantly from 6-7(.)9 mm to 8-9(.)9 nun, before filling to F3/2; levels then rose sharply (3-fold) to F1 levels. Granulosa betaglycan niRNA expression rose 3-fold from 4-5(.)9 min to 8-9(.)9 mm, before falling 4-fold to F3/2; levels then rose sharply (4-fold) to F1 levels. ActRIIB levels did not vary significantly during follicular development. Thecal ActRI mRNA expression was similar from 4-7(.)9 mm then decreased significantly to a nadir at the F4 position, before increasing 2-fold to the F1 (P < 0(.)05). Although thecal ActRIIB and -IIA expression did not vary significantly from 4 nim to F3, ActRIIB expression increased significantly (2-fold) from F3 to F1 and ActIIA, increased 22-fold from F2 to F1 (P < 0(.)05). Thecal betaglycan fell to a nadir at F6 after follicle selection; levels then increased significantly to F2, before filling similar to 50% in the F I. In all follicles studied expression of betaglycan and ActRI (granulosa: 1-0(.)65, P < 0-001, n=144/group; theca: r=0(.)49, P < 0-001, n=144/group) was well correlated. No significant correlations were identified between betaglycan and ActRIIA or -IIB. Considering all follicles analysed, granulosa mRNA expression of betaglycan, ActRI ActRIIA and ActRIIB were all significantly lower than in corresponding thecal tissue (betaglycan, 11(.)4-fold; ActRIIB, 5(.)1-fold; ActR(.) 3-8-fold: ActRIIA, 2(.)8-fold). The co-localization of type-I and -II activin receptors and betaglycan on granulosa and thecal cells are consistent with a local auto/paracrine role of inhibins and activins in modulating ovarian follicle development, selection and progression in the domestic fowl.

Multiplexed quantitative proteomics using differential metabolic 15N-labeling

There are several advantages of using metabolic labeling in quantitative proteomics. The early pooling of samples compared to post-labeling methods eliminates errors from different sample processing, protein extraction and enzymatic digestion. Metabolic labeling is also highly efficient and relatively inexpensive compared to commercial labeling reagents. However, methods for multiplexed quantitation in the MS-domain (or ‘non-isobaric’ methods), suffer from signal dilution at higher degrees of multiplexing, as the MS/MS signal for peptide identification is lower given the same amount of peptide loaded onto the column or injected into the mass spectrometer. This may partly be overcome by mixing the samples at non-uniform ratios, for instance by increasing the fraction of unlabeled proteins. We have developed an algorithm for arbitrary degrees of nonisobaric multiplexing for relative protein abundance measurements. We have used metabolic labeling with different levels of 15N, but the algorithm is in principle applicable to any isotope or combination of isotopes. Ion trap mass spectrometers are fast and suitable for LC-MS/MS and peptide identification. However, they cannot resolve overlapping isotopic envelopes from different peptides, which makes them less suitable for MS-based quantitation. Fourier-transform ion cyclotron resonance (FTICR) mass spectrometry is less suitable for LC-MS/MS, but provides the resolving power required to resolve overlapping isotopic envelopes. We therefore combined ion trap LC-MS/MS for peptide identification with FTICR LC-MS for quantitation using chromatographic alignment. We applied the method in a heat shock study in a plant model system (A. thaliana) and compared the results with gene expression data from similar experiments in literature.

Transcription profiling in human platelets reveals LRRFIP1 as a novel protein regulating platelet function

Within the healthy population, there is substantial, heritable, and interindividual variability in the platelet response. We explored whether a proportion of this variability could be accounted for by interindividual variation in gene expression. Through a correlative analysis of genome-wide platelet RNA expression data from 37 subjects representing the normal range of platelet responsiveness within a cohort of 500 subjects, we identified 63 genes in which transcript levels correlated with variation in the platelet response to adenosine diphosphate and/or the collagen-mimetic peptide, cross-linked collagen-related peptide. Many of these encode proteins with no reported function in platelets. An association study of 6 of the 63 genes in 4235 cases and 6379 controls showed a putative association with myocardial infarction for COMMD7 (COMM domain-containing protein 7) and a major deviation from the null hypo thesis for LRRFIP1 [leucine-rich repeat (in FLII) interacting protein 1]. Morpholino-based silencing in Danio rerio identified a modest role for commd7 and a significant effect for lrrfip1 as positive regulators of thrombus formation. Proteomic analysis of human platelet LRRFIP1-interacting proteins indicated that LRRFIP1 functions as a component of the platelet cytoskeleton, where it interacts with the actin-remodeling proteins Flightless-1 and Drebrin. Taken together, these data reveal novel proteins regulating the platelet response.

Genome-wide mapping of 5-hydroxymethylcytosine in three rice cultivars reveals its preferential localization in transcriptionally silent transposable element genes

5-Hydroxymethylcytosine (5hmC), a modified form of cytosine that is considered the sixth nucleobase in DNA, has been detected in mammals and is believed to play an important role in gene regulation. In this study, 5hmC modification was detected in rice by employing a dot-blot assay, and its levels was further quantified in DNA from different rice tissues using liquid chromatography-multistage mass spectrometry (LC-MS/MS/MS). The results showed large intertissue variation in 5hmC levels. The genome-wide profiles of 5hmC modification in three different rice cultivars were also obtained using a sensitive chemical labelling followed by a next-generation sequencing method. Thousands of 5hmC peaks were identified, and a comparison of the distributions of 5hmC among different rice cultivars revealed the specificity and conservation of 5hmC modification. The identified 5hmC peaks were significantly enriched in heterochromatin regions,and mainly located in transposable element (TE) genes, especially around retrotransposons. The correlation analysis of 5hmC and gene expression data revealed a close association between 5hmC and silent TEs. These findings provide a resource for plant DNA 5hmC epigenetic studies and expand our knowledge of 5hmC modification.

Adaptation of Rhizobium leguminosarum to pea, alfalfa and sugar beet rhizospheres investigated by comparative transcriptomics

Background The rhizosphere is the microbe-rich zone around plant roots and is a key determinant of the biosphere's productivity. Comparative transcriptomics was used to investigate general and plant-specific adaptations during rhizosphere colonization. Rhizobium leguminosarum biovar viciae was grown in the rhizospheres of pea (its legume nodulation host), alfalfa (a non-host legume) and sugar beet (non-legume). Gene expression data were compared to metabolic and transportome maps to understand adaptation to the rhizosphere. Results Carbon metabolism was dominated by organic acids, with a strong bias towards aromatic amino acids, C1 and C2 compounds. This was confirmed by induction of the glyoxylate cycle required for C2 metabolism and gluconeogenesis in all rhizospheres. Gluconeogenesis is repressed in R. leguminosarum by sugars, suggesting that although numerous sugar and putative complex carbohydrate transport systems are induced in the rhizosphere, they are less important carbon sources than organic acids. A common core of rhizosphere-induced genes was identified, of which 66% are of unknown function. Many genes were induced in the rhizosphere of the legumes, but not sugar beet, and several were plant specific. The plasmid pRL8 can be considered pea rhizosphere specific, enabling adaptation of R. leguminosarum to its host. Mutation of many of the up-regulated genes reduced competitiveness for pea rhizosphere colonization, while two genes specifically up-regulated in the pea rhizosphere reduced colonization of the pea but not alfalfa rhizosphere. Conclusions Comparative transcriptome analysis has enabled differentiation between factors conserved across plants for rhizosphere colonization as well as identification of exquisite specific adaptation to host plants.

Tumour necrosis factor alpha induces rapid reduction in AMPA receptor-mediated calcium entry in motor neurones by increasing cell surface expression of the GluR2 subunit: relevance to neurodegeneration

The AMPA receptor (AMPAR) subunit GluR2, which regulates excitotoxicity and the inflammatory cytokine tumour necrosis factor alpha (TNF alpha) have both been implicated in motor neurone vulnerability in Amyotrophic Lateral Sclerosis/Motor Neurone Disease. TNF alpha has been reported to increase cell surface expression of AMPAR subunits to increase synaptic strength and enhance excitotoxicity, but whether this mechanism occurs in motor neurones is unknown. We used primary cultures of mouse motor neurones and cortical neurones to examine the interaction between TNF alpha receptor activation, GluR2 availability, AMPAR-mediated calcium entry and susceptibility to excitotoxicity. Short exposure to a physiologically relevant concentration of TNFalpha (10 ng/ml, 15 min) caused a marked redistribution of both GluR1 and GluR2 to the cell surface as determined by cell surface biotinylation and immunofluorescence. Using Fura-2 AM microfluorimetry we showed that exposure to TNFalpha caused a rapid reduction in the peak amplitude of AMPA-mediated calcium entry in a PI3-kinase and p38 kinase-dependent manner, consistent with increased insertion of GluR2-containing AMPAR into the plasma membrane. This resulted in a protection of motor neurones against kainate-induced cell death. Our data therefore, suggests that TNF alpha acts primarily as a physiological regulator of synaptic activity in motor neurones rather than a pathological drive in ALS

Toll-like receptor expression in murine DC subsets: lack of TLR7 expression by CD8α+ DC correlates with unresponsiveness to imidazoquinolines

Toll-like receptors (TLR) recognize microbial and viral patterns and activate dendritic cells (DC). TLR distribution among human DC subsets is heterogeneous: plasmacytoid DC (PDC) express TLR1, 7 and 9, while other DC types do not express TLR9 but express other TLR. Here, we report that mRNA for most TLR is expressed at similar levels by murine splenic DC sub-types, including PDC, but that TLR3 is preferentially expressed by CD8α+ DC while TLR5 and TLR7 are selectively absent from the same subset. Consistent with the latter, TLR7 ligand activates CD8α– DC and PDC, but not CD8α+ DC as measured by survival ex vivo, up-regulation of surface markers and production of IL-12p40. These data suggest that the dichotomy in TLR expression between plasmacytoid and non-plasmacytoid DC is not conserved between species. However, lack of TLR7 expression could restrict the involvement of CD8α+ DC in recognition of certain mouse pathogens.

Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Dietary alpha-tocopherol affects differential gene expression in rat testes

Gene-chip technology was employed to study the effect of dietary vitamin E (VE) on gene expression in rat testes. Male albino rats were fed with either a diet deficient in VE or a standard diet containing VE. Differential gene expression was monitored at five individual time-points over a period of 14 months with all animals individually pro. led. Low VE intake resulted in the consistent upregulation of 7-dehydrocholesterol reductase and GATA binding protein 4, both involved in testosterone synthesis. Cyclin D3, important in cell cycle progression and Wilms tumor 1, related to cancer development, were also up-regulated in the vitamin E deficient animals. This study demonstrates that low dietary VE intake has long-term effects on gene expression in the testes. Our data provides insights into the possible molecular mechanisms underlying the beneficial effects of vitamin E on the male reproductive organ.

Effects of apoE genotype on macrophage inflammation and heme oxygenase-1 expression

In order to gain a more comprehensive understanding of the aetiology of apolipoprotein E4 genotype-cardiovascular disease (CVD) associations, the impact of the apoE genotype on the macrophage inflammatory response was examined. The murine monocyte-macrophage cell line (RAW 264.7) stably transfected to produce equal amounts of human apoE3 or apoE4 was used. Following LPS stimulation, apoE4-macrophages showed higher and lower concentrations of tumour necrosis factor alpha (pro-inflammatory) and interleukin 10 (anti-inflammatory), respectively, both at mRNA and protein levels. In addition, increased expression of heme oxygenase-1 (a stress-induced anti-inflammatory protein) was observed in the apoE4-cells. Furthermore, in apoE4-macrophages, an enhanced transactivation of the key redox sensitive transcription factor NF-kappa B was shown. Current data indicate that apoE4 macrophages have an altered inflammatory response, which may contribute to the higher CVD risk observed in apoE4 carriers. (c) 2007 Elsevier Inc. All rights reserved.

Inhibition of p38/CREB phosphorylation and COX-2 expression by olive oil polyphenols underlies their anti-proliferative effects

We investigated the anti-proliferative effects of an olive oil polyphenolic extract on human colon adenocarcinoma cells. Analysis indicated that the extract contained hydroxytyrosol, tyrosol and the various secoiridoid derivatives, including oleuropein. This extract exerted a strong inhibitory effect on cancer cell proliferation, which was linked to the induction of a G2/M phase cell cycle block. Following treatment with the extract (50 mu g/ml) the number of cells in the G2/M phase increased to 51.82 +/- 2.69% relative to control cells (15.1 +/- 2.5%). This G2/M block was mediated by the ability of olive oil polyphenols (50 mu g/ml) to exert rapid inhibition of p38 (38.7 +/- 4.7%) and CREB (28.6 +/- 5.5%) phosphorylation which led to a downstream reduction in COX-2 expression (56.9 +/- 9.3%). Our data suggest that olive oil polyphenols may exert chemo preventative effects in the large intestine by interacting with signalling pathways responsible for colorectal cancer development. (c) 2007 Elsevier Inc. All rights reserved.

Identification of hepatic molecular mechanisms of action of alpha-tocopherol using global gene expression profile analysis in rats

The recent discovery that vitamin E (VE) regulates gene activity at the transcriptional level indicates that VE may exert part of its biological effects by mechanisms which may be independent of its well-recognised antioxidant function. The objective of this study was the identification of hepatic vitamin E-sensitive genes and examination of the effects of VE on their corresponding biological endpoints. Two groups of male rats were randomly assigned to either a VE-sufficient diet or to a control diet deficient in VE for 290 days. High-density oligonucleotide microarrays comprising over 7000 genes were used to assess the transcriptional response of the liver. Differential gene expression was monitored over a period of 9 months, at four different time-points, and rats were individually profiled. This experimental strategy identified several VE-sensitive genes, which were chronically altered by dietary VE. VE supplementation down-regulated scavenger receptor CD36, coagulation factor IX and 5-alpha-steroid reductase type 1 mRNA levels while hepatic gamma glutamyl-cysteinyl synthetase was significantly up-regulated. Measurement of the corresponding biological endpoints such as activated partial thromboplastin time, plasma dihydrotestosterone and hepatic glutathione substantiated the gene chip data which indicated that dietary VE plays an important role in a range of metabolic processes within the liver. (C) 2004 Elsevier B.V. All rights reserved.