Optimal pooling for genome re-sequencing with ultra-high-throughput short-read technologies

New generation sequencing technologies offer unique opportunities and challenges for re-sequencing studies. In this article, we focus on re-sequencing experiments using the Solexa technology, based on bacterial artificial chromosome (BAC) clones, and address an experimental design problem. In these specific experiments, approximate coordinates of the BACs on a reference genome are known, and fine-scale differences between the BAC sequences and the reference are of interest. The high-throughput characteristics of the sequencing technology makes it possible to multiplex BAC sequencing experiments by pooling BACs for a cost-effective operation. However, the way BACs are pooled in such re-sequencing experiments has an effect on the downstream analysis of the generated data, mostly due to subsequences common to multiple BACs. The experimental design strategy we develop in this article offers combinatorial solutions based on approximation algorithms for the well-known max n-cut problem and the related max n-section problem on hypergraphs. Our algorithms, when applied to a number of sample cases give more than a 2-fold performance improvement over random partitioning

Intra- and interspecific polymorphisms ofLeishmania donovani andL. tropica minicircle DNA

A pair of degenerate polymerase chain reaction (PCR) primers (LEI-1, TCG GAT CC[C,T] [G,C]TG GGT AGG GGC GT; LEI-2, ACG GAT CC[G,C] [G,C][A,C]C TAT [A,T]TT ACA CC) defining a 0.15-kb segment ofLeishmania minicircle DNA was constructed. These primers amplified not only inter- but also intraspecifically polymorphic sequences. Individual sequences revealed a higher intraspecific than interspecific divergence. It is concluded that individual sequences are of limited relevance for species determination. In contrast, when a data base of 19 different sequences was analyzed in a dendrographic plot, an accurate species differentiation was feasible

Tissue-specific expression of high-voltage-activated dihydropyridine-sensitive L-type calcium channels

The cloning of the cDNA for the α1 subunit of L-type calcium channels revealed that at least two genes (CaCh1 and CaCh2) exist which give rise to several splice variants. The expression of mRNA for these α1 subunits and the skeletal muscle α2/δ, β and γ subunits was studied in rabbit tissues and BC3H1 cells. Nucleic-acid-hybridization studies showed that the mRNA of all subunits are expressed in skeletal muscle, brain, heart and aorta. However, the α1-, β- and γ-specific transcripts had different sizes in these tissues. Smooth muscle and heart contain different splice variants of the CaCh2 gene. The α1, β and γ mRNA are expressed together in differentiated but not in proliferating BC3H1 cells. A probe specific for the skeletal muscle α2/δ subunit did not hybridize to poly(A)-rich RNA from BC3H1 cells. These results suggest that different splice variants of the genes for the α1, β and γ subunits exist in tissues containing L-type calcium channels, and that their expression is regulated in a coordinate manner

Clustering exact matches of pairwise sequence alignments by weighted linear regression

<p>Abstract</p> <p>Background</p> <p>At intermediate stages of genome assembly projects, when a number of contigs have been generated and their validity needs to be verified, it is desirable to align these contigs to a reference genome when it is available. The interest is not to analyze a detailed alignment between a contig and the reference genome at the base level, but rather to have a rough estimate of where the contig aligns to the reference genome, specifically, by identifying the starting and ending positions of such a region. This information is very useful in ordering the contigs, facilitating post-assembly analysis such as gap closure and resolving repeats. There exist programs, such as BLAST and MUMmer, that can quickly align and identify high similarity segments between two sequences, which, when seen in a dot plot, tend to agglomerate along a diagonal but can also be disrupted by gaps or shifted away from the main diagonal due to mismatches between the contig and the reference. It is a tedious and practically impossible task to visually inspect the dot plot to identify the regions covered by a large number of contigs from sequence assembly projects. A forced global alignment between a contig and the reference is not only time consuming but often meaningless.</p> <p>Results</p> <p>We have developed an algorithm that uses the coordinates of all the exact matches or high similarity local alignments, clusters them with respect to the main diagonal in the dot plot using a weighted linear regression technique, and identifies the starting and ending coordinates of the region of interest.</p> <p>Conclusion</p> <p>This algorithm complements existing pairwise sequence alignment packages by replacing the time-consuming seed extension phase with a weighted linear regression for the alignment seeds. It was experimentally shown that the gain in execution time can be outstanding without compromising the accuracy. This method should be of great utility to sequence assembly and genome comparison projects.</p

Repositioning the Catalytic Triad Aspartic Acid of Haloalkane Dehalogenase: Effects on Stability, Kinetics, and Structure

