Enriched Gene Ontology Categories Research Articles

Array probe density and pathobiological relevant CpG calling bias in human disease and physiological DNA methylation profiling.

The HumanMethylation450 BeadChip array (450K; Infinium) is a widely used tool in epigenomics. A recognized concern in the 450K platform is the potential effect of the number of probes/gene (PG) on ranking differentially methylated (DM) CpGs (DM-CpGs) before testing for enrichment of gene ontology categories. We previously showed in a fatty acid (FA)-induced DNA methylation profiling study that when DM-CpGs are ranked by the number of called DM-CpGs-to-PG ratio, the 150 top-ranking gene list is enriched in pathways that overlap with the corresponding Affymetrix array-based expression data. In this study, a comparative analysis of thirteen 450K-based studies representing FA-stimulated cellular models, aging, diseased and normal tissues, revealed that the 150 top-ranking DM-CpGs are in high PG genes. This points to a significant false-negative rate in the low PG gene set when delta-beta-based ranking is performed. We show that PG is not related to the density of methylation-prone sites, as it does not follow gene length or GC content. Conversely, ranking genes by the number of DM-CpGs-to-PG ratio and analysing the 150 top-ranking entries yields significantly enriched gene disease- or tissue-specific function categories that are increased both in number and in the degree of overlap with expression data compared with delta-beta-only ranking or to the previously published gometh-based pipeline. The 15 top-ranking loci list is also significantly enriched in non-coding RNAs, a greatly underrepresented transcript type in 450K. In summary, the proposed simple normalization method yields pathobiologically relevant DM-CpGs. This method is relevant for the newly developed MethylationEPIC (Infinium) microarray.

Analysis of transcriptome in hickory (Carya cathayensis), and uncover the dynamics in the hormonal signaling pathway during graft process.

BackgroundHickory (Carya cathayensis), a woody plant with high nutritional and economic value, is widely planted in China. Due to its long juvenile phase, grafting is a useful technique for large-scale cultivation of hickory. To reveal the molecular mechanism during the graft process, we sequenced the transcriptomes of graft union in hickory.ResultsIn our study, six RNA-seq libraries yielded a total of 83,676,860 clean short reads comprising 4.19 Gb of sequence data. A large number of differentially expressed genes (DEGs) at three time points during the graft process were identified. In detail, 777 DEGs in the 7 d vs 0 d (day after grafting) comparison were classified into 11 enriched Gene Ontology (GO) categories, and 262 DEGs in the 14 d vs 0 d comparison were classified into 15 enriched GO categories. Furthermore, an overview of the PPI network was constructed by these DEGs. In addition, 20 genes related to the auxin-and cytokinin-signaling pathways were identified, and some were validated by qRT-PCR analysis.ConclusionsOur comprehensive analysis provides basic information on the candidate genes and hormone signaling pathways involved in the graft process in hickory and other woody plants.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3182-4) contains supplementary material, which is available to authorized users.

Comparative proteomics of cucurbit phloem indicates both unique and shared sets of proteins

Cucurbits are well-studied models for phloem biology but unusually possess both fascicular phloem (FP) within vascular bundles and additional extrafascicular phloem (EFP). Although the functional differences between the two systems are not yet clear, sugar analysis and limited protein profiling have established that FP and EFP have divergent compositions. Here we report a detailed comparative proteomics study of FP and EFP in two cucurbits, pumpkin and cucumber. We re-examined the sites of exudation by video microscopy, and confirmed that in both species, the spontaneous exudate following tissue cutting derives almost exclusively from EFP. Comparative gel electrophoresis and mass spectrometry-based proteomics of exudates, sieve element contents and microdissected stem tissues established that EFP and FP profiles are highly dissimilar, and that there are also species differences. Searches against cucurbit databases enabled identification of more than 300 FP proteins from each species. Few of the detected proteins (about 10%) were shared between the sieve element contents of FP and EFP, and enriched Gene Ontology categories also differed. To explore quantitative differences in the proteomes, we developed multiple reaction monitoring methods for cucumber proteins that are representative markers for FP or EFP and assessed exudate composition at different times after tissue cutting. Based on failure to detect FP markers in exudate samples, we conclude that FP is blocked very rapidly and therefore makes a minimal contribution to the exudates. Overall, the highly divergent contents of FP and EFP indicate that they are substantially independent vascular compartments.

