Screening differentially expressed genes in salivary exosomes of oral lichen planus by lncRNA microarray

doi:10.3877/cma.j.issn.1674-1366.2022.04.002

Abstract

Abstract:

Objective

To screen the abnormal expression of long non-coding RNAs (lncRNAs) and mRNA by lncRNA expression microarray in oral lichen planus (OLP) salivary exosomes, and to analyze and explore the possible molecular mechanism of lncRNA and mRNA in the occurrence and development of OLP.

Methods

Saliva samples from 9 OLP patients and 3 healthy controls were collected to isolate exosomes. Then, exosomes were detected by nanoparticle tracking analysis (NTA) , transmission electron microscopy and western blot analysis of exosome-specific biomarkers. The expression profiles of lncRNA and mRNA in the salivary exosomes of erosive OLP patients (EOLP group) and reticulated OLP patients (ROLP group) were compared with those of healthy controls by lncRNA expression microarray, and the differentially expressed genes were screened. The differential genes were analyzed by Gene Ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment analysis.

Results

NTA, transmission electron microscopy and Western blot confirmed the isolation of exosomes. Compared with the control group, there were 267 differentially expressed lncRNAs in the EOLP group, including 99 up-regulated and 168 down-regulated lncRNAs, and there were 122 differentially expressed mRNAs, of which 38 were up-regulated and 84 were down-regulated. There were 201 differentially expressed lncRNAs in the ROLP group, of which 83 were up-regulated and 118 were down-regulated, and there were 86 differentially expressed mRNAs, of which 32 were up-regulated and 54 were down-regulated. There were 50 identical differentially expressed mRNAs and 128 identical differentially expressed lncRNAs in the two groups. GO and KEGG analysis showed that differentially expressed genes involved in gene transcription, protein translation, immune response and other biological processes.

Conclusions

In this study, we determined the expression profile of lncRNA and mRNA in salivary exosomes of OLP patients, and identified differentially expressed lncRNAs and mRNAs associated with OLP. These lncRNAs and mRNAs may serve as important candidates for diagnosis and elucidate the pathogenesis of OLP.

Key words: Lichen planus, oral, Saliva, Exosomes, Long Noncoding RNA, RNA, messenger

Wenxin Zuo, Li Yuan, Tianhui Yang, Jianming Tang, Zhiyi Zhou, Fei He. Screening differentially expressed genes in salivary exosomes of oral lichen planus by lncRNA microarray[J]. Chinese Journal of Stomatological Research(Electronic Edition), 2022, 16(04): 208-218.

/ / Recommend

Add to citation manager EndNote|Ris|BibTeX

URL: https://zhkqyxyjzz.cma-cmc.com.cn/EN/10.3877/cma.j.issn.1674-1366.2022.04.002

https://zhkqyxyjzz.cma-cmc.com.cn/EN/Y2022/V16/I04/208

Figures/Tables 17

References 34

