Aperçu
L’épissage alternatif est un mécanisme fondamental par lequel un seul gène peut produire plusieurs isoformes d’ARNm, augmentant ainsi considérablement la capacité de codage du génome. Près de 95 % des gènes multi-exons humains subissent un épissage alternatif, générant des produits protéiques distincts qui peuvent avoir des fonctions différentes, voire opposées. Le choix des exons à inclure ou à exclure est étroitement réglementé et spécifique au tissu. L’analyse d’épissage alternatif vise à identifier et à quantifier ces différences au niveau des isoformes selon les conditions, les stades de développement ou les états pathologiques, révélant une couche supplémentaire de régulation génique au-delà des simples changements d’expression.
Méthodes
La détection d’événements d’épissage alternatifs à partir de données de séquençage d’ARN nécessite des approches informatiques spécialisées. Les outils de quantification des isoformes (tels que Salmon, Kallisto ou RSEM) estiment l’abondance de chaque isoforme de transcription connue. Les outils basés sur les événements (rMATS, MAJIQ ou SUPPA2) classent les événements d’épissage en catégories : saut d’exon de cassette, exons mutuellement exclusifs, sites d’épissage alternatifs 5’ ou 3’ et rétention d’intron. Ces outils utilisent des matrices de comptage de lectures cartographiées pour épisser les jonctions par rapport aux exons et testent les différences significatives à l’aide de modèles statistiques appariés. Détection d’épissage de novo (par StringTie ou Cufflinks) assemble des transcriptions sans annotation pour découvrir de nouvelles isoformes. La visualisation des différences d’épissage avec des tracés sashimi (de MISO ou IGV) permet d’interpréter des modèles complexes.
### Candidatures
La dérégulation de l’épissage alternatif est une caractéristique de nombreuses maladies. Dans le cancer, les facteurs d’épissage sont fréquemment mutés, produisant des isoformes spécifiques à la tumeur qui conduisent à la prolifération et aux métastases. Les troubles neurologiques tels que l’amyotrophie spinale et la démence frontotemporale résultent directement de défauts d’épissage. Comprendre les modèles d’épissage facilite l’interprétation de [la structure et les types d’ARN] (/guides/rna-structure-and-types.html) et s’appuie sur les concepts de [la transcription et le traitement de l’ARN] (/guides/transcription-and-rna-processing.html). L’analyse d’épissage recoupe également [la régulation des gènes et l’épigénétique] (/guides/gene-regulation-and-epigenetics.html), car l’état de la chromatine influence la sélection du site d’épissage cotranscriptionnel et la rétroaction d’épissage peut réguler l’élongation de la transcription.
Protocole pratique
Pour l’épissage différentiel basé sur des événements avec rMATS, commencez par les fichiers BAM de séquence d’ARN alignés par STAR à l’aide de l’indicateur « –outSAMstrandField intronMotif » pour un alignement tenant compte de l’épissage. Préparez un fichier de liste BAM et une annotation GTF, puis exécutez rmats.py --b1 sample1.bam,sample2.bam --b2 control1.bam,control2.bam --gtf annotation.gtf --od output_dir -t paired --readLength 150. rMATS détecte cinq types d’événements d’épissage : exon sauté (SE), exons mutuellement exclusifs (MXE), site d’épissage alternatif en 5’ (A5SS), site d’épissage alternatif en 3’ (A3SS) et intron retenu (RI). Les résultats incluent la différence de niveau d’inclusion (IncLevelDifference) et un taux de fausses découvertes (FDR) pour chaque événement. La métrique clé est la valeur du pourcentage d’épissage (PSI ou Ψ), qui quantifie la proportion de transcriptions incluant un exon. Le PSI va de 0 (exclu) à 1 (inclus), et un ΔΨ > 0,1 avec un FDR < 0,05 est un seuil de signification courant. Pour la quantification au niveau des isoformes sans classification d’événements, « Salmon » en mode sans alignement quantifie les abondances de transcription : « salmon quant -i transcriptome_index -l A -1 R1.fastq -2 R2.fastq -o isoform_output ». En aval, « tximport » importe des estimations au niveau de la transcription dans R, et « DRIMSeq » ou « satuRn » teste pour l’utilisation différentielle de la transcription. Visualisez les événements d’épissage avec rmats2sashimiplot, qui génère des tracés sashimi montrant la couverture de lecture sur les jonctions exon-exon. Intégrez l’analyse des QTL d’épissage (sQTL) à l’aide de « FastQTL » lorsque des données de génotype sont disponibles, reliant les variantes génétiques à la variation d’épissage entre les individus.
