Se requiere una suscripción a JoVE para ver este contenido. Inicie sesión o comience su prueba gratuita.

Artículo de método

Descubrimiento de genes conductores en tumores timos similares al cáncer de 29 años de origen colorrectal

6.2K visualizaciones

⸱

DOI:

10.3791/61077

⸱

22 de julio de 2020

 ,  ,  , 

En este artículo

Resumen

Aquí se presenta un protocolo para descubrir los genes conductores sobreexpresados que mantienen las células similares al cáncer derivadas de las células HT29 colorrectales. RNAseq con bioinformática disponible se realizó para investigar y examinar las redes de expresión génica para esclarecer un mecanismo potencial implicado en la supervivencia de las células tumorales dirigidas.

Resumen

Las células madre cancerosas desempeñan un papel vital contra las terapias clínicas, contribuyendo a la recaída tumoral. Hay muchos oncogenes involucrados en la tumorigenesis y el inicio de las propiedades del tallo del cáncer. Dado que la expresión génica en la formación de tumores tirectales derivados del cáncer no está clara, se necesita tiempo para descubrir los mecanismos que trabajan en un gen a la vez. Este estudio demuestra un método para descubrir rápidamente los genes conductores implicados en la supervivencia de las células similares al cáncer colorrectal in vitro. En este estudio se seleccionaron y utilizaron células cancerosas HT29 colorrectales que expresan el LGR5 cuando se cultivan como esferoides y acompañan a un aumento de los marcadores de vástago CD133. El protocolo presentado se utiliza para realizar RNAseq con bioinformática disponible para descubrir rápidamente los genes controladores sobreexpresados en la formación de tumores tirásticos derivados del tallo colorrectal. La metodología puede detectar y descubrir rápidamente posibles genes conductores en otros modelos de enfermedades.

Introducción

El cáncer colorrectal (CRC) es una de las principales causas de muerte con alta prevalencia y mortalidad en todo el mundo1,,2. Debido a mutaciones genéticas y amplificaciones, las células cancerosas crecen sin control proliferativo, lo que contribuye a la supervivencia celular3,la antipoptosis4y la íltividad del cáncer5,,6,,7. Dentro de un tejido tumoral, la heterogeneidad tumoral permite que las células tumorales se adapten y sobrevivan durante los tratamientos terapéuticos8. Las células madre cancerosas (CSC), con una mayor tasa de auto-renovación y pluripotencia que los tipos diferenciales de cáncer, son principalmente responsables de la recurrencia tumoral9,,10 y la CRC metastásica11. Los CCC presentan más resistencia a los medicamentos12,13,14 y propiedades anti-apoptosis15,16, sobreviviendo así a las quimioterapias tumorales.

Aquí, con el fin de investigar el mecanismo potencial para el tallo en las células madre CRC seleccionadas, RNAseq se realizó para examinar los genes expresados diferencialmente en los esferoides tumorales. Las células cancerosas pueden formar esferoides (también llamados tumoresferas) cuando se cultivan en condiciones de baja adherencia y estimuladas por factores de crecimiento añadidos al medio cultivado, incluyendo EGF, bFGF, HGF, y IL6. Por lo tanto, seleccionamos células tumorales CRC HT29 que resisten las quimioterapias con un aumento de STAT3 fosforilado cuando se tratan con oxaliplatina e irinotecon17. Además, HT29 expresó marcadores de tallo más altos cuando se cultiva en las condiciones de cultivo descritas. El modelo CSC derivado de HT29 expresaba mayores cantidades de receptor acoplado por proteína G rico en leucina 5 (LGR5)18, un marcador específico de células madre CRC19,20. Además, CD133, considerado un biomarcador general para las células madre cancerosas, también está muy expresado en la línea celular HT2921. El propósito de este protocolo es descubrir grupos de genes conductores en las tumoresferas similares al cáncer establecidas basadas en conjuntos de datos bioinformáticos en lugar de investigar oncogenesindividuales 22. Investiga posibles mecanismos moleculares a través del análisis RNAseq seguido de análisis bioinformáticos disponibles.

