Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Correlation of meta 1 expression with culture stage, cell morphology and infectivity in Leishmania (Leishmania) amazonensis promastigotes

The parasitic protozoan Leishmania (Leishmania) amazonensis alternates between mammalian and insect hosts. In the insect host, the parasites proliferate as procyclic promastigotes andthen differentiate into metacyclic infective forms. The meta 1 gene is preferentially expressed during metacyclogenesis. Meta 1 expression profile determination along parasite growth curves revealed that the meta 1 mRNA level peaked at the early stationary phase then decreased to an intermediate level. No correlation was observed between meta 1 expression and infectivity. Conversely, infectivity correlated with the increase of apoptotic cells in the late stationary phase.

Chemokines and chemokine receptors expression in the lesions of patients with American cutaneous leishmaniasis

American cutaneous leishmaniasis (ACL) presents distinct active clinical forms with different grades of severity, known as localised (LCL), intermediate (ICL) and diffuse (DCL) cutaneous leishmaniasis. LCL and DCL are associated with a polarised T-helper (Th)1 and Th2 immune response, respectively, whereas ICL, or chronic cutaneous leishmaniasis, is associated with an exacerbated immune response and a mixed cytokine expression profile. Chemokines and chemokine receptors are involved in cellular migration and are critical in the inflammatory response. Therefore, we evaluated the expression of the chemokines CXCL10, CCL4, CCL8, CCL11 and CXCL8 and the chemokine receptors CCR3, CXCR3, CCR5 and CCR7 in the lesions of patients with different clinical forms of ACL using immunohistochemistry. LCL patients exhibited a high density of CXCL10+, CCL4+ and CCL8+ cells, indicating an important role for these chemokines in the local Th1 immune response and the migration of CXCR3+ cells. LCL patients showed a higher density of CCR7+ cells than ICL or DCL patients, suggesting major dendritic cell (DC) migration to lymph nodes. Furthermore, DCL was associated with low expression levels of Th1-associated chemokines and CCL11+ epidermal DCs, which contribute to the recruitment of CCR3+ cells. Our findings also suggest an important role for epidermal cells in the induction of skin immune responses through the production of chemokines, such as CXCL10, by keratinocytes.

Conservation and developmental expression of ubiquitin isopeptidases in Schistosoma mansoni

Several genes related to the ubiquitin (Ub)-proteasome pathway, including those coding for proteasome subunits and conjugation enzymes, are differentially expressed during the Schistosoma mansoni life cycle. Although deubiquitinating enzymes have been reported to be negative regulators of protein ubiquitination and shown to play an important role in Ub-dependent processes, little is known about their role in S. mansoni . In this study, we analysed the Ub carboxyl-terminal hydrolase (UCHs) proteins found in the database of the parasite’s genome. An in silicoana- lysis (GeneDB and MEROPS) identified three different UCH family members in the genome, Sm UCH-L3, Sm UCH-L5 and SmBAP-1 and a phylogenetic analysis confirmed the evolutionary conservation of the proteins. We performed quantitative reverse transcription-polymerase chain reaction and observed a differential expression profile for all of the investigated transcripts between the cercariae and adult worm stages. These results were corroborated by low rates of Z-Arg-Leu-Arg-Gly-Gly-AMC hydrolysis in a crude extract obtained from cercariae in parallel with high Ub conjugate levels in the same extracts. We suggest that the accumulation of ubiquitinated proteins in the cercaria and early schistosomulum stages is related to a decrease in 26S proteasome activity. Taken together, our data suggest that UCH family members contribute to regulating the activity of the Ub-proteasome system during the life cycle of this parasite.

The expression level of BAALC-associated microRNA miR-3151 is an independent prognostic factor in younger patients with cytogenetic intermediate-risk acute myeloid leukemia.

Acute myeloid leukemia (AML) is a heterogeneous disease whose prognosis is mainly related to the biological risk conferred by cytogenetics and molecular profiling. In elderly patients (60 years) with normal karyotype AML miR-3151 have been identified as a prognostic factor. However, miR-3151 prognostic value has not been examined in younger AML patients. In the present work, we have studied miR-3151 alone and in combination with BAALC, its host gene, in a cohort of 181 younger intermediate-risk AML (IR-AML) patients. Patients with higher expression of miR-3151 had shorter overall survival (P=0.0025), shorter leukemia-free survival (P=0.026) and higher cumulative incidence of relapse (P=0.082). Moreover, in the multivariate analysis miR-3151 emerged as independent prognostic marker in both the overall series and within the unfavorable molecular prognostic category. Interestingly, the combined determination of both miR-3151 and BAALC improved this prognostic stratification, with patients with low levels of both parameters showing a better outcome compared with those patients harboring increased levels of one or both markers (P=0.003). In addition, we studied the microRNA expression profile associated with miR-3151 identifying a six-microRNA signature. In conclusion, the analysis of miR-3151 and BAALC expression may well contribute to an improved prognostic stratification of younger patients with IR-AML.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

