跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
分子细胞生物学。2008年9月;28(17): 5147–5161.
2008年6月23日在线发布。 数字对象标识:10.1128/MCB.00331-08号
预防性维修识别码:项目经理2519722
PMID:18573883

推测MicroRNA对S1结肠癌细胞株转录稳定性和蛋白翻译的调控及其mRNA 3′非翻译区ABCG2表达

关联数据

补充资料

摘要

ABCG2在临床药理学中被认为是一种重要的外排转运体,在化疗药物耐药中具有潜在的重要作用。为了确定调控ABCG2 mRNA在其3′非翻译区(3′UTR)表达的表观遗传机制,我们对S1亲代结肠癌细胞系及其耐药的ABCG2过度表达对应物进行了cDNA末端的3′快速扩增。我们发现亲代细胞中3′UTR的长度大于1500 bp,并且使用miRBase TARGETs数据库,在耐药细胞中ABCG2 mRNA缺失的3′UTR3′部分中确定了一个推测的microRNA(miRNA)结合位点,与最近报道的hsa-miR520h位点不同。我们假设ABCG2的3′UTR处假定的miRNA结合抑制ABCG2表达。在耐药S1MI80细胞中,由于3′UTR较短,miRNA不能与ABCG2 mRNA结合,因此mRNA降解和/或对蛋白翻译的抑制得到缓解,导致ABCG2过度表达。通过报告基因分析、miRNA结合位点的突变分析以及miRNA抑制剂或模拟物的强制表达,严格验证了这一假设。miRNA去除这种表观遗传调控可能与耐药癌细胞中ABCG2的过度表达有关。

ABCG2是一种普遍存在的ATP-结合盒(ABC)转运体,除了在底物药物的吸收、分布和消除中发挥重要作用外,还可能导致癌细胞产生多药耐药性(1).ABCG2公司在使用各种抗癌药物筛选的人类肿瘤细胞系中经常观察到过度表达(12,34,38,44,52).

调节的分子机制ABCG2公司表达不太清楚。与大多数TATA-less基因启动子一样ABCG2公司启动子包含许多Sp1、AP1和AP2位点和CCAAT盒。ABCG2公司也有一个假定的CpG岛位于基因上游(50). 迄今为止,大多数研究都在研究ABCG2公司集中在转录水平。两个功能顺式中的元素ABCG2公司启动子,即雌激素(15)和缺氧(28)反应元件和过氧化物酶体增殖物激活的受体反应元件ABCG2公司基因(48)已报告。已提出一种芳香烃受体反应元件,但其确切序列尚未确定(13). 在耐药MCF-7细胞中,由于ABCG2 mRNA 5′非翻译区(5′UTR)剪接变异体的差异表达,已观察到启动子的替代使用(39). 我们最近报道了DNA甲基化和组蛋白修饰在调节人类ABCG2公司在肾癌细胞系中(50)在一些药物选择耐药细胞系中(49)分别是。

ABCG2的3′UTR知之甚少。在真核mRNAs中,3′UTR通过调节核质mRNA转运、多聚腺苷酸化状态、亚细胞靶向性、翻译效率、稳定性和降解速率,在转录后水平上调节基因表达中发挥重要作用(11,40,45,47,54). 在GenBank(登录号:。NM_004827号)在UTR数据库(UTR数据库条目UTR:3HSA117529)中,长度约为2 kb,比人类mRNA平均770 bp长得多(24,30)这表明它在基因表达调控中可能具有一个或多个重要作用。

MicroRNAs(miRNAs)是控制动植物基本细胞过程的一大类基因调节分子(2,16,21,36). miRNAs的基因调控通常是通过与靶mRNAs 3′UTR序列形成不完美的杂种来介导的,从而诱导翻译抑制和/或mRNA降解(10,29,42).

在本研究中,我们证明ABCG2 mRNA在亲代S1结肠癌细胞系中采用的3′UTR比其耐药对应物更长,并且miRNA(hsa-miR-519c)通过推测的hsa-miR-519c结合位点降低内源性ABCG2 mRNA和蛋白水平,该结合位点位于仅在亲代细胞中发现的较长的3′UTR区域内。这些发现表明,逃避miRNA介导的翻译抑制和mRNA降解可能导致耐药癌细胞中ABCG2过度表达。

材料和方法

组织培养。

人类结肠癌细胞株S1及其耐药亚系S1MI80已在前面描述过(38). 本研究从NCI肿瘤药物筛选中选择SW620、MCF-7、H460和SF295细胞。细胞系保存在Iscove改良的Eagle培养基(S1、S1MI80和MCF-7细胞)或RPMI培养基(SW620、SF295和H460细胞)中,补充有10%胎牛血清、100单位/ml硫酸链霉素和100单位/ml硫酸青霉素G,并在37°C和5%CO中培养2.

RNA降解分析。

将父母S1和耐药S1MI80细胞在10 cm组织培养皿中培养至亚融合状态,然后将放线菌素D(西格玛,圣路易斯,密苏里州)添加至最终浓度5μg/ml,以阻止从头开始的RNA合成。放线菌素D处理后4、8和16小时,收集细胞,并用如下所述的半定量逆转录-PCR(RT-PCR)定量ABCG2 mRNA。作为对照,还监测c-myc和甘油醛-3-磷酸(GAPDH)mRNA水平。对于转染pcDNA3.1载体的HEK293细胞,在转染48小时后添加放线菌素D,然后在相同的时间点采集细胞。为了标准化目的,评估共转染的绿色荧光蛋白(GFP)(pEGFP-C1;BD Bioscience Clontech)的表达水平。记录的值是与添加放线菌素D之前的量相比,mRNA的剩余百分比。所有实验重复三次。

半定量RT-PCR。

使用Trizol试剂(Invitrogen,Carlsbad,CA)分离总RNA。使用商用cDNA合成试剂盒(马萨诸塞州伦道夫Bioline)逆转录RNA(1μg)。cDNA扩增是使用针对ABCG2公司(5′-CAATGGGATCATAAACCTG-3′[向前]和5′-GAGGCTGATAGAA-3′[向后])和GAPDH公司(5′ACCCAGTCACATGCACATCATCAC-3′[正向]和5′TCCACCACCTGTGCTGTA-3′[反向])。引物的设计是这样的,它们的内含子跨度为3到6英寸ABCG2公司和内含子8GAPDH公司GAPDH cDNA的扩增作为内部对照。PCR扩增在55°C的退火温度下进行28个周期(ABCG2)或25个周期(GAPDH),分别产生584-或330-bp产物。当cDNA合成混合物中省略逆转录酶时,未获得PCR产物。PCR产物在2%琼脂糖凝胶中溶解并用溴化乙锭染色。采集凝胶图像,并使用APP_Collage PPC4.0分析软件程序量化条带强度。

实时PCR。

对于mRNA降解试验,根据制造商的建议,使用LightCycler Fast Start DNA Master Taqman试剂盒(Universal ProbeLibrary;印第安纳波利斯罗氏分子生物化学公司)和LightCyscler 2.0系统,通过实时定量PCR方法重复测量ABCG2 mRNA水平。使用Roche UPL定制设计服务为ABCG2和GAPDH设计引物(用于归一化),并选择合适的荧光探针(分别为ABCG1和GAPDH的探针56和60)。用于分析ABCG2和GAPDH表达的引物序列如下:ABCG2正向,5′-TGGCTTAGACTCAAGCACAGC-3′;ABCG2反面,5′-TCGTCCTGCTTAGACATCC-3′;GAPDH正向,5′-AGCCACATCGCTCAGACC-3′;和GAPDH反向,5′-GCCAATACACACAAATCC-3′。相对ABCG2和GADPH表达水平通过读取MCF-7细胞cDNA五倍四倍连续稀释产生的标准曲线来确定,该曲线绘制了阈值周期与模板的对数输入量。

使用LightCycler Fast Start DNA Master Plus Sybr green I试剂盒(Roche Applied Science),通过定量实时RT-PCR重复分析ABCG2 3′UTR的相对稳定性。使用LightCycler系统在90°C下进行35个周期的定量PCR,持续10 s,引物特定退火温度持续10 s和72°C持续18 s。每个周期结束后,将温度升高至引物各自的采集温度,并在530 nm(Light Cycler荧光通道F1)处测量与双链DNA结合的Sybr绿的荧光。样品荧光超过背景荧光的循环数由LightCycler软件(3.5版)使用二阶导数方法确定。扩增阶段后进行熔融曲线分析,以消除非特异性扩增或引物二聚体形成的可能性。

3′RACE分析。

