A molecular roadmap of reprogramming somatic cells into iPS cells.

Factor-induced reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) is inefficient, complicating mechanistic studies. Here, we examined defined intermediate cell populations poised to becoming iPSCs by genome-wide analyses. We show that induced pluripotency elicits two transcriptional waves, which are driven by c-Myc/Klf4 (first wave) and Oct4/Sox2/Klf4 (second wave). Cells that become refractory to reprogramming activate the first but fail to initiate the second transcriptional wave and can be rescued by elevated expression of all four factors. The establishment of bivalent domains occurs gradually after the first wave, whereas changes in DNA methylation take place after the second wave when cells acquire stable pluripotency. This integrative analysis allowed us to identify genes that act as roadblocks during reprogramming and surface markers that further enrich for cells prone to forming iPSCs. Collectively, our data offer new mechanistic insights into the nature and sequence of molecular events inherent to cellular reprogramming.

Identifying predictive signalling networks for Vedolizumab response in ulcerative colitis.

BACKGROUND: In ulcerative colitis (UC), the molecular mechanisms that drive disease development and patient response to therapy are not well understood. A significant proportion of patients with UC fail to respond adequately to biologic therapy. Therefore, there is an unmet need for biomarkers that can predict patients' responsiveness to the available UC therapies as well as ascertain the most effective individualised therapy. Our study focused on identifying predictive signalling pathways that predict anti-integrin therapy response in patients with UC. METHODS: We retrieved and pre-processed two publicly accessible gene expression datasets (GSE73661 and GSE72819) of UC patients treated with anti-integrin therapies: (1) 12 non-IBD controls and 41 UC patients treated with Vedolizumab therapy, and (2) 70 samples with 58 non-responder and 12 responder UC patient samples treated with Etrolizumab therapy without non-IBD controls. We used a diffusion-based signalling model which is mainly focused on the T-cell receptor signalling network. The diffusion model uses network connectivity between receptors and transcription factors. RESULTS: The network diffusion scores were able to separate VDZ responder and non-responder patients before treatment better than the original gene expression. On both anti-integrin treatment datasets, the diffusion model demonstrated high predictive performance for discriminating responders from non-responders in comparison with 'nnet'. We have found 48 receptor-TF pairs identified as the best predictors for VDZ therapy response with AUC ≥ 0.76. Among these receptor-TF predictors pairs, FFAR2-NRF1, FFAR2-RELB, FFAR2-EGR1, and FFAR2-NFKB1 are the top best predictors. For Etrolizumab, we have identified 40 best receptor-TF pairs and CD40-NFKB2 as the best predictor receptor-TF pair (AUC = 0.72). We also identified subnetworks that highlight the network interactions, connecting receptors and transcription factors involved in cytokine and fatty acid signalling. The findings suggest that anti-integrin therapy responses in cytokine and fatty acid signalling can stratify UC patient subgroups. CONCLUSIONS: We identified signalling pathways that may predict the efficacy of anti-integrin therapy in UC patients and personalised therapy alternatives. Our results may lead to the advancement of a promising clinical decision-making tool for the stratification of UC patients.

Mucosal gene transcription of ulcerative colitis in endoscopic remission.

Aim/Objective: Ulcerative colitis (UC) is a chronic inflammatory bowel disease. In UC, a wide range of criteria are used for disease remission, with few studies investigating the differences between disease remission and normal control groups. This paper compares known inflammatory and healing mediators in the mucosa of UC in clinical remission and normal controls, in order to better describe the remission state.Method: Mucosal biopsies from 72 study participants (48 UC and 24 normal controls) were included from the Advanced Study of Inflammatory Bowel Disease (ASIB Study), Arctic University of Norway, Norway. Clinical remission was defined as Mayo clinical score ≤ 2, with endoscopic subscores of ≤ 1. Targeted gene transcription analyses were performed using hydrolysis probes and SYBR-green.Results: Among the mucosal transcripts examined, 10 genes were regulated in remission versus normal controls, 8 upregulated pro-inflammatory transcripts (IL1B, IL33, TNF, TRAF1, CLDN2, STAT1, STAT3 and IL13Ra2) and 2 downregulated (pro-inflammatory TBX21 and anti-inflammatory TGFB1). In total, 14 transcripts were regulated between the investigated groups. Several master transcription factors for T-cell development were upregulated in patients with Mayo endoscopic score of 1 in comparison to 0.Conclusions: The mucosa of UC in clinical and endoscopic remission differs from normal mucosa, suggesting a remaining dysregulation of inflammatory and wound healing mechanisms.

