多组织转录组分析揭示末端外显子使用耦合
探究mRNA亚型多样性是否由替代末端外显子使用驱动，以及第一和最后外显子使用是否相关。
使用HIT index框架分析GTEx 17,350个RNA-seq样本，量化替代第一外显子（AFE）和替代最后外显子（ALE）的使用，计算表达AFE和ALE数量的相关性以及基于位置顺序的使用相关性。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
mRNA起始与终止在空间上协调
mRNA initiation and termination are spatially coordinated.
实验完整度
高包含多组织转录组分析、长读长测序、CRISPR激活/干扰功能验证、RNAPII突变体及4sU-DRB-LRS动力学验证等多个独立证据层级。
主要模型
人类多组织RNA-seq样本（GTEx） PacBio长读长测序样本（ENCODE） HEK293T-A2细胞系 K562细胞系
重点核对
长读长测序中TSS和PAS坐标的Spearman相关系数（ρ）阈值：ρ≥0.2为PITA，ρ≤-0.2为anti-PITA CRISPR激活/干扰实验中sgRNA靶向第一外显子上游或下游400 bp区域 RNAPII延伸速率通过TT-seq与mNET-seq信号比值计算 4sU-DRB-LRS中DRB处理3小时，4sU标记10分钟 PITA富集使用ΔAUC量化，仅考虑读数≥10的基因和≥2的亚型
摘要
转录起始和终止决定驱动mRNA亚型多样性，但它们之间的关系仍知之甚少。通过系统分析转录起始位点和终止位点的联合使用，我们观察到使用上游起始位点的mRNA优先使用上游终止位点，下游位点的使用也类似耦合。我们的结果提示存在一个“位置性起始-终止轴”（PITA），其中替代末端位点的使用基于其基因组顺序而耦合。PITA在具有独特染色质特征的较长基因中富集。我们发现mRNA 5'起始选择直接影响3'末端，这取决于RNA聚合酶II的转运速度。我们的结果表明，空间组织和转录动力学将转录起始和mRNA 3'末端决定耦合起来，从而定义mRNA亚型表达。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 转录起始和终止是否受到共同调控，以及哪些机制驱动协调的RNA加工？
- 核心机制
- PITA轴基于末端外显子的位置顺序耦合TSS和PAS的使用；RNAPII延伸速率取决于起始位点，下游TSS起始的RNAPII延伸更快，倾向于使用下游更强的PAS。
- 主要证据
- 多组织短读长和长读长测序显示TSS和PAS使用基于位置顺序正相关；CRISPR激活/干扰改变AFE表达引起同 ordinal ALE表达一致变化；RNAPII慢速突变破坏耦合，快速突变增强耦合；4sU-DRB-LRS显示下游TSS起始的RNAPII延伸距离更长。
- 研究意义
- 提出PITA轴作为mRNA亚型调控的新范式，强调空间组织和转录动力学在耦合转录起始和终止决定中的作用，有助于理解细胞如何精细调控亚型表达和身份定义。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
长读长测序验证分子内TSS-PAS耦合
区分末端外显子使用是直接耦合（同一mRNA分子）还是间接耦合，并量化PITA耦合程度。
分析ENCODE PacBio ISO-seq数据（人类47个组织/细胞类型，109个样本），计算每个基因内读数起始和终止坐标的Spearman ρ，使用ΔAUC量化PITA富集。
CRISPR激活/干扰验证AFE对ALE的因果调控
直接测试改变AFE使用是否引起同ordinal ALE表达变化。
在HEK293T-A2细胞中使用dCas9-CRISPR激活或干扰扰动16个基因的27个特定AFE，通过qPCR或RNA-seq检测AFE和ALE表达变化。
RNAPII延伸速率突变体分析
测试改变RNAPII延伸速率是否影响PITA耦合。
分析过表达野生型、快速（R749H）或慢速（E1126C）RNAPII突变体的HEK293T-A2细胞的RNA-seq数据，计算AFE和ALE使用的相关性。
4sU-DRB-LRS测量TSS特异性延伸速率
估算每个TSS起始的RNAPII延伸距离，检验下游TSS是否具有更快的延伸速率。
在K562细胞中用DRB处理3小时，4sU标记10分钟，提取新生RNA，进行长读长测序（直接RNA测序和cDNA测序），计算读数从起始到终止的基因组距离。
RNAPII延伸速度计算与PITA基因特征分析
计算基因特异性RNAPII延伸速度，并分析PITA基因的延伸速度分布及启动子特征。
使用TT-seq/mNET-seq比值计算K562细胞中RNAPII延伸速度，分析PITA和非PITA基因在TSS和PAS周围的延伸速度，以及启动子预测和CpG岛分析。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
验证样本或模型特征
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | Dulbecco改良Eagle培养基 | Gibco | 11965118 |
| 胎牛血清 | Gibco | A31406–02 | |
| CRISPR转染 | 脂质体3000 | Invitrogen | L3000001 |
| CRISPR激活 | UniSAM载体 | Addgene | 99866 |
| CRISPR干扰 | dCas9载体 | Addgene | 159086 |
| RNA提取 | RNeasy迷你试剂盒 | Qiagen | 74104 |
| 逆转录 | Maxima H Minus第一链cDNA合成试剂盒 | Thermo Fisher Scientific | K1652 |
| 定量PCR | Maxima SYBR Green/ROX qPCR预混液（2X） | Thermo Fisher Scientific | K0222 |
| QuantStudio 5实时PCR系统 | Applied Biosystems | A28574 | |
| RNA-seq文库制备 | NEBNext® Poly(A) mRNA磁珠分离模块 | NEB | E7490L |
| NEBNext® Ultra™ II定向RNA文库制备试剂盒（Illumina） | NEB | E7760L | |
| PAS删除实验 | Neon转染系统 | Thermo-Fisher | MPK5000 |
| 细胞培养 | RPMI-1640培养基（含谷氨酰胺） | Gibco | 11875119 |
| 热灭活胎牛血清 | Thermo Fisher Scientific | A5256801 | |
| 青霉素/链霉素 | Gibco | 15140122 | |
| 转录抑制 | 5,6-二氯苯并咪唑1-β-D-呋喃核糖苷 | Sigma Aldrich | D1916 |
| 新生RNA标记 | 4-硫代尿苷 | ThermoFisher | J60679.MD |
| RNA提取 | Trizol | Thermo Fisher Scientific | 15596018 |
| RNA纯化 | RNA Clean & Concentrator-100柱 | Zymo Research | R1019 |
| rRNA去除 | RNAseH | Biosearch Technologies | H39500 |
| 直接RNA测序 | SQK-RNA004 | Oxford Nanopore Technologies | -- |
| cDNA测序 | SMARTer® PCR cDNA合成试剂盒 | Takara Bio | 634926 |
| Advantage 2 PCR试剂盒 | Takara Bio | 639207 | |
| SQK-LSK114 | Oxford Nanopore Technologies | -- | |
| FLO-MIN114流动池（R10） | Oxford Nanopore Technologies | FLO-MIN114 |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 长读长测序PITA鉴定 | PacBio读数过滤条件：5'端在经验第一外显子20 nt内，3'端在验证PAS 50 nt内；基因≥10条读数，亚型≥2条读数；PITA阈值ρ≥0.2 阅读提示：Materials and Methods中PacBio long read sequencing analysis和Identifying PITA genes小节 |
| CRISPR激活/干扰实验 | HEK293T-A2细胞，6孔板~80%汇合度，lipofectamine3000转染，3个生物学重复，转染后48小时收样；sgRNA靶向上游或下游400 bp 阅读提示：Materials and Methods中CRISPR-a and CRISPR-i experiments小节 |
| RNAPII延伸速率突变体分析 | HEK293T-A2细胞表达WT RNAPII、快速R749H或慢速E1126C突变体；分析恰好3个AFE和3个ALE的基因 阅读提示：Results中RNAPII elongation rates mediate the relationship between mRNA ends小节及Fig. 5A |
| 4sU-DRB-LRS | K562细胞，100 μM DRB处理3小时，最后5分钟加1 mM 4sU，洗涤后1 mM 4sU标记10分钟；直接RNA测序或cDNA测序 阅读提示：Materials and Methods中Estimation of TSS-specific elongation rates小节 |
| 延伸速度计算 | TT-seq和mNET-seq数据；延伸速度=TT-seq覆盖度/mNET-seq读数密度；排除平均归一化读数<0.05或>0.95的基因，排除TSS或PAS间距<5 kb的基因 阅读提示：Materials and Methods中Calculation of RNAPII elongation velocities小节 |