基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序（RNA-seq）数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件（其中新识别的可变剪接事件3922个）,包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.     

Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer

Background

The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)_nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.

Results

Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.

Conclusion

We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1143) contains supplementary material, which is available to authorized users.     

可变剪切在植物成花转换中的作用

精确调控成花转换, 确保植物在适宜环境下开花, 对于植物的成功繁殖和物种繁衍至关重要。开花由多种分子机制在转录、转录后和蛋白质水平进行调控。可变剪切(AS)是一种普遍的转录后水平调控过程, 可从单个基因产生多个转录本, 从而丰富转录组和蛋白质组的多样性。大量研究表明, 可变剪切在成花转换过程中发挥重要作用。根据发育和环境条件, AS能够影响mRNA的稳定性和/或蛋白亚型的功能, 从而调控开花相关基因的功能转录本和/或功能蛋白水平。揭示成花相关pre-mRNA的AS作用将进一步增进人们对开花相关基因功能以及整个成花转换调控网络的认识。该文归纳了涉及成花转换的AS研究进展, 并针对各个调控途径进行总结, 以期为进一步研究植物AS和成花转换调控机制提供参考。     

中华猕猴桃基因组可变剪接事件鉴定及分析

可变剪接使一个基因能产生多种m RNA成熟体,极大地增加蛋白多样性.采用中华猕猴桃基因组数据做参考数据,利用中华猕猴桃叶片和果实3个不同发育时期(未成熟、半成熟和成熟期)的转录组数据,从中华猕猴桃基因组(39040个基因)中共鉴定出11651个基因(占总基因数的29%)对应的32180个可变剪接事件.在可变剪接不同类型中,内含子保留类型的发生频率最高,占50%以上;3′可变位点类型频率约为5′端可变类型的2倍.GO富集分析结果表明,可变剪接的基因主要富集于酶调控及核苷酸结合相关功能的GO类别中,而组织特有可变剪接基因功能富集热点与组织的重要功能关联,叶片多为肌动蛋白及微管相关;未成熟果实与双组分信号系统相关;半成熟果实多与磷脂合成过程相关;成熟果实多与信号传递过程相关.另外,55.6%的维生素合成相关基因发生可变剪接事件,显著高于基因组水平的29.6%,暗示着可变剪接参与维生素合成相关基因代谢过程中的重要作用.通过对中华猕猴桃全基因组可变剪接的分析,为解析中华猕猴桃基因组及进一步开展相关分子育种工作提供依据.