本研究由美国南加州大学赫尔曼奥斯特罗牙科学院牙周病学、诊断科学与牙科卫生学系的黄宇（音译）、陈凯西（Casey Chen）等人，联合华盛顿大学微生物学系、泰国国王科技大学系统生物学与生物信息学研究组以及巴西圣保罗大学微生物学系的科研人员共同完成。研究成果以“Comparative genomic hybridization and transcriptome analysis with a pan-genome microarray reveal distinctions between JP2 and non-JP2 genotypes of *Aggregatibacter actinomycetemcomitans*”为题，于2013年发表于期刊 *Molecular Oral Microbiology*（第28卷，第1-17页）。

一、 学术背景

本研究的科学领域聚焦于口腔微生物学与病原菌基因组学，具体研究对象是伴放线聚集杆菌（Aggregatibacter actinomycetemcomitans, Aa）。Aa是一种革兰氏阴性兼性厌氧菌，是局限性侵袭性牙周炎等多种牙周病的主要致病菌，但也常存在于健康个体口腔中。该菌种包含不同的血清型和克隆谱系，菌株间的基因组内容存在显著差异（可达19.5%）。其中，血清型b型菌株中，若其白细胞毒素（leukotoxin）操纵子启动子区存在一个530碱基对（bp）的特征性缺失，则被定义为高毒力的JP2基因型克隆。大量流行病学和前瞻性研究表明，JP2克隆与非洲裔人群的侵袭性牙周炎密切相关，其致病风险远高于非JP2克隆。虽然JP2型菌株的高白细胞毒性已被证实，但其高毒力是否仅归因于此单一因素尚不明确。研究者推测，JP2基因型可能拥有一系列独特的、非JP2基因型所不具备的毒力决定簇。

本研究的核心目的在于验证一个假设：JP2与非JP2基因型在基因组内容（即基因的有无）和基因表达模式上存在差异。通过构建一个定制的Aa泛基因组（pan-genome）微阵列，研究者旨在系统性比较两种基因型菌株的基因组构成和转录组特征，以期发现与JP2型高毒力相关的特定基因组岛（genomic islands）或差异表达的基因通路，为深入解析Aa的致病机制提供关键线索。

二、 详细工作流程

本研究包含一系列严谨且环环相扣的实验与分析流程，主要分为五个关键步骤：

1. Aa泛基因组微阵列的设计、构建与验证

   • 目的与对象：开发一个能够覆盖Aa物种遗传多样性的高密度寡核苷酸微阵列。
   • 方法：
     • 数据基础：研究基于此前已测序的14株Aa基因组数据，并结合新测序的4株菌株（S23A， SCC4092， Aas4a， A160）的全基因组数据，共计18株Aa菌株的基因组信息。
     • 基因簇定义：对这18个基因组中的42,668个预测基因进行同源性聚类，最终定义了3426个同源基因簇（homologous gene clusters），这些基因簇共同构成了Aa的泛基因组。
     • 探针设计与芯片制作：从每个基因簇中选择最长的序列作为代表序列，使用eArray和ArrayOligoSelector软件设计长度为60bp的寡核苷酸探针。最终，设计了15,208个探针，代表3126个基因簇，外加536个对照探针，由安捷伦（Agilent）公司合成制作成8x15K规格的定制微阵列。
     • 探针过滤：为避免序列变异导致杂交信号不准确，进行了一系列严格的生物信息学过滤：要求探针必须与目标基因簇内所有成员的序列相似度≥80%，且不能与非目标基因簇匹配；同时，剔除代表基因数少于3个探针的基因簇。最终，用于分析的探针集包含10,934个探针，覆盖2676个基因簇（包括1762个核心基因和914个附属基因）。
     • 验证：使用11株已完成基因组测序的参考菌株（包括4株完整基因组和7株草图基因组）的基因组DNA与微阵列进行杂交，以验证微阵列检测基因存在与否的灵敏度和特异性。验证结果显示，灵敏度（被检出的基因数/已知存在于该基因组的基因数）在0.948–0.989之间，特异性（未被检出的基因数/已知不存在于该基因组的基因数）在0.867–0.986之间，表明该微阵列性能良好，可可靠地用于后续研究。

2. 菌株收集与基因分型

   • 研究对象：研究共纳入20余株Aa临床分离株，重点关注血清型b的菌株。菌株来源多样，包括不同种族背景（高加索人、非洲裔、西班牙裔等）和牙周健康状况（健康、牙龈炎、慢性牙周炎、侵袭性牙周炎）的个体。
   • 方法：通过PCR确定每株菌的血清型。同时，使用特异性引物扩增白细胞毒素操纵子启动子区，根据是否存在530-bp的缺失，将血清型b菌株明确划分为JP2基因型和非JP2基因型。研究中“非JP2”特指同样为血清型b但具有完整长启动子的菌株，以聚焦于JP2特异的基因组特征。

