Au total, 3290 DEG ont été identifiés à partir des 23348 gènes de GSE84884 (pSS), dont 2659 gènes régulés à la hausse et 631 gènes régulés à la baisse (Figure 1A). Pour GSE51092 (pSS), un total de 3290 DEG ont été identifiés à partir des 11409 gènes, dont 667 gènes régulés à la hausse et 587 gènes régulés à la baisse (Figure 1B). La base de données GeneCards a obtenu 102 DEG liés au pSS ovarien, et le score de corrélation des critères de dépistage était de ≥20. L’union et la déduplication des DEG liés au pSS obtenus à partir de GEO et de GeneCards ont abouti à 453 DEG (pSS-DEG). Au total, 7154 DEG ont été identifiés à partir des 25440 gènes de GSE32863 (LUAD), dont 3934 gènes régulés à la hausse et 3220 gènes régulés à la baisse (Figure 1C). En GSE75037 (LUAD), un total de 11235 DEG ont été identifiés à partir des 25440 gènes, dont 6125 gènes régulés à la hausse et 5110 gènes régulés à la baisse (Figure 1D). La base de données GeneCards a obtenu 600 DEG ovariens liés au LUAD, et le score de corrélation des critères de dépistage était de ≥20. L’union et la déduplication des DEG liés au LUAD obtenues à partir de GEO et de GeneCards ont abouti à 7478 DEG (LUAD-DEG). Au total, 233 DEG partagés (pSS-LUAD-DEG) ont été sélectionnés entre les pSS-DEG et les LUAD-DEG et visualisés à l’aide de diagrammes de Venn (tableau supplémentaire 1 et figure 1E).
En ce qui concerne l’analyse de l’enrichissement de la voie KEGG, les 10 principales voies de signalisation significatives des pSS-LUAD-DEG identifiées étaient liées aux voies métaboliques après suppression des voies de signalisation liées à la maladie. Plus précisément, ces voies comprennent la voie de signalisation PI3K-Akt, la voie de signalisation MAPK, l’interaction cytokine-récepteur des cytokines, la nécroptose, la régulation du cytosquelette d’actine, la voie de signalisation FoxO, la voie de signalisation du récepteur de type NOD, la sénescence cellulaire, la différenciation des cellules Th17 et la voie de signalisation JAK-STAT (Figure 2A). Les 10 principales voies de signalisation KEGG et les pSS-LUAD-DEG ont été visualisés, comme le montre la figure 2B.
L’analyse d’enrichissement GO des pSS-LUAD-DEG a révélé les 10 principaux processus biologiques suivants (Figure 3), et a été enrichi en composants cellulaires (CC) : réponse au virus, réponse immunitaire innée, réponse de défense au virus, réponse de défense au symbiote, régulation du processus viral, régulation du cycle de vie viral, régulation négative du processus viral, régulation négative de la réplication du génome viral, régulation de la réplication du génome viral, réponse immunitaire innée antivirale ; processus biologique (BP) : côté de la membrane, face externe de la membrane plasmique, matrice extracellulaire, structure d’encapsulation externe, matrice extracellulaire contenant du collagène, trimère de collagène, radeau de membrane, microdomaine de membrane, radeau de membrane plasmique, cavéole ; fonction moléculaire (MF) : liaison du récepteur des cytokines, activité des cytokines, activité régulatrice du récepteur de signalisation, liaison de la superfamille du récepteur du facteur de nécrose tumorale, activité activatrice du récepteur de signalisation, liaison du récepteur du facteur de nécrose tumorale, activité du ligand du récepteur, activité d’homodimérisation des protéines, activité de la protéine kinase du récepteur transmembranaire, activité de la protéine kinase.
Le réseau PPI comprenait 99 nœuds et 466 arêtes (Figure 4A). Les 20 gènes les plus élevés du réseau PPI sont STAT3, STAT1, TP53, TNF, IL6, IFNG, EGFR, ISG15, CCNA2, IL10, JAK2, MX1, IL1B, IFIT1, AKT1, SMAD3, ASPM, CCNB2, NUSAP1 et CEP55, qui servent de gènes pivots (figure 4B).
Une analyse de la courbe ROC a été effectuée pour évaluer la précision diagnostique des 20 gènes centraux (Figure 5). Les valeurs de l’aire sous la courbe (AUC) pour les courbes ROC de STAT3, STAT1, TP53, TNF, IL6, IFNG, EGFR, ISG15, CCNA2, IL10, JAK2, MX1, IL1B, IFIT1, AKT1, SMAD3, ASPM, CCNB2, NUSAP1 et CEP55 étaient de 0,543, 0,840, 0,724, 0,892, 0,965, 0,529, 0,721, 0,763, 0,933, 0,784, 0,836, 0,648, 0,936, 0,689, 0,662, 0,548, 0,960, 0,939, 0,934 et 0,953 pour GSE32863 (Figure 5A, D), 0,620, 0,811, 0,676, 0,878, 0,955, 0,546, 0,642, 0,765, 0,945, 0,759, 0,804, 0,625, 0,938, 0,753, 0,641, 0,540, 0,958, 0,945, 0,908 et 0,941 pour GSE75037 (figures 5B et E). De plus, pour l’ensemble de validation GSE31210 (Figure 5C, F), les valeurs AUC étaient de 0,505, 0,676, 0,819, 0,688, 0,863, 0,540, 0,720, 0,665, 0,814, 0,562, 0,766, 0,661, 0,806, 0,528, 0,607, 0,652, 0,906, 0,927, 0,910 et 0,898. Les valeurs de l’ASC pour IL6, CCNA2, JAK2, IL1B, ASPM, CCNB2, NUSAP1 et CEP55 étaient de >0,7, ce qui suggère que les huit ont une signification diagnostique.
Dans les tissus pulmonaires des souris BALB/c et NOD/Ltj, seule une petite quantité de facteurs pro-inflammatoires IL-6 et IL-1β a été détectée. Après stimulation avec des PM2,5, l’expression de l’IL-6 et de l’IL-1β a augmenté de manière significative, avec une élévation nettement plus élevée observée chez les souris NOD/Ltj par rapport aux souris BALB/c (Figure 6).
Après stimulation par les PM2,5, les niveaux d’expression de JAK2/STAT3 n’ont pas changé de manière significative dans les tissus pulmonaires des souris BALB/c et NOD/Ltj. Cependant, les niveaux de p-JAK2 et de p-STAT3 ont augmenté de manière significative. Les souris NOD/Ltj ont montré une expression significativement plus élevée de p-JAK2 et de p-STAT3 par rapport aux souris BALB/c après stimulation PM2,5 (Figure 7).
Après la stimulation par les PM2,5, l’expression de l’ARNm d’ASPM et de NUSAP1 n’a pas changé de manière significative chez les souris BALB/c et les souris NOD/Ltj, mais l’expression de l’ARNm de CCNA2, CCNB2 et CEP55 a augmenté (Figure 8). Notamment, l’expression de l’ARNm de CCNA2, CCNB2 et CEP55 était significativement plus élevée chez les souris NOD/Ltj par rapport aux souris BALB/c après stimulation PM2.5.

