Surrogate variable analysis using partial least squares (SVA-PLS) in gene expression studies

MOTIVATION: In a typical gene expression profiling study, our prime objective is to identify the genes that are differentially expressed between the samples from two different tissue types. Commonly, standard analysis of variance (ANOVA)/regression is implemented to identify the relative effects of these genes over the two types of samples from their respective arrays of expression levels. But, this technique becomes fundamentally flawed when there are unaccounted sources of variability in these arrays (latent variables attributable to different biological, environmental or other factors relevant in the context). These factors distort the true picture of differential gene expression between the two tissue types and introduce spurious signals of expression heterogeneity. As a result, many genes which are actually differentially expressed are not detected, whereas many others are falsely identified as positives. Moreover, these distortions can be different for different genes. Thus, it is also not possible to get rid of these variations by simple array normalizations. This both-way error can lead to a serious loss in sensitivity and specificity, thereby causing a severe inefficiency in the underlying multiple testing problem. In this work, we attempt to identify the hidden effects of the underlying latent factors in a gene expression profiling study by partial least squares (PLS) and apply ANCOVA technique with the PLS-identified signatures of these hidden effects as covariates, in order to identify the genes that are truly differentially expressed between the two concerned tissue types. RESULTS: We compare the performance of our method SVA-PLS with standard ANOVA and a relatively recent technique of surrogate variable analysis (SVA), on a wide variety of simulation settings (incorporating different effects of the hidden variable, under situations with varying signal intensities and gene groupings). In all settings, our method yields the highest sensitivity while maintaining relatively reasonable values for the specificity, false discovery rate and false non-discovery rate. Application of our method to gene expression profiling for acute megakaryoblastic leukemia shows that our method detects an additional six genes, that are missed by both the standard ANOVA method as well as SVA, but may be relevant to this disease, as can be seen from mining the existing literature.     

Characterization of a human 12/15-lipoxygenase promoter variant associated with atherosclerosis identifies vimentin as a promoter binding protein

Background

Sequence variation in the human 12/15 lipoxygenase (ALOX15) has been associated with atherosclerotic disease. We functionally characterized an ALOX15 promoter polymorphism, rs2255888, previously associated with carotid plaque burden.

Methodology/Principal Findings

We demonstrate specific in vitro and in vivo binding of the cytoskeletal protein, vimentin, to the ALOX15 promoter. We show that the two promoter haplotypes carrying alternate alleles at rs2255888 exhibit significant differences in promoter activity by luciferase reporter assay in two cell lines. Differences in in-vitro vimentin-binding to and formation of DNA secondary structures in the polymorphic promoter sequence are also detected by electrophoretic mobility shift assay and biophysical analysis, respectively. We show regulation of ALOX15 protein by vimentin.

Conclusions/Significance

This study suggests that vimentin binds the ALOX15 promoter and regulates its promoter activity and protein expression. Sequence variation that results in changes in DNA conformation and vimentin binding to the promoter may be relevant to ALOX15 gene regulation.     

Evolution of tryptophan biosynthetic pathway in microbial genomes: a comparative genetic study

Biosynthetic pathway evolution needs to consider the evolution of a group of genes that code for enzymes catalysing the multiple chemical reaction steps leading to the final end product. Tryptophan biosynthetic pathway has five chemical reaction steps that are highly conserved in diverse microbial genomes, though the genes of the pathway enzymes show considerable variations in arrangements, operon structure (gene fusion and splitting) and regulation. We use a combined bioinformatic and statistical analyses approach to address the question if the pathway genes from different microbial genomes, belonging to a wide range of groups, show similar evolutionary relationships within and between them. Our analyses involved detailed study of gene organization (fusion/splitting events), base composition, relative synonymous codon usage pattern of the genes, gene expressivity, amino acid usage, etc. to assess inter- and intra-genic variations, between and within the pathway genes, in diverse group of microorganisms. We describe these genetic and genomic variations in the tryptophan pathway genes in different microorganisms to show the similarities across organisms, and compare the same genes across different organisms to find the possible variability arising possibly due to horizontal gene transfers. Such studies form the basis for moving from single gene evolution to pathway evolutionary studies that are important steps towards understanding the systems biology of intracellular pathways.     

Aromatic-aromatic interactions: analysis of π-π interactions in interleukins and TNF proteins

Aromatic-aromatic hydrogen bonds are important in many areas of chemistry, biology and materials science. In this study we have analyzed the roles played by the π-π interactions in interleukins (ILs) and tumor necrosis factor (TNF) proteins. Majority of π-π interacting residues are conserved in ILs and TNF proteins. The accessible surface area calculations in these proteins reveal that these interactions might be important in stabilizing the inner core regions of these proteins. In addition to π-π interactions, the aromatic residues also form π-networks in ILs and TNF proteins. The results obtained in the present study indicate that π-π interactions and π-π networks play important roles in the structural stability of ILs and TNF proteins.     