Comparative Proteomic Analysis of the Response of Maize (Zea mays L.) Leaves to Long Photoperiod Condition.

Maize (Zea mays L.), an important industrial material and food source, shows an astonishing environmental adaptation. A remarkable feature of its post-domestication adaptation from tropical to temperate environments is adaptation to a long photoperiod (LP). Many photoperiod-related genes have been identified in previous transcriptomics analysis, but proteomics shows less evidence for this mechanism of photoperiod response. In this study, we sampled newly expanded leaves of maize at the three- and six-leaf stages from an LP-sensitive introgression line H496, the donor CML288, LP-insensitive inbred line, and recurrent parent Huangzao4 (HZ4) grown under long days (15 h light and 9 h dark). To characterize the proteomic changes in response to LP, the iTRAQ-labeling method was used to determine the proteome profiles of plants exposed to LP. A total of 943 proteins differentially expressed at the three- and six-leaf stages in HZ4 and H496 were identified. Functional analysis was performed by which the proteins were classified into stress defense, signal transduction, carbohydrate metabolism, protein metabolism, energy production, and transport functional groups using the WEGO online tool. The enriched gene ontology categories among the identified proteins were identified statistically with the Cytoscape plugin ClueGO + Cluepedia. Twenty Gene Ontology terms showed the highest significance, including those associated with protein processing in the endoplasmic reticulum, splicesome, ribosome, glyoxylate, dicarboxylate metabolism, L-malate dehydrogenase activity, and RNA transport. In addition, for subcellular location, all proteins showed significant enrichment of the mitochondrial outer membrane. The sugars producted by photosynthesis in plants are also a pivotal metabolic output in the circadian regulation. The results permit the prediction of several crucial proteins to photoperiod response and provide a foundation for further study of the influence of LP treatments on the circadian response in short-day plants.

Transcriptomic network analyses of leaf dehydration responses identify highly connected ABA and ethylene signaling hubs in three grapevine species differing in drought tolerance.

BackgroundGrapevine is a major food crop that is affected by global climate change. Consistent with field studies, dehydration assays of grapevine leaves can reveal valuable information of the plant’s response at physiological, transcript, and protein levels. There are well-known differences in grapevine rootstocks responses to dehydration. We used time-series transcriptomic approaches combined with network analyses to elucidate and identify important physiological processes and network hubs that responded to dehydration in three different grapevine species differing in their drought tolerance.ResultsTranscriptomic analyses of the leaves of Cabernet Sauvignon, Riparia Gloire, and Ramsey were evaluated at different times during a 24-h controlled dehydration. Analysis of variance (ANOVA) revealed that approximately 11,000 transcripts changed significantly with respect to the genotype x treatment interaction term and approximately 6000 transcripts changed significantly according to the genotype x treatment x time interaction term indicating massive differential changes in gene expression over time. Standard analyses determined substantial effects on the transcript abundance of genes involved in the metabolism and signaling of two known plant stress hormones, abscisic acid (ABA) and ethylene. ABA and ethylene signaling maps were constructed and revealed specific changes in transcript abundance that were associated with the known drought tolerance of the genotypes including genes such as VviABI5, VviABF2, VviACS2, and VviWRKY22. Weighted-gene coexpression network analysis (WGCNA) confirmed these results. In particular, WGCNA identified 30 different modules, some of which had highly enriched gene ontology (GO) categories for photosynthesis, phenylpropanoid metabolism, ABA and ethylene signaling. The ABA signaling transcription factors, VviABI5 and VviABF2, were highly connected hubs in two modules, one being enriched in gaseous transport and the other in ethylene signaling. VviABI5 was distinctly correlated with an early response and high expression for the drought tolerant Ramsey and with little response from the drought sensitive Riparia Gloire. These ABA signaling transcription factors were highly connected to VviSnRK1 and other gene hubs associated with sugar, ethylene and ABA signaling.ConclusionA leaf dehydration assay provided transcriptomic evidence for differential leaf responses to dehydration between genotypes differing in their drought tolerance. WGCNA proved to be a powerful network analysis approach; it identified 30 distinct modules (networks) with highly enriched GO categories and enabled the identification of gene hubs in these modules. Some of these genes were highly connected hubs in both the ABA and ethylene signaling pathways, supporting the hypothesis that there is substantial crosstalk between the two hormone pathways. This study identifies solid gene candidates for future investigations of drought tolerance in grapevine.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0804-6) contains supplementary material, which is available to authorized users.

