布鲁氏菌病的外泌体miRNA标志物及应用的制作方法

布鲁氏菌病的外泌体mirna标志物及应用
技术领域
1.本发明属于生物检测技术领域，尤其涉及布鲁氏菌病的外泌体mirna标志物及应用。


背景技术：

2.布鲁氏菌病(brucellosis，简称“布病”)是由布鲁氏菌属细菌侵入机体，引起人畜共患的传染性变态反应性疾病，多因接触患病的牛、羊等牲畜及其制品感染，发病隐匿、诊断延迟和致残率高，对患者的健康尤其危险，不具有典型的临床症状并易形成慢性感染性疾病。目前已经有12种布鲁氏菌被报道，其中引起人类感染的最常见的布鲁氏菌包括羊种、牛种和猪种，其中羊种布鲁氏菌毒力最强。
3.随着畜牧业的迅速发展，我国牲畜饲养量逐渐增加，畜牧产品的流通速度加快，使布鲁氏菌病扩散速度增快，流行范围也不断扩大，布鲁氏菌病感染已成为我国重要的公共卫生问题之一。
4.微小rna(microrna,mirna)是一类约为22个核苷酸大小的非编码单链rna，对多种生物功能具有重要的调节作用，可作为疾病的分期、分级、疾病预测以及临床干预措施监测等生物标志物。目前，布鲁氏菌病患者的精确临床诊断面临着许多挑战。血清mirnas的表达水平往往与机体状态存在密切关系，mirnas有望成为布鲁氏菌病精确诊断和预后判断的良好标志物。近年来，国内外研究学者已发表了许多有关mirnas作为生物标志的研究报道。但是相比于在肿瘤中已知的重要作用，mirnas在感染性疾病尤其是细菌感染所致疾病中的作用研究还较少。


技术实现要素：

