Genome Assemblies: Difference between revisions

From EchinoWiki
Jump to navigation Jump to search
imported>Echinobase
imported>Echinobase
(7 intermediate revisions by the same user not shown)
Line 19: Line 19:




Pm v1.0 Pm v2.0
Scaffold number 60,183 57,698
Scaffold N50 52,6141 76,341
Contig number 179,756 131,779
Contig N50 9,466 18,676


<table style="width:100%">
{| class="wikitable"
  <tr>
!
    <th>Firstname</th>
! Pm v1.0
    <th>Lastname</th>
! Pm v2.0
    <th>Age</th>
|-
  </tr>
| Scaffold number
  <tr>
| 60,183
    <td>Jill</td>
| 57,698
    <td>Smith</td>
|-
    <td>50</td>
| Scaffold N50
  </tr>
| 52,6141
  <tr>
| 76,341
    <td>Eve</td>
|-
    <td>Jackson</td>
| Contig number
    <td>94</td>
| 179,756
  </tr>
| 131,779
  <tr>
|-
    <td>John</td>
| Contig N50
    <td>Doe</td>
| 9,466
    <td>80</td>
| 18,676
  </tr>
|}
</table>


=== '''V1.0 Assembly'''===
=== '''V1.0 Assembly'''===
<u>What's New</u>
Pmin_1.0 is the latest (as of Apr 11, 2012) assembly of the genome of Patiria Miniata. The assembly tools CABOG (Celera Assembler), Newbler, ATLAS-Link, and ATLAS-GapFill were used to assemble a combination of 454 reads (fragment and 2.5kb insert paired ends;~15x coverage) and Illumina reads (300bp insert and 2.5kb insert paired ends;~70x coverage).
<u>Introduction</u>
This information is for the first release (Pmin_1.0) of the draft genome sequence of the Patiria miniata . This is a draft sequence and may contain errors so users should exercise caution.Typical errors in draft genome sequences include misassemblies of repeated sequences, collapses of repeated regions, and unmerged overlaps(e.g. due to polymorphisms) creating artificial duplications.
With a goal of solving the polymorphism issues of the data while maintaining the sequence continuity, The Pmin_1.0 assembly was generated in the following steps:
1) 454 reads were assembled by CABOG using settings less strignent than the default (unitigger=bog utgErrorRate=0.03 ovlErrorRate=0.08 cnsErrorRate=0.08 cgwErrorRate=0.14 doExtendClearRanges=0)
2) Both contig and degenerate sequences from the previous step were chopped into fake reads with ~11x coverage (500bp long; 460bp overlap; 80bp minimal length) for ctgs and 8x coverage(450bp long; 400bp overlap; 80bp minimal length) for degs. The fake reads were then assembled by Newbler with the option of -large.
3) Both 454 and iIlumina pair end reads were mapped to the contigs from the previous step. We used BLAT to map the 454 data and bwa(aln+samse) to map the Illumina data, both with the default options. Based on the mapping locations of the paired ends, contigs were then ordered and oriented into scaffolds using ATLAS-Link.
4) ATLAS-GapFill was then used to assemble the reads locally in an attempt to fill the gaps among the contigs within the scaffolds.This final step produced 770.5Mb sequences with contig N50 size of 9.5kb and scaffold N50 size of 50.3kb.
<u>Conditions for use</u>
These data are made available before scientific publication with the following understanding:
*The data may be freely downloaded, used in analyses, and repackaged in databases.
*Users are free to use the data in scientific papers analyzing particular genes and regions if the providers of this data (Baylor College of Medicine Human Genome Sequencing Center) are properly acknowledged. Please cite the BCM-HGSC web site or publications from BCM-HGSC referring to the genome sequence.
*The BCM-HGSC plans to publish the assembly and genomic annotation of the dataset, including large-scale identification of regions of evolutionary conservation and other features.
*This is in accordance with, and with the understandings in the Fort Lauderdale meeting discussing Community Resource Projects and the resulting NHGRI policy statement (http://www.genome.gov/page.cfm?pageID=10506537).
*Any redistribution of the data should carry this notice.


== '''''Lytechinus variegatus''''' ==
== '''''Lytechinus variegatus''''' ==

Revision as of 14:41, 4 December 2019

Echinoderm Genome Assemblies by Species


Strongylocentrotus purpuratus

Assembly_3.1 (Spur_3.1)

Assembly 2.6(Spur 2.6)

Assembly_2.5(Spur_2.5)

Assembly_2.1(Spur_2.1)

Assembly_0.5(Spur_0.5)

Patiria miniata

V2.0 Assembly

We sought to improve the Patiria miniata genome assembly with additional PacBio sequences. We generated a new PacBio read dataset at the Duke University Sequencing Center using our reference individual DNA. The read dataset contains 2 million reads and 15.8 billion bp. The read N50 is 10.4 Kb. We used PBJelly2 to combine the PacBio reads with the previously assembled contigs. The results were an improvement in contig size and number with only a small reduction in the number of scaffolds (Table). The P. miniata Gene v2.0 set was generated using MAKER2 pipeline from v2.0 genome assembly.


Pm v1.0 Pm v2.0
Scaffold number 60,183 57,698
Scaffold N50 52,6141 76,341
Contig number 179,756 131,779
Contig N50 9,466 18,676

V1.0 Assembly

What's New

Pmin_1.0 is the latest (as of Apr 11, 2012) assembly of the genome of Patiria Miniata. The assembly tools CABOG (Celera Assembler), Newbler, ATLAS-Link, and ATLAS-GapFill were used to assemble a combination of 454 reads (fragment and 2.5kb insert paired ends;~15x coverage) and Illumina reads (300bp insert and 2.5kb insert paired ends;~70x coverage).

Introduction

This information is for the first release (Pmin_1.0) of the draft genome sequence of the Patiria miniata . This is a draft sequence and may contain errors so users should exercise caution.Typical errors in draft genome sequences include misassemblies of repeated sequences, collapses of repeated regions, and unmerged overlaps(e.g. due to polymorphisms) creating artificial duplications.

With a goal of solving the polymorphism issues of the data while maintaining the sequence continuity, The Pmin_1.0 assembly was generated in the following steps:

1) 454 reads were assembled by CABOG using settings less strignent than the default (unitigger=bog utgErrorRate=0.03 ovlErrorRate=0.08 cnsErrorRate=0.08 cgwErrorRate=0.14 doExtendClearRanges=0)

