DNA signature-based approaches for bacterial detection and identification

Assigning the unassigned: A signature-based classification of rDNA metabarcodes reveals new deep-sea diversity

Environmental DNA metabarcoding reveals a vast genetic diversity of marine eukaryotes. Yet, most of the metabarcoding data remain unassigned due to the paucity of reference databases. This is particularly true for the deep-sea meiofauna and eukaryotic microbiota, whose hidden diversity is largely unexplored. Here, we tackle this issue by using unique DNA signatures to classify unknown metabarcodes assigned to deep-sea foraminifera. We analyzed metabarcoding data obtained from 311 deep-sea sediment samples collected in the Clarion-Clipperton Fracture Zone, an area of potential polymetallic nodule exploitation in the Eastern Pacific Ocean. Using the signatures designed in the 37F hypervariable region of the 18S rRNA gene, we were able to classify 802 unassigned metabarcodes into 61 novel lineages, which have been placed in 27 phylogenetic clades. The comparison of new lineages with other foraminiferal datasets shows that most novel lineages are widely distributed in the deep sea. Five lineages are also present in the shallow-water datasets; however, phylogenetic analysis of these lineages separates deep-sea and shallow-water metabarcodes except in one case. While the signature-based classification does not solve the problem of gaps in reference databases, this taxonomy-free approach provides insight into the distribution and ecology of deep-sea species represented by unassigned metabarcodes, which could be useful in future applications of metabarcoding for environmental monitoring.

Development of a database and standardized approach for rpoB sequence-based subtyping and identification of aerobic spore-forming Bacillales

Aerobic spore-forming Bacillales are a highly diverse and ubiquitous group that includes organisms that cause foodborne illnesses and food spoilage. Classical microbiological and biochemical identification of members of the order Bacillales represents a challenge due to the diversity of organisms in this group as well as the fact that the phenotypic-based taxonomic assignment of some named species in this group is not consistent with their phylogenomic characteristics. DNA-sequencing-based tools, on the other hand, can be fast and cost-effective, and can provide for a more reliable identification and characterization of Bacillales isolates. In comparison to 16S rDNA, rpoB was shown to better discriminate between Bacillales isolates and to allow for improved taxonomic assignment to the species level. However, the lack of a publicly accessible rpoB database, as well as the lack of standardized protocols for rpoB-based typing and strain identification, is a major challenge. Here, we report (i) the curation of a DNA sequence database for rpoB-based subtype classification of Bacillales isolates; (ii) the development of standardized protocols for generating rpoB sequence data, and a scheme for rpoB-based initial taxonomic identification of Bacillales isolates at the species level; and (iii) the integration of the database in a publicly accessible online platform that allows for the analysis of rpoB sequence data from uncharacterized Bacillales isolates. Specifically, we curated a database of DNA sequences for a 632-nt internal variable region within the rpoB gene from representative Bacillales reference type strains and a large number of isolates that we have previously isolated and characterized through multiple projects. As of May 21, 2021, the rpoB database contained more than 8350 rpoB sequences representing 1902 distinct rpoB allelic types that can be classified into 160 different genera. The database also includes 1129 rpoB sequences for representative Bacillales reference type strains as available on May 21, 2021 in the NCBI database. The rpoB database is integrated into the online Food Microbe Tracker platform (www.foodmicrobetracker.com) and can be queried using the integrated BLAST tool to initially subtype and taxonomically identify aerobic and facultative anaerobic spore-formers. While whole-genome sequencing is increasingly used in bacterial taxonomy, the rpoB sequence-based identification scheme described here provides a valuable tool as it allows for rapid and cost-effective initial isolate characterization, which can help to identify and characterize foodborne pathogens and food spoilage bacteria. In addition, the database and primers described here can also be adopted for metagenomics approaches that include rpoB as a target, improving discriminatory power and identification over what can be achieved using 16S rDNA as a target.

Diversity and plant growth promoting properties of rhizobia isolated from root nodules of Ononis arvensis

