轮状病毒A组 G9P[8]
轮状病毒A组G9P[8]基因型由VP7糖蛋白G9与VP4蛋白酶敏感蛋白P[8]组合而成，拥有经典的Wa样全基因组骨架（G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1）。其VP7基因至少可分化为四个与人类疾病相关的系统发育谱系：lineage I（1980年代早期株）、lineage II（印度新生儿株116E，即Rotavac疫苗株）、lineage III（1990年代中期以来的全球优势谱系）及lineage VI（中国t203样株，2010年后快速扩散）。G9P[8]是继G1–G4之后第五个具有全球流行病学重要性的轮状病毒基因型，已在亚洲、欧洲、非洲和美洲多国成为优势流行株。其VP7基因持续积累抗原表位区氨基酸替换，近年来DS-1样骨架重配株的出现进一步丰富了该基因型的遗传多样性。
发现与早期鉴定
G9血清型轮状病毒的发现始于1983年2月。当时费城儿童医院从一名18月龄急性胃肠炎患儿粪便中分离出一株与已知四种血清型（G1–G4）均无交叉中和反应的病毒，命名为WI61。Clark等人于1987年正式将其描述为「推定的人类新血清型」，这标志着G9作为第五个具有全球流行病学意义的轮状病毒血清型登台。后续回顾性研究在华盛顿特区1980年的存档粪便样本中检出了更早的G9P[8]病毒，将G9的流行时间线至少前推了三年。
WI61的VP7基因与当时已知的G1–G4血清型序列同源性不足85%，其RNA电泳型呈「长型」特征。随后日本（F45，1985年）和印度（116E，1986年从无症状新生儿中分离，属G9P[11]）也相继报告了G9毒株。
全球传播与基因型格局
G9在1980年代仅零星检出，但自1990年代中期开始在各大洲迅速扩散。1995–1996流行季，美国费城地区约50%的轮状病毒腹泻由G9引起，这是该血清型在西半球首次大规模暴发。此后G9P[8]在欧洲（匈牙利、比利时、意大利）、亚洲（泰国、中国、日本、印度）、非洲和南美洲持续被报告，并在多个国家成为优势基因型。
2012年后，G9P[8]取代G3P[8]成为中国大陆最主要的轮状病毒基因型。Liu等（2022）对2016–2018年中国115株G9P[8]的全基因组分析显示，绝大多数毒株拥有Wa样骨架：G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1。其中2018年有48株出现NSP4基因由E1向E2（DS-1样）的替换，提示同一基因型内正在发生骨架层面的重配。2020年泰国报告了完全DS-1样骨架的G9P[8]重配株（G9-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2），进一步拓展了该基因型的遗传多样性。
VP7谱系多样性
G9的VP7基因至少可分为六个系统发育谱系（lineage I–VI），其中与人类疾病密切相关的有四个。
- Lineage I：以WI61、F45、AU32为代表，均为1980年代分离的早期毒株。
- Lineage II：唯一代表为印度新生儿株116E（G9P[11]），即Rotavac疫苗株。该谱系VP7在氨基酸水平上与lineage I和III存在多处差异，尤其在与抗体结合相关的表面暴露位点上。
- Lineage III：自1990年代中期起成为全球优势谱系，涵盖了绝大多数当代G9P[8]流行株。与lineage I毒株相比，lineage III的VP7在7-1a、7-1b和7-2中和表位区域积累了多个氨基酸替换。
- Lineage VI：以中国t203样毒株为代表，2010年前后在北京重新出现后迅速扩散，并在2012年取代lineage III成为中国大部分地区的优势G9亚型。该谱系的VP7与lineage III存在四个特征性氨基酸替换及多个同义密码子替换。
贝叶斯进化分析估算G9 VP7基因的最晚共同祖先可追溯至约1935年，核苷酸进化速率约为1.38×10⁻³替换/位点/年。
关键氨基酸突变
G9P[8]流行株的VP7基因在与现有疫苗株的比较中暴露出若干值得关注的氨基酸差异，尤其在三个已鉴定的中和表位（7-1a、7-1b、7-2）区域。
| 基因 | 突变 | 类型 | 功能影响 |
|---|---|---|---|
| VP7 | N69-T70-S71（糖基化基序） | N-糖基化位点插入 | 中国G9-III及G9-VI流行株中普遍存在该N-X-T/S基序；Rotavac疫苗株116E（lineage II）及Rotasill疫苗株中缺失。该基序可在VP7表面添加N-聚糖侧链，空间遮蔽邻近中和表位，构成潜在的免疫逃逸机制（Mao et al., 2022）。 |
| VP7 | K143M | 氨基酸替换（表位7-2） | 北印度G9流行株中检出K143M替换；116E疫苗株为K143。M143位于VP7三聚体表面突出的7-2表位环区，侧链性质的改变可能削弱疫苗诱导抗体对流行株的识别（Gupta et al., 2021）。 |
| VP7 | D100N（lineage III）/ D100G（lineage II） | 氨基酸替换（表位7-1a） | Lineage I毒株（WI61等）为D100，lineage II（116E）为G100，lineage III全球流行株为N100。该位点位于7-1a表位核心区域，不同谱系的残基差异直接影响单克隆抗体中和效价（Hoshino et al., 2004; Cuffia et al., 2016）。 |
| VP7 | T87I | 氨基酸替换 | Lineage II（116E）特征性替换；lineage I及lineage III流行株为T87。I87位于VP7表面β-折叠环，靠近7-1a表位和抗体结合区域，可能影响VP7的抗原构象（Rippinger et al., 2010）。 |
疫苗交叉保护
当前全球广泛使用的两种口服轮状病毒疫苗对G9P[8]的保护机制不同。五价RotaTeq含五种人-牛重配株，分别为G1、G2、G3、G4及P1A[8]（后者携带牛G6 VP7与人P[8] VP4）。该疫苗不含G9 VP7组分，其对G9P[8]的保护来自P1A[8]重配株——因G9P[8]流行株共享P[8] VP4骨架，P1A[8]组分可诱导针对P[8]的交叉中和抗体。FDA于2017年基于临床试验数据将G9正式纳入RotaTeq的适应证。
单价Rotarix（G1P[8]）对G9P[8]提供异型保护，其诱导的交叉中和抗体可部分覆盖G9 VP7表位。在印度获批的Rotavac以116E毒株（G9P[11]，VP7属lineage II）为基础，是全球唯一直接靶向G9 VP7的轮状病毒疫苗。然而lineage II的VP7与当前主流流行株（lineage III/VI）在数个中和表位上存在差异。俄罗斯的研究显示当地G9-III流行株与Rotavac疫苗株在VP7表位区域有0–3个氨基酸差异；中国G9流行株中普遍存在的N69-T70-S71糖基化位点在116E中缺失。
真实世界数据显示RotaTeq在低儿童死亡率国家对轮状病毒胃肠炎的总疫苗有效性约为85%，对G9P[8]的有效性与对G1–G4相当。一项尼加拉瓜的病例对照研究未发现RotaTeq引入后特定毒株的持续替代。
临床与公共卫生影响
G9P[8]是全球五种主要轮状病毒基因型（G1P[8]、G2P[4]、G3P[8]、G4P[8]、G9P[8]）之一，合计占人类轮状病毒感染的约90%–95%。在疫苗引入前，G9的区域流行率在7%至16%之间波动。2016年全球约有128,500名5岁以下儿童死于轮状病毒腹泻，G9在亚洲部分地区（如中国、泰国、印度）的检出率已超过G1和G3。
G9P[8]呈现出两种值得警惕的进化趋势。一是VP7抗原表位的持续漂移——lineage VI在中国快速取代lineage III便是一个信号。二是骨架重配的活跃——DS-1样G9P[8]重配株在泰国和马来西亚的检出表明，该基因型可能正在突破Wa样骨架的限制。两种趋势叠加对抗原匹配性提出了更高要求。在轮状病毒疫苗覆盖率持续提升的背景下，对G9P[8]的分子流行病学监测不可松懈。
参考文献
- Clark HF, Hoshino Y, Bell LM, Groff J, Hess G, Bachman P, Offit PA. Rotavirus isolate WI61 representing a presumptive new human serotype. J Clin Microbiol. 1987;25(9):1757-1762. DOI: 10.1128/jcm.25.9.1757-1762.1987
- Hoshino Y, Jones RW, Ross J, Honma S, Santos N, Gentsch JR, Kapikian AZ. Rotavirus serotype G9 strains belonging to VP7 gene phylogenetic sequence lineage 1 may be more suitable for serotype G9 vaccine candidates than those belonging to lineage 2 or 3. J Virol. 2004;78(14):7795-7802. DOI: 10.1128/JVI.78.14.7795-7802.2004
- Laird AR, Gentsch JR, Nakagomi T, Nakagomi O, Glass RI. Characterization of serotype G9 rotavirus strains isolated in the United States and India from 1993 to 2001. J Clin Microbiol. 2003;41(7):3100-3111. DOI: 10.1128/jcm.41.7.3100-3111.2003
- Liu X, Wang M, Li S, Li J, Xiao J, Li H, Zhang Q, Kong X, Wang H, Li D, Duan Z. Genomic and evolutionary characteristics of G9P[8], the dominant group A rotavirus in China (2016–2018). Front Microbiol. 2022;13:997957. DOI: 10.3389/fmicb.2022.997957
- Mao T, Wang M, Wang J, Ma Y, Liu X, Wang M, Sun X, Li L, Li H, Zhang Q, Li D, Duan Z. Phylogenetic analysis of the viral proteins VP4/VP7 of circulating human rotavirus strains in China from 2016 to 2019 and comparison of their antigenic epitopes with those of vaccine strains. Front Cell Infect Microbiol. 2022;12:927490. DOI: 10.3389/fcimb.2022.927490
- Vesikari T, Matson DO, Dennehy P, Van Damme P, Santosham M, Rodriguez Z, et al. Safety and efficacy of a pentavalent human-bovine (WC3) reassortant rotavirus vaccine. N Engl J Med. 2006;354(1):23-33. DOI: 10.1056/NEJMoa052664
- Bhandari N, Rongsen-Chandola T, Bavdekar A, John J, Antony K, Taneja S, et al. Efficacy of a monovalent human-bovine (116E) rotavirus vaccine in Indian infants: a randomised, double-blind, placebo-controlled trial. Lancet. 2014;383(9935):2136-2143. DOI: 10.1016/S0140-6736(13)62630-6