使用Trizol试剂(Invitrogen)从亲本S1和耐药S1MI80细胞中提取总RNA。根据制造商的协议,使用3′快速扩增cDNA末端(3′RACE)cDNA扩增试剂盒(Bioline)进行第一链cDNA合成。3′RACE使用适配器底漆(AP)(表(表1;1; 参见补充材料中的图S3A),一个寡核苷酸(dT)引物,在其5′端具有定义的(非dT)序列,用于启动第一链cDNA合成。AP的寡核苷酸(dT)部分与mRNA poly(A)尾部杂交以引发RT反应。因此,AP 5′端的确定序列被纳入所有合成的cDNA分子中。然后用Bio-X-ACT长DNA聚合酶(Bioline)对第一链cDNA进行PCR。外部反向底漆(表(表1;1; 参见补充材料中的图S3A),然后可以与与已知ABCG2序列结合的上游引物结合使用,以从上游引物特异性扩增到mRNA多腺苷酸化位点。底漆ABCG2 ex16(表(表1;1; 参见补充材料中的图S3A)被用作我们研究中报告的3′RACE反应的上游引物。为了提高第一次PCR的特异性,使用嵌套的N1正向引物、外反向引物和第一次PCR混合物作为DNA模板进行嵌套PCR(PCR 1)。将PCR中的DNA作为混合物亚克隆到pCR2.1载体(Invitrogen)中。然后由NCI实验致癌实验室的DNA测序核心设施对多个克隆进行测序。所有序列与报告的mRNA序列进行比对ABCG2公司(GenBank登录号。纳米_004827). 为了证实具有长3′UTR的克隆的RACE发现,使用嵌套N1正向引物进行嵌套PCR 2(ABCG2公司核苷酸[nt]2401至2429)和anABCG2公司-将特定的嵌套N2反向引物退火到报告的mRNA序列的3′端附近(nt 4240至4263)(表(表1;1; 参见补充材料中的图S3A)。对于S1细胞,使用hsa-miR519c-特异性抑制剂(Dharmacon)转染后从细胞中获取的总RNA重复3′RACE分析。

表1。

本研究中使用的寡核苷酸引物

使用和底漆序列(5′-3′)(nt)
3英尺赛道
    ABCG2 ex16AAAGCAGGGCATCGATCTCTCAC(2338-2360)
    外部反向底漆CAGTCGTCTCTCGCAGGGTTCAGGCAGGCACTCTGG
    嵌套N1正向底漆GATTGTATTTCCTCACAATTGCCTACC(2401-2429)
    嵌套N2反向底漆TATACCAAGAACAGAAAAACTGA(4263-4240)
    3分赛道APCCCTGTTCAAGCGCATCTGAGGTAACCATGACCGTGC(T)18
3′UTR RT-PCR
    A-Fw(A-Fw)TTCCCCTTAATTCAGTATGATTTCACA(2464-3030)
    A-R公司TACCACTAGACAGGAATATCAGTAAGTTC公司
    B-Fw(B-Fw)GAGAATTGAATTCTGGAACTCCTGA(2731-3291)
    B-R公司在加拿大
    C-FW公司GAACTTACTGATTATTCCTGTCTAGTGGTA(3002-3558)
    C-R公司CTGAGTAGGCTTTTGTGTGCTACAATCTCTAT公司
    D-Fw(D-Fw)TCATAAGGAATCTGACCACTCTGGAT(3264-3812)
    德黑兰CATTAAGGGAAGCTGAGAGAGTAGAGACAGAGAGAGACATTTT公司
    电子Fwatagaagatttagcacacaaaagctactcag(3526-4024)
    电子-无线电CTCAGTTACTCCTGTAAGTGCTGA公司
    F-Fw(飞行-飞行)AAAGTTCCTGTCTACTCTTGCTCAGCTTCCTTAATG(3775-4388)
    F-R公司GTGGCTGTTTTTCTTTTGAAGTATATTTTGTG公司
荧光素酶报告子构建物的诱变制备
    3′UTR(3820/3841)
        miR519c站点FCGGAAAGATCGCCGTGTGTATTTATATAGAGATGTAGTCGGGGCGGC
        miR519c站点RGCCGCCCCGACTCTAGACAAATGCTATCTCTATAAAAATTACGCGCTACTTTCCG
    3′UTR加扰(3820/3841)
        扰码(3820-3841)FCGGAAAGATCGCGTAATTACTGTCGTAAAGTTTAAGATTCTAGAGTCGGGGCGGC
        加扰(3820-3841)RGCCGCCCCGACTCTAGAAATCTTAAACTTTAACGACAGTAATTAACGCGCGCTCTTCCG
    3′UTR FLΔ(3820/3841)
        删除(Δ3820-3841)FCTACTCTTGCTCTCTCTTCCCTTAATGTTAACATTCAAGCTAACATTAAGTTAACATTGGTAC
        删除(Δ3820-3841)Rgtaccaatgtaacttgttagctttgaatgtaacattaagggaagctgaggagagaggtag

生物信息学。

使用miRBase TARGETS对hsa-miR-519c的miRNA靶点进行预测(19)和RNAHybrid(43)软件。使用UTRScan程序分析ABCG2 3′UTR(http://bighost.area.ba.cnr.it/BIG/UTR主页/)搜索功能元素的存在(41).

RNA二级结构分析。

使用mfold程序进行RNA二级结构分析(56)(可在http://www.bioinfo.rpi.edu/applications/fold网站). 折叠温度固定在37°C。使用miR-519c结合序列(带有或不带有相邻核苷酸±100 bp)进行分析。序列用mfold以本地自动方式折叠。最稳定的热力学结构如补充材料中的图S8所示。

用miRNA抑制剂或模拟物转染细胞。

从Dharmacon购买了一种miRNA抑制剂和模拟物。使用寡效胺(Invitrogen)将高达120 nM的特异性hsa-miR519c抑制剂或总体积为10 ml的模拟物转染细胞。模拟物是包含成熟miRNA序列的双链RNA。当模拟物被转染到细胞中时,其结果类似于添加一个与内源性miRNA功能相似的外源性加工miRNA。抑制剂包含一个与成熟miRNA序列互补的序列,并充当诱饵。内源性miRNA与抑制剂紧密结合并被隔离,因此内源性靶点的结合较少。如下所述,通过干环RT-PCR评估miRNAs的表达。

miRNAs的干环RT-PCR。

如前所述,对成熟的hsa-miR-519c和hsa-miR-520h进行干环RT(6,46). 所有用于miRNA干环RT的试剂均来自Applied Biosystems(加利福尼亚州福斯特市),使用hsa-miR-519c(P/N:4373251)或hsa-miR 520h(P/N:4373258)特异性干环RT引物。组装RT反应混合物,并在25℃下培养10 min,在42℃下培养60 min,在70℃下培养5 min,然后在4℃下保持。PCR产物用以下引物扩增:miR-519c正向,5′-GCGCAAAGTGCTTTTAGGAT-3′;miR-519c反面,5′-GTGCAGTCGAGGT-3′;miR-520h正向,5′-GCGGCAAAGTGCTTCCCTTTAGAGT-3′;miR-520h反向,5′-GTGCAGGGTCCGAGGT-3′;U6正向,5′-CTCGCTCGGCAGCACA-3′;和U6反面,5′-AACGCTCACGAATTTGCGT-3′。PCR产物在2%琼脂糖凝胶中进行分析。U6小核糖核酸被用作内部对照。

荧光素酶报告分析。

荧光素酶报告基因构建物是使用携带猴病毒40(SV40)启动子荧光素酶表达单元的pGL3-控制载体(Promega,Madison,WI)制备的。5′UTR片段ABCG2公司(nt−1289 to+396)通过PCR从先前描述的启动子结构中扩增(50),并通过RT-PCR从H460细胞总RNA中产生一个3′UTR(nt+2462-+4430)片段。将其插入萤火虫萤光素酶基因的下游,位于pGL3对照载体中萤光素酶编码序列和poly(a)信号之间的XbaI位点。通过限制性映射和序列分析验证了结构的正确定位。根据全长(FL)3′UTR结构(碱基+2462到+4430),我们构建了一系列3′缺失结构,3′端终止于碱基+4060、+3028和+2761。碱基的编号是相对于GenBank登录号中定义的核苷酸1指定的。NM_004827号另一个结构体(3′UTR FLΔ3820/3841)是通过从FL结构体中删除假定的hsa-miR519c结合位点而制备的。在荧光素酶基因下游克隆了一个与假定的hsa-miR519c结合位点相对应的短序列,得到了3′UTR(nt3820-3841)结构。作为对照,以类似的方式制备带有从nt3820到3841的加扰序列的向量(见图。图4A;4A级; 也可参见补充材料中的图S9)。

保存图片、插图等的外部文件。对象名称为zmb017087649004.jpg

ABCG2荧光素酶报告分析。(A) 显示用于研究ABCG2 mRNA的3′UTR的荧光素酶报告子结构的示意图。亲本载体pGL3-control包含萤火虫荧光素酶的编码区,其表达由SV40启动子(Promega)驱动。核苷酸“1”是指ABCG2转录本中的第一个核苷酸,如序列NM_004827所定义。生成了以下结构:3′UTR FL(2462/4430);三个截断的3′UTR结构(基于FL 3′UTR-结构,从3′UTR3′端进行渐进删除,其3′端标记在FL结构的顶部,即nt 2761、3028和4060);3′UTR(3820/3841),仅包含hsa-miR-519c的22-nt假定miRNA结合位点;3′UTR FL(Δ3820-3841),包含缺乏hsa-miR-519c结合位点的FL 3′UTR;3′UTR(scramb 3820/3841),包含相同的22nt miRNA结合位点,但顺序随机;和5′UTR构建物(−1289/+396),通过插入荧光素酶编码序列下游ABCG2 5′UTRs的片段(−128至+396)制备(用作对照,以证明3′UTR对报告分析的荧光素瘤酶活性的特定影响)。图表未按比例绘制。补充材料中的图S9中提供了更详细的报告者结构图。(B) 如图所示,各种报告结构的荧光素酶报告活动。平均报告活动±SD(萤火虫/雷尼利亚荧光素酶单位)。pGL3-control表示没有ABCG2 3′UTR的矢量主干。(C) 与miRNA抑制剂(左)或模拟物(右)共同转染的细胞中FL ABCG2 3′UTR FL(2462/4430)构建物的荧光素酶报告活性。根据B组所述,获得了平均报告活动±SD。学生进行了统计分析t吨测试。

按照制造商的方案,使用Fugene 6转染试剂(Roche)将萤火虫荧光素酶构建物(200 ng)转染到12孔板中的S1细胞中。为了监测转染效率,用50 ng phRG-Basic质粒(Promega)共同转染细胞,编码雷尼利亚荧光素酶。转染48小时后,使用双荧光素酶报告分析系统(Promega)测量荧光。所有转染均重复进行,通过将萤火虫荧光素酶活性标准化为雷尼利亚每个样本的荧光素酶活性。每个构建体在三个独立的转染中进行测试。为了测试miRNA抑制剂或模拟物对报告构建物的影响,使用DharmacFECT Duo转染试剂(Dharmacon)进行转染。结果表示为三个独立实验的平均值±标准偏差(SD)。数据集由未配对的学生进行比较t吨测试,显著差异被认为是P(P)值<0.05。

蛋白质印迹分析。

用十二烷基硫酸钠-聚丙烯酰胺凝胶电泳分离转染HEK293或A549细胞制备的全细胞裂解物,并用ABCG2(Kamiya Biomedical Company,Seattle,WA)、GAPDH(American Research Products,Belmont,MA)、,和GFP(细胞信号技术,丹弗斯,马萨诸塞州)。用1:10000稀释的山羊抗鼠(IRDye800CW)或抗兔(IRDy680)二级抗体(Li-Cor)孵育后,用Odyssey(Li-Cor,Licolin,NE)红外成像系统分析印迹。

结果

亲本S1和耐药S1MI80细胞系中ABCG2 mRNA的稳定性。

用放线菌素D抑制RNA合成后,测定了亲本S1结肠癌细胞系及其耐药亚系S1MI80中ABCG2 mRNA的半衰期。放线菌素D处理后不同时间间隔分离总RNA,并通过RT-PCR分析ABCG2 mRNA水平。结果显示为用GAPDH标准化后ABCG2 mRNA的剩余百分比(图。(图1)。1). 放线菌素D处理后,S1MI80细胞中的ABCG2 mRNA在16小时以上相当稳定,而亲代S1细胞中的ABCG2 mRNA的半衰期约为6小时。亲代S1细胞中romidepsin上调ABCG2(49)没有改变其mRNA的稳定性(图。(图1)。1). 为了确保放线菌素D的转录阻断在我们的实验中有效,我们进行了RT-PCR,以检测c-myc和GAPDH mRNA,分别作为短半衰期(15至40分钟)和长半衰期的对照(超过24小时)(参见补充材料中的图S1)。重要的是,GAPDH和c-myc mRNA在未经处理和经romidepsin处理的S1和S1MI80细胞中的降解情况相似(参见补充材料中的图S1)。通过定量实时PCR(通用探针库基因检测;印第安纳波利斯罗氏应用科学公司)(参见补充材料中的图S2A)对mRNA稳定性研究中的ABCG2和GAPDH水平进行重复分析,结果与上述结果相当。我们推测ABCG2 mRNA降解速率的降低可能有助于耐药S1MI80细胞中ABCG2的上调。在其他几对亲本和ABCG2过度表达的耐药细胞系中也观察到不同的ABCG2 mRNA稳定性(数据未显示)。

保存图片、插图等的外部文件。对象名称为zmb017087649001.jpg

(Top)ABCG2 mRNA在耐药S1MI80细胞中比在亲本S1细胞中更稳定。虽然ABCG2也被罗美肽治疗上调(2 ng/ml,24 h),但其mRNA稳定性不受影响。放线菌素D(5μg/ml)抑制了新生RNA的合成,并在处理后0、4、8和16小时采集RNA。对ABCG2 mRNA进行RT-PCR分析,以追踪ABCG2mRNA随时间的剩余量。c-myc和GAPDH mRNA水平也分别作为快速降解和稳定mRNA的对照进行监测。记录的值是与添加放线菌素D前、用GAPDH标准化后的mRNA残留量相比的mRNA百分比。所示数据代表三个独立实验的平均值±SD。(底部)代表性凝胶图像,显示用放线菌素D预处理的细胞中ABCG2、c-myc和GAPDH的RT-PCR分析。

通过3′RACE分析鉴定人ABCG2 mRNA的3′端。

进行3′RACE反应以检测ABCG2 mRNA的3′端并确定poly(A)位点的位置(17). 使用AP与从亲本S1和耐药S1MI80细胞分离的模板mRNA杂交,进行第一链cDNA合成。然后对第一链cDNA进行PCR(第一PCR),使用与AP中唯一序列互补的外反向引物,结合ABCG2特异性正向引物ABCG2 ex16,在终止密码子上游的ABCG2外显子16处与nt 2338至2360互补(参见补充材料中的图S3A)。该PCR为抗性S1MI80亚系(第一个PCR)产生了多个小产物(200至400 bp)(参见补充材料中的图S3B)。然而,无法从亲本S1细胞中获得PCR产物,可能是因为该细胞系中ABCG2的表达水平极低。未经RT预处理的样品未发生显著反应(数据未显示)。为了提高第一次PCR的特异性,使用ABCG2特异的巢式N1正向引物(在ABCG2的第16外显子nt 2401至2429的补充,仍在终止密码子上游)和外反向引物(数据未显示)进行巢式PCR 1。嵌套PCR 1产生了更清晰的PCR条带,但基本上与第一个PCR中的模式相同。由于从我们的3′RACE分析中获得的ABGC2 3′UTR的长度(S1MI80细胞约200至400 bp长)比GenBank和UTR数据库中报告的长度短得多,因此我们评估了其他敏感的亲本细胞系(SW620、MCF7、SF295和H460),以确定这是否是一种细胞类型特异性观察。令我们惊讶的是,3′RACE分析也为其他亲本细胞系产生了长产物(~2 kb)和多个小产物(200至400 bp)(参见补充材料中的图S3B)(可能是因为该细胞系中ABCG2表达水平低,所以没有从SW620亲本细胞获得PCR产物)。

值得注意的是,在含有多个模板的PCR中,引物结合位点之间的序列长度不同(即来自不同长度的3′UTR的第一链cDNA),较短的模板通常优先扩增(4). 这可能导致通过3′RACE分析测定3′UTR长度时出现偏差,因为它低估了长3′UTR-PCR产物。为了克服这个问题并排除ABCG2 3′UTR中的选择性剪接,我们使用来自亲本和抗性细胞系对的总RNA以及与ABCG2 3′UTR cDNA中不同区域互补的引物进行了RT-PCR实验(图。(图2A;2安培; 表1)。1). 六个独立重叠的PCR片段(A至F),覆盖了报告的ABCG2 3′UTR序列(NM_004827号)被放大(图。(图2A)。2安培). 针对所有检测的RT样本,获得了所有六个PCR片段,其大小符合预测值,消除了ABCG2 3′UTR处选择性剪接的可能性。当cDNA合成混合物中省略逆转录酶时,未获得PCR产物。使用来自S1细胞的一系列输入基因组DNA构建校准曲线(参见补充材料中的图S4)。在用GAPDH归一化后,通过读取校准曲线来估计不同3′UTR片段与最上游区域A相比的相对丰度。对于亲本S1细胞,发现6个3′UTR区域的相对丰度相似(图。(图2B)。2B型). 然而,对于耐药S1MI80细胞,从A区到F区,下游更多的3′UTR区域的相对丰度越来越低,这表明耐药细胞中截断的3′UTR的比例越来越大(图。(图2B2B型).

保存图片、插图等的外部文件。对象名为zmb017087649002.jpg

RT-PCR分析研究亲本S1和耐药S1MI80细胞中ABCG2 3′UTR不同片段的相对丰度。(A) 用RT-PCR分析ABCG2 3′UTR内六个不同重叠片段的大致位置示意图。每个片段约500 bp长。(B) 用GAPDH归一化后,绘制了不同3′UTR片段相对丰度与A区数量的关系。开条,亲代S1细胞;封闭条,抗S1MI80细胞。所示数据代表三个独立实验的平均值±SD。由于hsa-miR-519c的结合位点位于ABCG2 3′UTR内的nt 3820至3841,因此只有可由hsa-miR 519c调节的长ABCG2转录物才能产生片段F的PCR产物(nt 3775至4388)。然而,长和短ABCG2转录物都为片段A(nt 2462至3030)产生PCR产物。通过比较片段A和F的PCR产物,我们估计亲本S1细胞中约77%的ABCG2转录物可以被hsa-miR-519c调节。(C) 六种不同3′UTR片段的相对mRNA稳定性。在放线菌素D(5μg/ml)处理后的不同时间间隔分离总RNA,并通过RT-PCR分析A组6个3′UTR区域的mRNA水平。代表性实验的数据绘制为每个3′UTR片段的mRNA剩余百分比,与添加放线菌素D之前的量相比,使用GAPDH标准化后的量。

接下来,使用放线菌素D处理后不同时间间隔从S1和S1MI80细胞中分离的总RNA,通过RT-PCR评估六个不同3′UTR区域的相对mRNA稳定性(图。(图2C;2摄氏度; 参见补充材料中的图S5)。虽然在亲本S1细胞中发现所有六个区域的降解速度相似,但在抗性S1MI80细胞中,上游区域(即A到C)比下游区域(即D到F)更稳定(图。(图2C;2摄氏度; 参见补充材料中的图S5)。使用SYBR绿色检测(Roche Applied Sciences)用定量实时PCR重复此3′UTR PCR分析,并获得了可比较的结果(参见补充材料中的图S6)。结合3′RACE结果,这表明耐药S1MI80细胞中存在较短的3′UTR转录物,并且它们比亲本S1细胞中较长的转录物更稳定。

不同长度ABCG2 3′UTR序列分析。

首先使用亲本H460细胞的3′RACE产物,然后使用亲本S1细胞和S1MI80细胞的产物,我们将3′RACE-PCR混合物作为一个整体(第一个PCR)克隆到TOPO载体中,并对得到的克隆进行DNA测序。图3A3A级说明了在GenBank中报道的ABCG2 mRNA转录物的3′UTR内的一些推定功能元件的相对位置(登录号:。NM_004827号)[包括多聚(A)信号、假定的富含AU的元素(ARE)和几个假定的miRNAs结合位点(见下文)]。从我们的测序结果中,发现了五个不同长度的ABCG2 3′UTR转录本。与3′RACE分析一致,虽然在亲代细胞中有一个长的(~2 kb;在nt 4419处结束)和四个短得多的(100至400 bp;在nt 2471、2525、2644和3013处结束)3′UTR转录物,但在耐药细胞中只发现了四个短的转录物。3′UTR处核苷酸的编号与ABCG2转录本中按顺序指定的第一个核苷酸有关NM_004827号,在整个报告中使用。Clustal X 2.0多序列比对软件的序列比对(ftp://ftp.ebi.ac.uk/pub/software/clustalw2)表示五个3′UTR序列非常相似,除了3′端的长度不同(参见补充材料中的图S7)。在我们的研究中,还报道了ABCG2(GenBank登录号:。BC092408号; UTR数据库条目UTR:3HSA117529)。该观察结果表明,不同长度是3′UTR的多聚腺苷化的结果,而不是选择性剪接。

保存图片、插图等的外部文件。对象名为zmb017087649003.jpg

(A) 报道的ABCG2 mRNA 3′UTR的GenBank序列(NM_004827)中的假定功能元件。终止密码子(TAA)、规范(AAUAAA)和非规范(AUUAAA)多(A)信号、三个假定ARE(UAAUUAUA、AUUAUUAAUA和GAGAAUUAAUU)的位置,标记了ABCG2的3′UTR内推测的hsa-miR-519c和miR-520h结合位点。(B) miRBase TARGETS第4版预测的四组假定miRNA靶点(阴影区域)示意图(19),在ABCG2的3′UTR范围内。关于这些推测的miRNAs及其在ABCG2 3′UTR上的结合位点序列的详细信息总结在表中表2。2.最上游的簇可以被13种不同的人类miRNAs(hsa-miR-518c、hsa-miR-519d、hsa-micR-302a、hsa-milR-302d和hsa-miR-520a-h)结合,第二个位点可以被hsa-miR142-5p结合,第三个位点可以由hsa-miR.100结合,最下游的位点可以被hsa-miR-512c结合。RNA杂交自由能(ΔG公司)还显示了hsa-miR-519c的计算和理论miRNA-mRNA双链配对。此外,还标记了hsa-miR-519c序列中的“种子”区域,该区域决定了假定miRNA的特异性,以及可能破坏miRNA-mRNA相互作用的假定结合位点中G-U摆动的位置。在hsa-miR-519c序列5′端第4位发现的G-U摆动可能对miRNA靶点识别有害,但已经证明这种G-U碱基配对通常在种子区可以耐受(9).

仔细观察补充材料中图S7所示的序列,发现在五个备选聚(A)端之前有两个聚(A。显然,位于nt 2626的poly(A)信号(AUUAAA;一个非经典信号),位于终止密码子TAA(nt 2459)下游,产生较短的3′UTR,而另一个位于nt 4405的poly。(图3A;3A级; 参见补充材料中的图S7)。此外,除了3′UTR较短的转录本(终止于nt 2471、2525和2644)外,ABGC2的所有其他转录本都有一个或多个假定ARE或其变体的副本(UAAUUUAUA、AUUAUUAAUA和GAGAAUUAAUU,分别发现于nt 2710、2857和4412)。AREs的作用是降解靶mRNA,因此有无AREs存在不同长度的ABCG2 3′UTR可能会区分其mRNA稳定性。

ABCG2 3′UTR包含四个推测的miRNA结合位点簇。

为了寻找可能参与ABCG2 mRNA稳定性的miRNA结合位点,使用miRBase TARGETS数据库进行了电子分析(网址:http://microrna.sanger.ac.uk) (19)和RNAHybrid(43)在ABCG2 3′UTR上寻找推测的miRNA结合位点。RNAHybrid程序使用自由能(ΔG公司)计算miRNA与相关mRNA内潜在位点之间有利结合相互作用的算法。这些分析确定了四个假定的miRNA结合位点(nt 2491至2512、2823至2841、3052至3073和3820至3841)(图。(图3B)3B公司)和几种潜在的miRNA(总结见表表2)。2). 在ABCG2 3′UTR上的四个假定的miRNA结合位点中,hsa-miR-519c特异性的结合序列位于较长的3′UTR的远端3′端,并且是抗性细胞中较短的3′UTR中唯一缺失的一个(位置3820)。请注意,这个位点只在人类身上发现。

表2。

ABCG2 3′UTR上推测的miRNAs及其结合位点序列

职位3′UTR的miRNA结合位点序列推测miRNAΔG公司(千卡/摩尔)
2492-2513CACAUAAAAAGAGCACUU公司hsa-miR-302a−13.24
UCACAUAAAAAGAGACACUUU公司hsa-miR-302d−18.53
2493-2514CACAUAAAAAGAGCACUUG公司hsa-miR-518c−17.27
2492-2514ACAUAAAAAGAGCACUUUG公司hsa-miR-519d−16.28
2492-2513ACAUAAAAAGAGCACUU公司hsa-miR-520a−15.10
2493-2513UCACAUAAAAAGAGACACUUU公司hsa-miR-520b−15.26
2491-2513年中非共和国hsa-miR-520c−16.59
CACAUAAAAAAAGAGCACUUU公司hsa-miR-520d−23.26
2493-2513ACAUAAAAAGAAGCACUUhsa-miR-520e基因−17.96
2491-2512中卡卡阿亚阿加加卡库hsa-miR-520f−15.56
2492-2515ACAUAAAAAGAGACACUUGA公司hsa-miR-520g−20.43
2494-2515AUAAAAAAGAAGACUUGA公司hsa-miR-520小时−17.43
2493-2513CACAUAAAAAAAGAGCACUUU公司hsa-miR-526b−15.43
2824-2842UUCAGUCUUAAUACUUAUG公司hsa-miR-142-5p−15.43
2953-2974GAGGUUGGAACUGGUGGU公司hsa-miR-100−24.71
3820-3842AGUAAAUACUUUUCAGCACUUAC公司hsa-miR-519c−14.48磅
核苷酸“1”是指ABCG2转录本中的第一个核苷酸,如NM_004827序列中所定义。

我们询问hsa-miR-519c是否参与ABCG2抑制。miRNA阻遏的一个重要决定因素是miRNA-mRNA双链在“种子”序列中的稳定性,即miRNA 5′端的nt2到8(5,10,20,23). 所谓的“7-mer-m8位点”种子匹配(包含与miRNA核苷酸8匹配的种子匹配)(20)预测为hsa-miR-519c(图。(图3B),3B公司)具有良好的杂交自由能(−14.4 kcal/mol)。

长期以来,靶点可及性被认为是有效反义寡核苷酸和siRNA-介导沉默的重要因素(31). 因此,我们还分析了预测的miRNA结合位点序列的稳定性,以确定我们是否可以找到可能影响可及性的RNA结构特征,并增强靶点预测的特异性(55). 因此,这些稳定或不稳定元件的存在可能有助于支持对假定的miRNA靶点的计算机筛选。使用mFold(56)根据自由能预测(ΔG公司)(见上文)和二级RNA结构。预测hsa-miR-519c结合位点序列的二级结构有一个不稳定的发夹环(补充材料中的图S8显示了mFold单独或在±100-bp相邻核苷酸的背景下为假定的hsa-miR-519c连接位点创建的二级RNA折叠结构),从而加强我们对miRNA靶向的电子预测。

验证S1细胞ABCG2 mRNA 3′UTR中假定的hsa-miR-519c靶点。

先前的研究表明,miRNA结合位点足以使报告基因结构中的miRNA依赖性基因沉默(10,53). 我们将不同长度的ABCG2 3′UTR直接插入荧光素酶开放阅读框下游,在组成活性SV40启动子的控制下(图。(图4A;4A级; 所有报告结构的全部细节可以在补充材料的图S9中找到)。在本试验中,荧光素酶活性水平与荧光素素酶mRNA的翻译效率和/或稳定性相关。FL ABCG2 3′UTR(nt 2462至4430)与不插入任何3′UTR的pGL3对照组相比,荧光素酶活性显著降低了90%(图。(图4B)。4B类). 通过在荧光素酶开放阅读框下游插入ABCG2(nt−1289 to+396)的5′UTR构建物对报告酶活性的最小抑制(减少<20%)也证明了3′UTR对荧光素素酶活性的特异性抑制。从ABCG2 3′UTR的3′端进行缺失的Reporter结构导致荧光素酶活性受到越来越少的抑制,这表明3′UTR3′端存在抑制因子(图。(图4B)。4B类). 重要的是,FL ABCG2 3′UTR构建物删除了假定的hsa-miR-519c结合位点(Δ3820-3841),将荧光素酶活性恢复到pGL3对照水平的40%左右。此外,仅hsa-miR-519c结合位点也可以将荧光素酶活性抑制到约20%的pGL3对照水平。然而,假定的hsa-miR-519c结合位点的扰乱序列对荧光素酶活性没有任何影响。

为了证实hsa-miR-519c是抑制ABCG2 3′UTR荧光素酶报告子结构活性的特异性miRNA,我们进行了实验,在S1细胞中,通过使用合成的mi RNA,改变了hsa-miR 519c的内源性池里迪安hsa-miR-519c抑制剂或模拟物(Dharmacon,Lafayette,CO)。我们研究中使用的特异性hsa-miR-519c抑制剂是一种化学修饰的反义寡核苷酸,与成熟的miRNA完全互补,与内源性miRNA紧密结合,导致其隔离(8). 特定的hsa-miR-519c抑制剂在最高剂量120 nM时将FL 3′UTR报告子结构中的荧光素酶活性恢复了约四倍(图。(图4C)。4摄氏度). 另一方面,hsa-miR-519c模拟物是一种含有成熟miRNA序列的双链RNA。当添加到细胞中时,模拟物的功能与内源性miRNA相同。当特定的hsa-miR-519c模拟物与含有ABCG2的FL3′UTR序列的报告质粒在S1细胞中共同转染时,与单独的报告结构相比,荧光素酶活性进一步呈剂量依赖性降低(图。(图4C)。4摄氏度). 作为对照,miRNA模拟物和抑制剂的通用阴性对照基于秀丽隐杆线虫miRNA序列(cel-miR-67[Dharmacon];该miRNA已被证实与人类、小鼠和大鼠的miRNA具有最小的序列一致性)对荧光素酶活性没有任何影响(图。(图4C)。4摄氏度). 此外,hsa-miR-519c模拟物和抑制剂均未对hsa-miR519c位点缺失的FL ABCG2 3′UTR报告子结构产生影响(参见补充材料中的图S10)。这些结果表明,hsa-miR-519c与人ABCG2 3′UTR中鉴定的靶结合位点特异性相互作用,从而抑制ABCG2 mRNA和蛋白质的表达。

由于推测的hsa-miR-519c位点已在S1细胞中验证为ABCG2 3′UTR,我们使用转染了特定hsa-miR519c抑制剂的S1细胞重复了3′RACE分析。我们使用特异性hsa-miR-519c抑制剂来增加S1细胞中ABCG2 3′UTR的水平(参见补充材料中的图S11A),随后通过巢式PCR 1和巢式PCR 2证实了S1细胞的长序列及其在S1MI80细胞中的缺失(参见补充材料中的图S11B)。

hsa-miR-519c通过翻译抑制和mRNA降解调节ABCG2的表达。

miRNAs的调节可能是由于翻译抑制和/或mRNA降解。为了评估hsa-miR-519c是否通过这些机制中的任何一种在下调ABCG2表达中发挥作用,需要一个ABCG2表达式载体(pcDNA3-ABCG2)及其两个衍生物,其中一个具有三个hsa-miR519c结合位点序列(pcDNA2-ABCG2-3×miRNA结合位点)的串联重复序列另一个串联重复三个hsa-miR-519c位点(pcDNA3-ABCG2-3×xscramble结合位点)的干扰序列,在载体中ABCG2编码序列的终止密码子后立即融合,转染到HEK293细胞中,并测量ABCG2的表达水平。作为内部对照,还评估了共转染GFP(pEGFP-C1;BD Bioscience Clontech,Mountain View,CA)的水平,并用于使ABCG2 mRNA水平正常化。在mRNA(−24%)和蛋白质(甚至更低,−49%)水平下,经3x×miRNA结合位点载体转染的细胞中ABCG2的表达低于未插入任何3′UTR的原始ABCG2表达载体转染细胞(参见补充材料中的图S12A)。加扰结合位点载体给出了与普通ABCG2表达载体相似的ABCG2表达。将hsa-miR-519c抑制剂与pcDNA3-ABCG2-3×miR519c结合位点载体共转染,以验证hsa-miR-519c是否特异性抑制ABCG2的表达。通过特定的hsa-miR-519c抑制剂而非通用阴性对照miRNA抑制剂,ABCG2 mRNA和蛋白表达水平恢复到用普通pcDNA3-ABCG2载体转染的细胞中的水平(参见补充材料中的图S12A)。

hsa-miR-519c抑制剂对ABCG2 mRNA和蛋白表达的差异作用是可重复的,这使我们假设在hsa-miR519c抑制ABCG2的过程中,蛋白翻译阻断可能比mRNA降解发挥更重要的作用。miRNA-mRNA配对的程度被认为会影响阻遏机制,完美配对会导致mRNA断裂,但不完美配对会造成蛋白质翻译阻滞(23). 我们如上所述制备了ABCG2表达载体,但紧接着在终止密码子之后,我们放置了与miR519c序列(pcDNA3-ABCG2-3×miRNA互补序列)完全互补的22-nt序列的三个拷贝的串联重复序列或三个加扰序列副本的串联重复(pcDNA3-ABCG2-3×加扰互补序列)。然后用这些新载体转染HEK293细胞,检测ABCG2的表达。与之前的载体一样,发现用3×miRNA互补序列载体转染的细胞中ABCG2的表达低于用普通ABCG2表达载体转染细胞中的表达,但在mRNA(−61%)和蛋白质(−70%)水平上都相似(参见补充材料中的图S12B),表明在影响这两个水平的同时,hsa-miR-519c对蛋白质翻译的影响在野生型3′UTR中占主导地位。3×加扰互补序列载体与原ABCG2表达载体表达的ABCG2相似。有趣的是,当载体仅包含hsa-miR-519c结合位点或互补序列的单个拷贝时,也获得了类似的观察结果(图。5A和B).

保存图片、插图等的外部文件。对象名为zmb0170876490005.jpg

hsa-miR-519c结合位点通过mRNA降解和蛋白翻译阻断促进ABCG2表达的调节。(A) 用ABCG2表达载体(pcDNA3-ABCG2)或其两种衍生物中的一种转染HEK293细胞,一种带有hsa-miR-519c结合序列的拷贝(pcDNA2-ABCG2-1×miRNA结合位点),另一种带有hsa-miR-5129c位点的打乱序列(pcDNA1-ABCG2-可填充结合位点)在载体中ABCG2编码序列的终止密码子后立即融合,并在转染48 h后评估ABCG2 mRNA和蛋白的表达。一种特定的hsa-miR-519c或阴性对照miRNA抑制剂也被共同转染到细胞中,以观察其是否可以增加ABCG2的表达水平。ABCG2表达水平用GFP(pEGFP-C1;DB Bioscience Clontech)进行标准化,并绘制在条形图中,与用pcDNA3-ABCG2载体转染的细胞的水平相关。所示数据代表三个独立实验的平均值±SD。(B) hsa-miR-519c和3′UTR序列之间的完美Watson-Crick配对主要促进ABCG2抑制中的mRNA降解。用与A组相似的ABCG2表达载体转染HEK293细胞。然而,在两个衍生载体中,不是hsa-miR-519c结合位点,而是与hsa-miR519c连接位点互补的22-nt序列的一个拷贝(pcDNA3-ABCG2-1×miRNA互补序列)或其扰码序列(pcDNA3-ABCG2-可压缩互补序列)在ABCG2编码序列的下游立即融合。发现用1×miRNA互补序列载体转染的细胞中ABCG2的表达低于用普通ABCG2表达载体转染而在终止密码子后没有任何东西的细胞,但在mRNA(−65%)和蛋白质(−70%)水平上的表达相似。该加扰互补序列向量给出了与原始ABCG2表达向量相似的ABCG2表达式。hsa-miR-519c也被发现负责阻遏,因为hsa-miR519c抑制剂,而非干扰控制miRNA抑制剂,可以恢复转染pcDNA3-ABCG2-1×miRNA互补序列的细胞中ABCG2的表达。ABCG2的表达水平用GFP标准化,并与转染pcDNA3-ABCG2载体的细胞中获得的水平相关。

在转染了各种载体(具有hsa-miR-519c结合位点或互补序列的单个拷贝)的HEK293细胞中,也测量了ABCG2 mRNA的稳定性。在转染后48小时加入放线菌素D以阻断转录,在不同时间间隔分离总RNA,并使用RT-PCR分析ABCG2 mRNA水平。作为一项内部控制,对共转染GFP的水平进行测量以进行归一化。如图所示。图6,6,ABCG2编码序列下游具有hsa-miR-519c互补序列的载体表达的ABCG2 mRNA是最不稳定的(半衰期[t吨1/2],~6 h),然后是带有野生型hsa-miR-519c结合位点的载体(t吨1/2,>16小时)。相反,在放线菌素D处理后24小时内,无任何3′UTR的ABCG2相当稳定。使用通用探针库基因检测,通过实时PCR重复检测,获得了类似的结果。这些结果强调了ABCG2 3′UTR(nt 3820至3841)中hsa-miR-519c结合位点的体内失稳作用,但仅部分通过促进mRNA降解。有趣的是,hsa-miR-519c与hsa-miR-519c互补序列的完美配对导致ABCG2 mRNA的最大降解。

保存图片、插图等的外部文件。对象名为zmb0170876490006.jpg

hsa-miR-519c结合位点通过促进mRNA降解降低ABCG2的表达。(顶部)ABCG2表达载体(pcDNA3-ABCG2)及其两种衍生物,一种具有hsa-miR-519c结合序列的拷贝(pcDNA3-ABCG2-1×miRNA结合位点),另一种具有与hsa-miR-519c互补的序列(pcDNA3-ABCG2-1×miRNA互补序列)在载体中ABCG2编码序列的终止密码子后立即融合,转染到HEK293细胞中,并监测ABCG2 mRNA的稳定性。在转染后48小时加入放线菌素D以阻断转录,在不同时间间隔分离总RNA,并使用RT-PCR分析ABCG2 mRNA水平。作为内部对照,评估GFP mRNA水平(来自pEGFP-C1[BD Bioscience Clontech]的共转染),并用于使ABCG2 mRNA水平正常化。该图显示了三个独立且可重复的实验的平均结果。ABCG2编码序列下游带有hsa-miR-519c互补序列的载体表达的ABCG2 mRNA最不稳定(t吨1/2,~6 h),然后是带有hsa-miR-519c结合位点的载体(t吨1/2,>16小时)。相反,在放线菌素D处理后24小时内,没有任何3′UTR的ABCG2 mRNA相当稳定。(底部)代表性凝胶图像,显示顶部面板中总结的RT-PCR分析。

hsa-miR-519c调节内源性ABCG2 mRNA和蛋白表达。

最后,我们评估了hsa-miR-519c在下调内源性ABCG2表达中是否具有功能作用。选择A549细胞进行我们的研究是因为通过Western blot分析可以在这些细胞中检测到内源性ABCG2蛋白的表达,并且与我们研究中使用的其他亲本细胞系类似,A549细胞具有长和短ABCG2 3′UTR,如3′RACE分析所示(数据未显示)。用不同浓度的hsa-miR-519c抑制剂(一种与成熟miRNA具有互补序列的单链修饰RNA)转染A549细胞,导致ABCG2在mRNA和蛋白质水平上的剂量依赖性诱导。使用120 nM miR519c抑制剂,ABCG2 mRNA和蛋白分别上调了两倍和六倍,而GAPDH中没有发现变化(图。(图7,7,顶部面板)。阴性对照miRNA抑制剂没有改变ABCG2、GAPDH或hsa-miR-519c水平。尽管miRNA抑制剂主要通过隔离内源性miRNA发挥作用,但较高剂量的miRNA抑制剂也会降低内源性miRNAs的水平(图。(图7,7,顶部面板)。事实上,反义寡核苷酸已被证明可以降低同源miRNAs的细胞浓度(8,29).

保存图片、插图等的外部文件。对象名为zmb0170876490007.jpg

hsa-miR-519c调节A549细胞内源性ABCG2的表达。将0到120 nM的特异性hsa-miR-519c抑制剂、hsa-miR-519c模拟物或相应的阴性对照物(基于cel-miR-67序列[Dharmacon];该miRNA已被证实与人类、小鼠和大鼠的miRNA具有最小序列同源性)转染A549细胞48小时。采集总RNA和全细胞裂解物,用于随后的RT-PCR和ABCG2表达的Western blot分析。ABCG2的表达用GAPDH标准化,并与模拟处理细胞的表达相对应。hsa-miR-519c的转录物也通过干环RT-PCR测定。该图显示了三个独立且可重复的实验的代表性结果。通用miRNA阴性对照抑制剂或模拟物不影响ABCG2、hsa-miR-519c或GAPDH的表达。

在这些研究的同时,A549细胞被转染特定的hsa-miR-519c模拟物,这是一种模拟内源性hsa-miR519c的合成RNA双链。120 nM的特定hsa-miR-519c模拟物分别抑制ABCG2 mRNA和蛋白水平约70%和90%,但GAPDH未受影响(图。(图7,7,底部面板)。阴性对照miRNA模拟物也没有改变ABCG2、GAPDH或hsa-miR-519c水平。

另一个抑制性miRNA位于ABCG2的3′UTR。

如荧光素酶报告物分析所示(图。(图4B),4B类)最短的ABCG2 3′UTR结构(nt 2462至2761)对报告活性的实质性抑制表明,在ABCG2 3'UTR的上游区域存在其他未知的抑制因子。在编写本报告时,Liao等人(34)据报道,hsa-miR-520h通过抑制ABCG2的表达来促进造血干细胞向祖细胞的分化。根据miRBase TARGETS数据库对miRNA靶点的预测(19)和RNAHybrid(43),假定的hsa-miR-520h结合序列位于更上游(nt 2494至2514,相对于序列中定义的第一个核苷酸NM_004827号)(图。(图3A;3A级; 参见图。图9)9)ABCG2 3′UTR比hsa-miR-519c高。更重要的是,hsa-miR-520h与亲本S1和耐药S1MI80细胞中的长短ABCG2 3′UTR结合。因此,我们通过S1和S1MI80细胞的干环RT-PCR检测hsa-miR-520h的表达水平,试图阐明其在调节ABCG2中的作用。有趣的是,S1MI80中的hsa-miR-520h水平仅为S1细胞的一半(43%±13%)(图。(图8A8安).

保存图片、插图等的外部文件。对象名为zmb0170876490008.jpg

(A) 亲本S1和耐药S1MI80细胞之间hsa-miR-520h表达水平的差异。通过干环RT-PCR测定hsa-miR-519c和hsa-miR520h的表达(7,46). U6小核糖核酸作为负荷对照。图中显示了三个独立且可重复的实验的代表性结果,在下面板中绘制了平均值(±SD)。(B) 推测的hsa-miR-520h结合位点在ABCG2 3′UTR,RNA杂交自由能(ΔG公司)计算,以及hsa-miR-520h的理论miRNA-mRNA双链配对。

保存图片、插图等的外部文件。对象名称为zmb017087649009.jpg

提出了miRNA参与3′UTR ABCG2调控的模型。一种推测的miRNA,hsa-miR-519c,与ABCG2的3′UTR结合,并抑制ABCG2在亲本细胞系中的表达。在耐药细胞系中,由于3′UTR较短,hsa-miR-519c不能与ABCG2 mRNA结合,因此mRNA降解和/或蛋白翻译抑制得到缓解,通过其他机制,如启动子的组蛋白修饰等,导致ABCG2过度表达(49). hsa-miR-520h是另一个假定的靶向ABCG2 mRNA的miRNA,这是在手稿准备期间报道的(34). 与我们研究中确定的hsa-miR-519c不同,在S1和S1MI80细胞中,hsa-miR-520h与ABCG2的3′UTR结合。有趣的是,S1MI80细胞比S1细胞表达更少的hsa-miR-520h,推测前者对ABCG2的抑制更少。

讨论

的过度表达ABCG2公司在用多种化疗药物筛选的癌细胞系中经常观察到该基因(1,12,35,38,44,52). 然而,对其上调的机制知之甚少。研究发现,基因扩增和染色体易位在基因表达增加中起作用ABCG2公司一些耐药细胞系中的基因(27). ABCG2公司药物选择细胞中的基因可能提供另一种新的机制ABCG2公司上调(39),这一发现类似于对MDR-1型,其中的5′区重排MDR-1型导致其他启动子捕获该基因(22,26,37). 最近,我们报道了一组允许的组蛋白修饰和染色质重塑因子Brg-1与ABCG2公司近端启动子是上调ABCG2公司在抗性细胞中(49). 为了进一步了解ABCG2公司在耐药癌细胞中,我们发现ABCG2 mRNA在一些敏感的亲代癌细胞系中的稳定性低于其药物选择和ABCG2过度表达的耐药细胞系。我们还观察到,亲代细胞产生ABCG2 mRNA和长3′UTR和短3′UTRs,而耐药细胞只产生短的UTRs。

由于真核生物mRNA的3′UTR含有影响mRNA稳定性、翻译和转运的调节元件,因此有理由怀疑,与长3′UTR-相比,短3′UTRs中缺少一些调节元件可能导致敏感和耐药细胞中ABCG2的差异调节。我们假设一个假定的miRNA与亲代细胞中ABCG2的3′UTR结合,从而抑制ABCG2表达,但该miRNA的结合位点在耐药细胞中较短的3′UTR中丢失。miRNAs在调节发育和分化过程中起着直接作用,在动植物物种中已鉴定出4000多个miRNAss(miRBase Registry[http://microrna.sanger.ac.uk/sequences网站]). miRNA失调与几种癌症的进展和侵袭性相关(,25). 电子分析显示,hsa-miR-519c的假定靶点位于敏感细胞中鉴定的长ABCG2 3′UTR转录物的远3′端。

事实上,最近有人提出了miRNA对不同长度的3′UTR替代转录物的差异抑制(32). 当miRNA靶点位于多聚腺苷酸化的两个变异位点之间时,较短的转录物将是无靶点的,可以逃避miRNA介导的抑制,而较长的转录物则会受到抑制。ABCG2 mRNA既包含典型的AAUAAA poly(a)信号,也包含其非经典变体AUUAAA。规范和非规范聚(A)信号通常处理方式不同。这为具有多个poly(a)位点的mRNA通过合成特定mRNA形式进行调节提供了一种方法,从而调节基因表达(14,18). 根据我们的3′RACE分析和ABCG2的后续DNA测序,似乎非经典AUUAAA poly(A)信号优先用于耐药细胞,仅提供短的3′UTR转录物。在亲代细胞中,使用标准AAUAAA和非标准AUUAAA信号(图。(图3;; 参见补充材料中的图S7),从而生成长和短3′UTR转录本。

我们对ABCG2是hsa-miR-519c真正靶点的假设进行了严格测试。基于当前预测算法(19),这种靶向性很可能是因为hsa-miR-519c与ABCG2 mRNA的3′UTR结合,其“7-mer-m8位点”与hsa-miR-519c种子序列互补(图。(图3B)。3B公司). 靶ABCG2 3′UTR区域中存在不稳定元素(发夹环),进一步增强了其成为直接靶的可能性(参见补充材料中的图S8)。利用报告基因分析,FL ABCG2 3′UTR显著抑制了报告基因的活性,而3′截断报告基因结构的抑制程度较小。此外,当从FL ABCG2 3′UTR中删除hsa-miR-519c结合位点时,对报告基因活性的抑制作用也显著降低。融合在荧光素酶报告基因编码序列下游的hsa-miR-519c结合序列仅能降低报告基因的活性,尽管这类结果通常需要三个miRNA结合序列拷贝和一个短连接序列(20,33,55). 此外,FL ABCG2 3′UTR荧光素酶结构体(而不是ABCG2 3'UTR 3′端截短的结构体)对hsa-miR-519c抑制剂或模拟物具有特异性反应(图。(图4C;4摄氏度; 参见补充材料中的图S10)。在瞬时转染实验中,通过比较ABCG2编码序列下游含有hsa-miR-519c结合位点的表达载体与含有扰码序列的表达载体,证明了hsa-miR519c的抑制作用。表达载体中hsa-miR-519c位点的突变消除了抑制作用。最后,在A549人肺癌细胞中,内源性ABCG2被特定的hsa-miR-519c抑制剂上调,而被特定的hsa-miR-5129c模拟物下调(图。(图7)。7). 因此,多项测试证实ABCG2是hsa-miR-519c的真正靶点。这些测试也被用于确认多种其他miRNA靶对,包括let-7-RAS(25),miR-1-手2(55),miR-127-BCL6(46),miR-130a-GAX(7)和miR-196-HOXB8(53).

尽管miRNAs抑制基因表达可以通过mRNA降解和/或抑制现有mRNA的翻译来介导(51)我们的结果表明,hsa-miR-519c抑制ABCG2蛋白翻译的效果比促进mRNA降解的效果更为显著,这可以通过hsa-miR519c与人类ABCG2 3′UTR的次优配对程度来预测。当同一ABCG2表达载体中的hsa-miR-519c结合位点被与hsa-miR-519c互补的序列替换时,ABCG2主要受到mRNA降解的抑制(图。5A和B和6)。6). miRNA和靶mRNA之间完美的Watson-Crick互补性允许mRNA按照实验使用小干扰RNA(siRNA)时发生基因沉默的相同机制进行切割(23).

本研究中使用的耐药S1MI80亚系是从LS174人结肠癌细胞系的S1克隆中通过米托蒽醌筛选获得的(38). 有趣的是,在S1MI80细胞中观察到染色体带4q21-q22没有扩增,但ABCG2基因下游的t(4,7)平衡易位(27). 在本研究中,我们观察到S1MI80细胞中没有长的3′UTR,只留下短的UTR,这与这种易位是一致的。由于这种易位的断点被认为是ABCG2公司基因,易位与抗性细胞中采用的差异3′UTR长度之间是否存在任何相关性尚待确定。

总之,我们已经证明ABCG2在亲代和耐药细胞的3′UTR处采用不同的多聚腺苷酸化位点(因此3′UTRs长度不同),并且选择性多聚腺苷化与hsa-miR-519c靶向之间存在相互作用。图99提出了一种模型,其中hsa-miR-519c仅在亲代细胞而非耐药细胞中靶向ABCG2的3′UTR,从而允许ABCG2通过基因启动子处的组蛋白修饰和染色质重塑事件在后者中过度表达(49). ABCG2 mRNA在hsa-miR-519c存在下通过其他机制高度表达,如基因扩增、易位或不同组蛋白修饰,面临选择压力,以避免这种miRNA的结合位点。当手稿正在准备时,另一种推测的miRNA(hsa-miR-520h)据报道,在造血干细胞分化为祖细胞的过程中,靶向ABCG2(34). 与我们研究中确定的hsa-miR-519c不同,在S1和S1MI80细胞中,hsa-miR-520h与ABCG2的3′UTR结合。有趣的是,S1MI80细胞表达的hsa-miR-520h(~40%)少于S1细胞。

可以考虑这项工作的两个影响。如果可以在3′UTR中识别出其他抑制性元件或其他抑制性miRNA结合位点,则可以生成含有这种miRNA的载体,以传递给患者肿瘤,以克服ABCG2过度表达对药物积累的限制,前提是可以克服将基因治疗转化为临床的问题。第二,药物转运蛋白基因的遗传变异(尤其是单核苷酸多态性[SNP])与改变癌症药物药理学和药物效应有关。迄今为止,大多数SNP在ABCG2公司位于其编码区域内。由于我们已经在ABCG2 3′UTR上确定了一个miRNA结合位点,因此SNP变异(如果在3′UTRs上发现)可以与miRNA结合相互作用,从而影响ABCG2的整体表达。因此,这些SNP的影响是可以预测的,并为研究转运蛋白活性的个体间变异开辟了另一条途径。毫无疑问,继续努力描述miRNA在调节ABCG2中的功能,最终可能提供一种改变ABCG2表达的新方法,以操纵ABCG2底物药物的处置。

补充材料

[补充材料]

致谢

这项研究得到了NIH、国家癌症研究所、癌症研究中心的内部研究项目的支持。

脚注

2008年6月23日提前出版。

本文的补充材料可以在http://mcb.asm.org/.

参考文献

1Allikmets,R.、L.M.Schriml、A.Hutchinson、V.Romano-Spica和M.Dean。染色体4q22上与多药耐药性有关的人类胎盘特异性ATP-结合盒基因(ABCP)。癌症研究。 585337-5339. [公共医学][谷歌学者]
2Ambros,V.2004年。动物microRNA的功能。自然 431350-355. [公共医学][谷歌学者]
3Bandres,E.、E.Cubedo、X.Agirre、R.Malumbres、R.Zarate、N.Ramirez、A.Abajo、A.Navarro、I.Moreno、M.Monzo和J.Garcia Foncillas。2006.通过实时PCR鉴定在结直肠癌和非肿瘤组织中差异表达的13种成熟微小RNA。摩尔癌症 529[PMC免费文章][公共医学][谷歌学者]
4C.F.Basler和M.S.Horwitz。B亚组35型腺病毒早期区域3 mRNA与C亚组腺病毒不同。病毒学 215165-177. [公共医学][谷歌学者]
5Brennecke,J.、A.Stark、R.B.Russell和S.M.Cohen。2005.microRNA-目标识别原理。《公共科学图书馆·生物学》。 e85。[PMC免费文章][公共医学][谷歌学者]
6Chen,C.,D.A.Ridzon,A.J.Broomer,Z.Zhou,D.H.Lee,J.T.Nguyen,M.Barbisin,N.L.Xu,V.R.Mahuvakar,M.R.Andersen,K.Q.Lao,K.J.Livak,K.J Guegler。2005.通过干环RT-PCR实时定量microRNA。核酸研究。 33e179。[PMC免费文章][公共医学][谷歌学者]
7Chen,Y.和D.H.Gorski。2008.通过下调抗血管生成同源框基因GAX和HOXA5的微RNA(miR-130a)调节血管生成。血液 1111217-1226.[PMC免费文章][公共医学][谷歌学者]
8Davis,S.、B.Lollo、S.Freier和C.Esau。2006.用反义寡核苷酸改进miRNA的靶向性。核酸研究。 342294-2304.[PMC免费文章][公共医学][谷歌学者]
9Didiano,D.和O.Hobert。2006.完美的种子配对通常不是miRNA-靶相互作用的可靠预测因子。自然结构。分子生物学。 13849-851. [公共医学][谷歌学者]
10Doench,J.G.和P.A.Sharp。2004.翻译抑制中microRNA靶点选择的特异性。基因发育。 18504-511.[PMC免费文章][公共医学][谷歌学者]
11Donnini,M.、A.Lapucci、L.Papucci,E.Witort、A.Jacquier、G.Brewer、A.Nicolin、S.Capaccioli和N.Schiavone。2004年。TINO的鉴定:一种新的进化保守的BCL-2 AU-富元件RNA-结合蛋白。生物学杂志。化学。 27920154-20166. [公共医学][谷歌学者]
12Doyle,L.A.、W.Yang、L.V.Abruzzo、T.Krogmann、Y.Gao、A.K.Rishi和D.D.Ross。人类MCF-7乳腺癌细胞的多药耐药转运蛋白。程序。国家。阿卡德。科学。美国 9515665-15670.[PMC免费文章][公共医学][谷歌学者]
13Ebert,B.、A.Seidel和A.Lampen。2005.BCRP作为在Caco-2细胞中代谢形成的苯并[a]芘结合物的转运体的鉴定及其由Ah-受体激动剂的诱导。致癌作用 261754-1763. [公共医学][谷歌学者]
14Edwards-Gilbert、G.、K.L.Veraldi和C.Milcarek。1997.复杂转录单位中的替代性聚(A)位点选择:意味着结束?核酸研究。 252547-2561.[PMC免费文章][公共医学][谷歌学者]
15Ee、P.L.、S.Kamalakaran、D.Tonetti、X.He、D.D.Ross和W.T.Beck。乳腺癌耐药蛋白(ABCG2)基因中一种新的雌激素反应元件的鉴定。癌症研究。 641247-1251. [公共医学][谷歌学者]
16Farh,K.K.,A.Grimson,C.Jan,B.P.Lewis,W.K.Johnston,L.P.Lim,C.B.Burge和D.P.Bartel。2005.哺乳动物微小RNA对mRNA抑制和进化的广泛影响。科学类 3101817-1821. [公共医学][谷歌学者]
17Frohman,M.A.1990年。RACE:cDNA末端的快速扩增,第28-38页。M.A.Inns、D.H.Gelfand、J.J.Sninsky和T.J.White(编辑),PCR协议:方法和应用指南。学术出版社,加利福尼亚州圣地亚哥。
18Graber,J.H.、C.R.Cantor、S.C.Mohr和T.F.Smith。1999生物信息学控制信号检测:不同物种的mRNA 3′末端处理序列。程序。国家。阿卡德。科学。美国 9614055-14060.[PMC免费文章][公共医学][谷歌学者]
19Griffiths-Jones,S.2006年。miRBase:微RNA序列数据库。方法分子生物学。 342129-138. [公共医学][谷歌学者]
20Grimson,A.、K.K.Farh、W.K.Johnston、P.Garrett-Engele、L.P.Lim和D.P.Bartel。2007.哺乳动物中微RNA靶向特异性:种子配对以外的决定因素。分子电池 2791-105.[PMC免费文章][公共医学][谷歌学者]
21他、L.和G.J.汉农。2004.微小RNA:在基因调控中发挥重要作用的小RNA。Nat.Rev.基因。 5522-531. [公共医学][谷歌学者]
22哈夫、L.M.、J.S.李、R.W.罗比和T.福乔。2006.导致MDR-1激活的基因重排特征。生物学杂志。化学。 28136501-36509. [公共医学][谷歌学者]
23Hutvagner,G.和P.D.Zamore。2002年,多重翻转RNAi酶复合物中的microRNA。科学类 2972056-2060. [公共医学][谷歌学者]
24Jensen、L.E.和A.S.Whitehead。2004年。膜结合IL-1R辅助蛋白mRNA的3′非翻译区导致组织特异性失稳。免疫学杂志。 1736248-6258. [公共医学][谷歌学者]
25Johnson,S.M.、H.Groshans、J.Shingara、M.Byrom、R.Jarvis、A.Cheng、E.Labourier、K.L.Reinert、D.Brown和F.J.Slack。2005年。RAS受let-7 microRNA家族调控。单元格 120635-647. [公共医学][谷歌学者]
26Knutsen,T.、L.A.Mickley、T.Ried、E.D.Green、S.du Manoir、E.Schrock、M.Macville、Y.Ming、R.Robey、M.Polymeropous、R.Torres和T.Fojo。1998年,随机染色体重排激活药物选择细胞系和药物难治性ALL患者中的耐药基因MDR1/P-糖蛋白的细胞遗传学和分子特征。基因染色体癌 2344-54. [公共医学][谷歌学者]
27Knutsen,T.、V.K.Rao、T.Ried、L.Mickey、E.Schneider、K.Miyake、B.M.Ghadimi、H.Padilla-Nash、S.Pack、L.Greenberger、K.Cowan、M.Dean、T.Fojo和S.Bates。2000.在独立衍生的米托蒽醌耐药细胞系中扩增4q21-q22和MXR基因。基因染色体癌 27110-116. [公共医学][谷歌学者]
28Krishnamurthy,P.、D.D.Ross、T.Nakanishi、K.Bailey Dell、S.Zhou、K.E.Mercer、B.Sarkadi、B.P.Sorrentino和J.D.Schuetz。干细胞标志物Bcrp/ABCG2通过与血红素的相互作用增强缺氧细胞的存活。生物学杂志。化学。 27924218-24225. [公共医学][谷歌学者]
29Krutzfeldt,J.、N.Rajewky、R.Braich、K.G.Rajeev、T.Tuschl、M.Manoharan和M.Stoffel。2005年,使用“antagomirs”在体内沉默microRNAs。自然 438685-689. [公共医学][谷歌学者]
30Lander,E.S.,L.M.Linton,B.Birren等人,2001年。人类基因组的初步测序和分析。自然 409860-921. [公共医学][谷歌学者]
31Lee,N.S.、T.Dihjima、G.Bauer、H.Li、M.J.Li,A.Ehsani、P.Salvatera和J.Rossi。针对HIV-1 rev转录物的小干扰RNA在人类细胞中的表达。自然生物技术。 20500-505. [公共医学][谷歌学者]
32Legendre,M.、W.Ritchie、F.Lopez和D.Gautheret。2006.替代转录物的差异抑制:miRNA靶点筛选。公共科学图书馆计算。生物。 2e43。[PMC免费文章][公共医学][谷歌学者]
33刘易斯、B.P.、C.B.伯格和D.P.巴特尔。2005年。通常两侧有腺苷的保守种子配对表明,数千人类基因是microRNA的靶标。单元格 12015-20. [公共医学][谷歌学者]
34Liao,R.、J.Sun、L.Zhang、G.Lou、M.Chen、D.Zhou、Z.Chen和S.Zhang。2008.微RNA在人类造血干细胞的发育中发挥作用。细胞生物化学杂志。 104805-817. [公共医学][谷歌学者]
35Maliepard,M.、M.A.van Gastelen、L.A.de Jang、D.Pluim、R.C.van Waardenburg、M.C.Ruevekamp-Helmers、B.G.Floot和J.H.Schellens。1999.拓扑替康选择的肿瘤细胞系中BCRP/MXR/ABCP基因的过度表达。癌症研究。 594559-4563. [公共医学][谷歌学者]
36Meister,G.、M.Landthaler、Y.Dorsett和T.Tuschi。2004.序列特异性抑制microRNA和siRNA诱导的RNA沉默。核糖核酸 10544-550.[PMC免费文章][公共医学][谷歌学者]
37Mickley,L.A.、B.A.Spengler、T.A.Knutsen、J.L.Biedler和T.Fojo。基因重排:MDR-1基因激活的新机制。临床杂志。投资。 991947年至1957年。[PMC免费文章][公共医学][谷歌学者]
38Miyake,K.、L.Mickley、T.Litman、Z.Zan、R.Robey、B.Cristensen、M.Brangi、L.Greenburger、M.Dean、T.Fojo和S.E.Bates。米托蒽醌耐药细胞中高过表达cDNA的分子克隆:与ABC转运基因同源性的证明。癌症研究。 598-13. [公共医学][谷歌学者]
39Nakanishi,T.、K.J.Bailer-Dell、B.A.Hassel、K.Shiozawa、D.M.Sullivan、J.Turner和D.D.Ross。BCRP mRNA新的5′非翻译区变异体在药物选择的癌细胞和正常人体组织中差异表达:对耐药性、组织特异性表达和替代启动子使用的影响。癌症研究。 665007-5011. [公共医学][谷歌学者]
40Nishimori,T.、H.Inoue和Y.Hirata。环氧合酶-2基因3′-非翻译区通过糖皮质激素受体参与其转录后调控。生命科学。 742505-2513之间。[公共医学][谷歌学者]
41Pesole,G.和S.Liuni。1999.用于真核生物mRNA 5′和3′非翻译区功能分析的互联网资源。趋势Genet。 15378. [公共医学][谷歌学者]
42Pillai,R.S.2005年。微小RNA的功能:微小RNA的多种机制?核糖核酸 111753-1761.[PMC免费文章][公共医学][谷歌学者]
43Rehmsmeier,M.、P.Steffen、M.Hochsmann和R.Giegerich。2004.快速有效地预测microRNA/目标双工体。核糖核酸 101507-1517.[PMC免费文章][公共医学][谷歌学者]
44Robey,R.、W.Y.Medina-Perez、K.Nishiyama、T.Lahusen、K.Miyake、T.Litman、A.M.Senderowicz、D.D.Ross和S.E.Bates。2001.ATP-结合盒半转运体ABCG2(Mxr/BCrp/ABCP1)在黄哌啶耐药的人类乳腺癌细胞中的过度表达。临床。癌症研究。 7145-152. [公共医学][谷歌学者]
45Roe、D.F.、G.L.Craviso和J.C.Waymire。2004年。烟碱刺激以需要转录的方式调节酪氨酸羟化酶mRNA半衰期和与3′UTR结合的蛋白质。大脑研究摩尔大脑研究。 12091-102. [公共医学][谷歌学者]
46Saito,Y.、G.Liang、G.Egger、J.M.Friedman、J.C.Chuang、G.A.Coetzee和P.A.Jones。2006.通过染色质修饰药物下调人类癌细胞中原癌基因BCL6,特异性激活microRNA-127。癌细胞 8435-443. [公共医学][谷歌学者]
47Schiavone,N.、P.Rosini、A.Quattrone、M.Donnini、A.Lapucci、L.Citti、A.Bevilacqua、A.Nicolin和S.Capaccioli。bcl-2 mRNA 3′非翻译区中的一个保守的AU-rich元件具有不稳定功能,参与凋亡过程中bcl-2的下调。美国财务会计准则委员会J。 14174-184. [公共医学][谷歌学者]
48Szatmari,I.、G.Vamosi、P.Brazda、B.L.Balint、S.Benko、L.Szeles、V.Jeney、C.Ozvegy-Laczka、A.Szanto、E.Barta、J.Balla、B.Sarkadi和L.Nagy。过氧化物酶体增殖物激活受体γ-调节的ABCG2表达对人类树突状细胞具有细胞保护作用。生物学杂志。化学。 28123812-23823. [公共医学][谷歌学者]
49To、K.K.W.、O.Polgar、L.Huff、M.Kuniaki和S.E.Bates。2008年。组蛋白修改ABCG2公司HDAC抑制剂治疗后的启动子反映了多药耐药细胞的情况。摩尔癌症研究。 6151-164.[PMC免费文章][公共医学][谷歌学者]
50收件人:K.K.W.、Z.Zhan和S.E.Bates。肾癌中ABCG2基因启动子甲基化异常。分子细胞。生物。 268572-8585.[PMC免费文章][公共医学][谷歌学者]
51Valencia-Sanchez,M.A.、J.Liu、G.J.Hannon和R.Parker。2006.通过miRNAs和siRNAs控制翻译和mRNA降解。基因发育。 20515-524. [公共医学][谷歌学者]
52Volk,E.L.、K.M.Farley、Y.Wu、F.Li、R.W.Robey和E.Schneider。野生型乳腺癌耐药蛋白的过度表达介导了甲氨蝶呤耐药。癌症研究。 625035-5040. [公共医学][谷歌学者]
53Yekta,S.、I.H.Shih和D.P.Bartel。2004.MicroRNA指导的HOXB8mRNA切割。科学类 304594-596. [公共医学][谷歌学者]
54Yu、H.、S.Stasinopoulos、P.Leedman和R.L.Medcalf。纤溶酶原激活物抑制剂2型mRNA的固有不稳定性由三四脯氨酸调节。生物学杂志。化学。 27813912-13918. [公共医学][谷歌学者]
55Zhao,Y.、E.Samal和D.Srivastava。2005.血清反应因子调节心肌发生过程中靶向Hand2的肌肉特异性microRNA。自然 436214-220. [公共医学][谷歌学者]
56Zuker,M.2003年。用于核酸折叠和杂交预测的Mfold web服务器。核酸研究。 313406-3415.[PMC免费文章][公共医学][谷歌学者]

文章来自分子和细胞生物学由以下人员提供泰勒和弗朗西斯