This is the first report describing isolates from root nodules of Ononis arvensis (field restharrow). The aim of this investigation was to describe the diversity, phylogeny, and plant growth promoting features of microsymbionts of O. arvensis, i.e., a legume plant growing in different places of the southern part of Poland. Twenty-nine bacterial isolates were characterized in terms of their phenotypic properties, genome fingerprinting, and comparative analysis of their 16S rRNA, nodC and acdS gene sequences. Based on the nodC and 16S rRNA gene phylogenies, the O. arvensis symbionts were grouped close to bacteria of the genera Rhizobium and Mesorhizobium, which formed monophyletic clusters. The acdS gene sequences of all the isolates tested exhibited the highest similarities to the corresponding gene sequences of genus Mesorhizobium strains. The presence of the acdS genes in the genomes of rhizobia specific for O. arvensis implies that these bacteria may promote the growth and development of their host plant in stress conditions. The isolated bacteria showed a high genomic diversity and, in the BOX-PCR reaction, all of them (except three) exhibited DNA fingerprints specific only for them. Our studies showed that restharrow isolates formed effective symbiotic interactions with their native host (O. arvensis) and Ononis spinosa but not with Trifolium repens and Medicago sativa belonging to the same tribe Trifolieae as Ononis species and not with Lotus corniculatus, representing the tribe Loteae.

Comparative analyses of phenotypic methods and 16S rRNA, khe, rpoB genes sequencing for identification of clinical isolates of Klebsiella pneumoniae

The present work aimed to evaluate 16S rRNA, khe and rpoB gene sequencing for the identification of Klebsiella pneumoniae in comparison with phenotypic methods. Fifteen clinical isolates were examined, which were initially identified as K. pneumoniae subsp. pneumoniae using the automated VITEK 32 system in two hospitals in Enshi City, China. Their identity was further supported by conventional phenotypic methods on the basis of morphological and biochemical characteristics. Using Bayesian phylogenetic analyses and haplotypes network reconstruction, 13 isolates were identified as K. pneumoniae, whereas the other two isolates (K19, K24) were classified as Shigella sp. and Enterobacter sp., respectively. Of the three genes, 16S rRNA and khe gene could discriminate the clinical isolates at the genus level, whereas rpoB could discriminate Klebsiella at the species and even subspecies level. Overall, the gene tree based on rpoB is more compatible with the currently accepted classification of Klebsiella than those based on 16S rRNA and khe genes, showing that rpoB can be a powerful tool for identification of K. pneumoniae isolates. Above all, our study challenges the utility of khe as a species-specific marker for identification of K. pneumoniae.

Molecular diagnosis of infectious diseases using cytological specimens

Pathologists have an important role in the diagnosis of infectious disease (ID). In many cases, a definitive diagnosis can be made using cytopathology alone. However, several ancillary techniques can be used on cytological material to reach a specific diagnosis by identifying the causative agent and consequently defining the management of the patient. This review aims to present the effectiveness of the application of molecular studies on cytological material to diagnose IDs and discuss the advantages and disadvantages of the various molecular techniques according to the type of cytological specimen and the infectious agents.

Whole-genome sequencing targets drug-resistant bacterial infections

During the past two decades, the technological progress of whole-genome sequencing (WGS) had changed the fields of Environmental Microbiology and Biotechnology, and, currently, is changing the underlying principles, approaches, and fundamentals of Public Health, Epidemiology, Health Economics, and national productivity. Today’s WGS technologies are able to compete with conventional techniques in cost, speed, accuracy, and resolution for day-to-day control of infectious diseases and outbreaks in clinical laboratories and in long-term epidemiological investigations. WGS gives rise to an exciting future direction for personalized Genomic Epidemiology. One of the most vital and growing public health problems is the emerging and re-emerging of multidrug-resistant (MDR) bacterial infections in the communities and healthcare settings, reinforced by a decline in antimicrobial drug discovery. In recent years, retrospective analysis provided by WGS has had a great impact on the identification and tracking of MDR microorganisms in hospitals and communities. The obtained genomic data are also important for developing novel easy-to-use diagnostic assays for clinics, as well as for antibiotic and therapeutic development at both the personal and population levels. At present, this technology has been successfully applied as an addendum to the real-time diagnostic methods currently used in clinical laboratories. However, the significance of WGS for public health may increase if: (a) unified and user-friendly bioinformatics toolsets for easy data interpretation and management are established, and (b) standards for data validation and verification are developed. Herein, we review the current and future impact of this technology on diagnosis, prevention, treatment, and control of MDR infectious bacteria in clinics and on the global scale.

DNA Sequence Signatures for Rapid Detection of Six Target Bacterial Pathogens Using PCR Assays

