Molecular Epidemiology of Mycoplasma pneumoniae: Genotyping Using Single Nucleotide Polymorphisms and SNaPshot Technology.

Molecular typing of Mycoplasma pneumoniae is an important tool for identifying grouped cases and investigating outbreaks. In the present study, we developed a new genotyping method based on single nucleotide polymorphisms (SNPs) selected from the whole-genome sequencing of eight M. pneumoniae strains, using the SNaPshot minisequencing assay. Eight SNPs, localized in housekeeping genes, predicted lipoproteins, and adhesin P1 genes were selected for genotyping. These SNPs were evaluated on 140 M. pneumoniae clinical isolates previously genotyped by multilocus variable-number tandem-repeat analysis (MLVA-5) and adhesin P1 typing. This method was also adapted for direct use with clinical samples and evaluated on 51 clinical specimens. The analysis of the clinical isolates using the SNP typing method showed nine distinct SNP types with a Hunter and Gaston diversity index (HGDI) of 0.836, which is higher than the HGDI of 0.583 retrieved for the MLVA-4 typing method, where the nonstable Mpn1 marker was removed. A strong correlation with the P1 adhesin gene typing results was observed. The congruence was poor between MLVA-5 and SNP typing, indicating distinct genotyping schemes. Combining the results increased the discriminatory power. This new typing method based on SNPs and the SNaPshot technology is a method for rapid M. pneumoniae typing directly from clinical specimens, which does not require any sequencing step. This method is based on stable markers and provides information distinct from but complementary to MLVA typing. The combined use of SNPs and MLVA typing provides powerful discrimination of strains.

Molecular characterisation of Brucella species.

The genus Brucella (Mayer and Shaw, 1920) currently consists often species with validly published names. Within most species further differentiation into biovars exists. Genetically, all Brucella species are highly related to each other, exhibiting sequence similarity values of 98% to 100% in aligned regions (core genome). The population structure is clonal. Despite this close genetic relatedness, the various species can be clearly distinguished from each other by application of high-resolution molecular typing tools, in addition to assessment of phenotype and host preference. Accurate species delineation can be achieved by conventional multiplex polymerase chain reaction (PCR), single nucleotide polymorphism (SNP) analysis and multilocus sequence typing (MLST) or multilocus sequence analysis (MLSA). The last is also suitable for phylogenetic reconstructions, owing to the highly clonal evolution of the different species. Highly discriminatory multilocus variable number of tandem repeats (VNTR) analysis (MLVA) allows both species delineation and differentiation of individual isolates and thus represents a perfect first-line toolfor molecular epidemiological studies within outbreak investigations. More recently,whole genome sequencing (WGS)and the resulting global genome-wide SNP analysis have become available. These novel approaches should help in further understanding the evolution, host specificity and pathogenicity of the genus Brucella.

Use of multilocus variable-number tandem repeat analysis (MLVA) in eight European countries, 2012.

Genotyping of important medical or veterinary prokaryotes has become a very important tool during the last decades. Rapid development of fragment-separation and sequencing technologies has made many new genotyping strategies possible. Among these new methods is multilocus variable-number tandem repeat analysis (MLVA). Here we present an update on the use of MLVA in eight European countries (Denmark, France, Germany, Ireland, Italy, the Netherlands, Norway and Sweden). Researchers in Europe have been active in developing and implementing a large array of different assays. MLVA has been used as a typing tool in several contexts, from aiding in resolving outbreaks of foodborne bacteria to typing organisms that may pose a bioterrorist threat, as well as in scientific studies.

A new highly discriminatory multiplex capillary-based MLVA assay as a tool for the epidemiological survey of Pseudomonas aeruginosa in cystic fibrosis patients.

Multiple locus variable number of tandem repeats (VNTR) analysis (MLVA) has been shown to provide a high level of information for epidemiological investigations and the follow-up of Pseudomonas aeruginosa chronic infection. In the present study, an automatized MLVA assay has been developed for the analysis of 16 VNTRs in two multiplex polymerase chain reactions (PCRs), followed by capillary electrophoresis. The result in the form of a code is directly usable for clustering analyses. This MLVA-16(Orsay) scheme was applied to the genotyping of 83 isolates from eight cystic fibrosis patients, demonstrating that the same genotype persisted during eight years of chronic infection in the majority of cases. Comparison with pulsed-field gel electrophoresis (PFGE) analysis showed that both methods were congruent, MLVA providing, in some cases, additional informativity. The evolution of strains during long-term infection was revealed by the presence of VNTR variants.

Characterization of genetic diversity of Bacillus anthracis in France by using high-resolution melting assays and multilocus variable-number tandem-repeat analysis.

