Differential Gene Expression Profiles May Differentiate Responder and Nonresponder Patients with Rheumatoid Arthritis for Methotrexate (MTX) Monotherapy and MTX plus Tumor Necrosis Factor Inhibitor Combined Therapy

Objective. We aimed to evaluate whether the differential gene expression profiles of patients with rheumatoid arthritis (RA) could distinguish responders from nonresponders to methotrexate (MTX) and, in the case of MTX nonresponders, responsiveness to MTX plus anti-tumor necrosis factor-alpha (anti-TNF) combined therapy. Methods. We evaluated 25 patients with RA taking MTX 15-20 mg/week as a monotherapy (8 responders and 17 nonresponders). All MTX nonresponders received intliximab and were reassessed after 20 weeks to evaluate their anti-TNF responsiveness using the European League Against Rheumatism response criteria. A differential gene expression analysis from peripheral blood mononuclear cells was performed in terms of hierarchical gene clustering, and an evaluation of differentially expressed genes was performed using the significance analysis of microarrays program. Results. Hierarchical gene expression clustering discriminated MTX responders from nonresponders, and MTX plus anti-TNF responders from nonresponders. The evaluation of only highly modulated genes (fold change > 1.3 or < 0.7) yielded 5 induced (4 antiapoptotic and CCL4) and 4 repressed (4 proapoptotic) genes in MTX nonresponders compared to responders. In MTX plus anti-TNF nonresponders, the CCL4, CD83, and BCL2A1 genes were induced in relation to responders. Conclusion. Study of the gene expression profiles of RA peripheral blood cells permitted differentiation of responders from nonresponders to MTX and anti-TNF. Several candidate genes in MTX non-responders (CCL4, HTRA2, PRKCD, BCL2A1, CAV1, TNIP1 CASP8AP2, MXD1, and BTG2) and 3 genes in MTX plus anti-TNF nonresponders (CCL4, CD83, and BCL2A1) were identified for further study. (First Release July 1 2012; J Rheumatol 2012;39:1524-32; doi:10.3899/jrheum.120092)

Potential costs of bacterial infection on storage protein gene expression and reproduction in queenless Apis mellifera worker bees on distinct dietary regimes

Insects are able to combat infection by initiating an efficient immune response that involves synthesizing antimicrobial peptides and a range of other defense molecules. These responses may be costly to the organism, resulting in it exploiting endogenous resources to maintain homeostasis or support defense to the detriment of other physiological needs. We used queenless worker bees on distinct dietary regimes that may alter hemolymph protein storage and ovary activation to investigate the physiological costs of infection with Serratia marcescens. The expression of the genes encoding the storage proteins vitellogenin and hexamerin 70a, the vitellogenin receptor, and vasa (which has a putative role in reproduction), was impaired in the infected bees. This impairment was mainly evident in the bees fed beebread, which caused significantly higher expression of these genes than did royal jelly or syrup, and this was confirmed at the vitellogenin and hexamerin 70a protein levels. Beebread was also the only diet that promoted ovary activation in the queenless bees, but this activation was significantly impaired by the infection. The expression of the genes encoding the storage proteins apolipophorins-I and -III and the lipophorin receptor was not altered by infection regardless the diet provided to the bees. Similarly, the storage of apolipophorin-I in the hemolymph was only slightly impaired by the infection, independently of the supplied diet. Taken together these results indicate that, infection demands a physiological cost from the transcription of specific protein storage-related genes and from the reproductive capacity. (C) 2012 Elsevier Ltd. All rights reserved.

Trypanosoma cruzi Gene Expression in Response to Gamma Radiation

Trypanosoma cruzi is an organism highly resistant to ionizing radiation. Following a dose of 500 Gy of gamma radiation, the fragmented genomic DNA is gradually reconstructed and the pattern of chromosomal bands is restored in less than 48 hours. Cell growth arrests after irradiation but, while DNA is completely fragmented, RNA maintains its integrity. In this work we compared the transcriptional profiles of irradiated and non-irradiated epimastigotes at different time points after irradiation using microarray. In total, 273 genes were differentially expressed; from these, 160 were up-regulated and 113 down-regulated. We found that genes with predicted functions are the most prevalent in the down-regulated gene category. Translation and protein metabolic processes, as well as generation of precursor of metabolites and energy pathways were affected. In contrast, the up-regulated category was mainly composed of obsolete sequences (which included some genes of the kinetoplast DNA), genes coding for hypothetical proteins, and Retrotransposon Hot Spot genes. Finally, the tyrosyl-DNA phosphodiesterase 1, a gene involved in double-strand DNA break repair process, was up-regulated. Our study demonstrated the peculiar response to ionizing radiation, raising questions about how this organism changes its gene expression to manage such a harmful stress.

