Articulo de referencia

Secuenciación de ARN con resolución temporal

Los métodos de secuenciación de ARN con resolución temporal son aplicaciones de la secuenciación de ARN que permiten observar la abundancia de ARN a lo largo del tiempo en una o...

Los métodos de secuenciación de ARN con resolución temporal son aplicaciones de la secuenciación de ARN que permiten observar la abundancia de ARN a lo largo del tiempo en una o varias muestras biológicas. La secuenciación de ADN de segunda generación ha posibilitado un análisis del transcriptoma rentable, de alto rendimiento e imparcial . [ 1 ] Normalmente, la secuenciación de ARN solo permite capturar una instantánea del transcriptoma en el momento de la toma de la muestra. [ 1 ] Esto requiere múltiples muestreos en diferentes momentos, lo que aumenta los costes económicos y de tiempo de los experimentos. Las innovaciones metodológicas y tecnológicas han permitido analizar el transcriptoma de ARN a lo largo del tiempo sin necesidad de realizar múltiples muestreos en distintos momentos.

Fondo

Si bien el ADN codifica todos los elementos funcionales de la vida, la información codificada debe convertirse en una forma funcional. Siguiendo el dogma central de la biología molecular , el ARN mensajero codifica la información genética para la producción de proteínas, las cuales, junto con el ARN funcional, llevan a cabo la mayoría de los procesos celulares necesarios para la vida. [ 2 ] Los cambios en la abundancia de ARN pueden utilizarse como una medida de cambios en el comportamiento celular, como el estrés térmico , la infección por virus o la oncogénesis . [ 3 ] El conocimiento de cómo cambia el transcriptoma durante los procesos celulares permite una mayor comprensión de los mecanismos exactos que subyacen a estos procesos.

Originalmente, la abundancia de ARN en todo el transcriptoma solo podía evaluarse utilizando métodos como microarrays de ADN o análisis serial de expresión génica (SAGE). [ 4 ] [ 5 ] Estos métodos son prohibitivos en diferentes aspectos; los microarrays, si bien son baratos, proporcionan resultados inconsistentes [ 6 ] y SAGE se basa en la secuenciación de Sanger , que proporciona un rendimiento limitado. Utilizando la secuenciación de segunda generación, en lugar de medir la hibridación relativa de secuencias a sondas en el caso de los microarrays o secuenciar segmentos cortos en el caso de SAGE, un investigador puede simplemente secuenciar el ARN total dentro de una muestra y medir las abundancias relativas de tipos específicos de ARN comparando la cantidad de veces que cada molécula de ARN fue secuenciada en una muestra dada.

Normalmente, en un experimento tradicional de RNA-seq, microarrays o SAGE, el ARN se extrae de una muestra biológica, como células cultivadas , y se analiza mediante el método elegido. Los datos obtenidos corresponden a la abundancia de ARN en las condiciones experimentales dadas en el momento de la extracción. Para muchas aplicaciones, como la comparación de la abundancia de moléculas de ARNm entre células expuestas a un fármaco y las que no lo están, este tipo de enfoque experimental es suficiente. Sin embargo, muchos procesos celulares de interés científico y médico son procesos que ocurren a lo largo del tiempo, como la diferenciación celular o la fagocitosis . [ 7 ] [ 8 ] El estudio de estos procesos requiere el análisis de la abundancia de ARN en una serie de puntos temporales.

Métodos

Comparación de métodos de secuenciación de ARN con resolución temporal. Las series temporales requieren muestras tanto anteriores como posteriores a todos los puntos temporales, así como la secuenciación por separado de todas las muestras biológicas. La purificación por afinidad reduce el número de muestras biológicas necesarias, pero aumenta el número de ejecuciones de secuenciación. La conversión de nucleótidos requiere el menor número de muestras biológicas y ejecuciones de secuenciación en general. h = horas.

Muestras de series temporales

Preparación de muestras y procesamiento de datos

El método más sencillo para evaluar la abundancia de ARN a lo largo del tiempo consiste en utilizar múltiples muestras tratadas exactamente igual, con la única diferencia de la duración del tratamiento. Por ejemplo, para investigar un proceso biológico que se estima que dura una hora, un investigador podría diseñar un experimento en el que el proceso se active durante cinco minutos, 15 minutos, 30 minutos, 45 minutos, una hora y dos horas en muestras de cultivo celular separadas antes de recolectar las células para el análisis de secuenciación de ARN (RNA-seq). De esta forma, el investigador obtendría mediciones del transcriptoma en cada uno de estos puntos temporales, y la comparación entre estas muestras indicaría qué procesos celulares se activan y desactivan con el tiempo.