Hypothalamus transcriptome profile suggests an anorexia-cachexia syndrome in the anx/anx mouse model

The anx/anx mouse displays poor appetite and lean appearance and is considered a good model for the study of anorexia nervosa. To identify new genes involved in feeding behavior and body weight regulation we performed an expression profiling in the hypothalamus of the anx/anx mice. Using commercial microarrays we detected 156 differentially expressed genes and validated 92 of those using TaqMan low-density arrays. The expression of a set of 87 candidate genes selected based on literature evidences was also quantified by TaqMan low-density arrays. Our results showed enrichment in deregulated genes involved in cell death, cell morphology and cancer as well as an alteration of several signaling circuits involved in energy balance including neuropeptide Y and melanocortin signaling. The expression profile along with the phenotype led us to conclude that anx/anx mice resemble the anorexia-cachexia syndrome typically observed in cancer, infection with human immunodeficiency virus or chronic diseases, rather than starvation, and that anx/anx mice could be considered a good model for the treatment and investigation of this condition.

Repeated exposure to Lutzomyia intermedia sand fly saliva induces local expression of interferon-inducible genes both at the site of injection in mice and in human blood.

During a blood meal, Lutzomyia intermedia sand flies transmit Leishmania braziliensis, a parasite causing tegumentary leishmaniasis. In experimental leishmaniasis, pre-exposure to saliva of most blood-feeding sand flies results in parasite establishment in absence of any skin damages in mice challenged with dermotropic Leishmania species together with saliva. In contrast, pre-immunization with Lu. intermedia salivary gland sonicate (SGS) results in enhanced skin inflammatory exacerbation upon co-inoculation of Lu. intermedia SGS and L. braziliensis. These data highlight potential unique features of both L. braziliensis and Lu. intermedia. In this study, we investigated the genes modulated by Lu. intermedia SGS immunization to understand their potential impact on the subsequent cutaneous immune response following inoculation of both SGS and L. braziliensis. The cellular recruitment and global gene expression profile was analyzed in mice repeatedly inoculated or not with Lu. intermedia. Microarray gene analysis revealed the upregulation of a distinct set of IFN-inducible genes, an immune signature not seen to the same extent in control animals. Of note this INF-inducible gene set was not induced in SGS pre-immunized mice subsequently co-inoculated with SGS and L. braziliensis. These data suggest the parasite prevented the upregulation of this Lu. intermedia saliva-related immune signature. The presence of these IFN-inducible genes was further analyzed in peripheral blood mononuclear cells (PBMCs) sampled from uninfected human individuals living in a L. braziliensis-endemic region of Brazil thus regularly exposed to Lu. intermedia bites. PBMCs were cultured in presence or absence of Lu. intermedia SGS. Using qRT-PCR we established that the IFN-inducible genes induced in the skin of SGS pre-immunized mice, were also upregulated by SGS in PBMCs from human individuals regularly exposed to Lu. intermedia bites, but not in PBMCs of control subjects. These data demonstrate that repeated exposure to Lu. intermedia SGS induces the expression of potentially host-protective IFN-inducible genes.

Regional variations in ABC transporter expression along the mouse intestinal tract.

The ATP-binding cassette (ABC) family of proteins comprise a group of membrane transporters involved in the transport of a wide variety of compounds, such as xenobiotics, vitamins, lipids, amino acids, and carbohydrates. Determining their regional expression patterns along the intestinal tract will further characterize their transport functions in the gut. The mRNA expression levels of murine ABC transporters in the duodenum, jejunum, ileum, and colon were examined using the Affymetrix MuU74v2 GeneChip set. Eight ABC transporters (Abcb2, Abcb3, Abcb9, Abcc3, Abcc6, Abcd1, Abcg5, and Abcg8) displayed significant differential gene expression along the intestinal tract, as determined by two statistical models (a global error assessment model and a classic ANOVA, both with a P < 0.01). Concordance with semiquantitative real-time PCR was high. Analyzing the promoters of the differentially expressed ABC transporters did not identify common transcriptional motifs between family members or with other genes; however, the expression profile for Abcb9 was highly correlated with fibulin-1, and both genes share a common complex promoter model involving the NFkappaB, zinc binding protein factor (ZBPF), GC-box factors SP1/GC (SP1F), and early growth response factor (EGRF) transcription binding motifs. The cellular location of another of the differentially expressed ABC transporters, Abcc3, was examined by immunohistochemistry. Staining revealed that the protein is consistently expressed in the basolateral compartment of enterocytes along the anterior-posterior axis of the intestine. Furthermore, the intensity of the staining pattern is concordant with the expression profile. This agrees with previous findings in which the mRNA, protein, and transport function of Abcc3 were increased in the rat distal intestine. These data reveal regional differences in gene expression profiles along the intestinal tract and demonstrate that a complete understanding of intestinal ABC transporter function can only be achieved by examining the physiologically distinct regions of the gut.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

