Helen Lockstone

IRF5 promotes inflammatory macrophage polarization and TH1-TH17 responses.

Polymorphisms in the gene encoding the transcription factor IRF5 that lead to higher mRNA expression are associated with many autoimmune diseases. Here we show that IRF5 expression in macrophages was reversibly induced by inflammatory stimuli and contributed to the plasticity of macrophage polarization. High expression of IRF5 was characteristic of M1 macrophages, in which it directly activated transcription of the genes encoding interleukin 12 subunit p40 (IL-12p40), IL-12p35 and IL-23p19 and repressed the gene encoding IL-10. Consequently, those macrophages set up the environment for a potent T helper type 1 (TH1)-TH17 response. Global gene expression analysis demonstrated that exogenous IRF5 upregulated or downregulated expression of established phenotypic markers of M1 or M2 macrophages, respectively. Our data suggest a critical role for IRF5 in M1 macrophage polarization and define a previously unknown function for IRF5 as a transcriptional repressor.

Renal Cyst Formation in Fh1-Deficient Mice Is Independent of the Hif/Phd Pathway: Roles for Fumarate in KEAP1 Succination and Nrf2 Signaling

Summary The Krebs cycle enzyme fumarate hydratase (FH) is a human tumor suppressor whose inactivation is associated with the development of leiomyomata, renal cysts, and tumors. It has been proposed that activation of hypoxia inducible factor (HIF) by fumarate-mediated inhibition of HIF prolyl hydroxylases drives oncogenesis. Using a mouse model, we provide genetic evidence that Fh1-associated cyst formation is Hif independent, as is striking upregulation of antioxidant signaling pathways revealed by gene expression profiling. Mechanistic analysis revealed that fumarate modifies cysteine residues within the Kelch-like ECH-associated protein 1 (KEAP1), abrogating its ability to repress the Nuclear factor (erythroid-derived 2)-like 2 (Nrf2)-mediated antioxidant response pathway, suggesting a role for Nrf2 dysregulation in FH-associated cysts and tumors.

Analysis of CD161 expression on human CD8+ T cells defines a distinct functional subset with tissue-homing properties.

CD8+ T lymphocytes play a key role in host defense, in particular against important persistent viruses, although the critical functional properties of such cells in tissue are not fully defined. We have previously observed that CD8+ T cells specific for tissue-localized viruses such as hepatitis C virus express high levels of the C-type lectin CD161. To explore the significance of this, we examined CD8+CD161+ T cells in healthy donors and those with hepatitis C virus and defined a population of CD8+ T cells with distinct homing and functional properties. These cells express high levels of CD161 and a pattern of molecules consistent with type 17 differentiation, including cytokines (e.g., IL-17, IL-22), transcription factors (e.g., retinoic acid-related orphan receptor γ-t, P = 6 × 10−9; RUNX2, P = 0.004), cytokine receptors (e.g., IL-23R, P = 2 × 10−7; IL-18 receptor, P = 4 × 10−6), and chemokine receptors (e.g., CCR6, P = 3 × 10−8; CXCR6, P = 3 × 10−7; CCR2, P = 4 × 10−7). CD161+CD8+ T cells were markedly enriched in tissue samples and coexpressed IL-17 with high levels of IFN-γ and/or IL-22. The levels of polyfunctional cells in tissue was most marked in those with mild disease (P = 0.0006). These data define a T cell lineage that is present already in cord blood and represents as many as one in six circulating CD8+ T cells in normal humans and a substantial fraction of tissue-infiltrating CD8+ T cells in chronic inflammation. Such cells play a role in the pathogenesis of chronic hepatitis and arthritis and potentially in other infectious and inflammatory diseases of man.

Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes

