La secuenciación del genoma y epigenoma de células individuales mediante transposasas (scGET-seq) es un método de secuenciación de ADN para perfilar la cromatina abierta y cerrada . A diferencia del ensayo de células individuales para la cromatina accesible a transposasas con secuenciación (scATAC-seq), que solo se dirige a la eucromatina activa , [ 1 ] scGET-seq también es capaz de sondear la heterocromatina inactiva . [ 2 ]
Esto se logra mediante el uso de TnH, que se crea mediante la unión del cromodominio (CD) de la proteína de heterocromatina-1-alfa (HP-1) a la transposasa Tn5 . TnH puede entonces dirigirse a la trimetilación de la lisina 9 de la histona 3 ( H3K9me3 ), un marcador de heterocromatina. [ 3 ]
De forma similar a la velocidad del ARN, que utiliza la relación entre el ARN empalmado y el no empalmado para inferir la cinética de los cambios en la expresión génica a lo largo del desarrollo celular, [ 4 ] la relación entre las señales de TnH y Tn5 obtenidas mediante scGET-seq puede utilizarse para calcular la velocidad de la cromatina, que mide la dinámica de la accesibilidad de la cromatina a lo largo de las vías de desarrollo celular. [ 2 ]
Historia
La regulación transcripcional está estrechamente ligada a los estados de la cromatina. La cromatina abierta, o permisiva a la transcripción, constituye solo el 2-3% del genoma , pero abarca el 94,4% de los sitios de unión de factores de transcripción . [ 5 ] [ 6 ] Por el contrario, el ADN más compactado , o heterocromatina , es responsable de la organización y estabilidad del genoma. [ 7 ] La densidad de la cromatina también cambia a lo largo de los procesos de diferenciación celular , [ 8 ] pero hay una falta de métodos de secuenciación de alto rendimiento para analizar directamente la heterocromatina.
Muchas enfermedades relacionadas con el genoma, como el cáncer, están estrechamente vinculadas a cambios en su epigenoma . Los cánceres, en particular, se caracterizan por la heterogeneidad unicelular, que puede impulsar la metástasis y la resistencia al tratamiento. [ 9 ] [ 10 ] Los mecanismos subyacentes a estos procesos aún se desconocen en gran medida, aunque el advenimiento de las tecnologías unicelulares , incluida la epigenómica unicelular , ha contribuido enormemente a su elucidación. [ 11 ]
En 2015, ATAC-seq , que utiliza la transposasa Tn5 para fragmentar y marcar la cromatina accesible, o eucromatina , para su secuenciación, se volvió factible con resolución de célula individual. [ 12 ] scGET-seq se basa en esta tecnología al proporcionar también información sobre la heterocromatina , ofreciendo una visión más completa de la estructura y la dinámica de la cromatina dentro de cada célula. [ 13 ]
Métodos

