A new approach for the analysis of bacterial microarray-based Comparative Genomic Hybridization: insights from an empirical study

BACKGROUND: Microarray-based Comparative Genomic Hybridization (M-CGH) has been used to characterize the extensive intraspecies genetic diversity found in bacteria at the whole-genome level. Although conventional microarray analytical procedures have proved adequate in handling M-CGH data, data interpretation using these methods is based on a continuous character model in which gene divergence and gene absence form a spectrum of decreasing gene conservation levels. However, whereas gene divergence may yet be accompanied by retention in gene function, gene absence invariably leads to loss of function. This distinction, if ignored, leads to a loss in the information to be gained from M-CGH data. We present here results from experiments in which two genome-sequenced strains of C. jejuni were compared against each other using M-CGH. Because the gene content of both strains was known a priori, we were able to closely examine the effects of sequence divergence and gene absence on M-CGH data in order to define analytical parameters for M-CGH data interpretation. This would facilitate the examination of the relative effects of sequence divergence or gene absence in comparative genomics analyses of multiple strains of any species for which genome sequence data and a DNA microarray are available. RESULTS: As a first step towards improving the analysis of M-CGH data, we estimated the degree of experimental error in a series of experiments in which identical samples were compared against each other by M-CGH. This variance estimate was used to validate a Log Ratio-based methodology for identification of outliers in M-CGH data. We compared two genome strains by M-CGH to examine the effect of probe/target identity on the Log Ratios of signal intensities using prior knowledge of gene divergence and gene absence to establish Log Ratio thresholds for the identification of absent and conserved genes. CONCLUSION: The results from this empirical study validate the Log Ratio thresholds that have been used in other studies to establish gene divergence/absence. Moreover, the analytical framework presented here enhances the information content derived from M-CGH data by shifting the focus from divergent/absent gene detection to accurate detection of conserved and absent genes. This approach closely aligns the technical limitations of M-CGH analysis with practical limitations on the biological interpretation of comparative genomics data

Extensive characterization of Campylobacter jejuni chicken isolates to uncover genes involved in the ability to compete for gut colonization

