Experimental Diabetes ResearchVolume 2012 (2012), Article ID 695214, 7 pagesdoi:10.1155/2012/695214

Research Article

Regulation of Pancreatic microRNA-7 Expression

Sharon Kredo-Russo,¹ Avital Ness,1,2 Amitai D. Mandelbaum,¹ Michael D. Walker,² and Eran Hornstein¹

¹Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 76100, Israel²Department of Biological Chemistry, Weizmann Institute of Science, Rehovot 76100, Israel

Received 12 January 2012; Accepted 8 March 2012

Academic Editor: Francis Lynn

Copyright © 2012 Sharon Kredo-Russo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Genome-encoded microRNAs (miRNAs) provide a posttranscriptional regulatory layer, which is important for pancreas development. Differentiation of endocrine cells is controlled by a network of pancreatic transcription factors including Ngn3 and NeuroD/Beta2. However, how specific miRNAs are intertwined into this transcriptional network is not well understood. Here, we characterize the regulation of microRNA-7 (miR-7) by endocrine-specific transcription factors. Our data reveal that three independent miR-7 genes are coexpressed in the pancreas. We have identified conserved blocks upstream of pre-miR-7a-2 and pre-miR-7b and demonstrated by functional assays that they possess promoter activity, which is increased by the expression of NeuroD/Beta2. These data suggest that the endocrine specificity of miR-7 expression is governed by transcriptional mechanisms and involves members of the pancreatic endocrine network of transcription factors.

1. Introduction

The development of the endocrine pancreas is governed by a network of transcription factors that specify different endocrine cell types, including insulin-producing beta cells, glucagon-producing alpha cells, somatostatin-producing delta cells, pancreatic polypeptide-producing PP cells, and ghrelin-producing epsilon cells [1–3]. The endocrine differentiation program is initiated by neurogenin3 (Ngn3) [4, 5]. Next, a complex network of transcription factors is activated to differentially specify the endocrine lineages (reviewed in [3, 6]).

Genome-encoded miRNAs act in concert with transcription factors to refine gene expression and confer robustness to developmental transitions [7–10]. Many miRNA genes are nested within introns of protein-coding genes and are subjected to transcriptional control with their host gene [11]. However, other miRNA genes are located in intergenic regions and are expressed autonomously. For example, in a previous study, we characterized the pancreas-enriched miR-375 and demonstrated that cell specificity is controlled transcriptionally through well-defined cis-regulatory elements [12].

Like miR-375, miR-7 is highly and selectively expressed in the endocrine pancreas of zebrafish, mouse, and human [13–16]. miR-7 is an evolutionarily conserved miRNA, encoded by a single gene in flies and by three different genomic loci in mammals. In mice two miR-7 genes are located in intergenic regions of Chr. 7 (mmu-mir-7a-2) and Chr. 17 (mmu-mir-7b), whereas a third miR-7 gene, mmu-mir-7a-1, is embedded within an intron of the gene encoding for the RNA-binding protein, Hnrnpk (MGI: 99894, on Chr. 13).

The two miR-7a genes generate an identical 22nt mature sequence, whereas miR-7b differs by a single nucleotide. However, functionally the three genes are identical, as they harbor the same “seed” sequence. Hence, all miR-7 genes coregulate the very same target set. While a defined set of targets is suggested for miR-7 [17] some of which have been experimentally validated [18], little is known about miR-7 promoter structure or the mechanisms controlling miR-7 expression.

In this study we characterized elements within the mmu-mir-7a-2 and mmu-mir-7b upstream regulatory sequences. We show that miR-7 responds to Ngn3 directly. However, our data suggest that NeuroD/Beta2, the primary effector of Ngn3, controls miR-7 and is probably responsible for maintenance of miR-7 expression in differentiated endocrine cells.

2. Materials and Methods

2.1. Quantitative PCR for Precursor and Mature miRNA

Extraction of total RNA was carried out by the miRNeasy Mini Kit (Qiagen). For precursor quantification, synthesis of cDNA was performed using miScript system (Qiagen). cDNA was synthesized from miRNAs using Taqman MicroRNA qPCR Assays (Applied Biosystems). qPCR analysis of mRNA was performed on LightCycler 480 System (Roche) using Kapa SYBR Green qPCR kit (Finnzymes). miRNA levels were normalized to the expression of small RNAs (sno234 and U6) and mRNA normalized to GAPDH and HPRT. (primer sequences are described in Supplementary Table 1 available online at doi: 10.1155/2012/695214).

