Expressão diferencial de mRNA de genes relacionados a m6A em câncer de mama
A análise dos dados TCGA e GTEx (1.034 tumores, 104 NATs, 92 tecidos saudáveis) revelou que a maioria dos genes relacionados ao m6A foi expressa diferencialmente (Figura 1). Em relação aos tecidos mamários saudáveis, 10 genes foram significativamente regulados para cima: KIAA1429, RBM15, RBM15B, HNRNPC, YTHDF1, YTHDF2, YTHDF3, IGF2BP1, IGF2BP3 e ALKBH5 (todos P < 0,0001, exceto RBM15 P < 0,001). Oito genes foram regulados para baixo: WTAP, METTL3, METTL14, HNRNPA2B1, YTHDC1, YTHDC2, IGF2BP2 e FTO (P < 0,05 a P < 0,0001). METTL16 e eIF3A não mostraram diferença significativa. Ao comparar tumores com NATs, HNRNPA2B1 foi regulado em alta, enquanto YTHDF3, eIF3A e ALKBH5 foram regulados para baixo; RBM15, RBM15B e YTHDF2 não apresentaram diferença (Figura 1). O uso de NATs como controles foi questionado, possivelmente explicando as discrepâncias nas comparações de tecidossaudáveis 18,19. Curiosamente, uma tendência de maior abundância de proteína YTHDF3 em todos os subtipos de câncer de mama foi observada no conjunto de dados CPTAC em comparação com tecidos normais (P < 0,0001). Além disso, uma expressão significativamente maior de mRNA de YTHDF3 foi observada apenas em subtipos luminais de câncer de mama em comparação com tecidos normais (P = 0,004, < 0,01) (Figura Suplementar S1). No geral, uma desregulação significativa dos genes relacionados ao m6A é evidente no câncer de mama.
Metilação por promotores e amplificação gênica de genes relacionados a m6A em câncer de mama
A metilação do DNA desempenha um papel importante na regulação epigenética da transcrição gênica, e a hipermetilação do promotor inibe a expressão domRNA 20. A análise MethHC demonstrou diferenças significativas de metilação dos promotores apenas para IGF2BP1, IGF2BP2 e IGF2BP3 em tumores versus controles (P < 0,005). A regulação negativa de IGF2BP2 mRNA pode ser parcialmente atribuída à sua hipermetilação. A análise do cBioPortal mostrou que WTAP e ALKBH5 apresentaram regulação negativa do mRNA, enquanto KIAA1429 (VIRMA), YTHDF1, YTHDF3 e IGF2BP1 apresentaram amplificação gênica associada à sua regulação para cima (Figura Suplementar S2, Figura Suplementar S3 e Figura Suplementar S4). Esses achados sugerem que os níveis de mRNA dos genes relacionados ao m6A não são predominantemente impulsionados pela metilação do promotor ou amplificação gênica.
Análise de sobrevivência de genes relacionados ao m6A em câncer de mama
A análise de sobrevivência de 947 pacientes revelou que baixa expressão de RBM15B e alta expressão de METTL16, HNRNPC, YTHDF1, YTHDF3 e IGF2BP1 estavam significativamente correlacionadas com a redução da sobrevivência geral (Figura 2 e Figura 3). A regressão multivariada de Cox, ajustando por idade, estágio e raça, confirmou que a regulação para cima do YTHDF3 era um preditor independente de baixa sobrevivência geral (HR: 1,024, IC 95%: 1,003–1,046, P = 0,024) (Tabela 1). Para validação externa, a análise usando a ferramenta online KMplotter mostrou que uma alta expressão de YTHDF3 prevê um prognóstico ruim da doença (Figura Suplementar S5). Outros genes não foram significativos de forma independente (Tabela Suplementar S1).
Rede molecular potencial e infiltração imune do gene relacionado a m6A YTHDF3
A análise STRING do YTHDF3 (com interatores limitados a 50) identificou uma rede com 224 arestas e enriquecimento significativo (P < 1e-16). YTHDF3 está intimamente associado a METTL3, ALKBH5, WTAP, METTL14 e outros (Figura 4A). A análise de enriquecimento do Metascape identificou 13 agrupamentos funcionais significativos (P < 0,01, contagem mínima de 3, fator de enriquecimento > 1,5; Figura 4B). O cluster mais bem classificado foi a regulação metabólica do mRNA (GO1903311, 22 genes). Outros grupos relevantes para o câncer incluíram reparação de dealquilação do DNA (GO0006307), processamento de pri-miRNA (GO0031053) e estabilização de telomerase (GO1904356), entre outros. A lista completa de nomes de clusters, genes associados e estatísticas de enriquecimento está fornecida na Tabela Suplementar S2. Além disso, os resultados revelaram correlações positivas significativas entre a expressão de YTHDF3 e as células imunes (todas P < 0,01) (Figura Suplementar S6). Assim, a regulação para cima do YTHDF3 pode estar associada ao metabolismo do mRNA, reparo do DNA, manutenção dos telômeros e infiltração imunológica.
Expressão aumentada de mRNA de YTHDF3 em câncer de mama
Tecidos de câncer de mama e tecidos normais adjacentes foram coletados de 20 pacientes que passaram por cirurgia no Departamento de Cirurgia de Mama, o Hospital Afiliado da Southwest Medical University, entre julho de 2023 e outubro de 2023. A análise qRT-PCR dessas amostras pareadas mostrou expressão significativamente maior de mRNA YTHDF3 em tecidos tumorais do que em tecidos normais adjacentes (P = 0,001) (Figura Suplementar S7 e Tabela Suplementar S3).
Disponibilidade de Dados
Todos os conjuntos de dados analisados neste estudo estão disponíveis publicamente nos seguintes repositórios de acesso aberto. Expressão gênica e dados clínicos: Dados de expressão de RNA-seq do câncer de mama e informações clínicas correspondentes foram obtidos do The Cancer Genome Atlas (TCGA-BRCA) via The Human Protein Atlas (https://www.proteinatlas.org/). Dados adicionais de expressão normal de tecido mamário foram obtidos do projeto Genotype-Tissue Expression (GTEx) (acesso dbGaP phs000424.v8.p2) através do Breast Cancer Gene-Expression Miner v4.7 (http://bcgenex.ico.unicancer.fr/). A validação da sobrevivência foi realizada usando o Plotter Kaplan-Meier (https://kmplot.com). Os dados de expressão de proteínas e mRNAs foram obtidos via UALCAN (http://ualcan.path.uab.edu/). Alterações genéticas: Dados de mutação somática e alteração de número de cópia para o conjunto de dados invasivo de câncer de mama (TCGA, Provisório, 1108 amostras) foram acessados via cBioPortal (https://www.cbioportal.org/). Metilação do promotor de DNA: Os dados de metilação foram obtidos do MethHC (http://methhc.mbc.nctu.edu.tw/), que contém 839 amostras de câncer de mama e controles normais adjacentes parcialmente pareados. Interação proteína-proteína e enriquecimento de vias de vias: A análise funcional de redes foi realizada usando STRING v12.0 (https://string-db.org/, escore de confiança 0,4) e Metascape (http://metascape.org/). A análise de infiltração imune foi realizada usando TIMER (https://cistrome.shinyapps.io/timer/). Quaisquer dados adicionais subjacentes aos números estão disponíveis junto ao autor correspondente, mediante solicitação razoável.

Figura 1. Expressão diferencial de mRNA de genes relacionados à N6-metiladenosina (m6A) em câncer de mama. Vermelho indica regulação para cima; Azul indica regulação para baixo, preto indica ausência de expressão e branco indica ausência de dados de expressão gênica. A codificação por cores reflete significância estatística. O teste Dunnett-Tukey-Kramer foi usado para comparações par a par. Nota: "Normal" compreende a comparação entre tumores e NATs no TCGA, bem como a comparação entre tumores e tecidos saudáveis no GTEx. Abreviações: NATs = tecidos adjacentes normais; TCGA = Atlas do Genoma do Câncer; GTEx = Expressão de Genótipo-Tecido. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 2. Associação entre a expressão aberrante do mRNA dos "escritores" m6A RBM15B e METTL16 e a sobrevivência geral em pacientes com câncer de mama. (A) Baixa expressão de RBM15B esteve associada a uma sobrevivência geral menor (teste log-rank, P = 0,008, <0,01); (B) Alta expressão de METTL16 estava associada a prognóstico ruim (teste log-rank, P = 0,013, <0,05). Os pacientes foram dicotomiados em grupos de alta e baixa expressão usando o valor ótimo de corte. n = 947 (coorte TCGA-BRCA). As curvas azul e verde representam, respectivamente, os grupos de baixa e alta expressão. A unidade de expressão do mRNA foi Fragmentos por Quilobase de transcrito por Milhão de leituras mapeadas. Abreviação: TCGA-BRCA = Carcinoma Invasivo de Mama do Atlas do Genoma do Câncer; m6A = N6-metiladenosina. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 3. Correlação entre a expressão aberrante de mRNA dos "leitores" m6A HNRNPC, YTHDF1, YTHDF3 e IGF2BP1 com a sobrevivência geral em pacientes com câncer de mama. (A) Alta expressão de mRNA de HNRNPC esteve associada a uma sobrevivência geral reduzida (teste log-rank, P = 0,01, <0,05). (B-D) Alta expressão de YTHDF1, YTHDF3 e IGF2BP1 esteve significativamente associada a um prognóstico pior (teste log-rank, P = 0,002, P = 0,014 e P < 0,001, respectivamente). Os pacientes foram dicotomiados em grupos de alta e baixa expressão usando o valor ótimo de corte. n = 947 (coorte TCGA-BRCA). As curvas azul e verde representam, respectivamente, os grupos de baixa e alta expressão. A unidade de expressão do mRNA do gênio foi Fragmentos por Quilobase de transcrito por Milhão de leituras mapeadas. Abreviações: TCGA-BRCA = Carcinoma Invasivo de Mama do Atlas do Genoma do Câncer; m6A = N6-metiladenosina. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 4. A potencial rede molecular reguladora do YTHDF3. (A) Rede de interação proteica construída consultando o banco de dados STRING com YTHDF3 como isca (escore de confiança ≥ 0,4, interatores máximos = 50). Larguras das bordas indicam a força das evidências de apoio. (B)A análise Metascape identificou 13 clusters funcionais significativos (P < 0,01, contagem mínima = 3, fator de enriquecimento > 1,5). As fontes ontológicas incluíram a Enciclopédia de Genes e Genomas de Kyoto e os Processos Biológicos de Ontologia Genética. Por favor, clique aqui para ver uma versão ampliada desta figura.
| Parâmetros | RH | IC 95% | valor-p |
| YTHDF3 | 1.024 | 1.003-1.046 | 0.024 |
| Idade | 1.036 | 1.023-1.050 | <0,001 |
| Corrida | 1.275 | 0.892-1.821 | 0.182 |
| Palco | 2.263 | 1.796-2.851 | <0,001 |
Tabela 1: Resultados da análise de regressão multivariada de sobrevivência de Cox do gene relacionado a m6A YTHDF3. Abreviações: HR = razão de risco; IC 95% = intervalo de confiança.
Figura Suplementar S1. Os níveis de proteína e mRNA de YTHDF3 em diferentes subtipos de câncer de mama e tecidos normais dos conjuntos de dados CPTAC e TCGA. A expressão proteica de YTHDF3 em (A) tumores primários e tecidos normais e (B) subtipos de carcinoma invasivo da mama. Valores Z indicam o número de desvios padrão em relação ao valor mediano entre amostras para o tipo correspondente de câncer. Os valores da razão espectral de contagem do CPTAC log2 foram normalizados primeiro dentro de cada perfil amostral e depois entre as amostras. (C) A expressão de mRNA (RNA-seq, Fragmentos por Quilobase de transcrito por Milhão de leituras mapeadas) de YTHDF3 em diferentes subtipos de carcinoma invasivo de mama. Abreviações: CPTAC = Consórcio Clínico de Análise Proteômica de Tumores; TCGA = Atlas do Genoma do Câncer. Por favor, clique aqui para baixar este arquivo.
Figura Suplementar S2. Status de metilação do promotor de genes relacionados a m6A recuperado do banco de dados MethHC. (A) Metilação por promotores de "escritores" de m6A em câncer de mama e controles normais. (B) Metilação promotora de "borrachas" m6A em câncer de mama e controles normais. (C) Metilação promotora dos "leitores" de m6A em câncer de mama e controles normais. Vermelho, câncer de mama; Verde, controles normais. *P < 0,05, **P < 0,005. Abreviação: m6A = N6-metiladenosina. Por favor, clique aqui para baixar este arquivo.
Figura Suplementar S3. Mutação gênica e alterações no número de cópias analisadas pelo cBioPortal para "escritores" e "borrachas" m6A (TCGA, Provisório, n = 1.108). (A) Alterações em "escritores" m6A. (B) Alterações nas "borrachas" m6A. Gráficos de barras indicam a frequência de cada tipo de alteração: mutações missentido, amplificações, deleções profundas e múltiplas alterações. Apenas genes com frequência de alteração ≥ 1% são mostrados. A frequência geral de alteração por gene é exibida à direita. Abreviações: TCGA = Atlas do Genoma do Câncer; m6A = N6-metiladenosina. Por favor, clique aqui para baixar este arquivo.
Figura Suplementar S4. Mutação genética e alterações no número de cópias analisadas pelo cBioPortal para "leitores" m6A (TCGA, Provisório, n = 1.108). Gráficos de barras indicam a frequência de cada tipo de alteração: mutações missentido, amplificações, deleções profundas e múltiplas alterações. Apenas genes com frequência de alteração ≥ 1% são mostrados. A frequência geral de alteração por gene é exibida à direita. Abreviações: TCGA = Atlas do Genoma do Câncer; m6A = N6-metiladenosina. Por favor, clique aqui para baixar este arquivo.
Figura Suplementar S5. Correlação de YTHDF3 com a sobrevivência de pacientes com câncer de mama no KMplotter. (A) Alta expressão de mRNA YTHDF3 esteve associada a uma sobrevivência geral reduzida (P = 0,011, <0,05). (B-D) Alta expressão de YTHDF3 esteve significativamente associada a menor sobrevivência sem recidiva, sobrevivência livre de metástase à distância e sobrevivência pós-progressão (P = 1,9e-09, P = 0,06 e P = 0,029, respectivamente). Todos os valores P foram calculados pelo teste log-rank. O número de pacientes em risco é mostrado abaixo de cada gráfico. O corte ideal era determinado pelo recurso de melhor corte de seleção automática da ferramenta. Abreviação: KMplotter = Plotter de Kaplan-Meier. Por favor, clique aqui para baixar este arquivo.
Figura Suplementar S6. Correlação da expressão de YTHDF3 com infiltração imune no câncer de mama. A expressão de YTHDF3 correlacionou-se positivamente com os níveis de infiltração de células B, células T CD4+ , células T CD8+ , macrófagos, neutrófilos e DCs no conjunto de dados TIMER. A pureza tumoral foi usada como covariável na análise de correlação parcial de Spearman. Abreviações: DCs = células dendríticas; TIMER = Recurso de Estimativa do Imunismo Tumoral. Por favor, clique aqui para baixar este arquivo.
Figura suplementar S7. Validação quantitativa da reação em cadeia da polimerase por transcrição reversa da expressão do mRNA YTHDF3 em 20 tecidos pareados de câncer de mama e tecidos normais adjacentes. Os níveis de mRNA do YTHDF3 foram normalizados para β-actina. Os dados são apresentados como valores de 2⁻ΔΔCt. A comparação estatística foi realizada usando um teste t de Student em duas caudas pareadas nos valores deΔC t (P = 0,001). Abreviações: NC = tecidos mamários normais adjacentes; BrCa = tecidos do câncer de mama. Por favor, clique aqui para baixar este arquivo.
Tabela Suplementar S1. Resultados resumidos da análise multivariada de sobrevivência em genes relacionados a m6A. Abreviações : HR = razão de risco; IC 95% = intervalo de confiança. Por favor, clique aqui para baixar este arquivo.
Tabela Suplementar S2. Os resultados de anotação e enriquecimento de 13 genes relacionados ao YTHDF3 baseados no Metascape.Por favor, clique aqui para baixar este arquivo.
Tabela Suplementar S3. Informações clínicas de 20 amostras pareadas de câncer de mama.Por favor, clique aqui para baixar este arquivo.