RNA-Seq Profiling of Intact and Enucleated Oocyte SCNT Embryos Reveals the Role of Pig Oocyte Nucleus in Somatic Reprogramming.

The specific molecular mechanisms involved in somatic reprogramming remain unidentified. Removal of the oocyte genome is one of the primary causes of developmental failure in cloned embryos, whereas intact oocyte shows stronger reprogramming capability than enucleated oocyte. To identify the reason for the low efficiency of cloning and elucidate the mechanisms involved in somatic reprogramming by the oocyte nucleus, we injected pig cumulus cells into 539 intact MII oocytes and 461 enucleated MII oocytes. Following activation, 260 polyploidy embryos developed to the blastocyst stage whereas only 93 traditionally cloned embryos (48.2% vs. 20.2%, P < 0.01) reached blastocyst stage. Blastocysts generated from intact oocytes also had more cells than those generated from enucleated oocytes (60.70 vs. 46.65, P < 0.01). To identify the genes that contribute to this phenomenon, two early embryos in 2-cell and 4-cell stages were collected for single-cell RNA sequencing. The two kinds of embryos were found to have dramatically different transcriptome profiles. Intact oocyte nuclear transfer embryos showed 1,738 transcripts that were up-regulated relative to enucleated cloned embryos at the 2-cell stage and 728 transcripts that were down-regulated (|log2Ratio| ≥ 5). They showed 2,941 transcripts that were up-regulated during the 4-cell stage and 1,682 that were down-regulated (|log2Ratio| ≥ 5). The most significantly enriched gene ontology categories were those involved in the regulation of binding, catalytic activity, and molecular transducer activity. Other genes that were notably up-regulated and expressed in intact oocyte nuclear transfer embryos were metabolic process. This study provides a comprehensive profile of the differences in gene expression between intact oocyte nuclear transfer embryos and traditional nuclear transfer embryos. This work thus paves the way for further research on the mechanisms underlying somatic reprogramming by oocytes.

Identification of Differentially Expressed Genes Related to Dehydration Resistance in a Highly Drought-Tolerant Pear, Pyrus betulaefolia, as through RNA-Seq.

Drought is a major abiotic stress that affects plant growth, development and productivity. Pear is one of the most important deciduous fruit trees in the world, but the mechanisms of drought tolerance in this plant are still unclear. To better understand the molecular basis regarding drought stress response, RNA-seq was performed on samples collected before and after dehydration in Pyrus betulaefolia. In total, 19,532 differentially expressed genes (DEGs) were identified. These genes were annotated into 144 Gene Ontology (GO) terms and 18 clusters of orthologous groups (COG) involved in 129 Kyoto Encyclopedia of Genes and Genomes (KEGG) defined pathways. These DEGs comprised 49 (26 up-regulated, 23 down-regulated), 248 (166 up-regulated, 82 down-regulated), 3483 (1295 up-regulated, 2188 down-regulated), 1455 (1065 up-regulated, 390 down-regulated) genes from the 1 h, 3 h and 6 h dehydration-treated samples and a 24 h recovery samples, respectively. RNA-seq was validated by analyzing the expresson patterns of randomly selected 16 DEGs by quantitative real-time PCR. Photosynthesis, signal transduction, innate immune response, protein phosphorylation, response to water, response to biotic stimulus, and plant hormone signal transduction were the most significantly enriched GO categories amongst the DEGs. A total of 637 transcription factors were shown to be dehydration responsive. In addition, a number of genes involved in the metabolism and signaling of hormones were significantly affected by the dehydration stress. This dataset provides valuable information regarding the Pyrus betulaefolia transcriptome changes in response to dehydration and may promote identification and functional analysis of potential genes that could be used for improving drought tolerance via genetic engineering of non-model, but economically-important, perennial species.

