Chris:LabNotes/sci-Methyl Seq/Calendar/2017/2017-7-3: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Cjwei
>Cjwei
Line 169: Line 169:
  <u>GGTTGGTTGTGGTGGGTGATTGTGGTAGGGAAG      CTTCCCTACCACAATCACCCACCACAACCAACC      0.55    65.15    66        Pass</u> (Adpt1_v5)
  <u>GGTTGGTTGTGGTGGGTGATTGTGGTAGGGAAG      CTTCCCTACCACAATCACCCACCACAACCAACC      0.55    65.15    66        Pass</u> (Adpt1_v5)
  <u>TGGTGGGATGGGTGAAGGTTGTTGGGTGTGATT      AATCACACCCAACAACCTTCACCCATCCCACCA      0.52    65.90    67        Pass</u> (Adpt1_v5)
  <u>TGGTGGGATGGGTGAAGGTTGTTGGGTGTGATT      AATCACACCCAACAACCTTCACCCATCCCACCA      0.52    65.90    67        Pass</u> (Adpt1_v5)
  GATGTGGTGGTGAAGGTGATGGGTAGTGAGGGA      TCCCTCACTACCCATCACCTTCACCACCACATC      0.55    65.07    66.1      Pass
  <u>GATGTGGTGGTGAAGGTGATGGGTAGTGAGGGA      TCCCTCACTACCCATCACCTTCACCACCACATC      0.55    65.07    66.1      Pass</u> (Adpt2_v5 Barcode Seq)
  TGGGAAGGGAAGGTGGAGTGGAAGTTGGAGAGT      ACTCTCCAACTTCCACTCCACCTTCCCTTCCCA      0.55    65.94    67.2      Pass
  <u>TGGGAAGGGAAGGTGGAGTGGAAGTTGGAGAGT      ACTCTCCAACTTCCACTCCACCTTCCCTTCCCA      0.55    65.94    67.2      Pass</u> (Adpt2_v5 Barcode Seq)
  ATGTTGTGTGGAGGGTGTGGGTGGGATGAGAAA      TTTCTCATCCCACCCACACCCTCCACACAACAT      0.52    65.54    66.6      Pass
  <u>ATGTTGTGTGGAGGGTGTGGGTGGGATGAGAAA      TTTCTCATCCCACCCACACCCTCCACACAACAT      0.52    65.54    66.6      Pass</u> (Adpt2_v5 Barcode Seq)
  GTAGGTGTTGGGTTTGGAAGGTGAGGTGGTGGT      ACCACCACCTCACCTTCCAAACCCAACACCTAC      0.55    65.85    66.9      Pass
  GTAGGTGTTGGGTTTGGAAGGTGAGGTGGTGGT      ACCACCACCTCACCTTCCAAACCCAACACCTAC      0.55    65.85    66.9      Pass
  ATGGGAATGGGAGGGTGGGTGAGGTGGT<u>AATGA</u>      TCATTACCACCTCACCCACCCTCCCATTCCCAT      0.55    66.11    67.5      Pass
  ATGGGAATGGGAGGGTGGGTGAGGTGGT<u>AATGA</u>      TCATTACCACCTCACCCACCCTCCCATTCCCAT      0.55    66.11    67.5      Pass