Using Streptococcus pyogenes as a model, we previously established a stepwise computational workflow to effectively identify species-specific DNA signatures that could be used as PCR primer sets to detect target bacteria with high specificity and sensitivity. In this study, we extended the workflow for the rapid development of PCR assays targeting Enterococcus faecalis, Enterococcus faecium, Clostridium perfringens, Clostridium difficile, Clostridium tetani, and Staphylococcus aureus, which are of safety concern for human tissue intended for transplantation. Twenty-one primer sets that had sensitivity of detecting 5-50 fg DNA from target bacteria with high specificity were selected. These selected primer sets can be used in a PCR array for detecting target bacteria with high sensitivity and specificity. The workflow could be widely applicable for the rapid development of PCR-based assays for a wide range of target bacteria, including those of biothreat agents.

Improved PCR assay for the species-specific identification and quantitation of Legionella pneumophila in water

Legionellosis outbreak is a major global health care problem. However, current Legionella risk assessments may be compromised by uncertainties in Legionella detection methods, infectious dose, and strain infectivity. These limitations may place public health at significant risk, leading to significant monetary losses in health care. However, there are still unmet needs for its rapid identification and monitoring of legionellae in water systems. Therefore, in the present study, a primer set was designed based on a LysR-type transcriptional regulator (LTTR) family protein gene of Legionella pneumophila subsp. pneumophila str. Philadelphia 1 because it was found that this gene is structurally diverse among species through BLAST searches. The specificity of the primer set was evaluated using genomic DNA from 6 strains of L. pneumophila, 5 type strains of other related Legionella species, and other 29 reference pathogenic bacteria. The primer set used in the PCR assay amplified a 264-bp product for only targeted six strains of L. pneumophila. The assay was also able to detect at least 1.39 × 10(3) copies/μl of cloned amplified target DNA using purified DNA or 7.4 × 10(0) colony-forming unit per reaction when using calibrated cell suspension. In addition, the sensitivity and specificity of this assay were confirmed by successful detection of Legionella pneumophila in environmental water samples.

Validation and Application of a Real-time PCR Protocol for the Specific Detection and Quantification of Clavibacter michiganensis subsp. sepedonicus in Potato

Clavibacter michiganensis subsp. sepedonicus (Cms) multiplies very rapidly, passing through the vascular strands and into the stems and petioles of a diseased potato. Therefore, the rapid and specific detection of this pathogen is highly important for the effective control of the pathogen. Although several PCR assays have been developed for detection, they cannot afford specific detection of Cms. Therefore, in this study, a computational genome analysis was performed to compare the sequenced genomes of the C. michiganensis subspecies and to identify an appropriate gene for the development of a subspecies-specific PCR primer set (Cms89F/R). The specificity of the primer set based on the putative phage-related protein was evaluated using genomic DNA from seven isolates of Cms and 27 other reference strains. The Cms89F/R primer set was more specific and sensitive than the existing assays in detecting Cms in in vitro using Cms cells and its genomic DNA. This assay was also able to detect at least 1.47×10(2) copies/μl of cloned-amplified target DNA, 5 fg of DNA using genomic DNA or 10(-6) dilution point of 0.12 at OD600 units of cells per reaction using a calibrated cell suspension.

Improved PCR assay for the specific detection and quantitation of Escherichia coli serotype O157 in water

Escherichia coli serotype O157 is still a major global healthcare problem. However, only limited information is now available on the molecular and serological detection of pathogenic bacteria. Therefore, the development of appropriate strategies for their rapid identification and monitoring is still needed. In general, the sequence analysis based on stx, slt, eae, hlyA, rfb, and fliCh7 genes is widely employed for the identification of E. coli serotype O157; but there have been critical defects in the diagnosis and identification of E. coli serotype O157, in that they are also present in other E. coli serogroups. In this study, NCBI-BLAST searches using the nucleotide sequences of the putative regulatory protein gene from E. coli O157:H7 str. Sakai found sequence difference at the serotype level. The specific primers from the putative regulatory protein gene were designed and investigated for their sensitivity and specificity for detecting the pathogen in environment water samples. The specificity of the primer set was evaluated using genomic DNA from 8 isolates of E. coli serotype O157 and 32 other reference strains. In addition, the sensitivity and specificity of this assay were confirmed by successful identification of E. coli serotype O157 in environmental water samples. In conclusion, this study showed that the newly developed quantitative serotype-specific PCR method is a highly specific and efficient tool for the surveillance and rapid detection of high-risk E. coli serotype O157.

