Este estudo foi conduzido de acordo com a Declaração de Helsinque. A aprovação foi concedida pelo Comitê de Ética da Universidade Médica da Mongólia Interior em 21 de outubro de 2021. O consentimento informado por escrito foi obtido de todos os participantes incluídos no estudo.
Coleta de espécimes clínicos
Tecidos de CCR pareados e tecidos não cancerosos adjacentes foram coletados de 45 pacientes submetidos à ressecção radical entre dezembro de 2017 e fevereiro de 2020. Nenhum dos pacientes recebeu radioterapia ou quimioterapia pré-operatória, e pacientes com outras neoplasias foram excluídos. O tecido não canceroso adjacente foi coletado a pelo menos 5 cm da margem tumoral. As amostras foram congeladas imediatamente após a remoção cirúrgica e armazenadas em nitrogênio líquido até a extração de RNA. A coorte incluiu 27 pacientes do sexo masculino e 18 do sexo feminino; 13 pacientes tinham menos de 60 anos, e 32 tinham 60 anos ou mais. As variáveis clínico-patológicas incluíram diâmetro do tumor, estágio TNM, grau de diferenciação, trombo tumoral vascular, metástase linfonodal e nível pré-operatório de CEA.
Aquisição e pré-processamento de dados públicos
Dados de miRNA-seq, RNA-seq e dados clínicos de acesso aberto para TCGA-COAD e TCGA-READ foram obtidos no portal Genomic Data Commons. Tumores primários e amostras de tecido sólido normal foram identificados de acordo com o tipo de amostra do GDC. Amostras de tumores metastáticos e recorrentes foram excluídas da análise de expressão diferencial entre tumor e tecido normal. Para os arquivos de miRNA-seq do TCGA, a entrada ao nível da base de dados detectada nesta análise foi hsa-mir-192; os valores de expressão foram analisados como log2(RPM + 1), em que RPM indica leituras por milhão de leituras de miRNA mapeadas. Foram realizadas análises não pareadas e pareadas entre tumor e tecido normal quando amostras emparelhadas estavam disponíveis. Os conjuntos de dados GEO GSE89076 e GSE156355 foram utilizados para triagem de genes-alvo candidatos. Genes diferencialmente expressos provenientes do GEO foram triados utilizando valores ajustados p < 0,05 e mudança absoluta no log2 de pelo menos 2. Os nomes dos conjuntos de dados foram verificados quanto à consistência ao longo do manuscrito.
Análises diagnósticas, clínico-patológicas e de sobrevida
A análise da curva ROC (característica operacional do receptor) foi utilizada para avaliar a separação exploratória entre tumor e tecido normal nos conjuntos de dados TCGA-COAD, TCGA-READ e no coorte combinado COADREAD. Como o número de amostras normais era pequeno, foram realizadas análises de sensibilidade adicionais, incluindo reamostragem balanceada entre tumor e tecido normal com 5.000 iterações e análise leave-one-normal-out. Para as análises de associação clinicopatológica, as amostras de tumor primário foram divididas em grupos de alta e baixa expressão com base no nível mediano de expressão da miR-192 no TCGA. As análises foram repetidas para o coorte combinado COADREAD e separadamente para COAD e READ, sempre que o tamanho da amostra permitiu. A sobrevida global foi avaliada por meio da análise de Kaplan-Meier com o teste log-rank e regressão de Cox univariável. As análises diagnósticas, clinicopatológicas e de sobrevida foram relatadas separadamente para evitar a superestimação de conclusões prognósticas.
Análises de genes-alvo candidatos e genes centrais
Genes-alvo candidatos foram identificados integrando a predição do miRWalk com genes diferencialmente expressos de GSE89076 e GSE156355. A interseção entre essas três fontes foi mantida para análises exploratórias subsequentes. Uma rede de interação proteína-proteína foi construída utilizando o STRING, e os genes centrais foram selecionados com base no grau do nó e na conectividade do módulo por meio das análises cytoHubba e MCODE baseadas no Cytoscape. Os genes centrais mantidos (KIF20A, TPX2, CDCA5, CCNB1, CDK1, PLP1, NRXN1, GRIK3 e KIF5C) foram avaliados posteriormente utilizando dados de RNA-seq do TCGA. As diferenças de expressão entre tumor e tecido normal, as correlações de Spearman com a expressão do miR-192 do TCGA e as associações com sobrevida geral foram resumidas como análises exploratórias suplementares. Esses genes foram considerados apenas como alvos candidatos; nenhuma ligação direta, efeito de resgate ou regulação em nível proteico downstream foi inferida a partir dos dados bioinformáticos.
Análise de enriquecimento de GO e KEGG
Foram realizadas análises de enriquecimento de Gene Ontology (GO) e Kyoto Encyclopedia of Genes and Genomes (KEGG) para o conjunto de genes candidatos e os módulos significativos da rede. Os termos de GO foram agrupados em categorias de processo biológico, componente celular e função molecular. Resultados de enriquecimento com p < 0,05 foram considerados estatisticamente relevantes e foram interpretados como sinais de vias geradores de hipóteses, e não como mecanismos validados experimentalmente.
Extração de RNA e análise de qRT-PCR
O RNA total foi extraído de 45 amostras pareadas de tecido de CCR e tecido não canceroso adjacente utilizando o reagente TRIzol, conforme o protocolo do fabricante. A pureza e a concentração do RNA foram avaliadas espectrofotometricamente, e amostras com valores de A260/A280 entre 1,8 e 2,1 foram utilizadas para a transcrição reversa. A qRT-PCR foi realizada utilizando um sistema de transcrição reversa para miRNA e reagentes de qPCR para miRNA baseados em SYBR Green da TransGen Biotech. A transcrição reversa e a amplificação por qPCR foram realizadas conforme as instruções do fabricante dos reagentes. Cada reação de qRT-PCR foi executada em triplicata técnica. O RNU6B foi utilizado como referência interna. Os primers para miR-192-5p foram 5'-GGGGCTGACCTATGAATTGA-3' (sentido) e 5'-CAGTGCAGGGTCCGAGGT-3' (anti-sentido). Os primers para RNU6B foram 5'-ATTGGAACGATACAGAGAAGATT-3' (sentido) e 5'-GGAACGCTTCACGAATTTG-3' (anti-sentido). A análise da curva de dissociação foi utilizada para avaliar a especificidade da amplificação, quando suportada pela plataforma de qPCR. A expressão relativa foi calculada utilizando o método 2-ΔΔCt. Os resultados de tecidos tumorais pareados e tecidos adjacentes foram analisados como medições clínicas pareadas.
Cultura celular e transfeção
A linhagem celular de CCR humano HT29 foi obtida do estoque celular do laboratório dos autores e havia sido autenticada e confirmada como livre de contaminação por micoplasma. As células HT29 foram cultivadas em meio McCoy's 5A suplementado com 10% de soro bovino fetal, 100 U/mL de penicilina e 100 µg/mL de estreptomicina a 37 oC em atmosfera com 5% de CO2. As células foram repicadas quando atingiram cerca de 80% de confluência. Os miméticos de miR-192-5p e o oligonucleotídeo controle negativo con238 foram adquiridos da Sangon Biotech (Xangai, China) e utilizados para transfeção transitória. Células HT29 foram semeadas em placas de 6 poços e transfectadas quando a densidade celular atingiu aproximadamente 20–30%, utilizando uma concentração final de oligonucleotídeo de 50 nM e um protocolo de transfeção baseado em lipídeos, conforme as instruções do fabricante do reagente. A proporção exata entre reagente e oligonucleotídeo, bem como a quantidade de reagente, seguiu o protocolo recomendado pelo fabricante, conforme registrado nos arquivos do laboratório dos autores. As células foram incubadas com os complexos de transfeção em condições padrão de cultivo, e análises downstream de imagem da transfeção e ensaios de cicatrização (wound-healing) foram realizados após a transfeção, de acordo com o cronograma experimental. Apenas os dados obtidos com células HT29 foram incluídos no presente estudo.
Ensaio de cicatrização de ferida
Células HT29 transfectadas foram semeadas em placas de 6 poços e cultivadas até a formação de uma monocamada confluente. Um arranhão linear foi criado utilizando uma ponteira estéril de 200 µL, e as células destacadas foram removidas por lavagem com PBS. O meio foi então substituído por meio sem soro para reduzir os efeitos relacionados à proliferação durante a avaliação da migração. Imagens da ferida foram capturadas em 0 h, 24 h e 48 h. O ensaio de cicatrização da ferida foi realizado em três réplicas biológicas independentes, e múltiplos campos microscópicos foram medidos para cada grupo e ponto temporal. A largura da ferida foi quantificada a partir das imagens de microscopia utilizando software de análise de imagem; a versão exata do software utilizada nas medições originais não estava disponível nos registros fonte. A taxa de cicatrização da ferida foi calculada como [(largura da ferida em 0 h – largura da ferida no ponto temporal indicado) / largura da ferida em 0 h] × 100%. Como os dados disponíveis do manuscrito não incluem um experimento controle de proliferação separado, os resultados do ensaio de cicatrização da ferida foram interpretados como relacionados à migração, mas sem exclusão completa dos efeitos de proliferação.
Análise de correlação de assinatura imune
A associação entre a expressão de miR-192 do TCGA e assinaturas relacionadas ao sistema imunológico foi reavaliada utilizando dados de RNA-seq de tumores primários. Um total de 425 amostras tumorais com dados pareados de miRNA e RNA-seq foi incluído na análise de correlação de escores de assinatura. A expressão de miR-192 mostrou correlações negativas estatisticamente significativas, mas fracas, com várias assinaturas imunes após correção da taxa de falsa descoberta, incluindo células Tfh (rho = -0.205, FDR = 0.00049), células dendríticas ativadas (rho = -0.186, FDR = 0.00098), células Treg (rho = -0.185, FDR = 0.00098), células dendríticas plasmocitoides (rho = -0.169, FDR = 0.0027) e macrófagos (rho = -0.149, FDR = 0.0092). Correlações positivas foram observadas principalmente para assinaturas de eosinófilos (rho = 0.139, FDR = 0.0155) e de células Th17 (rho = 0.121, FDR = 0.0296). Esses achados indicam associações imunes heterogêneas e geralmente fracas, ao invés de um padrão uniforme de infiltração imune positiva.
Análise estatística
As análises estatísticas foram realizadas utilizando o GraphPad Prism 9.0, SPSS 25.0 e software R 3.6.3. As variáveis contínuas foram resumidas como média ± desvio padrão ou mediana com intervalo interquartílico, conforme apropriado. Comparações pareadas de tecidos foram analisadas usando testes pareados quando as amostras de tumor e tecido adjacente eram emparelhadas. Testes não paramétricos foram utilizados quando as suposições de normalidade não foram atendidas. Curvas ROC, reamostragem balanceada e análise leave-one-normal-out foram utilizadas para discriminação diagnóstica exploratória; testes qui-quadrado foram utilizados para associações clínico-patológicas categóricas; modelos log-rank e de Cox foram utilizados para análises de sobrevida; e a correlação de Spearman foi utilizada para análises imunes e de genes centrais. Valores de p bilaterais < 0,05 foram considerados estatisticamente significativos, com correção da taxa de falsas descobertas relatada para correlações de assinaturas imunes.