1.本发明涉及生物
技术领域：
：，特别是涉及一种用于检测或辅助检测质子射线辐射的生物标志物、试剂盒及其检测方法和应用。
背景技术：
：：2.随着我国载人航天和空间实验室建设的快速发展，登月和深空探测计划的实施将推动我国尖端科学技术的发展和国防实力的提升，还将奠定中国在未来世界新资源开发与利用、月球资源控制、太空探索与研究格局中的地位和作用。飞船交互对接、出舱活动、建立空间站、长期载人航天任务、载人登月等工程实施中，航天员在航天器舱内尤其是舱室外活动将面临暴露于复杂的空间辐射环境的健康危害。虽然空间辐射环境极为复杂，但质子仍是空间辐射中的主要成分且在空间辐射中占有重要地位。3.质子辐射能够使受作用物质的分子或原子产生电离，作用于人体或生物体时，可引起组织细胞中原子或分子的变化，使受照细胞被杀伤或发生变异，从而导致对人体的各种健康危害。然而现有技术中缺乏简单、快速的质子辐射的检测方法。因此，发展一种快速、安全监测航天员受照剂量的检测技术就显得越来越迫切。技术实现要素：4.本发明的目的在于提供一种用于检测或辅助检测质子射线辐射的微生物菌群生物标志物，通过检测所述微生物菌群生物标志物dna表达量或相对表达量，能准确检测到待测对象是否遭受质子射线辐射。5.为解决上述技术问题，本发明提供了以下技术方案：本发明提供了一种用于检测或辅助检测质子射线辐射的微生物菌群生物标志物，所述微生物菌群包括蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter中的至少一种。6.本发明提供了上述微生物菌群作为生物标志物在制备检测或辅助检测是否遭受质子射线辐射的产品中的应用。7.优选的，通过检测样本中所述微生物菌群dna的表达量或相对表达量以确定是否遭受质子射线辐射。8.更优选的，所述样本为粪便。9.本发明提供了一种试剂盒，包括：检测所述微生物菌群dna表达量或相对表达量的试剂。10.优选的，所述试剂包括检测所述微生物菌群dna表达量或相对表达量的引物或探针。11.更优选的，所述引物为引物组1、引物组2、引物组3和引物组4中的任意一组或几组；所述引物组1中，上游引物序列如seqidno.1所示，下游引物序列如seqidno.2所示；所述引物组2中，上游引物序列如seqidno.3所示，下游引物序列如seqidno.4所示；所述引物组3中，上游引物序列如seqidno.5所示，下游引物序列如seqidno.6所示；所述引物组4中，上游引物序列如seqidno.7所示，下游引物序列如seqidno.8所示。12.本发明还提供了上述的试剂盒在检测或辅助检测是否遭受质子射线辐射中的应用。13.本发明还提供了一种检测或辅助检测是否遭受质子射线辐射的方法，所述方法包括：提取样本中上述的微生物菌群的总dna进行pcr扩增，所得的pcr扩增产物经测序、数据分析获得所述微生物菌群dna的表达量或相对表达量，确定样本是否遭受质子射线辐射。14.优选的，遭受质子射线辐射的样本中微生物菌群dna表达量或相对表达量大于正常微生物菌群dna的表达量或相对表达量。15.本发明提供了一种用于检测或辅助检测质子射线辐射的微生物菌群生物标志物，对balb/c小鼠和c57bl/6小鼠辐射后，可检测到蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter的表达量增加，对待测对象的蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter表达量或相对表达量与未遭受质子射线辐射的balb/c小鼠和c57bl/6小鼠进行比较，能检测到待测对象是否遭受质子射线辐射，及时进行事故救援、伤员处理和临床治疗应用，对辐射处理具有重要的意义。附图说明16.图1为balb/c小鼠和c57bl/6小鼠质子照射与未照射后收集粪便进行的菌群主成分pcoa分析。17.图2为balb/c小鼠和c57bl/6小鼠质子照射与未照射后收集粪便进行的菌群在门一级的成分分析。18.图3为balb/c小鼠和c57bl/6小鼠质子照射与未照射后收集粪便进行的菌门中显著性差异分析；其中，a为balb/c小鼠质子照射后，表达量会发生变化的菌门；b为c57bl/6小鼠小鼠质子照射后表达量升高的菌门。19.图4为balb/c小鼠和c57bl/6小鼠质子照射与未照射后收集粪便进行的菌属中显著性差异分析；其中，a为balb/c小鼠质子照射后，表达量发生变化的菌属；b为c57bl/6小鼠小鼠质子照射后表达量发生变化的菌属。具体实施方式20.本发明提供了一种用于检测或辅助检测质子射线辐射的微生物菌群生物标志物，所述微生物菌群包括蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter中的至少一种。本发明中，所述蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter的dna表达量或相对表达量在遭受质子射线辐射的样本中升高。21.本发明提供了上述微生物菌群作为生物标志物在制备检测或辅助检测是否遭受质子射线辐射的产品中的应用。本发明中，所述产品通过检测样本中所述微生物菌群dna的表达量或相对表达量以确定是否遭受质子射线辐射。所述样本优选为待测样品的粪便。在本发明的具体实施例中，所述微生物菌群dna的表达量或相对表达量优选为c57bl/6小鼠和balb/c小鼠结直肠粪便中微生物菌群dna的表达量或相对表达量。22.本发明提供了一种试剂盒，包括：检测所述微生物菌群dna表达量或相对表达量的试剂。本发明中，所述试剂优选的包括检测所述微生物菌群dna表达量或相对表达量的引物或探针。所述引物为引物组1、引物组2、引物组3和引物组4中的任意一组或几组；所述引物组1中，上游引物序列如seqidno.1所示，下游引物序列如seqidno.2所示；所述引物组2中，上游引物序列如seqidno.3所示，下游引物序列如seqidno.4所示；所述引物组3中，上游引物序列如seqidno.5所示，下游引物序列如seqidno.6所示；所述引物组4中，上游引物序列如seqidno.7所示，下游引物序列如seqidno.8所示。本发明中，通过利用所述引物扩增微生物菌群dna不同的区域，检测该区域的基因表达量或相对表达量。所述引物扩增的基因片段如下表1所示：本发明还提供了上述的试剂盒在检测或辅助检测是否遭受质子射线辐射中的应用。本发明中，所述应用优选为利用上述试剂盒对所述微生物菌群dna进行扩增，确定所述微生物菌群dna的表达量或相对表达量，从而确定待测样品是否遭受质子射线辐射。23.本发明还提供了一种检测或辅助检测是否遭受质子射线辐射的方法，所述方法包括：提取样本中上述的微生物菌群的总dna进行pcr扩增，所得的pcr扩增产物经测序、数据分析获得所述微生物菌群dna的表达量或相对表达量，确定样本是否遭受质子射线辐射。24.本发明中，提取得到微生物菌群的总dna后优选的还包括对dna的提取质量进行检测；所述检测的方法优选为琼脂糖凝胶电泳。本发明中，对dna的提取质量进行检测后选优的包括对dna进行定量；所述定量的方法优选为紫外分光光度计。本发明对所述提取样本中微生物菌群总dna的方法并没有特殊限定，在本发明的具体实施例中，所述方法优选为ctab法。本发明对所述测序的平台并没有特殊限定，在本发明的具体实施例中，所述测序平台优选为novaseq6000测序仪。25.本发明中，所述数据分析能够用于将来自待测样品的所述微生物菌群dna表达量或相对表达量转换为所述待测样品的检测结果。本发明中，所述数据分析包括数据处理模块、多样性分析模块、物种注释并差异分析和高级分析模块；其中，所述数据拆分模块用于将对测序获得的双端数据进行拆分，并去除接头和barcode序列，再通过数据的拼接和过滤，dada2降噪；多样性分析模块用于将待测对象的所述微生物菌群dna表达量或相对表达量值与未遭受质子射线辐射的待测对象的微生物菌群dna表达量或相对表达量值进行比较；物种注释并差异分析和高级分析模块用于输出结论。本发明中，所述数据分析所用的软件和版本具体如表2所示。26.表2软件版本及分析要求分析内容软件名称版本序列拼接flash(http://ccb.jhu.edu/software/flash/)v1.2.8低质量序列过滤fqtrim(http://ccb.jhu.edu/software/fqtrim/)v0.94嵌合体过滤vsearch(https://github.com/torognes/vsearch)v2.3.4去噪（dada2）qiime2(https://qiime2.org/)2019.7去除测序接头cutadapt(http://cutadapt.readthedocs.org/en/stable/guide.html)v1.9去除低质量序列fqtrim(http://ccb.jhu.edu/software/fqtrim/)v0.94pcoa联川生物云平台（https://www.omicstudio.cn/tool/24）v8.4.0柱状图联川生物云平台（https://www.omicstudio.cn/tool/17）v8.4.0箱线图联川生物云平台（https://www.omicstudio.cn/tool/1）v8.4.0本发明中，若待测样品的微生物菌群dna表达量或相对表达量大于未遭受质子射线辐射待测样品的微生物菌群dna表达量或相对表达量，则结论是待测样品遭受了质子射线辐射；若待测样品的微生物菌群dna表达量或相对表达量小于或等于未遭受质子射线辐射待测样品的微生物菌群dna表达量或相对表达量，则结论是待测样品未遭受质子射线辐射。27.本发明中，所用原料、试剂与设备均为已知产品，采用常规市售产品即可。28.为了进一步说明本发明，下面结合实施例对本发明提供的技术方案进行详细地描述，但不能将它们理解为对本发明保护范围的限定。29.实施例1利用微生物菌群生物标志物蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter中的至少一种检测样本是否遭受质子射线辐射的方法：1.实验流程1.1微生物组总dna提取利用ctab法对小鼠粪便样本进行微生物组总dna的提取，并通过琼脂糖凝胶电泳检测dna提取质量，同时采用紫外分光光度计对dna进行定量。30.1.2pcr扩增利用表1中的引物组对步骤1.1中获得的dna进行扩增，其中，pcr反应的体系和条件分别如表3和表4所示。31.表3pcr反应体系pcr反应成分pcr反应体积phusionhotstartflex2xmastermix12.5微升forwardprimer2.5微升reverseprimer2.5微升templatedna50纳克addddh2oto25微升表4pcr反应条件pcr反应温度pcr反应时间循环数98℃30秒98℃10秒54℃30秒35循环72℃45秒72℃10分钟4℃∞1.3pcr产物定量步骤1.2得到的pcr产物利用ampurextbeads(beckmancoultergenomics,danvers,ma,usa)纯化，然后通过qubit(invitrogen,usa)定量，再按照样本的测序量要求，进行相应比例的混合。32.1.4扩增产物回收纯化pcr扩增产物通过2%琼脂糖凝胶电泳进行检测，回收采用ampurextbeads回收试剂盒。33.1.5扩增产物定量混样上机测序对纯化后的pcr产物使用agilent2100生物分析仪(agilent，美国)和illumina(kapabiosciences,woburn,ma，美国)的文库定量试剂盒进行评估，合格的文库浓度应在2nm以上。将合格的各上机测序文库（index序列不可重复）梯度稀释后，根据所需测序量按相应比例混合，并经naoh变性为单链进行上机测序；使用novaseq6000测序仪进行2×250bp的双端测序，相应试剂为novaseq6000spreagentkit(500cycles)。34.实施例2生物信息学分析流程1.数据拆分对测序获得的双端数据，首先需要根据barcode信息对样品进行数据拆分，并去除接头和barcode序列。35.2.数据拼接和过滤具体步骤有：1)去除rawdata的引物序列和平衡碱基序列；其中，使用的软件为cutadapt(v1.9)，参数：'-gr1-gr2-n1-o17-m100'；2)将每一对paired-endreads根据overlap区拼接合并成一条更长的tag；其中，使用的软件为flash（v1.2.8），参数：'-m10ꢀ‑m100ꢀ‑x0.25ꢀ‑t1ꢀ‑z'；3)对测序reads进行窗口法质量扫描，扫描窗口默认为100bp，当窗口内平均质量值低于20时，将read从窗口起始到3’终止的部分截掉；其中，使用的软件为fqtrim，参数：'-p33ꢀ‑w100ꢀ‑q20ꢀ‑l100ꢀ‑m5ꢀ‑p1ꢀ‑vꢀ‑otrim.fastq.gz'；4)去除截短后长度小于100bp的序列；5)去除截短后n（不确定模糊碱基）的含量在5%以上的序列；6)去除嵌合体序列；其中，使用的软件为vsearch（v2.3.4），参数默认；3.dada2去噪通过qiimedada2denoise-paired调用dada2进行长度过滤和去噪，获得asv(feature）特征序列和asv(feature）丰度表格，并去除singletonsasvs(即在全体样本中，序列总数仅为1的asv(feature），默认操作)。36.4.多样性分析基于得到的asv(feature）特征序列和asv(feature）丰度表格进行alpha多样性分析和beta多样性分析。其中，alpha多样性分析主要通过observed_species，shannon，simpson，chao1，goods_coverage，pielou_e七大指数对生境内的多样性进行评估。beta多样性分析主要通过计算四种距离（weighted_unifrac，unweighted_unifrac，jaccard，bray_curtis）采用六大分析对生境间（样本/分组间）的多样性进行评估分析。37.5.物种注释根据asv(feature）序列文件采用silva（release138，https://www.arbsilva.de/documentation/release138/）数据库以nt-16s数据库进行物种注释，并根据asv(feature）丰度表对各物种在样本中的丰度进行统计，注释的置信度阈值为：0.76.差异分析和高级分析基于得到的物种丰度统计信息，进行各比较组之间的差异分析。根据样本情况选择不同的统计学方法：采用fisher'sexacttest，对无生物学重复样品差异比较；mann-whitneyutest，对有生物学重复两组样品差异比较；kruskal-wallistest，对有生物学重复样品的多组之间比较，筛选阈值：p《0.05。38.实施例31.选择6-8周龄雄性小鼠，照射质子射线前在spf级动物房饲养1周以适应环境，给予足量的水和饲料，光照时间为12：12小时，动物伦理符合军事医学研究院辐射医学研究所相关规定。其中，小鼠为balb/c小鼠和c57bl/6小鼠，购自于斯贝福公司。39.将小鼠放进装置中依次进行5gy的全身照射，待3天后，在无菌条件下取小鼠粪便5颗左右放入冻存管后，立即用液氮速冻处理1小时，-80℃保存，待用。其中，质子辐射源来自中国原子能科学研究院100mev质子回旋加速器产生的90mev单能质子。40.2.利用实施例1所述方法检测小鼠在照射质子射线前后的dna表达量或相对表达量。41.3.经联川生物云平台（https://www.omicstudio.cn/tool/24）（版本v8.4.0）得到图1，经联川生物云平台（https://www.omicstudio.cn/tool/17）（版本v8.4.0）得到图2，由联川生物云平台（https://www.omicstudio.cn/tool/1）（版本v8.4.0）得到图3-4。42.可以看出，图1中c57bl/6小鼠质子照射后菌群与未照射组有明显的分离，证明差异显著，balb/c小鼠质子照射后菌群与未照射组有分离的趋势。图2中发生变化的菌群排名前十的是：拟杆菌门bacteroidetes，厚壁菌门firmicutes，变形菌门proteobacteria，疣微菌门verrucomicrobia，脱铁杆菌门deferribacteres，蓝藻门cyanobacteria，未分类unclassified，ε变形菌epsilonbacteraeota，放线菌门actinobacteria和patescibacteria。43.图3-4表明蓝藻门cyanobacteria和ε变形菌epsilonbacteraeota在balb/c和c57bl/6小鼠肠道菌群中均升高，anaerovorax和变形菌门螺杆菌属helicobacter在balb/c和c57bl/6小鼠肠道菌群中均升高，厚壁菌门梭菌属clostridium和firmicutes_unclassified在balb/c和c57bl/6小鼠肠道菌群中均降低。44.由以上实施例可知，在遭受了质子射线辐射的小鼠样本中，微生物菌群生物标志物蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacterdna表达量或相对表达量大于未遭受质子射线辐射的微生物菌群生物标志物的dna表达量或相对表达量。表明微生物菌群蓝藻门cyanobacteria、ε变形菌epsilonbacteraeota、anaerovorax、螺杆菌属helicobacter能够作为生物标志物用于检测或辅助检测样本是否遭受质子射线辐射，对辐射处理具有重要的意义。45.以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的
技术领域：
：，均同理包括在本发明的专利保护范围内。当前第1页12当前第1页12