Genotyping of Genetically Monomorphic Bacteria: DNA Sequencing in Mycobacterium tuberculosis Highlights the Limitations of Current Methodologies

Because genetically monomorphic bacterial pathogens harbour little DNA sequence diversity, most current genotyping techniques used to study the epidemiology of these organisms are based on mobile or repetitive genetic elements. Molecular markers commonly used in these bacteria include Clustered Regulatory Short Palindromic Repeats (CRISPR) and Variable Number Tandem Repeats (VNTR). These methods are also increasingly being applied to phylogenetic and population genetic studies. Using the Mycobacterium tuberculosis complex (MTBC) as a model, we evaluated the phylogenetic accuracy of CRISPR- and VNTR-based genotyping, which in MTBC are known as spoligotyping and Mycobacterial Interspersed Repetitive Units (MIRU)-VNTR-typing, respectively. We used as a gold standard the complete DNA sequences of 89 coding genes from a global strain collection. Our results showed that phylogenetic trees derived from these multilocus sequence data were highly congruent and statistically robust, irrespective of the phylogenetic methods used. By contrast, corresponding phylogenies inferred from spoligotyping or 15-loci-MIRU-VNTR were incongruent with respect to the sequence-based trees. Although 24-loci-MIRU-VNTR performed better, it was still unable to detect all strain lineages. The DNA sequence data showed virtually no homoplasy, but the opposite was true for spoligotyping and MIRU-VNTR, which was consistent with high rates of convergent evolution and the low statistical support obtained for phylogenetic groupings defined by these markers. Our results also revealed that the discriminatory power of the standard 24 MIRU-VNTR loci varied by strain lineage. Taken together, our findings suggest strain lineages in MTBC should be defined based on phylogenetically robust markers such as single nucleotide polymorphisms or large sequence polymorphisms, and that for epidemiological purposes, MIRU-VNTR loci should be used in a lineage-dependent manner. Our findings have implications for strain typing in other genetically monomorphic bacteria.     

A New Double Digestion Ligation Mediated Suppression PCR Method for Simultaneous Bacteria DNA-Typing and Confirmation of Species: An Acinetobacter sp. Model

We have designed a new ddLMS PCR (double digestion Ligation Mediated Suppression PCR) method based on restriction site polymorphism upstream from the specific target sequence for the simultaneous identification and differentiation of bacterial strains. The ddLMS PCR combines a simple PCR used for species or genus identification and the LM PCR strategy for strain differentiation. The bacterial identification is confirmed in the form of the PCR product(s), while the length of the PCR product makes it possible to differentiate between bacterial strains. If there is a single copy of the target sequence within genomic DNA, one specific PCR product is created (simplex ddLMS PCR), whereas for multiple copies of the gene the fingerprinting patterns can be obtained (multiplex ddLMS PCR). The described ddLMS PCR method is designed for rapid and specific strain differentiation in medical and microbiological studies. In comparison to other LM PCR it has substantial advantages: enables specific species'' DNA-typing without the need for pure bacterial culture selection, is not sensitive to contamination with other cells or genomic DNA, and gives univocal “band-based” results, which are easy to interpret. The utility of ddLMS PCR was shown for Acinetobacter calcoaceticus-baumannii (Acb) complex, the genetically closely related and phenotypically similar species and also important nosocomial pathogens, for which currently, there are no recommended methods for screening, typing and identification. In this article two models are proposed: 3′ recA-ddLMS PCR-MaeII/RsaI for Acb complex interspecific typing and 5′ rrn-ddLMS PCR-HindIII/ApaI for Acinetobacter baumannii intraspecific typing. ddLMS PCR allows not only for DNA-typing but also for confirmation of species in one reaction. Also, practical guidelines for designing a diagnostic test based on ddLMS PCR for genotyping different species of bacteria are provided.     

Multilocus sequence typing of Salmonella strains by high-throughput sequencing of selectively amplified target genes

Rapid development of next generation sequencing (NGS) technologies in recent years has made whole genome sequencing of bacterial genomes widely accessible. However, it is often unnecessary or not feasible to sequence the whole genome for most applications of genetic analyses in bacteria. Selectively capturing defined genomic regions followed by NGS analysis could be a promising approach for high-resolution molecular typing of a large set of strains. In this study, we describe a novel and straightforward PCR-based target-capturing method, hairpin-primed multiplex amplification (HPMA), which allows for simultaneous amplification of numerous target genes. To test the feasibility of NGS-based strain typing using HPMA, 20 target gene sequences were simultaneously amplified with barcode tagging in each of 41 Salmonella strains. The amplicons were then pooled and analyzed by 454 pyrosequencing. Analysis of the sequence data, as an extension of multilocus sequence typing (MLST), demonstrated the utility and potential of this novel typing method, MLST-seq, as a high-resolution strain typing method. With the rapidly increasing sequencing capacity of NGS, MLST-seq or its variations using different target enrichment methods can be expected to become a high-resolution typing method in the near future for high-throughput analysis of a large collection of bacterial strains.     

