Branches corresponding to partitions reproduced in less than 50% of bootstrap replicates were collapsed. The
MP tree was obtained using the Close-Neighbor-Interchange algorithm [17] with search level 3 [16, 17] in which the initial trees were obtained with the random addition of sequences (10 replicates). The tree is drawn to scale with branch lengths calculated by the average pathway method [17] and with the number of changes over the whole sequence as units. Estimates of Average Evolutionary Divergence over Sequence Pairs of stkP within penicillin susceptibility groups The number of amino acid and of nucleotide substitutions per site was averaged over all sequence pairs within each group by the Poisson correction JPH203 molecular weight method and the Maximum Composite Likelihood method, respectively, using
MEGA version 4 software [14]. Standard error estimates were obtained by the bootstrap procedure (1000 replicates). StkP modelling A 3D-model of the selleck chemicals kinase domain of the StkP protein (271 residues long) of strain R6 was obtained using the sequence (accession number NP_359169). BLASTP analysis indicated that the serine-threonine kinase Selumetinib from strain R6 has 63% sequence identity with serine-threonine kinase of Mycobacterium tuberculosis (PDB ID: 1o6yA). The following structure PDB ID: 1o6yA; 1mruA.pdb, 1mruB.pdb, 1y8gB.pdb and 1zmwB.pdb were used as a template for building a homology model for the kinase domain of StkP with the SWISS-MODEL server [18, 19]. Ramachandran plot analysis for phi and psi torsion angles indicated that 95.9% of residues were in the allowed region of IMP dehydrogenase the plot, which is
more than the average cut-off of 90% used in most reliable models [20]. The final alignment adjustments and visualisation were undertaken with Deep View/Swiss-PdbViewer version 3.7. Genotyping of pbp genes Genetic polymorphism of penA, pbpX and pbp1A genes (encoding PBP2B, PBP2X and PBP1A, respectively) of all clinical strains was investigated first by restriction fragment length polymorphism (RFLP) analysis. A number was given to each restriction pattern for each of the three pbp genes analysed, so the PBP profile has three numbers (for example: 4-9-7). The full genes were amplified by PCR using the primers described in Table 2 and 0.8 U of iProof Polymerase (Bio-Rad, Hercules, California) according to the manufacturer’s instructions, with 35 cycles at an annealing temperature of 56°C for 30 seconds. The amplification products of penA and pbpX were digested for 1 H with 5 U of both HaeIII and RsaI restriction endonucleases. The amplification product of pbp1A was similarly digested with HaeIII and DdeI (all restriction enzymes supplied by New England Biolabs, Beverly, Mas.). The digested products were separated on agarose gel. Dice coefficient of similarity was used for cluster analysis with the unweighted pair group method with arithmetic averages using BioNumerics software v3.5 (Applied Maths, Sint-Martens-Latem, Belgium). The position tolerance was set to 1.