TABLE 1

Genes most distinctly mutated between carriage- and disease-associated isolates

LocusProductDiseaseaCarriageaP valueb
GBSCOH1_RS00285Phosphoribosylformylglycinamidine synthase19 (11/8)9 (6/3)6.6900 × 10−4
GBSCOH1_RS01850cSerine/threonine protein kinase Stk117 (13/4)6 (3/3)4.7100 × 10−4
GBSCOH1_RS04925c23S rRNA methyltransferase12 (7/5)1 (1/0)0.0115
GBSCOH1_RS04975ABC transporter permease18 (15/3)4 (1/3)7.1600 × 10−7
GBSCOH1_RS05095Cell division protein FtsK28 (23/5)19 (12/7)4.5700 × 10−6
GBSCOH1_RS06600Cell wall anchor Srr218 (10/8)7 (6/1)2.3500 × 10−4
GBSCOH1_RS07650cAmidase14 (9/5)4 (2/2)0.0161
GBSCOH1_RS07705cTwo-component sensor histidine kinase CovS15 (11/4)5 (5/0)0.0088
GBSCOH1_RS08845cDNA polymerase III subunit alpha PolC17 (10/7)6 (4/2)4.7100 × 10−4
GBSCOH1_RS04480Type II CRISPR RNA-guided endonuclease Cas914 (8/6)19 (13/6)6.4900 × 10−5
GBSCOH1_RS05130Carbamoyl phosphate synthase large subunit10 (8/2)17 (15/2)4.5100 × 10−5
GBSCOH1_RS06095Peptidase C59 (7/2)21 (13/8)2.8400 × 10−12
GBSCOH1_RS08245X-prolyl-dipeptidyl aminopeptidase5 (3/2)14 (11/3)1.2200 × 10−4
GBSCOH1_RS08960cHypothetical protein7 (5/2)14 (12/2)2.3800 × 10−3
  • a Number of mutations exclusively acquired by strains associated with each clinical state. Values in parentheses correspond to the numbers of nonsynonymous and synonymous substitutions detected, respectively.

  • b Bonferroni-adjusted P values were obtained with the outlier test.

  • c Gene with a mutation frequency significantly biased toward either carriage or disease, independent of association with clade I or II (Fig. 1).