Differential gene expression and developmental competence in in vitro produced bovine embryos

The embryonic developmental block occurs at the 8-cell stage in cattle and is characterized by a lengthening of the cell cycle and an increased number of embryos that stop development. The maternal-embryonic transition arises at the same stage resulting in the transcription of many genes. Gene expression studies during this stage may contribute to the understanding of the physiological mechanisms involved in the maternal-embryonic transition. Herein we identified genes differentially expressed between embryos with high or low developmental competence to reach the blastocyst stage using differential display PCR. Embryos were analysed according to developmental kinetics: fast cleavage embryos showing 8 cells at 48 h post insemination (hpi) with high potential of development (F8), and embryos with slow cleavage presenting 4 cells at 48 hpi (54) and 8 cells at 90 hpi (S8), both with reduced rates of development to blastocyst. The fluorescence DDPCR method was applied and allowed the recovery of 176 differentially expressed bands with similar proportion between high and low development potential groups (52% to F8 and 48% in S4 and S8 groups). A total of 27 isolated fragments were cloned and sequenced, confirming the expected primer sequences and allowing the identification of 27 gene transcripts. PI3KCA and ITM2B were chosen for relative quantification of mRNA using real-time PCR and showed a kinetic and a time-related pattern of expression respectively. The observed results suggest the existence of two different embryonic genome activation mechanisms: fast-developing embryos activate genes related to embryonic development, and slow-developing embryos activate genes related to cellular survival and/or death.

Changes in spatial and temporal gene expression during incompatible interaction between common bean and anthracnose pathogen

Common bean, one of the most important legumes for human consumption, may have drastic reduction in yield due to anthracnose, a disease caused by the fungus Colletotrichum lindemuthianum. Rapid induction of the plant defense mechanisms is essential to establish an incompatible interaction with this pathogenic fungus. In this study, we evaluated spatial (leaves, epicotyls and hypocotyls) and temporal (24, 48, 72 and 96 hours after inoculation [HAI]) relative expression (RE) of 12 defense-related transcripts selected from previously developed ESTs libraries, during incompatible interaction between the resistant common bean genotype SEL 1308 and the avirulent anthracnose pathogen race 73, using real time quantitative RT-PCR (RT-qPCR) analysis. All selected transcripts, including the ones coding for pathogenesis-related (PR) proteins (PR1a, PR1b, PR2, and PR16a and PR16b) were differentially regulated upon pathogen inoculation. The expression levels of these transcripts were dependent on the tissue and time post inoculation. This study contributes to a better understanding of the kinetics of induced defenses against a fungal pathogen of common bean and may be used as a base line to study defenses against a broad range of pathogens including bacteria as well as non-host resistance. (C) 2012 Elsevier GmbH. All rights reserved.

Gene Expression Profile in Response to Doxorubicin-Rapamycin Combined Treatment of HER-2-Overexpressing Human Mammary Epithelial Cell Lines

HER-2-positive breast cancers frequently sustain elevated AKT/mTOR signaling, which has been associated with resistance to doxorubicin treatment. Here, we investigated whether rapamycin, an mTOR inhibitor, increased the sensitivity to doxorubicin therapy in two HER-2-overexpressing cell lines: C5.2, which was derived from the parental HB4a by transfection with HER-2 and SKBR3, which exhibits HER-2 amplification. The epithelial mammary cell line HB4a was also analyzed. The combined treatment using 20 nmol/L of rapamycin and 30 nmol/L of doxorubicin arrested HB4a and C5.2 cells in S to G(2)-M, whereas SKBR3 cells showed an increase in the G(0)-G(1) phase. Rapamycin increased the sensitivity to doxorubicin in HER-2-overexpressing cells by approximately 2-fold, suggesting that the combination displayed a more effective antiproliferative action. Gene expression profiling showed that these results might reflect alterations in genes involved in canonical pathways related to purine metabolism, oxidative phosphorylation, protein ubiquitination, and mitochondrial dysfunction. A set of 122 genes modulated by the combined treatment and specifically related to HER-2 overexpression was determined by finding genes commonly regulated in both C5.2 and SKBR3 that were not affected in HB4a cells. Network analysis of this particular set showed a smaller subgroup of genes in which coexpression pattern in HB4a cells was disrupted in C5.2 and SKBR3. Altogether, our data showed a subset of genes that might be more robust than individual markers in predicting the response of HER-2-overexpressing breast cancers to doxorubicin and rapamycin combination. Mol Cancer Ther; 11(2); 464-74. (C) 2011 AACR.

Gene expression of peripheral blood lymphocytes may discriminate patients with schizophrenia from controls