La secuenciación de próxima generación es un método de secuenciación de ADN de alto rendimiento, fácil de disponible y fiable basado en la ayuda computacional, utilizado para examinar exhaustivamente los genes del conductor para guiar las terapias tumorales23. La tecnología también se utiliza para detectar la expresión génica a partir de la transcripción inversa de una muestra de ARN aislada24. Sin embargo, al realizar el cribado con RNAseq, los genes más importantes a atacar con terapia pueden no tener el diferencial de expresión más alto entre las muestras experimentales y de control. Por lo tanto, se desarrollaron algunas bioinformáticas para clasificar e identificar genes basados en conjuntos de datos actuales como KEGG25,GO26,,27o PANTHER28,incluyendo Ingenuity Pathway Analysis (IPA)29 y NetworkAnalyst30. Este protocolo muestra la integración de RNAseq y NetworkAnalyst para descubrir rápidamente un grupo de genes en los esferoides derivados de HT29 seleccionados en comparación con las células HT29 parentales. También se sugiere la aplicación de este método a otros modelos de enfermedades para descubrir diferencias en genes importantes.

En comparación con la investigación de la expresión génica individual, una técnica de alto rendimiento proporciona ventajas para encontrar genes potenciales del conductor fácilmente para la medicina de precisión tumoral. Con conjuntos de datos útiles como KEGG, GO o PANTHER, se pueden identificar genes específicos basados en los modelos de la enfermedad, las vías de señalización o funciones específicas, y esto permite centrarse rápidamente en genes específicos e importantes, ahorrando tiempo y costos de investigación. Una aplicación similar se utiliza en estudios anteriores14,18,31. En particular, un tumor es más complicado porque diferentes tipos de tumores expresan genes y vías distintivas para la supervivencia y la proliferación. Por lo tanto, este protocolo puede detectar genes que distinguen diferentes tipos de tumores en diferentes circunstancias. Existe el potencial de encontrar estrategias eficaces contra el cáncer mediante la comprensión del mecanismo de expresión génica específica.

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Protocolo

1. Cultivo celular y formación de la tumoral

  1. Cultivo de células HT29 en un plato de 10 cm que contiene el medio de águila modificada de Dulbecco (DMEM) con 10% de suero bovino fetal (FBS) y 1% de antibiótico de penicilina-estreptomicina (P/S).
  2. Cultivar las células en una incubadora a 37oC con 5% de CO2 y 95% de humedad en condiciones asépticas, hasta que alcancen el 80% de confluencia.
  3. Trypsinize HT29 cells with 1 mL of 0.25% trypsin for 5 min at 37 oC and asly neutrally the trypsin by adding 2 mL of DMEM with 10% FBS and 1% P/S.
  4. Cuente las células HT29 usando un hemocicómetro.
  5. Añadir 2.000 células/pozo a una baja unión de 6 placas de pozo con 2 ml de DMEM sin suero con 1% P/S y complementado con 0.2% B27, 20 ng/ml del factor de crecimiento epidérmico (EGF), 20 ng/ml del factor de crecimiento de fibroblastos (bFGF), 20 ng/ml de factor de crecimiento de hepatocitos (HGF) y 20 ng/ml de interleucina 6 (IL6).
  6. Cultivar las células a 37oC con un 5% de CO2 y un 95% de humedad en condiciones asépticas.
  7. Añadir 0,5 ml de medio de células madre cancerosas cada 2 días hasta que las tumoresferas miden >100 m de diámetro durante al menos 7 días.
  8. Observar y medir el diámetro de la tumoressfera utilizando un microscopio invertido con un sistema de imágenes celulares digitales.
  9. Cuando las tumoresferas alcancen >100 m de diámetro, tripsan las tumoresferas con 0,25% de trippsina durante 5 min a 37oC y neutralicen la trippsina añadiendo 2 veces el volumen del medio de crecimiento. Cuente las células usando un hemocicómetro.
  10. Centrifugar a 1.200 rpm durante 10 min y retire el sobrenadante.
  11. Incubar 5 x 104 células con 2 l de anti-LGR5-PE y 2 l de anti-CD133-PE en 100 s de DMEM individualmente durante 30 min a temperatura ambiente, agitando a 200 rpm.
    NOTA: CD133 es un biomarcador general de células madre cancerosas que se expresa altamente en células HT29.
  12. Agregue 900 l de PBS y analice la expresión LGR5 y CD133 utilizando la citometría de flujo. Un cambio en la fluorescencia en el canal FL2-H indica la expresión génica.

