$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
È importante notare che questo studio si concentra sulla descrizione delle sfumature della revisione delle varianti in IGV ai fini della redazione di rapporti clinici. La redazione finale del rapporto clinico dovrebbe basarsi sulla pipeline validata, sui criteri di qualità e sulla politica di reporting di ciascun laboratorio. Una trattazione completa di questo argomento importante va oltre lo scopo del presente studio; tuttavia, standard e linee guida raccomandati sono stati elaborati dal College Americano di Genetica Medica e Genomica e dall'Associazione per la Patologia Molecolare7 e vengono regolarmente implementati nei laboratori clinici.
Interpretazione del caso clinico n. 1:
Data l'elevata VAF, la distribuzione equilibrata tra i filamenti e l'incidenza nota della mutazione KRAS c.34G>T (p.G12C) nel cancro del polmone (Figura 2A), la variante nel caso A può giustificare una terapia mirata ed è quindi fondamentale segnalarla a scopo clinico. Nonostante la positività rilevata per la stessa variante KRAS c.34 G>T nel secondo caso (sarcama) (Figura 2B), sono presenti diversi segnali di allarme che indicano una chiamata falsa positiva: (1) la frazione allica della variante (VAF) riportata è di circa l'1%, valore notevolmente inferiore a quanto ci si aspetterebbe per un evento driver in un campione con una cellularità tumorale stimata anche solo moderata (20%) e (2) le letture che supportano questa variante mostrano un bias di filamento, essendo quasi esclusivamente presenti sul filamento diretto. Unitamente alla mancanza di plausibilità biologica per questo tipo di tumore, le evidenze suggeriscono che il reperto nel caso B sia un artefatto.
Interpretazione del caso clinico n. 2:
Questo approfondimento illustra diverse tipologie di errori distinti di cui i clinici devono essere consapevoli. Innanzitutto, i chiamatori di varianti possono identificare varianti genetiche adiacenti correttamente o in modo improprio come due SNV oppure come un singolo MNV (Figura 2C). Ciò è dovuto alla differenza nell'implementazione dell'algoritmo di chiamata delle varianti consapevole degli aplotipi nei chiamatori di varianti somatiche. Solo alcuni chiamatori di varianti, pubblicati più di recente, incorporano questa funzionalità14,15. In secondo luogo, variazioni tecniche e basate sulla nomenclatura come queste vengono inserite nei database esterni da centri di sequenziamento che utilizzano diverse pipeline bioinformatiche. Di conseguenza, tali database possono contenere diverse descrizioni distinte di ciò che dovrebbe essere elencato come una singola variante genetica. A meno che non venga effettuata periodicamente una fase di normalizzazione del database, ciò può avere diversi effetti indesiderati sull'interpretazione clinica successiva. Nel complesso, queste possibili fonti di errore possono combinarsi e portare a un referto errato per il paziente. Questo scenario complesso sottolinea l'importanza della visualizzazione delle varianti mediante software di alta qualità come IGV.
Interpretazione del caso clinico n. 3:
Dato il supporto multiplo alla lettura per la duplicazione e il noto meccanismo di attivazione costitutiva di KIT con duplicazioni dell'esone 9 nei GIST, possiamo concludere che questa variante KIT c.1504_1509dup (p.A502_Y503dup) è reale e non un artefatto (Figura 3). Le mutazioni di KIT sono frequenti fattori oncogeni nei GIST, verificandosi in circa il 90% dei casi. Una caratterizzazione accurata della variante è fondamentale in questo contesto, poiché, sebbene le mutazioni attivanti dell'esone 11 siano le più comuni, i GIST con mutazioni nell'esone 9 presentano un comportamento diverso rispetto alla sensibilità agli inibitori della tirosina chinasi, modificando così l'approccio terapeutico. Data l'importanza clinica delle mutazioni di KIT nei GIST, questa variante dovrebbe essere inclusa nel rapporto finale.
Interpretazione del caso clinico n. 4:
Dato il tracciato del coverage, l'aspetto tipico di una breve delezione e il supporto dei read, oltre al noto meccanismo di attivazione costitutiva delle delezioni dell'esone 19 di EGFR nel carcinoma polmonare non a piccole cellule (NSCLC), si può concludere che la variante EGFR c.2235_2249del (p.E746_A750del) è reale e non un artefatto (Figura 4). Molte mutazioni di EGFR, inclusa questa ben nota delezione dell'esone 19, sono altamente sensibili agli inibitori della tirosina chinasi. Considerata l'importanza clinica delle mutazioni di EGFR nei NSCLC, questa variante dovrebbe essere inclusa nel rapporto finale.
Interpretazione del caso clinico n. 5 e del caso clinico n. 6:
Le mappe termiche confermano l'amplificazione di EGFR e il punto di rottura tra gli esoni 1 e 8, valutato con BLAT, che conferma entrambe le alterazioni in questo paziente (Figura 5 e Figura 6). L'amplificazione di EGFR è generalmente considerata un fattore prognostico di malattia più aggressiva, e la riorganizzazione genica di EGFRvIII è un biomarcatore specifico del GBM e potenzialmente bersagliabile16,17.
I flussi di lavoro selezionati possono essere riprodotti utilizzando i file salvati nel ramo principale del repository GitHub https://github.com/Eitan177/Demo_IGV (aggiornato il 6/9/2026), che contiene sei file di sessione, uno per ogni vignetta, cinque coppie di file bam e bai, uno per ogni vignetta che utilizza un file di allineamento, e 10 file seg per la vignetta cinque, che utilizza file di segmentazione anziché file di allineamento.