dREAM co-operates with insulator-binding proteins and regulates expression at divergently paired genes.

dREAM complexes represent the predominant form of E2F/RBF repressor complexes in Drosophila. dREAM associates with thousands of sites in the fly genome but its mechanism of action is unknown. To understand the genomic context in which dREAM acts we examined the distribution and localization of Drosophila E2F and dREAM proteins. Here we report a striking and unexpected overlap between dE2F2/dREAM sites and binding sites for the insulator-binding proteins CP190 and Beaf-32. Genetic assays show that these components functionally co-operate and chromatin immunoprecipitation experiments on mutant animals demonstrate that dE2F2 is important for association of CP190 with chromatin. dE2F2/dREAM binding sites are enriched at divergently transcribed genes, and the majority of genes upregulated by dE2F2 depletion represent the repressed half of a differentially expressed, divergently transcribed pair of genes. Analysis of mutant animals confirms that dREAM and CP190 are similarly required for transcriptional integrity at these gene pairs and suggest that dREAM functions in concert with CP190 to establish boundaries between repressed/activated genes. Consistent with the idea that dREAM co-operates with insulator-binding proteins, genomic regions bound by dREAM possess enhancer-blocking activity that depends on multiple dREAM components. These findings suggest that dREAM functions in the organization of transcriptional domains.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication.

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.

A shared role for RBF1 and dCAP-D3 in the regulation of transcription with consequences for innate immunity.

Previously, we discovered a conserved interaction between RB proteins and the Condensin II protein CAP-D3 that is important for ensuring uniform chromatin condensation during mitotic prophase. The Drosophila melanogaster homologs RBF1 and dCAP-D3 co-localize on non-dividing polytene chromatin, suggesting the existence of a shared, non-mitotic role for these two proteins. Here, we show that the absence of RBF1 and dCAP-D3 alters the expression of many of the same genes in larvae and adult flies. Strikingly, most of the genes affected by the loss of RBF1 and dCAP-D3 are not classic cell cycle genes but are developmentally regulated genes with tissue-specific functions and these genes tend to be located in gene clusters. Our data reveal that RBF1 and dCAP-D3 are needed in fat body cells to activate transcription of clusters of antimicrobial peptide (AMP) genes. AMPs are important for innate immunity, and loss of either dCAP-D3 or RBF1 regulation results in a decrease in the ability to clear bacteria. Interestingly, in the adult fat body, RBF1 and dCAP-D3 bind to regions flanking an AMP gene cluster both prior to and following bacterial infection. These results describe a novel, non-mitotic role for the RBF1 and dCAP-D3 proteins in activation of the Drosophila immune system and suggest dCAP-D3 has an important role at specific subsets of RBF1-dependent genes.

The duration of gastrin treatment affects global gene expression and molecular responses involved in ER stress and anti-apoptosis.