Figure 1 : La carte du volcan et le diagramme de Venn des DEG. (A) La carte du volcan de GSE84844. (B) La carte des volcans de GSE51092. (C) La carte des volcans de GSE32863. (D) La carte du volcan de GSE75037. (E) Diagramme de Venn des pSS-DEG et des LUAD-DEG. Les gènes régulés à la hausse sont colorés en rouge ; Les gènes régulés à la baisse sont colorés en vert. Abréviations : pSS = syndrome de Sjögren primaire ; LUAD = adénocarcinome pulmonaire ; DEG = gènes exprimés de manière différentielle. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 2 : Enrichissement de la voie KEGG des pSS-LUAD-DEG. (A) TOP10 Carte d’enrichissement KEGG des pSS-LUAD-DEG. (B) L’association entre les voies de signalisation TOP10 KEGG avec les pSS-LUAD-DEG. Abréviations : KEGG = Encyclopédie de Kyoto des gènes et des génomes ; pSS = syndrome de Sjögren primaire ; LUAD = adénocarcinome pulmonaire ; DEG = gènes exprimés de manière différentielle. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Analyse de l’enrichissement des GO des pSS-LUAD-DEG. Abréviations : GO = Ontologie génétique ; pSS = syndrome de Sjögren primaire ; LUAD = adénocarcinome pulmonaire ; DEG = gènes exprimés de manière différentielle. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Réseau PPI des pSS-LUAD-DEG. (A) Visualisation du réseau PPI. Les nuances plus rouges indiquent une connectivité plus forte du nœud. (B) Les 20 premiers gènes avec des degrés plus élevés dans le réseau PPI. Abréviations : IPP = interaction protéine-protéine ; pSS = syndrome de Sjögren primaire ; LUAD = adénocarcinome pulmonaire ; DEG = gènes exprimés de manière différentielle. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Analyse de la courbe ROC des gènes du moyeu. (A-F) Analyse de la courbe ROC des gènes du moyeu dans GSE32863, GSE75037 et GSE31210. Abréviations : ROC = Receiver Operating Characteristic. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : L’impact des PM2,5 sur les facteurs pro-inflammatoires dans les tissus pulmonaires de souris. (A) Expression de l’IL-6 dans les tissus pulmonaires de souris de chaque groupe. (B) Expression de l’IL-1β dans les tissus pulmonaires de souris de chaque groupe. a. groupe témoin à blanc ; b. groupe pSS ; c. Groupe PM2,5 ; d. Groupe pSS-PM2.5. * p < 0,05 par rapport au groupe témoin à blanc ; # p < 0,05 par rapport au groupe pSS ; Δ p < 0,05 par rapport au groupe PM2,5. Abréviations : PM2,5 = matières particulaires 2,5. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 7 : L’impact des PM2,5 sur la voie de signalisation JAK2/STAT3 dans les tissus pulmonaires de souris. (A) Expression de p-JAK2, JAK2, p-STAT3, STAT3 et β-actine dans les tissus pulmonaires de souris de chaque groupe. (B) Le rapport p-JAK2/JAK2 dans les tissus pulmonaires de souris de chaque groupe. (C) Le rapport p-STAT3/STAT3 dans les tissus pulmonaires de souris de chaque groupe. a. groupe témoin à blanc ; b. groupe pSS ; c. Groupe PM2,5 ; d. Groupe pSS-PM2.5. * p < 0,05 par rapport au groupe témoin à blanc ; # p < 0,05 par rapport au groupe pSS ; Δ p < 0,05 par rapport au groupe PM2,5. Abréviations : PM2,5 = matières particulaires 2,5. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 8 : L’impact des PM2,5 sur les gènes liés à LUAD dans les tissus pulmonaires de souris. (A) Expression de l’ARNm CCNB2 dans les tissus pulmonaires de souris de chaque groupe. (B) Expression de l’ARNm de CCNA2 dans les tissus pulmonaires de souris de chaque groupe. (C) Expression de l’ARNm ASPM dans les tissus pulmonaires de souris de chaque groupe. (D) Expression de l’ARNm de NUSAP1 dans les tissus pulmonaires de souris de chaque groupe. (E) Expression de l’ARNm CEP55 dans les tissus pulmonaires de souris de chaque groupe. a. groupe témoin à blanc ; b. groupe pSS ; c. Groupe PM2,5 ; d. Groupe pSS-PM2.5. * p < 0,05 par rapport au groupe témoin à blanc ; # p < 0,05 par rapport au groupe pSS ; Δ p < 0,05 par rapport au groupe PM2,5. Abréviations : PM2,5 = matières particulaires 2,5. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Tableau supplémentaire 1 : Liste des DEG partagées (pSS-LUAD-DEGs) sélectionnées entre les pSS-DEG et les LUAD-DEG. Veuillez cliquer ici pour télécharger ce fichier.
Fichier de codage supplémentaire 1 : Code R pour le filtrage des DEG. Veuillez cliquer ici pour télécharger ce fichier.