2.2. Cell Culture and Luciferase Reporter Assay

HEK-293T cells (American Type Culture Collection), betaTC-3 (a gift from Shimon Efrat), and MIN6 cells (a gift from Jun-ichi Miyazaki) were grown in Dulbecco’s Modified Eagle Medium (DMEM) with 10% FBS, 2 mM L-glutamine, 100 U/mL penicillin/streptomycin at 37°C at 5% CO₂ in a humidified incubator. Experiments on MIN6 cells were performed between passages 18 to 28.

For miR-7 promoter analyses, fragments of 256 bp, 206 bp, and 156 bp (representing miR-7a-2 “block 1”, miR-7a-2 “block 2”, and “block 3” of miR-7b, resp.) were subcloned into pGL3-basic, using restriction enzymes BglII and KpnI. Primer sequences are described in Table  S1. HEK-293T cells were transfected with 100 ng of the reporter, 20 ng of A20 Renilla reporter, and in addition 50 ng of pcDNA3 empty and/or NeuroD/Beta2 and Ngn3 expression vectors. Reporter activity was measured 48 h after transfection with the Dual-Luciferase Reporter Assay System (Promega).

For overexpression analysis, expression vectors for transcription factors, pcDNA3 empty vector, and CMV-GFP vectors were transfected to MIN6 cells using Lipofectamine 2000 reagent (Invitrogen), according to the manufacture’s instructions. miR-7 endogenous expression was analyzed by qPCR 48 h later.

2.3. Statistical Analysis

Analysis was performed using either Student’s  -test or two-way ANOVA by the JMP software. Results are given as mean ± SEM. The null hypothesis was rejected at the 0.05 level (**).

3. Results and Discussion

3.1. Three miR-7 Loci Are Expressed in Endocrine Cells

Since miR-7 gene has three genomic copies in mouse and human, we first determined which of them is expressed in endocrine cells. As previously shown [19], miR-7a and miR-7b were highly and specifically expressed in islets of Langerhans, relative to other organs such as heart and brain (Figures 1(a) and 1(b)). The precursors transcribed from the three different loci can be distinguished by quantitative real-time PCR (qPCR). Therefore, we designed specific primers for each of miR-7 precursors (for oligo sequences, see Table  S1) and performed a qPCR study in BetaTC3 cell line. This analysis revealed that all three miR-7 genes are expressed in cultured beta cells (Figure 1(c)). Consistent with these results, qPCR for the mature miR-7a and miR-7b revealed expression of both forms of miR-7 in beta cells (Figure 1(d)). Taken together, the analysis of miR-7 expression suggests that the three miR-7 loci are activated in beta cells and are responsible for the overall high expression of the mature miR-7.

(a)

(b)

(c)

(d)

Figure 1: Three miR-7 genes are expressed in beta cells. (a) and (b) Taqman qPCR analysis for mature miR-7a and miR-7b in several tissues: the aorta and pulmonary artery, heart, isolated islets, and brains from three months old mice. Data normalized to sno234. (c) qPCR analysis for miR-7 precursors reveals expression in beta-TC cultured cells. Data normalized to GAPDH and 18S. independent measurements, duplicates each. (d) Taqman qPCR analysis to mature miR-7a and miR-7b in beta-TC cultured cells. Data normalized to sno234. independent measurements in duplicates. Error bars represent ± SEM.

3.2. Sequence Analysis of Potential miR-7 Regulatory Regions

Comparative genomics is commonly used to identify conserved sequences of functional importance [12, 20]. Therefore, we compared the region upstream of the pre-miR-7 sequences among several vertebrate orthologs and identified three highly conserved sequences, which may be potentially functional promoter regions. Trimethylation of lysine 4 of histone 3 (H3K4) is involved in activation of transcription and marks the position of transcriptional start site of many genes, including insulin [21]. Accordingly, we identified typical H3K4 methylation pattern [22, 23], indicating possible transcriptional start sites, at the vicinity of miR-7 genomic regions.

For miR-7a-2 we identified two conserved sequences: the 256 bp long “block 1” and 206 bp “block 2,” which are located 1420 bp and 450 bp upstream of the miR-7a-2 hairpin, respectively (Figure 2(a)). For miR-7b, a single 156 bp long “block 3” conserved sequence was identified, positioned 1,235 bp upstream of the miR-7b hairpin (Figure 2(b)). miR-7a-1 was omitted from this analysis since it is embedded in an intron of the Hnrnpk gene and is likely regulated by the promoter of the host gene.

