Mutation Reporter Tool: An online tool to interrogate loci of interest, with its utility demonstrated using hepatitis B virus
© Bell and Kramvis; licensee BioMed Central Ltd. 2013
Received: 16 July 2012
Accepted: 28 January 2013
Published: 23 February 2013
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.
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.
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.
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
KeywordsHepatitis B virus Mutations Sequence analysis
Example organism: Hepatitis B virus
Hepatitis B virus (HBV) is one of the most important blood-borne pathogens and is endemic to the sub-Saharan African and southeast Asian regions. Worldwide, around 2 billion people have been exposed to the virus, 240 million are chronically infected, and more than half a million die annually from infection-related liver diseases . At approximately 3,200 nucleotides, the HBV genome is small and has been well-characterized. The genome codes for seven different proteins from four overlapping reading frames (ORFs). To date, nine different genotypes of HBV have been identified: A to D [2, 3], E and F [3–5], G , H , I [8–12] and genotype J has recently been proposed . Subgenotypes have been recogonized in genotypes A to D, F and I, and these are named numerically . Disease progression, clinical manifestation of illness and treatment response differ between these genotypes [15–17].
Mutations (single nucleotide polymorphisms, or “SNPs”) in the genetic sequence of HBV are common, as the virus polymerase lacks proof-reading ability . Patterns of mutations at various known loci have been used to characterize the virus . Certain patterns are characterstic of a particular genotype, or subgenotype [14, 19], and can therefore be used to identify, or “genotype”, a given sample. Patterns at other loci are characteristic of known drug-resistant mutants , or indicate other important characteristics, such as down-regulation of, for example, hepatitis B e antigen (HBeAg) . Therefore, the examination of nucleotides at one or more known loci, either together or individually, is routinely used to characterise HBV sequences. Identification of mutations of interest is not always straighforward, however, for a number of reasons. Firstly, the HBV genome is circular (numbered from position “1” at the EcoR1 restriction site), but sequence data is linear, and position “1” lies within a region of interest, which is typically sequenced both downstream and upstream of this position. Secondly, HBV genotypes are not the same length, ranging from 3182, for genotype D, to 3248 nucleotides for genotype G. Thirdly, insertions and/or deletions of varying length may be present in some isolates, or, fourthly, isolates may be recombinants of two or more known genotypes. Thus, automated analysis of the genome is complex and sequence data should be carefully curated.
Basic core promoter/precore (BCP/PC) mutations
HBeAg is a non-particulate secretory protein expressed by HBV. The pre-core/core open reading frame encodes for HBeAg . The basic core promoter (BCP), which covers the distal X region and the proximal pre-core (PC) region, directs transcription of PC mRNA, which is translated into the pre-core/core fusion protein that is the precursor of HBeAg. This protein has a signal peptide at its amino end that targets it to the endoplasmic reticulum, where it is post-translationally modified by truncation at a fixed site on its amino end and at variable sites on its carboxyl end . Various mutations within the BCP and PC regions affect the expression of HBeAg at the transcriptional, translational and post-translational levels [23, 24]. The BCP A1762T/G1764A mutations affect transcription of the PC mRNA . Mutations that affect HBeAg expression at the translational level include Kozak sequence (1809-1812) mutations and the G1896A stop codon mutation. Substitutions at 1809-1812 are found mainly in subgenotype A1. HBeAg expression is impaired by Kozak mutations by a leaky scanning mechanism . The classical G1896A transition leads to a tryptophan to stop codon mutation, which results in the truncation of HBeAg precursor and abrogation of HBeAg expression . The emergence of G1896A leads to the stabilization of the encapsidation signal (ε) on the pregenomic RNA in genotypes with 1858T, but is rarely found in strains which have 1858C . At the post-translational level, the G1862T mutation, characteteristic of subgenotype A1, introduces a phenylalanine, which interferes with signal peptide cleavage and maturation of HBeAg . Clinically, HBeAg is used as an index of viral replication, infectivity, severity of disease and response to antiviral treatment. Mutations that affect HBeAg expression are clinically relevant  and thus analysis of their distribution is important. We demonstrate the utility of the Mutation Reporter Tool using the BCP/PC mutations as an example.
Loci of interest and patterns of residues
Analysis of loci of interest, which may be dispersed across the genome, and the resulting patterns of these loci, has traditionally been a manual, interactive process, which is time-consuming and error-prone. A new online tool, the Mutation Reporter Tool, has been developed to rapidly and easily display loci of interest and patterns of residues for any sequence data (nucleotides or amino acids) submitted by the user. Feedback from members of the Hepatitis Virus Diversity Research Programme, who used development versions of the tool extensively to analyze HBV sequences, was incorporated into the present version.
Results and Discussion
The Mutation Reporter Tool is one component of a larger project currently in progress, and makes use of a common (shared) Python computer language module, consisting of a “Sequence” class, which contains several methods. The tool consists of a web-based front-end, with which the user interacts, and a CGI script, which runs the Python code and generates the output web-page. The tool has been developed to assist scientists with data analysis and does not require any specialist computer skills or installation. A detailed online tutorial is available. HBV sequence data will be used to demonstrate the utility of the tool.
For example, the basic core promoter/precore (BCP/PC) region of HBV is routinely sequenced. Within this sequence fragment of approximately 500 nucleotides, the highly-conserved motif “AGATTA” is found at co-ordinate 1750. A file containing aligned BCP/PC sequence data is submitted to the tool with “AGATTA” specified as the anchor motif and “1750” as the anchor position. Loci of interest downsteam of 1750 are then specified by their absolute (and known) co-ordinates in the genome. Loci, which are known to affect the expression of HBeAg, are found at 1762, 1764 and 1896. In subgenotype A1, the “Kozak” sequence, which modulates the translation of HBeAg, is located at position “1809-1812”. All these loci are therefore entered into the “Loci” field as “1762,1764,1809-1812,1896”. Only these loci are extracted from each sequence in the input file and included on the output page.
The loci of interest can optionally be grouped into columns according to the “Output grouping” field. The field accepts a comma-separated list of integers, which indicate the number of loci to group into one output column. If no output grouping is specified, the tool will output all loci into one output column. Using the previous example of loci and an output grouping of “2,4,1”, the output would place the nucleotides at 1762 and 1764 together into one column (specified by the output grouping of “2”), the Kozak sequence at 1809-1812 into another column, and nucleotide at 1896 into a third column.
If only some sequences from the input file are to be processed, a “regular expression” can be entered next to the “Include/Exclude” drop-down box. This will then either include (or exclude) sequences for which the FASTA ID matches (or does not match) the regular expression provided. A tutorial describing regular expressions is linked from the input page for reference. Subsets of sequence data stored in one FASTA file can therefore easily be analysed separately, without having to create additional files. FASTA IDs in the output are truncated to the number of characters specified on the input page. If “Output percentages” is not selected, absolute counts are given as output, instead of percentages.
Example Usage: HBV serotypes
In addition to genotypic classification, HBV strains can be classified into one of nine serological subtypes (serotypes) . This classification is determined by the amino acids present at either three or five known positions within the HBV surface antigen (HBsAg) [31–34]. HBV serotype is loosely correlated with genotype . A published decision tree summarizes the interpretation of the amino acid positions to detemine the HBV serotype . Translated (amino acid) sequence data covering the HBV surface gene can be submitted to the Mutation Reporter Tool with the five amino acid positions of interest (122, 160, 127, 159 and 140) specified. An output grouping of “1,1,1,1,1” should be specified to place each amino acid into its own column for easier reading. The amino acids at each position for each sequence can then be examined together with the decision tree to determine the HBV serotype.
Mutation distribution graph
A “Reference motif” can be specified on the input page. This motif should include the reference (“wild-type”) residue for each of the specified loci, in order. For example, if loci “1809-1812,1896” are specified and the input file consists of HBV subgenotype A1 sequences, the reference motif would be “TCATG”. If a reference motif is specified, the output page will include a graph, which indicates the percentage of non-reference (mutant) residues present at each locus. If the input sequence does not contain the ambiguous base “N”, then specifying a reference motif consisting only of “N” characters will result in the tool including all of the residues at each locus, as all residues will not match the reference residue of “N” at each locus. Additional parameters on the input page are used to customize the graph appearance. These include specifying the graph dimensions (in pixels). Loci at which all sequences contain only the reference residue can be suppressed by selecting the appropriate control on the input page. Selecting the “Y-Axis scaled to 100%” control will ensure that the Y-axis of the graph extends from 0% to 100%. This is useful when preparing several graphs which are to be compared with each other. If this control is not selected, the Y-axis will be scaled according to the input data. The raw data used to construct the graph can be downloaded in CSV format from a link on the output page.
Example analysis: Subgenotypes A1 and A2, and genotype D
When the “Discovery Mode” option on the input page is selected, the tool examines the distribution of residues at each of the specified loci and selects for processing only those loci which are not conserved across all input sequences. This mode can be used to “discover” loci of interest by specifying a range of loci, such as “1-100” for example, rather than specific, known loci. The tool will then examine the residues at loci 1 to 100, and will include for further processing and output only those loci at which two or more (different) residues are found. Loci at which only one residue is found will be excluded from the analysis entirely.
When “Discovery Mode” is selected, the “Output grouping”, “Reference motif” and graphing parameters are disabled, as the number and position of loci, which would be included in the final analysis is only known after the tool has processed the file. Also, as this is a “discovery mode”, it will not be known in advance which loci should logically be grouped together as a unit of interest.
A limitation, by design, is that the sequence data must be aligned. The position of the anchor motif in the first sequence is taken as the anchor position, and loci in all sequences are referenced according to this position. If the input sequence data are not aligned, or if the anchor motif is incorrect, the tool may return incorrect data. Whilst the number of loci, which can be specified, is not limited, it may not be feasible to enter more than a few dozen loci, as this generates a large amount of output data. All loci values must be greater than the anchor position. Updates to the tool will be made to address limitations as necessary.
As an online tool, available free of charge, no download or installation is required. As demonstrated, this tool can be used for both genotyping and serotyping of HBV without the requirement of computer skills or knowledge of phylogenetics. However, as the tool is genome-agnostic, it has a wide application and nucleotide or amino acid sequence data from any organism can be analysed. Loci of interest, which may be located many hundreds of residues apart, can easily be extracted and their distribution summarised. Rapid characterisation of many sequences, or subsets of sequences, can be achieved easily when the loci of interest are known. Using the “Discovery Mode”, conserved and therefore uninformative loci, are automatically ignored, and potential loci of interest can be found and identified.
The Mutation Reporter Tool consists of a web-based front-end (“client” interface) with which the user interacts, and a CGI (common gateway interface) script on a server, which runs Python language  code to generate the output web-page. The tool is one component of a larger project currently under development, which makes use of a common, shared Python library. The input FASTA file which the user specifies is saved locally (on the server) by the CGI script and then processed by the Python library. Methods within this library are responsible for loading sequence data from a FASTA file, processing the input parameters, and extracting the requested data from the FASTA file. The output HTML page is written to disk by the Python script. The optional output graphs are generated using the ggplot2 graphics library  in the R statistical programming language . If graphs are requested by the user, the Python script writes the relevant data to disk as a CSV (comma-separated value) file. A short R script, which is customized based on the input parameters specified, is also written to disk. The Python script then calls the R script, which generates the graph and writes it to disk. The images are then linked on the output HTML page. The tool is an online resource, which requires a client browser to connect to the tool’s web-server. As such, there is no stand-alone, offline version available for download.
The tool, which assumes that the submitted sequence data is aligned, finds the first occurrence of the anchor motif in the first sequence in the input file. The first character of the anchor motif is then considered to be at the position specified as the anchor position. Sequence data at each of the specified loci for all sequences in the file is then accessed and tabulated. Loci positions are mapped to positions in the sequence data using the anchor motif as an offset value. Data from the loci specified are grouped into columns according to the “output grouping” field. If this field is not specified, all loci are grouped into one output column. If a sequence ID pattern was specified, the tool executes the appropriate regular expression match on the FASTA IDs in the input file. In “Discovery Mode”, loci at which no variation is found are excluded.
TB received bursaries from the National Research Foundation (NRF), the Medical Research Foundation (MRC), the Poliomyelitis Research Foundation, National Bioinformatics Network and the University of the Witwatersrand. AK received funding from the NRF (GUN#65530) and the MRC. Mark Keyter and Mukhlid Yousif tested the tool extensively and provided valuable feedback and suggestions.
- WHO: World Health Organization: Hepatitis B fact sheet 204 (July 2012, Revision). 2012.http://www.who.int/mediacentre/factsheets/fs204/en/
- Okamoto H, Tsuda F, Sakugawa H, Sastrosoewignjo RI, Imai M, Miyakawa Y, Mayumi M: Typing hepatitis B virus by homology in nucleotide sequence: comparison of surface antigen subtypes. J Gen Virol 1988, 69:2575–2583. 10.1099/0022-1317-69-10-2575PubMedView Article
- Norder H, Couroucé AM, Magnius LO: Complete genomes, phylogenetic relatedness, and structural proteins of six strains of the hepatitis B virus, four of which represent two new genotypes. Virology 1994, 198:489–503. 10.1006/viro.1994.1060PubMedView Article
- Norder H, Hammas B, Lofdahl S, Courouce AM, Magnius LO: Comparison of the amino acid sequences of nine different serotypes of hepatitis B surface antigen and genomic classification of the corresponding hepatitis B virus strains. J Gen Virol 1992,73(Pt 5):1201–1208.PubMedView Article
- Naumann H, Schaefer S, Yoshida CF, Gaspar AM, Repp R, Gerlich WH: Identification of a new hepatitis B virus (HBV) genotype from Brazil that expresses HBV surface antigen subtype adw4. J Gen Virol 1993, 74:1627–1632. 10.1099/0022-1317-74-8-1627PubMedView Article
- Stuyver L, Gendt SD, Geyt CV, Zoulim F, Fried M, Schinazi RF, Rossau R: A new genotype of hepatitis B virus: complete genome and phylogenetic relatedness. J Gen Virol 2000, 81:67–74.PubMed
- Arauz-Ruiz P, Norder H, Robertson BH, Magnius LO: Genotype H: a new Amerindian genotype of hepatitis B virus revealed in Central America. J Gen Virol 2002, 83:2059–2073.PubMed
- Olinger CM, Jutavijittum P, Hübschen JM, Yousukh A, Samountry B, Thammavong T, Toriyama K, Muller CP: Possible new hepatitis B virus genotype, southeast Asia. Emerg Infect Dis 2008, 14:1777–1780. 10.3201/eid1411.080437PubMedView ArticlePubMed Central
- Tran TTH, Trinh TN, Abe K: New complex recombinant genotype of hepatitis B virus identified in Vietnam. J Virol 2008, 82:5657–5663. 10.1128/JVI.02556-07PubMedView Article
- Arankalle VA, Gandhe SS, Borkakoty BJ, Walimbe AM, Biswas D, Mahanta J: A novel HBV recombinant (genotype I) similar to Vietnam/Laos in a primitive tribe in eastern India. J Viral Hepatitis 2010, 17:501–510.
- Osiowy C, Kaita K, Solar K, Mendoza K: Molecular characterization of hepatitis B virus and a 9-year clinical profile in a patient infected with genotype I. J Med Virol 2010, 82:942–948. 10.1002/jmv.21758PubMedView Article
- Yu H, Yuan Q, Ge SX, Wang HY, Zhang YL, Chen QR, Zhang J, Chen PJ, Xia NS: Molecular and phylogenetic analyses suggest an additional Hepatitis B virus genotype “I”. PLOS One 2010, 5:e9297. 10.1371/journal.pone.0009297PubMedView ArticlePubMed Central
- Tatematsu K, Tanaka Y, Kurbanov F, Sugauchi F, Mano S, Maeshiro T, Nakayoshi T, Wakuta M, Miyakawa Y, Mizokami M: A genetic variant of hepatitis B virus divergent from known human and ape genotypes isolated from a Japanese patient and provisionally assigned to new genotype J. J Virol 2009, 83:10538–10547. 10.1128/JVI.00462-09PubMedView ArticlePubMed Central
- Kramvis A, Arakawa K, Yu MC, Nogueira R, Stram DO, Kew MC: Relationship of serological subtype, basic core promoter and precore mutations to genotypes/subgenotypes of hepatitis B virus. J Med Virol 2008, 80:27–46. 10.1002/jmv.21049PubMedView Article
- Mayerat C, Mantegani A, Frei PC: Does hepatitis B virus (HBV) genotype influence the clinical outcome of HBV infection? J Viral Hepatitis 1999, 6:299–304. 10.1046/j.1365-2893.1999.00174.xView Article
- Sumi H, Yokosuka O, Seki N, Arai M, Imazeki F, Kurihara T, Kanda T, Fukai K, Kato M, Saisho H: Influence of hepatitis B virus genotypes on the progression of chronic type B liver disease. Hepatology 2003, 37:19–26. 10.1053/jhep.2003.50036PubMedView Article
- Kramvis A, Kew MC: Relationship of genotypes of hepatitis B virus to mutations, disease progression and response to antiviral therapy. J Viral Hepatitis 2005, 12:456–464. 10.1111/j.1365-2893.2005.00624.xView Article
- Steinhauer DA, Holland JJ: Direct method for quantitation of extreme polymerase error frequencies at selected single base sites in viral RNA. J Virol 1986, 57:219–228.PubMedPubMed Central
- Kramvis A, Kew M, François G: Hepatitis B virus genotypes. Vaccine 2005, 23:2409–2423. 10.1016/j.vaccine.2004.10.045PubMedView Article
- Zoulim F, Locarnini S: Hepatitis B virus resistance to Nucleos(t)ide analogues. Gastroenterology 2009, 137:1593–1608. 10.1053/j.gastro.2009.08.063PubMedView Article
- Revill P, Yuen L, Walsh R, Perrault M, Locarnini S, Kramvis A: Bioinformatic analysis of the hepadnavirus e-antigen and its precursor identifies remarkable sequence conservation in all orthohepadnaviruses. J Med Virol 2010, 82:104–115. 10.1002/jmv.21645PubMedView Article
- Ou JH, Laub O, Rutter WJ: Hepatitis B virus gene function: the precore region targets the core antigen to cellular membranes and causes the secretion of the e antigen. Proc Nat Acad Sci USA 1986, 83:1578–1582. 10.1073/pnas.83.6.1578PubMedView ArticlePubMed Central
- Kramvis A, Kew MC: Epidemiology of hepatitis B virus in Africa, its genotypes and clinical associations of genotypes. Hepatology Res 2007, 37:S9-S19. 10.1111/j.1872-034X.2007.00098.xView Article
- Tong S: Impact of viral genotypes and naturally occurring mutations on biological properties of hepatitis B virus. Hepatology Res 2007, 37:S3-S8. 10.1111/j.1872-034X.2007.00097.xView Article
- Buckwold VE, Xu Z, Chen M, Yen TS, Ou JH: Effects of a naturally occurring mutation in the hepatitis B virus basal core promoter on precore gene expression and viral replication. J Virol 1996, 70:5845–5851.PubMedPubMed Central
- Ahn SH, Kramvis A, Kawai S, Spangenberg HC, Li J, Kimbi G, Kew M, Wands J, Tong S: Sequence variation upstream of precore translation initiation codon reduces hepatitis B virus e antigen production. Gastroenterology 2003, 125:1370–1378. 10.1016/j.gastro.2003.07.016PubMedView Article
- Carman WF, Hadziyannis S, Mcgarvey MJ, Jacyna MR, Karayiannis P, Makris A, Thomas HC: Mutation preventing formation of hepatitis B e antigen in patients with chronic hepatitis B infection. Lancet 1989, 334:588–591. 10.1016/S0140-6736(89)90713-7View Article
- Lok AS, Akarca U, Greene S: Mutations in the pre-core region of hepatitis B virus serve to enhance the stability of the secondary structure of the pre-genome encapsidation signal. Proc Nat Acad Sci USA 1994, 91:4077–4081. 10.1073/pnas.91.9.4077PubMedView ArticlePubMed Central
- Chen CY, Crowther C, Kew MC, Kramvis A: A valine to phenylalanine mutation in the precore region of hepatitis B virus causes intracellular retention and impaired secretion of HBe-antigen. Hepatology Res 2008, 38:580–592. 10.1111/j.1872-034X.2007.00315.xView Article
- Magnius LO, Norder H: Genotypes and molecular epidemiology of the Hepatitis B virus as reflected by sequence variability of the S-Gene. Intervirology 1995, 38:24–34.PubMed
- Wands JR, Wong MA, Shorey J, Brown RD, Marciniak RA, Isselbacher KJ: Hepatitis B viral antigenic structure: signature analysis by monoclonal radioimmunoassays. Proc Nat Acad Sci USA 1984, 81:2237–2241. 10.1073/pnas.81.7.2237PubMedView ArticlePubMed Central
- Mimms LT, Floreani M, Tyner J, Whitters E, Rosenlof R, Wray L, Goetze A, Sarin V, Eble K: Discrimination of hepatitis B virus (HBV) subtypes using monoclonal antibodies to the PreS1 and PreS2 domains of the viral envelope. Virology 1990, 176:604–619. 10.1016/0042-6822(90)90031-LPubMedView Article
- Swenson PD, Riess JT, Krueger LE: Determination of HBsAg subtypes in different high risk populations using monoclonal antibodies. J Virological Methods 1991, 33:27–38. 10.1016/0166-0934(91)90004-JView Article
- Purdy MA, Talekar G, Swenson P, Araujo A, Fields A: A new algorithm for deduction of Hepatitis B surface antigen subtype determinants from the amino acid sequence. Intervirology 2007, 50:45–51. 10.1159/000096312PubMedView Article
- Tanaka Y, Hasegawa I, Kato T, Orito E, Hirashima N, Acharya SK, Gish RG, Kramvis A, Kew MC, Yoshihara N, Shrestha SM, Khan M, Miyakawa Y, Mizokami M: A case-control study for differences among hepatitis B virus infections of genotypes A (subtypes Aa and Ae) and D. Hepatology 2004, 40:747–755. 10.1002/hep.20365PubMedView Article
- van Rossum G: Python. [http://www.python.org].
- Wickham H: ggplot2: Elegant Graphics for Data Analysis. New York: Springer; 2009.View Article
- R Core Team: R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2012. [http://www.R-project.org/]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.