2. Aislamiento de ARN

NOTA: Utilice un kit comercial (consulte La Tabla de materiales)con una columna rápida para el aislamiento del ARN siguiendo las instrucciones del fabricante.

  1. Añadir 50 l de PBS a las células cosechadas (2 x 105 células) y resuspenderlas con pipeteo.
  2. Añadir 200 l de tampón de liceo que contenga 2 l de beta-mercaptoetanol (-ME). Vórtice rápidamente y dejar reposar durante 5 minutos.
  3. Centrifugar la solución a 16.000 x g durante 10 min. Recoger el sobrenadante y mezclar con 200 l de etanol al 70%.
  4. Utilice la columna adjunta para retirar el disolvente mediante centrifugación a 14.000 x g durante 1 min.
  5. Lavar con la solución de lavado 1 y 2 para eliminar completamente los NON-RNA por centrifugación a 14.000 x g durante 1 min.
  6. Centrifugar de nuevo a 14.000 x g durante 2 min para eliminar el etanol residual.
  7. Añadir 50 l de agua destilada, centrífuga a 14.000 x g durante 1 min, y recoger la solución.
  8. Mida la concentración de ARN utilizando OD260 con un espectrofotómetro.
    Concentración de ARN (g/ml) á (OD260) x (40 g de ARN/ml) y OD260/OD280 > 2. Las muestras de ARN deben tener un número de integridad de ARN (RIN) > 7.

3. Perfilado RNAseq y análisis bioinformático

NOTA: El análisis de RNAseq se realizó comercialmente (ver Tabla de Materiales)para investigar los genes diferenciales en las tumoresferas derivadas de HT29 en comparación con las células HT29 parentales.

  1. Utilice servicios comerciales para los pasos de RNAseq, incluida la construcción de bibliotecas, el control de calidad de la biblioteca y la secuenciación de ADN.
  2. El informe de datos debe contener información importante, incluidos los recuentos de lectura, el cambio de plegado log2 y el valor p. Seleccione los genes diferenciales de acuerdo con los siguientes parámetros: los genes con un cambio de plegado de > 1 log2 con recuentos de lectura >100 en el grupo de tumoresfera HT29, y los genes <-1 log2 cambian con el recuento de lectura >100 en el grupo parental HT29. En este caso, se consideró aceptable un valor p < 0.05 y se utilizaron los datos (Tabla 1).
    NOTA: Aquí, se utilizó un recuento de genes >100 como umbral para continuar el estudio de un gen en particular y validar su expresión.
  3. Utilice software de análisis estadístico (consulte Tabla de materiales), para mostrar un mapa de calor e identificar genes sobreexpresados >1 y genes regulados hacia abajo < 1 en el cambio de plegado de log2.
  4. Utilice el software R para dibujar la gráfica del volcán con x: cambio de pliegue log2; y: -log10 (valor p) para mostrar los genes diferenciales.
    1. Instalar la biblioteca R
      install.packages(library(calibrate))
    2. Lea los datos en RStudio con el siguiente programa:
      res <-read.csv("/Users/xxx.csv", header-T)
      cabeza(s)
      with(res, plot(log2FoldChange, -log10(pvalue), pch-19, main-"HT29CSC vs HT29", xlim-c(-6,6), col-"#C0C0C0"))
      with(subset(res, pvalue<.05 & log2FoldChange>1), points(log2FoldChange, -log10(pvalue), pch-19, col-"red"))
      with(subset(res, pvalue<.05 & log2FoldChange<(-1)), points(log2FoldChange, -log10(pvalue), pch-19, col-"blue"))
      with(subset(res, pvalue>.05), points(log2FoldChange, -log10(pvalue), pch-19, col-"#444444"))
      ablina (h-1,3, lty-2)
      abline(v-1, lty-2)
      abline(v-(-1), lty-2)
    3. Ejecute la carrera para obtener la gráfica del volcán.

