Category Archives: LSU C.virginica Oil Spill MBD BS Sequencing

Sequencing Data – LSU C.virginica MBD BS-Seq

Our sequencing data (Illumina HiSeq2500, 100SE) for this project has completed by Univ. of Oregon Genomics Core Facility (order number 2112).

Samples sequenced/pooled for this run:

Sample Treatment Barcode
HB2 25,000ppm oil ATCACG
HB16 25,000ppm oil TTAGGC
HB30 25,000ppm oil TGACCA
NB3 No oil ACAGTG
NB6 No oil GCCAAT
NB11 No oil CAGATC

All code listed below was run on OS X 10.9.5

Downloaded all 15 fastq.gz files to Owl/web/nightingales/C_virginica:

$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_001.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_002.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_003.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_004.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_005.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_006.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_007.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_008.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_009.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_010.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_011.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_012.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_013.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_014.fastq.gz
$curl -O http://gcf.uoregon.edu:8080/job/download/2112?fileName=lane1_NoIndex_L001_R1_015.fastq.gz

 

Renamed all files by removing the beginning of each file name (2112?fileName=) and replacing that with 2112_:

$for file in 2112*lane1_NoIndex_L001_R1_0*; do mv "$file" "${file/#2112?fileName=/2112_}"; done

 

Created a directory readme.md (markdown) file to list & describe directory contents: readme.md

$ls *.gz >> readme.md

Note: In order for the readme file to appear in the web directory listing, the file cannot be all upper-case.

 

Created MD5 checksums for each fastq.gz file: checksums.md5

$md5 *.gz >> checksums.md5

Library Cleanup – LSU C.virginica MBD BS Library

I was contacted by the sequencing facility at the University of Oregon regarding a sample quality issue with our library.  As evidenced by the electropherogram below, there is a great deal of adaptor primer dimer (the peak at 128bp):

 

This is a problem because such a high quantity of adaptor sequence will result in the majority of reads coming off the Illumina being just adaptor sequences.

With the remainder of the library sample prepared earlier, I performed the recommended clean up procedure for removing adaptor sequences in the EpiNext Post-Bisulfite DNA Library Preparation Kit – Illumina (Epigentek).    Briefly:

  • Brought sample volume up to 20uL with NanoPure H2O (added 9.99uL)

  • Added equal volume of MQ Beads

  • Washed beads 3x w/80% EtOH

  • Eluted DNA w/12uL Buffer EB (Qiagen)

After clean up, quantified the sample via fluorescence using the Quant-iT DNA BR Kit (Life Technologies/Invitrogen).  Used 1uL of the sample and the standards.  All standards were run in duplicate and read on a FLx800 plate reader (BioTek).

Results are here: 20150122 – LSU_virginicaMBDlibraryCleanup

Library concentration = 2.46ng/uL

Brought the entire sample up to 20uL with Buffer EB (Qiagen) and a final concentration of 0.1% Tween-20 (required by the sequencing facility).

Sent sample to the University of Oregon to replace our previous submission.

Bisulfite NGS Library – LSU C.virginica Oil Spill MBD Bisulfite DNA Sequencing Submission

Combined the following libraries in equal quantities (17ng each) to create a single, multiplexed sample for sequencing (LSU_Oil_01):

  • HB2 – 1 (ATCACG)
  • HB16 – 3 (TTAGGC)
  • HB30 – 4 (TGACCA)
  • NB3 – 5 (ACAGTG)
  • NB6 – 6 (GCCAAT)
  • NB11 – 7 (CAGATC)

Quantified pooled libraries using the Quant-iT dsDNA BR Kit (Invitrogen) with a FLx800 plate reader (BioTek). Used 1μL of the pooled sample, run in duplicate. Used 1uL of standards, run in duplicate.

Results:

pooled libraries = 6.575ng/μL

Will submit to University of Oregon Genomics Core Facility for 100bp, single end Illumina HiSeq2500 sequencing. They need 10nM of sample. For a library with average size range of 300-400bp, this requires a sample volume of 20uL with a concentration of 2.28ng/μL in a solution of 0.1% Tween20 in Buffer EB (Qiagen).

Combined 6.94μL of pooled libraries with 13.06 of 0.1% Tween20/EB solution.

Submitted sample LSU_Oil_01 to University of Oregon Genomics Core Facility via O/N FedEx on dry ice. Sample was assigned order # 2112.

Bisulfite NGS Library Prep – LSU C.virginica Oil Spill MBD Bisulfite DNA and Emma’s C.gigas Larvae OA Bisulfite DNA (continued from yesterday)

Continued library prep from yesterday. Set up Library Amplification according to the protocol. The samples received the following Barcode Indices:

  • HB2 – 1 (ATCACG)
  • HB5 – 2 (CGATGT)
  • HB16 – 3 (TTAGGC)
  • HB30 – 4 (TGACCA)
  • NB3 – 5 (ACAGTG)
  • NB6 – 6 (GCCAAT)
  • NB11 – 7 (CAGATC)
  • NB21 – 12 (CTTGTA)
  • 1A1 – 2 (CGATGT)
  • 1A2 – 1 (ATCACG)
  • 6A1 – 4 (TGACCA)
  • 6A2 – 5 (ACAGTG)
  • 103B1 – 6 (GCCAAT)
  • 103B2 – 7 (CAGATC)
  • 105A4 – 12 (CTTGTA)
  • 105A5 – 11 (GGCTAC)

Due to differences in input DNA quantities, samples were run with different numbers of thermal cycles.

13 thermal cycles were run for the following samples:

  • 1A1
  • 105A4
  • 105A5