Figura 1: Standard dei file NGS e flusso di lavoro del sequenziamento. (A) Le sequenze corte generate su una flow cell vengono tracciate mediante misurazioni ottiche e memorizzate in formato FASTQ, che include sia i dati di sequenza sia i punteggi di qualità per ogni base. Queste sequenze vengono quindi ordinate e allineate al genoma di riferimento, producendo un file BAM in cui ogni lettura è "impilata" nel locus corrispondente. (B) Le discrepanze tra il genoma di riferimento e le sequenze (ad esempio, mutazioni puntiformi, piccoli indel) vengono identificate da algoritmi di chiamata delle varianti. (C) Infine, le varianti rilevate vengono raccolte in un file VCF, che riassume le loro posizioni, le basi di riferimento e alternative, e altri dati di annotazione. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 2: Visualizzazione basata su IGV della variante KRAS c.34G>T (p.G12C) e della variante BRAF c.1798_1799delinsAA (p.V600K). (A) Visualizzazione di una variante reale (KRAS c.34G>T (p.G12C)) in IGV, che mostra un’elevata profondità di lettura (migliaia di letture) e circa il 35–40% di letture con una sostituzione G→T. Le letture dirette (rosse) e inverse (blu) sono bilanciate, suggerendo una mutazione eterozigote vera. (B) Artifatto sospetto in un tessuto non rilevante, che mostra una variante con supporto minimo di letture (VAF 1%). La maggior parte delle letture della variante appare sul filamento diretto (rosso), indicando un bias di filamento. (C) Una variante complessa, BRAF c.1798_1799delinsAA (p.V600K), osservata in IGV, che rivela una sostituzione di due basi (GT→AA) in cis. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3: Visualizzazione di una piccola inserzione. (A) La vista "squished" dell'esone 9 di KIT in IGV aiuta a identificare la regione di interesse. (B) Regione di interesse (riquadro rosso) identificata visualizzando la barra di inserzione e le basi soft-clippate. (C) Impostazione della vista "Expanded" in IGV per consentire la visualizzazione dell'inserzione. (D) Cliccando sulla barra di inserzione (riquadro rosso) viene visualizzata la sequenza inserita. (E) La sequenza inserita viene confrontata con la sequenza di riferimento per identificare le basi duplicate. (F) Ordinamento delle letture in base al filamento, che mostra come le basi soft-clippate corrispondano all'inserzione nelle letture inverse. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 4: Visualizzazione di una piccola delezione. (A) Vista "schiacciata" dell'esone 19 di EGFR che mostra la delezione presente in numerose letture. (B) Utilizzo della traccia di copertura per esaminare il numero di letture nelle posizioni all'interno e adiacenti alla delezione. (C) Vista "espansa" che mostra la delezione, identificando il numero di basi eliminate. (D) Il confronto con la sequenza di riferimento (riquadro rosso) identifica le basi eliminate. Cliccare qui per visualizzare una versione più grande di questa figura.