Rapid TaqMan-based quantification of chlorophyll d-containing cyanobacteria in the genus Acaryochloris

Reports of the chlorophyll (Chl) d-containing cyanobacterium Acaryochloris have accumulated since its initial discovery in 1996. The majority of this evidence is based on amplification of the gene coding for the 16S rRNA, and due to the wide geographical distribution of these sequences, a global distribution of Acaryochloris species was suggested. Here, we present a rapid, reliable, and cost-effective TaqMan-based quantitative PCR (qPCR) assay that was developed for the specific detection of Acaryochloris species in complex environmental samples. The TaqMan probe showed detection limits of ~10 16S rRNA gene copy numbers based on standard curves consisting of plasmid inserts. DNA from five Acaryochloris strains, i.e., MBIC11017, CCMEE5410, HICR111A, CRS, and Awaji-1, exhibited amplification efficiencies of >94% when tested in the TaqMan assay. When used on complex natural communities, the TaqMan assay detected the presence of Acaryochloris species in four out of eight samples of crustose coralline algae (CCA), collected from temperate and tropical regions. In three out of these TaqMan-positive samples, the presence of Chl d was confirmed via high-performance liquid chromatography (HPLC), and corresponding cell estimates of Acaryochloris species amounted to 7.6 × 10(1) to 3.0 × 10(3) per mg of CCA. These numbers indicate a substantial contribution of Chl d-containing cyanobacteria to primary productivity in endolithic niches. The new TaqMan assay allows quick and easy screening of environmental samples for the presence of Acaryochloris species and is an important tool to further resolve the global distribution and significance of this unique oxyphototroph.

Sequencing, characterization, and gene expression analysis of the histidine decarboxylase gene cluster of Morganella morganii

The histidine decarboxylase gene cluster of Morganella morganii DSM30146(T) was sequenced, and four open reading frames, named hdcT1, hdc, hdcT2, and hisRS were identified. Two putative histidine/histamine antiporters (hdcT1 and hdcT2) were located upstream and downstream the hdc gene, codifying a pyridoxal-P dependent histidine decarboxylase, and followed by hisRS gene encoding a histidyl-tRNA synthetase. This organization was comparable with the gene cluster of other known Gram negative bacteria, particularly with that of Klebsiella oxytoca. Recombinant Escherichia coli strains harboring plasmids carrying the M. morganii hdc gene were shown to overproduce histidine decarboxylase, after IPTG induction at 37 °C for 4 h. Quantitative RT-PCR experiments revealed the hdc and hisRS genes were highly induced under acidic and histidine-rich conditions. This work represents the first description and identification of the hdc-related genes in M. morganii. Results support the hypothesis that the histidine decarboxylation reaction in this prolific histamine producing species may play a role in acid survival. The knowledge of the role and the regulation of genes involved in histidine decarboxylation should improve the design of rational strategies to avoid toxic histamine production in foods.

Improved PCR for identification of Pseudomonas aeruginosa

The aim of the present study was to develop a noble and specific marker for a quantitative polymerase chain reaction (PCR) assay for the species-specific detection of Pseudomonas aeruginosa based on the O-antigen acetylase gene. It is an important challenge to characterize populations of the bacterium P. aeruginosa, an opportunist by virtue of its physiological and genetic adaptability. However, molecular and serological methods currently available for sensitive and specific detection of P. aeruginosa are by no means satisfactory because there have been critical defects in the diagnosis and identification of P. aeruginosa strains in that these assays also detect other Pseudomonas species, or do not obtain amplified products from P. aeruginosa strains. Therefore, a primer set was designed based on the O-antigen acetylase gene of P. aeruginosa PA01 because it has been known that this gene is structurally diverse among species. The specificity of the primer set was evaluated using genomic DNA from six isolates of P. aeruginosa, 18 different species of Pseudomonas, and 23 other reference pathogenic bacteria. The primer set used in the PCR assay amplified a 232-bp amplicon for only six P. aeruginosa strains. The assay was also able to detect at least 1.41 × 10(3) copies/μl of cloned amplified target DNA using purified DNA, or 2.7 × 10(2) colony-forming unit per reaction when using calibrated cell suspension. In conclusion, this assay can be applied as a practical diagnostic method for epidemiological research and the sanitary management of water with a low level or latent infection of P. aeruginosa.

