1 / 56

DNA sequence analysis Gene prediction

DNA sequence analysis Gene prediction. Gene prediction methods Gene indices Mapping cDNA on genomic DNA Applications. DNA sequence analysis Gene prediction. Protein coding sequence. 3‘UTR. 5‘UTR. promotor. exon 1. exon 2. exon n. exon n-1.

ramona
Download Presentation

DNA sequence analysis Gene prediction

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. DNA sequence analysisGene prediction • Gene prediction methods • Gene indices • Mapping cDNA on genomic DNA • Applications Computational Molecular Biology MPI for Molecular Genetics

  2. DNA sequence analysisGene prediction Protein coding sequence 3‘UTR 5‘UTR promotor exon 1 exon 2 exon n exon n-1 Computational Molecular Biology MPI for Molecular Genetics

  3. Gene predictionStrategies for detecting ORFs / exons • Distribution of Stop-codons • Codon usage • Hexamer frequencies • Prediction of the coding frame • Splice site recognition (Eucaryotes only) Computational Molecular Biology MPI for Molecular Genetics

  4. Gene predictionCodon usage (single exon) coding Frame 1 non-coding Frame 2 Frame 3 Computational Molecular Biology MPI for Molecular Genetics

  5. coding sequence Gene predictionCodon usage (single exon) coding Frame 1 non-coding Frame 2 Frame 3 correct start Computational Molecular Biology MPI for Molecular Genetics

  6. Gene predictionCodon usage (multiple exons) Exons: 208. .295 1029. .1349 1500. .1688 2686. .2934 3326. .3444 3573. .3680 4135. .4309 4708. .4846 4993. .5096 7301. .7389 7860. .8013 8124. .8405 8553. .8713 9089. .9225 13841. .14244 coding Frame 1 non-coding Frame 2 Frame 3 Splice sites Computational Molecular Biology MPI for Molecular Genetics

  7. Gene predictionCodon usage (multiple exons) Exons: 208. .295 1029. .1349 1500. .1688 2686. .2934 3326. .3444 3573. .3680 4135. .4309 4708. .4846 4993. .5096 7301. .7389 7860. .8013 8124. .8405 8553. .8713 9089. .9225 13841. .14244 coding Frame 1 non-coding Frame 2 Frame 3 Splice sites Computational Molecular Biology MPI for Molecular Genetics

  8. Gene predictionAdditional criteria • Detection of start codons • Detection of potential promotor elements • Detection of repetitive sequences (mostly untranslated) • Homology to known genes of related organisms Computational Molecular Biology MPI for Molecular Genetics

  9. Gene predictionSoftware • GENSCAN (C.Burge & S.Karlin) • Grail (neural network; Ueberbacher et al.) • MZEF (M. Zhang,1997) • FGeneH, Hexon (V.Solovyev et al., 1994) • Genie, etc. All programs are using dynamic programming for detection of the optimal solution Computational Molecular Biology MPI for Molecular Genetics

  10. DNA sequences in public databases Human ~ 4 million ESTs + 130 000 RNAs Mouse ~ 2.7 million ESTs + 30 000 RNAs Computational Molecular Biology MPI for Molecular Genetics

  11. cDNA is usually oligo dT primed, or by random primers • Several cDNAs for the same mRNA may be generated AAAAAA... cDNA TTTTTT... Expressed sequence tags (EST) • Reverse transcriptase stops ‚randomly‘ mRNA Computational Molecular Biology MPI for Molecular Genetics

  12. Expressed sequence tags (EST) Dechiffered sequence (EST) Clone = mRNA fragment 3‘-primer <700 bp Vector (known sequence) Average: 1500 bp Computational Molecular Biology MPI for Molecular Genetics

  13. Expressed sequence tags (EST) • Isolation of mRNAs from tissue(s) • Generation of cDNAs reflecting parts of the RNAs • Cloning of cDNAs into a vector (often random orientation) • End sequencing of the clones Computational Molecular Biology MPI for Molecular Genetics

  14. Generation of ESTsbasecalling problems close to 5‘ end of EST close to 3‘ end of EST missing bases Computational Molecular Biology MPI for Molecular Genetics

  15. expressed sequence tags (ESTs) Coverage of an mRNA by ESTs putative mRNA AAAAAA... 5‘UTR exon 1 exon 2 3‘UTR Computational Molecular Biology MPI for Molecular Genetics

  16. Characteristics of ESTs • Highly redundant • Low sequence quality • (Cheap) • Reflect expressed genes • May be tissue/stage specific Computational Molecular Biology MPI for Molecular Genetics

  17. Gene indices Clustering of EST and mRNA sequences of an organism to reduce redundance in sequence data. Goal: Each cluster represents one gene or mRNA • UniGene (NCBI) • TIGR Gene Indices • STACK (SANBI) • GeneNest (DKFZ,MPI) Computational Molecular Biology MPI for Molecular Genetics

  18. Gene indicesGeneNest workflow EMBL database Unigene database Quality clipping Quality clipping BLAST/QUASAR search, clustering Assembly, Consensus sequences Visualization Computational Molecular Biology MPI for Molecular Genetics

  19. Gene indicesQuality clipping • Removal of vector sequence • Masking of repetitive sequences (e.g. Alu) • Removal of terminal sequences of low quality In order to cluster based on gene-specific sequence data the following steps have to be performed: Computational Molecular Biology MPI for Molecular Genetics

  20. Gene indicesClustering • Minimal % identity (e.g. > 95%) • Minimal length of match (e.g. >40 bp) • No internal matches (TIGR gene indices) • Same origin of tissue (only STACK) Sequences are usually clustered if the matching part between two sequences fullfills several (empirical) criteria: Computational Molecular Biology MPI for Molecular Genetics

  21. Gene indicesAssembly • Contigs, reflecting parts of different transcripts • One consensus sequence per contig • A relative order of the sequences (alignment) Sequences in a cluster are assembled to group those sequences which are globally similar, resulting in Computational Molecular Biology MPI for Molecular Genetics

  22. Gene indicesConsensus sequences • Reduced error rate • Consensus often longer than any single sequence contributing • Efficient database search • Detection of exon/intron boundaries and alternative splice variants Computational Molecular Biology MPI for Molecular Genetics

  23. Gene indices Alignment consensus Computational Molecular Biology MPI for Molecular Genetics

  24. Gene indices AlignmentSoftware • Phrap (Phil Green) • CAP3 (X. Huang) • TIGR assembler • GAP4 (R. Staden) Computational Molecular Biology MPI for Molecular Genetics

  25. GeneNest visualization(http://genenest.molgen.mpg.de) Computational Molecular Biology MPI for Molecular Genetics

  26. GeneNest visualization(http://genenest.molgen.mpg.de) Computational Molecular Biology MPI for Molecular Genetics

  27. TIGR Gene Indices(http://www.tigr.org/) Alignment scheme Computational Molecular Biology MPI for Molecular Genetics

  28. UniGene(http://www.ncbi.nih.nlm.gov/UniGene) Computational Molecular Biology MPI for Molecular Genetics

  29. UniGene(http://www.ncbi.nih.nlm.gov/UniGene) Computational Molecular Biology MPI for Molecular Genetics

  30. missing intron consensus sequence (  mRNA) exons Mapping of consensus sequences on genomic DNA genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  31. Mapping cDNA on genomic DNA Computational Molecular Biology MPI for Molecular Genetics

  32. Gene indicesApplications • Detection of exon/intron boundaries • Detection of alternative splicing • Detection of Single Nucleotide Polymorphisms • Genome annotation • Analysis of gene expression • Genome-genome comparison Computational Molecular Biology MPI for Molecular Genetics

  33. 5‘UTR exon 1 exon 2 exon 3 mRNA 1 5‘UTR exon 1 exon 3 mRNA 2 5‘UTR exon 1 exon 2 Alternative Splicing hnRNA Computational Molecular Biology MPI for Molecular Genetics

  34. Alignment of EST consensus sequences and genomic target genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  35. Detection of the appropriate genomic target sequence Local similarity of EST consensus and genomic DNA >96% identity genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  36. Cutting out genomic target sequence genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  37. splice variant consensus sequence (  mRNA) exons Alternative Splicing(mapping on genomic DNA) genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  38. putative exons alternative exon aligned GeneNest consensus genomic sequence SpliceNest(http://SpliceNest.molgen.mpg.de) Computational Molecular Biology MPI for Molecular Genetics

  39. Alternative Splicing(additional exon) Splice variants of adenylsuccinate lyase unspliced ? skipped exon gene prediction errors ? Computational Molecular Biology MPI for Molecular Genetics

  40. Alternative Splicing Splice variants of APECED gene alternative variants number of sequences genomic sequence Computational Molecular Biology MPI for Molecular Genetics

  41. Alternative splicing Computational Molecular Biology MPI for Molecular Genetics

  42. Alternative Splicing (alternative donor site) Computational Molecular Biology MPI for Molecular Genetics

  43. Alternative Splicing Computational Molecular Biology MPI for Molecular Genetics

  44. Alternative Splicing(alternative exons) Computational Molecular Biology MPI for Molecular Genetics

  45. SpliceNest(hypothetical gene Hs16936) Computational Molecular Biology MPI for Molecular Genetics

  46. Single Nucleotide Polymorphisms(SNP) • SNPs are single base differences within one species • Several million SNPs detected in Human • SNPs may be related to diseases Computational Molecular Biology MPI for Molecular Genetics

  47. Single Nucleotide Polymorphisms(SNP) SNP or basecalling error ? Computational Molecular Biology MPI for Molecular Genetics

  48. Genome Annotation / Ensembl(http://www.ensembl.org) Computational Molecular Biology MPI for Molecular Genetics

  49. Analysis of gene expressiontissue-specificity • Counting frequency of EST derived from a specific tissue within one sequence cluster • Searching for cluster/contigs which are tissue specific (e.g. tumor) • Searching for alternative splice variants which are potentially tissue specific Computational Molecular Biology MPI for Molecular Genetics

  50. liver tumor kidney Analysis of gene expressionPDZ-domain containing protein PDZK1 (Hs.15456) Computational Molecular Biology MPI for Molecular Genetics

More Related