Using high-resolution melting (HRM) analysis, we developed a cost-effective method to genotype a set of 13 phylogenetically informative single-nucleotide polymorphisms (SNPs) within the genome of Bacillus anthracis. SNP discrimination assays were performed in monoplex or duplex and applied to 100 B. anthracis isolates collected in France from 1953 to 2009 and a few reference strains. HRM provided a reliable and cheap alternative to subtype B. anthracis into one of the 12 major sublineages or subgroups. All strains could be correctly positioned on the canonical SNP (canSNP) phylogenetic tree, except the divergent Pasteur vaccine strain ATCC 4229. We detected the cooccurrence of three canSNP subgroups in France. The dominant B.Br.CNEVA sublineage was found to be prevalent in the Alps, the Pyrenees, the Auvergne region, and the Saône-et-Loire department. Strains affiliated with the A.Br.008/009 subgroup were observed throughout most of the country. The minor A.Br.001/002 subgroup was restricted to northeastern France. Multiple-locus variable-number tandem-repeat analysis using 24 markers further resolved French strains into 60 unique profiles and identified some regional patterns. Diversity found within the A.Br.008/009 and B.Br.CNEVA subgroups suggests that these represent old, ecologically established clades in France. Phylogenetic relationships with strains from other parts of the world are discussed.

High-throughput typing method to identify a non-outbreak-involved Legionella pneumophila strain colonizing the entire water supply system in the town of Rennes, France.

Two legionellosis outbreaks occurred in the city of Rennes, France, during the past decade, requiring in-depth monitoring of Legionella pneumophila in the water network and the cooling towers in the city. In order to characterize the resulting large collection of isolates, an automated low-cost typing method was developed. The multiplex capillary-based variable-number tandem repeat (VNTR) (multiple-locus VNTR analysis [MLVA]) assay requiring only one PCR amplification per isolate ensures a high level of discrimination and reduces hands-on and time requirements. In less than 2 days and using one 4-capillary apparatus, 217 environmental isolates collected between 2000 and 2009 and 5 clinical isolates obtained during outbreaks in 2000 and 2006 in Rennes were analyzed, and 15 different genotypes were identified. A large cluster of isolates with closely related genotypes and representing 77% of the population was composed exclusively of environmental isolates extracted from hot water supply systems. It was not responsible for the known Rennes epidemic cases, although strains showing a similar MLVA profile have regularly been involved in European outbreaks. The clinical isolates in Rennes had the same genotype as isolates contaminating a mall's cooling tower. This study further demonstrates that unknown environmental or genetic factors contribute to the pathogenicity of some strains. This work illustrates the potential of the high-throughput MLVA typing method to investigate the origin of legionellosis cases by allowing the systematic typing of any new isolate and inclusion of data in shared databases.

Molecular epidemiology of human and animal tuberculosis in Ibadan, Southwestern Nigeria.

From 2005 to 2007, Mycobacterium tuberculosis complex (MTC) strains were isolated from cattle, goats and pigs samples collected at the Bodija abattoir and from human samples from tuberculosis patients and livestock traders at the Akinyele cattle market in Ibadan, Southwestern Nigeria. Seventy four isolates obtained from humans (24) and livestock (50) were identified as MTC strains. Thirty two isolates were spoligotyped. Nineteen of these 32 isolates were identified as M. tuberculosis whilst 13 were identified as Mycobacterium bovis. M. bovis was isolated from two humans, whereas M. tuberculosis was isolated from a bovine, a pig and a goat. All the M. bovis isolates identified in this study belonged to the Africa 1 clonal complex. Multiple locus VNTR [variable number of tandem repeats] analysis (MLVA) was carried out on the 74 isolates. Three major clusters were defined. Group A consisted of 24 M. tuberculosis isolates (MLVA genotypes 1-18). One strain was isolated from a bovine and one from a pig. Group B consisted of 49 M. bovis strains (MLVA genotypes 19-48), mainly of cattle origin but also included four goat, nine pig and two human isolates. Group C consisted of a single M. tuberculosis isolate (MLVA genotype 49) obtained from a goat. Spoligotyping and MLVA confirmed it as clustering with the East Africa Indian clade found in humans in Sudan and the Republic of Djibouti. The isolation of three M. tuberculosis strains from livestock raises the question of their epidemiological importance as a source of infection for humans.

[Performances of the assay MTBDRplus(®) in the surveillance of rifampicin resistance in Mycobacterium tuberculosis]. / Performances du kit MTBDRplus(®) dans le cadre de la surveillance de la résistance à la rifampicine chez Mycobacterium tuberculosis.

