Este estudo foi aprovado pelo Comitê de Ética do Hospital Provincial de Shaanxi (Número de Aprovação 2024-R055) e foi conduzido de acordo com os princípios da Declaração de Helsinque. O consentimento informado por escrito foi obtido de todos os participantes antes da coleta de qualquer amostra sanguínea. Todas as amostras humanas de sangue e soro foram coletadas, manipuladas e armazenadas em conformidade com as diretrizes institucionais do Departamento de Medicina Laboratorial.
Todas as amostras de sangue venoso e soro foram tratadas como material biológico potencialmente infeccioso. O processamento das amostras foi realizado em laboratório de biossegurança nível 2 por pessoal treinado utilizando equipamento de proteção individual apropriado. Material plástico descartável, ponteiras de pipeta e microplacas usadas que tiveram contato com material biológico foram autoclavados ou coletados como resíduo biológico perigoso e descartados por meio do fluxo institucional de resíduos médicos licenciado, conforme as regulamentações locais.
Visão geral do desenho do estudo:
O fluxo de trabalho deste estudo é resumido em um esquema de desenho do estudo (Figura Suplementar S1) e é dividido em dois braços complementares. No braço in silico, conjuntos de dados transcriptômicos públicos, de célula única, proteicos, de modificações pós-traducionais, imunes e de sensibilidade a fármacos foram analisados para caracterizar o contexto biológico da hWAPL e gerar hipóteses mecanicistas. No braço experimental, a hWAPL sérica e o SCC-Ag foram medidos em uma coorte independente de centro único, e seu desempenho diagnóstico, isoladamente e em combinação, foi avaliado por análise da curva ROC e regressão logística. Os dois braços foram integrados apenas ao nível da interpretação, sendo o braço sérico a única validação experimental.
Perfilagem da expressão gênica em câncer panfocal e carcinoma de células escamosas do colo do útero
Dados de sequenciamento de RNA provenientes do The Cancer Genome Atlas (TCGA) foram obtidos por meio do UCSC Xena Browser. Todos os bancos de dados públicos e ferramentas web utilizados neste estudo foram acessados entre setembro e outubro de 2025, e as versões específicas e datas de acesso estão listadas na Tabela de Materiais. Os níveis de expressão de WAPL, quantificados como fragmentos por quilobase de transcrito por milhão de leituras mapeadas (FPKM), foram submetidos a uma transformação log₂ com a adição de um pseudocount de 1. Esses valores transformados foram então utilizados para realizar análises comparativas entre tecidos malignos e amostras normais adjacentes pareadas em múltiplos tipos de câncer. Para a investigação focada em malignidades cervicais, foram isolados os dados de expressão de WAPL específicos ao carcinoma de células escamosas do colo do útero e adenocarcinoma endocervical (CESC). O perfil de expressão nessas amostras tumorais foi posteriormente visualizado e comparado com os níveis de expressão em tecidos de controle epiteliais cervicais normais.
Análise dos isoformas do transcrito de WAPL
Os perfis de expressão dos isoformas do transcrito de WAPL foram examinados utilizando o banco de dados GEPIA2. Os isoformas específicos foram delimitados por seus identificadores de transcritos Ensembl correspondentes (ENST). Os níveis de expressão foram calculados como o logaritmo na base 2 dos transcritos por milhão mais um [log₂(TPM+1)] e são apresentados graficamente como gráficos de violino.
Análise de sequenciamento de RNA de célula única
O conjunto de dados de scRNA-seq GSE168652 foi obtido do repositório Gene Expression Omnibus (GEO). O pré-processamento dos dados, avaliação da qualidade, normalização e agrupamento foram realizados com a plataforma scCancerExplorer. Foram filtradas as células que expressavam menos de 200 genes e os genes detectados em menos de três células. A normalização foi realizada utilizando a abordagem de Log-Normalização com um fator de escala de 10.000.
Os 2.000 genes mais variáveis foram identificados utilizando transformação de estabilização de variância e posteriormente utilizados para análise de componentes principais (PCA). A Aproximação e Projeção Uniformes de Variedade (UMAP) para redução de dimensionalidade foi aplicada com base nos primeiros 30 componentes principais. O agrupamento baseado em grafos foi implementado por meio do algoritmo Louvain com uma resolução de 0,8. Os agrupamentos celulares foram anotados de acordo com genes marcadores estabelecidos. Os níveis de expressão de WAPL foram medidos em diferentes populações celulares, e as fases do ciclo celular foram determinadas utilizando o método CellCycleScoring.
Análise da expressão proteica e de modificações pós-traducionais
A expressão e distribuição subcelular de WAPL no CSCC foram avaliadas utilizando o banco de dados Human Protein Atlas. Os sítios de modificação pós-traducional (PTM) foram identificados usando o PhosphoSitePlus. Além disso, foram investigadas associações potenciais entre os níveis de expressão de WAPL e genes implicados em processos de modificação proteica, especificamente AKT1, AKT2, KAT2B, EP300, USP7 e KAT5, no CSCC, utilizando a ferramenta Xiantao Academic13,14,15,16,17.
Rede de coexpressão e análise de enriquecimento funcional
Os dados de sequenciamento de RNA para TCGA-CSCC, processados usando o alinhador STAR e normalizados em transcritos por milhão (TPM), foram obtidos do portal Genomic Data Commons (GDC). Para identificar genes coexpressos com WAPL, foi aplicado um limiar de triagem com base em um coeficiente de correlação de Pearson absoluto maior que 0,6 e um nível de significância estatística de P < 0,05.
Para investigar as funções biológicas associadas à WAPL, foi realizada inicialmente uma análise de expressão diferencial após estratificar os pacientes em grupos de alta e baixa expressão de WAPL, utilizando o valor mediano de expressão como ponto de corte. Os genes diferencialmente expressos foram identificados com a plataforma Xiantao Academic e posteriormente submetidos à análise de enriquecimento da Ontologia Genética, incluindo as categorias de processo biológico, componente celular e função molecular, juntamente com a análise de vias da Enciclopédia de Genes e Genomas de Kyoto. Independentemente, as interações proteicas envolvendo WAPL foram exploradas utilizando o STRING, e a rede de interação resultante foi reconstruída e analisada no Cytoscape.
A análise de enriquecimento de conjuntos de genes (GSEA) foi realizada utilizando as bases de dados Reactome, WikiPathways e KEGG, com critérios de triagem definidos em um escore de enriquecimento normalizado absoluto (|NES|) > 1 e uma taxa de falsa descoberta (FDR) < 0,25.
Análise da infiltração imune
As ferramentas computacionais TIMER 2.0 e xCell foram empregadas para analisar dados transcriptômicos do cohort TCGA-CESC, permitindo a avaliação quantitativa do estado de infiltração de 24 populações distintas de células imunes.
Previsão da sensibilidade a medicamentos
As previsões de sensibilidade a medicamentos foram obtidas a partir do banco de dados GDSC2 utilizando o pacote R oncoPredict. A concentração inibitória média (IC50) foi calculada para um painel de 86 compostos quimioterápicos.
Sobrevivência e análise diagnóstica
As curvas de sobrevida global (OS) e intervalo livre de progressão (PFI) foram construídas utilizando a plataforma analítica Xiantao, com as coortes de pacientes dicotomizadas no nível mediano de expressão de WAPL. A significância prognóstica das variáveis foi avaliada por meio de modelos de regressão de riscos proporcionais de Cox univariados e multivariados. Para avaliar a capacidade diagnóstica de WAPL e SCC (SERPINB3) em distinguir CSCC do tecido normal adjacente, foi aplicada análise da curva característica de operação do receptor (ROC) a dados transcriptômicos, expressos como transcritos por milhão (TPM), obtidos do The Cancer Genome Atlas.
Coleta de amostras de soro e ELISA
Este estudo foi aprovado pelo Comitê de Ética em Pesquisa do Hospital Provincial de Shaanxi (Número de Aprovação 2024-R055) e foi conduzido em conformidade com a Declaração de Helsinque. Um total de 89 espécimes de soro foram obtidos de pacientes com CSCC confirmado histopatologicamente, juntamente com 89 controles saudáveis do sexo feminino pareados por idade. A idade mediana foi de 54 anos, variando de 31 a 77 anos, no grupo com CSCC, e de 41 anos, variando de 30 a 52 anos, no grupo controle. Todos os casos de câncer cervical foram recém-diagnosticados e sem tratamento prévio no momento da coleta sanguínea; as amostras foram obtidas antes do início de qualquer cirurgia, quimioterapia ou radioterapia. Observamos que, embora os dois grupos tenham sido recrutados em um contexto clínico comparável, as idades medianas diferiram em 13 anos. Essa diferença é analisada diretamente na seção de Resultados e em uma análise de sensibilidade específica ajustada por idade, sendo reconhecida como uma limitação na seção de Discussão.
O sangue venoso foi coletado em tubos sem anticoagulante, deixado coagular à temperatura ambiente por 10–20 min e, posteriormente, centrifugado a 2.000–3.000 × g durante 20 min. O soro foi separado dentro de 1 h após a coleta, fracionado imediatamente e armazenado a −20 °C; o intervalo entre a punção venosa e o congelamento não excedeu 2 h para nenhuma amostra. As alíquotas passaram por um único ciclo de congelamento–descongelamento antes da análise, e todas as medidas foram concluídas dentro de três meses após o armazenamento. O soro resultante foi dividido em alíquotas e conservado a −20 °C. Quaisquer amostras hemolisadas ou lipêmicas foram excluídas da análise.
A quantificação dos níveis séricos da proteína hWAPL foi realizada utilizando um kit comercialmente disponível de ensaio imunoenzimático (ELISA) baseado no princípio do sanduíche com duplo anticorpo. Resumidamente, padrões de calibração e amostras-teste foram dispensados em duplicata em poços de microplaca previamente revestidos com um anticorpo de captura específico para hWAPL, seguidos por uma incubação de 30 min a 37 °C. Após uma série de etapas de lavagem para remover material não ligado, um anticorpo de detecção marcado com peroxidase do rabo-de-cavalo (HRP) e direcionado contra hWAPL foi adicionado e incubado em condições idênticas.
Após lavagem adicional, 50 µL de cada solução cromogênica de substrato A e B foram combinados nos poços e deixados desenvolver no escuro a 37 °C por 10 min. Solução de parada (50 µL) foi adicionada para encerrar a reação enzimática. A densidade óptica foi medida em um comprimento de onda primário de 450 nm, com leitura de referência em 630 nm, dentro de 15 min após a interrupção da reação. Uma curva padrão foi construída utilizando calibradores diluídos em série, abrangendo uma faixa de concentração de 10–450 pg/mL. A concentração de hWAPL para cada amostra desconhecida foi determinada por interpolação a partir dessa curva padrão e ajustada conforme qualquer diluição pré-analítica. O ensaio foi realizado rigorosamente de acordo com as instruções do fabricante. De acordo com as especificações do fabricante, a faixa de detecção funcional do kit é de 10–450 pg/mL. As amostras foram analisadas com uma diluição final de 5 vezes, e espécimes cuja densidade óptica excedesse a do calibrador mais alto foram diluídos novamente e reensaiados, aplicando-se o fator total de diluição durante o cálculo inverso. O desempenho analítico foi verificado na matriz sérica antes da análise das amostras do estudo: um pool sérico de alta concentração diluído em série apresentou linearidade aceitável ao longo da faixa de calibração, e a recuperação do calibrador recombinante em três pools séricos esteve dentro de uma janela aceitável, compatível com o design do kit de sanduíche com duplo anticorpo. O limite inferior de quantificação correspondeu ao calibrador mais baixo (10 pg/mL). Ressaltamos que este kit é rotulado pelo fabricante exclusivamente para uso em pesquisa e não constitui um dispositivo diagnóstico in vitro clinicamente validado; as concentrações absolutas e o ponto de corte derivado devem, portanto, ser interpretados como específicos da plataforma e provisórios, conforme discutido abaixo.
Concomitantemente, os níveis séricos de AgCC foram quantificados utilizando a plataforma Abbott Alinity i por meio de um ensaio imunométrico com micropartículas quimioluminescentes (CMIA), com limite superior de referência de 1,5 ng/mL. Todos os procedimentos analíticos foram realizados em duplicata por pessoal de laboratório que desconhecia o estado da doença e todas as variáveis clínicas, incluindo estágio clínico, grau histológico e status do HPV. Os coeficientes de variação interensaio e intraensaio para as metodologias empregadas foram mantidos abaixo de 10%.
Análise estatística
Como as concentrações séricas de hWAPL não atenderam ao pressuposto de normalidade de acordo com o teste de Shapiro–Wilk (todos os P < 0,05), as análises subsequentes foram realizadas utilizando métodos não paramétricos. Assim, as variáveis contínuas são apresentadas como medianas com intervalos interquartis (IIQs). As diferenças entre dois grupos independentes foram avaliadas utilizando o teste U de Mann–Whitney, enquanto comparações entre três ou mais grupos foram realizadas com o teste H de Kruskal–Wallis. As correlações entre variáveis foram avaliadas utilizando a análise de correlação de postos de Spearman para dados contínuos ou ordinais e o coeficiente de correlação pontual bisserial para variáveis binárias, como o status de infecção pelo HPV. Os desfechos de sobrevida foram comparados utilizando o teste log-rank, e as razões de risco com intervalos de confiança de 95% foram estimadas a partir de modelos de regressão de riscos proporcionais de Cox.
O desempenho diagnóstico foi avaliado por meio da análise da curva ROC, sendo o ponto de corte ótimo para cada biomarcador determinado pela maximização do índice de Youden. As diferenças na área sob a curva (AUC) foram testadas utilizando o método de DeLong. Para investigar a utilidade diagnóstica combinada, foi construído um modelo de regressão logística binária com o status da doença (CEEC = 1; controle = 0) como variável dependente e os níveis séricos de hWAPL e antígeno de CEC como preditores independentes. As probabilidades preditas resultantes foram salvas como um escore composto para análise subseqüente por curva ROC.
A sensibilidade diferencial aos fármacos entre os grupos com alta e baixa expressão de WAPL foi avaliada com o pacote limma, aplicando limiares de valor de P ajustado < 0,05 e |mudança no log₂| > 0,3. A significância estatística foi definida como valor de P bicaudal < 0,05. Todas as análises foram realizadas utilizando SPSS, R ou as ferramentas online especificadas. Os pacotes do R incluíram survival, survminer, pROC, ggplot2, oncoPredict e limma. Para todas as análises de alta dimensão nas quais muitas hipóteses foram testadas simultaneamente, ou seja, comparações pan-cancro, triagem de coexpressão, análises de expressão diferencial e enriquecimento funcional, e análise de sensibilidade a fármacos, os valores de P foram ajustados para comparações múltiplas utilizando a taxa de falsa descoberta de Benjamini–Hochberg, sendo os valores de P ajustados apresentados quando aplicável; as correlações par a par entre hWAPL e genes individuais ou frações de células imunes são apresentadas com valores de P nominais, e esses resultados correlacionais são considerados geradores de hipóteses, e não confirmatórios. Como alguns campos clinicopatológicos na coorte pública apresentavam preenchimento incompleto, os casos com dados ausentes foram excluídos da análise correspondente de forma par a par, em vez de serem imputados, e o número de casos avaliáveis é indicado para cada variável na Tabela 1. Para abordar a diferença de idade entre as coortes séricas, foi realizada uma análise de sensibilidade utilizando regressão logística binária ajustada para idade, e a correlação de Spearman entre hWAPL sérico e idade foi calculada dentro de cada grupo.