Mutation Reporter Tool: an online tool to interrogate loci of interest, with its utility demonstrated using hepatitis B virus.

BACKGROUND: An online tool, which extracts and summarises nucleotide or amino acid sequence data at specified loci of interest, was developed and tested using the basic core promoter/precore (BCP/PC) region of the hepatitis B virus (HBV). The tool is aimed at researchers without specialist computer skills. METHODS: The tool consists of a web-based front-end, with a CGI script, which runs Python code to generate an output web-page. The Python code searches the input sequence data for a specified anchor motif, after which it generates summary tables and graphs of residue and motif distributions. RESULTS: After the user provides an input file in FASTA format containing aligned sequence data (nucleotides or amino acids) and specifies an anchor motif at a known coordinate, the tool summarizes the nucleotides or amino acids at the specified loci, their frequency and analyzes motif patterns of the loci.The tool can output a graph that displays the frequency of mutations relative to a reference sequence. The tool was used to analyze the BCP/PC region of HBV belonging to subgenotypes A1, A2 and subgenotype D and to serotype HBV. The "Discovery Mode" ignores conserved loci and assists in identifying potential loci of interest. CONCLUSIONS: Although HBV was used to demonstrate the utility of the Mutation Reporter Tool, the tool has wide application as it is genome-agnostic: nucleotide or amino acid sequence data from any organism can be processed. Rapid characterisation of many sequences can be achieved easily when the loci of interest are known. The tool is available online, without charge, at http://hvdr.bioinf.wits.ac.za/tools.

Bioinformatic curation and alignment of genotyped hepatitis B virus (HBV) sequence data from the GenBank public database.

BACKGROUND: Hepatitis B virus (HBV) DNA sequence data from thousands of samples are present in the public sequence databases. No publicly available, up-to-date, multiple sequence alignments, containing full-length and subgenomic fragments per genotype, are available. Such alignments are useful in many analysis applications, including data-mining and phylogenetic analyses. RESULTS: By issuing a query, all HBV sequence data from the GenBank public database was downloaded (67,893 sequences). Full-length and subgenomic sequences, which were genotyped by the submitters (30,852 sequences), were placed into a multiple sequence alignment, for each genotype (genotype A: 5868 sequences, B: 4630, C: 7820, D: 8300, E: 2043, F: 985, G: 189, H: 108, I: 23), according to the results of offline BLAST searches against a custom reference library of full-length sequences. Further curation was performed to improve the alignment. CONCLUSIONS: The algorithm described in this paper generates, for each of the nine HBV genotypes, multiple sequence alignments, which contain full-length and subgenomic fragments. The alignments can be updated as new sequences become available in the online public sequence databases. The alignments are available at http://hvdr.bioinf.wits.ac.za/alignments.

Bioinformatics tools for small genomes, such as hepatitis B virus.

DNA sequence analysis is undertaken in many biological research laboratories. The workflow consists of several steps involving the bioinformatic processing of biological data. We have developed a suite of web-based online bioinformatic tools to assist with processing, analysis and curation of DNA sequence data. Most of these tools are genome-agnostic, with two tools specifically designed for hepatitis B virus sequence data. Tools in the suite are able to process sequence data from Sanger sequencing, ultra-deep amplicon resequencing (pyrosequencing) and chromatograph (trace files), as appropriate. The tools are available online at no cost and are aimed at researchers without specialist technical computer knowledge. The tools can be accessed at http://hvdr.bioinf.wits.ac.za/SmallGenomeTools, and the source code is available online at https://github.com/DrTrevorBell/SmallGenomeTools.

Analysis of ultra-deep pyrosequencing and cloning based sequencing of the basic core promoter/precore/core region of hepatitis B virus using newly developed bioinformatics tools.

