Genomic Diversity and Zoonotic Potential of Brucella neotomae.

After reports in 2017 of Brucella neotomae infections among humans in Costa Rica, we sequenced 12 strains isolated from rodents during 1955-1964 from Utah, USA. We observed an exact strain match between the human isolates and 1 Utah isolate. Independent confirmation is required to clarify B. neotomae zoonotic potential.

CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers.

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

Characterization of sixteen Achromobacter xylosoxidans phages from Abidjan, Côte d'Ivoire, isolated on a single clinical strain.

Sixteen bacteriophages of Achromobacter xylosoxidans distributed into four genera have been isolated from sewage water in Abidjan, Côte d'Ivoire, using a single clinical strain, and their genomes have been sequenced. Three podoviruses belonged to the genus Phikmvvirus, and these represent the first A. xylosoxidans phages of this genus. Seven podoviruses, distributed into three groups, belonged to the genus Jwalphavirus. Among the siphoviruses, three revealed similarities to Pseudomonas phage 73 and members of the genus Septimatrevirus, and three were YuA-like phages. The virulence of these phages toward a panel of 10 genetically diverse strains was tested, with the phiKMV-like phages showing the broadest host range.

CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins.

CRISPR (clustered regularly interspaced short palindromic repeats) arrays and their associated (Cas) proteins confer bacteria and archaea adaptive immunity against exogenous mobile genetic elements, such as phages or plasmids. CRISPRCasFinder allows the identification of both CRISPR arrays and Cas proteins. The program includes: (i) an improved CRISPR array detection tool facilitating expert validation based on a rating system, (ii) prediction of CRISPR orientation and (iii) a Cas protein detection and typing tool updated to match the latest classification scheme of these systems. CRISPRCasFinder can either be used online or as a standalone tool compatible with Linux operating system. All third-party software packages employed by the program are freely available. CRISPRCasFinder is available at https://crisprcas.i2bc.paris-saclay.fr.

A carrier state is established in Pseudomonas aeruginosa by phage LeviOr01, a newly isolated ssRNA levivirus.

ssRNA bacteriophages are very abundant but poorly studied, particularly in relation to their effect on bacterial evolution. We isolated a new Pseudomonas aeruginosa levivirus, vB_PaeL_PcyII-10_LeviOr01, from hospital waste water. Its genome comprises 3669 nucleotides and encodes four putative proteins. Following bacterial infection, a carrier state is established in a fraction of the cells, conferring superinfection immunity. Such cells also resist other phages that use type IV pili as a receptor. The carrier population is composed of a mixture of cells producing phage, and susceptible cells that are non-carriers. Carrier cells accumulate phage until they burst, releasing large quantities of virions. The continuous presence of phage favours the emergence of host variants bearing mutations in genes involved in type IV pilus biogenesis, but also in genes affecting lipopolysaccharide (LPS) synthesis. The establishment of a carrier state in which phage particles are continuously released was previously reported for some dsRNA phages, but has not previously been described for a levivirus. The present results highlight the importance of the carrier state, an association that benefits both phages and bacteria and plays a role in bacterial evolution.

Pseudolysogeny and sequential mutations build multiresistance to virulent bacteriophages in Pseudomonas aeruginosa.

Coevolution between bacteriophages (phages) and their prey is the result of mutualistic interactions. Here, we show that pseudolysogeny is a frequent outcome of infection by virulent phages of Pseudomonas aeruginosa and that selection of resistant bacterial mutants is favoured by continuous production of phages. We investigated the frequency and characteristics of P. aeruginosa strain PAO1 variants resisting infection by different combinations of virulent phages belonging to four genera. The frequency of resistant bacteria was 10- 5 for single phage infection and 10- 6 for infections with combinations of two or four phages. The genome of 27 variants was sequenced and the comparison with the genome of the parental PAO1 strain allowed the identification of point mutations or small indels. Four additional variants were characterized by a candidate gene approach. In total, 27 independent mutations were observed affecting 14 genes and a regulatory region. The mutations affected genes involved in biosynthesis of type IV pilus, alginate, LPS and O-antigen. Half of the variants possessed changes in homopolymer tracts responsible for frameshift mutations and these phase variation mutants were shown to be unstable. Eleven double mutants were detected. The presence of free phage DNA was observed in association with exclusion of superinfection in half of the variants and no chromosomal mutation could be found in three of them. Upon further growth of these pseudolysogens, some variants with new chromosomal mutations were recovered, presumably due to continuous evolutionary pressure.

