Questa pagina spiega come eseguire una pipeline di analisi genomica secondaria su Google Cloud utilizzando le best practice relative al Genoma Analysis Toolkit (GATK). Le best practice GATK sono fornite dal Broad Institute.
Il flusso di lavoro utilizzato in questo tutorial è un'implementazione delle best practice GATK per il rilevamento di varianti nei dati di sequenziamento dell'intero genoma (WGS). Il flusso di lavoro è scritto nel Broad Institute's Workflow Definition Language (WDL) e viene eseguito sul runner Cromwell WDL.
Obiettivi
Dopo aver completato questo tutorial, saprai come:
- Esegui una pipeline utilizzando le best practice GATK con dati della build 38 del genoma di riferimento umano
- Esegui una pipeline utilizzando le best practice di GATK usando i tuoi dati
Costi
In questo documento, utilizzi i seguenti componenti fatturabili di Google Cloud:
- Compute Engine
- Cloud Storage
Per generare una stima dei costi in base all'utilizzo previsto, utilizza il Calcolatore prezzi.
Prima di iniziare
- Accedi al tuo account Google Cloud. Se non conosci Google Cloud, crea un account per valutare le prestazioni dei nostri prodotti in scenari reali. I nuovi clienti ricevono anche 300 $di crediti gratuiti per l'esecuzione, il test e il deployment dei carichi di lavoro.
-
In the Google Cloud console, on the project selector page, select or create a Google Cloud project.
-
Assicurati che la fatturazione sia attivata per il tuo progetto Google Cloud.
-
Abilita le API Cloud Life Sciences, Compute Engine, and Cloud Storage.
- Installa Google Cloud CLI.
-
Per initialize gcloud CLI, esegui questo comando:
gcloud init
-
Aggiorna e installa i componenti di
gcloud
:gcloud components update
gcloud components install beta -
In the Google Cloud console, on the project selector page, select or create a Google Cloud project.
-
Assicurati che la fatturazione sia attivata per il tuo progetto Google Cloud.
-
Abilita le API Cloud Life Sciences, Compute Engine, and Cloud Storage.
- Installa Google Cloud CLI.
-
Per initialize gcloud CLI, esegui questo comando:
gcloud init
-
Aggiorna e installa i componenti di
gcloud
:gcloud components update
gcloud components install beta - Installa git per scaricare i file richiesti.
-
Per impostazione predefinita, Compute Engine applica quote delle risorse per prevenire utilizzi involontari. Aumentando le quote, puoi avviare
più macchine virtuali contemporaneamente, aumentando la velocità effettiva e riducendo
i tempi di risposta.
Per ottenere risultati ottimali in questo tutorial, devi richiedere una quota aggiuntiva superiore a quella predefinita del progetto. I suggerimenti per gli aumenti della quota sono forniti nel seguente elenco e le quote minime necessarie per eseguire il tutorial. Effettua le tue richieste di quota nella regione
us-central1
:- CPU: 101 (minimo 17)
- Persistent Disk Standard (GB): 10.500 (minimo 320)
- Indirizzi IP in uso: 51 (minimo 2)
Puoi lasciare vuoti gli altri campi di richiesta della quota per mantenere le quote attuali.
crea un bucket Cloud Storage
Crea un bucket Cloud Storage utilizzando il comando gsutil mb
. A causa di un requisito nel motore Cromwell, non utilizzare un trattino basso (_
) nel nome del bucket, altrimenti riscontrerai un errore.
gsutil mb gs://BUCKET
La pipeline restituisce risultati, log e file intermedi in questo bucket.
Scarica i file di esempio
Esegui questi comandi per scaricare lo script WDL e helper:
git clone https://github.com/broadinstitute/wdl-runner.git git clone https://github.com/gatk-workflows/broad-prod-wgs-germline-snps-indels.git
Il repository gatk-workflows/broad-prod-wgs-germline-snps-indels contiene i seguenti file necessari per eseguire la pipeline:
*.wdl
: definizione del flusso di lavoro*.inputs.json
: parametri di input, inclusi i percorsi dei file BAM e il genoma di riferimento*.options.json
: opzioni di runtime del flusso di lavoro
Puoi trovare il file di definizione della pipeline Cromwell utilizzato per eseguire le pipeline WDL nel repository broadinstitute/wdl-runner/wdl_runner/.
Esegui la pipeline utilizzando dati di esempio
Questa sezione mostra come eseguire la pipeline con dati WGS utilizzando la build 38 del genoma di riferimento umano. I file di input sono file BAM non allineati.
Per eseguire la pipeline, completa i seguenti passaggi:
Crea la variabile di ambiente
GATK_GOOGLE_DIR
che punta alla cartella contenente i file della pipeline Broad:export GATK_GOOGLE_DIR="${PWD}"/broad-prod-wgs-germline-snps-indels
Crea la variabile di ambiente
GATK_OUTPUT_DIR
che punta al bucket Cloud Storage e a una cartella per ioutput
del flusso di lavoro, i filework
intermedi elogging
:export GATK_OUTPUT_DIR=gs://BUCKET/FOLDER
Passa alla directory
/wdl_runner
nel repository che hai scaricato. Questa directory contiene il file di definizione per l'esecuzione di pipeline basate su WDL su Google Cloud:cd wdl-runner/wdl_runner/
Esegui la pipeline:
Scegli una delle seguenti opzioni a seconda che utilizzi un VPC predefinito o un VPC personalizzato:
VPC predefinito
gcloud beta lifesciences pipelines run \ --pipeline-file wdl_pipeline.yaml \ --location us-central1 \ --regions us-central1 \ --inputs-from-file WDL=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.wdl,\ WORKFLOW_INPUTS=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.hg38.inputs.json,\ WORKFLOW_OPTIONS=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.options.json \ --env-vars WORKSPACE=${GATK_OUTPUT_DIR}/work,\ OUTPUTS=${GATK_OUTPUT_DIR}/output \ --logging ${GATK_OUTPUT_DIR}/logging/
VPC personalizzato
Crea le variabili di ambiente
NETWORK
eSUBNETWORK
per specificare il nome della rete e della subnet VPC:export NETWORK=VPC_NETWORK export SUBNETWORK=VPC_SUBNET
Modifica il file
PairedEndSingleSampleWf.options.json
che si trova nella directorybroad-prod-wgs-germline-snps-indels
e modifica le zone in modo da includere solo quelle all'interno della regione della tua subnet. Ad esempio, se utilizzi una subnetus-central1
, il campozones
sarà simile a questo:"zones": "us-central1-a us-central1-b us-central1-c us-central1-f"
.gcloud beta lifesciences pipelines run \ --pipeline-file wdl_pipeline.yaml \ --location us-central1 \ --regions us-central1 \ --network ${NETWORK} \ --subnetwork ${SUBNETWORK} \ --inputs-from-file WDL=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.wdl,\ WORKFLOW_INPUTS=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.hg38.inputs.json,\ WORKFLOW_OPTIONS=${GATK_GOOGLE_DIR}/PairedEndSingleSampleWf.options.json \ --env-vars WORKSPACE=${GATK_OUTPUT_DIR}/work,\ OUTPUTS=${GATK_OUTPUT_DIR}/output,\ NETWORK=${NETWORK},\ SUBNETWORK=${SUBNETWORK} \ --logging ${GATK_OUTPUT_DIR}/logging/
Il comando restituisce un ID operazione nel formato
Running [operations/OPERATION_ID]
. Puoi utilizzare il comandogcloud beta lifesciences describe
per monitorare lo stato della pipeline eseguendo questo comando (assicurati che il valore del flag--location
corrisponda alla posizione specificata nel passaggio precedente):gcloud beta lifesciences operations describe OPERATION_ID \ --location=us-central1 \ --format='yaml(done, error, metadata.events)'
Il comando
operations describe
restituiscedone: true
al termine della pipeline.Puoi eseguire uno script incluso con
wdl_runner
per verificare ogni 300 secondi se il job è in esecuzione, è stato completato o ha restituito un errore:../monitoring_tools/monitor_wdl_pipeline.sh OPERATION_ID us-central1 300
Al termine della pipeline, esegui questo comando per elencare gli output nel bucket Cloud Storage:
gsutil ls gs://BUCKET/FOLDER/output/
Puoi visualizzare i file intermedi creati dalla pipeline e scegliere quelli da conservare oppure rimuoverli per ridurre i costi associati a Cloud Storage. Per rimuovere i file, consulta l'articolo sull'eliminazione dei file intermedi nel bucket Cloud Storage.
Esegui la pipeline di best practice GATK sui tuoi dati
Prima di eseguire la pipeline sui dati locali, devi copiarli in un bucket Cloud Storage.
Copia file di input
La pipeline può essere eseguita con file BAM non allineati archiviati in Cloud Storage. Se i file sono in un formato diverso, ad esempio BAM o FASTQ allineato, devi convertirli prima di poter essere caricati in Cloud Storage. Puoi convertirli localmente o utilizzare l'API Pipelines per convertirli nel cloud.
L'esempio seguente mostra come copiare un singolo file da un file system locale a un bucket Cloud Storage:
gsutil -m -o 'GSUtil:parallel_composite_upload_threshold=150M' cp FILE \ gs://BUCKET/FOLDER
Per altri esempi su come copiare file in un bucket Cloud Storage, consulta la sezione Copia di dati in Cloud Storage.
Lo strumento a riga di comando gsutil
verifica automaticamente i checksum, in modo che,
una volta avvenuto il trasferimento, i tuoi dati siano compatibili per l'utilizzo con le best
practice GATK.
Esegui la pipeline sui tuoi dati
Per eseguire le best practice GATK sui tuoi file BAM non allineati, crea una copia di PairedEndSingleSampleWf.hg38.inputs.json
, quindi aggiorna i percorsi in modo che puntino ai tuoi file in un bucket Cloud Storage.
Puoi quindi seguire i passaggi descritti in Eseguire la pipeline utilizzando dati di esempio,
utilizzando il file PairedEndSingleSampleWf.hg38.inputs.json
aggiornato.
Se i dati non sono costituiti da file BAM non allineati e contengono genomi di riferimento, sequenziamento degli esoma, riquadri mirati e dati somatici, devi utilizzare flussi di lavoro diversi. Per ulteriori informazioni, consulta il GATK Support Forum e il repository GitHub di Broad Institute.
Risoluzione dei problemi
La pipeline è configurata per utilizzare le istanze di Compute Engine in regioni e zone specifiche. Quando esegui l'interfaccia a riga di comando gcloud, vengono utilizzate automaticamente una regione e una zona predefinite in base alla località in cui è stato creato il progetto Google Cloud. Questo può comportare il seguente messaggio di errore durante l'esecuzione della pipeline:
"ERROR: (gcloud.beta.lifesciences.pipelines.run) INVALID_ARGUMENT: Error: validating pipeline: zones and regions cannot be specified together"
Per risolvere il problema, rimuovi la regione e la zona predefinite eseguendo questi comandi, quindi esegui di nuovo la pipeline:
gcloud config unset compute/zone gcloud config unset compute/region
Per ulteriori informazioni sull'impostazione della regione e della zona predefinite per il tuo progetto Google Cloud, consulta Modifica della zona o della regione predefinite.
Se riscontri problemi durante l'esecuzione della pipeline, consulta Risoluzione dei problemi dell'API Cloud Life Sciences.
GATK ha aspettative rigide in merito ai formati dei file di input. Per evitare problemi, puoi verificare che i file superino ValidateSamFile.
Se l'esecuzione GATK non va a buon fine, puoi controllare i log eseguendo questo comando:
gsutil ls gs://BUCKET/FOLDER/logging
Se si verificano errori di autorizzazione, verifica che il tuo account di servizio disponga dell'accesso in lettura ai file di input e in scrittura al percorso del bucket di output. Se stai scrivendo file di output in un bucket di un progetto Google Cloud non di tua proprietà, devi concedere all'account di servizio l'autorizzazione ad accedere al bucket.
esegui la pulizia
Eliminazione dei file intermedi nel bucket Cloud Storage
Quando esegui la pipeline, questa archivia i file intermedi in
gs://BUCKET/FOLDER/work
. Puoi rimuovere i file al termine del flusso di lavoro per ridurre gli addebiti di Cloud Storage.
Per visualizzare la quantità di spazio utilizzata nella directory work
, esegui questo comando. L'esecuzione del comando potrebbe richiedere diversi minuti a causa delle dimensioni dei file nella directory.
gsutil du -sh gs://BUCKET/FOLDER/work
Per rimuovere i file intermedi nella directory work
, esegui questo comando:
gsutil -m rm gs://BUCKET/FOLDER/work/**
Elimina il progetto
Il modo più semplice per eliminare la fatturazione è quello di eliminare il progetto utilizzato per il tutorial.
Per eliminare il progetto:
- Nella console Google Cloud, vai alla pagina Progetti.
-
Nell'elenco dei progetti, selezionare quello da eliminare e fai clic su Elimina progetto.
- Nella finestra di dialogo, digita l'ID del progetto e fai clic su Chiudi per eliminare il progetto.
Passaggi successivi
- Questo tutorial mostra come eseguire un flusso di lavoro predefinito in un caso d'uso limitato, ma non è destinato all'esecuzione in produzione. Per informazioni su come eseguire l'elaborazione dei dati genomici in un ambiente di produzione su Google Cloud, consulta Architettura di riferimento per l'elaborazione dei dati genomici.
- Il sito di GATK e i forum del Broad Institute forniscono ulteriori informazioni di base, documentazione e assistenza per gli strumenti GATK e WDL.