22 thermal cycles were run for the following samples:

  • HB2
  • HB5
  • HB16
  • HB30
  • NB3
  • NB6
  • NB11
  • NB21
  • 1A2
  • 6A1
  • 6A2
  • 103B1
  • 103B2

Samples were quantified with 1uL of each sample using the Quant-iT dsDNA BR Kit (Invitrogen). Used 5uL of each standard and standards were run in duplicate.

Results:

Bisulfite NGS Library Prep – LSU C.virginica Oil Spill Bisulfite DNA and Emma’s C.gigas Larvae OA Bisulfite DNA

Constructed next generation libraries (Illumina) using the bisulfite-treated DNA from yesterday using the EpiNext Post-Bisulfite DNA Library Preparation Kit – Illumina (Epigentek). Samples were processed according to the manufacturer’s protocol up to Section 8 (Library Amplification) with the following changes:

– Skipped Section 7.1 (recommended to do so in the protocol due to low quantity of input DNA)

Samples were stored O/N @ -20C.

dA Tailing Master Mix

10x Tailing Buffer 1.5uL x 17.6 = 26.4uL

Klenow 1uL x 17.6 = 17.6uL

H2O 0.5uL x 17.6 = 8.8uL

Add 3uL of master mix to each sample

Adaptor Ligation

2x Ligation Buffer 17uL x 17.6 – 299.2uL

T4 DNA Ligase 1uL x 17.6uL = 17.6uL

Adaptors 1uL x 17.6 = 17.6uL

Added 19uL of master mix to each sample

dsDNA Conversion Master Mix

5x Conversion Buffer 4uL x 17.6 = 70.4uL

C.P. 2uL x 17.6 = 35.2uL

H2O 3uL x 18.6 = 52.8uL

Add 9uL of master mix to each sample

End Repair

10x Buffer 2uL x 17.6 = 35.2uL

Enzyme 1uL x 17.6 = 17.6uL

H2O 5uL x 17.6 = 88uL

Added 8uL of master mix to each sample

Bisulfite Conversion – LSU C.virginica Oil Spill MBD DNA and Emma’s C.gigas Larvae OA DNA

Performed bisulfite conversion on MBD DNA samples from LSU C.virginica oil spill samples (see 201411202 and 20141126) and Emma’s C.gigas larvae OA DNA samples (see 20141121) with the Methylamp DNA Modification Kit (Epigentek).

Added 4uL of H2O to each of Emma’s DNA samples to bring them up to 24uL.

Samples were processed according to the manufacturer’s protocol.

Samples were eluted with 10uL of Solution R6 and stored @ -20C.

EtOH Precipitation – LSU C.virginica Oil Spill MBD Continued (from 20141126)

Precipitation was continued according to the MethylMiner Methylated DNA Enrichment Kit (Invitrogen). Since I will need sample volumes of 24uL for the subsequent bisulfite conversion, I resuspended the samples in 29uL of water (will use 2.5uL x 2 reps for quantification).

Samples to be quantified:

NC = non-captured (i.e. non-methylated)

E = eluted (i.e. methylated)

  • HB2 NC
  • HB5 NC
  • HB16 NC
  • HB30 NC
  • NB3 NC
  • NB6 NC
  • NB11 NC
  • NB21 NC
  • HB2 E
  • HB5 E
  • HB16 E
  • HB30 E
  • NB3 E
  • NB6 E
  • NB11 E
  • NB21 E
  • Control NC
  • Control E

Samples were quantified using the Quant-IT BS Kit (Invitrogen) with a plate reader (BioTek). All samples were run in duplicate. Used 2.5uL of each sample for quantification.

Samples were stored in @ -20C (FTR 209) in the bisulfite seq box created by Claire for this project.

Results:

20141202_LSU_Virginica_MBD:

https://docs.google.com/spreadsheets/d/1NrrVmYsUQcstnrt4583mYN2PeVav54luyFvVUEkcjWE/edit?usp=sharing

Methylated DNA Enrichment (MBD) – LSU C.virginica Oil Spill gDNA

Enrichment was performed using the MethylMiner Methylated DNA Enrichment Kit (Invitrogen) according to the manufacturer’s protocol with the following changes:

– Used 25uL of Dynabeads M-280 (10uL/ug of input DNA) and 15uL of MBD-Biotin Protein (7uL/ug of input DNA).

– Followed the corresponding instructions for the volumes listed above and for quantities of input DNA > 1ug – 10ug

– A single elution with 2000mM NaCl was performed

– EtOH precipitation: Samples were incubated over the long weekend at -80C.

Gel – Sheared LSU C.virginica Oil Spill gDNA (from yesterday)

Ran ~250ng of sheared C.virginica gDNA from yesterday’s shearing.

Results:

Ladder used: O’GeneRuler 100bp Ladder (ThermoFisher)

The shearing is, surprisingly, very inconsistent across the samples. The target average fragment size was ~350bp. However, most of these samples are <250bp. The MethylMiner Kit (Invitrogen) suggests that an average fragment length of 100 – 200bp is ideal for short-read high-throughput sequencing, but we’re going to perform a bisulfite conversion on these which will result in some additional fragmentation, further reducing the average fragment size. Will proceed with methylated DNA enrichment.

DNA Shearing – LSU C.virginica Oil Spill gDNA

Used the remainder of the “sheared” samples (see today’s earlier entry; ~2750ng). Brought the volumes up to 80uL and transferred to 0.5mL snap cap tubes. The volume of 80uL was selected because it’s above the minimum volume required for shearing in 0.5mL tubes (10uL according to the Biorupter 300 manual) and the MethylMiner Kit (Invitrogen) requires the input DNA volume to be <= 80uL.

DNA was sheared with the following parameters:

Low power

30 cycles of:

30s on

30s off

Target average fragment size is ~350bp.

See tomorrow (20141126) for the gel.