Base-calling algorithm with vocabulary (BCV) method for analyzing population sequencing chromatograms.

Sanger sequencing is a common method of reading DNA sequences. It is less expensive than high-throughput methods, and it is appropriate for numerous applications including molecular diagnostics. However, sequencing mixtures of similar DNA of pathogens with this method is challenging. This is important because most clinical samples contain such mixtures, rather than pure single strains. The traditional solution is to sequence selected clones of PCR products, a complicated, time-consuming, and expensive procedure. Here, we propose the base-calling with vocabulary (BCV) method that computationally deciphers Sanger chromatograms obtained from mixed DNA samples. The inputs to the BCV algorithm are a chromatogram and a dictionary of sequences that are similar to those we expect to obtain. We apply the base-calling function on a test dataset of chromatograms without ambiguous positions, as well as one with 3-14% sequence degeneracy. Furthermore, we use BCV to assemble a consensus sequence for an HIV genome fragment in a sample containing a mixture of viral DNA variants and to determine the positions of the indels. Finally, we detect drug-resistant Mycobacterium tuberculosis strains carrying frameshift mutations mixed with wild-type bacteria in the pncA gene, and roughly characterize bacterial communities in clinical samples by direct 16S rRNA sequencing.

A novel ATP-binding cassette transporter is responsible for resistance to viologen herbicides in the cyanobacterium Synechocystis sp. PCC 6803.

The charged quaternary ammonium compounds--methyl, ethyl and benzyl viologens--generate reactive oxygen species in photosynthetic cells. Three independent methyl viologen-resistant spontaneous mutants of Synechocystis sp. PCC 6803 were identified, in which the conserved R115 residue of the Slr1174 protein was replaced with G115, L115 and C115. The Slr1174 protein of the DUF990 family is related to the permease subunit of an ABC-2-type transporter and its R115 mutation was found to be solely responsible for the observed methyl viologen resistance. Bioinformatic analysis showed that in various bacterial genomes, two genes encoding another permease subunit and the ATPase component of an ATP-binding cassette transporter form putative operons with slr1174 orthologs, suggesting that the protein products of these genes may form functional transporters. The corresponding genes in Synechocystis sp. PCC 6803, i.e. slr0610 for the permease and slr1901 for the ATPase, did not form such an operon. However, insertional inactivation of any slr1174, slr0610 or slr1901 genes in both the wild-type and the R115-resistant mutant resulted in increased sensitivity to methyl, ethyl and benzyl viologens; moreover, single- and double-insertion mutants did not differ in their viologen sensitivity. Our data suggest that Slr1901, Slr1174 and Slr0610 form a heteromeric ATP-binding cassette-type viologen exporter, in which each component is critical for viologen extrusion. Because the greatest increase in mutant sensitivity was observed in the case of ethyl viologen, the three proteins have been named EvrA (Slr1901), EvrB (Slr1174) and EvrC (Slr0610). This is the first report of a function for a DUF990 family protein.