[1]	Tampa M，Caruntu C，Mitran M，et al. Markers of oral lichen planus malignant transformation[J]. Dis Markers，2018，2018：1959506. DOI：10.1155/2018/1959506.
[2]	González-Moles MÁ，Warnakulasuriya S，González-Ruiz I，et al. Worldwide prevalence of oral lichen planus：A systematic review and meta-analysis[J]. Oral Dis，2021，27（4）：813-828. DOI：10.1111/odi.13323.
[3]	Ghazi N，Khorasanchi M. Markers associated with malignant transformation of oral lichen planus：A review article[J]. Arch Oral Biol，2021，127：105158. DOI：10.1016/j.archoralbio.2021.105158.
[4]	Gupta S，Jawanda MK. Oral lichen planus：An update on etiology，pathogenesis，clinical presentation，diagnosis and management [J]. Indian J Dermatol，2015，60（3）：222-229. DOI：10.4103/0019-5154.156315.
[5]	Chiang CP，Chang JYF，Wang YP，et al. Oral lichen planus—Differential diagnoses，serum autoantibodies，hematinic deficiencies，and management[J]. J Formos Med Assoc，2018，117（9）：756-765. DOI：10.1016/j.jfma.2018.01.021.
[6]	Nosratzehi T. Oral lichen planus：An overview of potential risk factors，biomarkers and treatments[J]. Asian Pac J Cancer Prev，2018，19（5）：1161-1167. DOI：10.22034/APJCP.2018.19.5.1161.
[7]	Han Y，Jia L，Zheng Y，et al. Salivary exosomes：Emerging roles in systemic disease[J]. Int J Biol Sci，2018，14（6）：633-643. DOI：10.7150/ijbs.25018.
[8]	Cheshmi B，Cheshomi H. Salivary exosomes：Properties，medical applications，and isolation methods[J]. Mol Biol Rep，2020，47（8）：6295-6307. DOI：10.1007/s11033-020-05659-1.
[9]	Kowalczyk MS，Higgs DR，Gingeras TR. Molecular biology：RNA discrimination[J]. Nature，2012，482（7385）：310-311. DOI：10.1038/482310a.
[10]	Caley DP，Pink RC，Trujillano D，et al. Long noncoding RNAs，chromatin，and development[J]. Scientific World Journal，2010，10：90-102. DOI：10.1100/tsw.2010.7.
[11]	Bridges MC，Daulagala AC，Kourtidis A. LNCcation：lncRNA localization and function[J]. J Cell Biol，2021，220（2）：e202009045. DOI：10.1083/jcb.202009045.
[12]	Aune TM，Spurlock CF 3rd. Long non-coding RNAs in innate and adaptive immunity[J]. Virus Res，2016，212：146-160. DOI：10.1016/j.virusres.2015.07.003.
[13]	Ritchie ME，Phipson B，Wu D，et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies [J]. Nucleic Acids Res，2015，43（7）：e47. DOI：10.1093/nar/gkv007.
[14]	Alexa A，Rahnenfuhrer J，Depends R，et al. topGO：Enrichment analysis for gene ontology.Version 2.10.0[EB/OL]. [2010-10-01]. [2022-08-24]. URL
[15]	Yu G，Wang LG，Han Y，et al. clusterProfiler：An R package for comparing biological themes among gene clusters[J]. OMICS，2012，16（5）：284-287. DOI：10.1089/omi.2011.0118.
[16]	Zhou G，Zhang J，Ren XW，et al. Increased B7-H1 expression on peripheral blood T cells in oral lichen planus correlated with disease severity[J]. J Clin Immunol，2012，32（4）：794-801. DOI：10.1007/s10875-012-9683-2.
[17]	Li P，Kaslan M，Lee SH，et al. Progress in exosome isolation techniques[J]. Theranostics，2017，7（3）：789-804. DOI：10.7150/thno.18133.
[18]	Deng Y，Cao Y，Wang L，et al. The role and application of salivary exosomes in malignant neoplasms[J]. Cancer Manag Res，2021，13：5813-5820. DOI：10.2147/CMAR.S321225.
[19]	Zhan C，Yang X，Yin X，et al. Exosomes and other extracellular vesicles in oral and salivary gland cancers[J]. Oral Dis，2020，26（5）：865-875. DOI：10.1111/odi.13172.
[20]	Byun JS，Hong SH，Choi JK，et al. Diagnostic profiling of salivary exosomal microRNAs in oral lichen planus patients[J]. Oral Dis，2015，21（8）：987-993. DOI：10.1111/odi.12374.
[21]	Chen Y，Li Z，Chen X，et al. Long non-coding RNAs：From disease code to drug role[J]. Acta Pharm Sin B，2021，11（2）：340-354. DOI：10.1016/j.apsb.2020.10.001.
[22]	Wang J，Zhai X，Guo J，et al. Long non-coding RNA DQ786243 modulates the induction and function of CD4⁺ Treg cells through Foxp3-miR-146a-NF-κB axis：Implications for alleviating oral lichen planus[J]. Int Immunopharmacol，2019，75：105761. DOI：10.1016/j.intimp.2019.105761.
[23]	Salmena L，Poliseno L，Tay Y，et al. A ceRNA hypothesis：The Rosetta Stone of a hidden RNA language？[J]. Cell，2011，146（3）：353-358. DOI：10.1016/j.cell.2011.07.014.
[24]	Wang L，Cho KB，Li Y，et al. Long noncoding RNA（lncRNA）-mediated competing endogenous RNA networks provide novel potential biomarkers and therapeutic targets for colorectal cancer [J]. Int J Mol Sci，2019，20（22）：5758. DOI：10.3390/ijms20225758.
[25]	Giraldez MD，Spengler RM，Etheridge A，et al. Phospho-RNA-seq：A modified small RNA-seq method that reveals circulating mRNA and lncRNA fragments as potential biomarkers in human plasma[J]. EMBO J，2019，38（11）：e101695. DOI：10.15252/embj.2019101695.
[26]	Zhao H，Gao Y，Chen Q，et al. RAD51AP1 promotes progression of ovarian cancer via TGF-β/Smad signalling pathway[J]. J Cell Mol Med，2021，25（4）：1927-1938. DOI：10.1111/jcmm.15877.
[27]	Wang H，Dong H，Qiao L，et al. ZEB1 induces non-small cell lung cancer development by targeting microRNA-320a to increase the expression of RAD51AP1[J]. Exp Cell Res，2021，405（2）：112687. DOI：10.1016/j.yexcr.2021.112687.
[28]	Wu Y，Wang H，Qiao L，et al. Silencing of RAD51AP1 suppresses epithelial-mesenchymal transition and metastasis in non-small cell lung cancer[J]. Thorac Cancer，2019，10（9）：1748-1763. DOI：10.1111/1759-7714.13124.
[29]	Jing X，Chen Y，Chen Y，et al. Down-regulation of USP8 inhibits cholangiocarcinoma cell proliferation and invasion[J]. Cancer Manag Res，2020，12：2185-2194. DOI：10.2147/CMAR.S234586.
[30]	Liu YF，Yang A，Liu W，et al. NME2 reduces proliferation，migration and invasion of gastric cancer cells to limit metastasis[J]. PLoS One，2015，10（2）：e0115968. DOI：10.1371/journal.pone.0115968.
[31]	Kayama K，Watanabe S，Takafuji T，et al. GRWD1 negatively regulates p53 via the RPL11-MDM2 pathway and promotes tumorigenesis[J]. EMBO Rep，2017，18（1）：123-137. DOI：10.15252/embr.201642444.
[32]	Giannetti L，Dello Diago AM，Spinas E. Oral lichen planus[J]. J Biol Regul Homeost Agents，2018，32（2）：391-395.
[33]	Roopashree MR，Gondhalekar RV，Shashikanth MC，et al. Pathogenesis of oral lichen planus：A review[J]. J Oral Pathol Med，2010，39（10）：729-734. DOI：10.1111/j.1600-0714.2010.00946.x.
[34]	Bombeccari GP，Guzzi G，Tettamanti M，et al. Oral lichen planus and malignant transformation：A longitudinal cohort study [J]. Oral Surg Oral Med Oral Pathol Oral Radiol Endod，2011，112（3）：328-334. DOI：10.1016/j.tripleo.2011.04.009.

