Gli approcci di screening aggregato basati su CRISPR-Cas9 sono ben consolidati in diversi contesti cellulari. Tuttavia, lo screening nelle cellule chemoresistenti presenta sfide distinte e richiede ulteriori ottimizzazioni. Poiché le cellule chemoresistenti rappresentano modelli preziosi per studiare le vulnerabilità dei tumori chemoresistenti, è importante affrontare queste limitazioni nelle strategie di screening per identificare in modo affidabile i regolatori genetici o epigenetici della chemioresistenza.
La generazione di linee cellulari chemoresistenti può essere macchinosa e dipende fortemente dall'agente chemioterapico utilizzato. Sebbene le cellule possano adattarsi facilmente a certi farmaci, possono essere altamente suscettibili adaltri 25. Pertanto, è importante determinare i valori IC50 specifici per linee cellulari prima del protocollo di generazione delle linee cellulari resistenti. Nei casi in cui le cellule non tollerano l'esposizione al farmaco a livello IC50 su più passaggi, è possibile partire da una concentrazione di farmaco più bassa (ad esempio, IC10-20) (vedi Protocollo sezione 1)17. Poiché l'acquisizione della resistenza è un evento stocastico, diversi pozzi/piastre di cellule possono essere trattati contemporaneamente per aumentare la probabilità di ottenere cellule resistenti. È fondamentale caratterizzare in modo completo la popolazione cellulare ottenuta ai farmaci resistenti per descrivere i cambiamenti trascrizionali e fenotipici durante la procedura di aumento della dose rispetto al loro stato ingenuo. Approcci di profilazione globale come RNA-seq possono essere impiegati per identificare geni espressi differenzialmente rispetto alla linea parentale. Dal punto di vista funzionale, la resistenza acquisita ai farmaci dovrebbe essere validata con saggi di vitalità e proliferazione cellulare per dimostrare che le cellule resistenti mantengono la sopravvivenza in condizioni di trattamento farmacologico citotossiche per le cellule genitori. Inoltre, i potenziali cambiamenti nella loro capacità proliferativa e nel tasso di crescita dovrebbero essere misurati rispetto alla linea cellulare parentale per calcolare i rispettivi PDL. Poiché la resistenza acquisita al Taxol è nota per essere mediata, almeno in parte, da effluxi farmacologici mediati da ABC, l'espressione dei principali trasportatori ABC (ad esempio, ABCB1) dovrebbe essere valutata ai livelli di mRNA e proteina tramite qPCR e Western blotting, e può essere integrata da dati trascrictomici. Sulla base del profilo di espressione del trasportatore ABC, dovrebbe essere considerata la presenza del fenotipo MDR e la resistenza ad altri agenti chemioterapeutici dovrebbe essere valutata. Nel nostro studio precedente, abbiamo caratterizzato le linee cellulari resistenti al Taxol e identificato ABCB1 come uno dei geni più alzati al massimo, e abbiamo ritenuto che la resistenza multifarmacologica mediata da ABCB1 può causare fallimento nella selezionedella puromicina 17. Pertanto, dopo la generazione di cellule resistenti, l'espressione dei trasportatori ABC e del fenotipo MDR devono essere valutati prima di qualsiasi procedura di trasfezione.
Il metodo di aumento della dose offre vantaggi rispetto ai protocolli di esposizione ad alte dosi ai farmaci, poiché le cellule possono essere criopreservate ad ogni fase prima di aumentare la dose, offrendo flessibilità e tutelando contro la perdita di coltura. Inoltre, l'aumento graduale della dose permette l'acquisizione di fenotipi di resistenza nel tempo, consentendo l'evoluzione adattiva; a differenza del trattamento acuto ad alta dose, che comporta la selezione di sottopopolazioni resistenti preesistenti. La sua progressione graduale permette ai ricercatori di scegliere deliberatamente il livello di resistenza desiderato al quale fermare la selezione e valutare le proprietà del fenotipo risultante. Inoltre, i passaggi intermedi dell'aumento della dose potrebbero essere sfruttati per studiare i meccanismi di resistenza acquisita in modo temporale. Secondo la questione biologica sottostante, si potrebbero ottenere diversi livelli di resistenza per studiare sia i meccanismi adattivi precoci con livelli di resistenza più bassi sia fenotipi di resistenza stabili e completamente stabiliti in cellule più resistenti. Inoltre, le soglie di resistenza clinicamente rilevanti variano sostanzialmente a seconda del farmaco, della linea cellulare e del tasso di proliferazione intrinseca; pertanto, la tempistica di aumento della dose può variare di conseguenza (nel nostro caso, la generazione di cellule resistenti al Taxol richiedeva circa 6 mesi). In linea con questa logica, nel nostro studio precedente abbiamo stabilito due diversi livelli di resistenza con SUM159PT cellule e confrontato le loro caratteristiche distinte econdivise 17. Pertanto, raccomandiamo ai ricercatori di definire l'intervallo di IC50 target del farmaco di interesse in base alla rilevanza clinica e biologica.
È ben stabilito che l'astinenza della pressione del farmaco può portare alla perdita della chemioresistenza, poiché le cellule possono riprogrammare il loro epigenoma e iltrascrisimo 25. Per questo motivo, si suggerisce di valutare la stabilità del fenotipo di resistenza per diversi mesi. Se le cellule tendono a tornare a uno stato sensibile dopo la sospensione del farmaco, le cellule resistenti possono essere coltivate in presenza della dose di mantenimento, tipicamente corrispondente alla concentrazione finale del farmaco utilizzata durante la generazione della resistenza, per preservare il loro fenotipo di resistenza. Sebbene i modelli di linee cellulari chemoresistenti in vitro siano pratici, riproducibili e ampiamente utilizzati per studiare la chemioresistenza acquisita, potrebbero non ricapitolare completamente l'eterogeneità fenotipica osservata nei tumori dei pazienti e non avere il microambiente tumorale, le interazioni immunitarie e l'eterogeneità spaziale 3,4,25. Inoltre, il passaggio prolungato durante la generazione della resistenza può causare alterazioni genetiche ed epigenetiche che si discostano da quelle osservate nei campioni di pazienti, limitando potenzialmente la rilevanza traslazionale degli abbinamenti identificati. Pertanto, i geni e i meccanismi candidati identificati tramite questo protocollo dovrebbero essere ulteriormente validati in sistemi sperimentali più rilevanti dal punto di vista clinico, come modelli di co-coltura, organoidi derivati dal paziente o studi in vivo. EPIKOL offre vantaggi grazie alla sua natura focalizzata se si vogliono eseguire schermi complementari su questi modelli, dove il numero massimo di celle ottenute è intrinsecamente limitato.
Le cellule chemoresistenti spesso crescono più lentamente rispetto alle loro controparti genitori, il che può rappresentare una sfida pratica per gli screening genomici che richiedono grandi popolazioni cellulari per mantenere la copertura delle biblioteche. Librerie focalizzate come EPIKOL offrono il vantaggio di effettuare uno screening con meno materiale di partenza, fornendo al contempo una maggiore profondità dello sgRNA per gene, migliorando così sensibilità e riproducibilità all'interno dello spazio target definito. Tuttavia, per design, EPIKOL cattura solo i meccanismi di resistenza epigenetica e non copre i contributori non epigenetici alla chemioresistenza. La scelta della biblioteca dovrebbe quindi essere guidata da esperimenti precedenti e dalla specifica questione biologica. Gli utenti di queste librerie dovrebbero essere consapevoli dei limiti delle biblioteche focalizzate e interpretare i risultati all'interno dell'ambito della libreria.
Un ostacolo importante nello screening CRISPR dei modelli chemoresistenti è l'emergere del fenotipo MDR guidato dalla regolazione all'alza del trasportatoreABC 9,10. Diversi trasportatori ABC mediano l'efflusso di diversi substrati. Una volta generata una linea cellulare resistente ai farmaci, è importante valutare se la resistenza si estende ad altri agenti chemioterapeutici. La resistenza alla puromicina è causata principalmente dalla sovraespressione di ABCB1, e questo influisce sul processo di selezione degli approcci di screening CRISPR aggregati, poiché la maggior parte di queste librerie contiene cassette di selezione della puromicina per l'arricchimento delle cellule trasdotte. Qui abbiamo utilizzato verapamil, un inibitore di ABCB1 di prima generazione, per facilitare l'accumulo intracellulare di puromicina ed eliminare selettivamente le cellule resistenti non trasdotte che altrimenti sfuggerebbero alla selezione a causa dell'efflusso farmacologico mediato da ABCB1. Poiché il verapamilo si lega ad ABCB1 in modo reversibile, la funzione di ABCB1 viene ripristinata entro poche ore dopo la rimozione del verapamil. Pertanto, il verapamil non conferisce una vulnerabilità aggiuntiva alle cellule resistenti durante la fase successiva di trattamento farmacologico dello screening. Tuttavia, va riconosciuto diversi limiti del verapamil. Poiché il verapamilo è un inibitore dei canali calcici e di altri trasportatori, può influenzare la fisiologia cellulare oltre l'inibizione diABCB1 10. Inoltre, la concentrazione efficace di verapamillo varia tra le linee cellulari e i modelli MDR a seconda del livello di espressione di ABCB1 e dell'attività di altri trasportatori. Pertanto, la concentrazione di verapamillo dovrebbe essere ottimizzata per ciascuna linea cellulare, e potrebbe essere presa in considerazione una nuova generazione di inibitori di ABCB1 più specifici per ABCB1 se fossero preoccupati effetti fuori bersaglio. Come ulteriore passo fondamentale, la titolazione virale dovrebbe essere effettuata sulla cellula resistente stessa piuttosto che sulle controparti parentali, poiché l'infettività e la selezionabilità spesso differiscono per le cellule resistenti.
Un'altra considerazione importante è determinare la durata appropriata dello screening, poiché le cellule chemoresistenti spesso crescono più lentamente rispetto ai loro omologhi genitoriabbinati 26. Per tenere conto di questa differenza, abbiamo utilizzato calcoli di livello di raddoppio della popolazione (sezione 6 del protocollo) per valutare quando i gruppi di controllo e di trattamento hanno raggiunto la soglia desiderata di PDL. Gli screening possono quindi essere interrotti in giorni diversi, permettendo raddoppi della popolazione comparabili e consentendo un'interpretazione più sicura degli effetti degli sgRNA nel gruppo di trattamento.
Le PCR per la preparazione della libreria sono considerevolmente importanti per l'accuratezza dell'identificazione dei colpi a valle. Il bias introdotto da bassi input di gDNA e cicli PCR aumentati può portare all'amplificazione di alcune sequenze di sgRNA, causando arricchimenti o perdite di falsi positivi. Per mitigare questo, raccomandiamo di utilizzare la massima quantità possibile di gDNA, rispettare i numeri di ciclo PCR offerti da questo protocollo e eseguire reazioni parallele per il pool prima del sequenziamento. Quando possibile, la copertura bibliotecaria non dovrebbe scendere al di sotto delle soglie raccomandate (>300× copertura) per garantire una rilevazione affidabile degli sgRNA esauriti piuttosto che una caduta stocastica degli sgRNA a bassa abbondanza dovuta a una copertura bassa.
Dopo il sequenziamento, abbiamo utilizzato il framework di analisi MAGeCK e eseguito diverse analisi di controllo qualità a valle e passaggi di interpretazionedei dati 24. Innanzitutto, l'analisi delle componenti principali (PCA) è stata utilizzata per valutare la distribuzione e la riproducibilità delle repliche biologiche (dati non mostrati, output MAGeCK). Si prevede un clustering replicato in PCA all'interno di ciascun braccio; Tuttavia, una maggiore variabilità nel braccio trattato dal farmaco è una caratteristica comune degli screening per la chemioresistenza, poiché il farmaco può indurre cambiamenti stocastici. Se i controlli positivo e negativo si comportano come previsto in ciascun braccio, questa variabilità non impedisce l'identificazione dei veri positivi. Successivamente, è stata applicata la normalizzazione dei conteggi di sgRNA in letture per milione (RPM) nei momenti iniziali e finali per valutare se gli sgRNA che mirano ai geni essenziali fossero esauriti comeprevisto 14. Successivamente, la normalizzazione mediana è stata utilizzata per aggregare gli effetti a livello di sgRNA, permettendo l'identificazione di esaurimento a livellogenico 24. I punteggi di esaurimento a livello genico sono stati visualizzati utilizzando grafici a cascata, che classificano i geni dai più impoveriti a quelli arricchiti. Le alterazioni del numero di copie e la poliploidia derivanti da un'esposizione prolungata a farmaci possono confondere i risultati degli screening a causa di un taglio Cas9 inefficienteo eccessivo 27. Sebbene EPIKOL prenda di mira i regolatori della cromatina che generalmente non sono soggetti ad amplificazioni basate su Taxol, raccomandiamo di valutare lo stato della ploidia e i cambiamenti nel numero di copie nelle cellule resistenti prima dello screening. Se si sospettano cambiamenti nel numero di copia, durante l'analisi MAGeCK possono essere applicati strumenti di correzione computazionale come CRISPRcleanR o CERES.
Negli screening che coinvolgono due bracci sperimentali, gruppi di controllo e di trattamento, sono necessari confronti diversi per interpretare correttamente le prestazioni degli screening e le dipendenzebiologiche 17. I confronti iniziali e finali sono adatti a valutare le prestazioni dello screening, inclusi il comportamento atteso di geni essenziali e controlli non targeting. Il confronto tra l'endpoint del gruppo di controllo (DMSO) e il punto temporale iniziale consente di identificare i geni necessari per la crescita e la fitness delle cellule resistenti in assenza di pressione del farmaco. Al contrario, il confronto di due campioni endpoint, trattato con farmaco (Taxol) e controllo (DMSO), rivela i geni specifici necessari per la resistenza al farmaco piuttosto che per la fitness generale. Va notato che il confronto diretto dei due campioni finali porta a una distribuzione distorsionata dei controlli genici essenziali, poiché sono necessari per la sopravvivenza in entrambi i gruppi. Un metodo alternativo di visualizzazione potrebbe essere il grafico Volcano, che mostra le variazioni di piega log2 insieme ai valori di significatività, permettendo l'identificazione di geni significativamente impoveriti. Una volta identificati i geni candidati, i conteggi normalizzati degli sgRNA e i grafici di visualizzazione sono utili per decidere quali sono gli sgRNA con le migliori prestazioni per i saggi di validazione a valle. Per lo screening EPIKOL, raccomandiamo di dare priorità ai geni candidati supportati dal comportamento concordante di almeno 6 su 10 sgRNA, poiché ciò fornisce maggiore fiducia negli effetti a livello genico.
Qui forniamo un protocollo ottimizzato per lo screening EPIKOL su linee cellulari resistenti al Taxol che hanno aumentato la regolazione di ABCB1 e hanno mostrato il fenotipo MDR. Questo protocollo è ampiamente applicabile a qualsiasi modello chemoresistente che mostri resistenza incrociata alla puromicina o ad altri antibiotici di selezione, consentendo così lo screening CRISPR in queste cellule. Indicazioni pratiche per le sfide tecniche più comuni incontrate in tutte le fasi del protocollo sono fornite nella Tabella Supplementare 1.