Molecular approaches: advantages and artifacts in assessing bacterial diversity

Bacteria account for a major proportion of Earth’s biological diversity. They play essential roles in quite diverse environments and there has been an increasing interest in bacterial biodiversity. Research using novel and efficient tools to identify and characterize bacterial communities has been the key for elucidating biological activities with potential for industrial application. The current approach used for defining bacterial species is based on phenotypic and genomic properties. Traditional and novel DNA-based molecular methods are improving our knowledge of bacterial diversity in nature. Advances in molecular biology have been important for studies of diversity, considerably improving our knowledge of morphological, physiological, and ecological features of bacterial taxa. DNA–DNA hybridization, which has been used for many years, is still considered the golden standard for bacteria species identification. PCR-based methods investigating 16S rRNA gene sequences, and other approaches, such as the metagenome, have been used to study the physiology and diversity of bacteria and to identify novel genes with potential pharmaceutical and other biotechnological applications. We examined the advantages and limitations of molecular methods currently used to analyze bacterial diversity; these are mainly based on the 16S rRNA gene. These methods have allowed us to examine microorganisms that cannot be cultivated by routine methods and have also been useful for phylogenetic studies. We also considered the importance of improvements in microbe culture techniques and how we can combine different methods to allow a more appropriate assessment of bacterial diversity and to determine their real potential for industrial applications.     

DNA microarray technology: a new tool for the epidemiological typing of bacterial pathogens?

Genomic hybridization on whole genome arrays detects the presence or absence of similar DNA regions in sufficiently related microorganisms, allowing genome-wide comparison of their genetic contents. A whole genome array is based on a sequenced bacterial isolate, and is a collection of DNA probes fixed on a solid support. In a single hybridization experiment, the absence/presence status of all genes of the sequenced microbe in the queried isolate can be examined. The objective of this minireview is to summarize the past usage of DNA microarray technology for microbial strain characterizations, and to estimate its future utilization in epidemiological studies and molecular typing of bacterial pathogens. The studies reviewed here confirm the usefulness of microarray technology for the detection of genetic polymorphisms. However, the construction or purchase of DNA microarrays and the performance of strain to strain hybridization experiments are still prohibitively expensive for routine application. Future use of arrays in epidemiology is likely to depend on the development of more cost-effective protocols, more robust and simplified formats, and the adequate evaluation of their performance (efficacy) and convenience (efficiency) compared with other genotyping methods. It seems more likely that a more focused assay, concentrating on genomic regions of variability previously detected by genome-wide microarrays, will find broad application in routine bacterial epidemiology.     

Beyond bacteria: a study of the enteric microbial consortium in extremely low birth weight infants

Extremely low birth weight (ELBW) infants have high morbidity and mortality, frequently due to invasive infections from bacteria, fungi, and viruses. The microbial communities present in the gastrointestinal tracts of preterm infants may serve as a reservoir for invasive organisms and remain poorly characterized. We used deep pyrosequencing to examine the gut-associated microbiome of 11 ELBW infants in the first postnatal month, with a first time determination of the eukaryote microbiota such as fungi and nematodes, including bacteria and viruses that have not been previously described. Among the fungi observed, Candida sp. and Clavispora sp. dominated the sequences, but a range of environmental molds were also observed. Surprisingly, seventy-one percent of the infant fecal samples tested contained ribosomal sequences corresponding to the parasitic organism Trichinella. Ribosomal DNA sequences for the roundworm symbiont Xenorhabdus accompanied these sequences in the infant with the greatest proportion of Trichinella sequences. When examining ribosomal DNA sequences in aggregate, Enterobacteriales, Pseudomonas, Staphylococcus, and Enterococcus were the most abundant bacterial taxa in a low diversity bacterial community (mean Shannon-Weaver Index of 1.02 ± 0.69), with relatively little change within individual infants through time. To supplement the ribosomal sequence data, shotgun sequencing was performed on DNA from multiple displacement amplification (MDA) of total fecal genomic DNA from two infants. In addition to the organisms mentioned previously, the metagenome also revealed sequences for gram positive and gram negative bacteriophages, as well as human adenovirus C. Together, these data reveal surprising eukaryotic and viral microbial diversity in ELBW enteric microbiota dominated bytypes of bacteria known to cause invasive disease in these infants.     