Aims/hypothesisMicroRNAs regulate a broad range of biological mechanisms. To investigate the relationship between microRNA expression and type 2 diabetes, we compared global microRNA expression in insulin target tissues from three inbred rat strains that differ in diabetes susceptibility.MethodsUsing microarrays, we measured the expression of 283 microRNAs in adipose, liver and muscle tissue from hyperglycaemic (Goto–Kakizaki), intermediate glycaemic (Wistar Kyoto) and normoglycaemic (Brown Norway) rats (n = 5 for each strain). Expression was compared across strains and validated using quantitative RT-PCR. Furthermore, microRNA expression variation in adipose tissue was investigated in 3T3-L1 adipocytes exposed to hyperglycaemic conditions.ResultsWe found 29 significantly differentiated microRNAs (padjusted < 0.05): nine in adipose tissue, 18 in liver and two in muscle. Of these, five microRNAs had expression patterns that correlated with the strain-specific glycaemic phenotype. MiR-222 (padjusted = 0.0005) and miR-27a (padjusted = 0.006) were upregulated in adipose tissue; miR-195 (padjusted = 0.006) and miR-103 (padjusted = 0.04) were upregulated in liver; and miR-10b (padjusted = 0.004) was downregulated in muscle. Exposure of 3T3-L1 adipocytes to increased glucose concentration upregulated the expression of miR-222 (p = 0.008), miR-27a (p = 0.02) and the previously reported miR-29a (p = 0.02). Predicted target genes of these differentially expressed microRNAs are involved in pathways relevant to type 2 diabetes.ConclusionThe expression patterns of miR-222, miR-27a, miR-195, miR-103 and miR-10b varied with hyperglycaemia, suggesting a role for these microRNAs in the pathophysiology of type 2 diabetes, as modelled by the Gyoto–Kakizaki rat. We observed similar patterns of expression of miR-222, miR-27a and miR-29a in adipocytes as a response to increased glucose levels, which supports our hypothesis that altered expression of microRNAs accompanies primary events related to the pathogenesis of type 2 diabetes.

Foxp2 Regulates Gene Networks Implicated in Neurite Outgrowth in the Developing Brain

Forkhead-box protein P2 is a transcription factor that has been associated with intriguing aspects of cognitive function in humans, non-human mammals, and song-learning birds. Heterozygous mutations of the human FOXP2 gene cause a monogenic speech and language disorder. Reduced functional dosage of the mouse version (Foxp2) causes deficient cortico-striatal synaptic plasticity and impairs motor-skill learning. Moreover, the songbird orthologue appears critically important for vocal learning. Across diverse vertebrate species, this well-conserved transcription factor is highly expressed in the developing and adult central nervous system. Very little is known about the mechanisms regulated by Foxp2 during brain development. We used an integrated functional genomics strategy to robustly define Foxp2-dependent pathways, both direct and indirect targets, in the embryonic brain. Specifically, we performed genome-wide in vivo ChIP–chip screens for Foxp2-binding and thereby identified a set of 264 high-confidence neural targets under strict, empirically derived significance thresholds. The findings, coupled to expression profiling and in situ hybridization of brain tissue from wild-type and mutant mouse embryos, strongly highlighted gene networks linked to neurite development. We followed up our genomics data with functional experiments, showing that Foxp2 impacts on neurite outgrowth in primary neurons and in neuronal cell models. Our data indicate that Foxp2 modulates neuronal network formation, by directly and indirectly regulating mRNAs involved in the development and plasticity of neuronal connections.

MicroRNA-125a is over-expressed in insulin target tissues in a spontaneous rat model of Type 2 Diabetes