Uncoupling phosphate deficiency from its major effects on growth and transcriptome via PHO1 expression in Arabidopsis.

Inorganic phosphate (Pi) is one of the most limiting nutrients for plant growth in both natural and agricultural contexts. Pi-deficiency leads to a strong decrease in shoot growth, and triggers extensive changes at the developmental, biochemical and gene expression levels that are presumably aimed at improving the acquisition of this nutrient and sustaining growth. The Arabidopsis thaliana PHO1 gene has previously been shown to participate in the transport of Pi from roots to shoots, and the null pho1 mutant has all the hallmarks associated with shoot Pi deficiency. We show here that A. thaliana plants with a reduced expression of PHO1 in roots have shoot growth similar to Pi-sufficient plants, despite leaves being strongly Pi deficient. Furthermore, the gene expression profile normally triggered by Pi deficiency is suppressed in plants with low PHO1 expression. At comparable levels of shoot Pi supply, the wild type reduces shoot growth but maintains adequate shoot vacuolar Pi content, whereas the PHO1 underexpressor maintains maximal growth with strongly depleted Pi reserves. Expression of the Oryza sativa (rice) PHO1 ortholog in the pho1 null mutant also leads to plants that maintain normal growth and suppression of the Pi-deficiency response, despite the low shoot Pi. These data show that it is possible to unlink low shoot Pi content with the responses normally associated with Pi deficiency through the modulation of PHO1 expression or activity. These data also show that reduced shoot growth is not a direct consequence of Pi deficiency, but is more likely to be a result of extensive gene expression reprogramming triggered by Pi deficiency.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

A microRNA profile of human CD8.

BackgroundRecently, regulatory T (Treg) cells have gained interest in the fields of immunopathology, transplantation and oncoimmunology. Here, we investigated the microRNA expression profile of human natural CD8+CD25+ Treg cells and the impact of microRNAs on molecules associated with immune regulation.MethodsWe purified human natural CD8+ Treg cells and assessed the expression of FOXP3 and CTLA-4 by flow cytometry. We have also tested the ex vivo suppressive capacity of these cells in mixed leukocyte reactions. Using TaqMan low-density arrays and microRNA qPCR for validation, we could identify a microRNA `signature¿ for CD8+CD25+FOXP3+CTLA-4+ natural Treg cells. We used the `TargetScan¿ and `miRBase¿ bioinformatics programs to identify potential target sites for these microRNAs in the 3¿-UTR of important Treg cell-associated genes.ResultsThe human CD8+CD25+ natural Treg cell microRNA signature includes 10 differentially expressed microRNAs. We demonstrated an impact of this signature on Treg cell biology by showing specific regulation of FOXP3, CTLA-4 and GARP gene expression by microRNA using site-directed mutagenesis and a dual-luciferase reporter assay. Furthermore, we used microRNA transduction experiments to demonstrate that these microRNAs impacted their target genes in human primary Treg cells ex vivo.ConclusionsWe are examining the biological relevance of this `signature¿ by studying its impact on other important Treg cell-associated genes. These efforts could result in a better understanding of the regulation of Treg cell function and might reveal new targets for immunotherapy in immune disorders and cancer.

Gene expression profiling in human pathogenic dermatophytes in vitro and in vivo

Objectives: Dermatophytes are highly specialized fungi which are the most common agents of superficial mycoses in humans and animals. The particular ability of these microorganisms to invade and multiply within keratinized host structures is presumably linked to their secreted keratinolytic activity, which is therefore a major putative virulence attribute of these fungi. The overall adaptation and transcriptional response of dermatophytes during protein degradation and/or infection is largely unknown. Methods: A Trichophyton rubrum cDNA microarray was developed and used for the transcriptional analysis of T. rubrum and Arthroderma benhamiae cells during growth on protein substrates. Moreover, the gene expression profile in A. benhamiae cells was monitored during infection of guinea pigs. Results: T. rubrum and A. benhamiae cells activate a large set of genes encoding secreted endo- and exoproteases during growth on soy and keratin. In addition, other specifically induced factors with potential implication in protein utilization were identified, e.g. multiple transporters, metabolic enzymes, transcription factors and hypothetical proteins with unknown function. Notably however, the protease gene expression profile in the fungal cells during infection was significantly different from the pattern elicited during in vitro growth on keratin. Conclusions: Our results suggest specific functions of individual proteases during infection, which may not be restricted to the degradation of keratin. This first, broad in vivo transcriptional profiling approach in dermatophytes gives new molecular insights into pathogenicity associated adaptation mechanisms that make these microorganisms the most successful causitive agents of superficial mycoses.