Development of a miniaturised microarray-based assay for the rapid identification of antimicrobial resistance genes in Gram-negative bacteria

We describe the development of a miniaturised microarray for the detection of antimicrobial resistance genes in Gram-negative bacteria. Included on the array are genes encoding resistance to aminoglycosides, trimethoprim, sulphonamides, tetracyclines and beta-lactams, including extended-spectrum beta-lactamases. Validation of the array with control strains demonstrated a 99% correlation between polymerase chain reaction and array results. There was also good correlation between phenotypic and genotypic results for a large panel of Escherichia coli and Salmonella isolates. Some differences were also seen in the number and type of resistance genes harboured by E. coli and Salmonella strains. The array provides an effective, fast and simple method for detection of resistance genes in clinical isolates suitable for use in diagnostic laboratories, which in future will help to understand the epidemiology of isolates and to detect gene linkage in bacterial populations. (C) 2008 Published by Elsevier B.V. and the International Society of Chemotherapy.

Genetic diversity among Escherichia coli O157 : H7 isolates and identification of genes linked to human infections

An Escherichia coli oligonucleotide microarray based on three sequenced genomes was validated for comparative genomic microarray hybridization and used to study the diversity of E. coli O157 isolates from human infections and food and animal sources. Among 26 test strains, 24 (including both Shiga toxin [Stx]-positive and -negative strains) were found to be related to the two sequenced E. coli O157:117 strains, EDL933 and Sakai. However, these strains showed much greater genetic diversity than those reported previously, and most of them could not be categorized as either lineage I or H. Some genes were found more often in isolates from human than from nonhuman sources; e.g., ECs1202 and ECs2976, associated with stx2AB and stx1AB, were in all isolates from human sources but in only 40% of those from nonhuman sources. Some (but not all) lineage I-specific or -dominant genes were also more frequently associated with isolates from human. The results suggested that it might be more effective to concentrate our efforts on finding markers that are directly related to infection rather than those specific to certain lineages. In addition, two Stx-negative O157 cattle isolates (one confirmed to be 117) were significantly different from other Stx-positive and -negative E. coli O157:117 strains and were more similar to MG1655 in their gene content. This work demonstrates that not all E. coli O157:117 strains belong to the same clonal group, and those that were similar to E. coli K-12 might be less virulent.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background: Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity. Results: The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data. Conclusion: After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.

Epidemic multidrug-resistant (MDR-AmpC) Salmonella enterica serovar Newport strains contain three phage regions and a MDR resistance plasmid

Multidrug-resistant (MDR-AmpC) Salmonella enterica serovar Newport has caused serious disease in animals and humans in North America, whereas in the UK S. enterica serovar Newport is not associated with severe disease and usually sensitive to antibiotics; MDR S. Newport (not AmpC) strains have only been isolated from poultry. We found that UK poultry strains belonged to MLST type ST166 and were distinct from cattle isolates for being able to utilize D-tagotose and when compared by pulsed-field gel electrophoresis (PFGE), comparative genomic hybridization (CGH) and diversity arrays technology (DArT). Cattle strains belonged to the ST45 complex differing from ST166 at all seven loci. PFGE showed that 19 out of 27 cattle isolates were more than 85% similar to each other and some UK and US strains were indistinguishable. Both CGH and DArT identified genes (including phage-related ones) that were uniquely present in the US isolates and two such genes identified by DArT showed sequence similarities with the pertussis-like (artAB) toxin. This work demonstrates that MDR-AmpC S. Newport from the USA are genetically closely related to pan-susceptible strains from the UK, but contained three extra phage regions and a MDR plasmid.

Identification of genetic and phenotypic differences associated with prevalent and non-prevalent Salmonella Enteritidis phage types: analysis of variation in amino acid transport

In this study, differences at the genetic level of 37 Salmonella Enteritidis strains from five phage types (PTs) were compared using comparative genomic hybridization (CGH) to assess differences between PTs. There were approximately 400 genes that differentiated prevalent (4, 6, 8 and 13a) and sporadic (11) PTs, of which 35 were unique to prevalent PTs, including six plasmid-borne genes, pefA, B, C, D, srgC and rck, and four chromosomal genes encoding putative amino acid transporters. Phenotype array studies also demonstrated that strains from prevalent PTs were less susceptible to urea stress and utilized L-histidine, L-glutamine, L-proline, L-aspartic acid, gly-asn and gly-gln more efficiently than PT11 strains. Complementation of a PT11 strain with the transporter genes from PT4 resulted in a significant increase in utilization of the amino acids and reduced susceptibility to urea stress. In epithelial cell association assays, PT11 strains were less invasive than other prevalent PTs. Most strains from prevalent PTs were better biofilm formers at 37 degrees C than at 28 degrees C, whilst the converse was true for PT11 strains. Collectively, the results indicate that genetic and corresponding phenotypic differences exist between strains of the prevalent PTs 4, 6, 8 and 13a and non-prevalent PT11 strains that are likely to provide a selective advantage for strains from the former PTs and could help them to enter the food chain and cause salmonellosis.