BackgroundMicroRNAs (miRNAs) are non-coding RNA molecules involved in post-transcriptional control of gene expression of a wide number of genes, including those involved in glucose homeostasis. Type 2 diabetes (T2D) is characterized by hyperglycaemia and defects in insulin secretion and action at target tissues. We sought to establish differences in global miRNA expression in two insulin-target tissues from inbred rats of spontaneously diabetic and normoglycaemic strains.MethodsWe used a miRNA microarray platform to measure global miRNA expression in two insulin-target tissues: liver and adipose tissue from inbred rats of spontaneously diabetic (Goto-Kakizaki [GK]) and normoglycaemic (Brown-Norway [BN]) strains which are extensively used in genetic studies of T2D. MiRNA data were integrated with gene expression data from the same rats to investigate how differentially expressed miRNAs affect the expression of predicted target gene transcripts.ResultsThe expression of 170 miRNAs was measured in liver and adipose tissue of GK and BN rats. Based on a p-value for differential expression between GK and BN, the most significant change in expression was observed for miR-125a in liver (FC = 5.61, P = 0.001, Padjusted= 0.10); this overexpression was validated using quantitative RT-PCR (FC = 13.15, P = 0.0005). MiR-125a also showed over-expression in the GK vs. BN analysis within adipose tissue (FC = 1.97, P = 0.078, Padjusted= 0.99), as did the previously reported miR-29a (FC = 1.51, P = 0.05, Padjusted= 0.99). In-silico tools assessing the biological role of predicted miR-125a target genes suggest an over-representation of genes involved in the MAPK signaling pathway. Gene expression analysis identified 1308 genes with significantly different expression between GK and BN rats (Padjusted < 0.05): 233 in liver and 1075 in adipose tissue. Pathways related to glucose and lipid metabolism were significantly over-represented among these genes. Enrichment analysis suggested that differentially expressed genes in GK compared to BN included more predicted miR-125a target genes than would be expected by chance in adipose tissue (FDR = 0.006 for up-regulated genes; FDR = 0.036 for down-regulated genes) but not in liver (FDR = 0.074 for up-regulated genes; FDR = 0.248 for down-regulated genes).ConclusionMiR-125a is over-expressed in liver in hyperglycaemic GK rats relative to normoglycaemic BN rats, and our array data also suggest miR-125a is over-expressed in adipose tissue. We demonstrate the use of in-silico tools to provide the basis for further investigation of the potential role of miR-125a in T2D. In particular, the enrichment of predicted miR-125a target genes among differentially expressed genes has identified likely target genes and indicates that integrating global miRNA and mRNA expression data may give further insights into miRNA-mediated regulation of gene expression.

CD161 Defines a Transcriptional and Functional Phenotype across Distinct Human T Cell Lineages

Summary The C-type lectin CD161 is expressed by a large proportion of human T lymphocytes of all lineages, including a population known as mucosal-associated invariant T (MAIT) cells. To understand whether different T cell subsets expressing CD161 have similar properties, we examined these populations in parallel using mass cytometry and mRNA microarray approaches. The analysis identified a conserved CD161++/MAIT cell transcriptional signature enriched in CD161+CD8+ T cells, which can be extended to CD161+ CD4+ and CD161+TCRγδ+ T cells. Furthermore, this led to the identification of a shared innate-like, TCR-independent response to interleukin (IL)-12 plus IL-18 by different CD161-expressing T cell populations. This response was independent of regulation by CD161, which acted as a costimulatory molecule in the context of T cell receptor stimulation. Expression of CD161 hence identifies a transcriptional and functional phenotype, shared across human T lymphocytes and independent of both T cell receptor (TCR) expression and cell lineage.

Variability of gene expression profiles in human blood and lymphoblastoid cell lines.

BackgroundReadily accessible samples such as peripheral blood or cell lines are increasingly being used in large cohorts to characterise gene expression differences between a patient group and healthy controls. However, cell and RNA isolation procedures and the variety of cell types that make up whole blood can affect gene expression measurements. We therefore systematically investigated global gene expression profiles in peripheral blood from six individuals collected during two visits by comparing five of the following cell and RNA isolation methods: whole blood (PAXgene), peripheral blood mononuclear cells (PBMCs), lymphoblastoid cell lines (LCLs), CD19 and CD20 specific B-cell subsets.ResultsGene expression measurements were clearly discriminated by isolation method although the reproducibility was high for all methods (range ρ = 0.90-1.00). The PAXgene samples showed a decrease in the number of expressed genes (P < 1*10-16) with higher variability (P < 1*10-16) compared to the other methods. Differentially expressed probes between PAXgene and PBMCs were correlated with the number of monocytes, lymphocytes, neutrophils or erythrocytes. The correlations (ρ = 0.83; ρ = 0.79) of the expression levels of detected probes between LCLs and B-cell subsets were much lower compared to the two B-cell isolation methods (ρ = 0.98). Gene ontology analysis of detected genes showed that genes involved in inflammatory responses are enriched in B-cells CD19 and CD20 whereas genes involved in alcohol metabolic process and the cell cycle were enriched in LCLs.ConclusionGene expression profiles in blood-based samples are strongly dependent on the predominant constituent cell type(s) and RNA isolation method. It is crucial to understand the differences and variability of gene expression measurements between cell and RNA isolation procedures, and their relevance to disease processes, before application in large clinical studies.

Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex

The human major histocompatibility complex (MHC) on chromosome 6p21 is a paradigm for genomics, showing remarkable polymorphism and striking association with immune and non-immune diseases. The complex genomic landscape of the MHC, notably strong linkage disequilibrium, has made resolving causal variants very challenging. A promising approach is to investigate gene expression levels considered as tractable intermediate phenotypes in mapping complex diseases. However, how transcription varies across the MHC, notably relative to specific haplotypes, remains unknown. Here, using an original hybrid tiling and splice junction microarray that includes alternate allele probes, we draw the first high-resolution strand-specific transcription map for three common MHC haplotypes (HLA-A1-B8-Cw7-DR3, HLA-A3-B7-Cw7-DR15, and HLA-A26-B18-Cw5-DR3-DQ2) strongly associated with autoimmune diseases including type 1 diabetes, systemic lupus erythematosus, and multiple sclerosis. We find that haplotype-specific differences in gene expression are common across the MHC, affecting 96 genes (46.4%), most significantly the zing finger protein gene ZFP57. Differentially expressed probes are correlated with polymorphisms between haplotypes, consistent with cis effects that we directly demonstrate for ZFP57 in a cohort of healthy volunteers (P = 1.2 × 10(-14)). We establish that alternative splicing is significantly more frequent in the MHC than genome-wide (72.5% vs. 62.1% of genes, P ≤ 1 × 10(-4)) and shows marked haplotypic differences. We also unmask novel and abundant intergenic transcription involving 31% of transcribed blocks identified. Our study reveals that the renowned MHC polymorphism also manifests as transcript diversity, and our novel haplotype-based approach marks a new step toward identification of regulatory variants involved in the control of MHC-associated phenotypes and diseases.

Use of multivariate analysis to suggest a new molecular classification of colorectal cancer

Molecular classification of colorectal cancer (CRC) is currently based on microsatellite instability (MSI), KRAS or BRAF mutation and, occasionally, chromosomal instability (CIN). Whilst useful, these categories may not fully represent the underlying molecular subgroups. We screened 906 stage II/III CRCs from the VICTOR clinical trial for somatic mutations. Multivariate analyses (logistic regression, clustering, Bayesian networks) identified the primary molecular associations. Positive associations occurred between: CIN and TP53 mutation; MSI and BRAF mutation; and KRAS and PIK3CA mutations. Negative associations occurred between: MSI and CIN; MSI and NRAS mutation; and KRAS mutation, and each of NRAS, TP53 and BRAF mutations. Some complex relationships were elucidated: KRAS and TP53 mutations had both a direct negative association and a weaker, confounding, positive association via TP53–CIN–MSI–BRAF–KRAS. Our results suggested a new molecular classification of CRCs: (1) MSI+ and/or BRAF‐mutant; (2) CIN+ and/or TP53– mutant, with wild‐type KRAS and PIK3CA; (3) KRAS‐ and/or PIK3CA‐mutant, CIN+, TP53‐wild‐type; (4) KRAS– and/or PIK3CA‐mutant, CIN–, TP53‐wild‐type; (5) NRAS‐mutant; (6) no mutations; (7) others. As expected, group 1 cancers were mostly proximal and poorly differentiated, usually occurring in women. Unexpectedly, two different types of CIN+ CRC were found: group 2 cancers were usually distal and occurred in men, whereas group 3 showed neither of these associations but were of higher stage. CIN+ cancers have conventionally been associated with all three of these variables, because they have been tested en masse. Our classification also showed potentially improved prognostic capabilities, with group 3, and possibly group 1, independently predicting disease‐free survival. Copyright