BACKGROUND: How cells decipher the duration of an external signal into different transcriptional outcomes is poorly understood. The hormone gastrin can promote a variety of cellular responses including proliferation, differentiation, migration and anti-apoptosis. While gastrin in normal concentrations has important physiological functions in the gastrointestine, prolonged high levels of gastrin (hypergastrinemia) is related to pathophysiological processes. RESULTS: We have used genome-wide microarray time series analysis and molecular studies to identify genes that are affected by the duration of gastrin treatment in adenocarcinoma cells. Among 403 genes differentially regulated in transiently (gastrin removed after 1 h) versus sustained (gastrin present for 14 h) treated cells, 259 genes upregulated by sustained gastrin treatment compared to untreated controls were expressed at lower levels in the transient mode. The difference was subtle for early genes like Junb and c-Fos, but substantial for delayed and late genes. Inhibition of protein synthesis by cycloheximide was used to distinguish between primary and secondary gastrin regulated genes. The majority of gastrin upregulated genes lower expressed in transiently treated cells were primary genes induced independently of de novo protein synthesis. This indicates that the duration effect of gastrin treatment is mainly mediated via post-translational signalling events, while a smaller fraction of the differentially expressed genes are regulated downstream of primary transcriptional events. Indeed, sustained gastrin treatment specifically induced prolonged ERK1/2 activation and elevated levels of the AP-1 subunit protein JUNB. Enrichment analyses of the differentially expressed genes suggested that endoplasmic reticulum (ER) stress and survival is affected by the duration of gastrin treatment. Sustained treatment exerted an anti-apoptotic effect on serum starvation-induced apoptosis via a PKC-dependent mechanism. In accordance with this, only sustained treatment induced anti-apoptotic genes like Clu, Selm and Mcl1, while the pro-apoptotic gene Casp2 was more highly expressed in transiently treated cells. Knockdown studies showed that JUNB is involved in sustained gastrin induced expression of the UPR/ER stress related genes Atf4, Herpud1 and Chac1. CONCLUSION: The duration of gastrin treatment affects both intracellular signalling mechanisms and gene expression, and ERK1/2 and AP-1 seem to play a role in converting different durations of gastrin treatment into distinct cellular responses.

Associations of breast cancer related exposures and gene expression profiles in normal breast tissue-The Norwegian Women and Cancer normal breast tissue study.

BACKGROUND: Normal breast tissue is utilized in tissue-based studies of breast carcinogenesis. While gene expression in breast tumor tissue is well explored, our knowledge of transcriptomic signatures in normal breast tissue is still incomplete. The aim of this study was to investigate variability of gene expression in a large sample of normal breast tissue biopsies, according to breast cancer related exposures (obesity, smoking, alcohol, hormone therapy, and parity). METHODS: We analyzed gene expression profiles from 311 normal breast tissue biopsies from cancer-free, post-menopausal women, using Illumina bead chip arrays. Principal component analysis and K-means clustering was used for initial analysis of the dataset. The association of exposures and covariates with gene expression was determined using linear models for microarrays. RESULTS: Heterogeneity of the breast tissue and cell composition had the strongest influence on gene expression profiles. After adjusting for cell composition, obesity, smoking, and alcohol showed the highest numbers of associated genes and pathways, whereas hormone therapy and parity were associated with negligible gene expression differences. CONCLUSION: Our results provide insight into associations between major exposures and gene expression profiles and provide an informative baseline for improved understanding of exposure-related molecular events in normal breast tissue of cancer-free, post-menopausal women.

Endogenous IL-1 receptor antagonist restricts healthy and malignant myeloproliferation.

Here we explored the role of interleukin-1ß (IL-1ß) repressor cytokine, IL-1 receptor antagonist (IL-1rn), in both healthy and abnormal hematopoiesis. Low IL-1RN is frequent in acute myeloid leukemia (AML) patients and represents a prognostic marker of reduced survival. Treatments with IL-1RN and the IL-1ß monoclonal antibody canakinumab reduce the expansion of leukemic cells, including CD34+ progenitors, in AML xenografts. In vivo deletion of IL-1rn induces hematopoietic stem cell (HSC) differentiation into the myeloid lineage and hampers B cell development via transcriptional activation of myeloid differentiation pathways dependent on NFκB. Low IL-1rn is present in an experimental model of pre-leukemic myelopoiesis, and IL-1rn deletion promotes myeloproliferation, which relies on the bone marrow hematopoietic and stromal compartments. Conversely, IL-1rn protects against pre-leukemic myelopoiesis. Our data reveal that HSC differentiation is controlled by balanced IL-1ß/IL-1rn levels under steady-state, and that loss of repression of IL-1ß signaling may underlie pre-leukemic lesion and AML progression.

Lung cancer and DNA repair genes: multilevel association analysis from the International Lung Cancer Consortium.

