$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
È importante notare che questo studio si concentra sulla descrizione delle sfumature della revisione delle varianti in IGV per la redazione di rapporti clinici. La redazione finale del rapporto clinico dovrebbe dipendere dalla specifica pipeline validata, dai criteri di qualità e dalle politiche di reporting di ciascun singolo laboratorio. Una trattazione completa di questo argomento importante va oltre gli scopi del presente studio; tuttavia, standard e linee guida raccomandati sono stati elaborati dal Collegio Americano di Genetica Medica e Genomica e dall'Associazione di Patologia Molecolare7 e vengono regolarmente applicati 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ò determinare una terapia mirata e, pertanto, è fondamentale segnalarla a livello clinico. Nonostante una chiamata positiva 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 allettica della variante riportata (VAF) è di circa l'1%, valore notevolmente inferiore a quanto ci si aspetterebbe per un evento driver in un campione stimato avere anche una bassa cellularità tumorale (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 cancro, le evidenze suggeriscono che il reperto nel caso B sia un artefatto.
Interpretazione del caso clinico n. 2:
Questo approfondimento illustra diversi errori distinti di cui i clinici devono essere consapevoli. In primo luogo, i chiamatori di varianti possono identificare in modo appropriato o inappropriato varianti genetiche adiacenti come due SNV oppure come un singolo MNV (Figura 2C). Ciò è dovuto alla differenza nell'implementazione dell'algoritmo di chiamata di varianti consapevole dell'aplotipo nei chiamatori di varianti somatiche. Solo pochi chiamatori di varianti, più recentemente pubblicati, incorporano questa caratteristica14,15. In secondo luogo, variazioni tecniche e basate sulla nomenclatura come queste vengono inserite in database esterni da centri di sequenziamento che utilizzano diverse pipeline bioinformatiche. Di conseguenza, tali database possono contenere molteplici 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 negativi sull'interpretazione clinica successiva. Nel complesso, queste possibili fonti di errore possono combinarsi e portare a un referto paziente errato. 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:
Data la traccia di copertura, l'aspetto tipico di una piccola delezione, il supporto dato dalla lettura e il noto meccanismo di attivazione costitutiva delle delezioni dell'esone 19 di EGFR nel NSCLC, si può concludere che questa 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 brevi 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 a un genoma di riferimento, producendo un file BAM in cui ogni sequenza è "impilata" nel suo locus corrispondente. (B) Le discrepanze tra il genoma di riferimento e le sequenze (ad esempio, mutazioni puntiformi, piccole inserzioni/delezioni) vengono identificate mediante 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 delle varianti KRAS c.34G>T (p.G12C) e 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 variazione G→T. Le letture forward (rosse) e reverse (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 dell'1%). La maggior parte delle letture della variante appare sul filamento forward (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 "compattata" 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 con clipping morbido. (C) Impostazione della vista "espansa" 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 di lettura, che mostra come le basi con clipping morbido 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 "compattata" 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 ingrandita 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 multicanale 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 si vedono barre 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 di tale segnale 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 è riprodotta 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 alla prima giunzione candidata e impostare la visualizzazione della traccia di allineamento su «Compatta». (B) Spostarsi lungo il locus, applicando iterativamente la scorciatoia per ordinare per base, «control-s», finché non si osserva un accumulo di letture con lo stesso tipo di soft-clipping nella parte superiore della traccia di allineamento. Maggiore è la quantità di letture con lo stesso tipo di soft-clipping (rappresentato dalla casella rossa) nell'accumulo, maggiore sarà la certezza riguardo alla presenza di 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 specifica sequenza in un database di riferimento. La selezione di una riga particolare sposterà la visualizzazione in quella posizione. (D) L'utilizzo nuovamente della funzione di ordinamento per base «control-s» mostra un accumulo di soft-clipping, questa volta sul lato opposto delle letture. Per confermare che non siano presenti 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 il confronto BLAT delle basi soft-clippate. La tabella risultante mostra il primo punto di rottura della giunzione, ovvero il locus inizialmente visualizzato, la giunzione tra l'esone 1 e l'introne 1. (E) La selezione della prima riga riporta al primo punto di rottura. Le basi presenti 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.Fare clic qui per scaricare il file.