Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia

BACKGROUND: It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS: The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS: The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS: This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Novel markers of the human follicle-associated epithelium identified by genomic profiling and microdissection.

BACKGROUND & AIMS: Regulation of gene expression in the follicle-associated epithelium (FAE) over Peyer's patches is largely unknown. CCL20, a chemokine that recruits immature dendritic cells, is one of the few FAE-specific markers described so far. Lymphotoxin beta (LTalpha1beta2) expressed on the membrane of immune cells triggers CCL20 expression in enterocytes. In this study, we measured expression profiles of LTalpha1beta2-treated intestinal epithelial cells and selected CCL20 -coregulated genes to identify new FAE markers. METHODS: Genomic profiles of T84 and Caco-2 cell lines treated with either LTalpha1beta2, flagellin, or tumor necrosis factor alpha were measured using the Affymetrix GeneChip U133A. Clustering analysis was used to select CCL20 -coregulated genes, and laser dissection microscopy and real-time polymerase chain reaction on human biopsy specimens was used to assess the expression of the selected markers. RESULTS: Applying a 2-way analysis of variance, we identified regulated genes upon the different treatments. A subset of genes involved in inflammation and related to the nuclear factor kappaB pathway was coregulated with CCL20 . Among these genes, the antiapoptotic factor TNFAIP3 was highly expressed in the FAE. CCL23 , which was not coregulated in vitro with CCL20 , was also specifically expressed in the FAE. CONCLUSIONS: We have identified 2 novel human FAE specifically expressed genes. Most of the CCL20 -coregulated genes did not show FAE-specific expression, suggesting that other signaling pathways are critical to modulate FAE-specific gene expression.

IQRray, a new method for Affymetrix microarray quality control, and the homologous organ conservation score, a new benchmark method for quality control metrics.

MOTIVATION: Microarray results accumulated in public repositories are widely reused in meta-analytical studies and secondary databases. The quality of the data obtained with this technology varies from experiment to experiment, and an efficient method for quality assessment is necessary to ensure their reliability. RESULTS: The lack of a good benchmark has hampered evaluation of existing methods for quality control. In this study, we propose a new independent quality metric that is based on evolutionary conservation of expression profiles. We show, using 11 large organ-specific datasets, that IQRray, a new quality metrics developed by us, exhibits the highest correlation with this reference metric, among 14 metrics tested. IQRray outperforms other methods in identification of poor quality arrays in datasets composed of arrays from many independent experiments. In contrast, the performance of methods designed for detecting outliers in a single experiment like Normalized Unscaled Standard Error and Relative Log Expression was low because of the inability of these methods to detect datasets containing only low-quality arrays and because the scores cannot be directly compared between experiments. AVAILABILITY AND IMPLEMENTATION: The R implementation of IQRray is available at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R. CONTACT: Marta.Rosikiewicz@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Bone anabolic potential of soy isoflavones and plant extracts and genomic changes during differentiation of hPOB-tert osteoblast-like cells

Summary For the nutritional management of bone health and the prevention of osteoporosis it is important to identify nutrients that positively influence the bone remodeling process at the cellular level. Soy isoflavones show promising osteoprotective effects in animals and humans but their mechanism of action in bone cells is yet poorly understood. Firstly, soy tissue cultures were characterized as a new and optimized source of isoflavones. A large variability in the isoflavone content was observed and high-producing strains (46.3 mg/g dry wt isoflavones) were identified. In the Ishikawa cells bioassay, the estrogenicity of isoflavones was confirmed to be 1000 to 10000 less than 17Mestradiol and that of the malonyl forms was shown for the first time (EC50 of 350 nM and 1880 nM for malonylgenistin and malonyldaidzin, respectively). The estrogenic activity of soya tissue culture extracts correlated to their isoflavone content. Secondly, the effects of phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway, as key mediators of bone formation, were investigated. Dietary achievable concentrations of genistein and daidzein (10vM), and statins (4xM) but not 17M estradiol (10nM), induced BMP-2 gene expression (by up to 3-fold) and inhibited the cholesterol biosynthetic pathway (by up to 50%) in the human osteoblastic cell line hP0B¬tert. In addition, several plant extracts (Cyperus rotundus, Lindera benzoin and Cnidium monnieri) induced BMP-2 gene expression but this induction was not restricted to the inhibition of the cholesterol synthesis pathway neither to the estrogenicity. Finally, the gene expression profiles during hP0B-tert differentiation induced by vitamin D and dexamethasone were analyzed with the Affymetrix human GeneChip. 1665 different genes and 98 ESTs were significantly regulated. The expression profiles of bone-related genes was largely in agreement with previously documented patterns, supporting the physiological relevance of the genomic results and the hP0B-tert cell line as a valid model of human osteoblast differentiation. The expression of alternative differentiation markers during the osteogenic treatment of hP0B-tert cells indicated that the adipocyte and myoblast differentiation pathways were repressed, confirming that these culture conditions allowed only osteoblast differentiation. The gene ontology analysis identified further sub-groups of genes that may be involved in the bone formation process. Aims of the thesis In order to define new strategies for the nutritional management of bone health and for the prevention of osteoporosis the major goal of the present work was to investigate the potential of phytonutrients to positively modulate the bone formation process at the cellular level and, in particular: 1.To select and optimise alternative plant sources containing high levels of isoflavones with estrogenic activity (Chapter 3). 2.To compare the effects of statins and phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway and to select new plant extracts with a bone anabolic potential (Chapter 4). 3.To further characterize the new human periosteal cell line, hP0B-tert, as a bone- formation model, by elucidating its gene expression profile during differentiation induced by vitamin D and dexamethasone (Chapter 5).

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

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

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