一起沙门菌引起的食源性疾病爆发的溯源分析

目的:对一起沙门菌引起的食源性疾病爆发进行溯源分析。方法:采用GB4789法对采集的样品进行分离及鉴定,采用16S r RNA基因分型方法及PFGE分型方法对分离的菌株进行分子生物学分析,并对爆发进行溯源分析。结果:生化及血清学结果表明,该起爆发分离的菌型为伦敦沙门氏菌。16S r RNA基因分型表明爆发所分离的菌株均为肠道沙门菌肠道亚种,菌株12 sam与其他4个菌株分子发育距离较远,均为16S r RNA基因分型的TYPE1-11型;PFGE分型结果表明菌株10 sam、16 sam、27 sam及29sam的PFGE带型相似度为100%,菌株12sam跟其他菌株相似率为96%。结论:GB4789法结果表明该起爆发是由伦敦沙门氏菌引起的,16S r RNA基因分型及PFGE分型方法的结果均表明该起食源性疾病来源一致。     

Efficient oligonucleotide probe selection for pan-genomic tiling arrays

Background  

Array comparative genomic hybridization is a fast and cost-effective method for detecting, genotyping, and comparing the genomic sequence of unknown bacterial isolates. This method, as with all microarray applications, requires adequate coverage of probes targeting the regions of interest. An unbiased tiling of probes across the entire length of the genome is the most flexible design approach. However, such a whole-genome tiling requires that the genome sequence is known in advance. For the accurate analysis of uncharacterized bacteria, an array must query a fully representative set of sequences from the species' pan-genome. Prior microarrays have included only a single strain per array or the conserved sequences of gene families. These arrays omit potentially important genes and sequence variants from the pan-genome.     

Molecular epidemiology in current times

Motivated to find options for prevention or intervention, molecular epidemiology aims to identify the host and microbial factors that determine the transmission, manifestation and progression of infectious disease. The genotyping of cultivatable bacterial strains is performed by either anonymous fingerprinting techniques or sequence-based exploration of variable genomic sites. Multilocus sequence typing of housekeeping genes and allele profiling of the core genome have become standard techniques of bacterial strain typing that may be supplemented by whole genome sequencing to explore all single nucleotide variants and/or the composition of the accessory genome. Next, novel protocols to investigate host and microbiome based upon smart third generation sequencing technologies are being developed for an effective surveillance, rapid diagnosis and real-time tracking of infectious diseases.     

A system for the construction of targeted unmarked gene deletions in the genus Burkholderia

Burkholderia species are extremely multidrug resistant, environmental bacteria with extraordinary bioremediation and biocontrol properties. At the same time, these bacteria cause serious opportunistic infections in vulnerable patient populations while some species can potentially be used as bioweapons. The complete DNA sequence of more than 10 Burkholderia genomes provides an opportunity to apply functional genomics to a collection of widely adaptable environmental bacteria thriving in diverse niches and establishing both symbiotic and pathogenic associations with many different organisms. However, extreme multidrug resistance hampers genetic manipulations in Burkholderia. We have developed and evaluated a mutagenesis system based on the homing endonuclease I-SceI to construct targeted, non-polar unmarked gene deletions in Burkholderia. Using the cystic fibrosis pathogen Burkholderia cenocepacia K56-2 as a model strain, we demonstrate this system allows for clean deletions of one or more genes within an operon and also the introduction of multiple deletions in the same strain. We anticipate this tool will have widespread environmental and biomedical applications, facilitating functional genomic studies and construction of safe strains for bioremediation and biocontrol, as well as clinical applications such as live vaccines for Burkholderia and other Gram-negative bacterial species.     

Bacterial phylogenetic clusters revealed by genome structure.