4. Selección del gen del conductor

  1. Seleccione Entrada de un solo gen en NetworkAnalyst.
  2. Copie y pegue los genes sobreexpresados seleccionados de la Tabla 1 con "humano" especificado como el organismo y el tipo de identificación Símbolo genético oficial.
    NOTA: Alternativamente, utilice Ensembl Gene ID para copiar y pegar.
  3. Inserte los datos haciendo clic en Cargar y Proceder para analizarlos mediante la interacción proteína-proteína (IPP) después del IBP genético.
  4. Utilice la base de datos de interacto STRING con un límite de puntuación de confianza de 900 para mostrar los genes de semillas que cruzan los genes cargados. Los genes de semillas que se asocian con más genes individuales fueron seleccionados como genes conductores que pueden estar involucrados en el mantenimiento de la formación de tumores de origen derivado de HT29.
    NOTA: Hay tres conjuntos de datos interactome para su uso: IMEx, STRING y Rolland. STRING contiene evidencia experimental de mayor confianza. Con números de genes cargados más bajos, IMEx se puede seleccionar para predecir y recoger los genes del controlador en las redes interactome.
  5. Seleccione Continuar en la vista general de asignación.
  6. Seleccione Blanco en el Fondo y Forzar Atlas en el botón Diseño.
  7. Seleccione PANTHER BP para analizar el grupo genético de regulación ascendente.
    NOTA: Esto muestra que HSPA5 fue responsable de la antipoptosis en las tumoresferas derivadas de HT29 en este estudio (Figura 3A). Para reducir el campo funcional específico, la clasificación KEGG, GO o PANTHER se puede utilizar alternativamente para seleccionar los genes de controlador específicos.

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Resultados

Para establecer el modelo para investigar el mecanismo en las células madre cancerosas, se utilizaron células TI29 colorrectales para cultivar las tumoresferas similares al cáncer in vitro en una placa de baja unión que contiene B27, EGF, bFGF, HGF e IL6. Las tumoresferas >100 m de diámetro se formaron en 7 días(Figura 1A). Las tumoresferas se tripsinaronizaron a células individuales y se analizaron utilizando citometría de flujo para detectar la expresión LGR5 y CD133. LGR5 aumentó en la...

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Discusión

En este estudio, se utilizaron tumores de tipo tallo cultivado como modelo en el análisis de datos RNAseq con bioinformática disponible. Para un modelo de enfermedad, se utilizaron tumores timóferas derivadas de HT29. Debido a que las tumoresferas tienen resistencia a los medicamentos contra las terapias tumorales, el modelo establecido se puede utilizar para investigar los mecanismos detallados de resistencia mediante la investigación de las diferencias en la expresión génica. Además, la tecnología genómica que utiliza ...

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Divulgaciones

Los autores no tienen divulgaciones financieras relevantes.

Agradecimientos

Los autores agradecen el Laboratorio Básico de Biología radiológica del Instituto de Investigación Radiológica, Chang Gung Memorial Hospital, por su apoyo técnico. Este estudio fue apoyado por subvenciones del hospital Chang Gung Memorial (CMRPD1J0321), Cheng Hsin General Hospital (CHGH 106-06), y Mackay Memorial Hospital (MMH-CT-10605 y MMH-106-61). Los organismos de financiación no tuvieron ninguna influencia en el diseño del estudio y la recopilación, análisis e interpretación de datos o en la escritura del manuscrito.

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Materiales