Preparación de la muestra
La preparación de muestras para scGET-seq comienza con la obtención de una suspensión de núcleos de células utilizando un método apropiado para el material de partida. [ 14 ]
El siguiente paso es producir la transposasa TnH . Tn5 es una transposasa que corta y liga adaptadores a regiones genómicas no unidas por nucleosomas (cromatina abierta). [ 15 ] HP-1a es un miembro de la familia HP1 y es capaz de reconocer y unirse específicamente a H3K9me3 . [ 16 ] [ 17 ] Su cromodominio utiliza un mecanismo de ajuste inducido para reconocer esta modificación de la cromatina. [ 18 ] La unión de los primeros 112 aminoácidos de HP-1a que contienen el cromodominio a Tn5 mediante un enlazador de tres poli-tirosina-glicina-serina (TGS) conduce a la creación de la transposasa TnH , que es capaz de dirigirse a la heterocromatina marcada por H3K9me3 . [ 2 ]
La preparación de la biblioteca se realiza utilizando un protocolo modificado para ATAC-seq de célula única , [ 19 ] donde la suspensión de núcleos se incuba secuencialmente primero con la transposasa Tn5 y luego con TnH. [ 2 ]
Análisis de datos
Los objetivos del análisis de datos son: [ 2 ]
- Para identificar y caracterizar poblaciones celulares distintas mediante agrupamiento
- Para perfilar la accesibilidad de la cromatina en todo el genoma.
- Para predecir variantes de número de copias y variantes de un solo nucleótido.
Preprocesamiento
- Tras la secuenciación , las lecturas deben desmultiplexarse y mapearse al genoma de referencia correspondiente . Las lecturas duplicadas se identifican y eliminan.
- Se identifican "picos" o regiones en el ADN enriquecidas en el número de lecturas mapeadas. [ 20 ]
- Se realiza un control de calidad y se filtran las células con un número bajo de lecturas o con pocas características detectadas.
- Se generan cuatro matrices de conteo (matrices donde cada columna es una célula y cada fila es una característica): Tn5-dhs, Tn5-complemento, TnH-dhs y TnH-complemento, que representan la señal de la cromatina accesible y compactada. [ 2 ]
Análisis
Reducción de dimensionalidad, visualización y agrupamiento.
Cada una de las matrices se filtra para eliminar las regiones compartidas, se normaliza y se transforma mediante logaritmo en base 2. La reducción de dimensionalidad lineal se realiza mediante análisis de componentes principales (PCA). Los grupos de células se identifican mediante un algoritmo k-NN [ 21 ] y el algoritmo de Leiden [ 22 ] . Finalmente, las cuatro matrices se combinan mediante factorización matricial [ 23 ] y reducción UMAP [ 24 ] .
Anotación de identificación celular
Hay dos enfoques para la anotación de la identidad celular: Anotación basada en la anotación de características de los picos ATAC, [ 25 ] y anotación basada en la integración con datos de scRNA-seq de referencia. [ 26 ]
Aplicaciones