Current bacterial taxonomy is mostly based on phenotypic criteria, which may yield misleading interpretations in classification and identification. As a result, bacteria not closely related may be grouped together as a genus or species. For pathogenic bacteria, incorrect classification or misidentification could be disastrous. There is therefore an urgent need for appropriate methodologies to classify bacteria according to phylogeny and corresponding new approaches that permit their rapid and accurate identification. For this purpose, we have devised a strategy enabling us to resolve phylogenetic clusters of bacteria by comparing their genome structures. These structures were revealed by cleaving genomic DNA with the endonuclease I-CeuI, which cuts within the 23S ribosomal DNA (rDNA) sequences, and by mapping the resulting large DNA fragments with pulsed-field gel electrophoresis. We tested this experimental system on two representative bacterial genera: Salmonella and Pasteurella. Among Salmonella spp., I-CeuI mapping revealed virtually indistinguishable genome structures, demonstrating a high degree of structural conservation. Consistent with this, 16S rDNA sequences are also highly conserved among the Salmonella spp. In marked contrast, the Pasteurella strains have very different genome structures among and even within individual species. The divergence of Pasteurella was also reflected in 16S rDNA sequences and far exceeded that seen between Escherichia and Salmonella. Based on this diversity, the Pasteurella haemolytica strains we analyzed could be divided into 14 phylogenetic groups and the Pasteurella multocida strains could be divided into 9 groups. If criteria for defining bacterial species or genera similar to those used for Salmonella and Escherichia coli were applied, the striking phylogenetic diversity would allow bacteria in the currently recognized species of P. multocida and P. haemolytica to be divided into different species, genera, or even higher ranks. On the other hand, strains of Pasteurella ureae and Pasteurella pneumotropica are very similar to those of P. multocida in both genome structure and 16S rDNA sequence and should be regarded as strains within this species. We conclude that large-scale genome structure can be a sensitive indicator of phylogenetic relationships and that, therefore, I-CeuI-based genomic mapping is an efficient tool for probing the phylogenetic status of bacteria.     

A unique DNA repair and recombination gene (recN) sequence for identification and intraspecific molecular typing of bacterial wilt pathogen Ralstonia solanacearum and its comparative analysis with ribosomal DNA sequences

Ribosomal gene sequences are a popular choice for identification of bacterial species and, often, for making phylogenetic interpretations. Although very popular, the sequences of 16S rDNA and 16-23S intergenic sequences often fail to differentiate closely related species of bacteria. The availability of complete genome sequences of bacteria, in the recent years, has accelerated the search for new genome targets for phylogenetic interpretations. The recently published full genome data of nine strains of R. solanacearum, which causes bacterial wilt of crop plants, has provided enormous genomic choices for phylogenetic analysis in this globally important plant pathogen. We have compared a gene candidate recN, which codes for DNA repair and recombination function, with 16S rDNA/16-23S intergenic ribosomal gene sequences for identification and intraspecific phylogenetic interpretations in R. solanacearum. recN gene sequence analysis of R. solanacearum revealed subgroups within phylotypes (or newly proposed species within plant pathogenic genus, Ralstonia), indicating its usefulness for intraspecific genotyping. The taxonomic discriminatory power of recN gene sequence was found to be superior to ribosomal DNA sequences. In all, the recN-sequence-based phylogenetic tree generated with the Bayesian model depicted 21 haplotypes against 15 and 13 haplotypes obtained with 16S rDNA and 16-23S rDNA intergenic sequences, respectively. Besides this, we have observed high percentage of polymorphic sites (S 23.04%), high rate of mutations (Eta 276) and high codon bias index (CBI 0.60), which makes the recN an ideal gene candidate for intraspecific molecular typing of this important plant pathogen.     

Linkage analysis by two-dimensional DNA typing.

In two-dimensional (2-D) DNA typing, genomic DNA fragments are separated, first according to size by electrophoresis in a neutral polyacrylamide gel and second according to sequence by denaturing gradient gel electrophoresis, followed by hybridization analysis using micro- and minisatellite core probes. The 2-D DNA typing method generates a large amount of information on polymorphic loci per gel. Here we demonstrate the potential usefulness of 2-D DNA typing in an empirical linkage study on the red factor in cattle, and we show an example of the 2-D DNA typing analysis of a human pedigree. The power efficiency of 2-D DNA typing in general is compared with that of single-locus typing by simulation. The results indicate that, although 2-D DNA typing is very efficient in generating data on polymorphic loci, its power to detect linkage is lower than single-locus typing, because it is not obvious whether a spot represents the presence of one or two alleles. It is possible to compensate for this lower informativeness by increasing the sample size. Genome scanning by 2-D DNA typing has the potential to be more efficient than current genotyping methods in scoring polymorphic loci. Hence, it could become a method of choice in mapping genetic traits in humans and animals.     

HLA Typing for the Next Generation