The purpose of the survey was the routine assessment of the MTBDRplus(®) kit performance in the determination and characterization of Mycobacterium tuberculosis resistance to rifampicin. The survey was carried out on a collection of 144 strains (126 of which were resistant to rifampicin) isolated on patients from 15 countries. Sensitivity to antituberculosis drugs was determined by a liquid culture system and the reference method was the amplification and sequencing of a target region of the rpoB gene whose mutations are responsible for rifampicin resistance (codons 507 to 533). The assessed kit was based on a reverse hybridization technique using eight overlapping probes covering the target region and four probes representing the most-frequently observed mutations. The assay performance was found excellent, specificity: 100%, sensitivity: 99.2%; 17 mutations affecting 10 codons were reported, two of which were newly identified.

CRISPR elements in Yersinia pestis acquire new repeats by preferential uptake of bacteriophage DNA, and provide additional tools for evolutionary studies.

The remarkable repetitive elements called CRISPRs (clustered regularly interspaced short palindromic repeats) consist of repeats interspaced with non-repetitive elements or 'spacers'. CRISPRs are present in both archaea and bacteria, in association with genes involved in DNA recombination and repair. In the Yersinia pestis genome, three such elements are found at three distinct loci, one of them being highly polymorphic. The authors have sequenced a total of 109 alleles of the three Y. pestis CRISPRs and they describe 29 new spacers, most being specific to one isolate. In nine strains of Yersinia pseudotuberculosis, 132 spacers were found, of which only three are common to Y. pestis isolates. In Y. pestis of the Orientalis biovar investigated in detail here, deletion of motifs is observed but it appears that addition of new motifs to a common ancestral element is the most frequent event. This takes place at the three different loci, although at a higher rate in one of the loci, and the addition of new motifs is polarized. Interestingly, the most recently acquired spacers were found to have a homologue at another locus in the genome, the majority of these inside an inactive prophage. This is believed to be the first time that the origin of the spacers in CRISPR elements has been explained. The CRISPR structure provides a new and robust identification tool.

Tandem repeats analysis for the high resolution phylogenetic analysis of Yersinia pestis.

BACKGROUND: Yersinia pestis, the agent of plague, is a young and highly monomorphic species. Three biovars, each one thought to be associated with the last three Y. pestis pandemics, have been defined based on biochemical assays. More recently, DNA based assays, including DNA sequencing, IS typing, DNA arrays, have significantly improved current knowledge on the origin and phylogenetic evolution of Y. pestis. However, these methods suffer either from a lack of resolution or from the difficulty to compare data. Variable number of tandem repeats (VNTRs) provides valuable polymorphic markers for genotyping and performing phylogenetic analyses in a growing number of pathogens and have given promising results for Y. pestis as well. RESULTS: In this study we have genotyped 180 Y. pestis isolates by multiple locus VNTR analysis (MLVA) using 25 markers. Sixty-one different genotypes were observed. The three biovars were distributed into three main branches, with some exceptions. In particular, the Medievalis phenotype is clearly heterogeneous, resulting from different mutation events in the napA gene. Antiqua strains from Asia appear to hold a central position compared to Antiqua strains from Africa. A subset of 7 markers is proposed for the quick comparison of a new strain with the collection typed here. This can be easily achieved using a Web-based facility, specifically set-up for running such identifications. CONCLUSION: Tandem-repeat typing may prove to be a powerful complement to the existing phylogenetic tools for Y. pestis. Typing can be achieved quickly at a low cost in terms of consumables, technical expertise and equipment. The resulting data can be easily compared between different laboratories. The number and selection of markers will eventually depend upon the type and aim of investigations.

Characterization of a tandem repeat polymorphism in Legionella pneumophila and its use for genotyping.

We have analyzed the variability of minisatellite sequences (also called variable-number tandem repeats [VNTRs]) in the genome of Legionella pneumophila. Based upon the genome sequence of the Philadelphia-1 strain (serogroup 1), 25 minisatellites were selected and their polymorphisms were analyzed by PCR with the DNA of serogroup 1 to 14 reference strains. For 22 markers, a PCR product of the expected size was found with the DNA of the Philadelphia-1 strain. Most of these markers did not amplify the DNA of other Legionella species or other bacteria used as controls. A polymorphism was observed for seven markers among the L. pneumophila strains tested. To check whether these markers could be used to compare strains of L. pneumophila, we analyzed two groups of isolates from clinical and environmental samples which had been independently genotyped by other methods. The results showed that, for the isolates in these two sets of samples, VNTR typing is as informative as pulsed-field gel electrophoresis for comparison of strains. Sequencing of one minisatellite from 14 reference strains was performed. Comparison of the sequences allowed a classification and confirmed the existence of subspecies of L. pneumophila. We also tested the usefulness of one very polymorphic marker as a tool for the rapid screening of colonies grown from water samples. This allowed the rapid identification of the L. pneumophila colonies and gave a first hint as to the presence of several strains in a single sample.

