Differenzielle mRNA-Expression von m6A-bezogenen Genen bei Brustkrebs
Die Analyse der TCGA- und GTEx-Daten (1.034 Tumore, 104 NATs, 92 gesunde Gewebe) zeigte, dass die meisten m6A-bezogenen Gene differenziell exprimiert wurden (Abbildung 1). Im Vergleich zu gesunden Brustgeweben waren 10 Gene signifikant hochreguliert: KIAA1429, RBM15, RBM15B, HNRNPC, YTHDF1, YTHDF2, YTHDF3, IGF2BP1, IGF2BP3 und ALKBH5 (alle P < 0,0001, außer RBM15 P < 0,001). Acht Gene wurden herunterreguliert: WTAP, METTL3, METTL14, HNRNPA2B1, YTHDC1, YTHDC2, IGF2BP2 und FTO (P < 0,05 bis P < 0,0001). METTL16 und eIF3A zeigten keinen signifikanten Unterschied. Beim Vergleich von Tumoren mit NATs wurde HNRNPA2B1 hochreguliert, während YTHDF3, eIF3A und ALKBH5 herunterreguliert waren; RBM15, RBM15B und YTHDF2 zeigten keinen Unterschied (Abbildung 1). Die Verwendung von NATs als Kontrollgruppen wurde hinterfragt, was möglicherweise die Diskrepanzen mit gesunden Gewebevergleichen erklärt18,19. Interessanterweise wurde im CPTAC-Datensatz im Vergleich zu normalen Geweben ein Trend zu einer höheren YTHDF3-Proteinhäufigkeit in allen Brustkrebs-Subtypen beobachtet (P < 0,0001). Darüber hinaus wurde eine signifikant höhere mRNA-Expression von YTHDF3 nur in luminalen Subtypen von Brustkrebs im Vergleich zu normalen Geweben beobachtet (P = 0,004, < 0,01) (Ergänzende Abbildung S1). Insgesamt ist eine erhebliche Dysregulation der m6A-bezogenen Gene bei Brustkrebs erkennbar.
Promotormethylierung und Genamplifikation von m6A-bezogenen Genen bei Brustkrebs
Die DNA-Methylierung spielt eine wichtige Rolle bei der epigenetischen Regulation der Gentranskription, und die Promotor-Hypermethylierung hemmt die mRNA-Expression20. Die MethHC-Analyse zeigte signifikante Promoter-Methylierungsunterschiede nur für IGF2BP1, IGF2BP2 und IGF2BP3 bei Tumoren im Vergleich zu Kontrollgruppen (P < 0,005). Die Downregulation der IGF2BP2 mRNA kann teilweise auf ihre Hypermethylierung zurückzuführen sein. cBioPortal-Analyse zeigte, dass WTAP und ALKBH5 eine mRNA-Downregulation aufwiesen, während KIAA1429 (VIRMA), YTHDF1, YTHDF3 und IGF2BP1 eine mit ihrer Upregulation assoziierte Genamplifikation aufwiesen (ergänzende Abbildung S2, ergänzende Abbildung S3 und ergänzende Abbildung S4). Diese Ergebnisse deuten darauf hin, dass die mRNA-Werte von m6A-bezogenen Genen nicht überwiegend durch Promotormethylierung oder Genamplifikation bestimmt werden.
Überlebensanalyse von m6A-bezogenen Genen bei Brustkrebs
Die Überlebensanalyse von 947 Patienten zeigte, dass eine niedrige RBM15B-Expression und eine hohe Expression von METTL16, HNRNPC, YTHDF1, YTHDF3 und IGF2BP1 signifikant mit einer verringerten Gesamtüberlebensrate korreliert waren (Abbildung 2 und Abbildung 3). Die multivariate Cox-Regression, Anpassung für Alter, Stadium und Rasse, bestätigte, dass die Hochregulierung von YTHDF3 ein unabhängiger Prädiktor für ein schlechtes Gesamtüberleben ist (HR: 1,024, 95%-KI: 1,003–1,046, P = 0,024) (Tabelle 1). Zur externen Validierung zeigte die Analyse mit dem Online-Tool KMplotter, dass eine hohe YTHDF3-Expression eine schlechte Krankheitsprognose vorhersagt (ergänzende Abbildung S5). Andere Gene waren unabhängig davon nicht signifikant (Ergänzende Tabelle S1).
Potenzielles molekulares Netzwerk und Immuninfiltration des m6A-bezogenen Gens YTHDF3
Die STRING-Analyse von YTHDF3 (mit begrenzten Interaktoren auf 50) identifizierte ein Netzwerk mit 224 Kanten und signifikanter Anreicherung (P < 1e-16). YTHDF3 ist eng mit METTL3, ALKBH5, WTAP, METTL14 und anderen verbunden (Abbildung 4A). Die Metascape-Anreicherungsanalyse identifizierte 13 signifikante funktionale Cluster (P < 0,01, Mindestanzahl von 3, Anreicherungsfaktor > 1,5; Abbildung 4B). Der am höchsten bewertete Cluster war die mRNA-Stoffwechselregulation (GO1903311, 22 Gene). Weitere krebsrelevante Cluster waren unter anderem DNA-Dealkylationsreparatur (GO0006307), Pri-miRNA-Verarbeitung (GO0031053) und Telomerase-Stabilisierung (GO1904356). Die vollständige Liste der Clusternamen, assoziierten Gene und Anreicherungsstatistiken ist in der ergänzenden Tabelle S2 bereitgestellt. Darüber hinaus zeigten die Ergebnisse signifikante positive Korrelationen zwischen der YTHDF3-Expression und Immunzellen (alle P < 0,01) (Ergänzende Abbildung S6). Daher könnte die Upregulation von YTHDF3 mit mRNA-Stoffwechsel, DNA-Reparatur, Telomererhaltung und Immuninfiltration in Verbindung stehen.
Hochregulierte mRNA-Expression von YTHDF3 bei Brustkrebs
Brustkrebsgewebe und übereinstimmende angrenzende normale Brustgewebe wurden von 20 Patienten gesammelt, die zwischen Juli 2023 und Oktober 2023 in der Abteilung für Brustchirurgie, dem angeschlossenen Krankenhaus der Southwest Medical University, operiert wurden. qRT-PCR-Analyse dieser gepaarten Proben zeigte eine signifikant höhere YTHDF3-mRNA-Expression im Tumorgewebe als in benachbarten normalen Geweben (P = 0,001) (Ergänzende Abbildung S7 und Ergänzungstabelle S3).
Datenverfügbarkeit
Alle in dieser Studie analysierten Datensätze sind in den folgenden Open-Access-Repositorien öffentlich verfügbar. Genexpression und klinische Daten: Brustkrebs-RNA-seq-Expressionsdaten und entsprechende klinische Informationen wurden aus dem Cancer Genome Atlas (TCGA-BRCA) über den Human Protein Atlas (https://www.proteinatlas.org/) gewonnen. Weitere normale Brustgewebeexpressionsdaten wurden aus dem Genotype-Tissue Expression (GTEx)-Projekt (dbGaP-Accession phs000424.v8.p2) mit Breast Cancer Gene-Expression Miner v4.7 (http://bcgenex.ico.unicancer.fr/) gewonnen. Die Überlebensvalidierung erfolgte mit dem Kaplan-Meier-Plotter (https://kmplot.com). Protein- und mRNA-Expressionsdaten wurden über UALCAN (http://ualcan.path.uab.edu/) gewonnen. Genetische Veränderungen: Somatische Mutations- und Kopienzahl-Veränderungsdaten für den Datensatz für invasiven Brustkrebs (TCGA, Provisional, 1108 Proben) wurden über cBioPortal (https://www.cbioportal.org/) abgerufen. DNA-Promotormethylierung: Methylierungsdaten wurden aus MethHC (http://methhc.mbc.nctu.edu.tw/) gewonnen, das 839 Brustkrebsproben und teilweise gepaarte benachbarte normale Kontrollgruppen enthält. Protein-Protein-Interaktion und Weganreicherung: Die funktionelle Netzwerkanalyse wurde mit STRING v12.0 (https://string-db.org/, Konfidenzwert 0,4) und Metascape (http://metascape.org/) durchgeführt. Die Immuninfiltrationsanalyse wurde mit TIMER (https://cistrome.shinyapps.io/timer/) durchgeführt. Alle zusätzlichen Daten, die den Zahlen zugrunde liegen, sind auf angemessene Anfrage vom korrespondierenden Autor verfügbar.

Abbildung 1. Differenzielle mRNA-Expression von N6-Methyladenosin (m6A)-bezogenen Genen bei Brustkrebs. Rot steht für eine Hochregulierung; Blau zeigt eine Downregulation an, Schwarz keinen Unterschied in der Expression und Weiß keine Genexpressionsdaten. Die Farbcodierung spiegelt die statistische Signifikanz wider. Der Dunnett-Tukey-Kramer-Test wurde für paarweise Vergleiche verwendet. Hinweis: "Normal" umfasst den Vergleich zwischen Tumoren und NATs bei TCGA sowie den Vergleich zwischen Tumoren und gesundem Gewebe in GTEx. Abkürzungen: NATs = normale benachbarte Gewebe; TCGA = Der Cancer Genome Atlas; GTEx = Genotyp-Gewebeexpression. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 2. Zusammenhang zwischen abweichender mRNA-Expression der m6A-"Schriftsteller" RBM15B und METTL16 und dem Gesamtüberleben bei Brustkrebspatientinnen. (A) Eine niedrige Expression von RBM15B war mit einem kürzeren Gesamtüberleben verbunden (Log-Rank-Test, P = 0,008, <0,01); (B) Eine hohe Expression von METTL16 war mit einer schlechten Prognose assoziiert (Log-Rang-Test, P = 0,013, <0,05). Die Patienten wurden nach dem optimalen Grenzwert in Hoch- und Niedrigausdrucksgruppen eingeteilt. n = 947 (TCGA-BRCA-Kohorte). Die blaue und grüne Kurve repräsentieren jeweils die Gruppen mit niedriger und hoher Expression. Die Einheit der mRNA-Expression war Fragmente pro Kilobase des Transkripts pro Million abgebildete Reads. Abkürzung: TCGA-BRCA = The Cancer Genome Atlas Breast Invasive Carcinom; m6A = N6-Methyladenosin. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 3. Korrelation zwischen aberranter mRNA-Expression der m6A-"Leser" HNRNPC, YTHDF1, YTHDF3 und IGF2BP1 mit der Gesamtüberlebensrate bei Brustkrebspatientinnen. (A) Eine hohe HNRNPC-mRNA-Expression war mit einer verringerten Gesamtüberlebensrate assoziiert (Log-Rank-Test, P = 0,01, <0,05). (B-D) Eine hohe Expression von YTHDF1, YTHDF3 und IGF2BP1 war signifikant mit einer schlechteren Prognose assoziiert (logaritarischer Test, P = 0,002, P = 0,014 bzw. P < 0,001). Die Patienten wurden nach dem optimalen Grenzwert in Hoch- und Niedrigausdrucksgruppen eingeteilt. n = 947 (TCGA-BRCA-Kohorte). Die blaue und grüne Kurve repräsentieren jeweils die Gruppen mit niedriger und hoher Expression. Die Einheit der Gen-mRNA-Expression war Fragmente pro Kilobase des Transkripts pro Million kartierte Reads. Abkürzungen: TCGA-BRCA = The Cancer Genome Atlas Breast Invasive Carcinom; m6A = N6-Methyladenosin. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 4. Das potenzielle molekulare regulatorische Netzwerk von YTHDF3. (A) Protein-Interaktionsnetzwerk, das durch Abfrage der STRING-Datenbank mit YTHDF3 als Köder erstellt wird (Konfidenzwert ≥ 0,4, maximale Interaktoren = 50). Kantenbreiten geben die Stärke der unterstützenden Beweise an. (B)Die Metascape-Analyse identifizierte 13 signifikante funktionale Cluster (P < 0,01, Mindestanzahl = 3, Anreicherungsfaktor > 1,5). Ontologische Quellen umfassten die Kyoto Encyclopedia of Genes and Genomes Pathways und Gene Ontology Biological Processes. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
| Parameter | HR | 95%-KI | p-Wert |
| YTHDF3 | 1.024 | 1.003-1.046 | 0.024 |
| Alter | 1.036 | 1.023-1.050 | <0,001 |
| Rennen | 1.275 | 0.892-1.821 | 0.182 |
| Bühne | 2.263 | 1.796-2.851 | <0,001 |
Tabelle 1: Ergebnisse der multivariaten Cox-Überlebensregressionsanalyse des m6A-bezogenen Gens YTHDF3. Abkürzungen: HR = Hazard Ratio; 95 % KI = Konfidenzintervall.
Ergänzende Abbildung S1. Die Protein- und mRNA-Spiegel von YTHDF3 in verschiedenen Subtypen von Brustkrebs und normalen Geweben stammen aus den CPTAC- und TCGA-Datensätzen. Die Proteinexpression von YTHDF3 in (A) primären Tumoren und normalen Geweben sowie (B) Brustinvasiven Karzinom-Subtypen. Z-Werte geben die Anzahl der Standardabweichungen vom Medianwert über die Proben hinweg für den entsprechenden Krebstyp an. Die Spektralzahlverhältniswerte von CPTAC log2 wurden zunächst innerhalb jedes Stichprobenprofils und dann über die Proben verteilt. (C) Die mRNA-Expression (RNA-seq, Fragmente pro Kilobase des Transkripts pro Million kartierte Lesungen) von YTHDF3 in verschiedenen Subtypen des brustinvasiven Karzinoms. Abkürzungen: CPTAC = das Klinische Konsortium für proteomische Tumoranalyse; TCGA = Der Cancer Genome Atlas. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S2. Promotor-Methylierungsstatus von m6A-bezogenen Genen, abgerufen aus der MethHC-Datenbank. (A) Promotormethylierung von m6A-"Schreibern" bei Brustkrebs und normalen Kontrollen. (B) Promotormethylierung von m6A-"Radiergummin" bei Brustkrebs und normalen Kontrollen. (C) Promotormethylierung von m6A-"Lesern" bei Brustkrebs und normalen Kontrollen. Rot, Brustkrebs; Grün, normale Steuerung. *P < 0,05, **P < 0,005. Abkürzung: m6A = N6-Methyladenosin. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S3. Genmutationen und Änderungen der Kopienzahl, analysiert von cBioPortal für m6A-"Schreiber" und "Radiergummis" (TCGA, Provisional, n = 1.108). (A) Änderungen an m6A-"Schriftstellern". (B) Änderungen an m6A-"Radiergummin". Balkendiagramme zeigen die Häufigkeit jeder Veränderungsart an: Fehlsense-Mutationen, Amplifikationen, tiefe Löschungen und mehrere Veränderungen. Nur Gene mit einer Veränderungshäufigkeit ≥ 1 % werden gezeigt. Die Gesamthäufigkeit der Veränderungen pro Gen wird rechts angezeigt. Abkürzungen: TCGA = The Cancer Genome Atlas; m6A = N6-Methyladenosin. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S4. Genmutation und Kopienzahl-Veränderungen, analysiert von cBioPortal für m6A-"Leser" (TCGA, Provisional, n = 1.108). Balkendiagramme zeigen die Häufigkeit jeder Veränderungsart an: Fehlsense-Mutationen, Amplifikationen, tiefe Löschungen und mehrere Veränderungen. Nur Gene mit einer Veränderungshäufigkeit ≥ 1 % werden gezeigt. Die Gesamthäufigkeit der Veränderungen pro Gen wird rechts angezeigt. Abkürzungen: TCGA = The Cancer Genome Atlas; m6A = N6-Methyladenosin. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S5. Korrelation von YTHDF3 mit dem Überleben von Brustkrebspatientinnen in KMplotter. (A) Eine hohe YTHDF3-mRNA-Expression war mit einer verringerten Gesamtüberlebensrate assoziiert (P = 0,011, <0,05). (B-D) Die hohe Expression von YTHDF3 war signifikant mit einem schlechteren rezidivfreien Überleben, fernem metastasfreiem Überleben und postprogressionärem Überleben verbunden (P = 1,9e-09, P = 0,06 bzw. P = 0,029). Alle P-Werte wurden durch den Log-Rank-Test berechnet. Die Anzahl der gefährdeten Patienten wird unter jedem Diagramm dargestellt. Der optimale Cut-off wurde durch die Auto-Select-Best-Cutoff-Funktion des Werkzeugs bestimmt. Abkürzung: KMplotter = Kaplan-Meier Plotter. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S6. Korrelation der YTHDF3-Expression mit der Immuninfiltration bei Brustkrebs. Die YTHDF3-Expression korrelierte positiv mit den Infiltrationsniveaus von B-Zellen, CD4+- T-Zellen, CD8+- T-Zellen, Makrophagen, Neutrophilen und DCs im TIMER-Datensatz. Die Tumorreinheit wurde als Kovariaten in der partiellen Spearman-Korrelationsanalyse verwendet. Abkürzungen: DCs = dendritische Zellen; TIMER = Ressource zur Tumorimmunabschätzung. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Abbildung S7. Quantitative Reverse-Transkriptions-Polymerase-Kettenreaktionsvalidierung der YTHDF3-mRNA-Expression in 20 gepaarten Brustkrebsgeweben und angrenzenden normalen Geweben. Die YTHDF3-mRNA-Spiegel wurden auf β-Aktin normalisiert. Die Daten werden als 2⁻ΔΔCt-Werte präsentiert. Der statistische Vergleich wurde mit einem zweiseitigen gepaarten Studenten-t-Test zu ΔC-t-Werten (P = 0,001) durchgeführt. Abkürzungen: NC = angrenzende, normale Brustgewebe; BrCa = Brustkrebsgewebe. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle S1. Zusammengefasste Ergebnisse der multivariaten Überlebensanalyse in m6A-bezogenen Genen. Abkürzungen: HR = Hazard Ratio; 95 % KI = Konfidenzintervall. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle S2. Die Annotations- und Anreicherungsergebnisse von 13 YTHDF3-bezogenen Genen basieren auf Metascape.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle S3. Klinische Informationen von 20 gepaarten Brustkrebsproben.Bitte klicken Sie hier, um diese Datei herunterzuladen.