使用短读和长读测序技术生成尿细菌完整基因组的杂交从 头 基因组组装

作为机构审查委员会批准的研究19MR0011（UTD）和STU 032016-006（UTSW）的一部分，从同意妇女收集的尿液中培养细菌。 1. 改良增强型尿培养 注意：所有培养步骤必须在无菌条件下进行。对所有器械、溶液和培养基进行灭菌。用70%乙醇清洁工作区域，然后设置本生燃烧器，并在靠近火焰的地方小心地工作，以减少污染的机会。或者，可以使用II类生物安全柜来维持无菌环境。穿戴适当的个人防护装备（PPE），以避免接触潜在的致病微生物。 接种甘油储备尿液和集落分离 在室温（RT）下解冻储存甘油的尿液。解冻后，涡旋样品5秒以混合。在无菌微量离心管中，在无菌1x磷酸盐缓冲盐水（PBS）中制备1：3和1：30稀释的尿液，最终体积为100μL。注意：将500μL未稀释的尿液和500μL50%无菌甘油混合在冷冻管中并储存在-80°C下来制备甘油储备尿液。 使用前在37°C下预热琼脂平板15分钟。请参阅 图1， 了解适合常见尿细菌属的培养基类型和培养条件。在接种前通过移液将稀释的尿液混合均匀，将稀释的尿液平板100μL放在所需的琼脂平板上，并使用无菌玻璃珠传播样品。将100μL的1x PBS稀释剂的板放在单独的板上作为无生长对照。注意：如果试图培养常见的泌尿病理性物种（例如， 大肠杆菌，克雷伯氏菌属，粪肠球菌 等），建议使用显色琼脂（材料表），因为它可以轻松识别泌尿病源性细菌物种（图1）。粘菌素纳利迪酸 （CNA） 或 MRS 琼脂有助于从已知含有革兰氏阴性尿路病原体的尿液中分离出挑剔的革兰氏阳性种属（例如 乳酸杆菌属），这些菌株在非选择性琼脂中可能优于挑剔的种类。 将板倒置在35°C的所需大气条件下孵育24小时，用于泌尿道病原体，3-5天用于挑剔细菌（图1）。 孵育期过后，从培养箱中取出板。从每个板中，选择表现出独特颜色，形态或溶血模式的菌落。 使用无菌环将细菌菌落重新划到相应的琼脂上，并将平板倒置在所需气氛中孵育2-5天，以获得良好分离的菌落。注意：如果使用BAP进行原代培养，则在显色琼脂上修补菌落可能提供有关样品中细菌群异质性的有用信息。 在液体肉汤和甘油放养细菌分离物中培养 一旦获得与母体菌落形态相匹配的分离菌落，选择单个菌落并使用无菌接种回路接种到3mL液体肉汤中。参见 图1， 了解能够支持普通尿液微生物群属生长的肉汤。用石膏膜密封琼脂板，并将其储存在4°C下2-4天。将液体培养物在所需的大气条件下孵育1-5天，直到培养物明显浑浊。 观察生长后，涡旋培养物，然后将1mL过夜培养物加入500μL无菌50%甘油在2mL冷冻管中;密封并通过倒置轻轻混合。为每个菌落准备两个甘油储备（一个用作备用）并储存在-80°C。 2. 通过16S rRNA基因桑格测序鉴定细菌种类 注意：微生物身份可以使用基质辅助激光解吸电离飞行时间质谱（MALDI-TOF）20进行确认。 菌落-聚合酶链反应 （PCR） 通过加入12.5μL的2x Taq聚合酶预混液，0.5μL的10μM 8F引物，0.5μL的10μM 1492R引物（材料表）和11.5μL无核酸酶水21，在PCR管中制备25μLPCR反应。注意：如果对多个样品进行PCR，请制作Taq聚合酶混合物，引物和无菌无核酸酶水的反应预混液。然后将25μL等分到每个PCR管中。 要进行集落-PCR，请使用无菌牙签或移液器吸头从重新条纹中擦拭分离良好的集落。将菌落重悬于步骤2.1.1制备的PCR反应混合物中。轻轻混合。通过以2000 x g的快速旋转收集管底部的液体。注：确保样品没有气泡。包括单独含有PCR反应混合物的无模板对照（NTC）样品。 将样品管放入热循环器中，并运行以下程序：95°C3分钟;40次循环：95°C持续30秒，51°C持续30秒，72°C持续1分钟30秒;72°C持续10分钟;保持在10°C。 凝胶提取和物种鉴定 完成PCR运行后，在0.5x Tris-硼酸盐-EDTA（TBE）缓冲液中制备的1%琼脂糖凝胶上检查PCR产物。在浇注凝胶之前，加入溴化乙锭（EtBr）。然后，使用梳子将凝胶浇铸到至少20μL样品体积的孔中。注意：EtBr是一种怀疑致癌的插层剂。处理时始终戴上手套和PPE，并根据机构的指导方针处理含有EtBr的材料。 凝胶凝固后，将凝胶放入充满0.5x TBE缓冲液的电泳槽中，然后取出梳子。将第一孔中的1 kb分子量标准品和10-20μL的PCR反应加载到后续孔中。以 100-140 V 电压运行，直至解决。在紫外光下观察凝胶，并确认在~1.5 kb处存在NTC孔中不存在的明确定义的条带。注意：紫外线对皮肤和眼睛有害，在观察凝胶时使用适当的防护装置，并穿戴适当的个人防护用品。注意：对于某些细菌，集落 PCR 可能不成功。从分离的gDNA进行PCR是替代选择22。 使用剃须刀切除~1.5 kb的条带，并将凝胶插条转移到干净的微量离心管中。按照制造商的说明（材料表）继续进行凝胶提取方案。通过微量分光光度计测量纯化DNA的浓度。注：需要>10 ng/μL的浓度，1.7-2.0之间的A260/280是可以接受的。 根据任何所选 Sanger 测序服务的指南，为每个样品准备两个 Sanger 测序反应，一个使用 8F，另一个使用 1492R 引物在无核酸酶水中。 收到测序数据后，将DNA序列上传到NCBI基本局部比对搜索工具（BLAST）网站（blast.ncbi.nlm.nih.gov/Blast.cgi），选择核苷酸BLAST（blastn），选择rRNA / ITS数据库16S核糖体RNA序列（细菌和古菌），然后运行Megablast程序。隔离点可以通过对数据库中引用的最高质量命中来识别。注意：一些细菌物种在其16S rRNA序列中表现出高同一性，并且仅通过这种方法可能无法区分。物种形成将需要DNA同源性和生化分析，以自信地区分同一属23的成员。 3. 基因组DNA（gDNA）的提取 注：本节利用 材料表中 引用的gDNA提取试剂盒中提供的试剂和离心柱，从不同细菌物种中高产提取优质基因组DNA。下面提供了建议的修改和说明。 根据制造商的说明准备试剂盒试剂盒。 将来自良好分离的菌落的细菌接种到培养基中，并在图1中注明的温度和大气压下孵育直到观察到足够的生长，在适当的无菌肉汤（图1）中制备 3-10mL 培养物。 孵育后，使用分光光度计24测量培养物在600nm（OD 600）处的光密度。 通过以1：10的比例稀释过夜培养物来制备样品以进行定量。包括无菌培养基的空白，以便测量。通过从样品读数中减去空白读数并乘以稀释因子10来计算光密度。 使用OD600 测量和预先建立的OD600 与CFU / mL的物种比率，计算获得2 x 109 个细胞所需的培养物数量。 将所需的培养体积在5000×g下离心5分钟至沉淀。 吸出上清液并将沉淀重悬于200μL冷TE缓冲液中（在程序开始时在冰上预冷）。 将样品在5000×g下离心2分钟。 除去上清液，然后将沉淀重悬于180μL酶解缓冲液（ELB）中，并加入20μL预煮的RNase A（10mg / mL）。为了有效裂解革兰氏阳性细菌，加入18μL突变溶素（25 kU / mL）。涡旋均匀，然后将样品在37°C的旋转器上孵育2小时。注意：建议将制造商方案中描述的ELB用于革兰氏阳性和革兰氏阴性细菌。 按照制造商的说明继续。注意：如果需要，重复洗脱步骤一次或两次，以获得额外的gDNA产量。 按照第4节中的说明评估提取的gDNA的质量，如果将在1周内使用，则将gDNA储存在4°C。或者，将gDNA保持在-20°C以长期储存。 4. 评估提取的gDNA的质量 为了通过凝胶电泳评估质量，如第2.2小节所述制备1%琼脂糖凝胶。在干净的管中制备样品：在副膜上混合1-2μL提取的gDNA和3μL的2x上样染料。上样后运行凝胶，然后在紫外光下可视化。注意：成功的gDNA提取将通过凝胶顶部的离散条带和最小的涂片来证明（图2A）。拖尾表明剪切。如果没有明显的 gDNA 波段和/或大量涂片，则重复 gDNA 提取。考虑减少 RNase A 和 Proteinase K 的孵育时间。如果观察到1.5-3 kb左右的两个条带，这表明RNA污染（图2B）。准备新鲜的RNase A并重复提取。 要使用微量分光光度计评估质量，请使用微量分光光度计测量gDNA浓度和吸光度比A260/280。浓度>50 ng/μL，A260/280在1.7-2.0之间是可以接受的。注意：低gDNA产量可能是由于低输入，高输入，核酸酶污染，裂解不足。吸光度高于该范围表示RNA污染。如果 gDNA 质量较差，则重复提取。 要通过荧光计评估质量，请按照制造商的说明使用高灵敏度测定试剂盒和荧光计仪器（材料表）量化gDNA浓度。需要>50纳克/μL的浓度。 5. 配对端的下一代短读测序和文库制备 注：短读测序可以在不同的读取长度和方向上对各种仪器进行。建议对细菌 WGS 进行 150 bp（300 周期）配对测序。文库制备和测序均可外包给核心设施或商业实验室。 根据制造商的说明（材料表）准备测序库。遵循制造商推荐的最终加载库浓度;但是，建议的修改是在1.8 pM下加载池化库，以便在NextSeq仪器上生成最佳读取。 虽然是可选的，但使用生物分析仪（材料表）来评估合并文库片段分布，并确保片段大小平均为600 bp。 6. 纳米孔 MinION 测序文库制备 根据制造商的协议（材料表）准备测序库。使用两个条形码扩展试剂盒，可在单个流通池上多路复用多达 24 个样品。建议分两部分进行文库制备，一次12个样品，同时多路复用24个样品。所有24个样本可以如下所述合并。注意：样品可以在完成原生条形码连接后在4°C下储存过夜 – 如有必要，这可在方案中提供停止点。在文库制备方案的天然条形码连接部分结束时，建议将每个样品的等摩尔量汇集到最大DNA质量（ng）可能。 为此，请根据制造商的说明使用荧光计（材料表）在条形码连接后量化所有样品。估计具有最低dsDNA浓度的样品的体积，然后计算该样品中发现的总dsDNA。使用此数字可确定将汇集在一起的所有其他样品的等摩尔量。注意：由于等摩尔计算将最大化池化dsDNA的量，从而产生高容量池（>65μL），因此需要清理以浓缩池。 dsDNA池清理和浓缩 将2.5倍体积的顺磁性珠（材料表）加入DNA池中，然后轻轻轻地轻拂管子以混合内容物。将试管置于室温下旋转器5分钟，以2000×g旋转样品并在磁铁上沉淀。 加入250μL新鲜制备的70%乙醇（在无核酸酶的水中），注意不要干扰沉淀。吸出乙醇并重复乙醇洗涤一次。 第二次吸气后，以2000×g旋转样品并将其放回磁铁上。移走任何残留的乙醇，让样品干燥约30秒。 从磁铁中取出管并将沉淀重悬于60-70μL无核酸酶水中。在室温下孵育2分钟。将样品沉淀在磁铁上，直到洗脱物清晰，然后取出洗脱并转移到干净的1.5mL微量离心管中。 使用荧光计定量浓缩池，然后制备等分试样以进入适配器连接步骤：在65μL最终体积中制备700ng样品。将池的其余部分保留在4°C，以便在第一次运行完成后完成第二次运行。 按照制造商的指示继续进行适配器连接，并将样品加载到流通池上。启动排序运行。注意：在上样前，从流通池底液端口吸入空气和~200 μL储存缓冲液。这对于流池引物的成功和样品上样至关重要。通过流通池的底漆口抽取和浇注溶液时，使用p1000移液器和吸头。 根据制造商的说明对库进行排序。 打开用于排序的操作软件，然后单击”开始”。输入实验的名称，建议的命名法包括运行日期和用户名。单击”继续选择试剂盒”，选择使用的相应库制备试剂盒和条形码扩展包，然后单击”继续运行选项”。 如果计划为第二次运行准备足够的库，请将运行长度调整为 48 小时（否则，请保留默认的 72 小时）。单击”继续基本呼叫”。 检查基本调用选项 “配置：快速基本调用 “，并确保将 “条形码” 设置为 “已启用”， 以便输出 FASTQ 文件将修剪条形码序列，并根据条形码解复用到单独的目录中。单击 继续 输出。 选择保存输出排序数据的位置。如果仅保存 FASTQ 输出，则预计大约需要 30-50 Gb 的数据，如果还保存 FAST5 输出，则>500 Gb 的数据。取消选中筛选选项 Qscore： 7 |Readlength：如果未筛选 ，如果计划继续进行第 7.2 节中描述的过滤，否则请保留选中状态并将 Readlength 调整为 200。 单击” 继续运行安装程序”并查看所有设置。如果设置正确，请单击” 开始”，否则单击” 上一步 “并进行任何必要的调整。 如果需要，可以按照制造商的说明清洗流通池，并重新装入剩余的池中。对剩余的池重复 6.2 中的步骤，一旦第一次运行完成并且流通池已洗涤。注：设置第二次运行时，根据制造商对先前在48小时内运行中使用的流通池的建议，将偏置电压调整为-250 mV。 7. 评估和准备阅读材料 注意：推荐的目录结构如图 4 所示。在继续下面的计算步骤之前，创建 在桌面中找到的目录，即Long_Reads，Short_Reads和Trimmed_Reads。 短读 （图 3)注意：短读以 FASTQ 格式生成。这些文件包含每个 FASTQ 的最大读取次数 4000 次。这些通常被压缩（.gz存档）并组织成多个文件。根据平台的不同，条形码通常会被修剪。某些程序接受压缩格式的文件，其他程序可能需要在导入之前提取它们。读数必须通过质量控制（QC）步骤，以确保基因组组装过程中的数据准确性。如果CLC基因组学工作台不可用，则可以使用替代程序来修剪和QC短读，例如Trimomatic25 或 Trim Galore （https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/） 用于修整，FastQC （http://www.bioinformatics.babraham.ac.uk/projects/fastqc/） 用于评估读取质量。建议将平均短读覆盖率（通过将读取次数乘以平均读取长度并除以基因组大小来估计）>100x。打开基因组学工作台软件（材料表）并导入所有配对端的短读FASTQ文件。配对的文件将自动生成。 在CLC_Data下创建一个新文件夹，方法是单击顶部工具栏上的 “新建 “，然后选择 “文件夹…”以存储文件。根据需要命名文件夹，建议的约定是使用示例 ID。将以下步骤的所有输出保存到此文件夹中。 在顶部工具栏上，单击” 导入 “按钮，然后选择 “Illumina…导航到与示例对应的所有短读文件并选定。确保选中配对读取选项，然后取消选中 删除失败的读取 选项。单击” 下一步”，选择” 保存”，然后再次单击” 下一步 “。选择将导入的文件保存在上一步中创建的新文件夹中，然后单击 完成。 为隔离创建所有配对文件的序列列表;这会将读取的数据连接到单个文件中，以便于分析。 在顶部工具栏上，单击”新建”按钮，然后选择”序列列表…”在左侧的目录列表中，选择要连接的文件，然后使用箭头将这些文件移动到右侧的选定文件列表中。单击”下一步”，选择”保存”，然后再次单击”下一步”。选择保存序列列表，然后单击”完成”。 生成序列列表后，立即使用示例 ID 对其进行重命名。 在序列列表上运行 QC 测序读取 工具：此过程将评估短读 NGS 生成的读取的整体质量参数。 在工具箱菜单（左下角窗口）中搜索 “读取排序”工具的 QC。 双击该工具，然后选择要分析的序列列表，然后单击” 下一步”。 确保选中所有输出选项，然后选择结果处理下保存。单击”下一步”并指定以保存输出文件，然后单击”完成”。 在序列列表上运行 “修剪读取” 工具：将根据质量、长度和模糊性进行修剪。此过程假定在此步骤之前已修剪排序中使用的条形码。 在工具箱（左下角窗口）中搜索 “修剪读取” 工具。双击 修剪 读取，然后选择要分析的序列列表，然后单击 下一步。 质量修剪：将质量评分限制设置为0.01，并将模糊的核苷酸保留在2。单击”下一步”。注：参数可由用户自行决定调整;这些是建议的设置。 取消选中 自动直读适配器修整 （仅当适配器在导入到 CLC 之前已从读取中修整时才执行此操作）。单击” 下一步 “并选中 “放弃长度以下的读取”，使用默认值 15。 单击”下一步”，选中”创建报告”，然后选择”保存”。单击”下一步”并指定保存输出文件的位置。单击”完成”。 导出修剪的序列列表：后续的混合组装和分析将在 CLC 之外完成，并且需要导出修剪的短读文件。 从左上角的目录导航中，选择在步骤 7.1.4 中生成的已修剪文件，然后单击顶部工具栏上的” 导出 “。选择 Fastq 作为导出文件类型，然后单击 “下一步”。选中 将配对序列列表 导出到两个文件。然后，单击” 下一步 “并选择要将文件导出到的Trimmed_Reads目录。单击 “完成”。确保已修整的短读文件已成功导出为扩展名为 .fastq的两个文件（R1 和 R2）。注意：修剪后的序列列表必须导出为两个文件，CLC 通常将其指定为 R1 和 R2。这一点至关重要，因为下游混合组件需要这样设置短读数据输入。 重命名导出的文件，请不要在文件名中使用空格和特殊字符。为简单起见，建议的格式trimmed_short_file。R1.fastq. 长（小黄人）读取（图3）注意：以下用于准备混合组装的Long（MinION）测序读取的管道利用了由命令行执行的NanoFilt和Nanostat程序26。 在继续之前安装这些工具，并熟悉 UNIX 的基础知识以执行这些命令。建议使用默认终端和 Bash Shell。常见终端命令和用法的课程指南可在 软件木工27中找到。以下说明假定生成的文件将使用条形码命名法（NB01，NB02等）命名，并保存在Long_Reads目录中。或者，可以在设置排序运行时使用 MinKNOW 完成读取筛选。建议平均长读取覆盖率为 >100 倍。建议的平均读取长度为 >2000 bp;因此，所需的长读取数低于短读取数。 在Long_Reads目录中为运行中使用的每个条形码（条形码01，条形码02等）创建新目录（图4）。将与每个条形码对应的所有 .fastq 文件复制到相应的文件夹中。合并每次运行的每个条形码的所有 .fastq 文件。 打开 终端 并使用 cd 命令导航到Long_Reads目录中的条形码目录 ：cd Desktop/Long_Reads/barcode01 通过执行以下命令，将每个条形码的所有 .fastq 文件连接成一个 .fastq 文件 ：cat *.fastq > NB01.fastq注意：此命令将每个 FASTQ 文件中的所有读取内容合并到一个名为 NB01.fastq 的大型单 FASTQ 中。 使用 NanoStat 通过执行以下命令来评估样本的读取质量 ：NanoStat –fastq NB01.fastq 通过将输出复制到文本或 Word 文件中以供将来参考来记录结果。 使用 NanoFilt 通过执行以下命令来过滤 MinION 读取丢弃的读数，Q < 7 且长度< 200：NanoFilt -q 7 -l 200 bp NB01.fastq | gzip > NB01 _trimmed.fastq.gz 通过执行命令在步骤 7.2.6 中生成的修剪过的文件上运行NanoStat：NanoStat –fastq NB01 _trimmed.fastq.gz 通过将输出复制到文本或Word文件中来记录结果，并与步骤7.2.4中的结果进行比较，以确保筛选成功（表1）。 对排序运行中使用的每个条形码重复步骤 7.2.2 到 7.2.8。注意：在步骤 7.2.6 中生成的 NB01_trimmed.fastq.gz 文件将用于混合程序集。 8. 生成杂交基因组组装 注：以下装配管线利用独轮车19、28、29、30组合了第7.1节和第7.2节中准备的短读和长读（图3）。安装独轮车及其依赖项并执行以下命令。在步骤 7.1.5 中导出的短读文件假定trimmed_short_file命名。R1.fastq和trimmed_short_file。R2.fastq 用于简化。 将短读文件和长读文件组织到名为 Trimmed_Reads 的目录中。该目录必须包含以下内容： 用于修剪长读取 的 .fastq.gz 文件（在步骤 7.2.6 中生成）。 两个 .fastq 文件（R1 和 R2），用于修整的短读（在步骤 7.1.5 中生成）。 导航到使用终端中的 cd 命令存储已读文件的目录 Trimmed_Reads：cd 桌面/Trimmed_Reads 进入正确的目录后，通过执行以下命令，压缩两个短读取文件，使它们也采用 .fastq.gz 格式 ：gzip trimmed_short_file。R1.快餐 对 R1 和 R2 重复步骤 8.2。检查所有读取的文件现在是否都采用 .fastq.gz 格式，并验证所有文件是否都与同一隔离点匹配。 通过运行以下命令，使用独轮车开始混合程序集：独轮车-1 trimmed_short_file。 R1.fastq.gz -2 trimmed_short_file。R2.fastq.gz -l NB01 _trimmed.fastq.gz -o unicycler_output_directory注意：-o 指定将保存独轮车输出的目录，执行命令后，独轮车将创建此目录;不要事先生成目录。运行时间因所用计算机的计算能力以及基因组大小和读取次数而异。这可能需要4小时到1或2天不等。该协议是在具有250 Gb RAM的CentOS Linux 7机器，具有2.5 GHz 12个实用内核和48个虚拟内核的Intel Xeon （R）CPU上执行的。或者，具有 16 Gb RAM 和 2.6 GHz 6 核处理器的个人计算机可以在更长的处理时间内计算这些程序集。 运行完成后，查看独轮车.log文件以确保没有错误 – 记录生成的重叠群的数量、大小和状态（完整、不完整）。 如果标识了不完整的重叠群（在独轮车日志中表示为不完整），请通过在步骤 8.4 中的命令中添加以下标志，以粗体模式重新运行独轮车：–mode bold。注意：粗体模式将降低装配过程中接受的长读桥接的质量阈值;这可能会产生完整的装配，但装配质量可能会降低。建议仅在必要时使用粗体模式，并作为重叠群连接的初步证据，以便稍后通过PCR确认。 9. 评估装配质量 注意：以下协议使用Bandage31 和QUAST32，这两个程序必须在使用前设置（图2 和 图4）。下载后不需要安装绷带，QUAST需要熟悉基本的命令行用法。还建议使用基准通用单拷贝正交学（BUSCO）33评估基因组完整性。 绷带： 点击 文件.然后，选择” 荷载图 “并选择保存到由 Unicycler 在步骤 8.4 中生成的unicycler_output_directory的 assembly.gfa 文件。加载后，单击左侧工具栏上的” 绘制图形 “按钮，然后查看如何连接和组织重叠群（称为节点），以评估程序集是否完成（图 5）。注：完整的组件由两端连接的单个圆形重叠群表示（图5A，B）。不完整的组件具有多个连接在一起的重叠群或线性的（图5C）。小的线性重叠群可能并不完整，因为它们可能表明线性染色体外元素。覆盖，也称为深度，将在绷带中注明，并代表染色体重叠群的相对丰度，在Unicycler中归一化为1x。 夸斯特 在终端中，使用 cd 命令导航到存储独轮车输出的文件夹 ：cd Desktop/Trimmed_Reads/unicycler_output_directory注意：程序集所在位置的路径中不允许有空格，即，指向 Unicycler 输出的任何目录的名称中都不能包含空格。或者，将 assembly.fasta 文件复制到桌面以便于访问。 通过执行以下命令来运行 QUAST：quast assembly.fasta -o quast_output_directory 查看 QUAST 在输出目录中生成的报告quast_output_directory。 10. 基因组注释 注意：下面的注释管道使用 Prokka34，这是一个命令行工具，必须在使用前安装。或者，通过自动GUI K-Base（材料表）使用Prokka，或通过Web服务器RAST35注释基因组。如果将基因组沉积到NCBI中，它们将使用原核基因组注释管道（PGAP）36自动注释。 使用 cd 命令在终端内导航到存储独轮车输出的文件夹（请参阅步骤 9.2.1）。然后，通过执行以下命令运行 Prokka：prokka –前缀 sample_ID –outdir prokka_output_directory assembly.fasta注意：–前缀将根据指定的sample_ID命名所有输出文件。–outdir将创建一个具有指定名称的输出目录，其中将保存所有Prokka输出文件;不要事先为 Prokka 创建输出目录。 通过打开.tsv表和/或将生成的.gff文件上传到序列分析软件中以可视化和分析注释来查看注释（图 6）。 可以根据感兴趣的遗传因素生成特定类型的注释。建议从基因组流行病学中心（www.genomicepidemiology.org/）Web服务器上的用户友好工具开始进行初步分析37，38，39，40，41。其他用于检测CRISPR-cas系统和前噬菌体的工具是可用的（图3）42，43。 11. 建议的数据民主化做法 如果可能，将所有原始读取数据以及组装的基因组存放在公共存储库中，例如NCBI序列读取存档（SRA）和Genbank。在NCBI沉积过程中，基因组通过PGAP管道自动注释。