The Interleukin-1 receptor antagonist is a direct target gene of PPARalpha in liver.

BACKGROUND/AIMS: The Peroxisome Proliferator-Activated Receptor (PPAR) alpha belongs to the superfamily of Nuclear Receptors and plays an important role in numerous cellular processes, including lipid metabolism. It is known that PPARalpha also has an anti-inflammatory effect, which is mainly achieved by down-regulating pro-inflammatory genes. The objective of this study was to further characterize the role of PPARalpha in inflammatory gene regulation in liver. RESULTS: According to Affymetrix micro-array analysis, the expression of various inflammatory genes in liver was decreased by treatment of mice with the synthetic PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. In contrast, expression of Interleukin-1 receptor antagonist (IL-1ra), which was acutely stimulated by LPS treatment, was induced by PPARalpha. Up-regulation of IL-1ra by LPS was lower in PPARalpha -/- mice compared to Wt mice. Transactivation and chromatin immunoprecipitation studies identified IL-1ra as a direct positive target gene of PPARalpha with a functional PPRE present in the promoter. Up-regulation of IL-1ra by PPARalpha was conserved in human HepG2 hepatoma cells and the human monocyte/macrophage THP-1 cell line. CONCLUSIONS: In addition to down-regulating expression of pro-inflammatory genes, PPARalpha suppresses the inflammatory response by direct up-regulation of genes with anti-inflammatory properties.

Genome-wide association study of recurrent major depressive disorder in two European case-control cohorts.

Major depressive disorder (MDD) is a highly prevalent disorder with substantial heritability. Heritability has been shown to be substantial and higher in the variant of MDD characterized by recurrent episodes of depression. Genetic studies have thus far failed to identify clear and consistent evidence of genetic risk factors for MDD. We conducted a genome-wide association study (GWAS) in two independent datasets. The first GWAS was performed on 1022 recurrent MDD patients and 1000 controls genotyped on the Illumina 550 platform. The second was conducted on 492 recurrent MDD patients and 1052 controls selected from a population-based collection, genotyped on the Affymetrix 5.0 platform. Neither GWAS identified any SNP that achieved GWAS significance. We obtained imputed genotypes at the Illumina loci for the individuals genotyped on the Affymetrix platform, and performed a meta-analysis of the two GWASs for this common set of approximately half a million SNPs. The meta-analysis did not yield genome-wide significant results either. The results from our study suggest that SNPs with substantial odds ratio are unlikely to exist for MDD, at least in our datasets and among the relatively common SNPs genotyped or tagged by the half-million-loci arrays. Meta-analysis of larger datasets is warranted to identify SNPs with smaller effects or with rarer allele frequencies that contribute to the risk of MDD.

Uncovering hidden duplicated content in public transcriptomics data.

As part of the development of the database Bgee (a dataBase for Gene Expression Evolution), we annotate and analyse expression data from different types and different sources, notably Affymetrix data from GEO and ArrayExpress, and RNA-Seq data from SRA. During our quality control procedure, we have identified duplicated content in GEO and ArrayExpress, affecting ∼14% of our data: fully or partially duplicated experiments from independent data submissions, Affymetrix chips reused in several experiments, or reused within an experiment. We present here the procedure that we have established to filter such duplicates from Affymetrix data, and our procedure to identify future potential duplicates in RNA-Seq data. Database URL: http://bgee.unil.ch/

MIMAS 3.0 is a Multiomics Information Management and Annotation System.

