Descripción general
El empalme alternativo es un mecanismo fundamental mediante el cual un solo gen puede producir múltiples isoformas de ARNm, ampliando enormemente la capacidad de codificación del genoma. Casi el 95% de los genes multiexónicos humanos se someten a un empalme alternativo, generando productos proteicos distintos que pueden tener funciones diferentes, incluso opuestas. La elección de qué exones se incluyen o excluyen está estrictamente regulada y es específica de cada tejido. El análisis de empalme alternativo tiene como objetivo identificar y cuantificar estas diferencias a nivel de isoformas entre condiciones, etapas de desarrollo o estados patológicos, revelando una capa adicional de regulación genética más allá de los simples cambios de expresión.
Métodos
La detección de eventos de empalme alternativos a partir de datos de RNA-seq requiere enfoques computacionales especializados. Las herramientas de cuantificación de isoformas (como Salmon, Kallisto o RSEM) estiman la abundancia de cada isoforma de transcripción conocida. Las herramientas basadas en eventos (rMATS, MAJIQ o SUPPA2) clasifican los eventos de empalme en categorías: omisión de exones del casete, exones mutuamente excluyentes, sitios de empalme alternativos 5’ o 3’ y retención de intrones. Estas herramientas utilizan matrices de recuento de mapeo de lecturas para empalmar uniones versus exones y prueban diferencias significativas utilizando modelos estadísticos emparejados. Detección de empalme de novo (mediante StringTie o Cufflinks) reúne transcripciones sin anotaciones para descubrir nuevas isoformas. Visualizar las diferencias de empalme con gráficos de sashimi (de MISO o IGV) ayuda a interpretar patrones complejos.
Aplicaciones
La desregulación del empalme alternativo es un sello distintivo de muchas enfermedades. En el cáncer, los factores de empalme frecuentemente mutan, produciendo isoformas específicas de tumores que impulsan la proliferación y la metástasis. Los trastornos neurológicos como la atrofia muscular espinal y la demencia frontotemporal son el resultado directo de defectos de empalme. Comprender los patrones de empalme ayuda a interpretar la estructura y los tipos de ARN y se basa en conceptos de la transcripción y el procesamiento de ARN. El análisis de empalme también se cruza con la regulación genética y la epigenética, ya que el estado de la cromatina influye en la selección del sitio de empalme cotranscripcional, y la retroalimentación del empalme puede regular el alargamiento de la transcripción.
Protocolo práctico
Para el empalme diferencial basado en eventos con rMATS, comience con archivos BAM RNA-seq alineados por STAR usando el indicador --outSAMstrandField intronMotif para la alineación con reconocimiento de empalme. Prepare un archivo de lista BAM y una anotación GTF, luego ejecute rmats.py --b1 sample1.bam,sample2.bam --b2 control1.bam,control2.bam --gtf annotation.gtf --od output_dir -t paired --readLength 150. rMATS detecta cinco tipos de eventos de empalme: exón omitido (SE), exones mutuamente excluyentes (MXE), sitio de empalme 5’ alternativo (A5SS), sitio de empalme 3’ alternativo (A3SS) e intrón retenido (RI). Los resultados incluyen la diferencia de nivel de inclusión (IncLevelDifference) y una tasa de descubrimiento falso (FDR) para cada evento. La métrica clave es el valor de porcentaje empalmado (PSI o Ψ), que cuantifica la proporción de transcripciones que incluyen un exón. PSI varía de 0 (excluido) a 1 (incluido), y un ΔΨ > 0,1 con FDR < 0,05 es un umbral de significación común. Para la cuantificación a nivel de isoforma sin clasificación de eventos, Salmon en modo sin alineación cuantifica la abundancia de transcripciones: salmon quant -i transcriptome_index -l A -1 R1.fastq -2 R2.fastq -o isoform_output. Posteriormente, tximport importa estimaciones a nivel de transcripción en R, y DRIMSeq o satuRn prueban el uso diferencial de transcripciones. Visualice eventos de empalme con rmats2sashimiplot, que genera gráficos de sashimi que muestran la cobertura de lectura en las uniones exón-exón. Integre el análisis de empalme QTL (sQTL) utilizando “FastQTL” cuando haya datos de genotipo disponibles, vinculando variantes genéticas con la variación de empalme entre individuos.