In vitro maturation alters gene expression in bovine oocytes.

Gene expression profiling of in vivo- and in vitro-matured bovine oocytes can identify transcripts related to the developmental potential of oocytes. Nonetheless, the effects of in vitro culturing oocytes are yet to be fully understood. We tested the effects of in vitro maturation on the transcript profile of oocytes collected from Bos taurus indicus cows. We quantified the expression of 1488 genes in in vivo- and in vitro-matured oocytes. Of these, 51 genes were up-regulated, whereas 56 were down-regulated (≥2-fold) in in vivo-matured oocytes in comparison with in vitro-matured oocytes. Quantitative real-time polymerase chain reaction (PCR) of nine genes confirmed the microarray results of differential expression between in vivo- and in vitro-matured oocytes (EZR, EPN1, PSEN2, FST, IGFBP3, RBBP4, STAT3, FDPS and IRS1). We interrogated the results for enrichment of Gene Ontology categories and overlap with protein-protein interactions. The results revealed that the genes altered by in vitro maturation are mostly related to the regulation of oocyte metabolism. Additionally, analysis of protein-protein interactions uncovered two regulatory networks affected by the in vitro culture system. We propose that the differentially expressed genes are candidates for biomarkers of oocyte competence. In vitro oocyte maturation can affect the abundance of specific transcripts and are likely to deplete the developmental competence.

Abstract P4-09-23: Gene expression signatures of microcalcifications among Taiwanese breast cancers

Abstract Introduction: Microcalcification is one of the most common radiological and pathological features of breast ductal carcinoma in situ (DCIS), and to a lesser extent, invasive breast cancer. In current study, we evaluated the transcriptional profiles associated with the phenomenon of ectopic mammary mineralization and a gene expression signature is derived. Materials and methods: a total of 109 consecutive breast invasive cancers were prospectively collected and assayed with Affymetrix Human Genome U133 Plus 2.0 microarrays. The presence of microcalcification was confirmed by histopathological examinations as well as reviews of pre-operative mammography. The associations of gene expression profiles with microcalcifications and relevant clinical features such as DCIS including both comedo/high-grade and non-comedo subtypes were tested. Results: Microcalcifications were presented in 84 (80%) of the study population as confirmed by pathological examination. Of these 84 patients, 81 (96%) were grown with coexistent DCIS microscopically, while only 8 (38%) of the 21 patients without concurrent microcalcifications, the invasive tumors were accompanied with DCIS (Chi-square test, P&amp;lt;0.001). In addition, high-grade (comedo) type DCIS were presented in 44 (54%) of the 81 cancers with microcalcifications whereas only 15% (n=2) of tumors with DCIS but without microcalcificaitons were of high-grade (comedo) type. There were 69 genes differentially expressed between breast cancers with and without microcalcifications (nominal P&amp;lt;0.001 with 10,000 random permutations), and 11 were associated with high-grade (comedo) type DCIS including APOD, CCDC183, SLMO1, SLC6A5, FMO1, QPRT and CES4A. The enriched Gene Ontology categories encompasses glycosaminoglycan, aminoglycan metabolic processes, Golgi apparatus cellular component and protein ubiquitination, indicating an active secretory process. The intersect (18 probesets) of microcalcificaion and DCIS-associated genes provided the best predictive accuracy of 82% with Bayesian compound covariate predictor. Performance of compound covariate predictor classifier:ClassSensitivitySpecificityPPVNPVWith microcalcification0.8260.2170.7980.25Without microcalcification0.2170.8260.250.798 Compared with mammography alone, the diagnostic accuracy of gene expression-based signature is much improved (cross-validated ROC AUC: 0.738). Discussion: Our study suggested that mammary microcalcification is not only the earlier detectable radiological finding for disease screening but the phenomenon itself resulted from distinct biological processes that constituted the molecular heterogeneity of human breast cancers. Further studies to evaluate the prognostic significance of microcalcifications are warranted. Citation Format: Huang C-C, Huang C-S, Tu S-H, Tsai M-L. Gene expression signatures of microcalcifications among Taiwanese breast cancers. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P4-09-23.