BACKGROUND: DNA sequence integrity, mRNA concentrations and protein-DNA interactions have been subject to genome-wide analyses based on microarrays with ever increasing efficiency and reliability over the past fifteen years. However, very recently novel technologies for Ultra High-Throughput DNA Sequencing (UHTS) have been harnessed to study these phenomena with unprecedented precision. As a consequence, the extensive bioinformatics environment available for array data management, analysis, interpretation and publication must be extended to include these novel sequencing data types. DESCRIPTION: MIMAS was originally conceived as a simple, convenient and local Microarray Information Management and Annotation System focused on GeneChips for expression profiling studies. MIMAS 3.0 enables users to manage data from high-density oligonucleotide SNP Chips, expression arrays (both 3'UTR and tiling) and promoter arrays, BeadArrays as well as UHTS data using MIAME-compliant standardized vocabulary. Importantly, researchers can export data in MAGE-TAB format and upload them to the EBI's ArrayExpress certified data repository using a one-step procedure. CONCLUSION: We have vastly extended the capability of the system such that it processes the data output of six types of GeneChips (Affymetrix), two different BeadArrays for mRNA and miRNA (Illumina) and the Genome Analyzer (a popular Ultra-High Throughput DNA Sequencer, Illumina), without compromising on its flexibility and user-friendliness. MIMAS, appropriately renamed into Multiomics Information Management and Annotation System, is currently used by scientists working in approximately 50 academic laboratories and genomics platforms in Switzerland and France. MIMAS 3.0 is freely available via http://multiomics.sourceforge.net/.

Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization.

Glomalean fungi induce and colonize symbiotic tissue called arbuscular mycorrhiza on the roots of most land plants. Other fungi also colonize plants but cause disease not symbiosis. Whole-transcriptome analysis using a custom-designed Affymetrix Gene-Chip and confirmation with real-time RT-PCR revealed 224 genes affected during arbuscular mycorrhizal symbiosis. We compared these transcription profiles with those from rice roots that were colonized by pathogens (Magnaporthe grisea and Fusarium moniliforme). Over 40% of genes showed differential regulation caused by both the symbiotic and at least one of the pathogenic interactions. A set of genes was similarly expressed in all three associations, revealing a conserved response to fungal colonization. The responses that were shared between pathogen and symbiont infection may play a role in compatibility. Likewise, the responses that are different may cause disease. Some of the genes that respond to mycorrhizal colonization may be involved in the uptake of phosphate. Indeed, phosphate addition mimicked the effect of mycorrhiza on 8% of the tested genes. We found that 34% of the mycorrhiza-associated rice genes were also associated with mycorrhiza in dicots, revealing a conserved pattern of response between the two angiosperm classes.

Genome-wide linkage and association analyses to identify genes influencing adiponectin levels: the GEMS Study.

Adiponectin has a variety of metabolic effects on obesity, insulin sensitivity, and atherosclerosis. To identify genes influencing variation in plasma adiponectin levels, we performed genome-wide linkage and association scans of adiponectin in two cohorts of subjects recruited in the Genetic Epidemiology of Metabolic Syndrome Study. The genome-wide linkage scan was conducted in families of Turkish and southern European (TSE, n = 789) and Northern and Western European (NWE, N = 2,280) origin. A whole genome association (WGA) analysis (500K Affymetrix platform) was carried out in a set of unrelated NWE subjects consisting of approximately 1,000 subjects with dyslipidemia and 1,000 overweight subjects with normal lipids. Peak evidence for linkage occurred at chromosome 8p23 in NWE subjects (lod = 3.10) and at chromosome 3q28 near ADIPOQ, the adiponectin structural gene, in TSE subjects (lod = 1.70). In the WGA analysis, the single-nucleotide polymorphisms (SNPs) most strongly associated with adiponectin were rs3774261 and rs6773957 (P < 10(-7)). These two SNPs were in high linkage disequilibrium (r(2) = 0.98) and located within ADIPOQ. Interestingly, our fourth strongest region of association (P < 2 x 10(-5)) was to an SNP within CDH13, whose protein product is a newly identified receptor for high-molecular-weight species of adiponectin. Through WGA analysis, we confirmed previous studies showing SNPs within ADIPOQ to be strongly associated with variation in adiponectin levels and further observed these to have the strongest effects on adiponectin levels throughout the genome. We additionally identified a second gene (CDH13) possibly influencing variation in adiponectin levels. The impact of these SNPs on health and disease has yet to be determined.

A genome-wide association search for type 2 diabetes genes in African Americans.

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Pregnancy-induced chromatin remodeling in the breast of postmenopausal women.

Early pregnancy and multiparity are known to reduce the risk of women to develop breast cancer at menopause. Based on the knowledge that the differentiation of the breast induced by the hormones of pregnancy plays a major role in this protection, this work was performed with the purpose of identifying what differentiation-associated molecular changes persist in the breast until menopause. Core needle biopsies (CNB) obtained from the breast of 42 nulliparous (NP) and 71 parous (P) postmenopausal women were analyzed in morphology, immunocytochemistry and gene expression. Whereas in the NP breast, nuclei of epithelial cells were large and euchromatic, in the P breast they were small and hyperchromatic, showing strong methylation of histone 3 at lysine 9 and 27. Transcriptomic analysis performed using Affymetrix HG_U133 oligonucleotide arrays revealed that in CNB of the P breast, there were 267 upregulated probesets that comprised genes controlling chromatin organization, transcription regulation, splicing machinery, mRNA processing and noncoding elements including XIST. We concluded that the differentiation process induced by pregnancy is centered in chromatin remodeling and in the mRNA processing reactome, both of which emerge as important regulatory pathways. These are indicative of a safeguard step that maintains the fidelity of the transcription process, becoming the ultimate mechanism mediating the protection of the breast conferred by full-term pregnancy.