Brucella vulpis sp. nov., isolated from mandibular lymph nodes of red foxes (Vulpes vulpes).

Two slow-growing, Gram-negative, non-motile, non-spore-forming, coccoid bacteria (strains F60T and F965), isolated in Austria from mandibular lymph nodes of two red foxes (Vulpes vulpes), were subjected to a polyphasic taxonomic analysis. In a recent study, both isolates were assigned to the genus Brucella but could not be attributed to any of the existing species. Hence, we have analysed both strains in further detail to determine their exact taxonomic position and genetic relatedness to other members of the genus Brucella. The genome sizes of F60T and F965 were 3 236 779 and 3 237 765 bp, respectively. Each genome consisted of two chromosomes, with a DNA G+C content of 57.2 %. A genome-to-genome distance of >80 %, an average nucleotide identity (ANI) of 97 % and an average amino acid identity (AAI) of 98 % compared with the type species Brucella melitensis confirmed affiliation to the genus. Remarkably, 5 % of the entire genetic information of both strains was of non-Brucella origin, including as-yet uncharacterized bacteriophages and insertion sequences as well as ABC transporters and other genes of metabolic function from various soil-living bacteria. Core-genome-based phylogenetic reconstructions placed the novel species well separated from all hitherto-described species of the genus Brucella, forming a long-branched sister clade to the classical species of Brucella. In summary, based on phenotypic and molecular data, we conclude that strains F60T and F965 are members of a novel species of the genus Brucella, for which the name Brucella vulpis sp. nov. is proposed, with the type strain F60T ( = BCCN 09-2T = DSM 101715T).

Generation of a CRISPR database for Yersinia pseudotuberculosis complex and role of CRISPR-based immunity in conjugation.

The clustered regularly interspaced short palindromic repeat - CRISPR-associated genes (CRISPR-Cas) system is used by bacteria and archaea against invading conjugative plasmids or bacteriophages. Central to this immunity system are genomic CRISPR loci that contain fragments of invading DNA. These are maintained as spacers in the CRISPR loci between direct repeats and the spacer composition in any bacterium reflects its evolutionary history. We analysed the CRISPR locus sequences of 335 Yersinia pseudotuberculosis complex strains. Altogether 1902 different spacer sequences were identified and these were used to generate a database for the spacer sequences. Only ∼10% of the spacer sequences found matching sequences. In addition, surprisingly few spacers were shared by Yersinia pestis and Y. pseudotuberculosis strains. Interestingly, 32 different protospacers were present in the conjugative plasmid pYptb32953. The corresponding spacers were identified from 35 different Y. pseudotuberculosis strains indicating that these strains had encountered pYptb32953 earlier. In conjugation experiments, pYptb32953-specific spacers generally prevented conjugation with spacer-positive and spacer-free strains. However, some strains with one to four spacers were invaded by pYptb32953 and some spacer-free strains were fully resistant. Also some spacer-positive strains were intermediate resistant to conjugation. This suggests that one or more other defence systems are determining conjugation efficiency independent of the CRISPR-Cas system.

Development of a multiple-locus variable-number tandem-repeat typing scheme for genetic fingerprinting of Burkholderia cenocepacia and application to nationwide epidemiological analysis.

Organisms of the Burkholderia cepacia complex are especially important pathogens in cystic fibrosis (CF), with a propensity for patient-to-patient spread and long-term respiratory colonization. B. cenocepacia and Burkholderia multivorans account for the majority of infections in CF, and major epidemic clones have been recognized throughout the world. The aim of the present study was to develop and evaluate a multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) scheme for B. cenocepacia. Potential VNTR loci were identified upon analysis of the annotated genome sequences of B. cenocepacia strains AU1054, J2315, and MCO-3, and 10 of them were selected on the basis of polymorphisms and size. A collection of 100 B. cenocepacia strains, including epidemiologically related and unrelated strains, as well as representatives of the major epidemic lineages, was used to evaluate typeability, epidemiological concordance, and the discriminatory power of MLVA-10 compared with those of pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). Longitudinal stability was assessed by testing 39 successive isolates from 14 patients. Typeability ranged from 0.91 to 1, except for that of one marker, which was not amplified in 53% of the B. cenocepacia IIIA strains. The MLVA types were shown to be stable in chronically colonized patients and within outbreak-related strains, with excellent epidemiological concordance. Epidemic and/or globally distributed lineages (epidemic Edinburgh-Toronto electrophoretic type 12 [ET-12], sequence type 32 [ST-32], ST-122, ST-234, and ST-241) were successfully identified. Conversely, the discriminatory power of MLVA was lower than that of PFGE or MLST, although PFGE variations within the epidemic lineages sometimes masked their genetic relatedness. In conclusion, MLVA represents a promising cost-effective first-line tool in B. cenocepacia surveillance.

