13 July 2020: Database Analysis
Prediction of Specific Subtypes and Common Markers of Non-Small Cell Lung Cancer Based on Competing Endogenous RNA Network
Yao Liu1ABCEF, Hao Wang2CDEF, Wenhan Yang3AFG, Youhui Qian1AFG*DOI: 10.12659/MSM.922280
Med Sci Monit 2020; 26:e922280
Abstract
BACKGROUND: There are various pathological types of lung cancer, including squamous cell carcinoma and adenocarcinoma. Although both of them are lung cancers, there are significant differences in diagnosis, pathogenesis, location, imaging, metastasis, and treatment. According to the competing endogenous RNA (ceRNA) theory, long non-coding RNAs (lncRNAs) compete with encoding protein genes (mRNAs) to connect with miRNAs, thus affecting the level of mRNA.d
MATERIAL AND METHODS: First, using the t test, we identified mRNAs and lncRNAs that have different expressions (fold change >2, P<0.01) in normal samples and in tumor samples. We calculated the significance of the shared miRNAs for mRNAs and lncRNAs by hypergeometric test (P<0.01). Further, mRNA and lncRNA pairs with co-expression relationships in cancer samples were used to establish ceRNA networks. Then, the random walk algorithm was used to optimize the specific ceRNA networks and identify potential prognostic markers of survival. Finally, we built a common ceRNA network to find markers of non-small-cell lung cancer.
RESULTS: We identified some potential key markers, such as PVT1, LINC00472, CDKN2A, and FAM83B, in LUSC and HOXA11-AS, HNF1A-AS1, LINC00511, and HOTAIR in LUAD by analyzing the ceRNA networks. Moreover, a number of common ceRNA pairs, such as CDC25C/CDK1/RRM2-LINC00355, have been found, and they are also significant markers for tumor survival and prognosis.
CONCLUSIONS: In summary, the present study provides a comparative analysis in 2 kinds of lung cancer ceRNA networks. Some specific and common markers we predicted that may be of great importance for clinical diagnosis and treatment.
Keywords: Neoplasms, Squamous Cell, RNA Stability, Tumor Markers, Biological, Carcinoma, Squamous Cell, Databases, Genetic, Gene Regulatory Networks, RNA, Messenger
Background
Lung cancer is a common cancer worldwide, with high morbidity and mortality [1]; 85% of lung cancer is non-small-cell lung cancer (NSCLC) [2,3]. Although significant progress has been made in treatment, such as chemotherapy and radiation therapy, the 5-year survival rate for NSCLC is still 15% [4]. Therefore, it is very important to understand the molecular mechanism of non-small cell lung cancer and to seek new therapeutic strategies. In addition, regional or remote transfers result in higher mortality [5]. Long non-coding RNAs (lncRNA) are a class of single-stranded RNA molecules with a length of more than 200 nt, located in the nucleus or cytoplasm and lacking protein-coding functions [6,7]. The expression of genes is regulated at many levels in the form of RNA (e.g., epigenetics, transcriptional regulation, and post-transcriptional regulation) [8,9]. lncRNA is only expressed at specific developmental stages, and it has tissue- or cell-specificity. Research on CeRNA (competing endogenous RNAs) and internal competitive RNAs reveals a new mechanism of interaction between RNAs. It is known that microRNAs can cause gene silencing by binding mRNAs, while ceRNAs can regulate gene expression by competitively competitive binding to microRNAs [5]. Recent studies have reported that lncRNAs are associated with other RNA transcripts by competing for internal RNAs (ceRNAs) mechanisms [10–12]. mRNAs compete with lncRNAs for miRNA molecules to alleviate the inhibition of miRNA targets [13–15]. This type of lncRNA-related competitive triplet is a subclass of ceRNA and occurs widely in humans and several other species [16–18]. Therefore, this study constructed and analyzed the specific ceRNA network of lung adenocarcinoma and squamous cell carcinoma. By constructing the specific ceRNA network and universal network of lung cancer, we found some important markers of cancer that could be used as a guide and reference for the further study of lung cancer.
Material and Methods
IDENTIFICATION OF DIFFERENTIAL EXPRESSION GENES:
We downloaded the gene expression data (FPKM) of LUAD and LUSC from TCGA project. First, we used PCA method to extract the classification feature information of 2 kinds of cancer expression profile data, and greatly reduced the dimension of expression profile data and visualized the classification features while maintaining high classification accuracy. Differentially expressed genes of LUAD and LUSC were obtained by unpaired
BUILDING THE SPECIFIC DE-CERNA NETWORK:
The miRNA targets were obtained from StarBase2.0, TargetScan, mirTarbase, LncBase, and MiRcode databases. mRNA and lncRNA pairs with more than zero shared miRNAs and a significant hypergeometric test P value (P<0.01) were retained. In total, we integrated more than human 800 000 lncRNA-PCG interactions. Further, we used expression profiles of lung cancer to calculate the expression correlations of the differential mRNA-lncRNA pairs by Pearson correlation analysis. Finally, we defined the consistent expression pairs as potential ceRNA relationships (cor >0; P<0.01). We used Cytoscape software to analyze the topological properties of the network.
OPTIMIZING THE DE-CERNA NETWORK AND BUILDING THE SUB-NET:
The random walk (RW) algorithm, which is based on a network
where r represents the weight of a node in each migration with the probability assigned to the neighbor node, where the default value of 0.7, w is used as a standardized adjacency matrix representing the network, and p 0 means the initial weight of the node. The weight vector, p t, represents the new weight of the node after t-time migration of the network. We obtained gene-disease associations from DisGeNET, Genetic Association Database, Lnc2Cancer, and LncRNADisease. Then, we mapped these known pathogenic genes of lung cancer to the ceRNA network as seed nodes. We used the random walk algorithm according to the topological structure of the ceRNA network to score each node. Finally, we selected the genes with the top random walk score (10 times the number of seed nodes for each cancer) as the nodes to extract the sub-net from the ceRNA network. We used Cytoscape software to visualize the sub-network.
BUILDING A COMMON CERNA NETWORK BETWEEN LUAD AND LUSC:
To obtain the common markers of lung cancer, we extracted the same ceRNA relationship between squamous cell carcinoma and adenocarcinoma and we defined this as a common ceRNA network. We extracted the simultaneous lncRNA-mRNA pairs in the 2 cancers to build a common CeRNA network. Finally, we input the genes list of each module into the DAVID (
SURVIVAL ANALYSIS:
We downloaded gene expression and clinical data from the TCGA project. The relationship between gene expression and prognosis of patients was assessed by Kaplan-Meier analysis. The mean gene expression was used as the cut-off to classify patients into high- and low-risk groups. For a gene in the ceRNA network, we calculated an average expression to evaluate the survival relation of the gene using the log-rank test (P<0.05). We used R package (“Survival”) to evaluate the differences between different groups in the survival analysis experiments.
Results
THE FUNCTIONS OF DIFFERENTIALLY EXPRESSED GENES:
We used PCA to reduce the expression profile of the 2 cancers (LUAD, LUSC) and showed the classification characteristics, finding that both sets of data maintain high classification accuracy at the same time (Figure 1). After the data quality evaluation, we used the non-paired t test method to find genes that are differentially expressed (|Log 2 FoldChange| >1 and P<0.01) between normal and cancer samples (Figure 2A–2D). The statistical significance P_value from low to high sequencing of 1000 differentially expressed genes have been annotated on some important functional GO terms cell multiplication (Figure 2E, 2F) such as “cell division”, “DNA replication”, “G1/S transition of mitotic cell cycle”, “regulation of cell cycle”, “sister chromatid cohesion”, “spindle organization”, “mitotic sister chromatid segregation”, “G2/M transition of mitotic cell cycle”, and pathways associated with cancer such as “cell cycle”. These differences in mRNAs and lncRNAs were used in the following analyses.
LUNG CANCER-SPECIFIC DE-CERNA NETWORK ANALYSIS:
We used the shared miRNA and expression correlation between mRNAs and lncRNAs to build the specific DE-ceRNA network of 2 kinds of lung cancer (see Materials and Methods section above). To retain the highly expressed genes, we deleted the nodes in the network with an average expression of less than 0.1 in all samples. The LUSC DE-ceRNA network contains 16 436 edges, 1628 lncRNAs nodes, and 1580 mRNAs nodes. The LUAD DE-ceRNA network contains 9953 edges, 1475 lncRNAs nodes, and 1011 mRNAs nodes. By analyzing the topological properties of the network with Cytoscape software, we found that only 1 frequency of shared neighbor node was the highest, with increasing frequency of shared neighbors (Figure 3A, 3E). The distribution of the node degrees of the 2 networks follow the power law distribution (Figure 3B, 3F). Most of the nodes have a small number of connections, and very few nodes have a very large number of connections, indicating that the network construction is a non-random process. The closeness centrality and betweenness centrality of the network show a positive correlation with the number of shared neighbors (Figure 3C, 3D, 3G, 3H). Through the topology analysis of the network, we determined that the network has strong connectivity, and since the size of the network is large, the next step is to optimize the network.
OPTIMIZING THE DE-CERNA NETWORK AND FINDING MARKER GENES:
To optimize the network, we defined known pathogenic genes obtained from 4 gene-associated disease databases as seed nodes to score all nodes in the network by random walk [19]. Then, we reinvested the top genes with random walk scores from high to low ranking into the DE-ceRNA network to get edge information and build a sub-network. Each sub-network is automatically divided into up (orange color) and down (blue color) parts (Figure 4A, 4B). There are nodes in the network that are significantly higher than other nodes, and the contribution of these nodes to the network may be greater. Therefore, we extracted 1 step of neighbor nodes of 2 networks with degree of nodes greater than 20 and maintained the original connection. These nodes that participate in multiple connections in the network are defined as markers such as oncogene PVT1 [20,21], LINC00472 [22], multiple tumor suppressor CDKN2A [23], prognostic Biomarkers in NSCLC FAM83B [24], cancer susceptibility gene BRCA2 in LUSC and HOXA11-AS [25], HNF1A-AS1 [26], LINC00511 [27,28], and HOTAIR [29–32] in LUAD (Figure 4A, 4B). There are many genes in this network that are significantly associated with survival of patients with 2 types of cancer (Figure 5). The expression of the genes PVT1, KC877982.1, LINC01133, MAGI2-AS3, DSCR9, SFTA1P, LINC00968, and SOX2-OTcan significantly distinguish the survival rate of LUSC patients (Figure 6A), and the expression of the genes AC064807.2, LINC00987, LINC01238, TMPO-AS1, OGFRP1, DRAIC, AC023421.1, and RNF144A-AS1 can significantly distinguish the survival rate of LUAD patients (Figure 6B).
COMMON CERNA NETWORK ANALYSIS:
To find common markers, we obtained the intersection of lncRNA-mRNA pairs between 2 cancers to build a common ceRNA network. The network contains 2680 lncRNA-mRNA pairs, 633 lncRNAs and 596 mRNAs (Figure 7A, 7B). Genes that occur at high frequencies in 2 types of cancer are thought to be more important in the treatment and detection of lung cancer. In the common network, mRNAs are also enriched in important pathways and functional terms of cancer development, such as “hsa04110: Cell cycle”, “hsa04115: p53 signaling pathway”, “GO: 0008283~cell proliferation”, “GO: 0051301~cell division”, “h_rbPathway: RB Tumor Suppressor/Checkpoint Signaling in response to DNA damage”, and “GAD_DISEASE: lung cancer” (Figure 8A). We selected cell cycle and p53 pathway modules to display in this common ceRNA network (Figure 8B). For the cell cycle pathway module, we speculated that lncRNAs connected by the coding genes (CDC6, CDK1, E2F2, PKMYT1, TTK, CHEK1, ESPL1, PTTG1, CDC25C, CDC25A, CCNB1, CCNE2, MAD2L1, CCNB2, PLK1, BUB1, BUB1B, ORC6, ORC1, CCNA2) in the pathway may be an important molecule that affects the cell cycle and the increased number of cancer cells, and they exist in an endogenous competitive relationship (CDC25C-LINC00355; ORC6-DDX11-AS1; BUB1/CCNA2/PLK1-VPS9D1-AS1; ORC1/ORC6-LINC00337). In the same way, lncRNA connected by the coding genes (CCNB1, CCNE2, CDK1, CCNB2, SERPINB5, RRM2, CHEK1, GTSE1) may also be an important molecule in the P53 cancer pathway. These close regulatory relationships are more likely to be pathogenic factors of NSCLC.
Discussion
The CeRNA mechanism plays an important role in the development of lung cancer. For example, lncRNA PVT1 regulates expression of HIF1α via functioning as ceRNA for miR-199a-5p in non-small cell lung cancer under hypoxia. We identified some cancer markers by building and analyzing specific or common ceRNA networks [9,33,34]. The specific treatment of 2 subtypes of non-small cell lung cancer is a problem that both researchers and clinicians want to solve, but so far there is no good explanation or evidence [35–37]. We constructed 2 specific networks with mRNA and lncRNA using the specificity of lncRNA to find markers and survival-related genes of 2 different cancers. The representative markers we found in the LUSC network and LUAD network are PVT1, HOTAIR, FAM83B, LINC00472, HOXA11-AS, and CDKN2A, which have also been reported as potential risk factors for malignant tumors such as gastric cancer, glioma, and cervical cancer [20–22,24,25,28–31]. In addition, we analyzed the influence of gene expression on patient survival by combining the clinical information of lung cancer patients. We screened a few genes related to lung cancer patient survival and prognosis, such as PVT1, LINC01133, TMPO-AS1, and DRAIC.
To predict some of the markers in both cancers, we also built a common ceRNA network, and we selected 2 important cancer pathways – Cell cycle and P53 – to analyze the subnetworks associated with them. We found that the CDC25C/CDK1/RRM2-LINC00355 ceRNA relationship can affect non-small cell lung cancer. It is reported in recent studies that exosomes-mediated LINC00355 regulates bladder cancer cell proliferation and invasion [38]. Identification of these markers will further improve understanding of lncRNA-mediated pathway regulation and ceRNA functions. Our study validated the conclusions of some previous experiments and also produced novel and innovative results awaiting verification by biologists. Continued investigation of the ceRNA networks analyzed here will aid in the development of targeted therapy for lung cancer. We will plan to collect more cancer expression data and clinical information to identify molecular markers more precisely using our network-based approach in future studies.
Conclusions
Accumulating evidence suggests that lncRNAs play important roles in lung cancer. However, only a few specific biomarkers have been identified in subtypes of non-small cell lung cancers so far. In this study, the network analysis and functional genes prediction of non-small cell lung cancer were assessed from 2 aspects – specificity and universality. These results suggest that our method identifies different cancer subtypes with remarkable molecular and cancer-related pathways, which would help to improve personalized cancer management. Our results provide a basis and evidence for future research and clinical treatment, and our experimental data will be of value to researchers.
Figures
Figure 1. Classification of samples by lung cancer expression profile. (A) LUSC mRNA expression profile. (B) LUSC lncRNA expression profile. (C) LUAD mRNA expression profile. (D) LUAD lncRNA expression profile. Figure 2. Enrichment analysis of differentially expressed genes. (A) Differentially expressed mRNAs for LUSC. (B) Differentially expressed lncRNAs for LUSC. (C) Differentially expressed mRNAs for LUAD. (D) Differentially expressed lncRNAs for LUAD. (E) Differentially expressed mRNAs enrichment analysis for LUSC. (F) Differentially expressed mRNAs enrichment analysis for LUAD. Figure 3. Topology attribute of specific DE-ceRNA network. (A–D) Attributes of DE-ceRNA network for LUSC. (E–H) Attributes of DE-ceRNA network for LUAD. Figure 4. Visualization of optimized DE-ceRNA Network. (A) Optimized DE-ceRNA Network for LUSC. (B) Optimized DE-ceRNA Network for LUAD. The border colour of the seed nodes is red. Figure 5. Visualization of sub-DE-ceRNA Network. (A) Sub-DE-ceRNA Network for LUSC. (B) Sub-DE-ceRNA Network for LUAD. The border colour of the seed nodes is red. Figure 6. KM survival curve of marker genes. (A) KM survival curves of lncRNAs and mRNAs for LUSC. (B) KM survival curves of lncRNAs and mRNAs for LUAD. Figure 7. Common ceRNA network for non-small cell lung cancer. (A) The intersection of the ceRNA pairs, lncRNAs, or mRNAs for 2 kinds of non-small cell lung cancer. (B) Visualization of common ceRNA network. Figure 8. Common markers for non-small cell lung cancer. (A) Enrichment analysis of mRNAs in common ceRNA network. (B) Network modules related to cancer-related pathway.References
1. Bray F, Ferlay J, Soerjomataram I, Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries: Cancer J Clin, 2018; 68; 394-424
2. Wang XW, Guo QQ, Wei Y, Construction of a competing endogenous RNA network using differentially expressed lncRNAs, miRNAs and mRNAs in nonsmall cell lung cancer: Oncol Rep, 2019; 42; 2402-15
3. Qiao F, Li N, Li W, Integrative bioinformatics analysis reveals potential long non-coding RNA biomarkers and analysis of function in non-smoking females with lung cancer: Med Sci Monit, 2018; 24; 5771-78
4. Lei Y, Shi Y, Duan J, Identification of alternative splicing and lncRNA genes in pathogenesis of small cell lung cancer based on their RNA sequencing: Adv Clin Exp Med, 2019; 28; 1043-50
5. Wang X, Yin H, Zhang L, The construction and analysis of the aberrant lncRNA-miRNA-mRNA network in non-small cell lung cancer: J Thorac Dis, 2019; 11; 1772-78
6. Amaral PP, Dinger ME, Mercer TR, Mattick JS, The eukaryotic genome as an RNA machine: Science, 2008; 319; 1787-89
7. Carninci P, Kasukawa T, Katayama S, The transcriptional landscape of the mammalian genome: Science, 2005; 309; 1559-63
8. Kapranov P, Cheng J, Dike S, RNA maps reveal new RNA classes and a possible function for pervasive transcription: Science, 2007; 316; 1484-88
9. Nakaya HI, Amaral PP, Louro R, Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription: Genome Biol, 2007; 8; R43
10. Ebert MS, Sharp PA, Emerging roles for natural microRNA sponges: Curr Biol, 2010; 20; R858-61
11. Arvey A, Larsson E, Sander C, Target mRNA abundance dilutes microRNA and siRNA activity: Mol Syst Biol, 2010; 6; 363
12. Calin GA, Liu CG, Ferracin M, Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas: Cancer Cell, 2007; 12; 215-29
13. Xiong DD, Li ZY, Liang L, The LncRNA NEAT1 accelerates lung adenocarcinoma deterioration and binds to mir-193a-3p as a competitive endogenous RNA: Cell Physiol Biochem, 2018; 48; 905-18
14. Zhang L, Fang F, He X, Long noncoding RNA TP73-AS1 promotes non-small cell lung cancer progression by competitively sponging miR-449a/EZH2: Biomed Pharmacother, 2018; 104; 705-11
15. Ning P, Wu Z, Hu A, Integrated genomic analyses of lung squamous cell carcinoma for identification of a possible competitive endogenous RNA network by means of TCGA datasets: Peer J, 2018; 6; e4254
16. Salmena L, Poliseno L, Tay Y, A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language?: Cell, 2011; 146; 353-58
17. Panzitt K, Tschernatsch MM, Guelly C, Characterization of HULC, a novel gene with striking up-regulation in hepatocellular carcinoma, as noncoding RNA: Gastroenterology, 2007; 132; 330-42
18. Wang J, Liu X, Wu H, CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer: Nucleic Acids Res, 2010; 38; 5366-83
19. Abdullahi A, Shohaimi S, Kilicman A, Ibrahim MH, Stochastic models in seed dispersals: random walks and birth-death processes: J Biol Dyn, 2019; 13; 345-61
20. Wang Y, Chen W, Lian J, The lncRNA PVT1 regulates nasopharyngeal carcinoma cell proliferation via activating the KAT2A acetyltransferase and stabilizing HIF-1alpha: Cell Death Differ, 2020; 27(2); 695-710
21. Wang C, Han C, Zhang Y, Liu F, LncRNA PVT1 regulate expression of HIF1alpha via functioning as ceRNA for miR199a5p in nonsmall cell lung cancer under hypoxia: Mol Med Rep, 2018; 17; 1105-10
22. Chen L, Zhang W, Li DY, Regulatory network analysis of LINC00472, a long noncoding RNA downregulated by DNA hypermethylation in colorectal cancer: Clin Genet, 2018; 93; 1189-98
23. Zeng H, Jorapur A, Shain AH, Bi-allelic Loss of CDKN2A initiates melanoma invasion via BRN2 activation: Cancer Cell, 2018; 34; 56-68.e9
24. Richtmann S, Wilkens D, Warth A, FAM83A and FAM83B as prognostic biomarkers and potential new therapeutic targets in NSCLC: Cancers (Basel), 2019; 11(5) pii: E652
25. Sun M, Nie F, Wang Y, LncRNA HOXA11-AS promotes proliferation and invasion of gastric cancer by scaffolding the chromatin modification factors PRC2, LSD1, and DNMT1: Cancer Res, 2016; 76; 6299-310
26. Liu L, Chen Y, Li Q, Duan P, lncRNA HNF1A-AS1 modulates non-small cell lung cancer progression by targeting miR-149-5p/Cdk6: J Cell Biochem, 2019; 120; 18736-50
27. Ding J, Yang C, Yang S, LINC00511 interacts with miR-765 and modulates tongue squamous cell carcinoma progression by targeting LAMC2: J Oral Pathol Med, 2018; 47; 468-76
28. Yan L, Wu X, Liu Y, Xian W, LncRNA Linc00511 promotes osteosarcoma cell proliferation and migration through sponging miR-765: J Cell Biochem, 2018 [Epub ahead of print]
29. Li J, Li S, Chen Z, miR-326 reverses chemoresistance in human lung adenocarcinoma cells by targeting specificity protein 1: Tumour Biol, 2016; 37; 13287-94
30. Loewen G, Jayawickramarajah J, Zhuo Y, Shan B, Functions of lncRNA HOTAIR in lung cancer: J Hematol Oncol, 2014; 7; 90
31. Malissovas N, Ninou E, Michail A, Politis PK, Targeting long non-coding RNAs in nervous system cancers: New insights in prognosis, diagnosis and therapy: Curr Med Chem, 2019; 26; 5649-63
32. Zhuang Y, Wang X, Nguyen HT, Induction of long intergenic non-coding RNA HOTAIR in lung cancer cells by type I collagen: J Hematol Oncol, 2013; 6; 35
33. Shao T, Wu A, Chen J, Identification of module biomarkers from the dysregulated ceRNA-ceRNA interaction network in lung adenocarcinoma: Mol Biosyst, 2015; 11; 3048-58
34. Sui J, Xu SY, Han J, Integrated analysis of competing endogenous RNA network revealing lncRNAs as potential prognostic biomarkers in human lung squamous cell carcinoma: Oncotarget, 2017; 8; 65997-6018
35. Zhu TG, Xiao X, Wei Q, Revealing potential long non-coding RNA biomarkers in lung adenocarcinoma using long non-coding RNA-mediated competitive endogenous RNA network: Braz J Med Biol Res, 2017; 50; e6297
36. Li J, Tian H, Pan J, Pecanex functions as a competitive endogenous RNA of S-phase kinase associated protein 2 in lung cancer: Cancer Lett, 2017; 406; 36-46
37. Sui J, Li YH, Zhang YQ, Integrated analysis of long non-coding RNA-associated ceRNA network reveals potential lncRNA biomarkers in human lung adenocarcinoma: Int J Oncol, 2016; 49; 2023-36
38. Yan L, Wang P, Fang W, Liang C, Cancer-associated fibroblasts-derived exosomes-mediated transfer of LINC00355 regulates bladder cancer cell proliferation and invasion: Cell Biochem Funct, 2020; 38(3); 257-65
Figures
In Press
Clinical Research
Evaluating Lipiodol Efficacy in Proximal Tubal Occlusion Treatment: Single-Center Experience and Literature...Med Sci Monit In Press; DOI: 10.12659/MSM.946266
Review article
Impact of Traditional Chinese Medicine Antioxidants on Oxidative Stress and Drug-Induced Liver Injury: A Re...Med Sci Monit In Press; DOI: 10.12659/MSM.945147
Meta-Analysis
Reliability of Extraoral Scanners in Capturing 3D Geometry for Dental Prostheses: A Systematic ReviewMed Sci Monit In Press; DOI: 10.12659/MSM.946470
Review article
Skin Barrier Dysfunction in Acne Vulgaris: Pathogenesis and Therapeutic ApproachesMed Sci Monit In Press; DOI: 10.12659/MSM.945336
Most Viewed Current Articles
17 Jan 2024 : Review article 6,957,034
Vaccination Guidelines for Pregnant Women: Addressing COVID-19 and the Omicron VariantDOI :10.12659/MSM.942799
Med Sci Monit 2024; 30:e942799
14 Dec 2022 : Clinical Research 1,969,335
Prevalence and Variability of Allergen-Specific Immunoglobulin E in Patients with Elevated Tryptase LevelsDOI :10.12659/MSM.937990
Med Sci Monit 2022; 28:e937990
16 May 2023 : Clinical Research 697,010
Electrophysiological Testing for an Auditory Processing Disorder and Reading Performance in 54 School Stude...DOI :10.12659/MSM.940387
Med Sci Monit 2023; 29:e940387
07 Jan 2022 : Meta-Analysis 263,050
Efficacy and Safety of Light Therapy as a Home Treatment for Motor and Non-Motor Symptoms of Parkinson Dise...DOI :10.12659/MSM.935074
Med Sci Monit 2022; 28:e935074