Use of a LightCycler gyrA mutation assay for rapid identification of mutations conferring decreased susceptibility to ciprofloxacin in multiresistant Salmonella enterica serotype typhimurium DT104 isolates

A LightCycler-based PCR-hybridization gyrA mutation assay (GAMA) was developed to rapidly detect gyrA point mutations in multiresistant (MR) Salmonella enterica serotype Typhimurium DT104 with decreased susceptibility to ciprofloxacin (MIC, 0.25 to 1.0 mg/liter). Ninety-two isolates (49 human, 43 animal) were tested with three individual oligonucleotide probes directed against an Asp-87-to-Asn (GAC --> AAC) mutation, an Asp-87-to-Gly (GAC --> GGC) mutation, and a Ser-83-to-Phe (TCC --> TTC) mutation. Strains homologous to the probes could be distinguished from strains that had different mutations by their probe-target melting temperatures. Thirty-seven human and 30 animal isolates had an Asp-87-to-Asn substitution, 6 human and 6 animal isolates had a Ser-83-to-Phe substitution, and 5 human and 2 animal isolates had an Asp-87-to-Gly substitution. The remaining six strains all had mismatches with the three probes and therefore different gyrA mutations. The sequencing of gyrA from these six isolates showed that one human strain and two animal strains had an Asp-87-to-Tyr (GAC --> TAC) substitution and two animal strains had a Ser-83-to-Tyr (TCC --> TAC) substitution. One animal strain had no gyrA mutation, suggesting that this isolate had a different mechanism of resistance. Fifty-eight of the strains tested were indistinguishable by several different typing methods including antibiograms, pulsed-field gel gel electrophoresis, and plasmid profiling, although they could be further subdivided according to gyrA mutation. This study confirmed that MR DT104 with decreased susceptibility to ciprofloxacin from humans and food animals in England and Wales may have arisen independently against a background of clonal spread of MR DT104.

The tagged microarray marker (TAM) method for high-throughput detection of single nucleotide and indel polymorphisms

The tagged microarray marker (TAM) method allows high-throughput differentiation between predicted alternative PCR products. Typically, the method is used as a molecular marker approach to determining the allelic states of single nucleotide polymorphisms (SNPs) or insertion-deletion (indel) alleles at genomic loci in multiple individuals. Biotin-labeled PCR products are spotted, unpurified, onto a streptavidin-coated glass slide and the alternative products are differentiated by hybridization to fluorescent detector oligonucleotides that recognize corresponding allele-specific tags on the PCR primers. The main attractions of this method are its high throughput (thousands of PCRs are analyzed per slide), flexibility of scoring (any combination, from a single marker in thousands of samples to thousands of markers in a single sample, can be analyzed) and flexibility of scale (any experimental scale, from a small lab setting up to a large project). This protocol describes an experiment involving 3,072 PCRs scored on a slide. The whole process from the start of PCR setup to receiving the data spreadsheet takes 2 d.

Identification of core and variable components of the Salmonella enterica subspecies I genome by microarray

We have performed microarray hybridization studies on 40 clinical isolates from 12 common serovars within Salmonella enterica subspecies I to identify the conserved chromosomal gene pool. We were able to separate the core invariant portion of the genome by a novel mathematical approach using a decision tree based on genes ranked by increasing variance. All genes within the core component were confirmed using available sequence and microarray information for S. enterica subspecies I strains. The majority of genes within the core component had conserved homologues in Escherichia coli K-12 strain MG1655. However, many genes present in the conserved set which were absent or highly divergent in K-12 had close homologues in pathogenic bacteria such as Shigella flexneri and Pseudomonas aeruginosa. Genes within previously established virulence determinants such as SPI1 to SPI5 were conserved. In addition several genes within SPI6, all of SPI9, and three fimbrial operons (fim, bcf, and stb) were conserved within all S. enterica strains included in this study. Although many phage and insertion sequence elements were missing from the core component, approximately half the pseudogenes present in S. enterica serovar Typhi were conserved. Furthermore, approximately half the genes conserved in the core set encoded hypothetical proteins. Separation of the core and variant gene sets within S. enterica subspecies I has offered fundamental biological insight into the genetic basis of phenotypic similarity and diversity across S. enterica subspecies I and shown how the core genome of these pathogens differs from the closely related E. coli K-12 laboratory strain.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