(a)

(b)

Figure 2: Conserved sequences upstream of miR-7 genes. Screen capture from the UCSC Genome Browser reveals fragments of Mus musculus chromosome 7 (a) and chromosome 17 (b). Tracks for conservation within mammals and the annotation of the pre-miR-7 sequences are depicted. Dashed squares demarcate the sequences that were further investigated. Scale bars represent 500 bp.

3.3. Activity of the miR-7 Promoters in Cultured Cells

To functionally characterize miR-7 promoters, the fragments consisting of blocks 1–3 were fused separately to a firefly luciferase reporter gene, on a promoterless plasmid (pGL3-basic), and transfected into the pancreatic beta cell line MIN6 and the embryonic kidney line HEK 293.

In HEK-293T cells, “block 2” produced a twofold increase in luciferase activity, relative to pGL3-basic reporter (“Ctrl”), indicating weak promoter activity (Figure 3(b), middle panel). “Block 3”, on the other hand, produced a large (200-fold) increase in luciferase activity, indicating the presence of a strong promoter in this region (Figure 3(b), right panel). Surprisingly “block 1” resulted in repression of the luciferase activation, suggesting that this conserved sequence may be involved in repression of gene transcription (Figure 3(b), left panel). We then examined whether the putative promoter sequences could drive luciferase activation in MIN6 beta cells. We observed that only the DNA sequences dubbed “block 3” showed significant promoter activity (Figure 3(c)). Taken together, these results show that “block 2” and “block 3” possess weak and strong promoter activity, respectively. The activity of these fragments is not restricted to pancreatic cells, and presumably these promoter fragments require additional elements to confer selectivity in vivo.

(a)

(b)

(c)

Figure 3: miR-7 putative promoters induce reporter activation. (a) Schematic representation of conserved sequences, upstream of pre-miR-7a-2 hairpins (“block 1” and “block 2”) and pre-miR-7b (“block 3”) that were subcloned upstream of pGL3 basic luciferase reporter. (b) Relative activation of luciferase firefly to renilla in HEK-239 cells, transfected with various reporters. pGL3-basic serves as control (“Ctrl”). (c) Relative luciferase activation in MIN6 beta cells transfected with the indicated reporters. experiments, in triplicates. Error bars represent ± SEM.

3.4. bHLH Transcription Factors Directly Induce miR-7 Promoter Activity

It has been previously shown that the bHLH transcription factors, Ngn3 and NeuroD1, play a central role in pancreas endocrine development and in mature beta-cell function, respectively. These proteins function through E-boxes (consensus sequence CAxxTG) located in target gene promoters [24]. Since two conserved E-box elements were identified in miR-7 promoter sequences (Figure 3(a)), we tested whether Ngn3 and NeuroD1 expression regulates miR-7 promoter activity.

For this, we performed transactivation experiments, in which luciferase reporter constructs were cotransfected into HEK-293T cells in the presence of expression vectors encoding the endocrine transcription factors. These analyses showed activation of “block 2” promoter in response to both NeuroD/Beta2a and Ngn3. In the presence of both factors, an additive effect was observed (Figure 4(b)). With “block 3”, NeuroD/Beta2 produced significant activation, whereas Ngn3 had little or no effect. The expression of block 1-containing promoter was not significantly affected by any of the co-transfected transacting factors.

(a)

(b)

(c)

Figure 4: Transactivation of miR-7 reporters in HEK 293T. (a) The luciferase activity of a reporter driven by “block 1” did not demonstrate response to Ngn3 or NeuroD/Beta2. However, reporters driven by “block 2” are induced by both Ngn3 and NeuroD/Beta2 (b). Highest luciferase response is induced by a combination of Ngn3 and NeuroD/Beta2. (c) “Block 3” transcription is transactivated only by NeuroD/Beta2 but not by Ngn3. Luciferase activity normalized to A20-Renilla expression and to transactivation by a control pcDNA3-empty vector. experiments, in triplicates. Error bars represent ±SEM (**).

Our study, therefore, suggests that the sequence upstream of miR-7b is able to activate transcription in cultured beta cells and can be activated by the endocrine transcription factor NeuroD/Beta2. The block 2 region upstream of miR-7a-2 can be activated both by NeuroD/Beta 2 and by Ngn3. However, this sequence shows very low promoter activity in the beta cell line MIN6 and therefore probably works in concert with other elements in order to transcribe the miR-7a-2 locus in beta cells.