Lung cancer (LC) is the leading cause of cancer-related death worldwide and tobacco smoking is the major associated risk factor. DNA repair is an important process, maintaining genome integrity and polymorphisms in DNA repair genes may contribute to susceptibility to LC. To explore the role of DNA repair genes in LC, we conducted a multilevel association study with 1655 single nucleotide polymorphisms (SNPs) in 211 DNA repair genes using 6911 individuals pooled from four genome-wide case-control studies. Single SNP association corroborates previous reports of association with rs3131379, located on the gene MSH5 (P = 3.57 × 10-5) and returns a similar risk estimate. The effect of this SNP is modulated by histological subtype. On the log-additive scale, the odds ratio per allele is 1.04 (0.84-1.30) for adenocarcinomas, 1.52 (1.28-1.80) for squamous cell carcinomas and 1.31 (1.09-1.57) for other histologies (heterogeneity test: P = 9.1 × 10(-)(3)). Gene-based association analysis identifies three repair genes associated with LC (P < 0.01): UBE2N, structural maintenance of chromosomes 1L2 and POLB. Two additional genes (RAD52 and POLN) are borderline significant. Pathway-based association analysis identifies five repair pathways associated with LC (P < 0.01): chromatin structure, DNA polymerases, homologous recombination, genes involved in human diseases with sensitivity to DNA-damaging agents and Rad6 pathway and ubiquitination. This first international pooled analysis of a large dataset unravels the role of specific DNA repair pathways in LC and highlights the importance of accounting for gene and pathway effects when studying LC.

Identifying anti-TNF response biomarkers in ulcerative colitis using a diffusion-based signalling model.

Motivation: Resistance to anti-TNF therapy in subgroups of ulcerative colitis (UC) patients is a major challenge and incurs significant treatment costs. Identification of patients at risk of nonresponse to anti-TNF is of major clinical importance. To date, no quantitative computational framework exists to develop a complex biomarker for the prognosis of UC treatment. Modelling patient-wise receptor to transcription factor (TF) network connectivity may enable personalized treatment. Results: We present an approach for quantitative diffusion analysis between receptors and TFs using gene expression data. Key TFs were identified using pandaR. Network connectivities between immune-specific receptor-TF pairs were quantified using network diffusion in UC patients and controls. The patient-specific network could be considered a complex biomarker that separates anti-TNF treatment-resistant and responder patients both in the gene expression dataset used for model development and separate independent test datasets. The model was further validated in rheumatoid arthritis where it successfully discriminated resistant and responder patients to tocilizumab treatment. Our model may contribute to prognostic biomarkers that may identify treatment-resistant and responder subpopulations of UC patients. Availability and implementation: Software is available at https://github.com/Amy3100/receptor2tfDiffusion. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

DNA hypo-methylation facilitates anti-inflammatory responses in severe ulcerative colitis.

Severe ulcerative colitis (UC) is a potentially life-threatening disease with a potential colorectal cancer (CRC) risk. The aim of this study was to explore the relationship between transcriptomic and genome-wide DNA methylation profiles in a well-stratified, treatment-naïve severe UC patient population in order to define specific epigenetic changes that could be responsible for the grade of disease severity. Mucosal biopsies from treatment-naïve severe UC patients (n = 8), treatment-naïve mild UC (n = 8), and healthy controls (n = 8) underwent both whole transcriptome RNA-Seq and genome-wide DNA bisulfite- sequencing, and principal component analysis (PCA), cell deconvolutions and diverse statistical methods were applied to obtain a dataset of significantly differentially expressed genes (DEGs) with correlation to DNA methylation for severe UC. DNA hypo-methylation correlated with approximately 80% of all DEGs in severe UC when compared to mild UC. Enriched pathways of annotated hypo-methylated genes revealed neutrophil degranulation, and immuno-regulatory interactions of the lymphoid system. Specifically, hypo-methylated anti-inflammatory genes found for severe UC were IL10, SIGLEC5, CD86, CLMP and members of inflammasomes NLRP3 and NLRC4. Hypo-methylation of anti-inflammatory genes during severe UC implies an interplay between the epithelium and lamina propria in order to mitigate inflammation in the gut. The specifically DNA hypo-methylated genes found for severe UC can potentially be useful biomarkers for determining disease severity and in the development of new targeted treatment strategies for severe UC patients.

