Lower bounds on multiple sequence alignment using exact 3-way alignment

<p>Abstract</p> <p>Background</p> <p>Multiple sequence alignment is fundamental. Exponential growth in computation time appears to be inevitable when an optimal alignment is required for many sequences. Exact costs of optimum alignments are therefore rarely computed. Consequently much effort has been invested in algorithms for alignment that are heuristic, or explore a restricted class of solutions. These give an upper bound on the alignment cost, but it is equally important to determine the quality of the solution obtained. In the absence of an optimal alignment with which to compare, lower bounds may be calculated to assess the quality of the alignment. As more effort is invested in improving upper bounds (alignment algorithms), it is therefore important to improve lower bounds as well. Although numerous cost metrics can be used to determine the quality of an alignment, many are based on sum-of-pairs (SP) measures and their generalizations.</p> <p>Results</p> <p>Two standard and two new methods are considered for using exact 2-way and 3-way alignments to compute lower bounds on total SP alignment cost; one new method fares well with respect to accuracy, while the other reduces the computation time. The first employs exhaustive computation of exact 3-way alignments, while the second employs an efficient heuristic to compute a much smaller number of exact 3-way alignments. Calculating all 3-way alignments exactly and computing their average improves lower bounds on sum of SP cost in <it>v</it>-way alignments. However judicious selection of a subset of all 3-way alignments can yield a further improvement with minimal additional effort. On the other hand, a simple heuristic to select a random subset of 3-way alignments (a random packing) yields accuracy comparable to averaging all 3-way alignments with substantially less computational effort.</p> <p>Conclusion</p> <p>Calculation of lower bounds on SP cost (and thus the quality of an alignment) can be improved by employing a mixture of 3-way and 2-way alignments.</p

Evaluation of a novel real-time PCR test based on the ssrA gene for the identification of group B streptococci in vaginal swabs

<p>Abstract</p> <p>Background</p> <p>Despite the implementation of prevention guidelines, early-onset group B streptococci (GBS) disease remains a cause of neonatal morbidity and mortality worldwide. Strategies to identify women who are at risk of transmitting GBS to their infant and the administration of intrapartum antibiotics have greatly reduced the incidence of neonatal GBS disease. However, there is a requirement for a rapid diagnostic test for GBS that can be carried out in a labour ward setting especially for women whose GBS colonisation status is unknown at the time of delivery. We report the design and evaluation of a real-time PCR test (<it>RiboSEQ </it>GBS test) for the identification of GBS in vaginal swabs from pregnant women.</p> <p>Methods</p> <p>The qualitative real-time PCR <it>RiboSEQ </it>GBS test was designed based on the bacterial <it>ssrA </it>gene and incorporates a competitive internal standard control. The analytical sensitivity of the test was established using crude lysate extracted from serial dilutions of overnight GBS culture using the IDI Lysis kit. Specificity studies were performed using DNA prepared from a panel of GBS strains, related streptococci and other species found in the genital tract environment. The <it>RiboSEQ </it>GBS test was evaluated on 159 vaginal swabs from pregnant women and compared with the GeneOhm™ StrepB Assay and culture for the identification of GBS.</p> <p>Results</p> <p>The <it>RiboSEQ </it>GBS test is specific and has an analytical sensitivity of 1-10 cell equivalents. The <it>RiboSEQ </it>GBS test was 96.4% sensitive and 95.8% specific compared to "gold standard" culture for the identification of GBS in vaginal swabs from pregnant women. In this study, the <it>RiboSEQ </it>GBS test performed slightly better than the commercial BD GeneOhm™ StrepB Assay which gave a sensitivity of 94.6% and a specificity of 89.6% compared to culture.</p> <p>Conclusion</p> <p>The <it>RiboSEQ </it>GBS test is a valuable method for the rapid, sensitive and specific detection of GBS in pregnant women. This study also validates the <it>ssrA </it>gene as a suitable and versatile target for nucleic acid-based diagnostic tests for bacterial pathogens.</p

Spanning forests and the q-state Potts model in the limit q \to 0

We study the q-state Potts model with nearest-neighbor coupling v=e^{\beta J}-1 in the limit q,v \to 0 with the ratio w = v/q held fixed. Combinatorially, this limit gives rise to the generating polynomial of spanning forests; physically, it provides information about the Potts-model phase diagram in the neighborhood of (q,v) = (0,0). We have studied this model on the square and triangular lattices, using a transfer-matrix approach at both real and complex values of w. For both lattices, we have computed the symbolic transfer matrices for cylindrical strips of widths 2 \le L \le 10, as well as the limiting curves of partition-function zeros in the complex w-plane. For real w, we find two distinct phases separated by a transition point w=w_0, where w_0 = -1/4 (resp. w_0 = -0.1753 \pm 0.0002) for the square (resp. triangular) lattice. For w > w_0 we find a non-critical disordered phase, while for w < w_0 our results are compatible with a massless Berker-Kadanoff phase with conformal charge c = -2 and leading thermal scaling dimension x_{T,1} = 2 (marginal operator). At w = w_0 we find a "first-order critical point": the first derivative of the free energy is discontinuous at w_0, while the correlation length diverges as w \downarrow w_0 (and is infinite at w = w_0). The critical behavior at w = w_0 seems to be the same for both lattices and it differs from that of the Berker-Kadanoff phase: our results suggest that the conformal charge is c = -1, the leading thermal scaling dimension is x_{T,1} = 0, and the critical exponents are \nu = 1/d = 1/2 and \alpha = 1.Comment: 131 pages (LaTeX2e). Includes tex file, three sty files, and 65 Postscript figures. Also included are Mathematica files forests_sq_2-9P.m and forests_tri_2-9P.m. Final journal versio