Deep sequencing-based characterization of transcriptome of trifoliate orange (Poncirus trifoliata (L.) Raf.) in response to cold stress.

BackgroundTrifoliate orange (Poncirus trifoliata (L.) Raf.) is extremely cold hardy after a full acclimation; however the underlying molecular mechanisms underlying this economically valuable trait remain poorly understood. In this study, global transcriptome profiles of trifoliate orange under cold conditions (4 °C) over a time course were generated by high-throughput sequencing.ResultsMore than 68 million high-quality reads were produced and assembled into a non-redundant data of 77,292 unigenes with an average length of 1112 bp (N50 = 1778 bp). Of these, 23,846 had significant sequence similarity to known genes and these were assigned to 61 gene ontology (GO) categories and 25 clusters of orthologous groups (COG) involved in 128 KEGG pathways. Sequences derived from cold-treated and control plants were mapped to the assembled transcriptome, resulting in the identification of 5549 differentially expressed genes (DEGs). These comprised 600 (462 up-regulated, 138 down-regulated), 2346 (1631 up-regulated, 715 down-regulated), and 5177 (2702 up-regulated, 2475 down-regulated) genes from the cold-treated samples at 6, 24 and 72 h, respectively. The accuracy of the RNA-seq derived transcript expression data was validated by analyzing the expression patterns of 17 DEGs by qPCR. Plant hormone signal transduction, plant-pathogen interaction, and secondary metabolism were the most significantly enriched GO categories amongst in the DEGs. A total of 60 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones, such as abscisic acid, ethylene and gibberellin, were affected by the cold stress. Meanwhile, levels of putrescine progressively increased under cold, which was consistent with up-regulation of an arginine decarboxylase gene.ConclusionsThis dataset provides valuable information regarding the trifoliate orange transcriptome changes in response to cold stress and may help guide future identification and functional analysis of genes that are importnatn for enhancing cold hardiness.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1629-7) contains supplementary material, which is available to authorized users.

BASiCS: Bayesian Analysis of Single-Cell Sequencing Data.

Single-cell mRNA sequencing can uncover novel cell-to-cell heterogeneity in gene expression levels in seemingly homogeneous populations of cells. However, these experiments are prone to high levels of unexplained technical noise, creating new challenges for identifying genes that show genuine heterogeneous expression within the population of cells under study. BASiCS (Bayesian Analysis of Single-Cell Sequencing data) is an integrated Bayesian hierarchical model where: (i) cell-specific normalisation constants are estimated as part of the model parameters, (ii) technical variability is quantified based on spike-in genes that are artificially introduced to each analysed cell’s lysate and (iii) the total variability of the expression counts is decomposed into technical and biological components. BASiCS also provides an intuitive detection criterion for highly (or lowly) variable genes within the population of cells under study. This is formalised by means of tail posterior probabilities associated to high (or low) biological cell-to-cell variance contributions, quantities that can be easily interpreted by users. We demonstrate our method using gene expression measurements from mouse Embryonic Stem Cells. Cross-validation and meaningful enrichment of gene ontology categories within genes classified as highly (or lowly) variable supports the efficacy of our approach.

Gene expression response of mice after a single dose of 137CS as an internal emitter.