[À l'origine dans / Was originally part of : Fac. Méd. vétérinaire - Chaire de recherche en salubrité des viandes]Background: Campylobacter jejuni is responsible for human foodborne enteritis. This bacterium is a remarkable colonizer of the chicken gut, with some strains outcompeting others for colonization. To better understand this phenomenon, the objective of this study was to extensively characterize the phenotypic performance of C. jejuni chicken strains and associate their gut colonizing ability with specific genes. Results: C. jejuni isolates (n = 45) previously analyzed for the presence of chicken colonization associated genes were further characterized for phenotypic properties influencing colonization: autoagglutination and chemotaxis as well as adhesion to and invasion of primary chicken caecal cells. This allowed strains to be ranked according to their in vitro performance. After their in vitro capacity to outcompete was demonstrated in vivo, strains were then typed by comparative genomic fingerprinting (CGF). In vitro phenotypical properties displayed a linear variability among the tested strains. Strains possessing higher scores for phenotypical properties were able to outcompete others during chicken colonization trials. When the gene content of strains was compared, some were associated with different phenotypical scores and thus with different outcompeting capacities. Use of CGF profiles showed an extensive genetic variability among the studied strains and suggested that the outcompeting capacity is not predictable by CGF profile. Conclusion: This study revealed a wide array of phenotypes present in C. jejuni strains, even though they were all recovered from chicken caecum. Each strain was classified according to its in vitro competitive potential and its capacity to compete for chicken gut colonization was associated with specific genes. This study also exposed the disparity existing between genetic typing and phenotypical behavior of C. jejuni strains

Comparative genomic assessment of Multi-Locus Sequence Typing: rapid accumulation of genomic heterogeneity among clonal isolates of Campylobacter jejuni

<p>Abstract</p> <p>Background</p> <p>Multi-Locus Sequence Typing (MLST) has emerged as a leading molecular typing method owing to its high ability to discriminate among bacterial isolates, the relative ease with which data acquisition and analysis can be standardized, and the high portability of the resulting sequence data. While MLST has been successfully applied to the study of the population structure for a number of different bacterial species, it has also provided compelling evidence for high rates of recombination in some species. We have analyzed a set of <it>Campylobacter jejuni </it>strains using MLST and Comparative Genomic Hybridization (CGH) on a full-genome microarray in order to determine whether recombination and high levels of genomic mosaicism adversely affect the inference of strain relationships based on the analysis of a restricted number of genetic loci.</p> <p>Results</p> <p>Our results indicate that, in general, there is significant concordance between strain relationships established by MLST and those based on shared gene content as established by CGH. While MLST has significant predictive power with respect to overall genome similarity of isolates, we also found evidence for significant differences in genomic content among strains that would otherwise appear to be highly related based on their MLST profiles.</p> <p>Conclusion</p> <p>The extensive genomic mosaicism between closely related strains has important implications in the context of establishing strain to strain relationships because it suggests that the exact gene content of strains, and by extension their phenotype, is less likely to be "predicted" based on a small number of typing loci. This in turn suggests that a greater emphasis should be placed on analyzing genes of clinical interest as we forge ahead with the next generation of molecular typing methods.</p

In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence

<p>Abstract</p> <p>Background</p> <p>Many approaches have been used to study the evolution, population structure and genetic diversity of <it>Escherichia coli </it>O157:H7; however, observations made with different genotyping systems are not easily relatable to each other. Three genetic lineages of <it>E. coli </it>O157:H7 designated I, II and I/II have been identified using octamer-based genome scanning and microarray comparative genomic hybridization (mCGH). Each lineage contains significant phenotypic differences, with lineage I strains being the most commonly associated with human infections. Similarly, a clade of hyper-virulent O157:H7 strains implicated in the 2006 spinach and lettuce outbreaks has been defined using single-nucleotide polymorphism (SNP) typing. In this study an <it>in silico </it>comparison of six different genotyping approaches was performed on 19 <it>E. coli </it>genome sequences from 17 O157:H7 strains and single O145:NM and K12 MG1655 strains to provide an overall picture of diversity of the <it>E. coli </it>O157:H7 population, and to compare genotyping methods for O157:H7 strains.</p> <p>Results</p> <p><it>In silico </it>determination of lineage, Shiga-toxin bacteriophage integration site, comparative genomic fingerprint, mCGH profile, novel region distribution profile, SNP type and multi-locus variable number tandem repeat analysis type was performed and a supernetwork based on the combination of these methods was produced. This supernetwork showed three distinct clusters of strains that were O157:H7 lineage-specific, with the SNP-based hyper-virulent clade 8 synonymous with O157:H7 lineage I/II. Lineage I/II/clade 8 strains clustered closest on the supernetwork to <it>E. coli </it>K12 and <it>E. coli </it>O55:H7, O145:NM and sorbitol-fermenting O157 strains.</p> <p>Conclusion</p> <p>The results of this study highlight the similarities in relationships derived from multi-locus genome sampling methods and suggest a "common genotyping language" may be devised for population genetics and epidemiological studies. Future genotyping methods should provide data that can be stored centrally and accessed locally in an easily transferable, informative and extensible format based on comparative genomic analyses.</p

Pan-genome sequence analysis using Panseq: an online tool for the rapid analysis of core and accessory genomic regions

<p>Abstract</p> <p>Background</p> <p>The pan-genome of a bacterial species consists of a core and an accessory gene pool. The accessory genome is thought to be an important source of genetic variability in bacterial populations and is gained through lateral gene transfer, allowing subpopulations of bacteria to better adapt to specific niches. Low-cost and high-throughput sequencing platforms have created an exponential increase in genome sequence data and an opportunity to study the pan-genomes of many bacterial species. In this study, we describe a new online pan-genome sequence analysis program, Panseq.</p> <p>Results</p> <p>Panseq was used to identify <it>Escherichia coli </it>O157:H7 and <it>E. coli </it>K-12 genomic islands. Within a population of 60 <it>E. coli </it>O157:H7 strains, the existence of 65 accessory genomic regions identified by Panseq analysis was confirmed by PCR. The accessory genome and binary presence/absence data, and core genome and single nucleotide polymorphisms (SNPs) of six <it>L. monocytogenes </it>strains were extracted with Panseq and hierarchically clustered and visualized. The nucleotide core and binary accessory data were also used to construct maximum parsimony (MP) trees, which were compared to the MP tree generated by multi-locus sequence typing (MLST). The topology of the accessory and core trees was identical but differed from the tree produced using seven MLST loci. The Loci Selector module found the most variable and discriminatory combinations of four loci within a 100 loci set among 10 strains in 1 s, compared to the 449 s required to exhaustively search for all possible combinations; it also found the most discriminatory 20 loci from a 96 loci <it>E. coli </it>O157:H7 SNP dataset.</p> <p>Conclusion</p> <p>Panseq determines the core and accessory regions among a collection of genomic sequences based on user-defined parameters. It readily extracts regions unique to a genome or group of genomes, identifies SNPs within shared core genomic regions, constructs files for use in phylogeny programs based on both the presence/absence of accessory regions and SNPs within core regions and produces a graphical overview of the output. Panseq also includes a loci selector that calculates the most variable and discriminatory loci among sets of accessory loci or core gene SNPs.</p> <p>Availability</p> <p>Panseq is freely available online at <url>http://76.70.11.198/panseq</url>. Panseq is written in Perl.</p

Genotypic Characterization and Prevalence of Virulence Factors among Canadian \u3ci\u3eEscherichia coli\u3c/i\u3e O157:H7 Strains

In this study, the association between genotypic and selected phenotypic characteristics was examined in a collection of Canadian Escherichia coli O157:H7 strains isolated from humans and cattle in the provinces of Alberta, Ontario, Saskatchewan, and Quebec. In a subset of 69 strains selected on the basis of specific phage types (PTs), a strong correlation between the lineage-specific polymorphism assay (LSPA6) genotype and PT was observed with all strains of PTs 4, 14, 21, 31, 33, and 87 belonging to the LSPA6 lineage I (LSPA6-LI) genotype, while those of PTs 23, 45, 67, and 74 belonged to LSPA6 lineage II (LSPA6-LII) genotypes. This correlation was maintained when additional strains of each PT were tested. E. coli O157:H7 strains with the LSPA6-LI genotype were much more common in the collection than were the LSPA6-LII or lineage I/II (LSPA6-LI/II)-related genotypes (82.6, 11.2, and 5.8%, respectively). Of the strains tested, proportionately more LSPA6-LI than LSPA6-LII genotype strains were isolated from humans (52.7% versus 19.7%) than from cattle (47.8% versus 80.2%). In addition, 96.7% of the LSPA6-LII strains carried the stx2c variant gene, while only 50.0% of LSPA6-LI/II and 2.7% of LSPA6-LI strains carried this gene. LSPA6-LII strains were also significantly more likely to possess the colicin D gene, cda (50.8% versus 23.2%), and have combined resistance to streptomycin, sulfisoxazole, and tetracycline (72.1% versus 0.9%) than were LSPA6-LI strains. The LSPA6 genotype- and PT-related characteristics identified may be important markers of specific ecotypes of E. coli O157:H7 that have unique epidemiological and virulence characteristics

Genotypic Characterization and Prevalence of Virulence Factors among Canadian \u3ci\u3eEscherichia coli\u3c/i\u3e O157:H7 Strains

In this study, the association between genotypic and selected phenotypic characteristics was examined in a collection of Canadian Escherichia coli O157:H7 strains isolated from humans and cattle in the provinces of Alberta, Ontario, Saskatchewan, and Quebec. In a subset of 69 strains selected on the basis of specific phage types (PTs), a strong correlation between the lineage-specific polymorphism assay (LSPA6) genotype and PT was observed with all strains of PTs 4, 14, 21, 31, 33, and 87 belonging to the LSPA6 lineage I (LSPA6-LI) genotype, while those of PTs 23, 45, 67, and 74 belonged to LSPA6 lineage II (LSPA6-LII) genotypes. This correlation was maintained when additional strains of each PT were tested. E. coli O157:H7 strains with the LSPA6-LI genotype were much more common in the collection than were the LSPA6-LII or lineage I/II (LSPA6-LI/II)-related genotypes (82.6, 11.2, and 5.8%, respectively). Of the strains tested, proportionately more LSPA6-LI than LSPA6-LII genotype strains were isolated from humans (52.7% versus 19.7%) than from cattle (47.8% versus 80.2%). In addition, 96.7% of the LSPA6-LII strains carried the stx2c variant gene, while only 50.0% of LSPA6-LI/II and 2.7% of LSPA6-LI strains carried this gene. LSPA6-LII strains were also significantly more likely to possess the colicin D gene, cda (50.8% versus 23.2%), and have combined resistance to streptomycin, sulfisoxazole, and tetracycline (72.1% versus 0.9%) than were LSPA6-LI strains. The LSPA6 genotype- and PT-related characteristics identified may be important markers of specific ecotypes of E. coli O157:H7 that have unique epidemiological and virulence characteristics

Development of a comparative genomic fingerprinting assay for rapid and high resolution genotyping of Arcobacter butzleri

Sherpa Romeo green journal. Open access, distributed under the terms of the Creative Commons Attribution (CC-BY) License.Background Molecular typing methods are critical for epidemiological investigations, facilitating disease outbreak detection and source identification. Study of the epidemiology of the emerging human pathogen Arcobacter butzleri is currently hampered by the lack of a subtyping method that is easily deployable in the context of routine epidemiological surveillance. In this study we describe a comparative genomic fingerprinting (CGF) method for high-resolution and high-throughput subtyping of A. butzleri. Comparative analysis of the genome sequences of eleven A. butzleri strains, including eight strains newly sequenced as part of this project, was employed to identify accessory genes suitable for generating unique genetic fingerprints for high-resolution subtyping based on gene presence or absence within a strain. Results A set of eighty-three accessory genes was used to examine the population structure of a dataset comprised of isolates from various sources, including human and non-human animals, sewage, and river water (n=156). A streamlined assay (CGF40) based on a subset of 40 genes was subsequently developed through marker optimization. High levels of profile diversity (121 distinct profiles) were observed among the 156 isolates in the dataset, and a high Simpson’s Index of Diversity (ID) observed (ID > 0.969) indicate that the CGF40 assay possesses high discriminatory power. At the same time, our observation that 115 isolates in this dataset could be assigned to 29 clades with a profile similarity of 90% or greater indicates that the method can be used to identify clades comprised of genetically similar isolates. Conclusions The CGF40 assay described herein combines high resolution and repeatability with high throughput for the rapid characterization of A. butzleri strains. This assay will facilitate the study of the population structure and epidemiology of A. butzleri.Ye

Diabetes Causes Dysfunctional Dopamine Neurotransmission Favoring Nigrostriatal Degeneration in Mice

This article also appears in: Special Collection: COVID-19 Resources.[Background]: Numerous studies indicate an association between neurodegenerative and metabolic diseases. Although still a matter of debate, growing evidence from epidemiological and animal studies indicate that preexisting diabetes increases the risk to develop Parkinson's disease. However, the mechanisms of such an association are unknown.[Objectives]: We investigated whether diabetes alters striatal dopamine neurotransmission and assessed the vulnerability of nigrostriatal neurons to neurodegeneration.[Methods]: We used streptozotocin‐treated and genetically diabetic db/db mice. Expression of oxidative stress and nigrostriatal neuronal markers and levels of dopamine and its metabolites were monitored. Dopamine release and uptake were assessed using fast‐scan cyclic voltammetry. 6‐Hydroxydopamine was unilaterally injected into the striatum using stereotaxic surgery. Motor performance was scored using specific tests.[Results]: Diabetes resulted in oxidative stress and decreased levels of dopamine and its metabolites in the striatum. Levels of proteins regulating dopamine release and uptake, including the dopamine transporter, the Girk2 potassium channel, the vesicular monoamine transporter 2, and the presynaptic vesicle protein synaptobrevin‐2, were decreased in diabetic mice. Electrically evoked levels of extracellular dopamine in the striatum were enhanced, and altered dopamine uptake was observed. Striatal microinjections of a subthreshold dose of the neurotoxin 6‐hydroxydopamine in diabetic mice, insufficient to cause motor alterations in nondiabetic animals, resulted in motor impairment, higher loss of striatal dopaminergic axons, and decreased neuronal cell bodies in the substantia nigra.[Conclusions]: Our results indicate that diabetes promotes striatal oxidative stress, alters dopamine neurotransmission, and increases vulnerability to neurodegenerative damage leading to motor impairment.Funded by the Spanish Ministries of Economy and Competitiveness (grants BFU2014‐52149‐R and BFU2017‐89336‐R to M.V., SAF2016‐78207‐R and PCIN‐2015‐098 to R.M., and BFU2017‐88393‐P to E.D.M.) and of Health, Social Services and Equality (PNSD‐2016I033 to R.M.) and by the Ramón Areces Foundation (ref. 172275 to R.M.). Supported by Medical Research Council UK iCASE award (to S.J.C., R.A.), a Biotechnology and Biological Sciences Research Council UK studentship (to S.V.M.), and Parkinson's UK (J‐1403 to S.J.C.). Partially supported by FEDER funds. CIBERDEM and CIBERNED are initiatives of the Instituto de Salud Carlos III. I.P.T. was supported by a fellowship from the Spanish Ministry of Education, Culture and Sports (FPU 14/04457).Peer reviewe

Comparative genomics profiling of clinical isolates of Aeromonas salmonicida using DNA microarrays

BACKGROUND: Aeromonas salmonicida has been isolated from numerous fish species and shows wide variation in virulence and pathogenicity. As part of a larger research program to identify virulence genes and candidates for vaccine development, a DNA microarray was constructed using a subset of 2024 genes from the draft genome sequence of A. salmonicida subsp. salmonicida strain A449. The microarray included genes encoding known virulence-associated factors in A. salmonicida and homologs of virulence genes of other pathogens. We used microarray-based comparative genomic hybridizations (M-CGH) to compare selected A. salmonicida sub-species and other Aeromonas species from different hosts and geographic locations. RESULTS: Results showed variable carriage of virulence-associated genes and generally increased variation in gene content across sub-species and species boundaries. The greatest variation was observed among genes associated with plasmids and transposons. There was little correlation between geographic region and degree of variation for all isolates tested. CONCLUSION: We have used the M-CGH technique to identify subsets of conserved genes from amongst this set of A. salmonicida virulence genes for further investigation as potential vaccine candidates. Unlike other bacterial characterization methods that use a small number of gene or DNA-based functions, M-CGH examines thousands of genes and/or whole genomes and thus is a more comprehensive analytical tool for veterinary or even human health research