Fortalezas

Este método es el más común para la medición de ARN a lo largo del tiempo en modelos de cultivo celular, principalmente debido a su simplicidad. Cada muestra biológica solo necesita procesarse exactamente de la misma manera, y el factor tiempo se ajusta fácilmente en la mayoría de los protocolos experimentales. Además, dado que cada punto temporal constituye una muestra independiente, se puede recolectar y secuenciar una mayor cantidad de ARN para un estudio.

Debilidades

La necesidad de múltiples muestras para la recopilación de datos con resolución temporal aumenta el costo del experimento e introduce un mayor potencial de errores técnicos. Si bien el precio de la secuenciación masivamente paralela ha disminuido considerablemente desde su introducción, sigue siendo prohibitivamente caro para muchos laboratorios realizar estudios de RNA-seq a gran escala. Este problema se agrava al añadir puntos temporales, lo que incrementa el número de muestras en un múltiplo del número de puntos temporales; usar dos puntos temporales en lugar de uno duplica el número de muestras necesarias en un experimento. En consecuencia, muchos estudios que utilizan RNA-seq de series temporales se ven limitados, ya sea por el tamaño de la muestra, lo que reduce la potencia estadística , [ 9 ] o por el número de puntos temporales, lo que reduce su resolución temporal, o por ambos. Finalmente, al requerir un mayor número de muestras biológicas, existe un mayor riesgo de que el error humano afecte los resultados, lo que puede llevar a conclusiones erróneas [ 10 ] [ 11 ].

Purificación por afinidad

Comparación de flujos de trabajo de marcaje metabólico.

Preparación de muestras y procesamiento de datos

En este enfoque, las muestras de cultivo celular se cultivan con nucleótidos marcados que permiten la purificación selectiva de moléculas de ARN recién sintetizadas. Un enfoque popular es el marcaje pulsado con 4-tiouridina (4-sU), un análogo del uracilo que se incorpora a las moléculas de ARN recién sintetizadas. [ 12 ] En este tipo de experimento, un investigador suplementaría las células con 4-sU en el momento del experimento o poco antes. Cuando el tratamiento experimental presumiblemente afecta la expresión del ARN, el ARN recién sintetizado se marcaría con 4-sU. El ARN recién sintetizado se marca con un grupo tiol reactivo , lo que permite unir moléculas útiles al ARN. [ 13 ] La biotina es una molécula popular para usar en este tipo de ensayo, ya que es económica y se une increíblemente fuerte y selectivamente a la estreptavidina . La incubación de ARN biotinilado con perlas que contienen estreptavidina permite la purificación selectiva del ARN recién sintetizado. A partir de aquí, el ARN recién sintetizado y el ARN total se secuencian por separado y se comparan para detectar diferencias.

Fortalezas

Affinity purification makes use of the incredibly popular biotin-streptavidin system for fractionation of biological materials. Binding of biotin to streptavidin is incredibly strong (Kd < 10−14 mol/L).[14] It is also highly specific, which results in minimal background signal from non-specific binding events. Furthermore, time resolution is obtained in a single biological sample, resulting in reduced biological variability compared to using separate samples for each time point.

Weaknesses

The weaknesses of this method are mainly centered around efficiency. One major difficulty is uptake of 4-sU into cultured cells. If 4-sU is given too early, then it will be incorporated into RNA that was not synthesized before the cell began responding to the experimental conditions. If it is given too late, then early stages of the cellular response are not captured by the experiment. The rate of uptake of 4-sU can be measured, but this requires additional experiments to determine optimal dosage and time. Furthermore, these parameters need to be measured in the specific cell lines of interest, as different cell lines may take up 4-sU more slowly than others. RNA is known to be prone to degradation in vitro. It is common for experimental protocols involving RNA to include a number of steps to reduce chances of Ribonuclease contamination or spontaneous degradation of samples, as RNA quality affects RNA-seq results.[15] Metabolic labeling involves a number of additional steps that must be performed in the laboratory on RNA that is in solution. Since metabolic labeling requires that the RNA be kept unfrozen in liquid solution, some level of spontaneous degradation is unavoidable, although it is usually not to such an extent that results are affected. Of greater risk is the chances of ribonuclease contamination, which would render a sample useless, wasting time and resources. It is important for researchers working with RNA in any capacity to minimize unnecessary handling of RNA due to these risks. One additional drawback of using this method is, given equivalent sample size, more sequencing runs are required compared to a time-series experiment. This is because multiple RNA samples corresponding to the initial time point must be sequenced.

