1 / 22

Double-Ended Shotgun Sequencing of PA14

Double-Ended Shotgun Sequencing of PA14. Daniel G. Lee 10/30/02. Determination of PA14 Genomic Sequence and Whole-Genome Alignment with PAO1. The complete genome of a related P. aeruginosa strain, PAO1, has been determined. The genome size is 6.2 Mb.

hmenzel
Download Presentation

Double-Ended Shotgun Sequencing of PA14

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Double-Ended Shotgun Sequencing of PA14 Daniel G. Lee 10/30/02

  2. Determination of PA14 Genomic Sequence and Whole-Genome Alignment with PAO1 • The complete genome of a related P. aeruginosa strain, PAO1, has been determined. The genome size is 6.2 Mb. • PAO1 is less virulent than PA14 in almost all of our model hosts. • PA14 contains additional DNA (sometimes large islands of DNA) not found in PAO1. Some of these additional genes may be responsible for the enhanced virulence of PA14. • A complete PA14 genomic sequence will allow us to: • identify all the (DNA) differences between PA14 and PAO1 (and later evaluate their contribution to virulence). • Simplify the bioinformatics component of the PA14 Unigene library. • Design a microarray (whole-genome or PA14-specific).

  3. PA14 Sequencing - Outline • Sequencing workflow. • Finishing. • Annotation and whole-genome alignment. • Integration with PA14 insertion library. • Requirements for publication.

  4. PA14 Genomic DNA Prep • Original PA14 RifR isolate from LGR. • 500 ml culture. • Alkaline lysis, with: • CTAB ppt. • 2 x Chloroform/Isoamyl Alcohol extraction. • 3 x Phenol extraction. • 1 x Phenol/Chloroform/Isoamyl Alcohol extraction. • 1 x Chloroform/Isoamyl Alcohol extraction. • Isoamyl alcohol ppt. • Resupended in 5 ml TE @1.174 mg/ml (5.87 mg total).

  5. PA14 Genomic DNA Prep Plasmid preps of PA14 library. Ligate PA14 fragments into vector and transform E. coli Sequencing of amplification products Linear amplification of inserts using dideoxy terminators Shear PA14 DNA and Size Fractionate Contig assembly (PHRED and PHRAP) • Genome-wide alignment with PAO1 • order contigs. • identify gaps for sequence finishing. • identify differences between PA14 and PAO1. Finishingand Annotation Workflow for Double-Ended Shotgun Sequencing of PA14

  6. Plasmid Library Construction • Shear DNA using nitrogen (cleavage more random than sonication). • Fill-in to produce blunt ends. • Size fractionate on low-melt agarose gel. • 1-3 kb fragments (700 bp). • 3-7 kb fragments • Ligate. • Transform. • Pick colonies.

  7. Plasmid Preps • O/N cultures in 96-well plates. • Freeze cell pellets. • Alkaline-lysis mini-preps in 96-well plates. • 604/650 plates done. • Dry DNA pellets O/N. • Resuspend DNA in H2O. • Transfer to 384-well plate. • QC by agarose gel.

  8. Sequencing Reactions • Set up reaction mix: • Labelled ddNTPs. • dNTPs • Buffer • Taq • Forward or reverse sequencing primer. • Aliquot rxn mix to 384-well PCR plate; freeze. • Add 3 ml DNA to each well (or 3 ml vector for “PCR control”). • “PCR”

  9. DNA Sequencing EtOH ppt. PCR reactions. Dry. Rssp. in H2O. Add previously characterized PCR reactions as “sequencing controls”. ABI Prism sequencers (liquid polymer capillary sequencer, 96 reactions at a time).

  10. ABI sequencer outputs electropherograms. • PHRED determines identity of base as well as quality score.

  11. Contig Assembly Electropherograms. PHRED - determines base identity and quality score for each position. PHRAP - aligns sequences to assemble contigs, determines consensus sequence and quality score for each position.

  12. CCG-AATTCCGGCTTTACGACGACTTGGGGTACCA ccg-aatt-cggctttacg aatgccggcattacg tt--ggc-ttacgaccctttg-ggt t--ggc-ttacg--gactaggggtacca Contig Assembly

  13. PA14 Sequencing:Current Status (as of 10/17/02) • Total amount of sequence: ~ 6 Mb (6.5 X coverage) • 72,000 sequences (36,000 clones). • 390 (out of ~650) 96-well plates sequenced. • 604 plates mini-prepped • Total number of “contigs”: < 2000? • 1 contig ~ 44 kb • 1 contig ~ 35 kb • ~ 10 contigs > 25 kb • ~ 12 contigs 20-25 kb • most contigs are 5-10 kb. • Library consists of ~1 kb inserts (current plans to introduce a library of 3-6 kb inserts). As of 10/21/02, one contig 73 kb, many > 50 kb.

  14. PA14 Genomic DNA Prep Plasmid preps of PA14 library. Ligate PA14 fragments into vector and transform E. coli Sequencing of amplification products Linear amplification of inserts using dideoxy terminators Shear PA14 DNA and Size Fractionate Contig assembly (PHRED and PHRAP) • Genome-wide alignment with PAO1 • order contigs. • identify gaps for sequence finishing. • identify differences between PA14 and PAO1. Finishingand Annotation Workflow for Double-Ended Shotgun Sequencing of PA14

  15. B A A: gaps in regions corresponding to PAO1 sequence. B: gaps in PA14-specific regions. Comparisons of PAO1 and PA14

  16. Tools for Genome-Wide Alignments of PAO1 and PA14 • Software Packages available from TIGR for Alignments. • MUMmer 2.1 - aligns MUMs (maximal unique matches) for two input sequences (two 3-4 Mb genomes aligned in under 30 seconds, using less than 100 Mb of memory, on a typical desktop computer running Unix/Linux). • NUCmer - alignments of highly similar sequences that may have large rearrangements (i.e. -- a group of assembly contigs vs. a complete genome). • PROmer - amino acid translation in all 6 frames for protein/peptide alignments. Useful for comparative genome annotation. • DisplayMUMs for graphical analysis of MUMmer output.

  17. Tools for Annotation of PA14 • PROmer - amino acid translation in all 6 frames for protein/peptide alignments. Useful for comparative genome annotation. • Jonathan’s automated suite of annotation tools (Hrp project)

  18. Approaches for Finishing • PCR amplification and directed sequencing of gapped regions. • Isolation of cosmid clones spanning gaps, subcloning, sequencing of subclones (using universal primers). • (Direct genomic sequencing). • (Altering sequencing reaction conditions for regions that are difficult to sequence through).

  19. B A Finishing • Methods: • PCR. • Cosmids • (Directed sequencing) • (Altered rxn. Conditions) • Considerations: • Type of gap. • Anticipated size of gap. • Quality/nature of sequence at junction.

  20. Integration with PA14 Unigene Library

  21. Requirements for Publication • “Finished” PA14 sequence. • Sufficient quality. • No gaps? • Annotation. • Comparison to PAO1. 4. What else? • Virulence data? • Proteomics? • Others?

  22. ACKNOWLEDGEMENTS • Harvard Partners • Genome Center: • K. Montgomery • G. Grills • L. Li • W. Brown • J. Decker • R. Elliot • L. Gendal • K. Osborn • A. Parerra • C. Xi • P. Juels • R. Kucherlapati MGH: N. Liberati S. Miyata J. Urbach F. Ausubel X. He M. Saucier L. Rahme

More Related