<u>ATGAATGTGTGTGGGAGGATGGTGTGGGAAGGG      CCCTTCCCACACCATCCTCCCACACACATTCAT      0.55    65.83    66.9      Pass</u> (Adpt2_v5 Universal Primer)
*Below is the output of the script for ACGT fragment in Adpt1 (will be digested away using DraIII).  For this filler sequence, we want to decrease the melting temperature and length of the oligo to a normal PCR primer as follows (based on IDT primer design <https://www.idtdna.com/pages/decoded/decoded-articles/pipet-tips/decoded/2013/10/21/designing-pcr-primers-and-probes>:
*Below is the output of the script for ACGT fragment in Adpt1 (will be digested away using DraIII).  For this filler sequence, we want to decrease the melting temperature and length of the oligo to a normal PCR primer as follows (based on IDT primer design <https://www.idtdna.com/pages/decoded/decoded-articles/pipet-tips/decoded/2013/10/21/designing-pcr-primers-and-probes>:
  Tm = 60-64C
  Tm = 60-64C

Revision as of 21:41, 5 July 2017

sci-Methyl Seq Barcode 1 Design v5; Barcode 2 Design v4

Background

  • We want to make to major changes to the adapter designs:
    • Adpt2 Y-Adapter Design:
      • We want to redesign the adapter sequences in order to avoid having to use a single primer PCR reaction, which seems to be giving us some trouble. In addition, we would be unable to perform PCR after bisulfite conversion to add appropriate sequencing adapters since the ends of the fragments with Adp2 would contain the same sequence. Instead, we would need to use commercial kits to perform library construction (such as Accel-NGS Methyl-Seq <https://swiftbiosci.com/products/accel-ngs-methyl-seq-dna-library-kit/>). This could prove problematic as further optimization/cost would be added in order to use these commercial kits.
      • Instead, we want to redesign Adpt2 to be Y-adapters with unique sequences on both strands that serve as primer binding sites. The general design would be as follows:
                                      3'
                                     /
                                    /
                                   /
         5'  -----[Barcode2]UMI----
                           3' <----
                                   \
                                    \
                                     \
                                      5'
    • Use 3-base sticky end for Adpt1/Adpt2 ligation:
      • Previously, I've tried a CC/GG sticky end, which seems to still allow Adpt2 to be annealing directly onto the original template fragment non-specifically. Consequently, instead, we will try using a three-base sticky end to help with increasing specificity of Adpt1/Adpt2 ligation.
      • To achieve this, we will use a restriction enzyme to cut the ends of Adpt1/Adpt2 in order to create the 3-base sticky end (use DraIII-HF from NEB).
Note: NoG sequences in red/H, NoC sequences in blue/D
HpyCH4III: ACN|GT
           TG|NCA
DraIII: CACNNN|GTG
        GTG|NNNCAC
Adpt1: (cut by HpyCH4III and DraIII)
    5' /5Phos/GT-----TTHH[Barcode1]-----CACNNN 3'
    3'       TCA-----AADD[Barcode1]-----GTG    5'
Adpt2: (cut by DraIII)
                                                3'
                                               /
                                              /
                                             /
    5' /5Phos/GTG-----[Barcode2]DDDDDDDD-----
    3'     NNNCAC-----[Barcode2]HHHHHHHH-----
                                             \
                                              \
                                               \
                                                5'

Overview

End Repair/dA-Tailing
     5'  -----A 3'
     3' A-----  5'

           |
           V
Adpt1_v5: (second strand synthesis; cut by HpyCH4III and DraIII)
    5' -----ACA|GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
                                 3' <-----GTG|NNNCAC----- 5'
                          | Second strand synthesis
                          V
    5' -----ACA|GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
    3' -----TG|TCA-----AADD[Barcode1]-----GTG|NNNCAC----- 5'
                          | Cut by HpyCH4III
                          V
    5' /5Phos/GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
    3'       TCA-----AADD[Barcode1]-----GTG|NNNCAC----- 5'
                          | Cut by DraIII
                          V
    5' /5Phos/GT-----TTHH[Barcode1]-----CACNNN 3'
    3'       TCA-----AADD[Barcode1]-----GTG/5Phos/    5'
Ligate Adpt1_v5 (using same optimized ligation protocol as before)
     5'    GTG-----[Barcode1]DDAA-----ACT|-----A|GT-----TTHH[Barcode1]-----CACNNN 3'
     3' NNNCAC-----[Barcode1]HHTT-----TG|A-----|TCA-----AADD[Barcode1]-----GTG    5'
           |
           V
Adpt2_v4: (second strand synthesis; cut by DraIII)
                                                     3'
                                                    /
                                                   /
                                                  /
    5' -----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
    3'                                       -----
                                                  \
                                                   \
                                                    \
                                                     5'
                          | Second strand synthesis
                          V
                                                     3'
                                                    /
                                                   /
                                                  /
    5' -----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
    3' -----GTG|NNNCAC-----[Barcode2]HHHHHHHH-----
                                                  \
                                                   \
                                                    \
                                                     5'
                          | Cut by DraIII
                          V
                                                3'
                                               /
                                              /
                                             /
    5' /5Phos/GTG-----[Barcode2]DDDDDDDD-----
    3'     NNNCAC-----[Barcode2]HHHHHHHH-----
                                             \
                                              \
                                               \
                                                5'
Ligate Adpt2_v4 (using same optimized ligation protocol as before)
    5'                                                                                                                                              3'
     \                                                                                                                                             /
      \                                                                                                                                           /
       \                                                                                                                                         /
        -----HHHHHHHH[Barcode2]-----CACNNN|GTG-----[Barcode1]DDAA-----ACT|-----A|GT-----TTHH[Barcode1]-----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
        -----DDDDDDDD[Barcode2]-----GTG|NNNCAC-----[Barcode1]HHTT-----TG|A-----|TCA-----AADD[Barcode1]-----GTG|NNNCAC-----[Barcode2]HHHHHHHH-----
       /                                                                                                                                         \
      /                                                                                                                                           \
     /                                                                                                                                             \
    3'                                                                                                                                              5'
           |
           V

Bisulfite conversion (C->U)

           |
           V

PCR (P7 will be added first followed by P5)
 P5: 5' AATGATACGGCGACCACCGA 3'
 P7: 5' CAAGCAGAAGACGGCATACGAGAT 3'
    5'                                                                                                                                              3'
    P5                                                                                                                                              P7
     \                                                                                                                                             /
      \                                                                                                                                           /
       \                                                                                                                                         /
        -----HHHHHHHH[Barcode2]-----CACNNN|GTG-----[Barcode1]DDAA-----ACT|-----A|GT-----TTHH[Barcode1]-----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
        -----DDDDDDDD[Barcode2]-----GTG|NNNCAC-----[Barcode1]HHTT-----TG|A-----|TCA-----AADD[Barcode1]-----GTG|NNNCAC-----[Barcode2]HHHHHHHH-----
       /                                                                                                                                         \
      /                                                                                                                                           \
     /                                                                                                                                             \
    P7                                                                                                                                              P5
    3'                                                                                                                                              5'
           |
           V

Sequencing (based on HiSeq 2500 Dual-Indexed on Single-Read Flow Cell File:Cw 20170705 Indexed Sequencing Overview.pdf.gz)
    P5                                                                                                                                              P7
     \                                                                                                                                             /
      \                                                                                                                                           /
       \                                                                                                                                         /
        -----HHHHHHHH[Barcode2]-----CACNNN|GTG-----[Barcode1]DDAA-----ACT|-----A|GT-----TTHH[Barcode1]-----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
                                              ----->             ----->                                              ----->
                                              Index2 (i5)        Read1                                               Index1 (i7)

Potential Filler Sequences

  • We want to redesign the potential filler sequences to be used for the adapters. The purpose of the filler sequences are two fold: (1) they will be used as sequencing primer binding sites on the Illumina flowcell, and (2) the filler sequences at the end of Adpt2 will be used as primer sites for adding P5/P7 sequences.
  • In order to satisfy the first purpose of the filler sequences, we want to ensure that they exhibit similar properties as standard Illumina sequencing primer sites. This means that they must closely match the following parameters. This is based on <http://nextgen.mgh.harvard.edu/CustomPrimer.html> and File:Cw 20170705 Illumina Adapter Sequences.pdf
Tm = 66C (recommended to use IDT oligoAnalyzer to calculate Tm using Neighbor Joining Method)
Length = 33bp
GC = 52%
  • We can also use truncated versions of these filler sequences as primer binding sites for adding the P5/P7 sequences during PCR (see the above Illumina Adapter Sequences document under Nextera adapters)
  • In order to design these adapters, I'll be using a new custom script to generate random oligos (File:Cw 20170705 primerGenerator.v2.txt). The script will generate random sequences using the given parameters:
Length = 33bp
Min Tm = 65C
Max Tm = 70C
GC content = 45-55%
No single base runs of C's or G's >3
No more than 3 C's or G's at either ends of the filler sequence
  • Below is the output of the script using the above parameters (noC):
Seq                                     revComp(Seq)                            GC      Tm       IDT Tm    Primer Stats Notes (http://www.bioinformatics.org/sms2/pcr_primer_stats.html)
GGTTGGTTGTGGTGGGTGATTGTGGTAGGGAAG       CTTCCCTACCACAATCACCCACCACAACCAACC       0.55    65.15    66        Pass (Adpt1_v5)
TGGTGGGATGGGTGAAGGTTGTTGGGTGTGATT       AATCACACCCAACAACCTTCACCCATCCCACCA       0.52    65.90    67        Pass (Adpt1_v5)
GATGTGGTGGTGAAGGTGATGGGTAGTGAGGGA       TCCCTCACTACCCATCACCTTCACCACCACATC       0.55    65.07    66.1      Pass (Adpt2_v5 Barcode Seq)
TGGGAAGGGAAGGTGGAGTGGAAGTTGGAGAGT       ACTCTCCAACTTCCACTCCACCTTCCCTTCCCA       0.55    65.94    67.2      Pass (Adpt2_v5 Barcode Seq)
ATGTTGTGTGGAGGGTGTGGGTGGGATGAGAAA       TTTCTCATCCCACCCACACCCTCCACACAACAT       0.52    65.54    66.6      Pass (Adpt2_v5 Barcode Seq)
GTAGGTGTTGGGTTTGGAAGGTGAGGTGGTGGT       ACCACCACCTCACCTTCCAAACCCAACACCTAC       0.55    65.85    66.9      Pass
ATGGGAATGGGAGGGTGGGTGAGGTGGTAATGA       TCATTACCACCTCACCCACCCTCCCATTCCCAT       0.55    66.11    67.5      Pass
ATGAATGTGTGTGGGAGGATGGTGTGGGAAGGG       CCCTTCCCACACCATCCTCCCACACACATTCAT       0.55    65.83    66.9      Pass (Adpt2_v5 Universal Primer)
Tm = 60-64C
Length = 22
Seq                     revComp(Seq)            GC      Tm       IDT Tm    Primer Stats Notes (http://www.bioinformatics.org/sms2/pcr_primer_stats.html)
TTACTGCCGTTGCGAACGCGAA  TTCGCGTTCGCAACGGCAGTAA  0.55    60.24    62        Pass (Adpt1_v5)
TGTGCGTTTGTGGCGCAACAAC  GTTGTTGCGCCACAAACGCACA  0.55    60.12    61.9      More than 3 self-annealing/hairpin bases
TGACCACGTGCGTTCGTTGCAA  TTGCAACGAACGCACGTGGTCA  0.55    60.36    62.3      More than 3 self-annealing bases
TGTGCGATGCGGCAATGACCAA  TTGGTCATTGCCGCATCGCACA  0.55    60.24    62.5      More than 50% of the bases are self annealing
TGCTATGCTGCGCGGAATGCAA  TTGCATTCCGCGCAGCATAGCA  0.55    60.37    62.5      More than 3 self-annealing bases

Barcode 1 Adapter Design (Adpt1_v5)

  • As described above, the following is how we would eventually create dsDNA Adpt1_v5:
Adpt1_v5: (second strand synthesis; cut by HpyCH4III and DraIII)
    5' -----ACA|GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
                                 3' <-----GTG|NNNCAC----- 5'
                          | Second strand synthesis
                          V
    5' -----ACA|GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
    3' -----TG|TCA-----AADD[Barcode1]-----GTG|NNNCAC----- 5'
                          | Cut by HpyCH4III
                          V
    5' /5Phos/GT-----TTHH[Barcode1]-----CACNNN|GTG----- 3'
    3'       TCA-----AADD[Barcode1]-----GTG|NNNCAC----- 5'
                          | Cut by DraIII
                          V
    5' /5Phos/GT-----TTHH[Barcode1]-----CACNNN 3'
    3'       TCA-----AADD[Barcode1]-----GTG    5'
  • Today, I'll be ordering two different oligos to perform second strand synthesis and restriction enzyme digestion: (using Barcode1 = ACCTC, HH = AC)
Adpt1_v5 Barcode Seq (This is top strand containing Barcode1 Sequence)
5' AATGAACAGTCCAACCAACACCACCCACTAACACCATCCCTTCTTAC[ACCTC]ACCACCCTACCCACTTCCAACAACCCACACTAACACAATGTGTTCGCGTTCGCAACGGCAGTAA 3'
Adpt1_v5 Universal Primer (Primer to start second strand synthesis)
5' TTACTGCCGTTGCGAACGCGAA 3'

Barcode 2 Adapter Design (Adpt2_v4)

  • As described above, the following is how we would eventually create dsDNA Adpt2_v4:
Adpt2_v4: (second strand synthesis; cut by DraIII)
                                                     3'
                                                    /
                                                   /
                                                  /
    5' -----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
    3'                                       -----
                                                  \
                                                   \
                                                    \
                                                     5'
                          | Second strand synthesis
                          V
                                                     3'
                                                    /
                                                   /
                                                  /
    5' -----CACNNN|GTG-----[Barcode2]DDDDDDDD-----
    3' -----GTG|NNNCAC-----[Barcode2]HHHHHHHH-----
                                                  \
                                                   \
                                                    \
                                                     5'
                          | Cut by DraIII
                          V
                                                3'
                                               /
                                              /
                                             /
    5' /5Phos/GTG-----[Barcode2]DDDDDDDD-----
    3'     NNNCAC-----[Barcode2]HHHHHHHH-----
                                             \
                                              \
                                               \
                                                5'
  • Below are the two oligos I'll be ordering today to perform second strand synthesis and restriction enzyme digestion: (using Barcode2 = GATAGG, DDDDDDDD = TGGAGAGG, NNN=CTT [such that won't anneal to the A-tailed sticky end of the original template with GAA in final sticky end])
Adpt2_v4 Barcode Seq (This is top strand containing Barcode2 Sequence)
5' AATGACACCTTGTGGATGTGGTGGTGAAGGTGATGGGTAGTGAGGGA[GATAGG]TGGAGAGGTGGGAAGGGAAGGTGGAGTGGAAGTTGGAGAGTATGTTGTGTGGAGGGTGTGGGTGGGATGAGAAA 3' (italicized portion is part of Y-tail)
Adpt2_v4 Universal Primer (Primer to start second strand synthesis and complete Y-adapter formation)
5' GTAGGTGTTGGGTTTGGAAGGTGAGGTGGTGGTACTCTCCAACTTCCACTCCACC 3' (italicized portion is part of Y-tail)