Differential expression analysis with global network adjustment

<p>Background: Large-scale chromosomal deletions or other non-specific perturbations of the transcriptome can alter the expression of hundreds or thousands of genes, and it is of biological interest to understand which genes are most profoundly affected. We present a method for predicting a gene’s expression as a function of other genes thereby accounting for the effect of transcriptional regulation that confounds the identification of genes differentially expressed relative to a regulatory network. The challenge in constructing such models is that the number of possible regulator transcripts within a global network is on the order of thousands, and the number of biological samples is typically on the order of 10. Nevertheless, there are large gene expression databases that can be used to construct networks that could be helpful in modeling transcriptional regulation in smaller experiments.</p> <p>Results: We demonstrate a type of penalized regression model that can be estimated from large gene expression databases, and then applied to smaller experiments. The ridge parameter is selected by minimizing the cross-validation error of the predictions in the independent out-sample. This tends to increase the model stability and leads to a much greater degree of parameter shrinkage, but the resulting biased estimation is mitigated by a second round of regression. Nevertheless, the proposed computationally efficient “over-shrinkage” method outperforms previously used LASSO-based techniques. In two independent datasets, we find that the median proportion of explained variability in expression is approximately 25%, and this results in a substantial increase in the signal-to-noise ratio allowing more powerful inferences on differential gene expression leading to biologically intuitive findings. We also show that a large proportion of gene dependencies are conditional on the biological state, which would be impossible with standard differential expression methods.</p> <p>Conclusions: By adjusting for the effects of the global network on individual genes, both the sensitivity and reliability of differential expression measures are greatly improved.</p&gt

EzArray: A web-based highly automated Affymetrix expression array data management and analysis system

<p>Abstract</p> <p>Background</p> <p>Though microarray experiments are very popular in life science research, managing and analyzing microarray data are still challenging tasks for many biologists. Most microarray programs require users to have sophisticated knowledge of mathematics, statistics and computer skills for usage. With accumulating microarray data deposited in public databases, easy-to-use programs to re-analyze previously published microarray data are in high demand.</p> <p>Results</p> <p>EzArray is a web-based Affymetrix expression array data management and analysis system for researchers who need to organize microarray data efficiently and get data analyzed instantly. EzArray organizes microarray data into projects that can be analyzed online with predefined or custom procedures. EzArray performs data preprocessing and detection of differentially expressed genes with statistical methods. All analysis procedures are optimized and highly automated so that even novice users with limited pre-knowledge of microarray data analysis can complete initial analysis quickly. Since all input files, analysis parameters, and executed scripts can be downloaded, EzArray provides maximum reproducibility for each analysis. In addition, EzArray integrates with Gene Expression Omnibus (GEO) and allows instantaneous re-analysis of published array data.</p> <p>Conclusion</p> <p>EzArray is a novel Affymetrix expression array data analysis and sharing system. EzArray provides easy-to-use tools for re-analyzing published microarray data and will help both novice and experienced users perform initial analysis of their microarray data from the location of data storage. We believe EzArray will be a useful system for facilities with microarray services and laboratories with multiple members involved in microarray data analysis. EzArray is freely available from <url>http://www.ezarray.com/</url>.</p

GENE-Counter: A Computational Pipeline for the Analysis of RNA-Seq Data for Gene Expression Differences

GENE-counter is a complete Perl-based computational pipeline for analyzing RNA-Sequencing (RNA-Seq) data for differential gene expression. In addition to its use in studying transcriptomes of eukaryotic model organisms, GENE-counter is applicable for prokaryotes and non-model organisms without an available genome reference sequence. For alignments, GENE-counter is configured for CASHX, Bowtie, and BWA, but an end user can use any Sequence Alignment/Map (SAM)-compliant program of preference. To analyze data for differential gene expression, GENE-counter can be run with any one of three statistics packages that are based on variations of the negative binomial distribution. The default method is a new and simple statistical test we developed based on an over-parameterized version of the negative binomial distribution. GENE-counter also includes three different methods for assessing differentially expressed features for enriched gene ontology (GO) terms. Results are transparent and data are systematically stored in a MySQL relational database to facilitate additional analyses as well as quality assessment. We used next generation sequencing to generate a small-scale RNA-Seq dataset derived from the heavily studied defense response of Arabidopsis thaliana and used GENE-counter to process the data. Collectively, the support from analysis of microarrays as well as the observed and substantial overlap in results from each of the three statistics packages demonstrates that GENE-counter is well suited for handling the unique characteristics of small sample sizes and high variability in gene counts