A rat model of intravesical delivery of small interfering RNA for studying urinary carcinoma.

PURPOSE: siRNA has been used successfully in loss-of-function studies in vitro, but neither in vivo nor in clinical applications. The aims of the present study were (1) to establish rat models for in vivo delivery of siRNA to bladder cancer, and (2) to identify potential targets for siRNA. METHODS: The rat models of human urinary carcinoma and rat urinary carcinoma cell line (AY-27) were induced by tobacco-related chemical carcinogens, either N-[4-(5-nitro-2-furyl)-2-thiazolyl]formamide (FANFT) or N-butyl-N-(4-hydroxybutyl) nitrosamine (BBN). A syngeneic orthotopic bladder cancer model was established by inoculation of AY-27 cells. A fluorescence-labelled negative control siRNA with cationic and neutral liposomes was tested both in vitro (AY-27 cells) and in vivo. RESULTS: siRNA was highly accumulated in the cancer cells as early as 12 h and remained at least for 24 h after a single dose in vivo. Numerous CD3+ T cells appeared mainly in the periphery area of the tumour. Bioinformatics analysis revealed a list of concordantly highly expressed genes, possible siRNA targets, in the animal models as well as human urinary carcinoma. Literature search on siRNA and bladder cancer provided a list of genes used as siRNA targets. CONCLUSION: The methodology and data presented in the present study provide a number of opportunities for basic research on urinary carcinogenesis and for translational research on evaluation of siRNA therapeutic strategies for urinary carcinoma in the native organ, where hormonal, neural and immunological processes more closely resemble the clinical situation.

Fibrosis Mediators in the Colonic Mucosa of Acute and Healed Ulcerative Colitis.

OBJECTIVES: A healed intestinal mucosa is the aim of therapy in acute ulcerative colitis (UC). Disruption of mucosal wound healing may lead to severe complications including intestinal fibrosis. This study examined mucosal gene expression in the healing process of acute UC with a special focus on known mediators of fibrosis. METHODS: Endoscopic biopsies from patients with acute, moderate to severe UC were analyzed with a quantitative polymerase chain reaction array for 84 genes involved in fibrosis pathways. All patients were treated with infliximab (anti- tumor necrosis factor). Biopsies were taken before therapy and when disease remission was reached, defined as a Mayo score of ≤2, with an endoscopic subscore of 0 or 1. A healthy control group was included. Immunostaining of matrix metallopeptidase 9 and smooth muscle actin was performed. RESULTS: Mucosal biopsies from acute UC (n = 28), remission UC (n = 28), and healthy controls (n = 13) were analyzed. Fibrosis and extracellular matrix-associated genes were upregulated in the endoscopically healed UC mucosa vs controls, with collagen type III alpha 1 chain, actin alpha 2, lysyl oxidase, TIMP metallopeptidase inhibitor 3, and caveolin 1 uniquely showing no overlap with acute disease. Pro- and antifibrotic mediators (interleukin [IL]13 receptor subunit alpha 2, IL1B, IL10, tumor necrosis factor, snail family transcriptional repressor 1, and C-C motif chemokine ligand 2) were upregulated in both acute and healed UC compared with controls. An attenuated pattern of the canonical transforming growth factor beta (TGFB) pathway was observed in acute UC and to a lesser extent in the healed mucosa, except for TGFB2, which was enhanced. DISCUSSION: The endoscopically healed mucosa of UC showed a persisting dysregulation of fibrosis-associated mediators compared with controls, including extracellular matrix remodeling, profibrotic cytokines, and TGFB signaling pathways.

Phylogenetic backgrounds and virulence profiles of atypical enteropathogenic Escherichia coli strains from a case-control study using multilocus sequence typing and DNA microarray analysis.