Research suggests that 4-sU labeling may result in transcriptional changes on its own, which would affect any results obtained using this method.[16]

Nucleotide conversion

Sample preparation and data processing

La conversión de nucleótidos funciona transformando algunos nucleótidos del ARN recién sintetizado en otros, que pueden detectarse mediante secuenciación. SLAMseq y Timelapse-seq son ejemplos de estos enfoques. [ 17 ] [ 18 ] Al igual que en la purificación por afinidad, las células se incuban con 4-sU. Tras la extracción del ARN de las muestras, se tratan con yodoacetamida (SLAMseq) o 2,2,2-trifluoroetilamina y periodato de sodio (Timelapse-seq), lo que convierte el 4-sU en un análogo de citosina que se secuencia como un nucleótido de citosina en lugar de uracilo. Durante el alineamiento de secuencias y el procesamiento de datos, las conversiones de U a C se utilizan para cuantificar el número de transcritos recién sintetizados en comparación con el ARN total. [ 19 ]

Fortalezas

Este método comparte muchas ventajas con la purificación por afinidad; en particular, el hecho de que no se requieren múltiples muestras para una serie temporal. Este método elimina la necesidad de múltiples ejecuciones de secuenciación para múltiples puntos temporales, ya que todo el ARN se procesa conjuntamente en el instrumento de secuenciación y el ARN marcado se separa del no marcado in silico . Esto reduce significativamente los costos de secuenciación, ya que ahora se puede obtener resolución temporal sin necesidad de muestras adicionales ni ejecuciones de secuenciación adicionales. Además, al secuenciar múltiples puntos temporales simultáneamente, se reduce aún más la variabilidad técnica introducida por el procesamiento de muestras, además de la reducida variabilidad biológica que proporciona la estrategia experimental 4-sU.

Debilidades

Al igual que sus ventajas, este método comparte muchas desventajas con los métodos de purificación por afinidad. En particular, la captación de 4-sU y el mayor manejo de las muestras. Dado que Timelapse-seq se basa en métodos de química sintética para convertir nucleótidos, las reacciones incompletas conllevan una subestimación de la abundancia de ARN recién sintetizado y pueden generar variabilidad entre muestras.

Secuenciación de transcripciones nacientes

Preparación de muestras y procesamiento de datos

A diferencia del marcaje metabólico, la secuenciación de transcripciones nacientes (NET-seq) secuencia directamente las transcripciones que aún están siendo transcritas por la ARN polimerasa II . [ 20 ] Este método permite estudiar la dinámica de la elongación de la transcripción , lo cual no es posible con las técnicas de marcaje metabólico. Para un experimento de NET-seq, las células se tratan como en un experimento estándar de RNA-seq hasta que se lisan . La lisis se realiza de manera que los complejos ARN-proteína permanezcan intactos, y la ARN polimerasa II se inmunoprecipita del lisado. El ARN que estaba siendo transcrito a partir del ADN todavía está unido a la ARN polimerasa y posteriormente se eluye de la polimerasa y se secuencia.

Fortalezas

Dado que NET-seq extrae transcritos que no han completado su transcripción, es posible obtener resolución de un solo nucleótido en el nucleótido sintetizado más recientemente. Esto es valioso para el estudio de fenómenos como la cinética transcripcional . Además, permite el estudio de transcritos inestables que se degradan poco después de la transcripción. El enfoque general de inmunoprecipitación de proteínas de unión a ARN es de gran utilidad para comprender otras áreas de la biología del ARN, como el empalme .

Debilidades

Este método se basa en la inmunoprecipitación de la ARN polimerasa II. La inmunoprecipitación presenta varios problemas, incluidas las interacciones de unión no específicas, que pueden provocar la inmunoprecipitación de moléculas de ARN no deseadas. La resolución temporal de NET-seq se limita a la elongación de la transcripción. Si bien es posible comparar las abundancias relativas entre transcritos mediante NET-seq, este no es el objetivo del método.

