MicroRNAs as tools to predict glucocorticoid response in inflammatory bowel diseases

Abstract

In spite of the introduction in therapy of highly effective biological agents, glucocorticoids (GCs) are still employed to induce remission in moderate to severe inflammatory bowel diseases (IBD), but considerable inter-individual differences in their efficacy and side effects have been reported. The effectiveness of these drugs is indeed very variable and side effects, particularly severe in pediatric patients, are common and often unpredictable: the understanding of the complex gene regulation mediated by GCs could shed light on the causes of this variability. In this context, microRNAs (miRNAs) represent a new and promising field of research. miRNAs are small non-coding RNA molecules that suppress gene expression at post-transcriptional level, and are fine-tuning regulators of diverse biological processes, including the development and function of the immune system, apoptosis, metabolism and inflammation. Emerging data have implicated the deregulated expression of certain miRNA networks in the pathogenesis of autoimmune and inflammatory diseases, such as IBD. There is a great interest in the identification of the role of miRNAs in the modulation of pharmacological response; however, the association between miRNA and GC response in patients with IBD has not yet been evaluated in a prospective clinical study. The identification of miRNAs differently expressed as a consequence of GC treatment in comparison to diagnosis, represents an important innovative approach that could be translated into clinical practice. In this review we highlight the altered regulation of proteins involved in GC molecular mechanism by miRNAs, and their potential role as molecular markers useful for predicting in advance GC response.

Key Words: Glucocorticoids, Inflammatory bowel diseases, MicroRNA, Molecular markers, Pharmacogenomics

Core tip: Studies on microRNAs (miRNAs) and pharmacogenomics represent a promising investigation topic that could increase the understanding of the pharmacology of steroids in inflammatory bowel diseases (IBDs) and possibly in other diseases. A number of studies have shown that glucocorticoids (GCs) can modify the expression profiles of different miRNAs, however, the obtained results have been highly variable, and to date it is not possible to recognize a specific miRNA pattern regulated by GCs. Moreover, existing studies employed techniques based on the use of reverse transcription quantitative polymerase chain reaction and microarrays, through the analysis and quantification of already known miRNAs. Using next generation sequencing technologies, it could be possible to detect novel, still unrecognised miRNAs, and identify new miRNA isoforms (iso-miRs) as well. This innovative approach could be a valuable tool for a better understanding of the role of miRNAs to predict steroid response in IBDs. In the future, the increased availability and the reduced costs of RNA profiling should enable the clinicians to stratify patients on specific miRNA biomarkers before starting GC treatment.

INTRODUCTION

To date, a curative pharmacological therapy for inflammatory bowel diseases (IBD) does not exist and the therapeutic approach is mainly aimed at controlling inflammation, with drugs capable of inducing and maintaining remission. Despite the introduction in therapy of highly effective biological agents, in IBD patients with moderate to severe disease glucocorticoids (GCs) are effective in inducing remission and are still considered the standard for treatment[1]. In spite of the large clinical use, the benefits of these agents are often narrowed by high inter-individual variability. Given the high incidence of suboptimal response, associated with a significant number of side effects, the identification of subjects that are most likely to respond poorly to these agents is extremely important. However, the mechanisms of this variability are scarcely understood and there is presently no means to predict the response in advance[2-5]; in this context, microRNAs (miRNAs) represent a new and promising field of research.