Atypical enteropathogenetic Escherichia coli (EPEC) strains are frequently detected in children with diarrhea but are also a common finding in healthy children. The aim of this study was to compare the phylogenetic ancestry and virulence characteristics of atypical (eae positive, stx and bfpA negative) EPEC strains from Norwegian children with (n = 37) or without (n = 19) diarrhea and to search for an association between phylogenetic ancestry and diarrhea. The strains were classified in phylogenetic groups by phylogenetic marker genes and in sequence types (STs) by multilocus sequence typing. Phylogenetic ancestry was compared to virulence characteristics based on DNA microarray analysis. Serotyping and pulsed-field gel electrophoresis (PFGE) were also performed. All four phylogenetic groups, 26 different STs, and 20 different clonal groups were represented among the 56 atypical EPEC strains. The strains were separated into three clusters by overall virulence gene profile; one large cluster with A, B1, and D strains and two clusters with group B2 strains. There was considerable heterogeneity in the PFGE profiles and serotypes, and almost half of the strains were O nontypeable. The efa1/lifA gene, previously shown to be statistically linked with diarrhea in this strain collection (J. E. Afset et al., J. Clin. Microbiol. 44:3703-3711, 2006), was present in 8 of 26 STs. The two phylogenetic groups B1 and D were weakly associated with diarrhea (P = 0.06 and P = 0.09, respectively). In contrast, group B2 was isolated most frequently from healthy controls (P = 0.05). In conclusion, the atypical EPEC strains were heterogeneous both phylogenetically and by virulence profile. Phylogenetic ancestry was less useful as a predictor of diarrhea than were specific virulence genes.

Genome-wide DNA Methylation in Treatment-naïve Ulcerative Colitis.

BACKGROUND AND AIMS: The aim of this study was to investigate the genome-wide DNA methylation status in treatment-naïve ulcerative colitis [UC], and to explore the relationship between DNA methylation patterns and gene expression levels in tissue biopsies from a well-stratified treatment-naïve UC patient group. METHODS: Mucosal biopsies from treatment-naïve patients [n = 10], and a healthy control group [n = 11] underwent genome-wide DNA bisulfite sequencing. Principal component analysis [PCA] and diverse statistical methods were applied to obtain a dataset of differentially methylated genes. DNA methylation annotation was investigated using the UCSC Genome Browser. Gene set enrichments were obtained using the Kyoto Encyclopaedia of Genes and Genomes [KEGG] and PANTHER. RESULTS: Of all significantly differentially expressed genes [DEGs], 25% correlated with DNA methylation patterns; 30% of these genes were methylated at CpG sites near their transcription start site [TSS]. Hyper-methylation was observed for genes involved in homeostasis and defence, whereas hypo-methylation was observed for genes playing a role in immune response [i.e. chemokines and interleukins]. Of the differentially DNA methylated genes, 25 were identified as inflammatory bowel disease [IBD] susceptibility genes. Four genes [DEFFA6, REG1B, BTNL3, OLFM4] showed DNA methylation in the absence of known CpG islands. CONCLUSIONS: Genome-wide DNA methylation analysis revealed distinctive functional patterns for hyper-and hypo-methylation in treatment-naïve UC. These distinct patterns could be of importance in the development and pathogenesis of UC. Further investigation of DNA methylation patterns may be useful in the development of the targeting of epigenetic processes, and may allow new treatment and target strategies for UC patients.

Transcriptomic Landscape of Treatment-Naïve Ulcerative Colitis.

BACKGROUND AND AIMS: Ulcerative colitis [UC] is a chronic inflammatory disease that effects the gastrointestinal tract and is considered one of the most prominent and common forms of inflammatory bowel disease [IBD]. This study aimed to define and describe the entire transcriptomic landscape in a well-stratified, treatment-naïve UC patient population compared with control patients by using next-generation technology, RNA-Seq. METHODS: Mucosal biopsies from treatment-naïve UC patients [n = 14], and healthy controls [n = 16] underwent RNA-Seq. Principal component analysis [PCA], cell deconvolution methods, and diverse statistical methods were applied to obtain and characterise a dataset of significantly differentially expressed genes [DEGs]. RESULTS: Analyses revealed 1480 significantly DEGs in treatment-naïve UC when compared with controls. Cell populations of monocytes, T cells, neutrophils, B cells/ lymphoid cells, and myeloid cells were increased during inflammation, whereas the fraction of epithelial cells were reduced in UC, which is reflected by the DEGs; 79 DEGs were identified as IBD susceptibility genes, and 58 DEGs were expressed in a gender-specific manner. MUC5B, REG3A, DEFA5, and IL33 might be considered as colorectal cancer [CRC] risk factors following UC in males. AQP9 together with CLDN2 may have a role regulating tissue-specific physiological properties in tight junctions in UC. An additional functional role for AQP9 in the synthesis and/or the function of mucus can be implied. CONCLUSIONS: This study reveals new potential players in UC pathogenesis in general, and provides evidence for a gender-dependent pathogenesis for UC. These results can be useful for the development of personalised treatment strategies for UC in the future.