A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis.

BACKGROUND: Some pathogenic bacteria are genetically very homogeneous, making strain discrimination difficult. In the last few years, tandem repeats have been increasingly recognized as markers of choice for genotyping a number of pathogens. The rapid evolution of these structures appears to contribute to the phenotypic flexibility of pathogens. The availability of whole-genome sequences has opened the way to the systematic evaluation of tandem repeats diversity and application to epidemiological studies. RESULTS: This report presents a database (http://minisatellites.u-psud.fr) of tandem repeats from publicly available bacterial genomes which facilitates the identification and selection of tandem repeats. We illustrate the use of this database by the characterization of minisatellites from two important human pathogens, Yersinia pestis and Bacillus anthracis. In order to avoid simple sequence contingency loci which may be of limited value as epidemiological markers, and to provide genotyping tools amenable to ordinary agarose gel electrophoresis, only tandem repeats with repeat units at least 9 bp long were evaluated. Yersinia pestis contains 64 such minisatellites in which the unit is repeated at least 7 times. An additional collection of 12 loci with at least 6 units, and a high internal conservation were also evaluated. Forty-nine are polymorphic among five Yersinia strains (twenty-five among three Y. pestis strains). Bacillus anthracis contains 30 comparable structures in which the unit is repeated at least 10 times. Half of these tandem repeats show polymorphism among the strains tested. CONCLUSIONS: Analysis of the currently available bacterial genome sequences classifies Bacillus anthracis and Yersinia pestis as having an average (approximately 30 per Mb) density of tandem repeat arrays longer than 100 bp when compared to the other bacterial genomes analysed to date. In both cases, testing a fraction of these sequences for polymorphism was sufficient to quickly develop a set of more than fifteen informative markers, some of which show a very high degree of polymorphism. In one instance, the polymorphism information content index reaches 0.82 with allele length covering a wide size range (600-1950 bp), and nine alleles resolved in the small number of independent Bacillus anthracis strains typed here.

A dominant modifier of transgene methylation is mapped by QTL analysis to mouse chromosome 13.

The single-copy hepatitis B virus transgene in the E36 transgenic mouse strain undergoes methylation changes in a parent-of-origin, tissue, and strain-specific fashion. In a C57BL/6 background, the paternally transmitted transgene is methylated in 30% of cells, whereas it is methylated in more than 80% of cells in (BALB/c x C57BL/6) F1 mice. We established previously that several genetic factors were likely to contribute to the transgene methylation profile, some with demethylating and some with de novo methylating activities. Using quantitative trait loci (QTL) mapping, we have now localized one major modifier locus on chromosome 13 (Mod13), which explains a 30% increase in the methylation level of this transgene with no effect on the flanking endogenous sequences. No other QTL could be identified, except for a demethylating activity of low significance located on chromosome 12. Recombinant inbred mice containing a BALB/c allele of Mod13 were then used to show that the presence of Mod13 is sufficient to induce de novo methylation. A segregation between de novo methylation and repression of transgene expression was uncovered, suggesting that this genetic system is also useful for the identification of factors that interpret methylation patterns in the genome.

Children of chernobyl cleanup workers do not show elevated rates of mutations in minisatellite alleles.

The disaster at the Chernobyl Nuclear Power Plant in April 1986 was accompanied by the release of large amounts of radioisotopes, resulting in the contamination of extensive regions of the Ukraine, Byelorus and the Russian Federation. Cleanup workers (liquidators) and people living on land contaminated with radioactive materials were most exposed. To assess the genetic effects of exposure to ionizing radiation after the Chernobyl accident, we have measured the frequency of inherited mutant alleles at seven hypermutable minisatellite loci in 183 children born to Chernobyl cleanup workers (liquidators) and 163 children born to control families living in nonirradiated areas of the Ukraine. There was no significant difference in the frequency of inherited mutant alleles between the exposed and control groups. The exposed group was then divided into two subgroups according to the time at which the children were conceived with respect to the fathers' work at the power plant. Eighty-eight children were conceived either while their fathers were working at the facility or up to 2 months later (Subgroup 1). The other 95 children were conceived at least 4 months after their fathers had stopped working at the Chernobyl site (Subgroup 2). The frequencies of mutant alleles were higher for the majority of loci (i.e. 1.44 times higher for CEB1) in Subgroup 1 than in Subgroup 2. This result, if confirmed, would reconcile the apparently conflicting results obtained in the chronically exposed Byelorus population and the Hiroshima-Nagasaki A-bomb survivors.

Minisatellites: mutability and genome architecture.

Minisatellites have been found in association with important features of human genome biology such as gene regulation, chromosomal fragile sites, and imprinting. Our knowledge of minisatellite biology has greatly increased in the past 10 years owing to the identification and careful analysis of human hypermutable minisatellites, experimental models in yeast, and recent in vitro studies of minisatellite recombination properties. In parallel, minisatellites have been put forward as potential biomarkers for the monitoring of genotoxic agents such as ionizing radiation. We summarize and discuss recent observations on minisatellites. In addition we take advantage of recent whole chromosome sequence data releases to provide a unifying view which may facilitate the annotation of tandem repeat sequences.

The SH2D2A gene encoding the T-cell-specific adapter protein (TSAd) is localized centromeric to the CD1 gene cluster on human Chromosome 1.

The SH2D2A gene encoding the T-cell-specific adapter protein (TSAd), was isolated from a human Chromosome (Chr) 1 cosmid library (LLNL, UK HGMP). The gene spans 11 kilobases and contains nine exons and eight introns. Four alternative transcript variants were observed in activated T cells. Three single-nucleotide polymorphisms were identified within intron 2. A variable number of GA repeats was found at position -340 from the first coding ATG. Linkage analysis using this marker in eight CEPH families showed that the SH2D2A gene is located close to the D1S2624 marker on Chr 1q21-1q22. Physical mapping of a PAC and BAC contig containing the CD1 gene cluster telomeric to D1S2624 failed to identify a clone harboring the SH2D2A gene. Thus the SH2D2A gene is located centromeric to the CD1 gene cluster on Chr 1.

Meiotic instability of human minisatellite CEB1 in yeast requires DNA double-strand breaks.

Minisatellites are tandemly repeated DNA sequences of 10-100-bp units. Some minisatellite loci are highly unstable in the human germ line, and structural analysis of mutant alleles has suggested that repeat instability results from a recombination-based process. To provide insights into the molecular mechanism of human minisatellite instability, we developed Saccharomyces cerevisiae strains carrying alleles of the most unstable human minisatellite locus, CEB1 (ref. 2). We observed that CEB1 is destabilized in meiosis, resulting in a variety of intra- and inter-allelic gains or losses of repeat units, similar to rearrangements described in humans. Using mutations affecting the initiation of recombination (spo11) or mismatch repair (msh2 pms1 ), we demonstrate that meiotic destabilization depends on the initiation of homologous recombination at nearby DNA double-strand break (DSBs) sites and involves a 'rearranged heteroduplex' intermediate. Most of the human and yeast data can be explained and unified in the context of DSB repair models.

Finding new human minisatellite sequences in the vicinity of long CA-rich sequences.

Microsatellites and minisatellites are two classes of tandem repeat sequences differing in their size, mutation processes, and chromosomal distribution. The boundary between the two classes is not defined. We have developed a convenient, hybridization-based human library screening procedure able to detect long CA-rich sequences. Analysis of cosmid clones derived from a chromosome 1 library show that cross-hybridizing sequences tested are imperfect CA-rich sequences, some of them showing a minisatellite organization. All but one of the 13 positive chromosome 1 clones studied are localized in chromosomal bands to which minisatellites have previously been assigned, such as the 1pter cluster. To test the applicability of the procedure to minisatellite detection on a larger scale, we then used a large-insert whole-genome PAC library. Altogether, 22 new minisatellites have been identified in positive PAC and cosmid clones and 20 of them are telomeric. Among the 42 positive PAC clones localized within the human genome by FISH and/or linkage analysis, 25 (60%) are assigned to a terminal band of the karyotype, 4 (9%) are juxtacentromeric, and 13 (31%) are interstitial. The localization of at least two of the interstitial PAC clones corresponds to previously characterized minisatellite-containing regions and/or ancestrally telomeric bands, in agreement with this minisatellite-like distribution. The data obtained are in close agreement with the parallel investigation of human genome sequence data and suggest that long human (CA)s are imperfect CA repeats belonging to the minisatellite class of sequences. This approach provides a new tool to efficiently target genomic clones originating from subtelomeric domains, from which minisatellite sequences can readily be obtained. [The sequence data described in this paper have been submitted to the EMBL data library under accession nos. AJ000377-AJ000383.]