Direcciones futuras

Aparte del muestreo de series temporales, actualmente no existen métodos para comparar más de dos puntos temporales. Los experimentos de marcaje metabólico solo permiten comparar la abundancia de ARN antes y después del marcaje por pulsos. Resulta interesante poder observar modificaciones en el transcriptoma a lo largo de varios puntos temporales en una misma muestra, ya que esto proporcionaría una mayor resolución temporal en los estudios. Los métodos de marcaje metabólico existentes son de interés para este fin; si se utilizaran diferentes marcadores metabólicos en distintos puntos temporales, se podrían investigar puntos temporales intermedios. Sin embargo, estos enfoques deben desarrollarse con cuidado, ya que los sesgos en los métodos de marcaje y en los pasos de procesamiento de las muestras podrían generar resultados erróneos si se comparan datos de diferentes métodos.

Se ha informado que el marcaje metabólico con 4-sU afecta el fenotipo celular . [ 16 ] En la práctica actual, esto es inevitable y se tolera ya que los datos obtenidos aún se ajustan a los modelos biológicos actuales, así como al hecho de que las muestras de 4-sU se comparan con muestras de 4-sU en la mayoría de los casos. Sin embargo, esto puede dar lugar a conclusiones espurias, especialmente si existe alguna interacción entre el efecto del 4-sU y la condición experimental elegida. No es posible distinguir las diferencias en los niveles de ARN como debidas a las condiciones experimentales estudiadas o como resultado del tratamiento con 4-sU. La identificación de sustancias químicas de marcaje que no afecten el fenotipo celular eliminaría por completo estos problemas.