名称	log₂FC	P值	表达情况
lnc-EVX1-2	1.662 166	0.000 024 2	上调
lnc-NHS-2	1.783 052	0.000 041 4	上调
lnc-GEMIN4-4	2.095 217	0.000 073 8	上调
linc-PDE6B	1.453 091	0.000 111	上调
lnc-NIPAL2-2	1.017 79	0.000 186	上调
lnc-HRH4-11	1.593 248	0.000 186	上调
lnc-ZMYM6NB-2	2.224 494	0.000 189	上调
lnc-TRA2B-2	1.333 532	0.000 195	上调
lnc-WNK1-4	2.782 025	0.000 211	上调
lnc-C12orf4-2	1.222 968	0.000 226	上调
lnc-RPS24-7	-2.682 63	0.000 000 556	下调
lnc-AL078585.1-14	-2.045 71	0.000 001 19	下调
lnc-PRKAA1-1	-2.280 28	0.000 002 82	下调
lnc-CYP26B1-4	-2.518 64	0.000 003 68	下调
lnc-COL9A1-2	-1.911 62	0.000 003 99	下调
lnc-USP12-7	-1.933 52	0.000 004 78	下调
lnc-C14orf39-2	-1.760 29	0.000 004 99	下调
lnc-CHM-7	-2.327 78	0.000 006 78	下调
lnc-C2orf43-5	-2.209	0.000 007 59	下调
lnc-RCOR3-1	-2.198 98	0.000 007 98	下调

名称	log₂FC	P值	表达情况
lnc-EVX1-2	1.662 166	0.000 024 2	上调
lnc-NHS-2	1.783 052	0.000 041 4	上调
lnc-GEMIN4-4	2.095 217	0.000 073 8	上调
linc-PDE6B	1.453 091	0.000 111	上调
lnc-NIPAL2-2	1.017 79	0.000 186	上调
lnc-HRH4-11	1.593 248	0.000 186	上调
lnc-ZMYM6NB-2	2.224 494	0.000 189	上调
lnc-TRA2B-2	1.333 532	0.000 195	上调
lnc-WNK1-4	2.782 025	0.000 211	上调
lnc-C12orf4-2	1.222 968	0.000 226	上调
lnc-RPS24-7	-2.682 63	0.000 000 556	下调
lnc-AL078585.1-14	-2.045 71	0.000 001 19	下调
lnc-PRKAA1-1	-2.280 28	0.000 002 82	下调
lnc-CYP26B1-4	-2.518 64	0.000 003 68	下调
lnc-COL9A1-2	-1.911 62	0.000 003 99	下调
lnc-USP12-7	-1.933 52	0.000 004 78	下调
lnc-C14orf39-2	-1.760 29	0.000 004 99	下调
lnc-CHM-7	-2.327 78	0.000 006 78	下调
lnc-C2orf43-5	-2.209	0.000 007 59	下调
lnc-RCOR3-1	-2.198 98	0.000 007 98	下调

