Visão geral
O splicing alternativo é um mecanismo fundamental pelo qual um único gene pode produzir múltiplas isoformas de mRNA, expandindo enormemente a capacidade de codificação do genoma. Quase 95% dos genes humanos multi-éxons passam por splicing alternativo, gerando produtos proteicos distintos que podem ter funções diferentes, até mesmo opostas. A escolha de quais exons são incluídos ou excluídos é rigorosamente regulada e específica do tecido. A análise de splicing alternativo visa identificar e quantificar essas diferenças no nível de isoformas entre condições, estágios de desenvolvimento ou estados de doença, revelando uma camada adicional de regulação genética além de simples alterações de expressão.
Métodos
A detecção de eventos de splicing alternativos a partir de dados de RNA-seq requer abordagens computacionais especializadas. Ferramentas de quantificação de isoformas (como Salmon, Kallisto ou RSEM) estimam a abundância de cada isoforma transcrita conhecida. Ferramentas baseadas em eventos (rMATS, MAJIQ ou SUPPA2) classificam os eventos de splicing em categorias: salto de exon de cassete, éxons mutuamente exclusivos, locais de splicing alternativos de 5’ ou 3’ e retenção de íntrons. Essas ferramentas usam matrizes de contagem de mapeamento de leituras para unir junções versus exons e testar diferenças significativas usando modelos estatísticos emparelhados. Detecção de emenda de novo (por StringTie ou Cufflinks) reúne transcrições sem anotação para descobrir novas isoformas. A visualização das diferenças de emenda com gráficos de sashimi (de MISO ou IGV) ajuda a interpretar padrões complexos.
Aplicativos
A desregulação do splicing alternativo é uma marca registrada de muitas doenças. No câncer, os fatores de splicing sofrem frequentemente mutação, produzindo isoformas específicas do tumor que impulsionam a proliferação e a metástase. Distúrbios neurológicos, como atrofia muscular espinhal e demência frontotemporal, resultam diretamente de defeitos de emenda. A compreensão dos padrões de splicing ajuda na interpretação de estrutura e tipos de RNA e se baseia em conceitos de transcrição e processamento de RNA. A análise de splicing também se cruza com regulação genética e epigenética, já que o estado da cromatina influencia a seleção do local de splice cotranscricional, e o feedback de splicing pode regular o alongamento da transcrição.
Protocolo Prático
Para emenda diferencial baseada em eventos com rMATS, comece com arquivos RNA-seq BAM alinhados por STAR usando o sinalizador --outSAMstrandField intronMotif para alinhamento com reconhecimento de emenda. Prepare um arquivo de lista BAM e uma anotação GTF e execute rmats.py --b1 sample1.bam,sample2.bam --b2 control1.bam,control2.bam --gtf annotation.gtf --od output_dir -t paired --readLength 150. O rMATS detecta cinco tipos de eventos de splicing: éxon ignorado (SE), éxons mutuamente exclusivos (MXE), local de splice alternativo 5’ (A5SS), local de splice alternativo 3’ (A3SS) e íntron retido (RI). Os resultados incluem diferença no nível de inclusão (IncLevelDifference) e uma taxa de descoberta falsa (FDR) para cada evento. A métrica principal é o valor Percent Spliced In (PSI ou Ψ), que quantifica a proporção de transcrições que incluem um éxon. O PSI varia de 0 (excluído) a 1 (incluído), e um ΔΨ > 0,1 com FDR < 0,05 é um limite de significância comum. Para quantificação em nível de isoforma sem classificação de eventos, Salmon no modo livre de alinhamento quantifica abundâncias de transcritos: salmon quant -i transcriptome_index -l A -1 R1.fastq -2 R2.fastq -o isoform_output. Downstream, tximport importa estimativas em nível de transcrição para R e teste DRIMSeq ou satuRn para uso diferencial de transcrição. Visualize eventos de splicing com rmats2sashimiplot, que gera gráficos de sashimi mostrando cobertura de leitura em junções éxon-éxon. Integre a análise de splicing QTL (sQTL) usando FastQTL quando dados genotípicos estiverem disponíveis, ligando variantes genéticas à variação de splicing entre indivíduos.
