BackgroundMany high throughput sequencing (HTS) approaches, such as the Roche/454 platform, produce sequences in which the quality of the sequence (as measured by a Phred-like quality scores) decreases linearly across a sequence read. Undertaking quality trimming of this data is essential to enable confidence in the results of subsequent downstream analysis. Here, we have developed a novel, highly sensitive and accurate approach (QTrim) for the quality trimming of sequence reads generated using the Roche/454 sequencing platform (or any platform with long reads that outputs Phred-like quality scores).ResultsThe performance of QTrim was evaluated against all other available quality trimming approaches on both poor and high quality 454 sequence data. In all cases, QTrim appears to perform equally as well as the best other approach (PRINSEQ) with these two methods significantly outperforming all other methods. Further analysis of the trimmed data revealed that the novel trimming approach implemented in QTrim ensures that the prevalence of low quality bases in the resulting trimmed data is substantially lower than PRINSEQ or any of the other approaches tested.ConclusionsQTrim is a novel, highly sensitive and accurate algorithm for the quality trimming of Roche/454 sequence reads. It is implemented both as an executable program that can be integrated with standalone sequence analysis pipelines and as a web-based application to enable individuals with little or no bioinformatics experience to quality trim their sequence data.

In this preliminary study we show that in 2008, 3 years after antiretroviral therapy was introduced into the Karonga District, Malawi, a greater than expected number of drug-naive individuals have been infected with HIV-1 subtype C virus harboring major and minor drug resistance mutations (DRMs). From a sample size of 40 reverse transcriptase (RT) consensus sequences from drug-naive individuals we found five showing NRTI and four showing NNRTI mutations with one individual showing both. From 29 protease consensus sequences, again from drug-naive individuals, we found evidence of minor DRMs in three. Additional major and minor DRMs were found in clonal sequences from a number of individuals that were not present in the original consensus sequences. This clearly illustrates the importance of sequencing multiple HIV-1 variants from individuals to fully assess drug resistance.

Erythropoietic protoporphyria (EPP) is caused by deficiency of ferrochelatase (FECH), which incorporates iron into protoporphyrin IX (PPIX) to form heme. Excitation of accumulated PPIX by light generates oxygen radicals that evoke excessive pain and, after longer light exposure, cause ulcerations in exposed skin areas of individuals with EPP. Moreover, ∼5% of the patients develop a liver dysfunction as a result of PPIX accumulation. Most patients (∼97%) have a severe FECH mutation (Mut) in trans to an intronic polymorphism (c.315-48C), which reduces ferrochelatase synthesis by stimulating the use of an aberrant 3′ splice site 63 nt upstream of the normal site for exon 4. In contrast, with the predominant c.315-48T allele, the correct splice site is mostly used, and individuals with a T/Mut genotype do not develop EPP symptoms. Thus, the C allele is a potential target for therapeutic approaches that modify this splicing decision. To provide a model for pre-clinical studies of such approaches, we engineered a mouse containing a partly humanized Fech gene with the c.315-48C polymorphism. F1 hybrids obtained by crossing these mice with another inbred line carrying a severe Fech mutation (named m1Pas) show a very strong EPP phenotype that includes elevated PPIX in the blood, enlargement of liver and spleen, anemia, as well as strong pain reactions and skin lesions after a short period of light exposure. In addition to the expected use of the aberrant splice site, the mice also show a strong skipping of the partly humanized exon 3. This will limit the use of this model for certain applications and illustrates that engineering of a hybrid gene may have unforeseeable consequences on its splicing.