A framework for significance analysis of gene expression data using dimension reduction methods.

BACKGROUND: The most popular methods for significance analysis on microarray data are well suited to find genes differentially expressed across predefined categories. However, identification of features that correlate with continuous dependent variables is more difficult using these methods, and long lists of significant genes returned are not easily probed for co-regulations and dependencies. Dimension reduction methods are much used in the microarray literature for classification or for obtaining low-dimensional representations of data sets. These methods have an additional interpretation strength that is often not fully exploited when expression data are analysed. In addition, significance analysis may be performed directly on the model parameters to find genes that are important for any number of categorical or continuous responses. We introduce a general scheme for analysis of expression data that combines significance testing with the interpretative advantages of the dimension reduction methods. This approach is applicable both for explorative analysis and for classification and regression problems. RESULTS: Three public data sets are analysed. One is used for classification, one contains spiked-in transcripts of known concentrations, and one represents a regression problem with several measured responses. Model-based significance analysis is performed using a modified version of Hotelling's T2-test, and a false discovery rate significance level is estimated by resampling. Our results show that underlying biological phenomena and unknown relationships in the data can be detected by a simple visual interpretation of the model parameters. It is also found that measured phenotypic responses may model the expression data more accurately than if the design-parameters are used as input. For the classification data, our method finds much the same genes as the standard methods, in addition to some extra which are shown to be biologically relevant. The list of spiked-in genes is also reproduced with high accuracy. CONCLUSION: The dimension reduction methods are versatile tools that may also be used for significance testing. Visual inspection of model components is useful for interpretation, and the methodology is the same whether the goal is classification, prediction of responses, feature selection or exploration of a data set. The presented framework is conceptually and algorithmically simple, and a Matlab toolbox (Mathworks Inc, USA) is supplemented.

Optimization of cDNA microarrays procedures using criteria that do not rely on external standards.

BACKGROUND: The measurement of gene expression using microarray technology is a complicated process in which a large number of factors can be varied. Due to the lack of standard calibration samples such as are used in traditional chemical analysis it may be a problem to evaluate whether changes done to the microarray procedure actually improve the identification of truly differentially expressed genes. The purpose of the present work is to report the optimization of several steps in the microarray process both in laboratory practices and in data processing using criteria that do not rely on external standards. RESULTS: We performed a cDNA microarry experiment including RNA from samples with high expected differential gene expression termed "high contrasts" (rat cell lines AR42J and NRK52E) compared to self-self hybridization, and optimized a pipeline to maximize the number of genes found to be differentially expressed in the "high contrasts" RNA samples by estimating the false discovery rate (FDR) using a null distribution obtained from the self-self experiment. The proposed high-contrast versus self-self method (HCSSM) requires only four microarrays per evaluation. The effects of blocking reagent dose, filtering, and background corrections methodologies were investigated. In our experiments a dose of 250 ng LNA (locked nucleic acid) dT blocker, no background correction and weight based filtering gave the largest number of differentially expressed genes. The choice of background correction method had a stronger impact on the estimated number of differentially expressed genes than the choice of filtering method. Cross platform microarray (Illumina) analysis was used to validate that the increase in the number of differentially expressed genes found by HCSSM was real. CONCLUSION: The results show that HCSSM can be a useful and simple approach to optimize microarray procedures without including external standards. Our optimizing method is highly applicable to both long oligo-probe microarrays which have become commonly used for well characterized organisms such as man, mouse and rat, as well as to cDNA microarrays which are still of importance for organisms with incomplete genome sequence information such as many bacteria, plants and fish.