AIMS: The aims of this study were to develop bioinformatics tools to explore ultra-deep pyrosequencing (UDPS) data, to test these tools, and to use them to determine the optimum error threshold, and to compare results from UDPS and cloning based sequencing (CBS). METHODS: Four serum samples, infected with either genotype D or E, from HBeAg-positive and HBeAg-negative patients were randomly selected. UDPS and CBS were used to sequence the basic core promoter/precore region of HBV. Two online bioinformatics tools, the "Deep Threshold Tool" and the "Rosetta Tool" (http://hvdr.bioinf.wits.ac.za/tools/), were built to test and analyze the generated data. RESULTS: A total of 10952 reads were generated by UDPS on the 454 GS Junior platform. In the four samples, substitutions, detected at 0.5% threshold or above, were identified at 39 unique positions, 25 of which were non-synonymous mutations. Sample #2 (HBeAg-negative, genotype D) had substitutions in 26 positions, followed by sample #1 (HBeAg-negative, genotype E) in 12 positions, sample #3 (HBeAg-positive, genotype D) in 7 positions and sample #4 (HBeAg-positive, genotype E) in only four positions. The ratio of nucleotide substitutions between isolates from HBeAg-negative and HBeAg-positive patients was 3.5 ∶ 1. Compared to genotype E isolates, genotype D isolates showed greater variation in the X, basic core promoter/precore and core regions. Only 18 of the 39 positions identified by UDPS were detected by CBS, which detected 14 of the 25 non-synonymous mutations detected by UDPS. CONCLUSION: UDPS data should be approached with caution. Appropriate curation of read data is required prior to analysis, in order to clean the data and eliminate artefacts. CBS detected fewer than 50% of the substitutions detected by UDPS. Furthermore it is important that the appropriate consensus (reference) sequence is used in order to identify variants correctly.

Fragment merger: an online tool to merge overlapping long sequence fragments.

While PCR amplicons extend to a few thousand bases, the length of sequences from direct Sanger sequencing is limited to 500-800 nucleotides. Therefore, several fragments may be required to cover an amplicon, a gene or an entire genome. These fragments are typically sequenced in an overlapping fashion and assembled by manually sliding and aligning the sequences visually. This is time-consuming, repetitive and error-prone, and further complicated by circular genomes. An online tool merging two to twelve long overlapping sequence fragments was developed. Either chromatograms or FASTA files are submitted to the tool, which trims poor quality ends of chromatograms according to user-specified parameters. Fragments are assembled into a single sequence by repeatedly calling the EMBOSS merger tool in a consecutive manner. Output includes the number of trimmed nucleotides, details of each merge, and an optional alignment to a reference sequence. The final merge sequence is displayed and can be downloaded in FASTA format. All output files can be downloaded as a ZIP archive. This tool allows for easy and automated assembly of overlapping sequences and is aimed at researchers without specialist computer skills. The tool is genome- and organism-agnostic and has been developed using hepatitis B virus sequence data.

Genotyping and molecular characterization of hepatitis B virus from human immunodeficiency virus-infected individuals in southern Africa.

Hepatitis B virus (HBV) and human immunodeficiency virus (HIV) are hyperendemic in sub-Saharan Africa. The HBV genotypes prevailing in HIV-infected Africans are unknown. Our aim was to determine the HBV genotypes in HIV-infected participants and to identify clinically significant HBV mutations. From 71 HBV DNA(+ve) HIV-infected participants, 49 basic core promoter/precore (BCP/PC) and 29 complete S regions were successfully sequenced. Following phylogenetic analysis of 29 specimens in the complete S region, 28 belonged to subgenotype A1 and one to D3. Mutations affecting HBeAg expression at the transcriptional (1762T1764A), translational (Kozak 1809-1812, initiation 1814-1816, G1896A with C1858T), or post translational levels (G1862T), were responsible for the high HBeAg-negativity observed. The G1862T mutation occurred only in subgenotype A1 isolates, which were found in one third (7/21) of HBsAg(-ve) participants, but in none of the 18 HBsAg(+ve) participants (p<0.05). Pre-S deletion mutants were detected in four HBsAg(+ve) and one HBsAg(-ve) participant/s. The following mutations occurred significantly more frequently in HBV isolated in this study than in strains of the same cluster of the phylogenetic tree: ps1F25L, ps1V88L/A; ps2Q10R, ps2 R48K/T, ps2A53V and sQ129R/H, sQ164A/V/G/D, sV168A and sS174N (p<0.05). ps1I48V/T occurred more frequently in females than males (p<0.05). Isolates with sV168A occurred more frequently in participants with viral loads >200 IU per ml (p<0.05) and only sS174N occurred more frequently in HBsAg(-ve) than in HBsAg(+ve) individuals (p<0.05). Prior to initiation of ART, ten percent, 3 of 29 isolates sequenced, had drug resistance mutations rtV173L, rtL180M+rtM204V and rtV214A, respectively. This study has provided important information on the molecular characteristics of HBV in HIV-infected southern Africans prior to ART initiation, which has important clinical relevance in the management of HBV/HIV co-infection in our unique setting.

Hepatitis B virus infection in human immunodeficiency virus infected southern African adults: occult or overt--that is the question.

Hepatitis B virus (HBV) and human immunodeficiency virus (HIV) share transmission routes and are endemic in sub-Saharan Africa. The objective of the present study was to use the Taormina definition of occult HBV infection, together with stringent amplification conditions, to determine the prevalence and characteristics of HBV infection in antiretroviral treatment (ART)-naïve HIV(+ve) adults in a rural cohort in South Africa. The presence of HBV serological markers was determined by enzyme linked immunoassay (ELISA) tests. HBV DNA-positivity was determined by polymerase chain reaction (PCR) of at least two of three different regions of the HBV genome. HBV viral loads were determined by real-time PCR. Liver fibrosis was determined using the aspartate aminotransferase-to-platelet ratio index. Of the 298 participants, 231 (77.5%) showed at least one HBV marker, with 53.7% HBV DNA(-ve) (resolved) and 23.8% HBV DNA(+ve) (current) [8.7% HBsAg(+ve): 15.1% HBsAg(-ve)]. Only the total number of sexual partners distinguished HBV DNA(+ve) and HBV DNA(-ve) participants, implicating sexual transmission of HBV and/or HIV. It is plausible that sexual transmission of HBV and/or HIV may result in a new HBV infection, superinfection and re-activation as a consequence of immunesuppression. Three HBsAg(-ve) HBV DNA(+ve) participants had HBV viral loads <200 IU/ml and were therefore true occult HBV infections. The majority of HBsAg(-ve) HBV DNA(+ve) participants did not differ from HBsAg(+ve) HBV DNA(+ve) (overt) participants in terms of HBV viral loads, ALT levels or frequency of liver fibrosis. Close to a quarter of HIV(+ve) participants were HBV DNA(+ve), of which the majority were HBsAg(-ve) and were only detected using nucleic acid testing. Detection of HBsAg(-ve) HBV DNA(+ve) subjects is advisable considering they were clinically indistinguishable from HBsAg(+ve) HBV DNA(+ve) individuals and should not be overlooked, especially if lamivudine is included in the ART.