Advances in multiplex PCR: balancing primer efficiencies and improving detection success

1. Multiplex PCR is a valuable tool in many biological studies but it is a multifaceted procedure that has to be planned and optimised thoroughly to achieve robust and meaningful results. In particular, primer concentrations have to be adjusted to assure an even amplification of all targeted DNA fragments. Until now, total DNA extracts were used for balancing primer efficiencies; however, the applicability for comparisons between taxa or different multiple-copy genes was limited owing to the unknown number of template molecules present per total DNA.

2. Based on a multiplex system developed to track trophic interactions in high Alpine arthropods, we demonstrate a fast and easy way of generating standardised DNA templates. These were then used to balance the amplification success for the different targets and to subsequently determine the sensitivity of each primer pair in the multiplex PCR.

3. In the current multiplex assay, this approach led to an even amplification success for all seven targeted DNA fragments. Using this balanced multiplex PCR, methodological bias owing to variation in primer efficiency will be avoided when analysing field-derived samples.

4. The approach outlined here allows comparing multiplex PCR sensitivity, independent of the investigated species, genome size or the targeted genes. The application of standardised DNA templates not only makes it possible to optimise primer efficiency within a given multiplex PCR, but it also offers to adjust and/or to compare the sensitivity between different assays. Along with other factors that influence the success of multiplex reactions, and which we discuss here in relation to the presented detection system, the adoption of this approach will allow for direct comparison of multiplex PCR data between systems and studies, enhancing the utility of this assay type.

Evolutionary and experimental assessment of novel markers for detection of Xanthomonas euvesicatoria in plant samples

BACKGROUND

Bacterial spot-causing xanthomonads (BSX) are quarantine phytopathogenic bacteria responsible for heavy losses in tomato and pepper production. Despite the research on improved plant spraying methods and resistant cultivars, the use of healthy plant material is still considered as the most effective bacterial spot control measure. Therefore, rapid and efficient detection methods are crucial for an early detection of these phytopathogens.

METHODOLOGY

In this work, we selected and validated novel DNA markers for reliable detection of the BSX Xanthomonas euvesicatoria (Xeu). Xeu-specific DNA regions were selected using two online applications, CUPID and Insignia. Furthermore, to facilitate the selection of putative DNA markers, a customized C program was designed to retrieve the regions outputted by both databases. The in silico validation was further extended in order to provide an insight on the origin of these Xeu-specific regions by assessing chromosomal location, GC content, codon usage and synteny analyses. Primer-pairs were designed for amplification of those regions and the PCR validation assays showed that most primers allowed for positive amplification with different Xeu strains. The obtained amplicons were labeled and used as probes in dot blot assays, which allowed testing the probes against a collection of 12 non-BSX Xanthomonas and 23 other phytopathogenic bacteria. These assays confirmed the specificity of the selected DNA markers. Finally, we designed and tested a duplex PCR assay and an inverted dot blot platform for culture-independent detection of Xeu in infected plants.

SIGNIFICANCE

This study details a selection strategy able to provide a large number of Xeu-specific DNA markers. As demonstrated, the selected markers can detect Xeu in infected plants both by PCR and by hybridization-based assays coupled with automatic data analysis. Furthermore, this work is a contribution to implement more efficient DNA-based methods of bacterial diagnostics.

Identification of DNA signatures suitable for use in development of real-time PCR assays by whole-genome sequence approaches: use of Streptococcus pyogenes in a pilot study

A stepwise computational approach using three layers of publicly available software was found to effectively identify DNA signatures for Streptococcus pyogenes. PCR testing validated that 9 out of 15 signature-derived primer sets could detect as low as 5 fg of target DNA with high specificity. The selected signature-derived primer sets were successfully evaluated against all 23 clinical isolates. The approach is readily applicable for designing molecular assays for rapid detection and characterization of various pathogenic bacteria.

Identification of Xanthomonas fragariae, Xanthomonas axonopodis pv. phaseoli, and Xanthomonas fuscans subsp. fuscans with novel markers and using a dot blot platform coupled with automatic data analysis