To identify a classifier in schizophrenia, blood gene expression profiling was applied to patients with schizophrenia under different treatments and to controls. Expression of six genes discriminated patients with sensitivity of 89.3% and specificity of 90%, supporting the use of peripheral blood as biological material for diagnosis in schizophrenia. (C) 2012 Elsevier Ireland Ltd. All rights reserved.

GATA3 and a dominant regulatory gene expression profile discriminate operational tolerance in human transplantation

Some organ-transplanted patients achieve a state of "operational tolerance" (01) in which graft function is maintained after the complete withdrawal of immunosuppressive drugs. We used a gene panel of regulatory/inflammatory molecules (FOXP3, GATA3, 100, TGFB1, TGFBR1/TBX21, TNF and IFNG) to investigate the gene expression profile in peripheral blood mononuclear cells of renal-transplanted individuals experiencing OT compared to transplanted individuals not displaying OT and healthy individuals (HI). OT subjects showed a predominant regulatory (REG) profile with higher gene expression of GATA3, FOXP3, TGFB1 and TGFB receptor 1 compared to the other groups. This predominant REG gene expression profile displayed stability over time. The significant GATA3 gene and protein expressions in OT individuals suggest that a Th2 deviation may be a relevant pathway to OT. Moreover, the capacity of the REG/INFLAMMA gene panel to discriminate OT by peripheral blood analysis indicates that this state has systemic repercussions. (C) 2011 Elsevier Inc. All rights reserved.

Bayesian model accounting for within-class biological variability in Serial Analysis of Gene Expression (SAGE)

Abstract Background An important challenge for transcript counting methods such as Serial Analysis of Gene Expression (SAGE), "Digital Northern" or Massively Parallel Signature Sequencing (MPSS), is to carry out statistical analyses that account for the within-class variability, i.e., variability due to the intrinsic biological differences among sampled individuals of the same class, and not only variability due to technical sampling error. Results We introduce a Bayesian model that accounts for the within-class variability by means of mixture distribution. We show that the previously available approaches of aggregation in pools ("pseudo-libraries") and the Beta-Binomial model, are particular cases of the mixture model. We illustrate our method with a brain tumor vs. normal comparison using SAGE data from public databases. We show examples of tags regarded as differentially expressed with high significance if the within-class variability is ignored, but clearly not so significant if one accounts for it. Conclusion Using available information about biological replicates, one can transform a list of candidate transcripts showing differential expression to a more reliable one. Our method is freely available, under GPL/GNU copyleft, through a user friendly web-based on-line tool or as R language scripts at supplemental web-site.

Simcluster: clustering enumeration gene expression data on the simplex space

Abstract Background Transcript enumeration methods such as SAGE, MPSS, and sequencing-by-synthesis EST "digital northern", are important high-throughput techniques for digital gene expression measurement. As other counting or voting processes, these measurements constitute compositional data exhibiting properties particular to the simplex space where the summation of the components is constrained. These properties are not present on regular Euclidean spaces, on which hybridization-based microarray data is often modeled. Therefore, pattern recognition methods commonly used for microarray data analysis may be non-informative for the data generated by transcript enumeration techniques since they ignore certain fundamental properties of this space. Results Here we present a software tool, Simcluster, designed to perform clustering analysis for data on the simplex space. We present Simcluster as a stand-alone command-line C package and as a user-friendly on-line tool. Both versions are available at: http://xerad.systemsbiology.net/simcluster. Conclusion Simcluster is designed in accordance with a well-established mathematical framework for compositional data analysis, which provides principled procedures for dealing with the simplex space, and is thus applicable in a number of contexts, including enumeration-based gene expression data.

GEDI: a user-friendly toolbox for analysis of large-scale gene expression data

