1 / 19

Quality and Error Control Coding for DNA Microarrays

Quality and Error Control Coding for DNA Microarrays. Olgica Milenkovic ECE Department University of Colorado, Boulder. IEEE Denver ComSoc. Outline. DNA Microarrays VLSIPS (Very Large Scale Immobilized Polymer Synthesis) Production of DNA Microarrays (http://www.affymetrix.com/)

mary-kirkland
Download Presentation

Quality and Error Control Coding for DNA Microarrays

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. Quality and Error Control Coding for DNA Microarrays Olgica Milenkovic ECE Department University of Colorado, Boulder IEEE Denver ComSoc

  2. Outline • DNA Microarrays • VLSIPS (Very Large Scale Immobilized Polymer Synthesis) • Production of DNA Microarrays (http://www.affymetrix.com/) • Base Scheduling • Mask Design • Quality-Control Coding • Error-Correcting DNA Microarrays (Multiplexed Arrays) • Production of Multiplexed DNA Microarrays • Base/Color Scheduling • Mask Design • Quality-Control Coding IEEE Denver ComSoc

  3. DNA microarrays I Slide #1 Goal: Determining which genes are expressed (active) and which are unexpressed (inactive) Comparative gene expression study of multiple cells Protein Coding Sequence Protein Coding Sequence Control of Transcription & Translation Transcription Translation Transcription Translation Protein Protein Gene expression and co-regulation IEEE Denver ComSoc

  4. DNA microarrays II Slide #2 Creating the `cell cultures’ to be compared… `Green’ Cell Culture : `Red’ Cell Culture: DNA Subsequence DNA Subsequence 3’- AATTT CGC… - 5’ 5’ - UUAAAGCG… - 3’ 3’ - AATTTCGC… - 5’ 3’ - AATTTCGC… - 5’ 3’- AATTT CGC… - 5’ 3’ - UUAAAGCG… - 5’ 3’- AATTT CGC… - 5’ 3’- AATTT CGC… - 5’ mRNA mRNA cDNA cDNA “Color Coding” “Color Coding” Creation of tagged cDNA sequences from first cell type Creation of tagged cDNA sequences from first cell type IEEE Denver ComSoc

  5. DNA microarrays III Slide #3 Gene Probes Spots Complementary sequences hybridizewith each other, forming stable double-helices Hybridization: 3’-AAGCT-5’ 5’-TTCGA-3’ DNA microarray is scanned by laser light of different wave-lengths IEEE Denver ComSoc

  6. Probe synthesis in microarrays I Slide #4 VLSIPS (Gene Chip, AFFYMETRIX, Array Manufacturing Manual) Linkers Mask Quartz Wafer Linker Activation IEEE Denver ComSoc

  7. Probe synthesis in microarrays II Slide #5 VLSIPS (Gene Chip, AFFYMETRIX, Array Manufacturing Manual) Solution of one DNA base (A or T or G or C) Solution of one DNA base (A) IEEE Denver ComSoc

  8. Slide #6 Base scheduling I Fixed probe length: N Production steps A T G C A T G C A T G C A T G C Spots 1 2 3 4 5 Synchronous schedule (length 4N) CTGA ACAA IEEE Denver ComSoc

  9. Slide #7 Base scheduling II Production steps A G G C T T G C T T G C C C G C Spots 1 2 3 4 5 Asynchronous schedule IEEE Denver ComSoc

  10. Base Scheduling III Slide #8 • Shortest asynchronous base schedule • Shortest common super-sequence of set of M sequences (NP-hard) ESN(M,k) – expected length of a longest common subsequence of M randomly chosen sequences of length N over an alphabet of size k No significant gain for N≈20-30 Periodic schedule used instead (length 4N) IEEE Denver ComSoc

  11. Mask Design Slide #9 Border-length minimization Feldman and Pevzner, 1994 Hannehalli et.al., 2002 Kahng et.al. 2003, 2004 Key idea: Arrange the probes on the array in such a way that the border-length of all masks is minimal Border-length graph: complete graph on M vertices, weight of edges equal to the Hamming distance between probes Greedy traveling salesman algorithm+ threading (discrete space-filing curve) IEEE Denver ComSoc

  12. Quality Control Slide #10 Hubbell and Pevzner, 1999 Sengupta and Tompa, 2002 Colbourn et.al., 2002 Manufacture identical probes at several quality-control spots in order to test precision of production steps Quality control (fidelity) spots IEEE Denver ComSoc

  13. Relevant coding-theoretic ideas Slide #11 • Balanced code (Sengupta and Tompa, 2002): • An b×v binary matrix of zeros and ones with • each row has weight k; • each column has weight bounded between l and b-l, for some constant l; • any pair of columns is at least at Hamming distance d apart; • Superimposed designs in Renyi’s search model (Kautz and Singleton, 1964, Dyachkov and Rykov, 1983): • An b×v binary matrix of zeros and ones with • all Boolean sums composed of no more than s columns are distinct; • each row has weight exactly t; • Additional constraints: the Boolean sums form an error-correcting code with prescribed minimum distance d; IEEE Denver ComSoc

  14. Error-correcting microarray design Slide #12 • Probe multiplexing (Khan et.al, 2003) Probes Excluding hybridization effects, spot formation quality and under iid measurement noise, s p o t s X – vector of RNA levels corresponding to N genes Y – total concentration of RNA at all spots S - hybridization affinity matrix, T - spot quality matrix Decoding algorithm: numerical optimization IEEE Denver ComSoc

  15. VLSIPS/analysis for multiplexed arrays Slide #13 Features: • Multiple polymer synthesis at one given spot (for simplicity, will consider only two probes per spot) • Can use two different classes of linkers sensitive to different wavelengths so to select probes for extension (say, `blue’ and `green’ and`cyan’) Spots A T G C A T G C A T G C A T G C 1 2 3 4 5 6 gb gbccbcgg cb b g bg

  16. VLSIPS/analysis for multiplexed arrays Slide #14 Scheduling: shortest schedule of bases/colors (Using results from V. Dancık, Expected Length of Longest Common Subsequences, 1994) Set-up: two identical sets of M `blue’ and M `green’ randomly and uniformly chosen sequences of length N over the alphabet of size four Length of shortest schedule Chvatal-Sankoff constants Synchronous schedule, no `cyan’ colored steps: 8N IEEE Denver ComSoc

  17. Slide #15 Mask design: A C G T A C G T AT,CA AC,CC GT,GA TT,TA S1 S2 S3 S4 bgc cgc cb s1 s4 s3 s2 L(M)=4, L(M)=4, L(M)=2, L(M)=2, L(M)=2, L(M)=2, L(M)=2 s1 s3 s2 s4 L(M)=2, L(M)=2, L(M)=2, L(M)=2, L(M)=2, L(M)=2, L(M)=2

  18. Slide #16 Mask Design / Scheduling Neighborhood graph: complete graph with M vertices labeled by two distinct sequences No `cyan’ steps: weight of edge between two vertices sums of Hamming distances Issues: For reasons of controlled hybridization, different probes (blue and green) at the same spot should have fairly large Hamming distance (Milenkovic and Kashyap, 2005) Border-length minimization becomes less effective With cyan colored steps involved, the distance measure also depends on the longest common subsequence of the probes at the same spot IEEE Denver ComSoc

  19. Slide #17 Quality Control Coding Theorem: Assume that there exists a linear error-control code with parameters [n,k,d] containing the all-ones codeword. Then one can construct a quality control array for a multiplexed DNA chip with 2(2k-2) disjoint blue and green production steps and M probes such that the length of each quality control probe is 2(k-1)-1, and that the weights w of the columns in the quality control array satisfy Furthermore, with such an array any collection of less than n/(n-d) failed blue or green steps, respectively, can be uniquely identified. Open question: how does one extend this result for schedules involving `cyan’ colored production steps, and under `spot’ failures. IEEE Denver ComSoc

More Related