Cesium-137 is a radionuclide of concern in fallout from reactor accidents or nuclear detonations. When ingested or inhaled, it can expose the entire body for an extended period of time, potentially contributing to serious health consequences ranging from acute radiation syndrome to increased cancer risks. To identify changes in gene expression that may be informative for detecting such exposure, and to begin examining the molecular responses involved, we have profiled global gene expression in blood of male C57BL/6 mice injected with 137CsCl. We extracted RNA from the blood of control or 137CsCl-injected mice at 2, 3, 5, 20 or 30 days after exposure. Gene expression was measured using Agilent Whole Mouse Genome Microarrays, and the data was analyzed using BRB-ArrayTools. Between 466-6,213 genes were differentially expressed, depending on the time after 137Cs administration. At early times (2-3 days), the majority of responsive genes were expressed above control levels, while at later times (20-30 days) most responding genes were expressed below control levels. Numerous genes were overexpressed by day 2 or 3, and then underexpressed by day 20 or 30, including many Tp53-regulated genes. The same pattern was seen among significantly enriched gene ontology categories, including those related to nucleotide binding, protein localization and modification, actin and the cytoskeleton, and in the integrin signaling canonical pathway. We compared the expression of several genes three days after 137CsCl injection and three days after an acute external gamma-ray exposure, and found that the internal exposure appeared to produce a more sustained response. Many common radiation-responsive genes are altered by internally administered 137Cs, but the gene expression pattern resulting from continued irradiation at a decreasing dose rate is extremely complex, and appears to involve a late reversal of much of the initial response.

Co-regulated gene expression by oestrogen receptor α and liver receptor homolog-1 is a feature of the oestrogen response in breast cancer cells

Oestrogen receptor α (ERα) is a nuclear receptor that is the driving transcription factor expressed in the majority of breast cancers. Recent studies have demonstrated that the liver receptor homolog-1 (LRH-1), another nuclear receptor, regulates breast cancer cell proliferation and promotes motility and invasion. To determine the mechanisms of LRH-1 action in breast cancer, we performed gene expression microarray analysis following RNA interference for LRH-1. Interestingly, gene ontology (GO) category enrichment analysis of LRH-1–regulated genes identified oestrogen-responsive genes as the most highly enriched GO categories. Remarkably, chromatin immunoprecipitation coupled to massively parallel sequencing (ChIP-seq) to identify genomic targets of LRH-1 showed LRH-1 binding at many ERα binding sites. Analysis of select binding sites confirmed regulation of ERα−regulated genes by LRH-1 through binding to oestrogen response elements, as exemplified by the TFF1/pS2 gene. Finally, LRH-1 overexpression stimulated ERα recruitment, while LRH-1 knockdown reduced ERα recruitment to ERα binding sites. Taken together, our findings establish a key role for LRH-1 in the regulation of ERα target genes in breast cancer cells and identify a mechanism in which co-operative binding of LRH-1 and ERα at oestrogen response elements controls the expression of oestrogen-responsive genes.

GeneMANIA Prediction Server 2013 Update