Abstract Background Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills. Results Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al. Conclusion GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.

Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Abstract Background Spotted cDNA microarrays generally employ co-hybridization of fluorescently-labeled RNA targets to produce gene expression ratios for subsequent analysis. Direct comparison of two RNA samples in the same microarray provides the highest level of accuracy; however, due to the number of combinatorial pair-wise comparisons, the direct method is impractical for studies including large number of individual samples (e.g., tumor classification studies). For such studies, indirect comparisons using a common reference standard have been the preferred method. Here we evaluated the precision and accuracy of reconstructed ratios from three indirect methods relative to ratios obtained from direct hybridizations, herein considered as the gold-standard. Results We performed hybridizations using a fixed amount of Cy3-labeled reference oligonucleotide (RefOligo) against distinct Cy5-labeled targets from prostate, breast and kidney tumor samples. Reconstructed ratios between all tissue pairs were derived from ratios between each tissue sample and RefOligo. Reconstructed ratios were compared to (i) ratios obtained in parallel from direct pair-wise hybridizations of tissue samples, and to (ii) reconstructed ratios derived from hybridization of each tissue against a reference RNA pool (RefPool). To evaluate the effect of the external references, reconstructed ratios were also calculated directly from intensity values of single-channel (One-Color) measurements derived from tissue sample data collected in the RefOligo experiments. We show that the average coefficient of variation of ratios between intra- and inter-slide replicates derived from RefOligo, RefPool and One-Color were similar and 2 to 4-fold higher than ratios obtained in direct hybridizations. Correlation coefficients calculated for all three tissue comparisons were also similar. In addition, the performance of all indirect methods in terms of their robustness to identify genes deemed as differentially expressed based on direct hybridizations, as well as false-positive and false-negative rates, were found to be comparable. Conclusion RefOligo produces ratios as precise and accurate as ratios reconstructed from a RNA pool, thus representing a reliable alternative in reference-based hybridization experiments. In addition, One-Color measurements alone can reconstruct expression ratios without loss in precision or accuracy. We conclude that both methods are adequate options in large-scale projects where the amount of a common reference RNA pool is usually restrictive.

Modeling gene expression regulatory networks with the sparse vector autoregressive model

Abstract Background To understand the molecular mechanisms underlying important biological processes, a detailed description of the gene products networks involved is required. In order to define and understand such molecular networks, some statistical methods are proposed in the literature to estimate gene regulatory networks from time-series microarray data. However, several problems still need to be overcome. Firstly, information flow need to be inferred, in addition to the correlation between genes. Secondly, we usually try to identify large networks from a large number of genes (parameters) originating from a smaller number of microarray experiments (samples). Due to this situation, which is rather frequent in Bioinformatics, it is difficult to perform statistical tests using methods that model large gene-gene networks. In addition, most of the models are based on dimension reduction using clustering techniques, therefore, the resulting network is not a gene-gene network but a module-module network. Here, we present the Sparse Vector Autoregressive model as a solution to these problems. Results We have applied the Sparse Vector Autoregressive model to estimate gene regulatory networks based on gene expression profiles obtained from time-series microarray experiments. Through extensive simulations, by applying the SVAR method to artificial regulatory networks, we show that SVAR can infer true positive edges even under conditions in which the number of samples is smaller than the number of genes. Moreover, it is possible to control for false positives, a significant advantage when compared to other methods described in the literature, which are based on ranks or score functions. By applying SVAR to actual HeLa cell cycle gene expression data, we were able to identify well known transcription factor targets. Conclusion The proposed SVAR method is able to model gene regulatory networks in frequent situations in which the number of samples is lower than the number of genes, making it possible to naturally infer partial Granger causalities without any a priori information. In addition, we present a statistical test to control the false discovery rate, which was not previously possible using other gene regulatory network models.

Multivariate gene expression analysis reveals functional connectivity changes between normal/tumoral prostates

Abstract Background Prostate cancer is a leading cause of death in the male population, therefore, a comprehensive study about the genes and the molecular networks involved in the tumoral prostate process becomes necessary. In order to understand the biological process behind potential biomarkers, we have analyzed a set of 57 cDNA microarrays containing ~25,000 genes. Results Principal Component Analysis (PCA) combined with the Maximum-entropy Linear Discriminant Analysis (MLDA) were applied in order to identify genes with the most discriminative information between normal and tumoral prostatic tissues. Data analysis was carried out using three different approaches, namely: (i) differences in gene expression levels between normal and tumoral conditions from an univariate point of view; (ii) in a multivariate fashion using MLDA; and (iii) with a dependence network approach. Our results show that malignant transformation in the prostatic tissue is more related to functional connectivity changes in their dependence networks than to differential gene expression. The MYLK, KLK2, KLK3, HAN11, LTF, CSRP1 and TGM4 genes presented significant changes in their functional connectivity between normal and tumoral conditions and were also classified as the top seven most informative genes for the prostate cancer genesis process by our discriminant analysis. Moreover, among the identified genes we found classically known biomarkers and genes which are closely related to tumoral prostate, such as KLK3 and KLK2 and several other potential ones. Conclusion We have demonstrated that changes in functional connectivity may be implicit in the biological process which renders some genes more informative to discriminate between normal and tumoral conditions. Using the proposed method, namely, MLDA, in order to analyze the multivariate characteristic of genes, it was possible to capture the changes in dependence networks which are related to cell transformation.

Identification of protein-coding and non-coding RNA expression profiles in CD34+ and in stromal cells in refractory anemia with ringed sideroblasts

Abstract Background Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods In this study, gene expression profiles of CD34+ cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results In CD34+ cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value ≤ 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value ≤ 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34+ cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.