2) Both contig and degenerate sequences from the previous step were chopped into fake reads with ~11x coverage (500bp long; 460bp overlap; 80bp minimal length) for ctgs and 8x coverage(450bp long; 400bp overlap; 80bp minimal length) for degs. The fake reads were then assembled by Newbler with the option of -large.

3) Both 454 and iIlumina pair end reads were mapped to the contigs from the previous step. We used BLAT to map the 454 data and bwa(aln+samse) to map the Illumina data, both with the default options. Based on the mapping locations of the paired ends, contigs were then ordered and oriented into scaffolds using ATLAS-Link.

4) ATLAS-GapFill was then used to assemble the reads locally in an attempt to fill the gaps among the contigs within the scaffolds.This final step produced 770.5Mb sequences with contig N50 size of 9.5kb and scaffold N50 size of 50.3kb.

Conditions for use

These data are made available before scientific publication with the following understanding:

  • The data may be freely downloaded, used in analyses, and repackaged in databases.
  • Users are free to use the data in scientific papers analyzing particular genes and regions if the providers of this data (Baylor College of Medicine Human Genome Sequencing Center) are properly acknowledged. Please cite the BCM-HGSC web site or publications from BCM-HGSC referring to the genome sequence.
  • The BCM-HGSC plans to publish the assembly and genomic annotation of the dataset, including large-scale identification of regions of evolutionary conservation and other features.
  • Any redistribution of the data should carry this notice.

Lytechinus variegatus

Assembly LvPtE5C

Assembly LvMSCB

Assembly 2.2 (Lvar_2.2)

Assembly 0.4 (Lvar_0.4)