Referencias

  1. 1 2 Wang Z, Gerstein M, Snyder M (enero de 2009). "RNA-Seq: una herramienta revolucionaria para la transcriptómica" . Nature Reviews. Genetics . 10 (1): 57– 63. doi : 10.1038/nrg2484 . PMC 2949280. PMID 19015660 .  
  2. Crick FH (1958). "Sobre la síntesis de proteínas". En Sanders FK (ed.). Simposios de la Sociedad de Biología Experimental, Número XII: La replicación biológica de macromoléculas . Cambridge University Press. págs. 138–163 . 
  3. Ozsolak F, Milos PM (febrero de 2011). " Secuenciación de ARN: avances, desafíos y oportunidades" . Nature Reviews. Genetics . 12 (2): 87– 98. doi : 10.1038/nrg2934 . PMC 3031867. PMID 21191423 .  
  4. Schena M, Shalon D, Davis RW, Brown PO (octubre de 1995). "Monitoreo cuantitativo de patrones de expresión génica con un microarreglo de ADN complementario". Science . 270 ( 5235): 467– 470. Bibcode : 1995Sci...270..467S . doi : 10.1126/science.270.5235.467 . PMID 7569999. S2CID 6720459 .  
  5. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW (octubre de 1995). "Análisis serial de la expresión génica". Science . 270 ( 5235): 484– 487. Bibcode : 1995Sci...270..484V . doi : 10.1126/science.270.5235.484 . PMID 7570003. S2CID 16281846 .  
  6. Okoniewski MJ, Miller CJ (junio de 2006). "Las interacciones de hibridación entre conjuntos de sondas en microarrays de oligos cortos conducen a correlaciones espurias" . BMC Bioinformatics . 7 : 276. doi : 10.1186/1471-2105-7-276 . PMC 1513401. PMID 16749918 .  
  7. Spies D, Ciaudo C (2015). " Dinámica en transcriptómica: avances en el análisis de la secuenciación de ARN y el análisis posterior" . Computational and Structural Biotechnology Journal . 13 : 469–477 . doi : 10.1016/j.csbj.2015.08.004 . PMC 4564389. PMID 26430493 .  
  8. Hejblum BP, Skinner J, Thiébaut R (junio de 2015). "Análisis de conjuntos de genes a lo largo del tiempo para datos de expresión génica longitudinal" . PLOS Computational Biology . 11 (6) e1004310. Bibcode : 2015PLSCB..11E4310H . doi : 10.1371/journal.pcbi.1004310 . PMC 4482329. PMID 26111374 .  
  9. Schurch NJ, Schofield P, Gierliński M, Cole C, Sherstnev A, Singh V, et al. (junio de 2016). "¿Cuántas réplicas biológicas se necesitan en un experimento de RNA-seq y qué herramienta de expresión diferencial se debe usar?" . RNA . 22 (6): 839– 851. doi : 10.1261/rna.053959.115 . PMC 4878611 . PMID 27022035 .   
  10. McIntyre LM, Lopiano KK, Morse AM, Amin V, Oberg AL, Young LJ, Nuzhdin SV (junio de 2011). " RNA-seq: variabilidad técnica y muestreo" . BMC Genomics . 12 : 293. doi : 10.1186/1471-2164-12-293 . PMC 3141664. PMID 21645359 .  
  11. Liu Y, Zhou J, White KP (febrero de 2014). "Estudios de expresión diferencial de RNA-seq: ¿más secuencia o más replicación?" . Bioinformatics . 30 ( 3): 301– 304. doi : 10.1093/bioinformatics/btt688 . PMC 3904521. PMID 24319002 .  
  12. Dölken L, Ruzsics Z, Rädle B, Friedel CC, Zimmer R, Mages J, et al. (septiembre de 2008). "Perfil de expresión génica de alta resolución para el análisis simultáneo de parámetros cinéticos de la síntesis y degradación del ARN" . RNA . 14 (9): 1959–1972 . doi : 10.1261/rna.1136108 . PMC 2525961. PMID 18658122 .   
  13. Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. (mayo de 2011). " Cuantificación global del control de la expresión génica en mamíferos" (PDF) . Nature . 473 (7347): 337– 342. Bibcode : 2011Natur.473..337S . doi : 10.1038/nature10098 . PMID 21593866. S2CID 205224972 .   
  14. Green NM (1975). "Avidina". Advances in Protein Chemistry . 29 : 85–133 . doi : 10.1016/S0065-3233(08)60411-8 . ISBN 978-0-12-034229-7. PMID 237414 . 
  15. Gallego Romero I, Pai AA, Tung J, Gilad Y (mayo de 2014). "RNA-seq: impacto de la degradación del ARN en la cuantificación de transcritos" . BMC Biology . 12 : 42. doi : 10.1186/1741-7007-12-42 . PMC 4071332. PMID 24885439 .  
  16. 1 2 Burger K, Mühl B, Kellner M, Rohrmoser M, Gruber-Eber A, Windhager L, et al. (octubre de 2013). "La 4-tiouridina inhibe la síntesis de ARNr y provoca una respuesta de estrés nucleolar" . RNA Biology . 10 (10): 1623– 1630. doi : 10.4161/rna.26214 . PMC 3866244. PMID 24025460 .   
  17. Herzog VA, Reichholf B, Neumann T, Rescheneder P, Bhat P, Burkard TR, et al. (diciembre de 2017). " Alquilación de ARN ligada a tiol para evaluar la dinámica de la expresión" . Nature Methods . 14 (12): 1198– 1204. doi : 10.1038/nmeth.4435 . PMC 5712218. PMID 28945705 .   
  18. Schofield JA, Duffy EE, Kiefer L, Sullivan MC, Simon MD (marzo de 2018). "TimeLapse-seq: agregando una dimensión temporal a la secuenciación de ARN a través de la recodificación de nucleósidos" . Nature Methods . 15 (3): 221– 225. doi : 10.1038/nmeth.4582 . PMC 5831505. PMID 29355846 .  
  19. Neumann T, Herzog VA, Muhar M, von Haeseler A, Zuber J, Ameres SL, Rescheneder P (mayo de 2019). "Cuantificación de conversiones de nucleótidos inducidas experimentalmente en conjuntos de datos de secuenciación de alto rendimiento" . BMC Bioinformatics . 20 (1): 258. doi : 10.1186/s12859-019-2849-7 . PMC 6528199. PMID 31109287 .  
  20. Churchman LS, Weissman JS (enero de 2011). " La secuenciación de transcripciones nacientes visualiza la transcripción con resolución de nucleótidos" . Nature . 469 (7330): 368–373 . Bibcode : 2011Natur.469..368C . doi : 10.1038/nature09652 . PMC 3880149. PMID 21248844 .