This Article OPEN ACCESS ARTICLE Full Text Full Text (PDF) Supplemental Research Data All Versions of this Article: gr.5894107v1 17/7/1101    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Kim, J. H. Articles by Li, L. M. Search for Related Content PubMed PubMed Citation Articles by Kim, J. H. Articles by Li, L. M. Social Bookmarking                   What's this?

Methods

Diploid genome reconstruction of Ciona intestinalis and comparative analysis with Ciona savignyi

¹ Department of Computer Science, Yonsei University, Seoul, 120-749, Republic of Korea; ² Molecular and Computational Biology Program, Department of Biological Sciences, University of Southern California, Los Angeles, California 90089, USA; ³ Department of Mathematics, University of Southern California, Los Angeles, California 90089, USA

One of the main goals in genome sequencing projects is to determine a haploid consensus sequence even when clone libraries are constructed from homologous chromosomes. However, it has been noticed that haplotypes can be inferred from genome assemblies by investigating phase conservation in sequenced reads. In this study, we seek to infer haplotypes, a diploid consensus sequence, from the genome assembly of an organism, Ciona intestinalis. The Ciona intestinalis genome is an ideal resource from which haplotypes can be inferred because of the high polymorphism rate (1.2%). The haplotype estimation scheme consists of polymorphism detection and phase estimation. The core step of our method is a Gibbs sampling procedure. The mate-pair information from two-end sequenced clone inserts is exploited to provide long-range continuity. We estimate the polymorphism rate of Ciona intestinalis to be 1.2% and 1.5%, according to two different polymorphism counting schemes. The distribution of heterozygosity number is well fit by a compound Poisson distribution. The N50 length of haplotype segments is 37.9 kb in our assembly, while the N50 scaffold length of the Ciona intestinalis assembly is 190 kb. We also infer diploid gene sequences from haplotype segments. According to our reconstruction, 85.4% of predicted gene sequences are continuously covered by single haplotype segments. Our results indicate 97% accuracy in haplotype estimation, based on a simulated data set. We conduct a comparative analysis with Ciona savignyi, and discover interesting patterns of conserved DNA elements in chordates.

E-mail jonghkim{at}usc.edu; fax (213) 740-2437.

[Supplemental material is available at www.genome.org. The diploid genome sequence of Ciona intestinalis is downloadable from http://www-rcf.usc.edu/~lilei/diploid.html. The software to reconstruct haplotypes, called hapBuild, is available on request.]

Article published online before print. Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.5894107

CiteULike    Connotea    Del.icio.us    Digg    Reddit    Technorati    What's this?

This article has been cited by other articles:

				V. Bansal and V. Bafna HapCUT: an efficient and accurate algorithm for the haplotype assembly problem Bioinformatics, August 15, 2008; 24(16): i153 - i159. [Abstract] [PDF]

				V. Bansal, A. L. Halpern, N. Axelrod, and V. Bafna An MCMC algorithm for haplotype assembly from whole-genome sequence data Genome Res., August 1, 2008; 18(8): 1336 - 1346. [Abstract] [Full Text] [PDF]

				Y. Harada, Y. Takagaki, M. Sunagawa, T. Saito, L. Yamada, H. Taniguchi, E. Shoguchi, and H. Sawada Mechanism of Self-Sterility in a Hermaphroditic Chordate Science, April 25, 2008; 320(5875): 548 - 550. [Abstract] [Full Text] [PDF]

				G. Denisov, B. Walenz, A. L. Halpern, J. Miller, N. Axelrod, S. Levy, and G. Sutton Consensus generation and variant detection by Celera Assembler Bioinformatics, April 15, 2008; 24(8): 1035 - 1040. [Abstract] [Full Text] [PDF]