Anchoring of palladium(II) in chemically modified mesoporous silica: An efficient heterogeneous catalyst for Suzuki cross-coupling reaction

The synthesis and characterization of a highly efficient and reusable catalyst, Pd(II) immobilized in mesoporous silica MCM-41, are described. Pd(II) Schiff-base moiety has been anchored onto mesoporous silica surface via silicon alkoxide chemistry. The catalyst has been characterized by small-angle X-ray diffraction (SAX), FTIR and electronic spectroscopy as well as elemental analysis. The catalyst is used in Suzuki cross-coupling reaction of various aryl halides, including less reactive chlorobenzene, and phenylboronic acid to give biaryls in excellent yields without any additive or ligand. High selectivity for the bi-aryl products containing both electron-donating and electron-withdrawing substituents, mild reaction conditions and possibility of easy recycle makes the catalyst highly desirable to address the industrial needs and environmental concerns.     

Constitutive activation of chaperone-mediated autophagy in cells with impaired macroautophagy

Three different types of autophagy-macroautophagy, microautophagy, and chaperone-mediated autophagy (CMA)-contribute to degradation of intracellular components in lysosomes in mammalian cells. Although some level of basal macroautophagy and CMA activities has been described in different cell types and tissues, these two pathways are maximally activated under stress conditions. Activation of these two pathways is often sequential, suggesting the existence of some level of cross-talk between both stress-related autophagic pathways. In this work, we analyze the consequences of blockage of macroautophagy on CMA activity. Using mouse embryonic fibroblasts deficient in Atg5, an autophagy-related protein required for autophagosome formation, we have found that blockage of macroautophagy leads to up-regulation of CMA, even under basal conditions. Interestingly, different mechanisms contribute to the observed changes in CMA-related proteins and the consequent activation of CMA during basal and stress conditions in these macroautophagy-deficient cells. This work supports a direct cross-talk between these two forms of autophagy, and it identifies changes in the lysosomal compartment that underlie the basis for the communication between both autophagic pathways.     

RankAggreg, an R package for weighted rank aggregation

Background  

Researchers in the field of bioinformatics often face a challenge of combining several ordered lists in a proper and efficient manner. Rank aggregation techniques offer a general and flexible framework that allows one to objectively perform the necessary aggregation. With the rapid growth of high-throughput genomic and proteomic studies, the potential utility of rank aggregation in the context of meta-analysis becomes even more apparent. One of the major strengths of rank-based aggregation is the ability to combine lists coming from different sources and platforms, for example different microarray chips, which may or may not be directly comparable otherwise.     

Induction of Autophagy by Cystatin C: A Mechanism That Protects Murine Primary Cortical Neurons and Neuronal Cell Lines

Cystatin C (CysC) expression in the brain is elevated in human patients with epilepsy, in animal models of neurodegenerative conditions, and in response to injury, but whether up-regulated CysC expression is a manifestation of neurodegeneration or a cellular repair response is not understood. This study demonstrates that human CysC is neuroprotective in cultures exposed to cytotoxic challenges, including nutritional-deprivation, colchicine, staurosporine, and oxidative stress. While CysC is a cysteine protease inhibitor, cathepsin B inhibition was not required for the neuroprotective action of CysC. Cells responded to CysC by inducing fully functional autophagy via the mTOR pathway, leading to enhanced proteolytic clearance of autophagy substrates by lysosomes. Neuroprotective effects of CysC were prevented by inhibiting autophagy with beclin 1 siRNA or 3-methyladenine. Our findings show that CysC plays a protective role under conditions of neuronal challenge by inducing autophagy via mTOR inhibition and are consistent with CysC being neuroprotective in neurodegenerative diseases. Thus, modulation of CysC expression has therapeutic implications for stroke, Alzheimer''s disease, and other neurodegenerative disorders.     

Chloroplasts and mitochondria have multiple heat tolerant isozymes of SOD and APX in leaf and inflorescence in Chenopodium album

Thermal stability of antioxidant defense enzymes superoxide dismutase (SOD, EC 1.15.1.1) and ascorbate peroxidase (APX, EC 1.11.1.11) was studied in chloroplasts and mitochondria of leaf and inflorescence in heat adaptive weed Chenopodium album. Leaf samples were taken in March (31 °C/14 °C) and young inflorescence (INF) was sampled at flowering in April (40 °C/21 °C). Leaf and INF chloroplast and mitochondrial fractions were subjected to elevated temperatures in vitro (5–100 °C) for 30′. SOD and APX showed activity even after boiling treatment in both chloroplast and mitochondria of leaf and INF. SOD was more heat stable than APX in both chloroplasts and mitochondria in both the tissues. Chloroplast contained more heat stable SOD and APX isozymes than mitochondria in both leaf and INF. To the best of our knowledge this is the first report showing presence of thermostable APX isozymes (100 °C for 30′) in chloroplasts and mitochondria in C. album. Heat stable isozymes of SOD and APX in chloroplasts and mitochondria in leaves and inflorescence may contribute to heat tolerance in C. album.