3. 比较基因组杂交（Comparative Genomic Hybridization, CGH）分析基因组内容

   • 目的：比较JP2与非JP2基因型菌株的基因存在/缺失模式，识别基因型特异的基因和基因组岛。
   • 对象与样本量：对12株血清型b菌株（6株JP2型，6株非JP2型）的基因组DNA进行CGH分析。每株菌设置生物学重复。
   • 实验流程：
     • 提取菌株基因组DNA，用Cy3荧光染料标记。
     • 将标记的DNA与定制微阵列杂交。
     • 使用安捷伦扫描仪扫描，并通过Feature Extractor软件提取信号数据。
     • 数据分析：对原始信号进行背景扣除、归一化（除以探针中的A碱基数）和log2转换。同一基因簇的探针信号取平均值。针对每次杂交，根据信号值的双峰分布，选择一个阈值来确定每个基因的“存在”或“缺失”状态。
     • 聚类分析：基于所有菌株的基因存在/缺失二进制数据，使用MEV软件进行层次聚类分析，以可视化菌株间的亲缘关系。
     • PCR验证：对CGH识别出的6个基因型特异基因，通过设计特异性引物进行PCR扩增，以验证CGH结果的可靠性。

4. 转录组（Transcriptome）分析

   • 目的：比较JP2与非JP2基因型菌株在指数生长期的基因表达差异。
   • 对象与样本量：选取3株JP2型菌株（HK1651， D28S-1， D41S-1）和2株非JP2型菌株（SCC1398， ANH9381）进行转录组分析。每株菌设置两个独立的生物学重复培养。
   • 实验流程：
     • 细菌在含0.6%酵母提取物的胰蛋白酶大豆肉汤中培养至对数生长中期。
     • 提取总RNA，确保无DNA污染。
     • 使用MessageAmp II-Bacteria试剂盒，通过体外转录（IVT）方法将RNA扩增并用Cy3标记。
     • 将标记的cRNA与同一套定制微阵列杂交。
     • 数据分析：提取处理后的信号值，对所有阵列数据进行分位数归一化。分析JP2型与非JP2型菌株之间共有基因（1952个）的表达差异。使用t检验（p < 0.05）筛选差异表达基因。利用DOOR数据库预测HK1651菌株的操纵子结构，重点关注那些至少包含两个基因、且在两个基因型间表达差异倍数≥2倍的操纵子。
     • qRT-PCR验证：选取5个差异表达基因（如lkta, tolr, *napd*等），使用同一批RNA样本进行实时定量PCR验证，以确认微阵列结果的准确性。

5. 生物信息学分析

   • 目的：深入解读CGH和转录组分析结果，定位基因型特异的基因组岛，分析差异表达操纵子的启动子结构。
   • 方法：
     • 基因组岛鉴定：结合已发表的全基因组序列数据（如HK1651和ANH9381），将CGH发现的基因型特异基因定位到具体的基因组坐标上，确定它们是否位于已知或新的基因组岛内，并分析这些岛的插入位点及对周边基因的影响。
     • 启动子序列比对：对差异表达的操纵子（除已知的leukotoxin操纵子外），在JP2代表株HK1651和非JP2代表株ANH9381之间进行上游启动子区域的序列比对，寻找可能导致表达差异的序列变异。

三、 主要研究结果

1. 菌株基因组内容聚类结果：基于2676个基因存在/缺失信息的层次聚类分析显示，Aa菌株清晰地分为两大分支。JP2基因型的6株血清型b菌株形成了一个独立且紧密的簇，与非JP2基因型的血清型b菌株以及血清型c菌株明显分开。这表明JP2克隆在基因组构成上具有高度保守性和独特性。

2. 基因型特异的基因与基因组岛：CGH分析共鉴定出21个在JP2型与非JP2型之间存在与否完全不同的基因。

   • JP2特异基因与基因组岛：发现9个仅存在于JP2型菌株的基因。其中，三个基因（P-cluster03231， 03273， 03132）在JP2代表株HK1651中位于一个23，394 bp的基因组岛IAAI-1上。该岛插入在一个截短的苹果酸脱氢酶基因和一个截短的黄素蛋白同源基因之间。而在非JP2菌株（如ANH9381）的相同基因座，则是完整的苹果酸脱氢酶和黄素蛋白基因。另外两个JP2特异基因（P-cluster03717， 03907）位于另一个7，244 bp的基因组岛IAAI-3上。
   • 非JP2特异基因与基因组岛：发现12个仅存在于非JP2型菌株的基因。其中四个基因位于非JP2代表株ANH9381的一个基因组岛IANH-1上。特别值得注意的是，在HK1651和ANH9381的同一个基因座（*guaA*和*aa01895*之间），分别插入了不同的基因组岛（IAAI-3和IANH-2），表明该位点是基因组岛插入的“热点”。此外，在另一株非JP2菌株I23c的相同位点，发现了第三个不同的基因组岛。这种“相同位点，不同岛屿”的现象强烈暗示了水平基因转移在此处的高度活跃性。PCR验证了部分基因型特异的基因存在模式。