Mitochondrial Physiology and Gene Expression Analyses Reveal Metabolic and Translational Dysregulation in Oocyte-Induced Somatic Nuclear Reprogramming

While reprogramming a foreign nucleus after somatic cell nuclear transfer (SCNT), the enucleated oocyte (ooplasm) must signal that biomass and cellular requirements changed compared to the nucleus donor cell. Using cells expressing nuclear-encoded but mitochondria-targeted EGFP, a strategy was developed to directly distinguish maternal and embryonic products, testing ooplasm demands on transcriptional and post-transcriptional activity during reprogramming. Specifically, we compared transcript and protein levels for EGFP and other products in pre-implantation SCNT embryos, side-by-side to fertilized controls (embryos produced from the same oocyte pool, by intracytoplasmic injection of sperm containing the EGFP transgene). We observed that while EGFP transcript abundance is not different, protein levels are significantly lower in SCNT compared to fertilized blastocysts. This was not observed for Gapdh and Actb, whose protein reflected mRNA. This transcript-protein relationship indicates that the somatic nucleus can keep up with ooplasm transcript demands, whilst transcription and translation mismatch occurs after SCNT for certain mRNAs. We further detected metabolic disturbances after SCNT, suggesting a place among forces regulating post-transcriptional changes during reprogramming. Our observations ascribe oocyte-induced reprogramming with previously unsuspected regulatory dimensions, in that presence of functional proteins may no longer be inferred from mRNA, but rather depend on post-transcriptional regulation possibly modulated through metabolism

A origem das parcerias público-privada na governança global da educação

Durante a última década, a globalização da governança educacional por meio de parcerias público-privadas (PPP) tem gerado considerável debate quanto ao seu significado, propósito, status e resultados. Este debate é particularmente aquecido no setor da educação por causa da ampla aceitação da educação como atividade complexa, social e política que deve permanecer, em grande parte, se não totalmente, no setor público, servindo a interesses públicos. O artigo analisa a rápida expansão das parcerias público-privadas em educação (PPPE) articulada à introdução de regras de mercado no setor. Neste estudo nos concentramos sobre o papel de uma rede de desenvolvimento global, fundamental na globalização de um tipo particular de PPPE, indicando que a ideia de PPP encaixa-se em um projeto mais amplo de reconstituição da educação pública no âmbito do setor de serviços, a ser governada como parte da construção de uma sociedade de mercado.Over the past decade, the globalization and governing of education through Public Private Partnerships (PPPs) have generated considerable debate as to their meaning, purpose, status and outcomes. This debate is particularly heated in the education sector because of the widely-held view that education is a complex social and political activity that should remain largely, if not wholly, in the public sector serving public interests. The article analyses the rapid expansion of Education Public Private Partnerships (EPPPs) and the associated introduction of market rules into the education sector. We focus on the role of a key global development network in globalizing a particular kind of ePPPs, and show that the EPPP idea fi ts into a wider project of reconstituting public education as an education services industry to be governed as part of the construction of a market society

Saltatory remodeling of Hox chromatin in response to rostrocaudal patterning signals

Hox genes controlling motor neuron subtype identity are expressed in rostrocaudal patterns that are spatially and temporally collinear with their chromosomal organization. Here we demonstrate that Hox chromatin is subdivided into discrete domains that are controlled by rostrocaudal patterning signals that trigger rapid, domain-wide clearance of repressive histone H3 Lys27 trimethylation (H3K27me3) polycomb modifications. Treatment of differentiating mouse neural progenitors with retinoic acid leads to activation and binding of retinoic acid receptors (RARs) to the Hox1–Hox5 chromatin domains, which is followed by a rapid domain-wide removal of H3K27me3 and acquisition of cervical spinal identity. Wnt and fibroblast growth factor (FGF) signals induce expression of the Cdx2 transcription factor that binds and clears H3K27me3 from the Hox1–Hox9 chromatin domains, leading to specification of brachial or thoracic spinal identity. We propose that rapid clearance of repressive modifications in response to transient patterning signals encodes global rostrocaudal neural identity and that maintenance of these chromatin domains ensures the transmission of positional identity to postmitotic motor neurons later in development.Leona M. and Harry B. Helmsley Charitable TrustNational Institutes of Health (U.S.) (Grant P01 NS055923)Smith Family Foundatio