High-throughput sequencing of Bacillus anthracis in France: investigating genome diversity and population structure using whole-genome SNP discovery.

BACKGROUND: Single nucleotide polymorphisms (SNPs) are ideal signatures for subtyping monomorphic pathogens such as Bacillus anthracis. Here we report the use of next-generation sequencing technology to investigate the historical, geographic and genetic diversity of Bacillus anthracis in France. 122 strains isolated over a 60-years period throughout the country were whole-genome sequenced and comparative analyses were carried out with a focus on SNPs discovery to discriminate regional sub-groups of strains. RESULTS: A total of 1581 chromosomal SNPs precisely establish the phylogenetic relationships existing between the French strains. Phylogeography patterns within the three canSNP sub-lineages present in France (i.e. B.Br.CNEVA, A.Br.011/009 and A.Br.001/002) were observed. One of the more remarkable findings was the identification of a variety of genotypes within the A.Br.011/009 sub-group that are persisting in the different regions of France. The 560 SNPs defining the A.Br.011/009- affiliated French strains split the Trans-Eurasian sub-group into six distinct branches without any intermediate nodes. Distinct sub-branches, with some geographic clustering, were resolved. The 345 SNPs defining the major B.Br CNEVA sub-lineage clustered three main phylogeographic clades, the Alps, the Pyrenees, and the Massif Central, with a small Saône-et-Loire sub-cluster nested within the latter group. The French strains affiliated to the minor A.Br.001/002 group were characterized by 226 SNPs. All recent isolates collected from the Doubs department were closely related. Identification of SNPs from whole-genome sequences facilitates high-resolution strain tracking and provides the level of discrimination required for outbreak investigations. Eight diagnostic SNPs, representative of the main French-specific phylogeographic clusters, were therefore selected and developed into high-resolution melting SNP discriminative assays. CONCLUSIONS: This work has established one of the most accurate phylogenetic reconstruction of B. anthracis population structure in a country. An extensive next-generation sequencing (NGS) dataset of 122 French strains have been created that allowed the identification of novel diagnostic SNPs useful to rapidly determine the geographic origin of any strain found in France.

Progenitor "Mycobacterium canettii" clone responsible for lymph node tuberculosis epidemic, Djibouti.

"Mycobacterium canettii," an opportunistic human pathogen living in an unknown environmental reservoir, is the progenitor species from which Mycobacterium tuberculosis emerged. Since its discovery in 1969, most of the ≈70 known M. canettii strains were isolated in the Republic of Djibouti, frequently from expatriate children and adults. We show here, by whole-genome sequencing, that most strains collected from February 2010 through March 2013, and associated with 2 outbreaks of lymph node tuberculosis in children, belong to a unique epidemic clone within M. canettii. Evolution of this clone, which has been recovered regularly since 1983, may mimic the birth of M. tuberculosis. Thus, recognizing this organism and identifying its reservoir are clinically important.

Brucella papionis sp. nov., isolated from baboons (Papio spp.).