3. 转录组差异分析结果：转录组分析鉴定出150个基因在JP2型与非JP2型菌株间的表达水平存在显著差异（p < 0.05）。通过操纵子水平分析，进一步聚焦到5个具有显著表达差异的操纵子（每个操纵子至少有两个基因表达差异≥2倍）：

   • Leukotoxin操纵子：如预期所料，JP2型菌株中该操纵子的表达水平显著高于非JP2型，这与JP2型启动子缺失导致的高转录活性一致。
   • nap操纵子：编码膜结合硝酸还原酶系统（membrane-bound nitrate reductase system）。该操纵子所有基因在JP2型菌株中的表达水平极低，而在非JP2型中表达很高，差异倍数最高可达近30倍。
   • tol-pal操纵子：编码革兰氏阴性菌中维持细胞包膜完整性的Tol-Pal系统。该操纵子在JP2型菌株中的表达水平显著高于非JP2型。
   • 一个包含*trkA*等基因的操纵子：在JP2型中表达更高。
   • sgab/tkt操纵子：在JP2型中表达较低。
   • 启动子分析：对上述除leukotoxin外的四个差异表达操纵子进行启动子序列比对，未发现JP2型与非JP2型之间存在明显的启动子结构差异（除sgab/tkt操纵子有一个非JP2特异的插入片段外）。这提示观察到的表达差异可能源于菌株间不同的基因调控机制，而非启动子序列本身。
   • qRT-PCR验证：对5个选定的差异表达基因进行qRT-PCR验证，结果与微阵列数据具有良好的一致性（三个基因相关性极佳，另外两个趋势一致），证实了转录组分析结果的可靠性。

四、 结论与意义

本研究首次系统地揭示了伴放线聚集杆菌JP2与非JP2基因型在全基因组内容和对数生长期基因表达谱上的显著区别。主要结论如下： 1. 基因组差异：JP2克隆是一个基因组内容高度保守且独特的群体。研究鉴定出多个JP2特异和非JP2特异的基因组岛，并发现同一基因组位点可作为不同基因组岛插入的“热点”，这为了解Aa的基因组可塑性和进化（特别是水平基因转移的作用）提供了新视角。 2. 表达谱差异：除了已知高表达的leukotoxin外，JP2型菌株还表现出其他关键基因通路的差异化调节。其中，nap操纵子（硝酸还原酶）在非JP2型中高表达而在JP2型中几乎沉默，以及tol-pal操纵子在JP2型中高表达，是两个最突出的发现。这些差异可能与菌株在不同微环境（如牙周袋）中的适应性、能量代谢、毒力或与宿主免疫系统的相互作用有关。 3. 研究价值： * 科学价值：本研究构建并验证的Aa泛基因组微阵列，是研究该物种遗传多样性和基因表达的有力工具。研究结果将JP2型高毒力的研究焦点，从单一的白细胞毒素，扩展到了特定的基因组获得/缺失事件和全局性的基因表达重编程。这为理解细菌克隆分化与毒力进化的关系提供了范例。 * 应用价值：鉴定出的基因型特异基因组岛和差异表达操纵子（如nap和tol-pal），为设计假设驱动的实验指明了方向。未来可以通过基因敲除、回补、功能分析等手段，直接验证这些遗传和表达差异是否以及在多大程度上贡献了JP2克隆的高致病性，从而可能揭示新的药物靶点或诊断标志物。

五、 研究亮点

1. 方法学创新：成功设计并验证了首个基于18株菌全基因组信息的Aa泛基因组微阵列，实现了对高遗传多样性物种进行高效、经济的基因组比较和转录组分析。
2. 系统性发现：首次在同一研究中整合了CGH和转录组分析，全面描绘了JP2与非JP2基因型在基因存在与否和基因表达水平两个层面的差异，提供了更完整的生物学图景。
3. 关键生物学发现：
   • 明确了JP2克隆在基因组聚类上的独立性。
   • 发现了基因组岛插入“热点” 这一有趣现象，为研究细菌基因组进化机制提供了具体案例。
   • 鉴定出nap和tol-pal这两个与JP2毒力可能高度相关的新候选通路，尤其是nap操纵子在两种基因型间的极端表达差异，是一个极具潜力的研究方向。
4. 清晰的科学逻辑：研究从假设出发，通过定制工具、系统比较、多角度验证（PCR， qRT-PCR）、深入生物信息学分析，最终将结果聚焦到少数但可能至关重要的生物学差异上，逻辑严谨，结论可靠。