@article {710, title = {Transcriptomic Analysis of a Diabetic Skin-Humanized Mouse Model Dissects Molecular Pathways Underlying the Delayed Wound Healing Response.}, journal = {Genes (Basel)}, volume = {12}, year = {2020}, month = {2020 12 31}, abstract = {

Defective healing leading to cutaneous ulcer formation is one of the most feared complications of diabetes due to its consequences on patients{\textquoteright} quality of life and on the healthcare system. A more in-depth analysis of the underlying molecular pathophysiology is required to develop effective healing-promoting therapies for those patients. Major architectural and functional differences with human epidermis limit extrapolation of results coming from rodents and other small mammal-healing models. Therefore, the search for reliable humanized models has become mandatory. Previously, we developed a diabetes-induced delayed humanized wound healing model that faithfully recapitulated the major histological features of such skin repair-deficient condition. Herein, we present the results of a transcriptomic and functional enrichment analysis followed by a mechanistic analysis performed in such humanized wound healing model. The deregulation of genes implicated in functions such as angiogenesis, apoptosis, and inflammatory signaling processes were evidenced, confirming published data in diabetic patients that in fact might also underlie some of the histological features previously reported in the delayed skin-humanized healing model. Altogether, these molecular findings support the utility of such preclinical model as a valuable tool to gain insight into the molecular basis of the delayed diabetic healing with potential impact in the translational medicine field.

}, keywords = {Animals, Diabetes Mellitus, Experimental, Gene Expression Profiling, Gene Expression Regulation, Gene ontology, Humans, Metabolic Networks and Pathways, Mice, Mice, Nude, Microarray Analysis, Molecular Sequence Annotation, Principal Component Analysis, Signal Transduction, Skin, Skin Transplantation, Skin Ulcer, Streptozocin, Tissue Engineering, Transcriptome, Transplantation, Heterologous, Wound Healing}, issn = {2073-4425}, doi = {10.3390/genes12010047}, author = {Le{\'o}n, Carlos and Garcia-Garcia, Francisco and Llames, Sara and Garc{\'\i}a-P{\'e}rez, Eva and Carretero, Marta and Arriba, Mar{\'\i}a Del Carmen and Dopazo, Joaquin and Del Rio, Marcela and Escamez, Maria Jos{\'e} and Mart{\'\i}nez-Santamar{\'\i}a, Luc{\'\i}a} } @article {457, title = {Global Transcriptome Analysis of Primary Cerebrocortical Cells: Identification of Genes Regulated by Triiodothyronine in Specific Cell Types.}, journal = {Cereb Cortex}, volume = {27}, year = {2017}, month = {2017 01 01}, pages = {706-717}, abstract = {

Thyroid hormones, thyroxine, and triiodothyronine (T3) are crucial for cerebral cortex development acting through regulation of gene expression. To define the transcriptional program under T3 regulation, we have performed RNA-Seq of T3-treated and untreated primary mouse cerebrocortical cells. The expression of 1145 genes or 7.7\% of expressed genes was changed upon T3 addition, of which 371 responded to T3 in the presence of cycloheximide indicating direct transcriptional regulation. The results were compared with available transcriptomic datasets of defined cellular types. In this way, we could identify targets of T3 within genes enriched in astrocytes and neurons, in specific layers including the subplate, and in specific neurons such as prepronociceptin, cholecystokinin, or cortistatin neurons. The subplate and the prepronociceptin neurons appear as potentially major targets of T3 action. T3 upregulates mostly genes related to cell membrane events, such as G-protein signaling, neurotransmission, and ion transport and downregulates genes involved in nuclear events associated with the M phase of cell cycle, such as chromosome organization and segregation. Remarkably, the transcriptomic changes induced by T3 sustain the transition from fetal to adult patterns of gene expression. The results allow defining in molecular terms the elusive role of thyroid hormones on neocortical development.

}, keywords = {Animals, Astrocytes, Cells, Cultured, Cerebral Cortex, Fluorescent Antibody Technique, Gene Expression Profiling, Mice, 129 Strain, Mice, Inbred BALB C, Mice, Inbred C57BL, Neurons, Piperazines, Transcriptome, Triiodothyronine}, issn = {1460-2199}, doi = {10.1093/cercor/bhv273}, author = {Gil-Iba{\~n}ez, Pilar and Garcia-Garcia, Francisco and Dopazo, Joaquin and Bernal, Juan and Morte, Beatriz} } @article {439, title = {Integrated gene set analysis for microRNA studies.}, journal = {Bioinformatics}, volume = {32}, year = {2016}, month = {2016 09 15}, pages = {2809-16}, abstract = {

MOTIVATION: Functional interpretation of miRNA expression data is currently done in a three step procedure: select differentially expressed miRNAs, find their target genes, and carry out gene set overrepresentation analysis Nevertheless, major limitations of this approach have already been described at the gene level, while some newer arise in the miRNA scenario.Here, we propose an enhanced methodology that builds on the well-established gene set analysis paradigm. Evidence for differential expression at the miRNA level is transferred to a gene differential inhibition score which is easily interpretable in terms of gene sets or pathways. Such transferred indexes account for the additive effect of several miRNAs targeting the same gene, and also incorporate cancellation effects between cases and controls. Together, these two desirable characteristics allow for more accurate modeling of regulatory processes.

RESULTS: We analyze high-throughput sequencing data from 20 different cancer types and provide exhaustive reports of gene and Gene Ontology-term deregulation by miRNA action.

AVAILABILITY AND IMPLEMENTATION: The proposed methodology was implemented in the Bioconductor library mdgsa http://bioconductor.org/packages/mdgsa For the purpose of reproducibility all of the scripts are available at https://github.com/dmontaner-papers/gsa4mirna

CONTACT: : david.montaner@gmail.com

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

}, keywords = {Computational Biology, Gene Expression Profiling, Gene ontology, Gene Regulatory Networks, High-Throughput Nucleotide Sequencing, Humans, MicroRNAs, Neoplasms, Reproducibility of Results}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btw334}, author = {Garcia-Garcia, Francisco and Panadero, Joaquin and Dopazo, Joaquin and Montaner, David} } @article {446, title = {Integrating transcriptomic and metabolomic analysis to understand natural leaf senescence in sunflower.}, journal = {Plant Biotechnol J}, volume = {14}, year = {2016}, month = {2016 Feb}, pages = {719-34}, abstract = {

Leaf senescence is a complex process, which has dramatic consequences on crop yield. In sunflower, gap between potential and actual yields reveals the economic impact of senescence. Indeed, sunflower plants are incapable of maintaining their green leaf area over sustained periods. This study characterizes the leaf senescence process in sunflower through a systems biology approach integrating transcriptomic and metabolomic analyses: plants being grown under both glasshouse and field conditions. Our results revealed a correspondence between profile changes detected at the molecular, biochemical and physiological level throughout the progression of leaf senescence measured at different plant developmental stages. Early metabolic changes were detected prior to anthesis and before the onset of the first senescence symptoms, with more pronounced changes observed when physiological and molecular variables were assessed under field conditions. During leaf development, photosynthetic activity and cell growth processes decreased, whereas sucrose, fatty acid, nucleotide and amino acid metabolisms increased. Pathways related to nutrient recycling processes were also up-regulated. Members of the NAC, AP2-EREBP, HB, bZIP and MYB transcription factor families showed high expression levels, and their expression level was highly correlated, suggesting their involvement in sunflower senescence. The results of this study thus contribute to the elucidation of the molecular mechanisms involved in the onset and progression of leaf senescence in sunflower leaves as well as to the identification of candidate genes involved in this process.

}, keywords = {Gas Chromatography-Mass Spectrometry, Gene Expression Profiling, Gene Expression Regulation, Plant, Gene ontology, Genes, Plant, Helianthus, Ions, metabolomics, Oligonucleotide Array Sequence Analysis, Plant Leaves, Principal Component Analysis, RNA, Messenger, Transcription Factors}, issn = {1467-7652}, doi = {10.1111/pbi.12422}, author = {Moschen, Sebasti{\'a}n and Bengoa Luoni, Sof{\'\i}a and Di Rienzo, Julio A and Caro, Mar{\'\i}a Del Pilar and Tohge, Takayuki and Watanabe, Mutsumi and Hollmann, Julien and Gonzalez, Sergio and Rivarola, M{\'a}ximo and Garcia-Garcia, Francisco and Dopazo, Joaquin and Hopp, Horacio Esteban and Hoefgen, Rainer and Fernie, Alisdair R and Paniego, Norma and Fernandez, Paula and Heinz, Ruth A} } @article {454, title = {The transcriptomics of an experimentally evolved plant-virus interaction.}, journal = {Sci Rep}, volume = {6}, year = {2016}, month = {2016 04 26}, pages = {24901}, abstract = {

Models of plant-virus interaction assume that the ability of a virus to infect a host genotype depends on the matching between virulence and resistance genes. Recently, we evolved tobacco etch potyvirus (TEV) lineages on different ecotypes of Arabidopsis thaliana, and found that some ecotypes selected for specialist viruses whereas others selected for generalists. Here we sought to evaluate the transcriptomic basis of such relationships. We have characterized the transcriptomic responses of five ecotypes infected with the ancestral and evolved viruses. Genes and functional categories differentially expressed by plants infected with local TEV isolates were identified, showing heterogeneous responses among ecotypes, although significant parallelism existed among lineages evolved in the same ecotype. Although genes involved in immune responses were altered upon infection, other functional groups were also pervasively over-represented, suggesting that plant resistance genes were not the only drivers of viral adaptation. Finally, the transcriptomic consequences of infection with the generalist and specialist lineages were compared. Whilst the generalist induced very similar perturbations in the transcriptomes of the different ecotypes, the perturbations induced by the specialist were divergent. Plant defense mechanisms were activated when the infecting virus was specialist but they were down-regulated when infecting with generalist.

}, keywords = {Arabidopsis, Ecotype, Gene Expression Profiling, Host-Pathogen Interactions, Potyvirus}, issn = {2045-2322}, doi = {10.1038/srep24901}, author = {Hillung, Julia and Garcia-Garcia, Francisco and Dopazo, Joaquin and Cuevas, Jos{\'e} M and Elena, Santiago F} } @article {492, title = {Permanent cardiac sarcomere changes in a rabbit model of intrauterine growth restriction.}, journal = {PLoS One}, volume = {9}, year = {2014}, month = {2014}, pages = {e113067}, abstract = {

BACKGROUND: Intrauterine growth restriction (IUGR) induces fetal cardiac remodelling and dysfunction, which persists postnatally and may explain the link between low birth weight and increased cardiovascular mortality in adulthood. However, the cellular and molecular bases for these changes are still not well understood. We tested the hypothesis that IUGR is associated with structural and functional gene expression changes in the fetal sarcomere cytoarchitecture, which remain present in adulthood.

METHODS AND RESULTS: IUGR was induced in New Zealand pregnant rabbits by selective ligation of the utero-placental vessels. Fetal echocardiography demonstrated more globular hearts and signs of cardiac dysfunction in IUGR. Second harmonic generation microscopy (SHGM) showed shorter sarcomere length and shorter A-band and thick-thin filament interaction lengths, that were already present in utero and persisted at 70 postnatal days (adulthood). Sarcomeric M-band (GO: 0031430) functional term was over-represented in IUGR fetal hearts.

CONCLUSION: The results suggest that IUGR induces cardiac dysfunction and permanent changes on the sarcomere.

}, keywords = {Animals, biomarkers, Blood Pressure, Body Weight, Disease Models, Animal, Echocardiography, Female, Fetal Growth Retardation, Fetal Heart, Fetus, Gene Expression Profiling, Organ Size, Placenta, Pregnancy, Rabbits, Sarcomeres}, issn = {1932-6203}, doi = {10.1371/journal.pone.0113067}, author = {Torre, Iratxe and Gonz{\'a}lez-Tendero, Anna and Garc{\'\i}a-Ca{\~n}adilla, Patricia and Crispi, F{\'a}tima and Garcia-Garcia, Francisco and Bijnens, Bart and Iruretagoyena, Igor and Dopazo, Joaquin and Amat-Rold{\'a}n, Ivan and Gratac{\'o}s, Eduard} } @article {505, title = {Mammosphere formation in breast carcinoma cell lines depends upon expression of E-cadherin.}, journal = {PLoS One}, volume = {8}, year = {2013}, month = {2013}, pages = {e77281}, abstract = {

Tumors are heterogeneous at the cellular level where the ability to maintain tumor growth resides in discrete cell populations. Floating sphere-forming assays are broadly used to test stem cell activity in tissues, tumors and cell lines. Spheroids are originated from a small population of cells with stem cell features able to grow in suspension culture and behaving as tumorigenic in mice. We tested the ability of eleven common breast cancer cell lines representing the major breast cancer subtypes to grow as mammospheres, measuring the ability to maintain cell viability upon serial non-adherent passage. Only MCF7, T47D, BT474, MDA-MB-436 and JIMT1 were successfully propagated as long-term mammosphere cultures, measured as the increase in the number of viable cells upon serial non-adherent passages. Other cell lines tested (SKBR3, MDA-MB-231, MDA-MB-468 and MDA-MB-435) formed cell clumps that can be disaggregated mechanically, but cell viability drops dramatically on their second passage. HCC1937 and HCC1569 cells formed typical mammospheres, although they could not be propagated as long-term mammosphere cultures. All the sphere forming lines but MDA-MB-436 express E-cadherin on their surface. Knock down of E-cadherin expression in MCF-7 cells abrogated its ability to grow as mammospheres, while re-expression of E-cadherin in SKBR3 cells allow them to form mammospheres. Therefore, the mammosphere assay is suitable to reveal stem like features in breast cancer cell lines that express E-cadherin.

}, keywords = {Breast Neoplasms, Cadherins, Cell Line, Tumor, Cell Proliferation, Cluster Analysis, Female, gene expression, Gene Expression Profiling, Gene Expression Regulation, Neoplastic, Gene Knockdown Techniques, Humans, MCF-7 Cells, Neoplastic Stem Cells, Spheroids, Cellular, Tumor Cells, Cultured}, issn = {1932-6203}, doi = {10.1371/journal.pone.0077281}, author = {Manuel Iglesias, Juan and Beloqui, Izaskun and Garcia-Garcia, Francisco and Leis, Olatz and Vazquez-Martin, Alejandro and Eguiara, Arrate and Cufi, Silvia and Pavon, Andres and Menendez, Javier A and Dopazo, Joaquin and Martin, Angel G} } @article {507, title = {Novel genes detected by transcriptional profiling from whole-blood cells in patients with early onset of acute coronary syndrome.}, journal = {Clin Chim Acta}, volume = {421}, year = {2013}, month = {2013 Jun 05}, pages = {184-90}, abstract = {

BACKGROUND: Genome-wide expression analysis using microarrays has been used as a research strategy to discovery new biomarkers and candidate genes for a number of diseases. We aim to find new biomarkers for the prediction of acute coronary syndrome (ACS) with a differentially expressed mRNA profiling approach using whole genomic expression analysis in a peripheral blood cell model from patients with early ACS.

METHODS AND RESULTS: This study was carried out in two phases. On phase 1 a restricted clinical criteria (ACS-Ph1, n=9 and CG-Ph1, n=6) was used in order to select potential mRNA biomarkers candidates. A subsequent phase 2 study was performed using selected phase 1 markers analyzed by RT-qPCR using a larger and independent casuistic (ACS-Ph2, n=74 and CG-Ph2, n=41). A total of 549 genes were found to be differentially expressed in the first 48 h after the ACS-Ph1. Technical and biological validation further confirmed that ALOX15, AREG, BCL2A1, BCL2L1, CA1, COX7B, ECHDC3, IL18R1, IRS2, KCNE1, MMP9, MYL4 and TREML4, are differentially expressed in both phases of this study.

CONCLUSIONS: Transcriptomic analysis by microarray technology demonstrated differential expression during a 48 h time course suggesting a potential use of some of these genes as biomarkers for very early stages of ACS, as well as for monitoring early cardiac ischemic recovery.

}, keywords = {Acute Coronary Syndrome, Acute-Phase Proteins, Adult, biomarkers, Blood Cells, Early Diagnosis, gene expression, Gene Expression Profiling, Humans, Male, Middle Aged, Oligonucleotide Array Sequence Analysis, RNA, Messenger, Transcriptome}, issn = {1873-3492}, doi = {10.1016/j.cca.2013.03.011}, author = {Silbiger, Vivian N and Luchessi, Andr{\'e} D and Hirata, Ros{\'a}rio D C and Lima-Neto, L{\'\i}dio G and Cavichioli, D{\'e}bora and Carracedo, {\'A}ngel and Bri{\'o}n, Maria and Dopazo, Joaquin and Garcia-Garcia, Francisco and Dos Santos, Elizabete S and Ramos, Rui F and Sampaio, Marcelo F and Armaganijan, Dikran and Sousa, Amanda G M R and Hirata, Mario H} } @article {501, title = {Role of CPI-17 in restoring skin homoeostasis in cutaneous field of cancerization: effects of topical application of a film-forming medical device containing photolyase and UV filters.}, journal = {Exp Dermatol}, volume = {22}, year = {2013}, month = {2013 Jul}, pages = {494-6}, abstract = {

Cutaneous field of cancerization (CFC) is caused in part by the carcinogenic effect of the cyclobutane pyrimidine dimers CPD and 6-4 photoproducts (6-4PPs). Photoreactivation is carried out by photolyases which specifically recognize and repair both photoproducts. The study evaluates the molecular effects of topical application of a film-forming medical device containing photolyase and UV filters on the precancerous field in AK from seven patients. Skin improvement after treatment was confirmed in all patients by histopathological and molecular assessment. A gene set analysis showed that skin recovery was associated with biological processes involved in tissue homoeostasis and cell maintenance. The CFC response was associated with over-expression of the CPI-17 gene, and a dependence on the initial expression level was observed (P~=~0.001). Low CPI-17 levels were directly associated with pro-inflammatory genes such as TNF (P~=~0.012) and IL-1B (P~=~0.07). Our results suggest a role for CPI-17 in restoring skin homoeostasis in CFC lesions.

}, keywords = {Administration, Topical, Adult, Aged, Aged, 80 and over, Biopsy, Deoxyribodipyrimidine Photo-Lyase, Female, Gene Expression Profiling, Gene Expression Regulation, Enzymologic, Gene Expression Regulation, Neoplastic, Homeostasis, Humans, Inflammation, Intracellular Signaling Peptides and Proteins, Liposomes, Male, Middle Aged, Muscle Proteins, Phenotype, Phosphoprotein Phosphatases, Reactive Oxygen Species, Skin, Skin Neoplasms, Ultraviolet Rays}, issn = {1600-0625}, doi = {10.1111/exd.12177}, author = {Puig-Butille, Joan Anton and Malvehy, Josep and Potrony, Miriam and Trullas, Carles and Garcia-Garcia, Francisco and Dopazo, Joaquin and Puig, Susana} } @article {516, title = {IL1β induces mesenchymal stem cells migration and leucocyte chemotaxis through NF-κB.}, journal = {Stem Cell Rev Rep}, volume = {8}, year = {2012}, month = {2012 Sep}, pages = {905-16}, abstract = {

Mesenchymal stem cells are often transplanted into inflammatory environments where they are able to survive and modulate host immune responses through a poorly understood mechanism. In this paper we analyzed the responses of MSC to IL-1β: a representative inflammatory mediator. Microarray analysis of MSC treated with IL-1β revealed that this cytokine activateds a set of genes related to biological processes such as cell survival, cell migration, cell adhesion, chemokine production, induction of angiogenesis and modulation of the immune response. Further more detailed analysis by real-time PCR and functional assays revealed that IL-1β mainly increaseds the production of chemokines such as CCL5, CCL20, CXCL1, CXCL3, CXCL5, CXCL6, CXCL10, CXCL11 and CX(3)CL1, interleukins IL-6, IL-8, IL23A, IL32, Toll-like receptors TLR2, TLR4, CLDN1, metalloproteins MMP1 and MMP3, growth factors CSF2 and TNF-α, together with adhesion molecules ICAM1 and ICAM4. Functional analysis of MSC proliferation, migration and adhesion to extracellular matrix components revealed that IL-1β did not affect proliferation but also served to induce the secretion of trophic factors and adhesion to ECM components such as collagen and laminin. IL-1β treatment enhanced the ability of MSC to recruit monocytes and granulocytes in vitro. Blockade of NF-κβ transcription factor activation with IκB kinase beta (IKKβ) shRNA impaired MSC migration, adhesion and leucocyte recruitment, induced by IL-1β demonstrating that NF-κB pathway is an important downstream regulator of these responses. These findings are relevant to understanding the biological responses of MSC to inflammatory environments.

}, keywords = {Cell Adhesion, Cell Movement, Cell Proliferation, Chemokines, Chemotaxis, Leukocyte, Collagen, Fibronectins, Gene Expression Profiling, Gene Knockdown Techniques, HEK293 Cells, Humans, I-kappa B Kinase, Inflammation Mediators, Intercellular Signaling Peptides and Proteins, Interleukin-1beta, Laminin, Leukocytes, Mesenchymal Stem Cells, NF-kappa B, Oligonucleotide Array Sequence Analysis, RNA Interference, Signal Transduction}, issn = {2629-3277}, doi = {10.1007/s12015-012-9364-9}, author = {Carrero, Rub{\'e}n and Cerrada, Inmaculada and Lled{\'o}, Elisa and Dopazo, Joaquin and Garcia-Garcia, Francisco and Rubio, Mari-Paz and Trigueros, C{\'e}sar and Dorronsoro, Akaitz and Ruiz-Sauri, Amparo and Montero, Jos{\'e} Anastasio and Sep{\'u}lveda, Pilar} } @article {532, title = {Early transcriptional defense responses in Arabidopsis cell suspension culture under high-light conditions.}, journal = {Plant Physiol}, volume = {156}, year = {2011}, month = {2011 Jul}, pages = {1439-56}, abstract = {

The early transcriptional defense responses and reactive oxygen species (ROS) production in Arabidopsis (Arabidopsis thaliana) cell suspension culture (ACSC), containing functional chloroplasts, were examined at high light (HL). The transcriptional analysis revealed that most of the ROS markers identified among the 449 transcripts with significant differential expression were transcripts specifically up-regulated by singlet oxygen ((1)O(2)). On the contrary, minimal correlation was established with transcripts specifically up-regulated by superoxide radical or hydrogen peroxide. The transcriptional analysis was supported by fluorescence microscopy experiments. The incubation of ACSC with the (1)O(2) sensor green reagent and 2{\textquoteright},7{\textquoteright}-dichlorofluorescein diacetate showed that the 30-min-HL-treated cultures emitted fluorescence that corresponded with the production of (1)O(2) but not of hydrogen peroxide. Furthermore, the in vivo photodamage of the D1 protein of photosystem II indicated that the photogeneration of (1)O(2) took place within the photosystem II reaction center. Functional enrichment analyses identified transcripts that are key components of the ROS signaling transduction pathway in plants as well as others encoding transcription factors that regulate both ROS scavenging and water deficit stress. A meta-analysis examining the transcriptional profiles of mutants and hormone treatments in Arabidopsis showed a high correlation between ACSC at HL and the fluorescent mutant family of Arabidopsis, a producer of (1)O(2) in plastids. Intriguingly, a high correlation was also observed with ABA deficient1 and more axillary growth4, two mutants with defects in the biosynthesis pathways of two key (apo)carotenoid-derived plant hormones (i.e. abscisic acid and strigolactones, respectively). ACSC has proven to be a valuable system for studying early transcriptional responses to HL stress.

}, keywords = {Arabidopsis, Blotting, Western, Cell Culture Techniques, Cells, Cultured, Chloroplasts, Cluster Analysis, Gene Expression Profiling, Gene Expression Regulation, Plant, Hydrogen Peroxide, Light, mutation, Oligonucleotide Array Sequence Analysis, Photosystem II Protein Complex, Plant Growth Regulators, Reproducibility of Results, Reverse Transcriptase Polymerase Chain Reaction, RNA, Messenger, Signal Transduction, Stress, Physiological, Transcription, Genetic}, issn = {1532-2548}, doi = {10.1104/pp.111.177766}, author = {Gonz{\'a}lez-P{\'e}rez, Sergio and Guti{\'e}rrez, Jorge and Garcia-Garcia, Francisco and Osuna, Daniel and Dopazo, Joaquin and Lorenzo, Oscar and Revuelta, Jos{\'e} L and Arellano, Juan B} } @article {580, title = {Functional assessment of time course microarray data.}, journal = {BMC Bioinformatics}, volume = {10 Suppl 6}, year = {2009}, month = {2009 Jun 16}, pages = {S9}, abstract = {

MOTIVATION: Time-course microarray experiments study the progress of gene expression along time across one or several experimental conditions. Most developed analysis methods focus on the clustering or the differential expression analysis of genes and do not integrate functional information. The assessment of the functional aspects of time-course transcriptomics data requires the use of approaches that exploit the activation dynamics of the functional categories to where genes are annotated.

METHODS: We present three novel methodologies for the functional assessment of time-course microarray data. i) maSigFun derives from the maSigPro method, a regression-based strategy to model time-dependent expression patterns and identify genes with differences across series. maSigFun fits a regression model for groups of genes labeled by a functional class and selects those categories which have a significant model. ii) PCA-maSigFun fits a PCA model of each functional class-defined expression matrix to extract orthogonal patterns of expression change, which are then assessed for their fit to a time-dependent regression model. iii) ASCA-functional uses the ASCA model to rank genes according to their correlation to principal time expression patterns and assess functional enrichment on a GSA fashion. We used simulated and experimental datasets to study these novel approaches. Results were compared to alternative methodologies.

RESULTS: Synthetic and experimental data showed that the different methods are able to capture different aspects of the relationship between genes, functions and co-expression that are biologically meaningful. The methods should not be considered as competitive but they provide different insights into the molecular and functional dynamic events taking place within the biological system under study.

}, keywords = {Computer Simulation, Gene Expression Profiling, Oligonucleotide Array Sequence Analysis, Time Factors}, issn = {1471-2105}, doi = {10.1186/1471-2105-10-S6-S9}, author = {Nueda, Maria Jos{\'e} and Sebasti{\'a}n, Patricia and Tarazona, Sonia and Garcia-Garcia, Francisco and Dopazo, Joaquin and Ferrer, Alberto and Conesa, Ana} } @article {590, title = {Direct functional assessment of the composite phenotype through multivariate projection strategies.}, journal = {Genomics}, volume = {92}, year = {2008}, month = {2008 Dec}, pages = {373-83}, abstract = {

We present a novel approach for the analysis of transcriptomics data that integrates functional annotation of gene sets with expression values in a multivariate fashion, and directly assesses the relation of functional features to a multivariate space of response phenotypical variables. Multivariate projection methods are used to obtain new correlated variables for a set of genes that share a given function. These new functional variables are then related to the response variables of interest. The analysis of the principal directions of the multivariate regression allows for the identification of gene function features correlated with the phenotype. Two different transcriptomics studies are used to illustrate the statistical and interpretative aspects of the methodology. We demonstrate the superiority of the proposed method over equivalent approaches.

}, keywords = {Breast Neoplasms, Computational Biology, Databases, Genetic, Female, Gene Expression Profiling, Humans, Mathematical Computing, Multivariate Analysis, Phenotype}, issn = {1089-8646}, doi = {10.1016/j.ygeno.2008.05.015}, author = {Conesa, Ana and Bro, Rasmus and Garcia-Garcia, Francisco and Prats, Jos{\'e} Manuel and G{\"o}tz, Stefan and Kjeldahl, Karin and Montaner, David and Dopazo, Joaquin} } @article {593, title = {GEPAS, a web-based tool for microarray data analysis and interpretation.}, journal = {Nucleic Acids Res}, volume = {36}, year = {2008}, month = {2008 Jul 01}, pages = {W308-14}, abstract = {

Gene Expression Profile Analysis Suite (GEPAS) is one of the most complete and extensively used web-based packages for microarray data analysis. During its more than 5 years of activity it has continuously been updated to keep pace with the state-of-the-art in the changing microarray data analysis arena. GEPAS offers diverse analysis options that include well established as well as novel algorithms for normalization, gene selection, class prediction, clustering and functional profiling of the experiment. New options for time-course (or dose-response) experiments, microarray-based class prediction, new clustering methods and new tests for differential expression have been included. The new pipeliner module allows automating the execution of sequential analysis steps by means of a simple but powerful graphic interface. An extensive re-engineering of GEPAS has been carried out which includes the use of web services and Web 2.0 technology features, a new user interface with persistent sessions and a new extended database of gene identifiers. GEPAS is nowadays the most quoted web tool in its field and it is extensively used by researchers of many countries and its records indicate an average usage rate of 500 experiments per day. GEPAS, is available at http://www.gepas.org.

}, keywords = {Computer Graphics, Dose-Response Relationship, Drug, Gene Expression Profiling, Internet, Kinetics, Oligonucleotide Array Sequence Analysis, Software}, issn = {1362-4962}, doi = {10.1093/nar/gkn303}, author = {T{\'a}rraga, Joaqu{\'\i}n and Medina, Ignacio and Carbonell, Jos{\'e} and Huerta-Cepas, Jaime and Minguez, Pablo and Alloza, Eva and Al-Shahrour, F{\'a}tima and Vegas-Azc{\'a}rate, Susana and Goetz, Stefan and Escobar, Pablo and Garcia-Garcia, Francisco and Conesa, Ana and Montaner, David and Dopazo, Joaquin} }