<title>Background</title> <p>A large number of genomes of different strains and closely related species of pathogens have been sequenced and many others are in the process. A detailed analysis of these genomic sequences can help us to decipher and establish genotype to phenotype relationship. The organisms evolve through a series of molecular changes reflected in genomic sequences and some of these are evolutionarily selected based on survival in a specific ecological niche. Characterization of sequence alterations in closely related organisms can help us to understand genome evolution at the molecular level in short time span, for example emergence of new endemic strains in a few decades. Therefore it is important to catalog all the sequence differences between any two organisms so that these can be a basis for designing experiments linking phenotype to genotype. In pathogenic organisms such a database can be useful in identifying species and strain-specific markers that can be a basis for designing diagnostic reagents.</p> <p>A number of molecular mechanisms have been described that are responsible for genomic changes [##REF##12433587##1##]. These contribute to single nucleotide polymorphisms (SNP), variable number of tandem repeats, insertion/deletion with or without involving transposable elements and recombination. Many of these have been used as markers for identification of strains and diagnosis of pathogens [##REF##7690367##2##, ####UREF##0##3##, ##REF##9157152##4####9157152##4##]. <italic>M. tuberculosis </italic>is a major cause of morbidity and mortality throughout the world. Genomic variations in this organism have been used to type pathogenic strains in a limited scale [##UREF##1##5##,##REF##9634230##6##]. A comprehensive database of all the genomic variations of <italic>M. tuberculosis </italic>is not currently available though some attempts have been made in this direction. For example, MTBreg (please see Availability &amp; requirements for more information) covers variations that are detected using spoligotyping, MycoDB (please see Availability &amp; requirements for more information) [##UREF##2##7##], MycoperonDB (please see Availability &amp; requirements for more information) [##REF##17254314##8##] and GenoMycDB (please see Availability &amp; requirements for more information) [##REF##16755503##9##] have some features that allow comparison between two genomes in a limited manner. In this report we describe a comprehensive database of genomic differences among strains and species of <italic>Mycobacteria </italic>belonging to the <italic>M. tuberculosis </italic>complex. The variations have been identified using ABWGC, a comparative genomic tool previously described by us [##UREF##3##10##]. We hope that this database will be useful to clinicians and basic scientists interested in understanding <italic>Mycobacterial </italic>diseases.</p>
<title>Utility and discussion</title> <p>MGDD has a web interface for the retrieval of genomic diversity information. A search can be initiated by first selecting genomes from the "Query" and "Subject" scroll down menu-bar. Currently there is information about six organisms and these can be selected in a pair-wise manner (Fig. ##FIG##0##1##). These organisms are:</p> <p><italic>M. tuberculosis</italic> CDC1551 (NC_002755.2)</p> <p><italic>M. tuberculosis</italic> F11 (NC_009565.1)</p> <p><italic>M. tuberculosis</italic> H37Ra (NC_009525.1)</p> <p><italic>M. tuberculosis</italic> strain H37Rv (NC_000962.2)</p> <p><italic>M. bovis</italic> AF2122/97 (NC_002945.3)</p> <p><italic>M. bovis</italic> BCG str Pasteur 1173P2 (NC_008769.1)</p> <p>Each pair of organisms can be analyzed in two different ways by choosing each one as query and the other one as target genome. We recommend that a pair should be analyzed in both ways in order to get a comprehensive list of variations, particularly indels. After selection of organisms the type of variation needs to be selected from the search page. Currently there are four options available and one of these to be chosen among:</p> <p>SNP, Insertion, Repeat expansion, Divergent regions</p> <p>After submission of the selected information a detailed query page appears. For example, if SNP is selected the new page will ask for choosing one of the 20 different possible transitions in a user-defined menu-bar and the search can be made restrictive by specifying genomic coordinates or gene name (Fig. ##FIG##0##1##). The output would show all the indicated SNPs in the selected region along with annotation of genes that contain the SNPs (Fig. ##FIG##0##1##). For insertions, divergent regions and repeat expansion the query page has also the option of selecting output on the basis of size in nucleotides. There are four options at present and these are &gt;10, 10–50, 50–100 and &lt;100. Since one can select only one query at a time, an error message is displayed if more then one query is selected.</p> <p>Table ##TAB##0##1## gives the total data present in the database. However, the distribution of these changes among strains and species are different. For example, the number of SNPs between the two <italic>M. tuberculosis </italic>strains H37Ra and Rv are 588 and that between the two <italic>M. bovis </italic>strains are 1271 (Table ##TAB##1##2##). In general the number of variants, observed between the two <italic>M. tuberculosis </italic>strains were much less compared to that between the two <italic>M. bovis </italic>strains. This is consistent with the fact that <italic>M. tuberculosis </italic>strains H37Ra and Rv have been recently derived from H37 [##UREF##4##12##]. These differences are a result of evolutionary history of the organisms and can be useful to map all the potential mutation hotspots in these organisms.</p>
<title>Conclusion</title> <p>In this report we describe MGDD, a database of genomic variants computed from fully sequenced organisms belonging to the <italic>M. tuberculosis </italic>complex. It contains data pertaining to SNP, insertions, repeat expansion and regions that show sequence divergence. Since MGDD is modular information regarding new genomes can be incorporated as and when the sequences become available. The search tool is simple and user friendly and allows one to locate a specific variation in any part of the genome or a gene.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Variation in genomes among different closely-related organisms can be linked to phenotypic differences. A number of mechanisms, such as replication error, repeat expansion and contraction, recombination and transposition can contribute to genomic differences. These processes lead to generation of SNPs, different types of repeat-based and transposons or IS-element-based polymorphisms, inversions and duplications and changes in synteny. A database of all the variations in a group of organisms is not only useful for understanding genotype-phenotype relationship but also in clinical applications. There is no database available at present that provides information about detailed genomic variations among different strains and species of <italic>Mycobacterium tuberculosis </italic>complex, organisms responsible for human diseases.</p> <title>Description</title> <p>MGDD is a free web-based database that allows quick user friendly search to find different types of genomic variations among a group of fully sequenced organisms belonging to <italic>M. tuberculosis </italic>complex. The searches are based on data generated by pair wise comparison using a tool that has already been described. Different types of variations that can be searched are SNPs, indels, tandem repeats and divergent regions. The searches can be designed to find specific variations either in a given gene or any given location of the query genome with respect to any other genome currently available.</p> <title>Conclusion</title> <p>Web-based database MGDD can help to find all the possible differences that exists between two strains or species of <italic>M. tuberculosis </italic>complex. The search tool is very user-friendly and can be used by anyone not familiar with computational methods and will be useful to both clinicians and researchers working on tuberculosis and other Mycobacterial diseases.</p>
<title>Construction and content</title> <p>The database contains pre-computed data derived from full genome sequences of <italic>M. tuberculosis </italic>strains H37Rv, CDC1551, H37Ra, F11, <italic>Mycobacterium bovis </italic>AF2122/97 and <italic>M. bovis </italic>BCG str. Pasteur 1173P2 using ABWGC [##UREF##3##10##]. The variations are categorized as SNPs, insertions, divergent regions (based on lack of sequence identity) and tandem repeats. All computations have been carried out in a pair-wise fashion. In some cases, such as SNPs the results differ depending upon the genome that has been used as a query in a pair of genomes. The database contains two sets of data pertaining to using each genome as query sequence. Insertion in one genome can be considered as a deletion in another genome, so the database contains only the insertions. If the insertions are due to known insertion elements and phage sequences these have been pointed out so that the information can be used for devising methods for better diagnosis and strain identification. Tandem repeats were identified using ABWGC and verified by "Tandem Repeat Finder" [##REF##9862982##11##].</p> <p>MGDD is implemented by using three- tier architecture. The web based application is created by using Apache web server which is connected to the database using MYSQL through an application layer written in Perl-CGI.</p> <p>The information from MGDD can be obtained by selecting a specific query using the "search option" given in the MGDD browser.</p> <title>Availability and requirements</title> <p>The web server can be accessed at <ext-link ext-link-type="uri" xlink:href="http://mirna.jnu.ac.in/mgdd/"/>.</p> <p>MTBreg: <ext-link ext-link-type="uri" xlink:href="http://www.doe-mbi.ucla.edu/Services/MTBreg/"/></p> <p>MycoDB: <ext-link ext-link-type="uri" xlink:href="http://xbase.bham.ac.uk/mycodb/about.pl"/></p> <p>MycoperonDB: <ext-link ext-link-type="uri" xlink:href="http://www.cdfd.org.in/mycoperondb/index.html"/></p> <p>GenoMycDB: <ext-link ext-link-type="uri" xlink:href="http://157.86.176.108/~catanho/genomycdb/"/></p> <title>Authors' contributions</title> <p>AV has implemented the ABWGC and identified the different divergent regions in two genomes and tested the database. AS has made the database and web based application. RR and AB helped in conceptualizing the database and writing the manuscript. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>The authors thank Department of Biotechnology (COE and TB informatics project) for support.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>A typical output of a query (SNP)</bold>. The transition selected was 'gt' and <italic>M. tuberculosis </italic>H37Rv was compared with <italic>M. tuberculosis </italic>CDC1551.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Statistics and data composition of MGDD</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Type of Data</bold></td><td align="left"><bold>Total number of entries.</bold></td></tr></thead><tbody><tr><td align="left">Divergence</td><td align="left">829</td></tr><tr><td align="left">Insertions</td><td align="left">7865</td></tr><tr><td align="left">Repeat expansion</td><td align="left">578</td></tr><tr><td align="left">SNP</td><td align="left">68768</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Genomic variants in <italic>M. tuberculosis </italic>and <italic>M. bovis </italic>strains</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left"><italic>M. tuberculosis </italic>H37Rv compared to H37Ra</td><td align="left"><italic>M. bovis </italic>compared to <italic>M. bovis </italic>strain BCG</td></tr></thead><tbody><tr><td align="left">SNP</td><td align="left">582</td><td align="left">1272</td></tr><tr><td align="left">Divergent Region</td><td align="left">1</td><td align="left">4</td></tr><tr><td align="left">Repeat expansion</td><td align="left">3</td><td align="left">12</td></tr><tr><td align="left">Insertions</td><td align="left">22</td><td align="left">114</td></tr></tbody></table></table-wrap>
<graphic xlink:href="1471-2164-9-373-1"/>
[{"surname": ["Jackson", "Athanassopoulos", "Tsiamis", "Mansfield", "Sesma", "Arnold", "Gibbon", "Murillo", "Taylor", "Vivian"], "given-names": ["RW", "E", "G", "JW", "A", "DL", "MJ", "J", "JD", "A"], "article-title": ["Identification of a pathogenicity island, which contains genes for virulence and avirulence, on a large native plasmid in the bean pathogen "], "italic": ["Pseudomonas syringae "], "source": ["Proceeding of the National Academy of Sciences USA"], "year": ["1999"], "volume": ["96"], "fpage": ["10875"], "lpage": ["10880"], "pub-id": ["10.1073/pnas.96.19.10875"]}, {"surname": ["Fomukong", "Tang", "Maamary", "Ibrahum", "Ramayah", "Yates", "Zainuddin", "Dale"], "given-names": ["NG", "TH", "al", "WA", "S", "M", "ZF", "JW"], "article-title": ["Insertion sequence typing of "], "italic": ["Mycobacterium tuberculosis"], "source": ["Tuber Lung Disease"], "year": ["1994"], "volume": ["75"], "fpage": ["435"], "lpage": ["40"], "pub-id": ["10.1016/0962-8479(94)90117-1"]}, {"surname": ["Chaudhuri", "Pallen"], "given-names": ["RR", "MJ"], "article-title": ["xBASE, a collection of online databases for bacterial comparative genomics"], "source": ["Nucleic Acid Research"], "year": ["2006"], "volume": ["34"], "fpage": ["D335"], "lpage": ["337"], "pub-id": ["10.1093/nar/gkj140"]}, {"surname": ["Vishnoi", "Roy", "Bhattacharya"], "given-names": ["A", "R", "A"], "article-title": ["Comparative analysis of bacterial genomes: identification of divergent regions in mycobacterial genomes using an anchor based approach"], "source": ["Nucleic Acid Research"], "year": ["2007"], "volume": ["35"], "fpage": ["3654"], "lpage": ["3667"], "pub-id": ["10.1093/nar/gkm209"]}, {"surname": ["Steenken", "Gardner"], "given-names": ["W", "LV"], "article-title": ["History of H37 strain of tubercle bacillus"], "source": ["Amer Rev Tuberc"], "year": ["1946"], "volume": ["54"], "fpage": ["52"], "lpage": ["66"]}]
BMC Genomics. 2008 Aug 5; 9:373