Haloalkane dehalogenase (DhlA) catalyzes the hydrolysis of haloalkanes via an alkyl-enzyme intermediate. The covalent intermediate, which is formed by nucleophilic substitution with Asp124, is hydrolyzed by a water molecule that is activated by His289. The role of Asp260, which is the third member of the catalytic triad, was studied by site-directed mutagenesis. Mutation of Asp260 to asparagine resulted in a catalytically inactive D260N mutant, which demonstrates that the triad acid Asp260 is essential for dehalogenase activity. Furthermore, Asp260 has an important structural role, since the D260N enzyme accumulated mainly in inclusion bodies during expression, and neither substrate nor product could bind in the active-site cavity. Activity for brominated substrates was restored to D260N by replacing Asn148 with an aspartic or glutamic acid. Both double mutants D260N+N148D and D260N+N148E had a 10-fold reduced kcat and 40-fold higher Km values for 1,2-dibromoethane compared to the wild-type enzyme. Pre-steady-state kinetic analysis of the D260N+N148E double mutant showed that the decrease in kcat was mainly caused by a 220-fold reduction of the rate of carbon-bromine bond cleavage and a 10-fold decrease in the rate of hydrolysis of the alkyl-enzyme intermediate. On the other hand, bromide was released 12-fold faster and via a different pathway than in the wild-type enzyme. Molecular modeling of the mutant showed that Glu148 indeed could take over the interaction with His289 and that there was a change in charge distribution in the tunnel region that connects the active site with the solvent. On the basis of primary structure similarity between DhlA and other α/β-hydrolase fold dehalogenases, we propose that a conserved acidic residue at the equivalent position of Asn148 in DhlA is the third catalytic triad residue in the latter enzymes.

Anti-nausea effects and pharmacokinetics of ondansetron, maropitant and metoclopramide in a low-dose cisplatin model of nausea and vomiting in the dog: a blinded crossover study

Nausea is a subjective sensation which is difficult to measure in non-verbal species. The aims of this study were to determine the efficacy of three classes of antiemetic drugs in a novel low dose cisplatin model of nausea and vomiting and measure change in potential nausea biomarkers arginine vasopressin (AVP) and cortisol. A four period cross-over blinded study was conducted in eight healthy beagle dogs of both genders. Dogs were administered 18 mg/m2 cisplatin intravenously, followed 45 min later by a 15 min infusion of either placebo (saline) or antiemetic treatment with ondansetron (0.5 mg/kg; 5-HT3 antagonist), maropitant (1 mg/kg; NK1 antagonist) or metoclopramide (0.5 mg/kg; D2 antagonist). The number of vomits and nausea associated behaviours, scored on a visual analogue scale, were recorded every 15 min for 8 h following cisplatin administration. Plasma samples were collected to measure AVP, cortisol and antiemetic drug concentrations

The nucleotide and partial amino acid sequences of rat fetuin

Fetuins are among the major plasma proteins, yet their biological role has remained elusive. Here we report the molecular cloning of rat fetuin and the sequence analysis of a full-length clone, RF619 of 1456 bp with an open reading frame of 1056 bp encoding 352 amino acid residues. The coding part of RF619 was identical with the cDNA sequence of the natural inhibitor of the insulin receptor tyrosine kinase from rat (pp63) except for four substitutions and a single base insertion causing divergence of the predicted protein sequences. Partial amino acid sequences of rat plasma fetuin were in agreement with the predictions based on the RF619 cDNA. Purified rat fetuin inhibited the insulin receptor tyrosine kinase in vitro. Therefore, we conclude that RF619 and pp63 cDNA encode the same protein, i.e. authentic rat fetuin which is a functional tyrosine kinase inhibitor

A direct comparison of the KB™ Basecaller and phred for identifying the bases from DNA sequencing using chain termination chemistry

<p>Abstract</p> <p>Background</p> <p>Relatively recently, the software KB™ Basecaller has replaced <it>phred </it>for identifying the bases from raw sequence data in DNA sequencing employing dideoxy chemistry. We have measured quantitatively the consequences of that change.</p> <p>Results</p> <p>The high quality sequence segment of reads derived from the KB™ Basecaller were, on average, 30-to-50 bases longer than reads derived from <it>phred</it>. However, microbe identification appeared to have been unaffected by the change in software.</p> <p>Conclusions</p> <p>We have demonstrated a modest, but statistically significant, superiority in high quality read length of the KB™ Basecaller compared to <it>phred</it>. We found no statistically significant difference between the numbers of microbial species identified from the sequence data.</p