Allele-level resolution data at primary HLA typing is the ideal for most histocompatibility testing laboratories. Many high-throughput molecular HLA typing approaches are unable to determine the phase of observed DNA sequence polymorphisms, leading to ambiguous results. The use of higher resolution methods is often restricted due to cost and time limitations. Here we report on the feasibility of using Pacific Biosciences’ Single Molecule Real-Time (SMRT) DNA sequencing technology for high-resolution and high-throughput HLA typing. Seven DNA samples were typed for HLA-A, -B and -C. The results showed that SMRT DNA sequencing technology was able to generate sequences that spanned entire HLA Class I genes that allowed for accurate allele calling. Eight novel genomic HLA class I sequences were identified, four were novel alleles, three were confirmed as genomic sequence extensions and one corrected an existing genomic reference sequence. This method has the potential to revolutionize the field of HLA typing. The clinical impact of achieving this level of resolution HLA typing data is likely to considerable, particularly in applications such as organ and blood stem cell transplantation where matching donors and recipients for their HLA is of utmost importance.     

SNP-specific extraction of haplotype-resolved targeted genomic regions

The availability of genotyping platforms for comprehensive genetic analysis of complex traits has resulted in a plethora of studies reporting the association of specific single-nucleotide polymorphisms (SNPs) with common diseases or drug responses. However, detailed genetic analysis of these associated regions that would correlate particular polymorphisms to phenotypes has lagged. This is primarily due to the lack of technologies that provide additional sequence information about genomic regions surrounding specific SNPs, preferably in haploid form. Enrichment methods for resequencing should have the specificity to provide DNA linked to SNPs of interest with sufficient quality to be used in a cost-effective and high-throughput manner. We describe a simple, automated method of targeting specific sequences of genomic DNA that can directly be used in downstream applications. The method isolates haploid chromosomal regions flanking targeted SNPs by hybridizing and enzymatically elongating oligonucleotides with biotinylated nucleotides based on their selective binding to unique sequence elements that differentiate one allele from any other differing sequence. The targeted genomic region is captured by streptavidin-coated magnetic particles and analyzed by standard genotyping, sequencing or microarray analysis. We applied this technology to determine contiguous molecular haplotypes across a ~150 kb genomic region of the major histocompatibility complex.     

Minisatellite polymorphism as a tool to distinguish closely related Lactococcus lactis strains

Genome plasticity is considered as a means for bacteria to adapt to their environment. Plasticity in tandem repeat sequences on bacterial genomes has been recently exploited to trace the epidemiology of pathogens. Here, we examine the utility of minisatellite (i.e., a repeat unit of six nucleotides or more) typing in non-pathogenic food bacteria of the species Lactococcus lactis. Thirty-four minisatellites identified on the sequenced L. lactis ssp. lactis strain IL1403 genome were first analyzed in 10 closely related ssp. lactis strains, as determined by randomly amplified polymorphic DNA (RAPD). The selected tandem repeats varied in length, percent identity between repeats, and locations. We showed that: (i) the greatest polymorphism was in orfs encoding exported proteins or in intergenic regions; (ii) two thirds of minisatellites were little- or non-variable, despite as much as 90% identity between tandem repeats; and (iii) dendrograms based on either RAPD or minisatellite analyses were similar. Seven minisatellites identified in this study are potentially useful for lactococcal typing. We then asked whether tandem repeats in L. lactis were stable upon very long-term (up to two years) storage. Despite large rearrangements previously reported in derivative strains, just one of 10 minisatellites tested underwent an alteration, suggesting that tandem repeat rearrangements probably occur during active DNA replication. We conclude that multiple locus minisatellite analysis can be a valuable tool to follow lactococcal strain diversity.     

DAJIN enables multiplex genotyping to simultaneously validate intended and unintended target genome editing outcomes

Genome editing can introduce designed mutations into a target genomic site. Recent research has revealed that it can also induce various unintended events such as structural variations, small indels, and substitutions at, and in some cases, away from the target site. These rearrangements may result in confounding phenotypes in biomedical research samples and cause a concern in clinical or agricultural applications. However, current genotyping methods do not allow a comprehensive analysis of diverse mutations for phasing and mosaic variant detection. Here, we developed a genotyping method with an on-target site analysis software named Determine Allele mutations and Judge Intended genotype by Nanopore sequencer (DAJIN) that can automatically identify and classify both intended and unintended diverse mutations, including point mutations, deletions, inversions, and cis double knock-in at single-nucleotide resolution. Our approach with DAJIN can handle approximately 100 samples under different editing conditions in a single run. With its high versatility, scalability, and convenience, DAJIN-assisted multiplex genotyping may become a new standard for validating genome editing outcomes.

Genome editing can introduce designed mutations into a target genomic site, but also into unintended off-target sites. DAJIN, a novel nanopore sequencing data analysis tool, identifies and quantifies allele numbers and their mutation patterns, reporting consensus sequences and visualizing mutations in alleles at single-nucleotide resolution.