A single ancestral miR-7 gene in invertebrates has undergone genomic duplication in the vertebrate clade. While duplication of genes often leads to functional divergence of each locus, our data suggest that miR-7 expression from three independent loci contributes primarily to higher levels of expression. Furthermore, coregulation of miR-7a-2 and miR-7b by NeuroD/Beta2 suggests that these genes respond to similar transacting factors in the endocrine pancreas.

3.5. Endogenous miR-7 Expression Is Activated by NeuroD/Beta2 in Cultured Beta Cells

Finally, we determined whether expression of endogenous miR-7 can be modified by introducing endocrine transcription factors into beta cell lines. To this end, we transfected MIN6 cells with expression vectors for Pdx1, Ngn3, NeuroD/Beta2, Hnf1b, Insm1, and E2A. Of these, miR-7 expression was induced in cultured beta cells only by NeuroD/Beta2 (Figure 5(a)). This is consistent with the results obtained from reporter experiments and suggests that selective expression of miR-7 is controlled by lineage-specific transcription factors, primarily, NeuroD/Beta2. This supports the hypothesis that miR-7 itself may be a component of this cascade with functional roles in controlling endocrine cell development at a posttranscriptional level (see model in Figure 5(b)).

(a)

(b)

Figure 5: NeuroD/Beta2 regulates miR-7. (a) miR-7 expression in MIN6 cells (passages 22–24), transfected with expression vectors for various transcription factors. Endogenous miR-7 expression is upregulated upon NeuroD/Beta2 introduction, relative to GFP-expressing control vector. (b) A Schema describing miR-7 regulators.

4. Conclusions

During development, Ngn3 induces endocrine cell differentiation by upregulating transcription factors such as NeuroD/Beta2, Pax4, Arx, and Pax6. Indeed, loss of Ngn3 causes blockade in endocrine differentiation [5]. At the same time, Ngn3 also upregulates miR-7 expression through its effector, the transcription factor NeuroD/Beta2, that is likely involved in maintaining miR-7 expression also in mature cells.

In summary, our analysis identifies miR-7 as a novel component downstream of Ngn3 and NeuroD, embedded within the transcription factor network regulating pancreas development. Further dissection of the transcriptional mechanisms controlling expression of miR-7 and other endocrine miRNAs will contribute substantially to our overall understanding of the role of miRNA in pancreas development and function.

Authors’ Contribution

S. K. Russo, M. D. Walker, and E. Hornstein designed the experiments. S. K. Russo, A. Ness and A. D. Mandelbaum conducted all experiments and analyses in this work. S. K. Russo, M. D. Walker, and E. Hornstein, wrote the paper.

Acknowledgments

The authors thank Tali Avnit-Sagi for protocols, reagents, and productive discussions. They thank Jun-ichi Miyazaki for MIN6 cells. This paper was supported by grants to E. Hornstein from Juvenile Diabetes Research Foundation (no. 99-2007-71), GIF Young Investigator Award, the EFSD/D-Cure Young Investigator Award and EFSD Lilly Program, the Israel Science Foundation, and the Wolfson Family Charitable Trust. E. Hornstein is the incumbent of the Helen and Milton A. Kimmelman Career Development Chair and the Dr. Sydney Brenner Chair. M. D. Walker is the incumbent of the Marvin Meyer and Jenny Cyker Chair of Diabetes Research at the Weizmann Institute. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the paper.

References

B. M. Lyttle, J. Li, M. Krishnamurthy et al., “Transcription factor expression in the developing human fetal endocrine pancreas,” Diabetologia, vol. 51, no. 7, pp. 1169–1180, 2008.

M. Martin, V. Hauer, M. Messmer, C. Orvain, and G. Gradwohl, “Transcription factors in pancreatic development: animal models,” Endocrine Development, vol. 12, pp. 24–32, 2007.

C. Bonal and P. L. Herrera, “Genes controlling pancreas ontogeny,” International Journal of Developmental Biology, vol. 52, no. 7, pp. 823–835, 2008.

G. Gu, J. Dubauskaite, and D. A. Melton, “Direct evidence for the pancreatic lineage: NGN3+ cells are islet progenitors and are distinct from duct progenitors,” Development, vol. 129, no. 10, pp. 2447–2457, 2002.

