Protocol Exchange

Author Information

Subject Terms

Keywords

Abstract

Introduction & Materials

Procedures

Outcomes

Conflicting Financial Interests

Associated Publications

Supplementary Files

Figures

Qiang Wang

Jian-Qun Chen

Dacheng Tian

Figure 1

Figure 2

Figure 3

AlignDB: a computational protocol for analyzing the relationship between indels and SNPs in genome sequences

AlignDB is developed to investigate the distribution of single-nucleotide changes around insertion/deletions \(indels) in genome comparisons. There are two set of analysis, two-way and three-way, in AlignDB.





  The two-way analysis indicates that nucleotide divergence \(D) is substantially elevated surrounding indels and decreases monotonically to near-background levels over several hundred bases. D is significantly correlated with both size and abundance of nearby indels. In a comparison of closely related species, the three-way analysis is available. We find that derived nucleotide substitutions surrounding indels occur in significantly greater numbers on the lineage containing the indel than on the one containing the ancestral \(non-indel) allele; the same holds within species for single-nucleotide mutations surrounding polymorphic indels.





  AlignDB is freely available on request from the authors. The parameters are fully modifiable. The following protocol describes how to use AlignDB to do two-way and three way genomic analyses for indels.

Hardware:


* A machine capable of running a modern OS. \(e.g. Window 2000 or newer, Linux kernel 2.6 or newer, OS X, etc.)
  * 1 Gb of RAM.
  * 80 Gb of hard drive space.
  


Software:


* Perl 5.8.x or better.
  * MySQL 5.0.x.
  * ClustalW alignment program.
  * Microsoft Excel or OpenOffice.org.
  * AlignDB itself.
  


Data:


* Whole genome alignment.
  * Genomic annotation. We use Ensembl as the genomic annotation source.

1. Install all AlignDB requirement software and dependent Perl modules. Read the README file in the AlignDB package root directory.
  


2. Install AlignDB itself. Follow the steps in doc/tutorial.pdf and use the example data in data/ directory coming with AlignDB for a trial run.
  


3. Prepare a set of input data in Blastz .axt file format. For example, you can download the completed genomic alignments from the UCSC Genome Bioinformatics Site.
  


4. Starting the GUI shell. Enter the base directory where you installed AlignDB. Type the follow command:
      &gt;perl gui\gui2.pl


    The GUI interface shows in Fig.1. DO NOT close the command line window, AlignDB need it.





5. Two-way comparison steps.
      Step 1: test your server. Click on the &#x201C;DB Server&#x201D; notebook tab, and check MySQL server setting. 


    Step 2: initiate your database. Then click the &#x201C;Init. alignDB&#x201D; button, which will call init\init_alignDB.pl to build the skeleton of the database.


    Step 3: generate two-way analysis database. Click the &#x201C;Gene. alignDB&#x201D; button which will call init\gene_alignDB.pl to analysis alignments and build the database.


    Step 4: update features using Ensembl annotations. Click the &#x201C;Feature&#x201D; tab, change nothing and click the &#x201C;Upd. feature&#x201D; button, which will run init\update_feature.pl.


    Step 5: update misc features. Click the &#x201C;Misc&#x201D; tab. Click &#x201C;Upd. slippage&#x201D; to run init\update_indel_slippage.pl and then click &#x201C;Upd. isw&#x201D; to run init\update_isw_indel_id.pl.


    Step 6: do statistics and get results. Then click &#x201C;Common stat&#x201D; button. An Excel workbook with 25 worksheets will be generated.


    Step 7: charting. If you run AlignDB on Windows, just like me, you could click the &#x201C;Common chart&#x201D; button, which will run stat\common_chart_factory.pl that call Excel to draw charts in .xls files by OLE. Fig.2 is an example workbook.





6. Three-way comparison steps.
  Three-way comparison needs three genome sequences to accomplish. So you should add a new genome sequence into the previous two-way comparison. This can be accomplished by joining two two-way comparisons, which share a common target sequence.


    Step 1: build another two-way comparison. You can just do steps 1-3 and omit step 4-7 if you don&#x2019;t need statistical analysis for the new two-way database.


    Step 2: generate three-way analysis database. Load two completed databases into &#x201C;first db&#x201D; and &#x201C;second db&#x201D;. Generate &#x201C;goal db&#x201D; name with the default naming rule. Leave the &#x201C;T/Q/R:&#x201D; options unchanged, and click the &#x201C;Ref outgroup&#x201D; button. See Fig.3 for an example.


    Step 3: update two-way features. Three-way comparison database is also a legal two-way database. Do step 4-5 of two-way comparison to update the new three-way database.


    Step 4: update three-way features. Go to &#x201C;Misc&#x201D; tab in the &#x201C;Three-way&#x201D; pane and click the &#x201C;Upd. CpG&#x201D; button to run extra\update_snp_cpg.pl.


    Step 5: do statistics and get results. Then click &#x201C;Three stat&#x201D; button which will run stat\three_stat_factory.pl. An Excel workbook with 11 worksheets will be generated.


    Step 6: charting. This step is almost the same as step 7 in Two-way comparison. Just click the &#x201C;Three chart&#x201D; button, which will run stat\three_chart_factory.pl.

None

Computational biology and bioinformatics

Biotechnology

Biochemistry

An open repository of community-contributed protocols sponsored by Nature Research.

Learn more about Protocol Exchange

Method Article

Qiang Wang, Jian-Qun Chen, Dacheng Tian

Qiang Wang

Department of Biology, Nanjing University

Jian-Qun Chen

Department of Biology, Nanjing University

Dacheng Tian

Department of Biology, Nanjing University

DOI:

10.1038/nprot.2008.208

Download PDF

Figures

Introduction

Equipment

Procedure

Comments

Version History

CURRENT STATUS: Posted

Version 1

No community comments so far

Metrics

Comments: 0

HTML Views: ...

Subject Areas

Computational biology and bioinformatics

Biotechnology

Biochemistry

Associated Publications

Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes

by D Tian, Q Wang, P Zhang, +7

Nature (24 June, 2008)

An open repository of community-contributed protocols sponsored by Nature Research.

Learn more about Protocol Exchange

Method Article

Qiang Wang, Jian-Qun Chen, Dacheng Tian

Qiang Wang

Department of Biology, Nanjing University

Jian-Qun Chen

Department of Biology, Nanjing University

Dacheng Tian

Department of Biology, Nanjing University

Badges

Peer Review Timeline

Version History

CURRENT STATUS: Posted

Version 1

No community comments so far

Metrics

Metrics

Comments: 0

HTML Views: ...

Subject Areas

Subject Areas

Computational biology and bioinformatics

Biotechnology

Biochemistry

DOI & PDF

DOI:

10.1038/nprot.2008.208

Download PDF

License

Declarations

Figures

Introduction

Equipment

Procedure

Comments

Related Articles

Associated Publications

Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes

by D Tian, Q Wang, P Zhang, +7

Nature (24 June, 2008)

Privacy Policy

Terms of Service