Procesa datos genómicos con Cloud Life Sciences
En esta página, se explica cómo ejecutar una canalización de genómica que utiliza la API de Cloud Life Sciences para crear un archivo de índice (archivo BAI) a partir de un archivo binario que contiene secuencias de ADN (archivo BAM).
Por lo general, los archivos BAM son grandes y pueden tardar mucho tiempo en leerse con un visor de genomas. Usas un archivo BAI para ubicar las partes del archivo BAM que contienen la posición del genoma que te interesa.
Antes de comenzar
- Sign in to your Google Cloud account. If you're new to Google Cloud, create an account to evaluate how our products perform in real-world scenarios. New customers also get $300 in free credits to run, test, and deploy workloads.
-
In the Google Cloud console, on the project selector page, select or create a Google Cloud project.
-
Make sure that billing is enabled for your Google Cloud project.
-
Enable the Cloud Life Sciences, Compute Engine, and Cloud Storage JSON APIs.
- Install the Google Cloud CLI.
-
To initialize the gcloud CLI, run the following command:
gcloud init
-
In the Google Cloud console, on the project selector page, select or create a Google Cloud project.
-
Make sure that billing is enabled for your Google Cloud project.
-
Enable the Cloud Life Sciences, Compute Engine, and Cloud Storage JSON APIs.
- Install the Google Cloud CLI.
-
To initialize the gcloud CLI, run the following command:
gcloud init
- Instala Python 3.8.
Si usas Windows y dejaste la casilla de verificación correspondiente seleccionada cuando instalaste Google Cloud CLI, esto se hizo de forma automática.
También puedes usar Cloud Shell, que ya tiene instalada gcloud CLI.
Ejecuta la canalización
Para ejecutar la canalización, completa los siguientes pasos:
Crea un bucket en el que almacenarás el archivo BAI. Los depósitos son los contenedores básicos que conservan tus datos en Cloud Storage. Para crear un bucket llamado
PROJECT_ID-life-sciences
, ejecuta el comandogcloud storage buckets create
:gcloud storage buckets create gs://PROJECT_ID-life-sciences
Reemplaza PROJECT_ID por el ID del proyecto de Google Cloud. Debes usar un nombre de bucket único a nivel global.
Si se completa correctamente, el comando muestra lo siguiente:
Creating gs://PROJECT_ID-life-sciences
Para iniciar la canalización, ejecuta el comando
gcloud beta lifesciences pipelines run
:gcloud beta lifesciences pipelines run \ --regions us-east1 \ --command-line 'samtools index ${BAM} ${BAI}' \ --docker-image "gcr.io/cloud-lifesciences/samtools" \ --inputs BAM=gs://genomics-public-data/NA12878.chr20.sample.bam \ --outputs BAI=gs://PROJECT_ID-life-sciences/NA12878.chr20.sample.bam.bai
Si se completa correctamente, el comando muestra lo siguiente:
Running [projects/PROJECT_ID/operations/OPERATION_ID]
Toma nota de OPERATION_ID, que usarás en el siguiente paso.
Para realizar un seguimiento del estado de la canalización, ejecuta el comando
gcloud beta lifesciences operations wait
. Reemplaza OPERATION_ID con el valor impreso en el paso anterior. La canalización tarda unos minutos en finalizar.gcloud beta lifesciences operations wait OPERATION_ID
Después de que la operación finaliza, muestra el mensaje siguiente:
Waiting for [projects/PROJECT_ID/operations/OPERATION_ID]...done.
Para verificar que se haya generado el archivo BAI, ejecuta el comando
gcloud storage ls
:gcloud storage ls gs://PROJECT_ID-life-sciences
Si se completa correctamente, el comando muestra lo siguiente:
gs://PROJECT_ID-life-sciences/NA12878.chr20.sample.bam.bai
Ejecutaste una canalización con la API de Cloud Life Sciences para crear un archivo BAI a partir de un archivo BAM. Usa un visor de genomas para examinar el archivo BAM NA12878.chr20.sample.bam
con el archivo de índice NA12878.chr20.sample.bam.bai
.
Realiza una limpieza
Sigue estos pasos para evitar que se apliquen cargos a tu cuenta de Google Cloud por los recursos que usaste en esta página.
Borra el archivo BAI
Para borrar el archivo BAI generado, pero conservar el proyecto y el bucket que creaste, ejecuta el comando gcloud storage rm
:
gcloud storage rm PROJECT_ID-life-sciences/NA12878.chr20.sample.bam.bai
Borra el bucket
Si creaste el bucket específicamente para esta guía de inicio rápido y ya no lo necesitas, pero quieres conservar tu proyecto, borra el bucket con el comando gcloud storage rm
. Cuando se borra el bucket, también se borra el archivo BAI generado.
gcloud storage rm gs://PROJECT_ID-life-sciences --recursive
Borra el proyecto
Si creaste el proyecto específicamente para esta guía de inicio rápido y ya no lo necesitas, puedes borrarlo. Cuando borras el proyecto, también se borra el archivo BAI y el bucket de Cloud Storage.
- In the Google Cloud console, go to the Manage resources page.
- In the project list, select the project that you want to delete, and then click Delete.
- In the dialog, type the project ID, and then click Shut down to delete the project.
¿Cómo fue?
¿Qué sigue?
- Obtén más información sobre los conjuntos de datos públicos de la API de Cloud Life Sciences.
- Obtén información para cargar datos sobre variantes en Cloud Storage o BigQuery.
- Obtén información para analizar variantes con BigQuery.