Two Gram-negative, non-motile, non-spore-forming coccoid bacteria (strains F8/08-60(T) and F8/08-61) isolated from clinical specimens obtained from baboons (Papio spp.) that had delivered stillborn offspring were subjected to a polyphasic taxonomic study. On the basis of 16S rRNA gene sequence similarities, both strains, which possessed identical sequences, were assigned to the genus Brucella. This placement was confirmed by extended multilocus sequence analysis (MLSA), where both strains possessed identical sequences, and whole-genome sequencing of a representative isolate. All of the above analyses suggested that the two strains represent a novel lineage within the genus Brucella. The strains also possessed a unique profile when subjected to the phenotyping approach classically used to separate species of the genus Brucella, reacting only with Brucella A monospecific antiserum, being sensitive to the dyes thionin and fuchsin, being lysed by bacteriophage Wb, Bk2 and Fi phage at routine test dilution (RTD) but only partially sensitive to bacteriophage Tb, and with no requirement for CO2 and no production of H2S but strong urease activity. Biochemical profiling revealed a pattern of enzyme activity and metabolic capabilities distinct from existing species of the genus Brucella. Molecular analysis of the omp2 locus genes showed that both strains had a novel combination of two highly similar omp2b gene copies. The two strains shared a unique fingerprint profile of the multiple-copy Brucella-specific element IS711. Like MLSA, a multilocus variable number of tandem repeat analysis (MLVA) showed that the isolates clustered together very closely, but represent a distinct group within the genus Brucella. Isolates F8/08-60(T) and F8/08-61 could be distinguished clearly from all known species of the genus Brucella and their biovars by both phenotypic and molecular properties. Therefore, by applying the species concept for the genus Brucella suggested by the ICSP Subcommittee on the Taxonomy of Brucella, they represent a novel species within the genus Brucella, for which the name Brucella papionis sp. nov. is proposed, with the type strain F8/08-60(T) ( = NCTC 13660(T) = CIRMBP 0958(T)).

Staphylococcus aureus from 152 cases of bovine, ovine and caprine mastitis investigated by Multiple-locus variable number of tandem repeat analysis (MLVA).

Staphylococcus aureus is one of the main etiological agents of mastitis in ruminants. In the present retrospective study, we evaluated the potential interest of a previously described automated multiple loci Variable Number of Tandem Repeats (VNTR) Assay (MLVA) comprising 16 loci as a first line tool to investigate the population structure of S. aureus from mastitis. We determined the genetic diversity of S. aureus strains from cases of clinical and subclinical mastitis in dairy cattle (n = 118, of which 16 were methicillin-resistant), sheep (n = 18) and goats (n = 16). The 152 strains could be subdivided into 115 MLVA genotypes (including 14 genotypes for the ovine strains and 15 genotypes for the caprine strains). This corresponds to a discriminatory index (D) value of 0.9936. Comparison with published MLVA data obtained using the same protocol applied to strains from diverse human and animal origins revealed a low number (8.5%) of human-related MLVA genotypes among the present collection. Eighteen percent of the S. aureus mastitis collection belonged to clonal complexes apparently not associated with other pathological conditions. Some of them displayed a relatively low level of diversity in agreement with a restricted ecological niche. These findings provide arguments suggesting that specific S. aureus lineages particularly adapted to ruminant mammary glands have emerged and that MLVA is a convenient tool to provide a broad overview of the population, owing to the availability via internet of databases compiling published MLVA genotypes.

Avirulence of a spontaneous Francisella tularensis subsp. mediasiatica prmA mutant.

Francisella tularensis, the causative agent of tularemia, is divided into three subspecies. Two of these, subspecies holarctica and tularensis, are highly pathogenic to humans and consequently relatively well studied. The third subspecies, mediasiatica, is rarely isolated and remains poorly studied. It is distributed in the sparsely populated regions of Central Asia and Siberia. Curently this subspecies is not known to have been responsible for human infections in spite of its high virulence in laboratory animals. Subspecies mediasiatica is currently divided into three subgroups-MI, present in Central Asia, MII, present in southern Siberia, and MIII represented by a unique strain, 60(B)57, isolated in Uzbekistan in 1960. We describe here the unexpected observation that MIII strain 60(B)57 is avirulent and immunogenic. We observed that infection with this strain protected mice from challenge 21 days later with a virulent subsp. mediasiatica strain. With an increase of this interval, the protection for mice was significantly reduced. In contrast, guinea pigs were protected from challenge with strains of the subspecies holarctica and mediasiatica (but not subsp. tularensis) 90 days after infection with 60(B)57. We performed genome assembly based on whole genome sequencing data obtained using the Nanopore MinION for strain 60(B)57 and two subsp. mediasiatica strains representing the Central Asian MI and Siberian MII phylogenetic subgroups. The prmA gene is truncated due to a nonsense mutation in strain 60(B)57. The deletion of gene prmA has previously been shown to induce a loss of virulence in Francisella novicida the closest model organism suggesting that the observed mutation might the cause of the avirulence of strain 60(B)57.