miRNAs are small (18-24 nucleotides) non-coding RNAs, which bind the 3’UTRs and the coding exons of their target genes and inhibit gene expression[6] either by messenger RNA (mRNA) cleavage (most common in plants) or by translational repression (most common in metazoan)[7,8]. According to the miRNA database miRBase, 1872 precursors and 2578 human mature miRNA sequences have been published (http://www.mirbase.org[9,10]) and we are only on the verge of understanding their physiological impact on gene regulation. A single miRNA can regulate a multitude of mRNAs (approximately 200), and each mRNA can be regulated by multiple miRNAs[11,12]; overall, it is predicted that protein production for at least 20% of all human genes is regulated by miRNAs[13,14].

By affecting gene regulation, miRNAs are likely to be implicated in the control of diverse biological processes, such as cellular proliferation and apoptosis[10,15-17], stem cell differentiation[15,18-20], and organ development and morphogenesis[21,22]; in addition a strong association between miRNA expression dysregulation and induction of cancer has been shown[23-26]. Moreover, miRNAs have important regulatory roles in the innate and adaptive immune system[27-29], and characteristic miRNA expression profiles have been demonstrated even in IBD[30-33].

There is a lot of interest in identifying the role of miRNAs in the modulation of drug response[34], but studies about this topic are still very limited, and the possible correlation between miRNAs expression and variability on GC response in IBD patients has not yet been examined. A better knowledge of miRNAs role could lead to their use as biomarkers for IBD, and consequently, to the development of new strategies for therapy personalization in these diseases.

This review tries to highlight the altered regulation of proteins involved in GC molecular mechanism by miRNAs in different diseases and in vitro models, and their potential role as molecular markers useful for predicting in advance GC response.

GLUCOCORTICOIDS IN INFLAMMATORY BOWEL DISEASES

GCs are effective inhibitors of cytokine secretion and T-cell activation, and are consequently largely employed in different inflammatory conditions, including IBD. Despite the introduction of novel therapies, these agents are still currently used for induction of remission in moderate to severe IBDs, however, a wide variability in response to these agents is evident and, in these diseases, GC resistance or dependence is particularly frequent. Among the adult IBD population, a prospective analysis has described the 1-year outcome in patients with Crohn’s disease (CD) treated with a first oral prednisone course (40-60 mg/d) and tapering to a maintenance dose of 10-15 mg/d[35]. Prolonged steroid response was obtained in 44% of patients, 36% of subjects were steroid dependent while 20% of subjects did not respond and were steroid resistant; a high frequency of surgery was reported within 1 mo after steroid treatment. Similar results have been obtained in a retrospective American study: immediate outcomes for CD and ulcerative colitis (UC), respectively, were complete remission in 58% and 54% of cases, partial remission in 26% and 30%, resistance in 16% of patients[36]. In paediatric IBD patients, clinical reports have shown that up to 90% of subjects has a rapid improvement of symptoms when prednisone treatment is given; however, after 1 year, only 55% of patients were still in remission and were considered steroid responsive. In around 38% of patients, steroid therapy could not be discontinued as patients experienced an increase of disease activity when the dose was reduced (steroid dependent)[37].

Demographic and/or clinical markers[36,38,39] have been evaluated and related with this variability in GC response, but results have not been consistently replicated. Genetic and epigenetic markers are likely to complement clinical and demographic predictors: phenotypes resulting from genetic changes and regulation can markedly influence drug pharmacokinetics or alter drug efficacy and/or toxicity profiles. The identification of genetic biomarkers that can be useful for classifying the disease and help to improve therapy is paramount.

MOLECULAR MECHANISM OF GC ACTION

The effects of GCs are mediated by the glucocorticoid receptor (GR)-α, a member of the nuclear receptor superfamily of ligand-dependent transcription factors[40,41]. The human GR gene is encoded on chromosome 5q31.3 and consists of nine coding exons[42]. Alternative splicing of exon 9 generates two receptor isoforms, GR-α and GR-β[43-46]. GR-β is not able to bind GCs, resides constitutively in the nucleus of cells, has a longer half-life than GR-α, and does not transactivate GC-inducible reporter genes[47]. It has been suggested[48,49] that cell specific expression and function of GR isoforms may explain the tissue and individual selective actions of GCs.

The function of GR is conditioned by chaperone and co-chaperone proteins that form a molecular heterocomplex with the GR itself[50,51], required for proper ligand binding, receptor activation and transcription: abnormalities in proteins that make up the heterocomplex may contribute to altered GC responsiveness[52,53]. Several studies have demonstrated differences in the heterocomplex gene expression profiles in steroid resistant in comparison with responder patients, but it is not clear if this different expression is the cause of the variability in response or the consequence of GC treatment[54-59]. After GC binding and dissociation from heterocomplex proteins, the GR translocates into the nucleus; translocation is mediated by specific nuclear transport factors that belong to the importin β family of nuclear transporters, and in particular by importin 13[60]. The activated receptor then binds as homodimer two palindromic DNA-binding sites, the so-called glucocorticoid responsive elements (GREs), localized in the promoter region of target genes[61-63]. As a consequence of DNA binding, GCs can induce trans-activation and trans-repression processes: binding to positive GREs leads to activation of the transcription of anti-inflammatory [e.g., interleukin 10 (IL-10), Annexin 1] as well as of regulator proteins involved in metabolic processes (e.g., enzymes of gluconeogenesis)[64-66]. The second mechanism of GC action is trans-repression[67], which leads to a reduced expression of immune-regulatory and proinflammatory proteins such as cytokines [IL-1, IL-2, IL-6, tumor necrosis factor-α (TNF-α)] and prostaglandins[68], and is believed to be responsible for the majority of beneficial anti-inflammatory effects.

Steroid hormones can regulate gene expression post-transcriptionally, by destabilizing mRNAs[69]. In addition, these hormones can induce rapid non genomic effects within the cytoplasm; for example, they induce the release of Src kinase from the GR heterocomplex, resulting in lipocortin activation and inhibition of arachidonic acid release[70,71], and alter cytoplasmic ion content[72,73].

miRNAS AND GC RESPONSE

miRNA regulation by GCs

It has been demonstrated that activation of GR by GCs might induce or repress specific miRNAs in various target genes. The majority of studies have evaluated the effect of GCs on miRNA expression levels in tumor leukemic cells, during GC induced apoptosis[74].

Rainer et al[75] have correlated miRNA levels with expression data of their host genes in cell lines and clinical samples of children with acute lymphoblastic leukemia (ALL) undergoing systemic GC monotherapy. At least 5 miRNAs were significantly regulated by GC therapy. Importantly, the miR-15/16 cluster, which induces cell cycle arrest, was up-regulated by GCs in a subset of ALL patients and cell lines, consistent with the known apoptotic effect of GCs in immature lymphoblasts. Indeed, overexpression of miR-15b/16 increased GC sensitivity in leukemia cell lines whereas silencing miR-15b/16 with inhibitors decreased GC sensitivity in vitro.

Another study in a T-cell lymphoma cell line has shown that GC treatment repressed the expression of the miRNA cluster miR-17-92, which results in elevated protein expression of Bim, a proapoptotic member of the B-cell lymphoma-2 family (Bcl-2). Overexpression of miRNA cluster miR-17-92 decreased Bim induction, and attenuated GC mediated apoptosis, while cluster knockdown increased Bim induction and GC mediated apoptosis[76]. These findings suggest a novel mechanism that could contribute to the induction of lymphocyte apoptosis by GCs.

Harada et al[77] demonstrated that in the leukemic cell line RS4; 11 dexamethasone down-regulated miRNA levels; miR17HG was rapidly down-regulated, and chromatin immunoprecipitation demonstrated that the promoter is a target of GC transcriptional repression; in particular, the miR-17-92 cluster was identified as a prime target for dexamethasone induced repression. In the sensitive leukemia cell line SUP-B15, but not in the resistant line REH, dexamethasone reduced the expression of the mi-R17 family and concomitantly increased its target protein Bim. Up-regulation or inhibition of miR-17 resulted in a decrease and increase, respectively in Bim protein levels and in dexamethasone induced cytotoxicity. Down-regulation of miR-17 levels was observed in ex vivo patients’ leukemia cells that underwent dexamethasone induced apoptosis[77].

Another recent study[78], by genome wide miRNA microarray on diagnostic bone marrow samples of ALL pediatric patients treated with GCs, identified a reduced expression of miR-355 as the most significant miRNA abnormality associated with poor outcome. Moreover, the authors demonstrated that exogenous expression of miR-355 in ALL cells increases sensitization to prednisolone-induced apoptosis. MAPK1 was identified as a target of miR-355, and the MEK/ERK inhibitor treatment increased GC induced cytotoxicity through the activation of Bim.

Smith et al[79] have demonstrated that miRNAs are repressed during GC induced apoptosis of primary rat thymocytes, and further demonstrated the repression of the miRNA processing enzymes Dicer, Drosha and DGCR8/Pasha. Silencing of Dicer expression in two human leukemic lines significantly enhanced GC induced apoptosis, while overexpression of the GC-repressed miR-17-92 polycistron reduced apoptosis.

Among the few studies that have considered the effect of GCs on miRNA expression in non tumor cells, Ledderhose et al[80] in native and CD3/CD28 stimulated cells from healthy volunteers, demonstrated that miR-24 is expressed in human T cells, and expression is increased 1.7 fold upon stimulation. Hydrocortisone significantly enhanced by 3 fold the miRNA induction[80].

In human corneal fibroblast treated for 16 h with dexamethasone, genome microarray and microRNA analyses were used to evaluate gene and miRNA expression. In response to treatment with the steroid, 261 genes were up-regulated and 123 were down-regulated more than threefold. Several miRNAs, including miR-16, miR-21 and miR-29C were up-regulated, whereas miR-100 was down-regulated by the steroid, suggesting a posttranscriptional control of gene expression through miRNAs[81].

Studies of the miRNAs profile on mucosal biopsies of patients with eosinophilic esophagitis, before and after successful treatment with GCs were conducted by Lu and collaborators[82]; of the 377 miRNA sequences examined, 32 miRNAs were significantly up-regulated and 4 down-regulated in the biopsies obtained before treatment compared to samples obtained after GC therapy. miR-214 was the most up-regulated (150 fold) and miR-146b-5b, 146a, 145, 142-3p and 21 were up-regulated by at least 10 fold.

Williams et al[83], using a highly sensitive reverse transcription-polymerase chain reaction, measured 277 miRNAs in airway biopsies obtained from normal subjects and mild asthmatic patients before and after one month twice dayly treatment with inhaled budesonide. No significant difference in miRNA expression was evident in the airway biopsies of normal and asthmatic subjects, and, despite improved lung function, no change in miRNAs expression was evident after one month budesonide treatment. However, a specific miRNA expression profile was observed in different cell types (alveolar epithelial cells, airway smooth muscle cells, alveolar macrophages, lung fibroblasts).

Finally, in a recent study[84], activated human CD4⁺ T cells from healthy donors were exposed in vitro to 1 µmol/L of methylprednisolone and changes in miRNA and mRNA expression profiles were analyzed by microarrays; a number of steroid responsive genes and miRNAs were identified. Further studies with qPCR, flow cytometry and ELISA, demonstrated that methylprednisolone increased the expression of miR-98 and suppressed the levels of predicted targets, including the pro-inflammatory cytokine IL-13 and three TNF receptors FAS, FASL, and TNF receptor superfamily member 1B (TNFRSF1B): these data suggest that methylprednisolone acts through miR-98 to inhibit specific pro-inflammatory targets[84].

GR as miRNA target

The role of miRNAs in the regulation of the GR has been examined, indeed, computational studies showed that the 3’ UTR of the GR is predicted to contain numerous seed regions recognized by a variety of miRNAs[85].

Using a combination of in silico prediction of miRNA binding sites, miRNA overexpression studies and mutagenesis of the GR 3’UTR, Vreugdenhil and collaborators[86] found that miR-18 and miR-124a bind GR mRNA and decrease GR activity in neuronal tissues. These miRNAs were tested for their ability to alter the translational activity of GR and reduce GR protein levels in cell cultures in vitro; miR-18 and miR-124a overexpression reduced GR protein levels and impaired the activation of the GR responsive gene glucocorticoid-induced leucine zipper (GILZ). In addition these authors have demonstrated by miRNA reporter assay that miR-124a is able to bind to the predicted seed region in the GR 3’UTR.

Ledderose et al[80] have investigated the role of miR-124 in the regulation of GR expression; these authors have studied the influence of the GR isoforms (the active isoform α, and the dominant negative non-ligand-binding isoform β) on GC effects in human T-cells, and found that, in patients with critical illness-related corticosteroid insufficiency, miR-124 specifically down-regulated GR-α: a slight increase of miR-124 and a reduction of GR-α was observed in patient T-cells compared to healthy controls. The authors suggested a novel miR-124-mediated mechanism in the down-regulation of GR-α in patients with critical illness-related corticosteroid insufficiency, that could explain, at least in part, GC resistance in this disease.

Tessel et al[87] have identified and characterized miR-130b as an important down-regulator of GR in GC-resistant multiple myeloma cell line: the overexpression of this miRNA was also associated with a decreased regulation of the downstream GC controlled gene GILZ, suggesting this mechanism as one of the possible causes of resistance to GCs.

miRNA involved in IBD

The pathophysiology of IBD is not yet clear, and genetic, epigenetic, infectious and immunological factors seem to play a role. It has been suggested that the gastrointestinal inflammation is the result of an altered activation of the immune system to a luminal factor, such as intestinal flora, in genetically predisposed subjects.

Among the many biological processes regulated by miRNAs, it is now accepted that these small non coding RNAs contribute to the maintenance of immunological homeostasis at mucosal sites[88,89]. The role of miRNAs in the pathogenesis of IBD has been thoroughly considered (see recent reviews[32,90,91]), and it has been suggested that these small non coding RNAs represent an important player in the complex interactions which results in IBD clinical features. Of particular interest is the observation that miRNA expression changes during tissue progression from normal to inflamed and varies according to the type and evolutionary stage of IBD[92]. Indeed, a number of studies have identified a specific differential expression of miRNAs in IBD and unique miRNA expression profiles for the different subtypes of IBDs, both in human tissues collected by colonoscopic biopsies and in peripheral blood samples, have been demonstrated[32,90,91].

It has been argued that genetic polymorphisms in miRNAs, as well as in miRNA target genes can affect their regulatory function and, consequently, the expression level of their target mRNAs. Most studies have described an association between SNPs in miRNA genes and human cancers[93-98], and only recently the association between mRNA related SNPs and the risk of IBD has been examined[99]. Bioinformatic approaches have been used to analyze the association between diseases-linked SNPs, miRNAs and mRNAs: SNP data derived from genome wide association studies that were correlated with miRNA, revealed a CD phenocode comprising rs11209026, rs7807268, rs254215, rs2542151 in miR-125, rs11805303 in miR-519, and rs6908425 in miR-181[30]. Of interest, miR-181, miR-519 and miR-119 could target mRNAs encoded by genes involved in the importin pathway, whereas miR-181 and miR-125 are potential regulators of components of inflammasome pathway. Both importin and inflammasome are involved also in GC molecular mechanism: importin is a nuclear transport protein responsible for the translocation of the complex GR-GC into the nucleus[2], and variants in inflammasome gene have been correlated with steroid resistance in pediatric IBD patients[100].

An association between rs3746444 in miR-499 and UC susceptibility has been observed in 170 Japanese patients: this SNP may alter the function or expression of miR-499, altering the regulation of target mRNAs related to inflammatory immune responses, and influencing the pathophysiological features of UC[101]. Of particular interest is the observation that the rs3746444 AG genotype was associated also with steroid dependence and refractory phenotype, whereas the rs3746444 AA genotype was inversely related to hospitalization time, steroid dependence, and refractory phenotype. In addition, the rs11614913 TT genotype held a significantly higher risk of refractory phenotype.

CONCLUSION

There is a lot of interest in identifying the role of miRNAs in the modulation of drug response, but studies about this topic are still very limited, and the possible correlation between miRNAs expression and variability in GC response in IBD patients has not yet been extensively examined. Studies about miRNAs and pharmacogenomics may represent a promising investigation topic that could increase the understanding of the pharmacology of steroids in IBDs and possibly in other diseases.

A number of studies have shown that GCs can modify the expression profile of different miRNAs, however, the obtained results have been highly variable. The differences observed can possibly be ascribed to the different tissues or cell lines analysed or different experimental protocols, and to date it is not possible to recognize a specific miRNA pattern regulated by GCs.

miRNA regulation by GCs in IBDs has never been analyzed in clinical prospective studies, in which patients are followed from diagnosis and throughout steroid therapy: the identification of miRNAs differently expressed as a consequence of GC treatment in comparison to diagnosis, could be an important innovative approach. This type of study design will reduce to the minimum the effect of confounding factors and results should be easier to translate into clinical practice.

Moreover, existing studies employ techniques based on the use of reverse transcription quantitative PCR and microarrays, based on the analysis and quantification of already known miRNAs. Using next generation sequencing technologies it should be possible to detect novel, still unrecognised miRNAs, and identify new miRNA isoforms (iso-miRs) as well.

In the future, the increased availability and the reduced costs of RNA profiling should enable the clinicians to stratify patients on specific miRNA biomarkers before starting GC treatment. This will allow the personalization of therapy, avoiding a treatment doomed to failure, increasing efficacy and reducing toxicity.