High throughput profile-profile based fold recognition for the entire human proteome

BACKGROUND: In order to maintain the most comprehensive structural annotation databases we must carry out regular updates for each proteome using the latest profile-profile fold recognition methods. The ability to carry out these updates on demand is necessary to keep pace with the regular updates of sequence and structure databases. Providing the highest quality structural models requires the most intensive profile-profile fold recognition methods running with the very latest available sequence databases and fold libraries. However, running these methods on such a regular basis for every sequenced proteome requires large amounts of processing power.In this paper we describe and benchmark the JYDE (Job Yield Distribution Environment) system, which is a meta-scheduler designed to work above cluster schedulers, such as Sun Grid Engine (SGE) or Condor. We demonstrate the ability of JYDE to distribute the load of genomic-scale fold recognition across multiple independent Grid domains. We use the most recent profile-profile version of our mGenTHREADER software in order to annotate the latest version of the Human proteome against the latest sequence and structure databases in as short a time as possible. RESULTS: We show that our JYDE system is able to scale to large numbers of intensive fold recognition jobs running across several independent computer clusters. Using our JYDE system we have been able to annotate 99.9% of the protein sequences within the Human proteome in less than 24 hours, by harnessing over 500 CPUs from 3 independent Grid domains. CONCLUSION: This study clearly demonstrates the feasibility of carrying out on demand high quality structural annotations for the proteomes of major eukaryotic organisms. Specifically, we have shown that it is now possible to provide complete regular updates of profile-profile based fold recognition models for entire eukaryotic proteomes, through the use of Grid middleware such as JYDE.

On the structural differences between markers and genomic AC microsatellites

AC microsatellites have proved particularly useful as genetic markers. For some purposes, such as in population biology, the inferences drawn depend on the quantitative values of their mutation rates. This, together with intrinsic biological interest, has led to widespread study of microsatellite mutational mechanisms. Now, however, inconsistencies are appearing in the results of marker-based versus non-marker-based studies of mutational mechanisms. The reasons for this have not been investigated, but one possibility, pursued here, is that the differences result from structural differences between markers and genomic microsatellites. Here we report a comparison between the CEPH AC marker microsatellites and the global population of AC microsatellites in the human genome. AC marker microsatellites are longer than the global average. Controlling for length, marker microsatellites contain on average fewer interruptions, and have longer segments, than their genomic counterparts. Related to this, marker microsatellites show a greater tendency to concentrate the majority of their repeats into one segment. These differences plausibly result from scientists selecting markers for their high polymorphism. In addition to the structural differences, there are differences in the base composition of flanking sequences, marker flanking regions being richer in C and G and poorer in A and T. Our results indicate that there are profound differences between marker and genomic microsatellites that almost certainly affect their mutation rates. There is a need for a unified model of mutational mechanisms that accounts for both marker-derived and genomic observations. A suggestion is made as to how this might be done.

The chromosomal assessment of salt tolerant substituted tritipyrum using genomic fluorescent in situ hybridiization (FISH)

Wheat, although moderately tolerant to salt, can not be cultivated in many areas. However, in the triticeae tribe, some of the wild wheat relatives are highly tolerant, e.g. Thinopyrum bessarabicum, which grows on the sea shore. Eight primary hexaploid tritipyrum lines, amphiploids between Triticum durum and Thinopyrum bessarabicum have been produced which can set seed in at least 250 mM NaCl. These tritipyrums (2n=6x=42, AABBEbEb) due to reasons such as brittle rachis, continuous production of tillers, late maturity, tall stature and meiotic instability will not fulfill the requirements of a successful commercial salt tolerant crop. To overcome such problems the substituted tritipyrum, in which selected Eb chromosomes are replaced by D genome chromosomes of 6x wheat, was produced from 6x tritipyrum x 6x wheat hybrids (F1: 2n=6x=42, AABBDEb) followed by selfing and backcrossing with 6x tritipyrum. The fertile plants among the above progenies were screened by the genomic fluorescent in situ hybridization technique to identify their Eb and D chromosome constitution. This study showed that producing tritiprum with variable numbers of Eb and D genome chromosomes is feasible and that FISH is a useful technique for determining the number of Eb chromosomes present.

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, P