Lo studio è stato esaminato e approvato dal Comitato Etico dell'Ospedale Popolare di Hangzhou affiliato alla Scuola di Medicina dell'Università Westlake. Il numero di approvazione era IIT-20230528-0108-01. Il consenso informato scritto è stato ottenuto dai partecipanti arruolati in modo prospettico prima del prelievo del sangue. Per le cartelle cliniche incluse in modo retrospettivo e i campioni residui archiviati raccolti durante le cure di routine, il consenso informato è stato esonerato dal comitato etico poiché sono stati utilizzati dati anonimizzati e non è stata effettuata alcuna ulteriore intervento.
Riproducibilità, controllo dei bias e integrità dei dati
I criteri di idoneità, i tempi di prelievo del sangue, le regole di valutazione dell'imaging, le definizioni dei sottogruppi, i predittori candidati, le definizioni delle soglie, le regole per i dati mancanti e le procedure di convalida sono stati specificati in anticipo prima dell'analisi statistica. I campioni di sangue sono stati raccolti prima della conferma tramite imaging con mezzo di contrasto e prima di qualsiasi trattamento antitumorale. I test di laboratorio hanno utilizzato identificatori dello studio anonimizzati e i campioni dei gruppi di addestramento e di convalida sono stati analizzati in ordine casuale. Il personale di laboratorio era all'oscuro della diagnosi finale, dell'assegnazione al gruppo, dello stadio secondo la classificazione del Barcelona Clinic Liver Cancer (BCLC) e dello stato di negatività per AFP. La diagnosi e lo stadio radiologico sono stati determinati indipendentemente rispetto ai test di laboratorio. L'inserimento dei dati, la conversione delle unità, l'assegnazione ai gruppi, la codifica delle esclusioni, la codifica dei valori mancanti e i set di dati statistici sono stati verificati da due ricercatori prima dell'adattamento del modello. I gruppi di addestramento e di convalida sono stati mantenuti separati durante tutta la selezione delle variabili, l'adattamento del modello, la derivazione delle soglie e la convalida.
Popolazione oggetto dello studio e disegno dello studio
Questo studio diagnostico a doppia coorte ha sviluppato e convalidato un modello multivariabile per il rilevamento del carcinoma epatocellulare (HCC) utilizzando la protrombina anomala/proteina indotta da PIVKA-II, AFP e variabili cliniche di routine. La coorte di addestramento assemblata in modo retrospettivo includeva partecipanti valutati tra il 1° gennaio 2021 e il 31 dicembre 2023. Questa coorte è stata utilizzata per la selezione delle variabili, l'adattamento del modello, la derivazione delle soglie e la convalida interna tramite bootstrap. La coorte di convalida arruolata in modo prospettico includeva partecipanti valutati tra il 1° gennaio 2024 e il 31 dicembre 2025. Questa coorte è stata utilizzata esclusivamente per test su campioni esterni. Non è stata effettuata alcuna nuova selezione di variabili, riadattamento del modello, aggiornamento dei coefficienti, riottimizzazione delle soglie o ricalibrazione basata su sottogruppi nella coorte di convalida. Il flusso dei partecipanti, l'assegnazione alle coorti, il prelievo di sangue, l'elaborazione del plasma, i test sui biomarcatori, la valutazione delle immagini e il flusso di lavoro statistico sono mostrati in Figura 1.
I partecipanti sono stati reclutati da un percorso di sorveglianza o diagnostico per epatopatie ad alto rischio presso l'Ospedale Affiliato Hangzhou First People's Hospital, Scuola di Medicina, Westlake University. Lo stato di alto rischio è stato definito come infezione da epatite B cronica, infezione da epatite C cronica, cirrosi di qualsiasi eziologia o steatosi epatica con fibrosi documentata in base a valutazione clinica, referti di laboratorio o precedenti esami di imaging. I partecipanti idonei erano adulti di età ≥18 anni che avevano completato il prelievo ematico pre-imaging e una tomografia computerizzata (TC) multiphasica con mezzo di contrasto e/o una risonanza magnetica (RM) dinamica con mezzo di contrasto entro 14 giorni dal prelievo del sangue. I criteri di esclusione comprendevano trattamenti precedenti per carcinoma epatocellulare (HCC), inclusi resezione, ablazione, terapia transarteriosa o terapia sistemica; somministrazione di vitamina K nei 7 giorni precedenti il prelievo; terapia con antagonisti della vitamina K o altre terapie anticoagulanti che potessero influenzare l'interpretazione della coagulazione; colangite acuta o ostruzione biliare al momento del prelievo; gravidanza; dati mancanti relativi ad AFP, PIVKA-II o diagnosi finale tramite imaging; e qualità inaccettabile del campione senza un campione ripetuto valido.
La popolazione finale dello studio comprendeva 389 partecipanti. Il gruppo di addestramento includeva 205 partecipanti, di cui 112 casi di carcinoma epatocellulare confermati mediante imaging e 93 controlli senza carcinoma epatocellulare. Il gruppo di validazione includeva 184 partecipanti, di cui 58 casi di carcinoma epatocellulare confermati mediante imaging e 126 controlli senza carcinoma epatocellulare. Nel gruppo di addestramento, i 93 controlli senza carcinoma epatocellulare comprendevano 30 controlli sani e 63 controlli con malattia epatica cronica. Nel gruppo di validazione, i 126 controlli senza carcinoma epatocellulare comprendevano 0 controlli sani e 126 controlli con malattia epatica cronica.
Lo stato sintomatico non è stato assegnato sulla base della sola anamnesi di epatite o di cirrosi compensata. Un partecipante è stato classificato come sintomatico solo quando il referto precedente all'imaging ha documentato un dolore nuovo o peggiorato nel quadrante superiore destro, perdita di peso inspiegabile, distensione addominale, ittero, urine scure o una scompensazione clinicamente documentata non spiegata da un'altra causa confermata. I partecipanti senza questi reperti sono stati classificati come asintomatici o con lesioni rilevate durante il monitoraggio.
Adjudicazione diagnostica e stadiazione
Lo stato finale di carcinoma epatocellulare (HCC) è stato determinato mediante tomografia computerizzata multiphasica con mezzo di contrasto e/o risonanza magnetica dinamica con mezzo di contrasto eseguite nell'ambito delle cure abituali. L'imaging doveva includere almeno una fase arteriosa e una fase porto-venosa; i reperti della fase ritardata sono stati inclusi quando disponibili. Una lesione è stata classificata come HCC quando, in un fegato a rischio, l'imaging ha mostrato un'iper-rappresentazione in fase arteriosa con washout in fase porto-venosa e/o ritardata, con o senza capsula di enhancement. I partecipanti senza caratteristiche radiologiche compatibili con HCC e senza una diagnosi maligna alternativa sono stati classificati come non-HCC.
Le informazioni di imaging sono state estratte mediante un modulo predeterminato che registrava la modalità di imaging, le fasi ottenute, l'iper-rappresentazione in fase arteriosa, il washout, l'aspetto capsulare, il numero di lesioni, il diametro massimo delle lesioni, l'invasione macrovascolare e la diffusione extraepatica. Due lettori esperti hanno estratto in modo indipendente le caratteristiche di imaging dai referti radiologici finali. Le discrepanze sono state risolte mediante consenso, utilizzando gli stessi criteri predeterminati. Se erano disponibili sia TC che RM, è stata utilizzata la modalità con la caratterizzazione multiphasica più completa. Se entrambe le modalità erano complete, è stata utilizzata quella indicata come base per la decisione clinica finale.
Lo stadio BCLC è stato assegnato dopo che lo stato del carcinoma epatocellulare (HCC) era stato definitivamente stabilito. L'HCC in stadio precoce è stato definito come stadio BCLC 0–A. L'HCC negativo per AFP è stato definito come AFP <20 ng/mL al prelievo ematico effettuato prima dell'imaging. Queste etichette di sottogruppo sono state utilizzate per la valutazione delle prestazioni nei sottogruppi, non per l'addestramento del modello o per l'idoneità dei partecipanti.
Raccolta del campione e lavorazione del plasma
Il sangue venoso periferico è stato prelevato al mattino dopo un digiuno notturno di almeno 8 h, prima della conferma mediante TC/RMI con mezzo di contrasto e prima dell'inizio della terapia antitumorale. Sono stati raccolti complessivamente 10 mL di sangue venoso in provette per prelievo ematico con anticoagulante K2-EDTA. Le provette sono state capovolte 8–10 volte immediatamente dopo il prelievo e trasportate in laboratorio a 4 °C. L'intervallo tra il prelievo e la centrifugazione è stato mantenuto entro 2 h ed è stato registrato nel registro del campione.
Il plasma è stato separato mediante centrifugazione a 1.600 × g per 10 min a 4 °C utilizzando una centrifuga refrigerata. Il sopranatante plasmatico è stato trasferito in provette microcentrifuga da 1,5 mL prive di nucleasi e suddiviso in aliquote da 0,5–1,0 mL. Le aliquote sono state conservate a -80 °C fino al momento dell'analisi. Prima della misurazione dei biomarcatori, un'aliquote è stata scongelata una sola volta a temperatura ambiente, 20–25 °C, per 20–30 min, mescolata delicatamente per inversione ed analizzata entro 2 h dal completo scongelamento. I campioni con evidente coagulazione, emolisi da moderata a grave o marcata lipemia sono stati scartati in base ai criteri predeterminati di accettazione del laboratorio. Per il coorte di validazione prospettica, è stata effettuata una ripetuta raccolta del sangue prima della conferma mediante imaging, quando clinicamente fattibile.
Misurazione di AFP e PIVKA-II
AFP e PIVKA-II sono stati misurati mediante saggi immunoenzimatici a chemiluminescenza elettrochimica su un analizzatore immunoenzimatico automatizzato. L'AFP è stato registrato in ng/mL e il PIVKA-II in mAU/mL. I risultati sono stati esportati direttamente dal sistema informativo di laboratorio.
Ogni lotto analitico includeva calibratori specifici per il saggio e almeno due livelli interni di controllo qualità. L'accettazione del lotto richiedeva che tutti i risultati del controllo qualità rientrassero nell'intervallo di controllo di laboratorio predefinito. Gli obiettivi di precisione analitica erano un coefficiente di variazione all'interno della corsa ≤10% e un coefficiente di variazione tra le corse ≤15% sia per AFP che per PIVKA-II. Se il controllo qualità falliva o veniva rilevata una deriva del lotto, il lotto veniva respinto, l'analizzatore veniva ricalibrato e il lotto veniva ripetuto qualora rimanesse un volume sufficiente di campione. Per un risultato singolo implausibile o fuori intervallo, quando possibile venivano effettuate due misurazioni ripetute dallo stesso aliquota scongelato, e si utilizzava la media dei due valori ripetuti accettati. Se la ripetizione dei test falliva o il volume del campione era insufficiente, il risultato veniva contrassegnato come non valido e gestito secondo la regola prevista per i dati mancanti.
Estrazione di variabili cliniche di routine
Le variabili cliniche di routine sono state estratte dallo stesso incontro clinico del prelievo ematico pre-imaging per AFP/PIVKA-II. Le variabili considerate includevano età, sesso, alanina aminotransferasi, aspartato aminotransferasi, albumina, bilirubina totale, rapporto internazionale normalizzato (INR), tempo di protrombina e conta piastrinica. Quando erano disponibili più valori all'interno della finestra pre-imaging, è stato selezionato il valore più vicino al momento del campionamento di AFP/PIVKA-II. Le unità sono state standardizzate prima dell'analisi. L'età è stata registrata in anni; il sesso è stato codificato come femmina = 0 e maschio = 1; la conta piastrinica è stata registrata come ×10⁹/L; l'albumina come g/L; la bilirubina totale come μmol/L; l'AFP come ng/mL; e il PIVKA-II come mAU/mL.
Dopo l'estrazione sono stati eseguiti controlli di intervallo, controlli dell'unità di misura e controlli sui valori implausibili. I valori segnalati sono stati verificati confrontandoli con i documenti originali. Un valore è stato escluso soltanto quando il documento originale confermava un problema di qualità del campione, un errore nella conversione dell'unità di misura o una deviazione nella procedura di elaborazione. I valori estremi verificati sono stati mantenuti. Il dizionario delle variabili ha definito la struttura del dataset, i nomi delle variabili, le unità di misura, le regole di codifica, le regole di trasformazione e i codici per i valori mancanti.
Definizione del predittore e selezione delle variabili
I predittori candidati sono stati definiti prima dello sviluppo del modello. I biomarcatori principali erano AFP e PIVKA-II. I candidati clinici di routine erano età, sesso, alanina aminotransferasi, aspartato aminotransferasi, albumina, bilirubina totale, INR o tempo di protrombina e conta piastrinica. AFP e PIVKA-II sono stati trasformati mediante logaritmo naturale. I valori al di sotto del limite di segnalazione del saggio sono stati sostituiti con la metà del limite inferiore di segnalazione prima della trasformazione.
I confronti univariati tra partecipanti con HCC e senza HCC nella coorte di addestramento hanno utilizzato test t per campioni indipendenti per variabili approssimativamente distribuite in modo normale, test U di Mann-Whitney per variabili continue asimmetriche e test χ2 o test esatto di Fisher per variabili categoriali. Le variabili con P < 0,10 nella selezione univariata sono state prese in considerazione per la modellizzazione multivariabile. La multicollinearità è stata valutata mediante il fattore di inflazione della varianza (VIF), e un valore di VIF >5 è stato considerato indicativo di una collinearità preoccupante. Quando due variabili risultavano fortemente collineari, è stata mantenuta quella con un'interpretazione clinica più chiara e una migliore performance incrementale del modello.
Il modello finale di regressione logistica multivariabile è stato sviluppato soltanto nella coorte di addestramento. ln(AFP) e ln(PIVKA-II) sono stati inseriti obbligatoriamente nel modello. I predittori aggiuntivi sono stati selezionati mediante eliminazione all'indietro guidata dal criterio di informazione di Akaike. Dopo la selezione all'indietro, le variabili mantenute sono state ln(AFP), ln(PIVKA-II), età, sesso, albumina, INR e conta piastrinica.
Sviluppo del modello e calcolo del punteggio di rischio
Il modello diagnostico è stato costruito utilizzando una regressione logistica multivariata con lo stato di carcinoma epatocellulare (HCC) come risultato binario. La probabilità predetta di HCC è stata calcolata a partire dal predittore lineare η mediante la seguente equazione:
p(HCC)=1/[1+exp(-η)]
η = -4.862+0.247×ln(AFP)+0.712×ln(PIVKA-II)+0.018×età+0.331×sesso-0.062×albumina+0.554×INR-0.004×conteggio_piastrinico
L'AFP è stata misurata in ng/mL, la PIVKA-II in mAU/mL, l'età in anni, il sesso è stato codificato come femmina = 0 e maschio = 1, l'albumina è stata misurata in g/L, l'INR era senza unità di misura e il conteggio delle piastrine è stato misurato come ×109/L. Sono stati utilizzati logaritmi naturali per AFP e PIVKA-II. Un nomogramma è stato generato a partire dagli stessi coefficienti stimati. Per l'analisi della curva ROC (receiver operating characteristic), l'analisi di calibrazione, l'analisi della curva decisionale e la validazione, è stata utilizzata la probabilità predetta, e non la scala dei punti del nomogramma.
Definizioni di cutoff e gestione delle soglie
Le soglie di PIVKA-II sono state etichettate in base allo scopo analitico. Il cutoff clinico stabilito per PIVKA-II era di 40,0 mAU/mL. Il cutoff di Youden per il gruppo di addestramento per PIVKA-II era di 45,0 mAU/mL. Il cutoff a specificità fissa è stato determinato selezionando la soglia più vicina al 90% di specificità nel gruppo di addestramento; questo valore è risultato pari a 38,0 mAU/mL per PIVKA-II. Una soglia apparente utilizzata solo nella validazione di 37,5 mAU/mL è stata considerata esplorativa e non è stata utilizzata come cutoff primario per la validazione.
Per il modello combinato, la soglia binaria principale è stata ricavata dalla coorte di addestramento utilizzando l'indice di Youden ed è stata applicata invariata alla coorte di validazione. Per l'AFP è stata utilizzata la soglia clinica consolidata di 20,0 ng/mL per la definizione del sottogruppo negativo all'AFP e per il confronto con la pratica routinaria. Le soglie ottimizzate di AFP utilizzate nell'analisi ROC sono state etichettate separatamente rispetto alla soglia clinica. Le soglie sono state pertanto riportate in base allo scopo analitico e non considerate soglie intercambiabili.
Prestazioni, calibrazione e validazione del modello
La discriminazione è stata valutata mediante curve caratteristiche operative del ricevitore e area sotto la curva. Le AUC sono state riportate con intervalli di confidenza al 95%. Il confronto tra AUC ha utilizzato il test di DeLong. Le prestazioni della classificazione sono state riportate come sensibilità, specificità, valore predittivo positivo, valore predittivo negativo e accuratezza a cutoff predeterminati.
La calibrazione è stata valutata mediante intercetta di calibrazione, pendenza della calibrazione, punteggio di Brier e test di adattamento di Hosmer–Lemeshow. I grafici di calibrazione hanno confrontato le probabilità previste e osservate di HCC tra i diversi gruppi di rischio. La validazione interna ha utilizzato 1.000 ricampionamenti bootstrap nella coorte di addestramento per stimare l'AUC corretta per l'ottimismo e la pendenza della calibrazione corretta per l'ottimismo. La coorte di validazione è stata utilizzata esclusivamente per la valutazione delle prestazioni esterne. Nessun aggiornamento dei coefficienti, riparametrazione del modello, riottimizzazione delle soglie o ricalibrazione basata su sottogruppi è stato effettuato nella coorte di validazione.
L'utilità clinica è stata valutata mediante l'analisi della curva decisionale in un intervallo di probabilità soglia compreso tra 0,05 e 0,50. Il beneficio netto è stato confrontato per AFP singola, PIVKA-II singola, il modello basato sulle variabili cliniche e il modello combinato.
Gestione dei dati mancanti e analisi di sensibilità
L'analisi primaria ha utilizzato un modello basato sui casi completi per i partecipanti con tutte le variabili incluse nel modello finale. I valori mancanti sono stati riassunti per ciascun predittore candidato. Se la percentuale di dati mancanti superava il 5% per qualsiasi predittore incluso, è stata eseguita un'analisi di sensibilità mediante imputazione multipla tramite equazioni concatenate, utilizzando l'esito, tutti i predittori inclusi e l'indicatore della coorte nel modello di imputazione. Le variabili con una percentuale di dati mancanti superiore al 20% non sono state considerate idonee per la modellizzazione multivariata. I coefficienti del modello, l'AUC, la pendenza di calibrazione e le metriche di classificazione ottenute dall'analisi con dati imputati sono stati confrontati con i risultati ottenuti dai casi completi.
Software e reporting riproducibile
Tutte le analisi sono state eseguite utilizzando R versione 4.3.2 e RStudio versione 2023.12.1. La regressione logistica è stata effettuata utilizzando il pacchetto stats base. L'analisi ROC e i confronti secondo DeLong sono stati eseguiti utilizzando pROC versione 1.18.5. L'analisi di calibrazione è stata condotta utilizzando rms versione 6.7-1 e ResourceSelection versione 0.3-6. La validazione con bootstrap ha utilizzato 1.000 ricampionamenti con un seme casuale fisso pari a 70558. L'analisi della curva decisionale è stata eseguita utilizzando rmda versione 1.6. L'imputazione multipla, se necessaria, è stata effettuata utilizzando mice versione 3.16.0. Lo script di analisi ha riprodotto il flusso dei partecipanti, i coefficienti del modello, le curve ROC, le metriche di calibrazione, la validazione bootstrap, l'analisi della curva decisionale, le analisi basate su soglie e le tabelle statistiche.