Figura 5: Visualizzazione di un guadagno del numero di copie. (A) Un file seg contiene sei colonne: ID del campione, cromosoma, inizio, fine, marcatori nel segmento e media del segmento. (B) Il caricamento diretto dei file seg in IGV consente una visualizzazione agevole. (C) Visualizzazione dei dati in una vista genomica globale. Ogni campione appare come una riga di dati, con colori che indicano i valori: rosso per >0 e blu per <0, dove colori più intensi indicano valori più distanti da 0. Dopo aver aperto il campione di interesse e nove campioni casuali dello stesso run, si ottiene una mappa termica multi-campione composta da dieci righe. Il campione di interesse è evidenziato nella parte inferiore all'interno di un riquadro rosso. (D) Restringendo la vista al locus EGFR, si osserva una barra rossa nel campione di interesse in basso, mentre nei campioni casuali le barre sono bianche o blu. Questa differenza cromatica tra il campione di interesse e i campioni casuali nelle altre righe dimostra che la chiamata di amplificazione nel campione di interesse è superiore al rumore di fondo in questo locus. (E) Allargando la vista al cromosoma 7, si nota un contrasto tra il locus contenente EGFR nel nostro campione e l'assenza dello stesso nei campioni casuali. L'amplificazione di EGFR appare come una sottile striscia rossa nel campione inferiore. (F) Restringendo la vista al locus MET, l'amplificazione di MET è replicata in misura variabile nei nove campioni casuali oltre che nel campione di interesse, suggerendo che questa chiamata di amplificazione sia un artefatto. Questo evidenzia la differenza nell'aspetto tra una chiamata reale del numero di copie al locus EGFR e un artefatto. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 6: Visualizzazione di una riorganizzazione genica di EGFRvIII. Prima di esaminare questo caso, assicurarsi che la visualizzazione delle basi soft-clippate sia stata configurata (vedere il File Supplementare 1 per ulteriori dettagli). (A) Passare al primo punto di giunzione candidato e impostare la visualizzazione della traccia di allineamento su "Squished". (B) Spostarsi lungo il locus, applicando iterativamente la scorciatoia di ordinamento per base, "control-s", finché non diventa visibile un accumulo di letture con clipping simile nella parte superiore della traccia di allineamento. Maggiore è la quantità di letture con clipping simile (riquadro rosso) nell'accumulo, maggiore sarà la certezza riguardo a una riorganizzazione strutturale. (C) Lo strumento BLAT recupera la sequenza delle basi soft-clippate e identifica dove nel genoma tale sequenza è presente, insieme a un punteggio quantitativo che misura quanto strettamente la sequenza si allinea con una sequenza specifica in un database di riferimento. La selezione di una riga particolare porterà alla posizione corrispondente. (D) L'utilizzo della funzione di ordinamento per base "control-s" mostra nuovamente un accumulo di basi soft-clippate, questa volta sul lato opposto delle letture. Per confermare l'assenza di basi aggiuntive nel campione rispetto al riferimento hg38, fare clic con il tasto destro del mouse su una lettura con una sequenza di basi soft-clippate ed eseguire nuovamente BLAT sulle basi soft-clippate. La tabella risultante mostra il primo punto di rottura della giunzione, il locus inizialmente visualizzato, ovvero la giunzione tra l'esone 1 e l'introne 1. (E) La selezione della prima riga riporta al primo punto di rottura. Le basi alla giunzione dell'allineamento di riferimento BLAT sono identiche alle basi soft-clippate utilizzate come query, indicando che l'intera sequenza di letture con soft-clipping è contenuta nell'allineamento della sequenza di riferimento. Cliccare qui per visualizzare una versione ingrandita di questa figura.
File supplementare 1: Configurazione consigliata per IGV.Cliccare qui per scaricare il file.
File supplementare 2: Flussi di lavoro condensati.Clicca qui per scaricare questo file.