Phytosanitary regulations and the provision of plant health certificates still rely mainly on long and laborious culture-based methods of diagnosis, which are frequently inconclusive. DNA-based methods of detection can circumvent many of the limitations of currently used screening methods, allowing a fast and accurate monitoring of samples. The genus Xanthomonas includes 13 phytopathogenic quarantine organisms for which improved methods of diagnosis are needed. In this work, we propose 21 new Xanthomonas-specific molecular markers, within loci coding for Xanthomonas-specific protein domains, useful for DNA-based methods of identification of xanthomonads. The specificity of these markers was assessed by a dot blot hybridization array using 23 non-Xanthomonas species, mostly soil dwelling and/or phytopathogens for the same host plants. In addition, the validation of these markers on 15 Xanthomonas spp. suggested species-specific hybridization patterns, which allowed discrimination among the different Xanthomonas species. Having in mind that DNA-based methods of diagnosis are particularly hampered for unsequenced species, namely, Xanthomonas fragariae, Xanthomonas axonopodis pv. phaseoli, and Xanthomonas fuscans subsp. fuscans, for which comparative genomics tools to search for DNA signatures are not yet applicable, emphasis was given to the selection of informative markers able to identify X. fragariae, X. axonopodis pv. phaseoli, and X. fuscans subsp. fuscans strains. In order to avoid inconsistencies due to operator-dependent interpretation of dot blot data, an image-processing algorithm was developed to analyze automatically the dot blot patterns. Ultimately, the proposed markers and the dot blot platform, coupled with automatic data analyses, have the potential to foster a thorough monitoring of phytopathogenic xanthomonads.

A GPU-based calculation using the three-dimensional FDTD method for electromagnetic field analysis

Numerical simulations with the numerical human model using the finite-difference time domain (FDTD) method have recently been performed frequently in a number of fields in biomedical engineering. However, the FDTD calculation runs too slowly. We focus, therefore, on general purpose programming on the graphics processing unit (GPGPU). The three-dimensional FDTD method was implemented on the GPU using Compute Unified Device Architecture (CUDA). In this study, we used the NVIDIA Tesla C1060 as a GPGPU board. The performance of the GPU is evaluated in comparison with the performance of a conventional CPU and a vector supercomputer. The results indicate that three-dimensional FDTD calculations using a GPU can significantly reduce run time in comparison with that using a conventional CPU, even a native GPU implementation of the three-dimensional FDTD method, while the GPU/CPU speed ratio varies with the calculation domain and thread block size.

Automatic analysis of macroarrays images

The analysis of dot blot (macroarray) images is currently based on the human identification of positive/negative dots, which is a subjective and time consuming process. This paper presents a system for the automatic analysis of dot blot images, using a pre-defined grid of markers, including a number of ON and OFF controls. The geometric deformations of the input image are corrected, and the individual markers detected, both tasks fully automatically. Based on a previous training stage, the probability for each marker to be ON is established. This information is provided together with quality parameters for training, noise and classification, allowing for a fully automatic evaluation of a dot blot image.

Characterisation of Lactobacillus helveticus strains producing antihypertensive peptides by RAPD and inverse-PCR of IS elements

Lactobacillus helveticus is used for the manufacture of cheeses and milk-based products. Although it is not considered a probiotic microorganism, some strains demonstrated beneficial effects through the production of antihypertensive peptides from the hydrolysis of casein during milk fermentation. Strain-specificity of bioactive peptide production by L. helveticus makes the availability of reliable typing methods essential for both legal and good manufacturing processes. Accordingly, RAPD and inverse-PCR of five insertion sequence elements were comparatively evaluated for the molecular characterisation of four L. helveticus dairy cultures producing antihypertensive peptides and fourteen reference strains. Calculation of discriminatory indices and cluster analysis of the DNA fingerprints confirmed the suitability of both approaches for acceptable strain differentiation. Although RAPD was more discriminating, for a few test strains a neat discrimination was only achieved through multiplex inverse-PCR, thus suggesting the suitability of a combined analytical approach for a finer strain discrimination.

Automatic detection of molecular markers in digital images

Labeled molecular markers are an important tool in molecular biology. This work presents a method for the automatic identification of molecular markers in dot blot images. The method detects the location of markers in the image and their size. An experiment was made with 6 test images, which were used to produce an additional set of 222 images with various rotation, translation, contrast and noise levels. Over 7500 markers were identified automatically and compared with reference values obtained manually. The RMS error for the marker positioning in the original test images were between 1.1 and 3.8 pixels, which is about 1/10 of the typical radius (26 pixels). The method proposed was found to be almost insensitive to grid rotation and translation, and reasonably robust to image contrast changes and presence of noise.