GeneMANIA (http://www.genemania.org) is a flexible user-friendly web interface for generating hypotheses about gene function, analyzing gene lists and prioritizing genes for functional assays. Given a query gene list, GeneMANIA extends the list with functionally similar genes that it identifies using available genomics and proteomics data. GeneMANIA also reports weights that indicate the predictive value of each selected data set for the query. GeneMANIA can also be used in a function prediction setting: given a query gene, GeneMANIA finds a small set of genes that are most likely to share function with that gene based on their interactions with it. Enriched Gene Ontology categories among this set can sometimes point to the function of the gene. Seven organisms are currently supported (Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster, Mus musculus, Homo sapiens, Rattus norvegicus and Saccharomyces cerevisiae), and hundreds of data sets have been collected from GEO, BioGRID, IRefIndex and I2D, as well as organism-specific functional genomics data sets. Users can customize their search by selecting specific data sets to query and by uploading their own data sets to analyze.

Integrated analysis of miRNA and mRNA reveals that acrolein modulates GPI anchor biosynthesis in human primary endothelial cells

Acrolein (ACR) is a highly reactive α, β-unsaturated carbonyl compound that is related to several diabetic vascular complications. ACR, a ubiquitous aldehyde pollutant, is a product of incomplete combustion reactions and a major component of cigarette smoke. ACR is also endogenously generated during lipid peroxidation. Glycosylphosphatidylinositol (GPI)-anchored proteins is a membrane bound protein, which have a variety of cellular functions. It is reported that some of them showed the reduced level under diabetic conditions. In this study, we explored whether ACR may affect GPI anchor related gene expression in human umbilical vein endothelial cells (HUVECs) by modulation of miRNA expression that play a significant role in the post-transcriptional regulation of genes. We performed global gene expression analysis, acquired altered expression levels of 27 miRNAs and 21 mRNAs that are relevant to GPI anchored proteins, and observed their anti-correlations. Gene Ontology (GO) terms that are closely associated with GPI anchor biosynthesis were highest when analysis was based on biological processes characterized by enriched GO categories. Our results revealed that ACR may control miRNAs to regulate genes involving in GPI anchor biosynthesis in human endothelial cells.

Dandruff/seborrhoeic dermatitis is characterized by an inflammatory genomic signature and possible immune dysfunction: transcriptional analysis of the condition and treatment effects of zinc pyrithione

Dandruff/seborrhoeic dermatitis is a common scalp condition that is characterized by flakes, pruritus and sometimes mild erythema. These symptoms reflect tissue level events that are poorly understood at the molecular level. The purpose of this work was: (i) to compare gene expression profiles in subjects with dandruff vs. those of subjects without dandruff to determine the key physiological disruptions manifest in the condition; and (ii) to determine the effect on this profile of treatment with a shampoo containing potentiated zinc pyrithione (ZPT). In study 1, scalp biopsies were taken from 16 normal subjects and from involved and uninvolved sites in 15 subjects with dandruff. In study 2, 30 subjects with dandruff were treated for 3 weeks with a commercial ZPT shampoo (n = 15) or a vehicle (n = 15), and scalp lesional biopsies were collected at baseline and end of study for transcriptomic analysis. RNA was extracted from all biopsies and Affymetrix gene chips were used to analyse transcriptomic profiles, followed by bioinformatic analysis. Analysis of study 1 biopsies revealed more than 7000 individual probes differentially regulated in dandruff lesional skin relative to normal. Enriched Gene Ontology categories included: lipid metabolism, immune response, response to stimulus, apoptosis, cell proliferation, and epidermal development. The most striking feature of lesional skin relative to normal was the reciprocal expression of induced inflammatory genes and repressed lipid metabolism genes. Induced inflammatory genes were also enriched in dandruff uninvolved skin, suggesting the existence of predisposing factors associated with inflammation. Many genes increased in lesional skin were increased at the level of protein in stratum corneum samples (e.g. IL-1RA, S100A8, S100A9, S100A11, IL-8). Under conditions known to improve overall scalp condition, the ZPT shampoo treatment in study 2 produced a transcriptomic profile resembling that of normal scalp skin. These data provide novel insights into the nature of dandruff and the therapeutic action of potentiated ZPT-containing shampoo, and provide a basis to explore many new mechanistic questions related to these topics.

Effects of acupuncture on the gene expression profile of lung tissue from normal rats

Acupuncture has been demonstrated to be an effective treatment for various diseases. However, little attention has been paid to its physiological influences, especially on the changes in protein and mRNA levels following acupuncture treatment under normal conditions. In this study, we investigated the gene expression profile of lung tissue from acupuncture-treated normal rats and attempted to characterize the underlying mechanisms of the changes in expression. Three common acupoints, Dazhui (GV14), fengmen (BL12) and feishu (BL13) were selected for analysis, and 2 serial analyses of gene expression (SAGE) tag libraries of the lung tissues that were derived from the normal and acupuncture-treated rats were established. Bioinformatic analyses were carried out using the functional annotation tools of the database for annotation, visualization and integrated discovery (DAVID), the gene ontology (GO) Tree Machine and the Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. In total, 144 tags were differentially expressed (P<0.05), and the DAVID functional classification of genes demonstrated that the genes were divided into 6 types. Furthermore, GO Tree Machine analysis of the gene categories indicated that 10 enriched GO categories had become enriched after acupuncture, and that 15 KEGG pathways matched the differentially expressed tags of the 2 SAGE libraries. Our results show that the essential effects of acupuncture on normal rats include the regulation of macromolecular biosynthesis, transportation and metabolism. Cellular biosynthesis and cellular lipid metabolism are the common biological processes that occur in response to acupuncture under normal and morbid conditions, which may be the general physiological effects of acupuncture.

DEFOG: discrete enrichment of functionally organized genes

High-throughput biological experiments commonly result in a list of genes or proteins of interest. In order to understand the observed changes of the genes and to generate new hypotheses, one needs to understand the functions and roles of the genes and how those functions relate to the experimental conditions. Typically, statistical tests are performed in order to detect enriched Gene Ontology categories or pathways, i.e. the categories are observed in the genes of interest more often than is expected by chance. Depending on the number of genes and the complexity and quantity of functions in which they are involved, such an analysis can easily result in hundreds of enriched terms. To this end we developed DEFOG, a web-based application that facilitates the functional analysis of gene sets by hierarchically organizing the genes into functionally related modules. Our computational pipeline utilizes three powerful tools to achieve this goal: (1) GeneMANIA creates a functional consensus network of the genes of interest based on gene-list-specific data fusion of hundreds of genomic networks from publicly available sources; (2) Transitivity Clustering organizes those genes into a clear hierarchy of functionally related groups, and (3) Ontologizer performs a Gene Ontology enrichment analysis on the resulting gene clusters. DEFOG integrates this computational pipeline within an easy-to-use web interface, thus allowing for a novel visual analysis of gene sets that aids in the discovery of potentially important biological mechanisms and facilitates the creation of new hypotheses. DEFOG is available at http://www.mooneygroup.org/defog.

Differentiation-Promoting Culture of Competent and Noncompetent Keratinocytes Identifies Biomarkers for Head and Neck Cancer

Aberrant contact-inhibited proliferation and differentiation induction couple with tumor severity, albeit with an imprecise association with prognosis. Assessment of contact inhibition and differentiation-promoting culture in this study of normal and immortalized oral keratinocytes (NOK and SVpgC2a, respectively) demonstrated elevated cloning ability and saturation density in the immortalized versus normal state, including consistent absence of differentiated morphological features. Transcriptomic analysis implicated 48 gene ontology categories, 8 molecular networks, and 10 key regulator genes in confluency-induced differentiation of NOK, all of which remained nonregulated in SVpgC2a. The SVpgC2a versus NOK transcriptome enriched 52 gene ontology categories altogether, 18 molecular networks, and 39 key regulator genes, several of which were associated with epithelial-mesenchymal transition. Assessment of the previously described gene sets relative to training data sets of head and neck squamous cell carcinoma samples, one including data on tumor differentiation and patient outcome and one present in the Human Gene Expression Map, identified four genes with association to poor survival (COX7A1, MFAP5, MPDU1, and POLD1). This gene set predicted poor outcome in an independent data set of 71 head and neck squamous cell carcinomas. The present study defines, for the first time to our knowledge, the broad gene spectrum that couples to induction, and loss, of oral keratinocyte differentiation. Bioinformatics assessments of the results relative to clinical data generated novel differentiation-related tumor biomarkers relevant to patient outcome.

GO Trimming: Systematically reducing redundancy in large Gene Ontology datasets.

BackgroundThe increased accessibility of gene expression tools has enabled a wide variety of experiments utilizing transcriptomic analyses. As these tools increase in prevalence, the need for improved standardization in processing and presentation of data increases, as does the need to guard against interpretation bias. Gene Ontology (GO) analysis is a powerful method of interpreting and summarizing biological functions. However, while there are many tools available to investigate GO enrichment, there remains a need for methods that directly remove redundant terms from enriched GO lists that often provide little, if any, additional information.FindingsHere we present a simple yet novel method called GO Trimming that utilizes an algorithm designed to reduce redundancy in lists of enriched GO categories. Depending on the needs of the user, this method can be performed with variable stringency. In the example presented here, an initial list of 90 terms was reduced to 54, eliminating 36 largely redundant terms. We also compare this method to existing methods and find that GO Trimming, while simple, performs well to eliminate redundant terms in a large dataset throughout the depth of the GO hierarchy.ConclusionsThe GO Trimming method provides an alternative to other procedures, some of which involve removing large numbers of terms prior to enrichment analysis. This method should free up the researcher from analyzing overly large, redundant lists, and instead enable the concise presentation of manageable, informative GO lists. The implementation of this tool is freely available at: http://lucy.ceh.uvic.ca/go_trimming/cbr_go_trimming.py