Brucella abortus in Kazakhstan, population structure and comparison with worldwide genetic diversity.

Brucella abortus is the main causative agent of brucellosis in cattle, leading to severe economic consequences in agriculture and affecting public health. The zoonotic nature of the infection increases the need to control the spread and dynamics of outbreaks in animals with the incorporation of high resolution genotyping techniques. Based on such methods, B. abortus is currently divided into three clades, A, B, and C. The latter includes subclades C1 and C2. This study presents the results of whole-genome sequencing of 49 B. abortus strains isolated in Kazakhstan between 1947 and 2015 and of 36 B. abortus strains of various geographic origins isolated from 1940 to 2004. In silico Multiple Locus Sequence Typing (MLST) allowed to assign strains from Kazakhstan to subclades C1 and to a much lower extend C2. Whole-genome Single-Nucleotide Polymorphism (wgSNP) analysis of the 46 strains of subclade C1 with strains of worldwide origins showed clustering with strains from neighboring countries, mostly North Caucasia, Western Russia, but also Siberia, China, and Mongolia. One of the three Kazakhstan strains assigned to subclade C2 matched the B. abortus S19 vaccine strain used in cattle, the other two were genetically close to the 104 M vaccine strain. Bayesian phylodynamic analysis dated the introduction of B. abortus subclade C1 into Kazakhstan to the 19th and early 20th centuries. We discuss this observation in view of the history of population migrations from Russia to the Kazakhstan steppes.

New Research on the Bacillus anthracis Genetic Diversity in Siberia.

Anthrax is a particularly dangerous infection of humans and ungulates caused by the Gram-positive spore-forming bacterium Bacillus anthracis. The highly monomorphic and clonal species B. anthracis is commonly divided into three main lineages, A, B, and C, which in turn are divided into several canSNP groups. We report here a phylogenetic analysis based on the whole-genome sequence (WGS) data of fifteen strains isolated predominantly in Siberia or Central and Southern Russia. We confirm the wide distribution of the cluster of strains of the B.Br.001/002 group, endemic to the Russian Arctic, which is also present in the steppe zone of Southern Siberia. We characterize additional branches within the major A.Br.001/002 polytomy comprising the A.Br.Ames and A.Br.Sterne lineages, one of which is identified in the Arctic.

Rapid identification of international multidrug-resistant Pseudomonas aeruginosa clones by multiple-locus variable number of tandem repeats analysis and investigation of their susceptibility to lytic bacteriophages.

The objective of this study was to determine the genetic diversity of multidrug-resistant (MDR) Pseudomonas aeruginosa strains isolated over a period of 12 months in two French hospitals and to test their susceptibility to bacteriophages. A total of 47 MDR isolates recovered from hospitalized patients were genotyped using multiple-locus variable number of tandem repeats analysis. The genotypes were distributed into five clones (including 19, 5, 5, 3, and 3 isolates, respectively) and 12 singletons. Comparison to 77 MDR strains from three other countries, and MLST analysis of selected isolates showed the predominance of international MDR clones. The larger clone, CC235, contained 59 isolates displaying different antibiotic resistance mechanisms, including the presence of the GES1, VIM-2, VIM-4, and IMP-1 ß-lactamases. Three newly isolated P. aeruginosa bacteriophages were found to lyse 42 of the 44 analyzed strains, distributed into the different clonal complexes. This pilot study suggests that systematic genotyping of P. aeruginosa MDR strains could improve our epidemiological understanding of transmission at both the local (hospital) and the national level and that phage therapy could be an alternative or a complementary treatment to antibiotics for treating MDR-infected patients.

Intraspecies biodiversity of the genetically homologous species Brucella microti.

