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.

Loss-of-function variants in GLMN are associated with generalized skin hyperpigmentation with or without glomuvenous malformation

Amyloidosis Cutis Dyschromica Cases Caused by GPNMB Mutations with Different Inheritance Patterns.

Transcriptomic analysis of mRNA expression in giant congenital melanocytic nevi

Identification of a Novel MLPH Missense Mutation in a Chinese Griscelli Syndrome 3 Patient

Genetic and phenotypic heterogeneity of multiple lentigines and precise diagnosis in four Chinese families with multiple lentigines

Validation of Susceptibility Loci for Vitiligo Identified by GWAS in the Chinese Han Population

Identification of 11 Potentially Relevant Gene Mutations Involved in Growth Retardation, Intellectual Disability, Joint Contracture, and Hepatopathy

LB1599 Novel microenvironmental defects revealed by large scale transcriptome sequencing of lesional and nonlesional skin of vitiligo patients

Novel Genetic Variants of ALG6 and GALNTL4 of the Glycosylation Pathway Predict Cutaneous Melanoma-Specific Survival

A Novel Nonsense GLI3 Variant is Associated with Polydactyly and Syndactyly in a Family by Blocking the Sonic Hedgehog Signaling Pathway

Association of novel MUC16, MAP3K15 and ABCA1 mutation with giant congenital melanocytic nevus

Genetic Polymorphisms of GZMB and Vitiligo: A Genetic Association Study Based on Chinese Han Population.

A postzygotic GNA13 variant upregulates the RHOA/ROCK pathway and alters melanocyte function in a mosaic skin hypopigmentation syndrome

Two nonsense GLI3 variants are associated with polydactyly and syndactyly in two families by affecting the sonic hedgehog signaling pathway

A gain-of-function variant in SREBF1 causes generalized skin hyperpigmentation with congenital cataracts

GGCX mutations in a patient with overlapping pseudoxanthoma elasticum/cutis laxa-like phenotype

Deficient Mismatch Repair is Detected in Large-to-giant Congenital Melanocytic Naevi: Providing New Insight into Aetiology and Diagnosis

Uncovering molecular mechanisms for amelanotic/hypopigmented primary cutaneous melanoma

Gain-of-Function Mutation of KIT Ligand on Melanin Synthesis Causes Familial Progressive Hyperpigmentation

Association of GZMB polymorphisms and susceptibility to non-segmental vitiligo in a Korean population

Two Novel Frameshift Mutations in the GLI3 Gene Underlie Non-Syndromic Polydactyly in Chinese Families