Lista de materiales utilizados en este artículo
NombreEmpresaNúmero de catálogoComentarios
imagen celular iRiSLogos Biosystems, IncI10999para observar la formación de tumoreras
Citometría de flujoBD biosciencesFACSCaliburpara la detección de LGR5 y CD133 en las esferas tumorales
anti-LGR5-PEBiolegend373803reactivo de detección LGR5
anti-CD133-PEBiolegend372803Detección de CD133 reactivo
EGFGenScriptZ00333para el cultivo de tumores
bFGFGenScriptZ03116para el cultivo de tumoresferas
HGFGenScriptZ03229para el cultivo de tumoresferas
IL6GenScriptZ03034para el cultivo de tumoresferas
PureLink Kit de extracción de ARNInvitrogen 12183025 ARN total para el análisis de RNAseq
Rendimiento de RNAseqBiotools, TaiwánEl análisis de RNAseq se realiza comercialmente por Biotools, Ttaiwan
NetworkAnalystInstitute of Parasitology, McGill University, Montreal, Quebec, Canadáhttp://www.networkanalyst.ca/
PrismGraphPad Softwareun software de análisis estadístico
Sistema digital de ,

Referencias

  1. Rawla, P., Sunkara, T., Barsouk, A. Epidemiology of colorectal cancer: incidence, mortality, survival, and risk factors. Przegląd Gastroenterologiczny. 14 (2), 89-103 (2019).
  2. Wong, M. C., Ding, H., Wang, J., Chan, P. S., Huang, J. Prevalence and risk factors of colorectal cancer in Asia. Intestinal Research. 17 (3), 317-329 (2019).
  3. Liu, Q., et al. Positive expression of basic transcription factor 3 predicts poor survival of colorectal cancer patients: possible mechanisms involved. Cell Death & Disease. 10 (7), 509(2019).
  4. Slattery, M. L., et al. Dysregulated genes and miRNAs in the apoptosis pathway in colorectal cancer patients. Apoptosis. 23 (3-4), 237-250 (2018).
  5. Arteaga, C. L., Engelman, J. A. ERBB receptors: from oncogene discovery to basic science to mechanism-based cancer therapeutics. Cancer Cell. 25 (3), 282-303 (2014).
  6. Yarden, Y., Pines, G. The ERBB network: at last, cancer therapy meets systems biology. Nature Reviews Cancer. 12 (8), 553-563 (2012).
  7. Cheng, C. C., et al. YM155 as an inhibitor of cancer stemness simultaneously inhibits autophosphorylation of epidermal growth factor receptor and G9a-mediated stemness in lung cancer cells. PLoS One. 12 (8), 0182149(2017).
  8. Prasetyanti, P. R., Medema, J. P. Intra-tumor heterogeneity from a cancer stem cell perspective. Molecular Cancer. 16 (1), 41(2017).
  9. Zhao, Y., et al. CD133 expression may be useful as a prognostic indicator in colorectal cancer, a tool for optimizing therapy and supportive evidence for the cancer stem cell hypothesis: a meta-analysis. Oncotarget. 7 (9), 10023-10036 (2016).
  10. Choi, J. E., et al. Expression of epithelial-mesenchymal transition and cancer stem cell markers in colorectal adenocarcinoma: Clinicopathological significance. Oncology Reports. 38 (3), 1695-1705 (2017).
  11. Massard, C., Deutsch, E., Soria, J. C. Tumour stem cell-targeted treatment: elimination or differentiation. Annals of Oncology. 17 (11), 1620-1624 (2006).
  12. Grillet, F., et al. Circulating tumour cells from patients with colorectal cancer have cancer stem cell hallmarks in ex vivo culture. Gut. 66 (10), 1802-1810 (2017).
  13. Dallas, N. A., et al. Chemoresistant colorectal cancer cells, the cancer stem cell phenotype, and increased sensitivity to insulin-like growth factor-I receptor inhibition. Cancer Research. 69 (5), 1951-1957 (2009).
  14. Chang, Y. F., et al. STAT3 induces G9a to exacerbate HER3 expression for the survival of epidermal growth factor receptor-tyrosine kinase inhibitors in lung cancers. BMC Cancer. 19 (1), 959(2019).
  15. Catalano, V., et al. Colorectal cancer stem cells and cell death. Cancers (Basel). 3 (2), 1929-1946 (2011).
  16. Piggott, L., et al. Suppression of apoptosis inhibitor c-FLIP selectively eliminates breast cancer stem cell activity in response to the anti-cancer agent, TRAIL. Breast Cancer Research. 13 (5), 88(2011).
  17. Chung, S. Y., et al. Two novel SHP-1 agonists, SC-43 and SC-78, are more potent than regorafenib in suppressing the in vitro stemness of human colorectal cancer cells. Cell Death Discovery. 4, 25(2018).
  18. Cheng, C. C., et al. STAT3 exacerbates survival of cancer stem-like tumorspheres in EGFR-positive colorectal cancers: RNAseq analysis and therapeutic screening. Journal of Biomedical Science. 25 (1), 60(2018).
  19. Kleist, B., Xu, L., Li, G., Kersten, C. Expression of the adult intestinal stem cell marker Lgr5 in the metastatic cascade of colorectal cancer. International Journal of Clinical and Experimental Pathology. 4 (4), 327-335 (2011).
  20. Medema, J. P. Targeting the Colorectal Cancer Stem Cell. New England Journal of Medicine. 377 (9), 888-890 (2017).
  21. Sahlberg, S. H., Spiegelberg, D., Glimelius, B., Stenerlow, B., Nestor, M. Evaluation of cancer stem cell markers CD133, CD44, CD24: association with AKT isoforms and radiation resistance in colon cancer cells. PLoS One. 9 (4), 94621(2014).
  22. Xia, J., Gill, E. E., Hancock, R. E. NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data. Nature Protocols. 10 (6), 823-844 (2015).
  23. Gagan, J., Van Allen, E. M. Next-generation sequencing to guide cancer therapy. Genome Medicine. 7 (1), 80(2015).
  24. Panichnantakul, P., Bourgey, M., Montpetit, A., Bourque, G., Riazalhosseini, Y. RNA-Seq as a Tool to Study the Tumor Microenvironment. Methods in Molecular Biology. 1458, 311-337 (2016).
  25. Kanehisa, M., Sato, Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Science. 29 (1), 28-35 (2020).
  26. Ashburner, M., et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature Genetics. 25 (1), (2000).
  27. The Gene Ontology Collective. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Research. 47, 330-338 (2019).
  28. Mi, H., Muruganujan, A., Thomas, P. D. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Research. 41, Database issue 377-386 (2013).
  29. Yu, F., Shen, X. Y., Fan, L., Yu, Z. C. Genome-wide analysis of genetic variations assisted by Ingenuity Pathway Analysis to comprehensively investigate potential genetic targets associated with the progression of hepatocellular carcinoma. European Review for Medical and Pharmacological Sciences. 18 (15), 2102-2108 (2014).
  30. Zhou, G., et al. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Research. 47 (1), 234-241 (2019).
  31. Cheng, C. C., et al. Epidermal growth factor induces STAT1 expression to exacerbate the IFNr-mediated PD-L1 axis in epidermal growth factor receptor-positive cancers. Molecular Carcinogenesis. 57 (11), 1588-1598 (2018).
  32. Gallego Romero, I., Pai, A. A., Tung, J., Gilad, Y. RNA-seq: impact of RNA degradation on transcript quantification. BMC Biology. 12, 42(2014).

Acceso restringido. Inicie sesión o comience una prueba gratuita para ver este contenido.

Reimpresiones y permisos

Explorar más artículos

Células madre del cáncer colorrectaltumoresferas HT29análisis de RNA-Seqcitometría de flujomarcadores LGR5 CD133descubrimiento de genes controladorescribado bioinformáticointeracción proteína-proteínavalidación mediante PCR cuantitativaformación de tumoresferas