Actual
Utilizando la relación entre las señales de Tn5 y TnH, se pueden calcular valores cuantitativos que describen la rapidez y la dirección de la remodelación de la cromatina (velocidad de la cromatina). [ 2 ] Al aislar las regiones más dinámicas e identificar los factores de transcripción que se unen a ellas, la velocidad de la cromatina puede utilizarse para inferir los procesos epigenéticos dinámicos que ocurren dentro de una célula determinada y las contribuciones de diversos factores de transcripción a dichos procesos. [ 2 ]
Futuro
La remodelación de la cromatina precede a los cambios en la expresión génica y mejora la comprensión de las trayectorias y los mecanismos de los cambios celulares. [ 27 ] [ 28 ] Por lo tanto, las plataformas y herramientas para la integración de datos multimodales son áreas de investigación activa. [ 29 ] [ 30 ] [ 31 ] Se ha propuesto incorporar elementos temporales y de direccionalidad mediante la integración de la velocidad de la cromatina con la velocidad del ARN para revelar aún más información sobre las vías de diferenciación. [ 32 ] [ 33 ]
Limitaciones
scGET-seq tiene algunas de las mismas limitaciones que scATAC-seq. Ambos procesos requieren muestras de núcleos de células viables y una alta viabilidad celular. [ 13 ] Una baja viabilidad celular produce una alta contaminación de ADN de fondo que no representa con precisión las señales biológicas auténticas. Además, la escasez y la naturaleza ruidosa de los datos de scATAC-seq y scGET-seq dificultan el análisis, y aún no existe consenso sobre la mejor manera de gestionar estos datos. [ 34 ]
Otra limitación es que scGET-seq aún necesita la validación de los resultados de SNV mediante secuenciación masiva del genoma. Aunque existe una alta correlación de mutaciones entre la secuenciación masiva del exoma y los resultados de scGET-seq, scGET-seq no logra detectar todos los SNV del exoma. [ 2 ]
Referencias
- ↑ Yan F, Powell DR, Curtis DJ, Wong NC (febrero de 2020). "De las lecturas a la comprensión: una guía para autoestopistas sobre el análisis de datos ATAC-seq" . Genome Biology . 21 (1) 22. doi : 10.1186/s13059-020-1929-3 . PMC 6996192. PMID 32014034 .
- 1 2 3 4 5 6 7 8 9 Tedesco M, Giannese F, Lazarević D, Giansanti V, Rosano D, Monzani S, et al. (febrero de 2022). "Chromatin Velocity reveals epigenetic dynamics by single-cell profiling of heterochromatin and euchromatin". Nature Biotechnology . 40 (2): 235– 244. doi : 10.1038/s41587-021-01031-1 . hdl : 11368/3007419 . PMID 34635836 . S2CID 238637962 .
- ↑ Kouzarides T (febrero de 2007). "Modificaciones de la cromatina y su función" . Cell . 128 ( 4): 693–705 . doi : 10.1016/j.cell.2007.02.005 . PMID 17320507. S2CID 11691263 .
- ^ La Manno G, Soldatov R, Zeisel A, Braun E, Hochgerner H, Petukhov V, et al. (Agosto de 2018). "Velocidad del ARN de células individuales" . Naturaleza . 560 (7719): 494– 498. Bibcode : 2018Natur.560..494L . doi : 10.1038/s41586-018-0414-6 . PMC 6130801 . PMID 30089906 .
- ↑ Klemm SL, Shipony Z, Greenleaf WJ (abril de 2019). "Accesibilidad de la cromatina y el epigenoma regulador". Nature Reviews. Genetics . 20 (4): 207– 220. doi : 10.1038/s41576-018-0089-8 . PMID 30675018 . S2CID 59159906 .
- ↑ Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, et al. (septiembre de 2012). " El paisaje de cromatina accesible del genoma humano" . Nature . 489 (7414): 75–82 . Bibcode : 2012Natur.489...75T . doi : 10.1038/ nature11232 . PMC 3721348. PMID 22955617. S2CID 4304439 .
- ↑ Penagos-Puig A, Furlan-Magaril M (2020). "La heterocromatina como un importante impulsor de la organización del genoma" . Frontiers in Cell and Developmental Biology . 8 579137. doi : 10.3389/fcell.2020.579137 . PMC 7530337. PMID 33072761 .
- ↑ Golkaram M, Jang J, Hellander S, Kosik KS, Petzold LR (octubre de 2017). "El papel de la densidad de la cromatina en la heterogeneidad de la población celular durante la diferenciación de células madre" . Scientific Reports . 7 (1) 13307. Bibcode : 2017NatSR...713307G . doi : 10.1038/ s41598-017-13731-3 . PMC 5645312. PMID 29042584 .
- ↑ Dagogo-Jack I, Shaw AT (febrero de 2018). "Heterogeneidad tumoral y resistencia a las terapias contra el cáncer". Nature Reviews. Oncología Clínica . 15 (2): 81– 94. doi : 10.1038/nrclinonc.2017.166 . PMID 29115304. S2CID 2194691 .
- ↑ Lawson DA, Kessenbrock K, Davis RT, Pervolarakis N, Werb Z (diciembre de 2018). "Heterogeneidad tumoral y metástasis a nivel de célula individual" . Nature Cell Biology . 20 (12): 1349– 1360. doi : 10.1038/s41556-018-0236-7 . PMC 6477686. PMID 30482943 .
- ↑ Dai Z, Gu XY, Xiang SY, Gong DD, Man CF, Fan Y (noviembre de 2020). "Investigación y aplicación de la secuenciación de células individuales en la heterogeneidad tumoral y la resistencia a fármacos de las células tumorales circulantes" . Biomarker Research . 8 (1) 60. doi : 10.1186/s40364-020-00240-1 . PMC 7653877. PMID 33292625 .
- ↑ Pott S, Lieb JD (agosto de 2015). "Secuenciación ATAC de células individuales: la fuerza reside en los números" . Genome Biology . 16 (1) 172. doi : 10.1186/s13059-015-0737-7 . PMC 4546161. PMID 26294014 .
- 1 2 Tang L (diciembre de 2021). "Esbozo de la cromatina abierta y cerrada". Nature Methods . 18 (12): 1448. doi : 10.1038/s41592-021-01351-9 . PMID 34862496 . S2CID 244871731 .
- ↑ "Aislamiento de núcleos para secuenciación de ARN de células individuales y tejidos para secuenciación de ARN de células individuales: protocolo demostrado, preparación de muestras, expresión génica de células individuales, soporte oficial de 10x Genomics" . support.10xgenomics.com . Consultado el 2 de marzo de 2022 .
- ↑ Hsu FM, Gohain M, Chang P, Lu JH, Chen PY (enero de 2018). «Capítulo 4 - Bioinformática de datos epigenómicos generados a partir de secuenciación de nueva generación». En Tollefsbol TO (ed.). Epigenética en enfermedades humanas . Epigenética traslacional. Vol. 6 (segunda ed.). Academic Press. págs. 65–106 . doi : 10.1016/B978-0-12-812215-0.00004-2 . ISBN 978-0-12-812215-0.
- ↑ Bannister AJ, Zegerman P, Partridge JF, Miska EA, Thomas JO, Allshire RC, Kouzarides T (marzo de 2001). "Reconocimiento selectivo de la lisina 9 metilada en la histona H3 por el dominio cromo HP1". Nature . 410 ( 6824): 120– 124. Bibcode : 2001Natur.410..120B . doi : 10.1038/35065138 . PMID 11242054. S2CID 4334447 .
- ↑ Watanabe S, Mishima Y, Shimizu M, Suetake I, Takada S (mayo de 2018). "Interacciones de HP1 unido al dinucleosoma H3K9me3 mediante simulaciones moleculares y ensayos bioquímicos" . Biophysical Journal . 114 (10): 2336– 2351. Bibcode : 2018BpJ...114.2336W . doi : 10.1016/j.bpj.2018.03.025 . PMC 6129468. PMID 29685391 .
- ↑ Nielsen PR, Nietlispach D, Mott HR, Callaghan J, Bannister A, Kouzarides T, et al. (marzo de 2002). "Estructura del cromodominio HP1 unido a la histona H3 metilada en la lisina 9". Nature . 416 ( 6876): 103– 107. Bibcode : 2002Natur.416..103N . doi : 10.1038/nature722 . PMID 11882902. S2CID 4423019 .
- ↑ "Guía del usuario de los kits de reactivos ATAC de células individuales de cromo (v1.1 Química) - Guía del usuario - Soporte oficial de 10x Genomics" . support.10xgenomics.com . Consultado el 2 de marzo de 2022 .
- ↑ Baek S, Lee I (2020-01-01). "Análisis de secuenciación ATAC de células individuales: del preprocesamiento de datos a la generación de hipótesis" . Computational and Structural Biotechnology Journal . 18 : 1429–1439 . doi : 10.1016/j.csbj.2020.06.012 . PMC 7327298. PMID 32637041 .
- ↑ Polański K, Young MD, Miao Z, Meyer KB, Teichmann SA, Park JE (febrero de 2020). "BBKNN: alineación rápida por lotes de transcriptomas de células individuales" . Bioinformatics . 36 ( 3): 964–965 . doi : 10.1093/bioinformatics/btz625 . PMC 9883685. PMID 31400197 .
- ^ Traag VA, Waltman L, van Eck Nueva Jersey (marzo de 2019). "De Lovaina a Leiden: garantizar comunidades bien conectadas" . Informes científicos . 9 (1) 5233. arXiv : 1810.08473 . Código Bib : 2019NatSR...9.5233T . doi : 10.1038/s41598-019-41695-z . PMC 6435756 . PMID 30914743 .
- ↑ Žitnik M, Zupan B (enero de 2015). "Fusión de datos mediante factorización matricial". IEEE Transactions on Pattern Analysis and Machine Intelligence . 37 (1): 41– 53. arXiv : 1307.0803 . Bibcode : 2015ITPAM..37...41I . doi : 10.1109 / TPAMI.2014.2343973 . PMID 26353207. S2CID 362295 .
- ↑ "UMAP: Aproximación y proyección de variedades uniformes para la reducción de dimensión — documentación de umap 0.5" . umap-learn.readthedocs.io . Consultado el 4 de marzo de 2022 .
- ↑ Cittaro D (2022-02-21), dawe/scatACC , consultado el 2022-03-04
- ↑ Abdelaal T, Michielsen L, Cats D, Hoogduin D, Mei H, Reinders MJ, Mahfouz A (septiembre de 2019). " Una comparación de métodos de identificación automática de células para datos de secuenciación de ARN de células individuales" . Genome Biology . 20 (1) 194. doi : 10.1186/s13059-019-1795-z . PMC 6734286. PMID 31500660 .
- ↑ Stadhouders R, Vidal E, Serra F, Di Stefano B, Le Dily F, Quilez J, et al. (febrero de 2018). " Los factores de transcripción orquestan la interacción dinámica entre la topología del genoma y la regulación génica durante la reprogramación celular" . Nature Genetics . 50 (2): 238– 249. doi : 10.1038/s41588-017-0030-7 . PMC 5810905. PMID 29335546 .
- ↑ Ranzoni AM, Tangherloni A, Berest I, Riva SG, Myers B, Strzelecka PM, et al. (marzo de 2021). "Análisis integrador de ARN-Seq de célula única y ATAC-Seq de la hematopoyesis del desarrollo humano" . Cell Stem Cell . 28 (3): 472–487.e7. doi : 10.1016/j.stem.2020.11.015 . PMC 7939551. PMID 33352111 .
- ↑ Lin Y, Wu TY, Wan S, Yang JY, Wong WH, Wang YX (enero de 2022). "scJoint integra datos de secuenciación de ARN de células individuales y ATAC-seq a escala de atlas con aprendizaje por transferencia" . Nature Biotechnology . 40 ( 5): 703–710 . doi : 10.1038/s41587-021-01161-6 . PMC 9186323. PMID 35058621. S2CID 246150572 .
- ↑ Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. (junio de 2019). " Integración integral de datos de células individuales" . Cell . 177 (7): 1888–1902.e21. doi : 10.1016/j.cell.2019.05.031 . PMC 6687398. PMID 31178118 .
- ↑ Wang C, Sun D, Huang X, Wan C, Li Z, Han Y, et al. (agosto de 2020). "Análisis integradores del transcriptoma y reguloma de células individuales utilizando MAESTRO" . Genome Biology . 21 (1) 198. doi : 10.1186/s13059-020-02116-x . PMC 7412809. PMID 32767996 .
- ↑ Xu Y, Begoli E, McCord RP (2022). "SciCAN: Integración de datos de accesibilidad de cromatina de células individuales y expresión génica mediante una red adversaria consistente en el ciclo" . npj Systems Biology and Applications . 8 (1) 33. bioRxiv 10.1101/2021.11.30.470677 . doi : 10.1038/s41540-022-00245-6 . PMC 9464763. PMID 36089620 .
- ↑ Chen Z, King WC, Gerstein M, Zhang J (2022). " DeepVelo : Aprendizaje profundo del campo de velocidad transcriptómico de células individuales con ecuaciones diferenciales ordinarias neuronales" . Science Advances . 8 (48) eabq3745. bioRxiv 10.1101/2022.02.15.480564 . doi : 10.1126/sciadv.abq3745 . PMC 9710871 .
- ↑ Baek S, Lee I (enero de 2020). "Análisis de secuenciación ATAC de células individuales: del preprocesamiento de datos a la generación de hipótesis" . Computational and Structural Biotechnology Journal . 18 : 1429–1439 . doi : 10.1016/j.csbj.2020.06.012 . PMC 7327298. PMID 32637041 .
- Técnicas de biología molecular