G. Gradwohl, A. Dierich, M. LeMeur, F. Guillemot, and F. Guillemot, “Neurogenin3 is required for the development of the four endocrine cell lineages of the pancreas,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 4, pp. 1607–1611, 2000.

A. M. Ackermann and M. Gannon, “Molecular regulation of pancreatic β-cell mass development, maintenance, and expansion,” Journal of Molecular Endocrinology, vol. 38, no. 1-2, pp. 193–206, 2007.

D. P. Bartel, “MicroRNAs: genomics, biogenesis, mechanism, and function,” Cell, vol. 116, no. 2, pp. 281–297, 2004.

E. Hornstein and N. Shomron, “Canalization of development by microRNAs,” Nature Genetics, vol. 38, pp. S20–S24, 2006.

A. Stark, J. Brennecke, N. Bushati, R. B. Russell, and S. M. Cohen, “Animal microRNAs confer robustness to gene expression and have a significant impact on 3′UTR evolution,” Cell, vol. 123, no. 6, pp. 1133–1146, 2005.

X. Li, J. J. Cassidy, C. A. Reinke, S. Fischboeck, and R. W. Carthew, “A microRNA imparts robustness against environmental fluctuation during development,” Cell, vol. 137, no. 2, pp. 273–282, 2009.

H. K. Saini, S. Griffiths-Jones, and A. J. Enright, “Genomic analysis of human microRNA transcripts,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 45, pp. 17719–17724, 2007.

T. Avnit-Sagi, L. Kantorovich, S. Kredo-Russo, E. Hornstein, and M. D. Walker, “The promoter of the pri-miR-375 gene directs expression selectively to the endocrine pancreas,” PLoS One, vol. 4, no. 4, Article ID e5033, 2009.

E. Wienholds, W. P. Kloosterman, E. Miska et al., “Cell biology: microRNA expression in zebrafish embryonic development,” Science, vol. 309, no. 5732, pp. 310–311, 2005.

P. Landgraf, M. Rusu, R. Sheridan et al., “A mammalian microRNA expression atlas based on small RNA library sequencing,” Cell, vol. 129, no. 7, pp. 1401–1414, 2007.

M. Correa-Medina, V. Bravo-Egana, S. Rosero et al., “MicroRNA miR-7 is preferentially expressed in endocrine cells of the developing and adult human pancreas,” Gene Expression Patterns, vol. 9, no. 4, pp. 193–199, 2009.

V. Bravo-Egana, S. Rosero, R. D. Molano et al., “Quantitative differential expression analysis reveals miR-7 as major islet microRNA,” Biochemical and Biophysical Research Communications, vol. 366, no. 4, pp. 922–926, 2008.

B. P. Lewis, I. H. Shih, M. W. Jones-Rhoades, D. P. Bartel, and C. B. Burge, “Prediction of microRNA targets,” Cell, vol. 115, no. 7, pp. 787–798, 2003.

B. Kefas, J. Godlewski, L. Comeau et al., “microRNA-7 inhibits the epidermal growth factor receptor and the akt pathway and is down-regulated in glioblastoma,” Cancer Research, vol. 68, no. 10, pp. 3566–3572, 2008.

M. V. Joglekar, V. M. Joglekar, and A. A. Hardikar, “Expression of islet-specific microRNAs during human pancreatic development,” Gene Expression Patterns, vol. 9, no. 2, pp. 109–113, 2009.

L. A. Pennacchio and E. M. Rubin, “Genomic strategies to identify mammalian regulatory sequences,” Nature Reviews Genetics, vol. 2, no. 2, pp. 100–109, 2001.

J. Francis, S. K. Chakrabarti, J. C. Garmey, and R. G. Mirmira, “Pdx-1 links histone H3-Lys-4 methylation to RNA polymerase II elongation during activation of insulin transcription,” The Journal of Biological Chemistry, vol. 280, no. 43, pp. 36244–36253, 2005.

J. Harrow, F. Denoeud, A. Frankish et al., “GENCODE: producing a reference annotation for ENCODE,” Genome Biology, vol. 7, 1, pp. S41–S49, 2006.

P. Flicek, M. R. Amode, D. Barrell et al., “Ensembl 2011,” Nucleic Acids Research, vol. 39, pp. D800–D806, 2011.

E. Glick, D. Leshkowitz, and M. D. Walker, “Transcription factor BETA2 acts cooperatively with E2A and PDX1 to activate the insulin gene promoter,” The Journal of Biological Chemistry, vol. 275, no. 3, pp. 2199–2204, 2000.