Aberration of the modulatory functions of intronic microRNA hsa-miR-933 on its host gene ATF2 results in type II diabetes mellitus and neurodegenerative disease development
Human Genomics volume 14, Article number: 34 (2020)
MicroRNAs are ~ 22-nucleotide-long biological modifiers that act as the post-transcriptional modulator of gene expression. Some of them are identified to be embedded within the introns of protein-coding genes, these miRNAs are called the intronic miRNAs. Previous findings state that these intronic miRNAs are co-expressed with their host genes. This co-expression is necessary to maintain the robustness of the biological system. Till to date, only a few experiments are performed discretely to elucidate the functional relationship between few co-expressed intronic miRNAs and their associated host genes.
In this study, we have interpreted the underlying modulatory mechanisms of intronic miRNA hsa-miR-933 on its target host gene ATF2 and found that aberration can lead to several disease conditions. A protein-protein interaction network-based approach was adopted, and functional enrichment analysis was performed to elucidate the significantly over-represented biological functions and pathways of the common targets. Our approach delineated that hsa-miR-933 might control the hyperglycemic condition and hyperinsulinism by regulating ATF2 target genes MAP4K4, PRKCE, PEA15, BDNF, PRKACB, and GNAS which can otherwise lead to the development of type II diabetes mellitus. Moreover, we showed that hsa-miR-933 can regulate a target of ATF2, brain-derived neurotrophic factor (BDNF), to modulate the optimal expression of ATF2 in neuron cells to render neuroprotection for the inhibition of neurodegenerative diseases.
Our in silico model provides interesting resources for experimentations in a model organism or cell line for further validation. These findings may extend the common perception of gene expression analysis with new regulatory functionality.
MicroRNAs (miRNA) are short, single-stranded ~ 22-nucleotides-long RNA molecules, which are partially complementary to one or more messenger RNA (mRNA) molecules known as target mRNAs . In humans, hundreds of miRNA genes are predicted to be present, and so, the potential regulatory circuitry afforded by miRNA is huge [2, 3]. They can either downregulate the gene expression  or can also upregulate the translation of mRNAs . The expression of almost 20–30% of all protein-encoding genes may be altered by miRNAs at this post-transcriptional level regulation . miRNAs may act as key regulators of processes as miscellaneous as embryonic development, cell proliferation, cell growth, tissue differentiation, and apoptosis. Recent studies of miRNA expression involve miRNAs in cellular signaling networks and co-regulation with transcription factors. Accordingly, a mutation in miRNAs, dysfunction of miRNA biogenesis, and dysregulation may result in a broad spectrum of diseases. In addition, components required for miRNA processing and/or function have also been implicated in various disorders. Currently, there have been reported ~ 378 diseases which are associated with miRNAs .
Intronic miRNAs can be defined based on two factors; first, they must share the same promoter with their encoded target genes, and second, they are spliced out of the transcript of such encoded genes and further processed into mature miRNAs . About 37% of the known human miRNAs are located within the introns of protein-coding genes preferably known as host genes . About 26% of the human intronic miRNAs are transcribed from their own promoters . But the majority of human intronic miRNAs are transcriptionally linked to their host gene expression and processed from the same primary transcript . In humans, most of the intronic miRNAs also show correlated expression with their host genes . Besides Drosha-processed miRNAs, the second type of intronic miRNAs, mirtrons are discovered that bypass the Drosha cleavage by splicing  but exhibit the same co-expression patterns with their host genes.
Intronic miRNAs can negatively regulate their host genes by targeting the 3′-UTR of their host gene, inhibiting the host gene’s targets, or inhibiting the transcription of their host genes by a negative feedback loop. On the other hand, some intronic miRNAs can act as a positive regulator of their host gene by forming a positive feedback loop that upregulates the function of its host gene, working in concert with the host gene’s targets, or silencing antagonistic genes to its host gene .
Previous findings [15, 16] showed that the coupled expression of intronic miRNA and host gene were observed where miRNA can modulate the function of its host gene. Moreover, Steiman-Shimony et al. showed that intronic miRNAs may target transcripts whose genes/proteins are targeted by the host gene which can code a transcription factor , but no functional consequence of those interplays was revealed. Experiments to elucidate the functional regulations of the host gene and intronic miRNA were done in minuscule to date as they are time-consuming and extravagant [15, 16].
Like many other regulatory miRNAs, an important association of intronic miR-933 in various diseases like dementia , hyperlipidemia and cardiovascular diseases , and gastric cancers  was reported. It was also found that ATF2 and miR-933 share a common promoter . It is evident from the previous studies that the ATF2 transcription factor is involved in many diseases ranging from inflammatory diseases, diabetes, multiple neurodegenerative pathologies, different forms of cancers, etc. . But the link between intronic miR-933 and host gene ATF2 still remained elusive.
Although there are several intronic miRNAs that exist in the human genome, systemic identification of intronic miRNAs and their role in normal physiology and disease pathobiology is only explored in few cases . In this study, after genome-wide identification of all probable intronic miRNAs, we explored the relationship between intronic miRNA hsa-miR-933 and its host genes with a combination of functional enrichment analysis to elucidate the functional relationship between them and the mechanism of these functional regulations. Furthermore, it was anticipated that the host gene’s and intronic miRNA’s common target gene sets for a particular biological process or pathway can act antagonistically or synergistically.
A total of 822 intronic miRNAs and their associated host genes were identified
From the human genome, we identified a total of 822 intronic miRNAs (data not shown). Among these, hsa-miR-933 was identified as embedded in the first intron of the host protein-coding gene activating transcription factor 2 (ATF2). ATF2 is a protein that acts as a transcriptional activator which regulates the transcription of a variety of genes. We also observed that host gene ATF2 and miR-933 transcripts share common promoters (Fig. 1a) and are expressed in many similar tissue types (Fig. 1b–d). H3K27ac, an activating histone modification mark, is found to be enriched around the shared promoter of ATF2 and miR-933 (Fig. 1a).
A unique set of hsa-miR-933 targets identified from various target prediction algorithms
To enumerate the targets of hsa-miR-933, we have utilized the targets extracted from various microRNA target databases. From all these databases, we have identified the commonly predicted unique set of twenty-two genes to be targeted by the hsa-miR-933 (Additional file 2).
Host (ATF2) target genes were determined from ChIP-seq experiments
ChIP-seq experiment in the cell line GM12878 from the ENCODE database  provides binding sites of ATF2 protein. For peak-to-target gene annotation, we utilized the well-known target-calling method “closestBed” implemented in BEDTools . A similar algorithm is also utilized by another target-calling program ChIPpeakAnno . We have identified candidate target genes of ATF2 by analyzing its ChIP-seq peaks and the complete dataset of Homo sapiens 149,604 transcripts of protein-coding genes as downloaded from Ensembl (release 70) . By this approach, we have obtained a set of 10,108 targets of ATF2 (Additional file 3).
Common targets of ATF2 and hsa-miR-933 reveal a connection between protein GNAS, PRKACB, and PRKCE
Next, to identify what are the genes that could be regulated by both host (ATF2) and intronic miRNA-933, we sought to determine the common target gene of both intronic miRNA-933 and its host gene ATF2, by Venn diagram overlap analysis. We have found that 15 common genes are targeted significantly by both hsa-miR-933 and ATF2 (χ2 = 37.9891, p value < 0.00001) (Additional file 5). These genes are PEA15, ZFHX4, DAB2IP, CYB5B, DLG1, CAND1, MAP4K4, KPNA1, CUX1, ZNF521, KCMF1, BDNF, GNAS, PRKACB, and PRKCE (Fig. 1e).
These 15 targets were subjected to protein-protein interaction (PPI) network analysis which revealed a connection between GNAS, PRKACB, and PRKCE proteins (Fig. 1f). This information was used to cluster GO annotations and KEGG pathway enrichment analysis.
Intronic miRNA hsa-miR-933 regulates hyperinsulinemia and hyperglycemia in type II diabetes mellitus that may be caused by overexpression of ATF2
In order to illuminate the roles of the genes targeted by intronic miRNA hsa-miR-933, we have performed the enrichment analysis using the target genes for both ATF2 and hsa-miR-933. Significant enrichment (FDR ≤ 0.05) of the pathways necessary for the development of type II diabetes mellitus in terms of GO biological process (Fig. 2a) and KEGG pathways (Fig. 2b) was observed. Functional enrichment analysis showed significant enrichment of biological processes that can disrupt the normal glucose, insulin, and glucagon homeostasis in the cell by positively and negatively regulating them leading to the development of type II diabetes mellitus (Fig. 2a, b). The target genes that were significantly enriched for the development of type II diabetes mellitus in both modules are GNAS, PRKACB, PRKCE, MAP4K4, PEA15, and BDNF (Fig. 2c). The functions of the enriched genes are listed in Additional file 4.
The network analysis using these common target genes shows ATF2-mediated activation of three genes PRKACB, PRKCE, and ATK1. PRKACB activates GNAS which promotes glucagon secretion upon glucose starvation. On the contrary, PRKCE and also MAP4K4 activate ATK1 (Fig. 3). ATK1 inhibits the expression of PEA15, an overexpressed protein in type II diabetes mellitus, where it may contribute to insulin resistance in glucose uptake. Hyperinsulinemia and hyperglycemia are both the pathophysiologies of type II diabetes mellitus (Fig. 4a). From this information, it can be considered for the analysis if hsa-miR-933 can regulate the overexpressed PEA15, downregulated BDNF, or other ATF2 target genes (PRKACB, GNAS, PRKCE, and MAP4K4) to control type II diabetes mellitus (Fig. 4b).
ATF2 transcription factor has been identified as an important component of the insulin signaling system and in maintaining glucose homeostasis in the cell . ATF2 targets genes are involved in insulin action, β cell function, and type II diabetes mellitus. The other functions of ATF2 target genes include adipocyte dysfunction, inhibition of insulin signaling, lipid metabolism, glucose metabolism, β cell dysfunction, and many others [28,29,30]. ATF2 can also regulate gluconeogenesis for glucose production from the liver in response to glucose starvation . Some ATF2 target genes, for example, BDNF and PEA15, are also linked to the development of insulin resistance under glucose stimulation [32, 33]. On the other hand, insulin itself activates ATF2 by phosphorylation of Thr69 and Thr71  by the c-Jun/MAPK pathway. Improper activation of ATF2 target genes and ATF2 itself under conditions of insulin resistance can contribute to the development of type II diabetes mellitus .
In recent times, microRNAs are considered as a possible biomarker or a potential therapeutic for treating type II diabetes mellitus . With this correspondence, the role of intronic miRNA hsa-miR-933’s functional regulation on its host gene ATF2 can provide valuable insight into the progression and control of type II diabetes mellitus.
hsa-miR-933 plays a role in neuronal regeneration and protection by suppressing the overexpressed targets of host gene ATF2 responsible for neuron death
In order to reveal whether the targets of ATF2 and hsa-miR-933 are involved in other significant pathways, we analyzed the functional enrichment using the GOBP module of the target genes. Interestingly, we have observed a significant enrichment of biological processes that contrast between the positive and negative regulatory cycles of neuron apoptotic process (Fig. 5a). Enrichment analysis with the KEGG pathway module also shows the significant pathways necessary for the development of different diseases related to neuron degradation, for example, Huntington’s disease, Parkinson’s disease, and long-term depression (Fig. 5b). The target genes that were significantly enriched in both modules are GNAS, PRKACB, PRKCE, MAP4K4, PEA15, and BDNF (Fig. 5c). The functions of the enriched genes are listed in Additional file 4.
From network analysis, we have tracked down the involvement of ATF2 in the activation of two genes: BDNF via JUN; PRKCE via FOS. PRKCE directly inhibits JUN thus repressing the expression of BDNF leading to neurodegeneration (Fig. 6). On the other hand, recently, it was discovered that hsa-miR-933 is a target of BDNF that may be involved in cell growth, apoptosis, cell proliferation, or the regulation of the cell cycle. It is thus to be analyzed if hsa-miR-933 can act as a potential repressor of neurodegeneration or promote neurodegeneration (Fig. 7).
ATF2 is highly expressed with large variations in the brain (Fig. 1d) and plays a role in both neurodegeneration and neurogenesis. ATF2 seems to play a fundamental role in neuronal viability and in neurological functions in the normal brain. ATF2 is downregulated in the hippocampus and the caudate nucleus in Alzheimer’s, Parkinson’s, and Huntington’s diseases .
ATF2/JUN heterodimers bind and activate CASP3, a key executor of neuronal apoptosis . Following death receptor stimulation, there is phosphorylation and binding of ATF2/JUN to death-inducing ligand promoters (FASLG, TNF, TNFSF10), which allows the spread of death signals . Neuronal apoptosis requires the simultaneous activation of ATF2/JUN and downregulation of FOS . This function is negatively regulated by phosphorylation of ATF2 by PRKCE, which dictates its nuclear localization .
Loss of functional ATF2 leads to hyperphosphorylated JNK and p38 which results in somatic and visceral motor neuron degeneration . On the other hand, activated ATF2 promotes apoptosis of various brain cells, which are cerebellar granule neurons . ATF2/JUN heterodimers also promote the death of sympathetic neurons . Sometimes, ATF2/JUN inhibits JNK-mediated phosphorylation of JUN and protects sympathetic neurons from apoptosis .
Another ATF2 target BDNF can play a neuroprotective role against neural apoptosis . Moreover, recently, brain-derived neurotrophic factor (BDNF) has been identified as a possible target sequence for hsa-miR-933 . So, exploiting this information, it can be analyzed if ATF2 and hsa-miR-933 can play a neuroprotective role in neurodegenerative diseases by regulating their common target gene BDNF.
In recent times, biological modifier microRNA (miRNA) has captured extensive attention for being a potential candidate for studying their role in the regulation of gene expression to delineate a variety of physiological processes related to the progression and pathogenesis of human diseases.
Almost half of the experimentally identified human miRNAs are encoded in the introns of annotated protein-coding genes which can be preferably denoted as intronic miRNA embedded within the host gene. Along with the breakthrough discovery of important functions for these intronic miRNAs, the fact that made research in this field of such significance is intronic miRNAs are co-expressed and share similar transcriptional regulatory mechanisms with their host genes in humans. This parallel expression pattern suggests that these intronic miRNAs may have functions either similar or opposite to that of their host genes in a cell type- or tissue type-specific spatial manner or expression time-specific temporal manner.
Despite the mounting evidence for their importance in normal physiology, very little is known about the regulatory cascade controlled by intronic miRNA and host gene in humans to date. In this study, a genome-wide network-based and data mining approach was adopted to elucidate the functional relationship between intronic miRNAs and host gene, and the effects of intronic miRNAs on host gene targets to control cellular homeostasis.
ATF2 is activated by stress kinases, including JNK and p38. In response to stimuli, ATF2 is phosphorylated on threonine 69 and/or 71 by JNK or by p38. Phosphorylation on Thr69 and Thr71 of ATF2 and its dimerization are required to activate ATF2 transcription factor activity . At the same time, transcriptionally active dimers of ATF2 protein are regulated by ubiquitylation and proteasomal degradation , and phosphorylation of ATF2 on Thr69 and Thr71 promotes its ubiquitylation and degradation . This phosphorylation and dephosphorylation cascade is an important factor for the positive regulation of human insulin gene expression .
Hyperinsulinemia is a condition when the amount of insulin circulating in the blood is higher than expected. Insulin itself can create a positive feedback loop by activating ATF2 in this pathway solely or mediated by MAPK8 (preceded by MAP4K4) or PRKCE which can eventually lead to hyperinsulinemia . Insulin overproduction also leads to the overexpression of other ATF2 target genes that lead to β-pancreatic cell dysfunction—pathophysiology of type II diabetes mellitus .
One of the ATF2 target gene—MAP4K4 is an attenuator for insulin signaling . MAPK4K4 also plays an important role in insulin resistance in response to glucose stimulus for the development of type II diabetes mellitus . GNAS and PRKACB are another two ATF2 target genes that are involved in the regulation of glucagon secretion [51, 52] and insulin secretion for maintaining glucose homeostasis in cells. Impairment in this regulation can lead to hyperglycemia, and the resultant excess glucose in the bloodstream leads to the development of type II diabetes mellitus . Other ATF2 target genes BDNF downregulation and PEA15 upregulation are also linked to the development of insulin resistance under glucose stimulation, pathophysiology of type II diabetes mellitus [32, 33].
All these genes show significant enrichment in the functional enrichment analysis performed in our study with ATF2 and hsa-miR-933. Thus, hsa-miR-933 can be considered as a potential regulator of type II diabetes mellitus. Recent researches are more focused on considering microRNAs as the new therapeutics to treat type II diabetes mellitus . Recently, it has also been discovered that microRNA-30d can induce insulin production by targeting MAP4K4 . microRNAs 103/107 can regulate insulin sensitivity . In this context, it can be considered if hsa-miR-933 can also regulate MAP4K4 or other protein kinases (PRKCE, PRKACB) to control type II diabetes mellitus . Also, recently, brain-derived neurotrophic factor (BDNF) has been identified as a possible target sequence for hsa-miR-933 . So, it can also be considered if hsa-miR-933 can reverse the effect of insulin resistance in response to glucose stimulus by upregulating BDNF.
The significant over-representation of genes associated with type II diabetes mellitus, related pathways, and associated biological processes identified from our study, indicates the firm possibility of considering hsa-miR-933 as a potential regulator of treating the complications linked with the disease.
ATF2 transcription factor is a ubiquitously expressed protein in humans with a more abundant expression in the brain. This profuse expression enables it to have a key role in both neurodegeneration and neurogenesis . This transcription factor also has a function in neuronal migration during development . But overexpression of ATF2 in neuronal-like cell culture promotes nerve cell death . In addition, ATF2 has an essential role in neuronal viability and in neurological functions in the normal brain. In neurodegenerative diseases like Alzheimer’s, Parkinson’s, and Huntington’s diseases, ATF2 is downregulated in the hippocampus and the caudate nucleus . So, the proper regulation of ATF2 is necessary for the neurons in the brain’s central nervous system.
Activating transcription factor 2 (ATF2) is a member of the activator protein-1 family of transcription factors that promote neuronal apoptosis by c-Jun and c-Fos cascade in the cytoplasm. Neuronal apoptosis requires the simultaneous activation of ATF2/c-JUN and downregulation of c-FOS . This function is negatively regulated by phosphorylation of ATF2 by PRKCE, which dictates its nuclear localization . Loss of functional ATF2 leads to hyperphosphorylated JNK and p38 and results in somatic and visceral motor neuron degeneration . Sometimes, activated ATF2 also promotes apoptosis of various brain cells, for example, cerebellar granule neurons .
ATF2/c-JUN heterodimers also promote the death of sympathetic neurons [39, 42]. In a negative feedback loop, ATF2/c-JUN also inhibits JNK-mediated phosphorylation of JUN and protects sympathetic neurons from apoptosis . Additionally, ATF2 can also play a role in neurogenesis by maintaining a subset of neural progenitor cells . Another ATF2 target BDNF has a neuroprotective role against neural apoptosis [43, 59].
Considering this fact, ATF2 can be presented as an important regulator of nervous system viability. The most recent research field, treating neurodegenerative diseases with microRNAs, can be used to modulate any aberrant expression of ATF2 target genes and subsequent abnormality leading to neuronal degradation or to provide sufficient neuronal protection [60,61,62]. As stated earlier, BDNF which is also over-represented in the enrichment analysis of this study of neurodegenerative disease has been identified as a possible target sequence for hsa-miR-933 . So, exploiting this information, it can be analyzed whether ATF2 and hsa-miR-933 can play a neuroprotective role in neurodegenerative diseases by regulating their common target gene BDNF.
Our study of functional enrichment analysis sheds light on the functional interactions between intronic miRNAs and host genes. Therefore, these findings can have potential applications in the development of diagnostic and treatment methods.
The outcome of this study shows that intronic miRNAs and their host genes can coincide with functional relations. Using a data-driven as well as a knowledge-based computational approach, common targets of intronic miRNA and their associated host genes were analyzed for functional co-relations. A further GO analysis predicted an intronic miRNA-host gene interaction network that confirms that the predicted target genes tend to be regulated simultaneously. Taken together, these results indicate either synergistic or antagonistic regulatory effects mediated by either downregulation of genes with an opposed function or fine-tuning of intronic miRNA targets, co-operative to the host gene.
Identification of host gene of hsa-miR-933
Genomic location position data of Homo sapiens microRNA hsa-miR-933 and protein-coding genes were retrieved from Ensembl (release 70)  database using the BioMart portal. We used the IntersectBed feature of the BEDTools  to find the host gene in which the hsa-miR-933 is embedded. We also cross-checked the promoter of the transcripts of hsa-miR-933 using the web portal http://fantom.gsc.riken.jp/5/suppl/De_Rie_et_al_2017/ .
Obtaining the tissue-wide expression profiles of hsa-miR-933 and ATF2
We have retrieved the tissue-wide expression data of the host gene ATF2 from NCBI Bioproject ID: PRJEB4337. Human miRNA tissue atlas  was used to obtain the gene expression data of hsa-miR-933 in different tissue types and bodily fluids.
Identification of the targets of the host gene (ATF2)
After identifying the host gene, we have identified the targets of the host gene. For this, ChIP-seq experiment datasets (Additional file 1) of the host gene were retrieved from the ENCODE database  for analysis. For target calling, we used the “closestBed” feature of the BEDTools  to identify the nearest transcript from the ChIP-seq peak as the target gene using Ensembl (release 70) transcripts.
Identification of the targets of hsa-miR-933
We have extracted the targets of hsa-miR-933 from several miRNA target gene databases which include both experimental and predicted targets, namely TargetScan v6.2 , TarBase v4 , PITA v6 , PicTar , experimentally validated miRTarBase v4.5 , miRBase , and miRecords v4 . A commonly predicted miRNA target gene set was prepared from the targets retrieved from different microRNA target databases for further analysis.
Identification of common targets of intronic miRNA and host gene by overlap analysis
To find out the common targets of hsa-miR-933 and its host gene, we have conducted an overlap analysis using the interactive tool Venny (v 2.0.2) . The result generated the list of targets that are both targeted by the host gene and hsa-miR-933.
Functional enrichment analysis
The functional annotation of target genes is based on Gene Ontology (GO)  as extracted from the Ensembl database (release 70)  and KEGG pathway database . Accordingly, all genes are classified into the ontology categories biological process (GOBP), cellular component (GOCC), molecular function (GOMF), and pathways when possible. The GO/pathway categories that have at least 10 genes annotated were only considered. Gitools  was used for functional over-representation analysis and heatmap generation. Resulting p values were adjusted for multiple testing using the Benjamin and Hochberg’s method of false discovery rate (FDR) .
Clustering of biological processes and pathways with a functional similarity of host genes and intronic miRNA target gene sets
The host gene confers regulatory control by transcriptional or translational inhibition/activation of the hsa-miR-933 target genes or vice versa in possibly related biological processes or pathways. To test this hypothesis for all hosts and target genes, their respective annotations were compared and clustered manually based on the literature. Functional gene annotations as provided by the GO  and KEGG pathway  classify genes according to their molecular function, associated biological processes, or appearance within defined cellular components and pathways they appear.
Construction of network with enriched gene sets derived from functional enrichment analysis
We have constructed the protein-protein interaction networks of enriched targets of the host gene and the miRNA hsa-miR-933 using tools STRING (version 10)  for possibly related biological processes or pathways. Based on protein-protein interaction, and data and knowledge-driven approach, the hypothesis for a particular enriched biological process or pathway was predicted.
The whole workflow of this study is summarized in Additional file 6: Figure S1.
Availability of data and materials
Datasets used in the analyses can be found in various publicly available databases.
Transcription start site
GO biological process
Kyoto Encyclopedia of Genes and Genomes
False discovery rate
Farazi TA, Juranek SA, Tuschl T. The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members. Development. 2008;135(7):1201–14.
Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH, Cuppen E. Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005;120(1):21–4.
Xie X, Lu J, Kulbokas E, Golub TR, Mootha V, Lindblad-Toh K, et al. Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals. Nature. 2005;434(7031):338.
Wu L, Fan J, Belasco JG. MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci. 2006;103(11):4034–9.
Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007;318(5858):1931–4.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120(1):15–20.
Lu M, Zhang Q, Deng M, Miao J, Guo Y, Gao W, et al. An analysis of human microRNA and disease associations. PLoS One. 2008;3(10):e3420.
Joshi PK, Gupta D, Nandal UK, Khan Y, Mukherjee SK, Sanan-Mishra N. Identification of mirtrons in rice using MirtronPred: a tool for predicting plant mirtrons. Genomics. 2012;99(6):370–5.
Griffiths-Jones S, Grocock RJ, Van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(suppl_1):D140-D1D4.
Corcoran DL, Pandit KV, Gordon B, Bhattacharjee A, Kaminski N, Benos PV. Features of mammalian microRNA promoters emerge from polymerase II chromatin immunoprecipitation data. PLoS One. 2009;4(4):e5279.
Kim YK, Kim VN. Processing of intronic microRNAs. EMBO J. 2007;26(3):775–83.
Wang D, Lu M, Miao J, Li T, Wang E, Cui Q. Cepred: predicting the co-expression patterns of the human intronic microRNAs with their host genes. PLoS One. 2009;4(2):e4421.
Ruby JG, Jan CH, Bartel DP. Intronic microRNA precursors that bypass Drosha processing. Nature. 2007;448(7149):83.
Gao X, Qiao Y, Han D, Zhang Y, Ma N. Enemy or partner: relationship between intronic micrornas and their host genes. IUBMB Life. 2012;64(10):835–40.
Truscott M, Islam AB, López-Bigas N, Frolov MV. mir-11 limits the proapoptotic function of its host gene, dE2f1. Genes Dev. 2011;25(17):1820–34.
Truscott M, Islam ABMMK, Lightfoot J, López-Bigas N, Frolov MV. An intronic microRNA links Rb/E2F and EGFR signaling. PLoS Genet. 2014;10(7):e1004493.
Steiman-Shimony A, Shtrikman O, Margalit H. Assessing the functional association of intronic miRNAs with their host genes. RNA. 2018;24(8):991–1004.
Dias IHK, Brown CL, Shabir K, Polidori MC, Griffiths HR. miRNA 933 expression by endothelial cells is increased by 27-hydroxycholesterol and is more prevalent in plasma from dementia patients. J Alzheimers Dis. 2018;64:1009–17.
Xu J, Chen Z, Wang Y, Wang X, Chen L, Yuan T, et al. Several circulating miRNAs related to hyperlipidemia and atherosclerotic cardiovascular diseases. Lipids Health Dis. 2019;18(1):104.
Zhang Y, Ma Y, Xu W, Li W, Min P, Qiu J, et al. Association of microRNA-933 variant with the susceptibility to gastric cancer. J BUON. 2017;22:390–5.
de Rie D, Abugessaisa I, Alam T, Arner E, Arner P, Ashoor H, et al. An integrated expression atlas of miRNAs and their promoters in human and mouse. Nat Biotechnol. 2017;35(9):872–8.
Watson G, ZeA R, Lau E. ATF2, a paradigm of the multifaceted regulation of transcription factors in biology and disease. Pharmacol Res. 2017;119:347–57.
Rosenbloom KR, Dreszer TR, Long JC, Malladi VS, Sloan CA, Raney BJ, et al. ENCODE whole-genome data in the UCSC Genome Browser: update 2012. Nucleic Acids Res. 2012;40(Database issue):D912–7.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.
Zhu LJ, Gazin C, Lawson ND, Pagès H, Lin SM, Lapointe DS, et al. ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data. BMC Bioinformatics. 2010;11(1):237.
Cunningham F, Achuthan P, Akanni W, Allen J, Amode MR, Armean IM, et al. Ensembl 2019. Nucleic Acids Res. 2018;47(D1):D745–D51.
Cao W, Daniel KW, Robidoux J, Puigserver P, Medvedev AV, Bai X, et al. p38 mitogen-activated protein kinase is the central regulator of cyclic AMP-dependent transcription of the brown fat uncoupling protein 1 gene. Mol Cell Biol. 2004;24(7):3057–67.
Baan B, van der Zon GCM, Maassen JA, Ouwens DM. The nuclear appearance of ERK1/2 and p38 determines the sequential induction of ATF2-Thr71 and ATF2-Thr69 phosphorylation by serum in JNK-deficient cells. Mol Cell Endocrinol. 2009;311(1):94–100.
Liu H, Deng X, Shyu Y, Li JJ, Taparowsky E, Hu C-D, Liu H, Deng X, Shyu YJ, Li JJ, Taparowsky EJ, Hu CD. Mutual regulation of c-Jun and ATF2 by transcriptional activation and subcellular localization. EMBO J 25: 1058-1069. EMBO J. 2006;25:1058–69.
Ouwens DM, de Ruiter ND, van der Zon GCM, Carter AP, Schouten J, van der Burgt C, et al. Growth factors can activate ATF2 via a two-step mechanism: phosphorylation of Thr71 through the Ras-MEK-ERK pathway and of Thr69 through RalGDS-Src-p38. EMBO J. 2002;21(14):3782–93.
Cheong J, Coligan JE, Shuman JD. Activating transcription factor-2 regulates phosphoenolpyruvate carboxykinase transcription through a stress-inducible mitogen-activated protein kinase pathway. J Biol Chem. 1998;273(35):22714–8.
Greig FH, Nixon GF. Phosphoprotein enriched in astrocytes (PEA)-15: a potential therapeutic target in multiple disease states. Pharmacol Ther. 2014;143(3):265–74.
Krabbe K, Nielsen A, Krogh-Madsen R, Plomgaard P, Rasmussen P, Erikstrup C, et al. Brain-derived neurotrophic factor (BDNF) and type 2 diabetes. Diabetologia. 2007;50(2):431–8.
Baan B, van Dam H, van der Zon GC, Maassen JA, Ouwens DM. The role of c-Jun N-terminal kinase, p38, and extracellular signal-regulated kinase in insulin-induced Thr69 and Thr71 phosphorylation of activating transcription factor 2. Mol Endocrinol. 2006;20(8):1786–95.
Kumar M, Nath S, Prasad HK, Sharma G, Li Y. MicroRNAs: a new ray of hope for diabetes mellitus. Protein & Cell. 2012;3(10):726–38.
Pearson A, Curtis M, Waldvogel H, Faull R, Dragunow M. Activating transcription factor 2 expression in the adult human brain: association with both neurodegeneration and neurogenesis. Neuroscience. 2005;133(2):437–51.
Song B, Xie B, Wang C, Li M. Caspase-3 is a target gene of c-Jun: ATF2 heterodimers during apoptosis induced by activity deprivation in cerebellar granule neurons. Neurosci Lett. 2011;505(2):76–81.
Herr I, Posovszky C, Di Marzio L, Cifone MG, Boehler T, Debatin K-M. Autoamplification of apoptosis following ligation of CD95-L, TRAIL and TNF-α. Oncogene. 2000;19(37):4255.
Ma C, Ying C, Yuan Z, Song B, Li D, Liu Y, et al. dp5/HRK is a c-Jun target gene and required for apoptosis induced by potassium deprivation in cerebellar granule neurons. J Biol Chem. 2007;282(42):30901–9.
Lopez-Bergami P, Lau E. Ronai Ze. Emerging roles of ATF2 and the dynamic AP1 network in cancer. Nat Rev Cancer. 2010;10(1):65–76.
Ackermann J, Ashton G, Lyons S, James D, Hornung J-P, Jones N, et al. Loss of ATF2 function leads to cranial motoneuron degeneration during embryonic mouse development. PLoS One. 2011;6(4):e19090.
Towers E, Gilley J, Randall R, Hughes R, Kristiansen M, Ham J. The proapoptotic dp5 gene is a direct target of the MLK-JNK-c-Jun pathway in sympathetic neurons. Nucleic Acids Res. 2009;37(9):3044–60.
Mao M, Zhiling W, Hui Z, Shengfu L, Dan Y, Jiping H. Cellular levels of TrkB and MAPK in the neuroprotective role of BDNF for embryonic rat cortical neurons against hypoxia in vitro. Int J Dev Neurosci. 2005;23(6):515–21.
Wu C, Lin H, Wang Q, Chen W, Luo H, Chen W, et al. Discrepant expression of microRNAs in transparent and cataractous human lenses. Invest Ophthalmol Vis Sci. 2012;53(7):3906–12.
Fuchs SY, Ze R. Ubiquitination and degradation of ATF2 are dimerization dependent. Mol Cell Biol. 1999;19(5):3289–98.
Firestein R, Feuerstein N. Association of activating transcription factor 2 (ATF2) with the ubiquitin-conjugating enzyme hUBC9 implication of the ubiquitin/proteasome pathway in regulation of ATF2 in T cells. J Biol Chem. 1998;273(10):5892–902.
Ban N, Yamada Y, Someya Y, Ihara Y, Adachi T, Kubota A, et al. Activating transcription factor-2 is a positive regulator in CaM kinase IV-induced human insulin gene expression. Diabetes. 2000;49(7):1142–8.
Thévenod F. Pathophysiology of diabetes mellitus type 2: roles of obesity, insulin resistance and β-cell dysfunction. Diabetes and Cancer. 19: Karger Publishers; 2008. p. 1-18.
Danai LV, Flach RJR, Virbasius JV, Menendez LG, Jung DY, Kim JH, et al. Inducible deletion of protein kinase Map4k4 in obese mice improves insulin sensitivity in liver and adipose tissues. Mol Cell Biol. 2015;35(13):2356–65.
Bouzakri K, Zierath JR. MAP4K4 gene silencing in human skeletal muscle prevents tumor necrosis factor-α-induced insulin resistance. J Biol Chem. 2007;282(11):7783–9.
Chen D, Zhang A, Fang M, Fang R, Ge J, Jiang Y, et al. Increased methylation at differentially methylated region of GNAS in infants born to gestational diabetes. BMC Med Genet. 2014;15(1):108.
Godoy-Matos AF. The role of glucagon on type 2 diabetes at a glance. Diabetol Metab Syndr. 2014;6(1):91.
Arner P, Kulyte A. MicroRNA regulatory networks in human adipose tissue and obesity. Nat Rev Endocrinol. 2015;11(5):276.
Zhao X, Mohan R, Özcan S, Tang X. MicroRNA-30d induces insulin transcription factor MafA and insulin production by targeting mitogen-activated protein 4 kinase 4 (MAP4K4) in pancreatic β-cells. J Biol Chem. 2012;287(37):31155–64.
Trajkovski M, Hausser J, Soutschek J, Bhat B, Akin A, Zavolan M, et al. MicroRNAs 103 and 107 regulate insulin sensitivity. Nature. 2011;474(7353):649.
Puri V, Virbasius J, Guilherme A, Czech M. RNAi screens reveal novel metabolic regulators: RIP140, MAP4k4 and the lipid droplet associated fat specific protein (FSP) 27. Acta Physiol. 2008;192(1):103–15.
Yu T, Li YJ, Bian AH, Zuo HB, Zhu TW, Ji SX, et al. The regulatory role of activating transcription factor 2 in inflammation. Mediat Inflamm. 2014;2014.
Sanalkumar R, Indulekha CL, Divya TS, Divya MS, Anto RJ, Vinod B, et al. ATF2 maintains a subset of neural progenitors through CBF1/Notch independent Hes-1 expression and synergistically activates the expression of Hes-1 in Notch-dependent neural progenitors. J Neurochem. 2010;113(4):807–18.
Almeida R, Manadas B, Melo C, Gomes J, Mendes C, Graos M, et al. Neuroprotection by BDNF against glutamate-induced apoptotic cell death is mediated by ERK and PI3-kinase pathways. Cell Death Differ. 2005;12(10):1329.
Maciotta Rolandin S, Meregalli M, Torrente Y. The involvement of microRNAs in neurodegenerative diseases. Front Cell Neurosci. 2013;7:265.
Majer A, Boese AS, Booth SA. The role of microRNAs in neurodegenerative diseases: implications for early detection and treatment. Regulatory RNAs: Springer; 2012. p. 443-73.
Wiedemann C. Neurodegenerative disease: microRNAs under threat from LRRK2. Nat Rev Neurosci. 2010;11(9):609.
Ludwig N, Leidinger P, Becker K, Backes C, Fehlmann T, Pallasch C, et al. Distribution of miRNA expression across human tissues. Nucleic Acids Res. 2016;44(8):3865–77.
Rosenbloom KR, Sloan CA, Malladi VS, Dreszer TR, Learned K, Kirkup VM, et al. ENCODE data in the UCSC Genome Browser: year 5 update. Nucleic Acids Res. 2012;41(D1):D56–63.
Vergoulis T, Vlachos IS, Alexiou P, Georgakilas G, Maragkakis M, Reczko M, et al. TarBase 6.0: capturing the exponential growth of miRNA targets with experimental support. Nucleic Acids Res. 2011;40(D1):D222–D9.
Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007;39(10):1278.
Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, et al. Combinatorial microRNA target predictions. Nat Genet. 2005;37(5):495.
Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39(Database issue):D163–9.
Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2018;47(D1):D155–D62.
Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA–target interactions. Nucleic Acids Res. 2008;37(suppl_1):D105-DD10.
Oliveros JC. VENNY. An interactive tool for comparing lists with Venn diagrams. 2015.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2007;36(suppl_1):D480-D4D4.
Perez-Llamas C, Lopez-Bigas N. Gitools: analysis and visualisation of genomic data using interactive heat-maps. PLoS One. 2011;6(5):e19541.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1-2):279–84.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2014;43(D1):D447–D52.
We thank Md. Wahid Murad for his valuable discussion.
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ChIP-seq Experiment data downloaded from ENCODE.
List of genes targeted by hsa-miR-933.
Partial list of 160 targets from a total of 10108 targets of the host gene ATF2.
Functions of the enriched genes for the development of hyperinsulinemia; and Neurodegeneration.
Contingency table of the Chi-square test for finding the statistical significance of the overlapped targets of miR-933 and ATF2.
Whole workflow of the study.
About this article
Cite this article
Islam, A.B.M.M.K., Mohammad, E. & Khan, M.AAK. Aberration of the modulatory functions of intronic microRNA hsa-miR-933 on its host gene ATF2 results in type II diabetes mellitus and neurodegenerative disease development. Hum Genomics 14, 34 (2020). https://doi.org/10.1186/s40246-020-00285-1
- Intronic microRNA
- diabetes mellitus
- Neurodegenerative diseases