RestrictionDigest: A powerful Perl module for simulating genomic restriction digests

Background: Reduced-representation sequencing technology iswidely used in genotyping for its economical and efficient features. A popular way to construct the reduced-representation sequencing libraries is to digest the genomic DNA with restriction enzymes. A key factor of this method is to determine the restriction enzyme(s). But there are few computer programs which can evaluate the usability of restriction enzymes in reduced-representation sequencing. SimRAD is an R package which can simulate the digestion of DNA sequence by restriction enzymes and return enzyme loci number as well as fragment number. But for linkage mapping analysis, enzyme loci distribution is also an important factor to evaluate the enzyme. For phylogenetic studies, comparison of the enzyme performance across multiple genomes is important. It is strongly needed to develop a simulation tool to implement these functions. Results: Here, we introduce a Perl module named RestrictionDigest with more functions and improved performance. It can analyze multiple genomes at one run and generate concise comparison of enzyme performance across the genomes. It can simulate single-enzyme digestion, double-enzyme digestion and size selection process and generate comprehensive information of the simulation including enzyme loci number, fragment number, sequences of the fragments, positions of restriction sites on the genome, the coverage of digested fragments on different genome regions and detailed fragment length distribution. Conclusions: RestrictionDigest is an easy-to-use Perl module with flexible parameter settings.With the help of the information produced by the module, researchers can easily determine the most appropriate enzymes to construct the reduced-representation libraries to meet their experimental requirements.

Polyphasic characterisation of Burkholderia cepacia complex species isolated from children with cystic fibrosis

Cystic fibrosis (CF) patients with Burkholderia cepacia complex (Bcc) pulmonary infections have high morbidity and mortality. The aim of this study was to compare different methods for identification of Bcc species isolated from paediatric CF patients. Oropharyngeal swabs from children with CF were used to obtain isolates of Bcc samples to evaluate six different tests for strain identification. Conventional (CPT) and automatised (APT) phenotypic tests, polymerase chain reaction (PCR)-recA, restriction fragment length polymorphism-recA, recA sequencing, and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) were applied. Bacterial isolates were also tested for antimicrobial susceptibility. PCR-recA analysis showed that 36 out of the 54 isolates were Bcc. Kappa index data indicated almost perfect agreement between CPT and APT, CPT and PCR-recA, and APT and PCR-recA to identify Bcc, and MALDI-TOF and recA sequencing to identify Bcc species. The recA sequencing data and the MALDI-TOF data agreed in 97.2% of the isolates. Based on recA sequencing, the most common species identified were Burkholderia cenocepacia IIIA (33.4%), Burkholderia vietnamiensis (30.6%), B. cenocepacia IIIB (27.8%), Burkholderia multivorans (5.5%), and B. cepacia (2.7%). MALDI-TOF proved to be a useful tool for identification of Bcc species obtained from CF patients, although it was not able to identify B. cenocepacia subtypes.