MIR181A2HG was down-regulated in psoriatic lesions
The expression of MIR181A2HG in psoriatic lesions and normal skins was explored using public available GEO datasets (GSE13355). MIR181A2HG was found to be down-regulated in psoriatic lesions (Figure 1A). qRT-PCR was performed to measure the expression level of MIR181A2HG in clinical specimens. Consistently, the results revealed that MIR181A2HG was lower in psoriatic lesions than in normal skin tissues (Figure 1B). GTEx database result showed that MIR181A2HG was highly expressed in skin tissues (Figure 1C). Taken together, these data suggested that MIR181A2HG may be involved in the pathogenesis of psoriasis.
Enhanced expression of MIR181A2HG inhibited keratinocytes proliferation
In order to determine the possible effects of MIR181A2HG on keratinocyte proliferation, we transfected the overexpression vector of MIR181A2HG into HaCaT keratinocytes. qRT-PCR assay revealed that MIR181A2HG overexpression plasmid transfection significantly increased MIR181A2HG level (Figure 2A). CCK8 assay revealed that MIR181A2HG up-regulation dramatically suppressed HaCaT keratinocytes growth in 72 and 96 hours (Figure 2B). Keratin 6 (KRT6) and keratin 16 (KRT16) were thought to be the hallmarks of psoriatic keratinocytes hyperproliferation [14–16]. qRT-PCR assay indicated that MIR181A2HG up-regulation resulted in significant reduction of the expression of KRT6 and KRT16 (Figure 2C). Consistent with the results of CCK8 assays, EdU incorporation rate was reduced following MIR181A2HG overexpression (Figure 2D). Meanwhile, we enhanced MIR181A2HG expression via a lentiviral vector in primary normal human epidermal keratinocytes (NHEKs) (Figure 2E). CCK8 assay revealed that MIR181A2HG up-regulation inhibited NHEKs growth (Figure 2F) and decreased EdU incorporation rate (Figure 2H). The expression of KRT6 and KRT16 was also reduced (Figure 2G). Collectively, these results suggested that enhanced expression of MIR181A2HG inhibited proliferation of keratinocytes.
Knockdown of MIR181A2HG promoted keratinocytes proliferation
We further explored the impact of MIR181A2HG knockdown on the proliferation of keratinocytes. qRT-PCR assay showed that MIR181A2HG expression was significantly suppressed by specific siRNA (Figure 3A and 3E). CCK8 assay revealed that MIR181A2HG knockdown could promote HaCaT keratinocytes and NHEKs growth (Figure 3B and 3F) and increased EdU incorporation rate (Figure 3D and 3H). qRT-PCR assay showed that MIR181A2HG knockdown resulted in significant up-regulation of KRT6 and KRT16 expression in HaCaT keratinocytes and NHEKs (Figure 3C and 3G). These results showed that MIR181A2HG knockdown promoted proliferation of keratinocytes.
As the host gene of miR-181a2/b2, the intron region of MIR181A2HG produces miR-181a/b (Figure 4A). Previous studies have shown that miR-181a/b plays an important role in the regulation of cell proliferation [17–20]. We found that miR-181a and miR-181b were down-regulated in psoriatic lesions based on GEO dataset (GSE145054) (Figure 4B). In order to explore whether the role of MIR181A2HG in regulating keratinocytes proliferation depends on miR-181a/b, we knocked down MIR181A2HG and detected the effect on expression of miR-181a/b. The results showed that knockdown of MIR181A2HG could not affect the expression of miR-181a and miR-181b (Figure 4C), suggesting that MIR181A2HG may regulate keratinocytes proliferation in a miR-181a/b-independent manner.
MIR181A2HG could interact with SRSF1 in keratinocytes
LncATLAS database analysis showed that the cellular localization of MIR181A2HG in different cells was not identical and showed cellular specificity (Figure 5A). Nuclear/cytoplasmic separation qRT-PCR result showed that MIR181A2HG transcript was localized both in the nucleus and cytoplasm, and slightly higher in the nucleus (Figure. 5B), implying that MIR181A2HG may play a regulatory role at pre-transcriptional, transcriptional and post-transcriptional level to regulate proliferation of keratinocytes.
As important regulatory RNAs, lncRNAs can play a variety of regulatory roles by binding to proteins. In order to screen the proteins interacted with MIR181A2HG, RNA pulldown-MS was used to identify the proteins that may bind to MIR181A2HG. Compared with the antisense strand probe, MIR181A2HG probe could pull down multiple proteins. 356 proteins were finally identified by MS (Figure 6A). KEGG pathway enrichment analysis showed that these proteins were mainly related to spliceosome, RNA transport, protein processing (Figure 6B). GO biological process analysis showed that these proteins were mainly involved in RNA processing and transport, mRNA metabolism (Figure 6C). These results suggested that MIR181A2HG may interact with RNA-binding proteins (RBPs) to regulate RNA processing, thereby affecting cell proliferation. We used Cytoscape software to construct the interaction network of these proteins (Figure 6D). Among them, 13 genes were considered as key genes (degree>50) (Figure 6E). The detailed information of these genes was shown in Table 1. In order to explore the possible RBPs interacting with MIR181A2HG, we took the intersection of these 13 key proteins with the interaction proteins predicted in catRAPID and ENCORI databases. The result showed that NOP56 and SRSF1 may be RBPs interacting with MIR181A2HG (Figure 6F).
Table 1
Detailed information of 13 key genes.
Gene symbols | Degrees | Full names | Gene function |
NOP56 | 54 | nucleolarprotein 56 | Nop56p is required for assembly of the 60S ribosomal subunit and is involved in pre-rRNA processing. The protein encoded by this gene is similar in sequence to Nop56p and is also found in the nucleolus. Expansion of a GGCCTG repeat from 3-8 copies to 1500-2500 copies in an intron of this gene results in spinocerebellar ataxia 36. |
SRSF1 | 53 | Serine and arginine rich splicing factor 1 | This gene encodes a member of the arginine/serine-rich splicing factor protein family. The encoded protein can either activate or repress splicing, depending on its phosphorylation state and its interaction partners. |
SNRPD2 | 63 | Small nuclear ribonucleoprotein D2 polypeptide | The protein encoded by this gene belongs to the small nuclear ribonucleoprotein core protein family. It is required for pre-mRNA splicing and small nuclear ribonucleoprotein biogenesis. |
POLR2A | 60 | RNA polymerase II subunit A | This gene encodes the largest subunit of RNA polymerase II, the polymerase responsible for synthesizing messenger RNA in eukaryotes. This subunit, in combination with several other polymerase subunits, forms the DNA binding domain of the polymerase, a groove in which the DNA template is transcribed into RNA. |
FBL | 57 | Fibrillarin | This gene product is a component of a nucleolar small nuclear ribonucleoprotein (snRNP) particle thought to participate in the first step in processing preribosomal RNA. It is associated with the U3, U8, and U13 small nuclear RNAs and is located in the dense fibrillar component (DFC) of the nucleolus. |
DDX39B | 56 | DExD-box helicase 39B | This gene encodes a member of the DEAD box family of RNA-dependent ATPases that mediate ATP hydrolysis during pre-mRNA splicing. The encoded protein is an essential splicing factor required for association of U2 small nuclear ribonucleoprotein with pre-mRNA, and it also plays an important role in mRNA export from the nucleus to the cytoplasm. |
SNRPD3 | 55 | Small nuclear ribonucleoprotein D3 polypeptide | This gene encodes a core component of the spliceosome, which is a nuclear ribonucleoprotein complex that functions in pre-mRNA splicing. |
RPL3 | 54 | Fibosomal protein L3 | This gene encodes a ribosomal protein that is a component of the 60S subunit. The protein belongs to the L3P family of ribosomal proteins and it is located in the cytoplasm. |
RPL5 | 53 | Fibosomal protein L5 | This gene encodes a member of the L18P family of ribosomal proteins and component of the 60S subunit. |
GNB2L1 | 53 | Guanine nucleotide binding protein beta polypeptide 2-like 1 | Component of the 40S ribosomal subunit involved in translational repression. |
EEF1A1 | 52 | Eukaryotic translation elongation factor 1 alpha 1 | This gene encodes an isoform of the alpha subunit of the elongation factor-1 complex, which is responsible for the enzymatic delivery of aminoacyl tRNAs to the ribosome. |
EIF4A1 | 51 | Eukaryotic translation initiation factor 4A1 | ATP-dependent RNA helicase which is a subunit of the eIF4F complex involved in cap recognition and is required for mRNA binding to ribosome. |
RPS9 | 50 | Ribosomal protein S9 | This gene encodes a ribosomal protein that is a component of the 40S subunit. The protein belongs to the S4P family of ribosomal proteins. |
In order to confirm whether MIR181A2HG interacts with NOP56 and SRSF1, RNA pulldown-Western blotting assay was performed. The results indicated that MIR181A2HG probe could bind to SRSF1 but not NOP56 (Figure 7A). RIP-qRT-PCR was carried out to further confirm the interaction between MIR181A2HG and SRSF1. The results showed that MIR181A2HG was enriched in SRSF1 but not GAPDH immunoprecipitates in NHEKs (Figure 7B). Notably, the expression level of SRSF1 was found to be higher in the psoriatic lesions than in normal skin tissues in GEO datasets (GSE13355) (Figure 7C). Interestingly, SRSF1 knockdown inhibited the expression of KRT16, which could be reversed by knockdown of MIR181A2HG (Figure 7D). Taken together, these results indicated MIR181A2HG may negatively regulate keratinocytes proliferation, at least in part, by interacting with SRSF1 (Figure 7E).