5.本发明的目的在于提供布鲁氏菌病的外泌体mirna标志物及应用，以解决上述技术问题。
6.本发明的目的之一在于提供布鲁氏菌病的外泌体mirna标志物，包括hsa-mir-11400、hsa-mir-199a-5p或hsa-mir-148a-5p中的一种或两种以上；hsa-mir-11400的核苷酸序列为seq id no:3，hsa-mir-199a-5p的核苷酸序列为seq id no:5，hsa-mir-148a-5p的核苷酸序列为seq id no:7。
7.本发明目的之二在于提供上述用于检测上述布鲁氏菌病的外泌体mirna标志物的试剂在制备检测试剂盒中的应用。
8.本发明的有益效果：本发明提供的外泌体hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p具有布鲁氏菌病诊断潜力，可以作为布鲁氏菌病精确诊断和预后判断的生物标志物。
附图说明
9.图1为血清外泌体的鉴定；a：血清来源外泌体中标记蛋白cd9和tsg101的蛋白质印
迹；b：使用zetasizer纳米zs90粒度仪测量的外泌体直径；c：从血清中纯化的外泌体的电子显微照片(比例尺500nm)；
10.图2为差异表达的血清外泌体mirna的热图；
11.图3为布鲁氏菌病患者中五种外泌体mirna的相对表达；a：hsa-mir-194-5p的表达水平；b：hsa-mir-11400的表达水平；c：hsa-mir-148a-5p的表达水平；d：hsa-mir-199a-5p的表达水平；e：hsa-mir-23b-3p的表达水平；
12.图4为通过roc曲线分析确定外泌体mirna的诊断能力；a：hsa-mir-11400；b：hsa-mir-199a-5p；c：hsa-mir-148a-5p；d：hsa-mir-11400和hsa-mir-199a-5p；e：hsa-mir-11400和hsa-mir-148a-5p；f：hsa-mir-148a-5p和hsa-mir-199a-5p；g：hsa-mir-11400,hsa-mir-148a-5p和hsa-mir-199a-5p；
13.图5为基因本体论(go)和京都基因和基因组百科全书(kegg)对血清衍生的外泌体mirnas(hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p)的免疫相关基因的分析；a：对明显富集的靶基因的go分析；在生物过程(bp)、细胞成分(cc)和分子功能(mf)术语中，列出了前5个显著富集的项目；p值、富集因子和mirna靶点数量被用来确定go的富集程度；b：通过kegg进行mirna靶点分析的途径富集度分析。allenrichment分数被评估为p值的负对数(-log p)；
14.图6为mirna-免疫相关基因调控网络。方块代表基因，椭圆代表mirnas；
15.图7为根据gse69597数据确定18个免疫基因的mrna表达；a：nfya；b：cyld；c：traf3；d：cdh1；e：plxna2；f：lepr；g：tgfb2；h：gnrhr；i：il-17ra；j：rora；k：cd48；l：prkca；m：tnf；n：actg1；o：ptgdr；p：cd22；q：acvr1b；r：cbl。
具体实施方式
16.下面通过具体实施方式进一步详细说明：
17.实施例1筛选布鲁氏菌病的外泌体mirna标志物
18.以下实验通过了伦理审批，并获得了所有参与者的知情同意。实验主要分为布鲁氏菌病组和健康对照组，参与者的信息参见表1。
19.布鲁氏菌病组入选标准：
20.根据中华人民共和国卫生行业标准《布鲁氏菌病诊断》(ws269-2019)，临床诊断为原发性布鲁氏菌症；年龄30-50岁。
21.布鲁氏菌可以从任何病理物质培养物中分离出来，如患者的血液、骨髓、其他体液和排泄物初筛中玫瑰红-孟加拉平板试验(rbpt)阳性，血清凝集试验(sat)，滴度≥1:100或病程持续一年以上仍有症状，滴度为1:50参与者没有并发糖尿病、艾滋病和其他自身免疫疾病。
22.排除标准：已接受药物治疗的布鲁氏菌病患者；患有传染病、基础疾病和肿瘤相关疾病者。
23.表1参与者信息及临床特征
[0024][0025]
以下实验的第一步和第二步仅涉及3名布鲁氏菌病患者和3名健康对照者外周血用于高通量测序。第三步验证血清外泌体mirna时，每组增加22名参与者，共涉及25名布鲁氏菌病患者和25名健康对照者。
[0026]
以下实验的所有数据均通过graphpad prism软件ver.8或r版本4.1.3进行统计分析。数据表示为平均值
±
标准差。基于方差的同质性，使用双尾student’s t检验或mann-whitney u检验对两组进行比较。使用roc曲线评估血清外泌体mirna对布鲁氏菌病的诊断能力。p《0.05被认为具有统计学意义。
[0027]
第一步：外泌体的分离和表征
[0028]
为了分离外泌体，3名布鲁氏菌病患者和3名健康者的血清样本预先通过0.22μm过滤器过滤。按照制造商的说明，使用qiagen exorneasy maxi kit试剂盒(qiagen，hilden，德国)分离外泌体。分离的外泌体用pbs稀释后，使用zetasizer纳米zs90粒度仪测量外泌体的直径。再通过透射电子显微镜(tem)观察外泌体的微观结构。
[0029]
为了验证外泌体是从布鲁氏菌病和健康对照组的血清中分离出来的，使用western blotting蛋白质印迹法分析评估外泌体特异性蛋白(tsg101和cd9)。
[0030]
将适量的ripa裂解物(ripa lysate)添加到分离的外泌体中以提取蛋白质。通过biosharp bca蛋白浓度测定试剂盒(biosharp，中国，货号bl521a)检测样品的蛋白质浓度为2680.07ng/ul。
[0031]
将等量的蛋白质样品加入到12％ sds聚丙烯酰胺凝胶中进行电泳。电泳后，将条带转移到pvdf膜上，并用脱脂奶粉密封1.5h；在4℃下将该膜与一抗抗兔多克隆cd9抗体(1:500，nbp2-67310，novus biologicals)和一抗抗小鼠多克隆tsg101抗体(1:5000，nb200-112，novus生物制品)孵育过夜。然后用二级抗体(1:5000，zb-2305或zb-2301，中山金桥生物技术)孵育1h。通过增强化学发光法(ecl)获得图像。
[0032]
图1为血清外泌体的鉴定结果；a：血清来源外泌体中标记蛋白cd9和tsg101的蛋白质印迹；b：使用zetasizer纳米zs90粒度仪测量的外泌体直径；c：从血清中纯化的外泌体的电子显微照片(比例尺500nm)。
[0033]
从图1a可以看出，tsg101和cd9在血清样本中表达，从而确定了外泌体的存在。
[0034]
从图1b可以看出，从血清中提取的外泌体的直径小于100nm。
[0035]
从图1c可以看出，外泌体的形态通过tem进一步观察，证明从血清样本中提取了外泌体。
[0036]
第二步、从外泌体中提取rna、高通量测序，以及分析已知mirna和新鉴定的mirna的表达水平
[0037]
根据mircute plus mirna qpcr试剂盒(tiangen)提取总rna。rna的浓度和纯度由nanodrop2000(thermo scientific,usa)测定，检测od260/od280的比值为1.96。od260代表核酸的吸光度，od280代表蛋白质的吸光度，od260/od280的比值在1.8-2.0之间认为rna较纯，rna样品被认为是高质量样品。当比值低于1.8时，考虑有蛋白质或酚污染。当比值高于2.0时，表明可能有异硫氰酸残存。
[0038]
使用高通量测序技术进行小rna测序。通过高通量测序获得原始图像数据文件，并通过碱基判读技术用将其转换为原始序列数据。使用cutadapt软件从原始序列读数中删除适配器序列，bowtie软件去除rrna、trna、miscrna、snrna和snorna；在mirbase 21.0数据库中(http://mirbase.org/)再次搜索过滤的序列以确认新的mirna。
[0039]
使用degseq软件分析已知mirna和新鉴定的mirna的表达水平。计算出的mirna表达指数为百万分之一(tpm)。将p值设置为《0.05，将折叠变化(fc)设置为≤0.5或≥2，以识别布鲁氏菌病患者和健康对照者之间血清外泌体中的差异表达mirnas(dems)。
[0040]
共在3名布鲁氏菌病患者和3名对照者的血清外泌体中鉴定出2006种mirna。两组之间共有41个成熟的mirna存在显著差异(p《0.05)。布鲁氏菌病组和对照组之间存在不同的mirna表达模式。其中，27个mirna下调，14个mirna上调。
[0041]
图2为差异表达的血清外泌体mirna的热图。热图显示了布鲁氏菌病和对照血清外泌体样品中mirna的差异表达(p《0.05)。与对照组相比，布鲁氏菌病患者样本中27个mirna表达下调，14个表达上调。与对照组相比，布鲁氏菌病患者中hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p的外泌体表达显著增加。
[0042]
第三步、qrt-pcr
[0043]
从25名布鲁氏菌病患者和25名健康对照者中提取血清外泌体，并进一步提取外泌体的总rna。使用mircute plus mirna第一链cdna试剂盒(kr211，tiangen)将总rna逆转录成cdna。使用mircute plus mirna qpcr试剂盒(fp411，tiangen)进行定量rt-pcr。所有mirna引物均购自上海sangon。u6被设定为内部控制，ct值被归一化。
[0044]
本实验基于测序结果，选择了前5个上调最大的dems进行验证，并在表2中列出了实验中用到的mirna引物。
[0045]
表2研究中使用的microrna引物
199a-5p；c：hsa-mir-148a-5p；d：hsa-mir-11400联合hsa-mir-199a-5p；e：hsa-mir-11400联合hsa-mir-148a-5p；f：hsa-mir-148a-5p联合hsa-mir-199a-5p；g：hsa-mir-11400联合hsa-mir-148a-5p和hsa-mir-199a-5p。
[0056]
从表5和图4可以得出，每一个单独的mirna都可以将布鲁氏菌病感染与健康对照区分开来：
[0057]
hsa-mir-11400的曲线下面积auc＝0.79(95％ cl:0.66
–
0.91)(p《0.05)；
[0058]
hsa-ir-199a-5p的曲线下面积auc＝0.81(95％ cl:0.69
–
0.92)(p《0.05)；
[0059]
hsa-mir-148a-5p的曲线下面积auc＝0.74(95％ cl:0.60
–
0.88)(p《0.05)。
[0060]
此外，mirna的组合通常会提高auc：
[0061]
组合hsa-mir-11400和hsa-mir-199a-5p的auc为0.90(95％cl:0.80
–
0.98)(p《0.05)；
[0062]
组合hsa-mir-138a-5p与hsa-mir-1240或hsa-mir-19a-5p，auc＞0.8(p《0.05)；
[0063]
同时组合hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p的auc为0.92(95％cl:0.84
–
1.00)(p《0.05)。
[0064]
表5.对血清外泌体单独和联合治疗布鲁氏菌病的预测价值
[0065][0066][0067]
第五步、游靶基因的预测
[0068]
通过mirwalk数据库预测mirna的靶基因。通过比较三个数据库的结果获得mirna的最终基因：targetscan、mirdb和mitarbase。
[0069]
预测了465个dems的靶基因，其中包括25个免疫相关的靶基因，这些免疫基因主要富集于癌症中的蛋白聚糖、nf-kappa b信号通路和il-17信号通路。
[0070]
通过dems-免疫基因网络构建，发现免疫基因可能主要受到hsa-mir-199a-5p和hsa-mir-148a-5p的调节。在这些免疫基因中，plxna2、il17ra、prkca、cd22、acvr1b和cbl的表达与gse69597数据集中一致。
[0071]
mirwalk预测了总共465个hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p的靶基因。此外，表6中列出了这些上调dems的靶基因数量。
[0072]
表6显著上调的dems的潜在目标基因
[0073]
上调的dems数量hsa-mir-1140026hsa-mir-199a-5p174hsa-mir-148a-5p265
靶基因总数量465
[0074]
第六步、筛选免疫基因，并对25个免疫基因进行分析
[0075]
使用immport数据库检索mirna下游靶基因中的免疫相关基因集。
[0076]
通过imimport数据库，本发明从dems的靶基因中筛选了25个免疫相关基因。
[0077]
为了解推定靶基因的生物学功能和涉及的信号通路，使用go(www.geneontology.org)和kegg(www.genome.jp/kegg)数据库对每个基因进行注释。
[0078]
使用“clusterprofiler”包进行go富集分析的结果如图5a所示。生物过程(bp)分析表明，免疫基因在造血、细胞连接组织、细胞-细胞连接组织等的正调节中显著富集；细胞成分(cc)显示靶免疫基因在膜区、膜微区和血清膜受体复合物中富集；分子功能(mf)分析表明，靶基因在受体配体活性和细胞因子受体活性方面显著丰富。
[0079]
此外，kegg通路分析表明，dem的免疫基因在癌症蛋白聚糖、nf-kappa b信号通路和il-17信号通路等中显著富集(图5b)。
[0080]
图5为基因本体论(go)和京都基因和基因组百科全书(kegg)对血清衍生的外泌体mirnas(hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p)的免疫相关基因的分析；a：对明显富集的靶基因的go分析；在生物过程(bp)、细胞成分(cc)和分子功能(mf)术语中，列出了前5个显著富集的项目；p值、富集因子和mirna靶点数量被用来确定go的富集程度；b：通过kegg进行mirna靶点分析的途径富集度分析。allenrichment分数被评估为p值的负对数(-log p)。
[0081]
从geo数据库下载的gse69597 mrna数据集来分析免疫基因的表达水平，数据集在gpl16791平台(illumina hiseq 2500)中获得，数据集内包含了46例布鲁氏菌病和46例健康者全血。确定布鲁氏菌病和对照全血之间差异表达免疫基因的students't检验。
[0082]
最后，有意义的免疫基因需要满足以下标准：首先，dems与其对应下游的靶基因呈相反关系，若dems上调，其靶基因则下调，其次，p《0.05有统计学差异。通过大数据筛选出的差异基因dems就是目标外泌体，通过目标外泌体基因预测下游的靶基因。一般来说，目标基因和其下游对应的靶基因是呈相反的关系，即目标基因上调，其对应下游的靶基因是下调的。
[0083]
为了更好地研究布鲁氏菌病血清外泌体中这些dems的分子免疫机制，使用cytoscape软件构建了dems免疫基因网络(图6)。图6为mirna-免疫相关基因调控网络。方块代表基因，椭圆代表mirnas。
[0084]
上调的dems与免疫基因之间的相互作用如下：hsa-mir-11400与两个免疫基因cd48和ptgdr相互作用(主要是研究目标基因对应下游哪些靶基因。比如hsa-mir-11400对应下游的靶基因有cd48和ptgdr，有助于后期筛选通路)。hsa-mir-148a-5p与10个免疫基因(tnf、cyld、plxna2、fgf5、pthlh、csf2rb、gnrhr、lepr、rora和cbl)相互作用；hsa-mir-199a-5p与13个免疫基因(nfya、traf3、kng1、nos1、actg1、tafa2、cdh1、cd22、kl、tgfb2、acvr1b、il-17ra和prkca)相互作用。
[0085]
图6为mirna-免疫相关基因调控网络。方块代表基因，椭圆代表mirnas。
[0086]
使用gse69597数据集鉴定了18种免疫基因(cd48、tnf、ptgdr、cyld、plxna2、gnrhr、lepr、rora、cbl、nfya、traf3、actg1、cdh1、cd22、tgfb2、acvr1b、il-17ra和prkca)的外周血水平。
[0087]
图7为根据gse69597数据确定18个免疫基因的mrna表达；a：nfya；b：cyld；c：traf3；d：cdh1；e：plxna2；f：lepr；g：tgfb2；h：gnrhr；i：il-17ra；j：rora；k：cd48；l：prkca；m：tnf；n：actg1；o：ptgdr；p：cd22；q：acvr1b；r：cbl。数据以平均值
±
sd表示。p＞0.05，*p＜0.05和**p＜0.01。对于上调的dems(hsa-mir-11400、hsa-mir-199a-5p和hsa-mir-148a-5p)，与正常对照组相比，布鲁氏菌病患者外周血中只有plxna2、il17ra、prkca、cd22、acvr1b和cbl的表达持续降低，其余表达增加或差异无统计学意义。因此，hsa-mir-199a-5p-cd22、hsa-mir-9a-5p-il17ra、hsa-ir-199a-5b-prkca、hsa-mir-199a-5p-acvr1bhsa-mir-148a-5p-plxna2和hsa-mir-138a-5p-cbl被鉴定为布鲁氏菌病的六种潜在调节途径。
[0088]
综上，本实施例从受试者的外周血中分离出血浆外泌体。通过透射电镜(tem)、纳米粒追踪分析(nta)和蛋白质印迹(wb)鉴定外显体的形态和蛋白质标记，证实了外泌体mirna在布鲁氏菌病中具有诊断潜力，为布鲁氏菌病提供潜在的诊断生物标志物，通过构建了潜在的外泌体mirna-mrna调控网络，为布鲁氏菌病的发病机制和治疗提供了新的思路。
[0089]
实施例2布鲁氏菌病的外泌体mirna标志物的检测试剂盒
[0090]
1、试剂盒组份
[0091]
表7试剂盒组份
[0092]
管号试剂名称主要组份规格1hsa-mir-11400引物引物seq id no:410-25μl
×
1管2hsa-mir-199a-5p引物引物seq id no:610-25μl
×
1管3hsa-mir-148a-5p引物引物seq id no:810-25μl
×
1管4miliq-water超纯水1ml
×
1管
[0093]
2、试剂盒使用方法
[0094]
(1)提取样本rna：提取待检样本中的rna；(2)基因扩增：对待检样本中的rna进行逆转录和dna扩增；(3)结果检测：对扩增结果进行检测。
[0095]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

技术特征：
1.布鲁氏菌病的外泌体mirna标志物，其特征在于，包括hsa-mir-11400、hsa-mir-199a-5p或hsa-mir-148a-5p中的一种或两种以上；hsa-mir-11400的核苷酸序列为seq id no:3，hsa-mir-199a-5p的核苷酸序列为seq id no:5，hsa-mir-148a-5p的核苷酸序列为seq id no:7。2.用于检测权利要求1所述的布鲁氏菌病的外泌体mirna标志物的试剂在制备检测试剂盒中的应用。

技术总结
本发明属于生物检测技术领域，具体公开了布鲁氏菌病的外泌体miRNA标志物，包括hsa-miR-11400、hsa-miR-199a-5p或hsa-miR-148a-5p中的一种或两种以上；hsa-miR-11400的核苷酸序列为SEQ ID NO:3，hsa-miR-199a-5p的核苷酸序列为SEQ ID NO:5，hsa-miR-148a-5p的核苷酸序列为SEQ ID NO:7。本发明还公开了上述用于检测上述布鲁氏菌病的外泌体miRNA标志物的试剂在制备检测试剂盒中的应用，以及用于检测布鲁氏菌病的外泌体miRNA标志物的试剂在制备检测试剂盒中的应用。本发明公开的在布鲁氏菌病中具有诊断潜力，为布鲁氏菌病提供潜在的诊断生物标志物。断生物标志物。断生物标志物。


技术研发人员：李智伟 陈志强 王倩 李书灵 丁剑冰 张峰波
受保护的技术使用者：新疆维吾尔自治区人民医院
技术研发日：2023.04.21
技术公布日：2023/8/22