Abstract:Background and objective: GCMN is a sporadic disease with an incidence ranging from 1/20,000 to 1/500000. So far, several studies have found that GCMN is related to somatic mutations, but most of them have focused on known pathogenic genes, and transcriptome sequencing based on large datasets is relatively uncommon. At present, the use of next-generation sequencing technologies and bioinformatics platforms makes genomic information study more comprehensive and efficient. In this study, the transcriptome differences between GCMN lesions and surrounding normal skin tissues were investigated using high-throughput transcriptome sequencing, and hub genes and pathways related to pathogenesis were identified, providing a theoretical foundation for further research into the pathogenesis of GCMN. Methods: Pathological skin tissue and surrounding normal skin tissue from GCMN patients, namely the pathological group (PG) and the control group (CG), were obtained. 1. All specimens were stained with HE to ensure that the samples met the experimental requirements. 2. Ten pairs of specimens were selected for high-throughput transcriptome sequencing, and the differentially expressed genes (DEGs) between the PG and the CG were obtained. The DEGs were analyzed by clusterProfiler R software for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The function of the subnetwork was analyzed and the hub genes were identified by the STRING database and Cytoscape software. 3. The expression differences of hub genes PTGS2, EGF, and SOX10 in pathological skin tissues and normal skin tissues were verified by qRT-PCR and immunofluorescence staining. Results: 1. HE staining revealed a lot of melanocytes in the dermis and subcutaneous tissues. They were found around the hair follicles, sweat glands, sebaceous glands, and blood vessel walls, or in a specific pattern. 2. The screening threshold was set at p < 0.01 and |log2fc|<1, and a total of 1163 DEGs were discovered between the PG and CG, with 519 genes up-regulated and 644 genes down-regulated in the pathological tissues. According to the GO functional analysis, 29 biological processes, 18 cell compositions, and 17 molecular functions were significantly enriched, with the majority of them being related to keratinocytes and the extracellular matrix. There were 779 nodes and 2359 interactions in the protein interaction network. Using the MCODE plug-in, the network was divided into 25 functional clusters. According to ClueGO results, Cluster5 was involved in melanin biosynthesis and melanocyte proliferation. Using 11 operation methods in the Cytohubba plug-in, PTGS2, EGF, and SOX10 in Cluster5 were chosen as hub genes. 3. qRT-PCR and immunofluorescent staining revealed that compared to normal skin tissue, the expression of SOX10 was significantly up-regulated, and the expression of PTGS2 and EGF was significantly down-regulated in pathological skin tissue(P < 0.001). Conclusions: In GCMN, keratinocytes and extracellular matrix may directly and indirectly affect melanocyte activity. PTGS2, EGF, and SOX10 are important genes and significantly differentially expressed in pathological and normal skin tissues. These findings may serve as a springboard for future research.

Genomic Structure of Mouse Tbx2 and Detection of Expression of Tbx2 in Normal and Malignance Melanophore by RT-PCR

Expression of mouse Tbx2 gene in normal and malignant melanophores by RT-PCR

Genome-wide Screening of Aberrant DNA Methylation Which Associated with Gene Expression in Mouse Skin Cancers.

Molecular cloning and characterization of a novel human BTB domain-containing gene, , which is down-regulated in glioma

Anti-tumor immune response varies among individuals: A gene expression profiling of mouse melanoma

Cloning and Characterization of Mouse Mterf Encoding a Mitochondrial Transcriptional Termination Factor.

The Expression, Purification, and Characterization of a Ras Oncogene (bras2) in Silkworm (bombyx Mori).

Effects of RNA Methylation N6-methyladenosine Regulators on Malignant Progression and Prognosis of Melanoma.

Tumor Associated Gene Expression of Human Leukemia Cell Line with High Tumorigenicity in Nude Mice

Transcriptomic analysis of mRNA expression in giant congenital melanocytic nevi

Bioinformatics Analysis to Identify the Critical Genes, Micrornas and Long Noncoding RNAs in Melanoma

Array Comparative Genome Hybridization for Tumor Classification and Gene Discovery in Mouse Models of Malignant Melanoma.

Transgenic mouse model for skin malignant melanoma

Microphthalmia-Associated Transcription Factor/T-Box Factor-2 Axis Acts Through Cyclin D1 to Regulate Melanocyte Proliferation

Induction Of Melanoma In Murine Macrophage Inflammatory Protein 2 Transgenic Mice Heterozygous For Inhibitor Of Kinase/Alternate Reading Frame

A Comparative Transcriptomic Analysis Of Uveal Melanoma And Normal Uveal Melanocyte

Comprehensive analysis and identification of key genes and signaling pathways in the occurrence and metastasis of cutaneous melanoma

The transcription factor TBX2 regulates melanogenesis in melanocytes by repressing Oca2

[Bioinformatics analysis to the differentially expressed genes of normal mucosa and carcinoma of colon].

Molecular characteristics and therapeutic implications of Toll-like receptor signaling pathway in melanoma

[A New Melanoma Antigen-Encoding Gene Subfamily in Human Chromosome X].