A Molecular Signature of Proteinuria in Glomerulonephritis

Proteinuria is the most important predictor of outcome in glomerulonephritis and experimental data suggest that the tubular cell response to proteinuria is an important determinant of progressive fibrosis in the kidney. However, it is unclear whether proteinuria is a marker of disease severity or has a direct effect on tubular cells in the kidneys of patients with glomerulonephritis. Accordingly we studied an in vitro model of proteinuria, and identified 231 “albumin-regulated genes” differentially expressed by primary human kidney tubular epithelial cells exposed to albumin. We translated these findings to human disease by studying mRNA levels of these genes in the tubulo-interstitial compartment of kidney biopsies from patients with IgA nephropathy using microarrays. Biopsies from patients with IgAN (n = 25) could be distinguished from those of control subjects (n = 6) based solely upon the expression of these 231 “albumin-regulated genes.” The expression of an 11-transcript subset related to the degree of proteinuria, and this 11-mRNA subset was also sufficient to distinguish biopsies of subjects with IgAN from control biopsies. We tested if these findings could be extrapolated to other proteinuric diseases beyond IgAN and found that all forms of primary glomerulonephritis (n = 33) can be distinguished from controls (n = 21) based solely on the expression levels of these 11 genes derived from our in vitro proteinuria model. Pathway analysis suggests common regulatory elements shared by these 11 transcripts. In conclusion, we have identified an albumin-regulated 11-gene signature shared between all forms of primary glomerulonephritis. Our findings support the hypothesis that albuminuria may directly promote injury in the tubulo-interstitial compartment of the kidney in patients with glomerulonephritis

Beta-carotene affects gene expression in lungs of male and female Bcmo1−/− mice in opposite directions

Molecular mechanisms triggered by high dietary beta-carotene (BC) intake in lung are largely unknown. We performed microarray gene expression analysis on lung tissue of BC supplemented beta-carotene 15,15′-monooxygenase 1 knockout (Bcmo1−/−) mice, which are—like humans—able to accumulate BC. Our main observation was that the genes were regulated in an opposite direction in male and female Bcmo1−/− mice by BC. The steroid biosynthetic pathway was overrepresented in BC-supplemented male Bcmo1−/− mice. Testosterone levels were higher after BC supplementation only in Bcmo1−/− mice, which had, unlike wild-type (Bcmo1+/+) mice, large variations. We hypothesize that BC possibly affects hormone synthesis or metabolism. Since sex hormones influence lung cancer risk, these data might contribute to an explanation for the previously found increased lung cancer risk after BC supplementation (ATBC and CARET studies). Moreover, effects of BC may depend on the presence of frequent human BCMO1 polymorphisms, since these effects were not found in wild-type mice

Molecular Characterization of Spontaneous Mesenchymal Stem Cell Transformation

Background. We previously reported the in vitro spontaneous transformation of human mesenchymal stem cells (MSC) generating a population with tumorigenic potential, that we termed transformed mesenchymal cells (TMC). Methodology/Principal Findings. Here we have characterized the molecular changes associated with TMC generation. Using microarrays techniques we identified a set of altered pathways and a greater number of downregulated than upregulated genes during MSC transformation, in part due to the expression of many untranslated RNAs in MSC. Microarray results were validated by qRT-PCR and protein detection. Conclusions/Significance. In our model, the transformation process takes place through two sequential steps; first MSC bypass senescence by upregulating c-myc and repressing p16 levels. The cells then bypass cell crisis with acquisition of telomerase activity, Ink4a/Arf locus deletion and Rb hyperphosphorylation. Other transformation-associated changes include modulation of mitochondrial metabolism, DNA damage-repair proteins and cell cycle regulators. In this work we have characterized the molecular mechanisms implicated in TMC generation and we propose a two-stage model by which a human MSC becomes a tumor cell