名称	log₂FC	P值	表达情况
FAM72A	1.477 287	0.000 107	上调
CLDN2	1.070 014	0.000 124	上调
SNORA79	1.844 52	0.000 334	上调
LOC101060385	2.476 236	0.000 341	上调
USP8	1.261 288	0.000 383	上调
MB	1.027 679	0.000 398	上调
DLD	1.003 385	0.000 441	上调
NARF	1.365 735	0.000 454	上调
PCDHB6	2.547 105	0.000 64	上调
RAD51AP1	1.710 067	0.000 797	上调
RPL11	-2.282 32	0.000 002 47	下调
ATP5PO	-1.747 77	0.000 003 24	下调
RPS12	-2.036 69	0.000 005 61	下调
RPS17	-1.718 66	0.000 022	下调
RPL8	-2.056 28	0.000 026	下调
RPL37A	-2.029 03	0.000 027 5	下调
NME2	-1.701 05	0.000 027 7	下调
VTRNA1-1	-2.335 26	0.000 032 4	下调
RPL27	-2.108 27	0.000 033 2	下调
RPS25	-2.299 16	0.000 033 3	下调

名称	log₂FC	P值	表达情况
FAM72A	1.477 287	0.000 107	上调
CLDN2	1.070 014	0.000 124	上调
SNORA79	1.844 52	0.000 334	上调
LOC101060385	2.476 236	0.000 341	上调
USP8	1.261 288	0.000 383	上调
MB	1.027 679	0.000 398	上调
DLD	1.003 385	0.000 441	上调
NARF	1.365 735	0.000 454	上调
PCDHB6	2.547 105	0.000 64	上调
RAD51AP1	1.710 067	0.000 797	上调
RPL11	-2.282 32	0.000 002 47	下调
ATP5PO	-1.747 77	0.000 003 24	下调
RPS12	-2.036 69	0.000 005 61	下调
RPS17	-1.718 66	0.000 022	下调
RPL8	-2.056 28	0.000 026	下调
RPL37A	-2.029 03	0.000 027 5	下调
NME2	-1.701 05	0.000 027 7	下调
VTRNA1-1	-2.335 26	0.000 032 4	下调
RPL27	-2.108 27	0.000 033 2	下调
RPS25	-2.299 16	0.000 033 3	下调

名称	log₂FC	P值	类型
lnc-C1orf124-3	2.689 691	0.000 034 6	上调
lnc-NHS-2	1.988 245	0.000 059 3	上调
lnc-FAM160A1-3	1.233 302	0.000 108	上调
AC012451.1	3.555 147	0.000 154	上调
lnc-ZMYM6NB-2	2.584 77	0.000 179	上调
lnc-NIPAL2-2	1.143 354	0.000 24	上调
lnc-EVX1-2	1.503 094	0.000 273	上调
lnc-GEMIN4-4	2.109 394	0.000 275	上调
lnc-HRH4-11	1.658 174	0.000 479	上调
lnc-WNK1-4	2.913 166	0.000 51	上调
lnc-RPS24-7	-2.737 15	0.000 002 43	下调
lnc-TRIM52-3	-2.638 4	0.000 012 1	下调
lnc-USP12-7	-1.887 67	0.000 030 1	下调
lnc-PRKAA1-1	-2.122 5	0.000 030 7	下调
lnc-AL078585.1-14	-1.768 69	0.000 030 8	下调
lnc-RPL8-1	-2.753 47	0.000 040 2	下调
lnc-C2orf43-5	-2.147 78	0.000 048 1	下调
lnc-RCOR3-1	-2.135 04	0.000 051	下调
lnc-HSPB3-2	-1.738 52	0.000 054 7	下调
lnc-CHCHD7-5	-1.657 17	0.000 064 3	下调

Please choose a citation manager

Content to export