Brucellosis is one of the major bacterial zoonoses worldwide. In the past decade, an increasing number of atypical Brucella strains and species have been described. Brucella microti in particular has attracted attention, because this species not only infects mammalian hosts but also persists in soil. An environmental reservoir may pose a new public health risk, leading to the reemergence of brucellosis. In a polyphasic approach, comprising conventional microbiological techniques and extensive biochemical and molecular techniques, all currently available Brucella microti strains were characterized. While differing in their natural habitats and host preferences, B. microti isolates were found to possess identical 16S rRNA, recA, omp2a, and omp2b gene sequences and identical multilocus sequence analysis (MLSA) profiles at 21 different genomic loci. Only highly variable microsatellite markers of multiple-locus variable-number tandem repeat (VNTR) analysis comprising 16 loci (MLVA-16) showed intraspecies discriminatory power. In contrast, biotyping demonstrated striking differences within the genetically homologous species. The majority of the mammalian isolates agglutinated only with monospecific anti-M serum, whereas soil isolates agglutinated with anti-A, anti-M, and anti-R sera. Bacteria isolated from animal sources were lysed by phages F1, F25, Tb, BK2, Iz, and Wb, whereas soil isolates usually were not. Rough strains of environmental origin were lysed only by phage R/C. B. microti exhibited high metabolic activities similar to those of closely related soil organisms, such as Ochrobactrum spp. Each strain was tested with 93 different substrates and showed an individual metabolic profile. In summary, the adaptation of Brucella microti to a specific habitat or host seems to be a matter of gene regulation rather than a matter of gene configuration.

Whole genome sequence analysis of Neisseria meningitidis strains circulating in Kazakhstan, 2017-2018.

Neisseria meningitidis (meningococcus) is a cosmopolitan bacterium that is often found in the upper respiratory tract of asymptomatic humans. However, N. meningitidis also causes meningeal inflammation and/or sepsis in humans with a periodic resurgence in incidence and high mortality rates. The pathogen is highly diverse genetically and antigenically, so that genotyping is considered important for vaccine matching to circulating strains. Annual incidence of meningococcal disease in Kazakhstan ranges between 0.2 and 2.5 cases per 100 thousand population. In total, 78 strains of N. meningitidis were isolated from clinical patients and contact persons during the years 2017-2018 in Kazakhstan. Of these, 41 strains including four from the patients and 37 from contacts, were sequenced using Illumina MiSeq. In silico typing was completed using the Neisseria pipeline 1.2 on the Galaxy Workflow Management System and PubMLST. Whole genome SNP (single nucleotide polymorphisms) trees were built using BioNumerics 8. Seven-gene multilocus sequence typing (MLST) identified ten sequence types (ST), two of which have not been previously described (ST-16025; ST-16027). ST-16025 was detected in two patients with invasive meningococcal disease in 2017 and 2018 in Akmola region and 16 contacts in 2017 in Turkistan region. This prevalent type ST-16025 demonstrates considerable intertypic diversity as it consists of three subcomplexes with a distance of more than 2000 SNPs. Invasive and carrier strains belong to different serogroups (MenB and MenC), PorA and FetA_VR. Two invasive strains were MenB, one MenC and one MenW (Hajj lineage). The strains from the contact persons were: MenC (n = 18), cnl (n = 9), MenY (n = 7), MenW (n = 1), MenB (n = 1) and one unidentifiable. Different numbers of alleles were present: 12, 11, 7, and 7 alleles for PorA, FetA, fHbp, and NHBA, respectively. This study is the first report of the genetic diversity of N. meningitidis strains in Kazakhstan. Despite limitations with the studied sample size, important conclusions can be drawn based on data produced. This study provides evidence for regulatory authorities with regard to changing routine diagnostic protocols to increase the collecting of samples for WGS.

Investigation of the population structure of Legionella pneumophila by analysis of tandem repeat copy number and internal sequence variation.

The population structure of the species Legionella pneumophila was investigated by multilocus variable number of tandem repeats (VNTR) analysis (MLVA) and sequencing of three VNTRs (Lpms01, Lpms04 and Lpms13) in selected strains. Of 150 isolates of diverse origins, 136 (86 %) were distributed into eight large MLVA clonal complexes (VACCs) and the rest were either unique or formed small clusters of up to two MLVA genotypes. In spite of the lower degree of genome-wide linkage disequilibrium of the MLVA loci compared with sequence-based typing, the clustering achieved by the two methods was highly congruent. The detailed analysis of VNTR Lpms04 alleles showed a very complex organization, with five different repeat unit lengths and a high level of internal variation. Within each MLVA-defined VACC, Lpms04 was endowed with a common recognizable pattern with some interesting exceptions. Evidence of recombination events was suggested by analysis of internal repeat variations at the two additional VNTR loci, Lpms01 and Lpms13. Sequence analysis of L. pneumophila VNTR locus Lpms04 alone provides a first-line assay for allocation of a new isolate within the L. pneumophila population structure and for epidemiological studies.