Illumina Methylation Library Construction – Oly/C.gigas Bisulfite-treated DNA

Took the bisulfite-treated DNA from 20151218 and made Illumina libraries using the TruSeq DNA Methylation Library Kit (Illumina).

Quantified the completed libraries using the Qubit 3.0 dsDNA BR Kit (ThermoFisher).

Evaluated the DNA with the Bioanalyzer 2100 (Agilent) using the DNA 12000 assay. Illumina recommended using the High Sensitivity assay, but we don’t have access to that so I figured I’d just give the DNA 12000 assay a go.

SampleName IndexNumber BarCode
1NF11 1 ATCACG
1NF15 2 CGATGT
1NF16 3 TTAGGC
1NF17 4 TGACCA
2NF5 5 ACAGTG
2NF6 6 GCCAAT
2NF7 7 CAGATC
2NF8 8 ACTTGA
M2 9 GATCAG
M3 10 TAGCTT
NF2_6 11 GGCTAC
NF_18 12 CTTGTA

 

Results:

Library Quantification (Google Sheet): 20151221_quantification_illumina_methylation_libraries

Test Name Concentration (ng/μL)
1NF11 Out of range
1NF15 2.14
1NF16 2.74
1NF17 2.64
2NF5 2.92
2NF6 Out of range
2NF7 2.42
2NF8 2.56
M2 Out of range
M3 2.1
NF2_6 2.38
NF2_18 Out of range

 

I used the Qubit’s BR (broad range) kit because I wasn’t sure what concentrations to expect. I need to use the high sensitivity kit to get a better evaluation of all the samples’ concentrations.

 

 

Bioanalyzer Data File (Bioanalyzer 2100): 2100_20expert_DNA_2012000_DE72902486_2015-12-21_16-58-43.xad

 

Ha! Well, looks like you definitely need to use the DNA High Sensitivty assay for the Bioanalyzer to pick up anything. Although, I guess you can see a slight hump in most of the samples at the appropriate sizes (~300bp); you just have to squint. ;)

Bioanalyzer – Bisulfite-treated Oly/C.gigas DNA

Following the guidelines of the TruSeq DNA Methylation Library Prep Guide (Illumina), I ran 1μL of each sample on an RNA Pico 6000 chip on the Seeb Lab’s Bioanalyzer 2100 (Agilent) to confirm that bisulfite conversion from earlier today worked.

Results:

Data File 1(Bioanlyzer 2100): 2100 expert_Eukaryote Total RNA Pico_DE72902486_2015-12-18_21-05-04.xad

Data File 1(Bioanlyzer 2100): 2100 expert_Eukaryote Total RNA Pico_DE72902486_2015-12-18_21-42-55.xad

 

 

Firstly, the ladder failed to produce any peaks. Not sure why this happened. Possibly not denatured? Seems unlikely, but next time I run the Pico assay, I’ll denature the ladder aliquot I use prior to running.

Overall, the samples look as they should (see image from TruSeq DNA Methylation Kit manual below), albeit some are a bit lumpy.

Bisulfite Treatment – Oly Reciprocal Transplant DNA & C.gigas Lotterhos DNA for BS-seq

After confirming that the DNA available for this project looked good, I performed bisulfite treatment on the following gDNA samples:

  • 1NF11
  • 1NF15
  • 1NF16
  • 1NF17
  • 2NF5
  • 2NF6
  • 2NF7
  • 2NF8
  • NF2_6
  • NF2_18
  • M2
  • M3

Sample names breakdown like this:

1NF#

1 = Fidalgo Bay outplants

NF = Fidalgo Bay broodstock origination

# = Sample number

2NF#

Same as above, but:

2 = Oyster Bay outplants

NF2_# (Oysters grown in Oyster Bay; DNA provided by Katherine Silliman)

NF2 = Fidalgo Bay broodstock origination, family #2

# = Sample number

M2/M3 = C.gigas from Katie Lotterhos

 

Followed the guidelines of the TruSeq DNA Methylation Library Prep Guide (Illumina).

Used the EZ DNA Methylation-Gold Kit (ZymoResearch) according to the manufacturer’s protocol with the following changes/notes:

  • Used 100ng DNA (per Illumina recs; Zymo recommends at least 200ng for “optimal results”).
  • Thermal cycling was performed in 0.5mL thin-wall tubes in a PTC-200 (MJ Research) using a heated lid
  • Centrifugations were performed at 13,000g
  • Desulphonation incubation for 20mins.

DNA quantity calculations are here (Google Sheet): 20151218_oly_bisulfite_calcs

Samples were stored @ -20C. Will check samples via Bioanalyzer before proceeding to library construction.

Agarose Gel – Oly gDNA for BS-seq Libraries, Take Two

The gel I ran earlier today looked real rough, due to the fact that I didn’t bother to equalize loading quantities of samples (I just loaded 1μL of all samples regardless of concentration). So, I’m repeating it using 100ng of DNA from all samples.

Additionally, this gel also includes C.gigas samples that Katie Lotterhos sent to us to see how they look.

Ran a 0.8% agarose, low-TAE gel, stained with ethidium bromide.

Results:

 

Look at that! The samples look MUCH nicer when they’re not overloaded and uniformly loaded!

Most have a prominent high molecular weight band (the band that’s closes to the top of the ladder, not the DNA visible in the wells). All exhibit smearing, but 2NF1 shows a weird accumulation of low molecular weight DNA.

Katie’s C.gigas samples (M1, M2, M3) look similar to the Olympia oyster gDNA, however her samples appear to have residual RNA in them (the fuzzy band ~500bp).

Will discuss with Steven which samples he wants to use for bisulfite treament and library construction.

Agarose Gel – Oly gDNA for BS-seq Libraries

Ran 1μL of each sample from yesterday’s DNA isolation on a 0.8% agarose, low-TAE gel, stained with ethidium bromide.

 

Results:

 

 

Since I didn’t load equal quantities of DNA, the intensities across the various samples is highly variable.

Those samples with high degree of smearing are also those with the highest concentrations. Thus, one would expect to be able to visualize a greater range of DNA sizes in a gel (because more DNA is present). Notice the samples with nice, high molecular weight bands and little smearing (1NF16, 1NF17). These are less than half the concentrations of all the samples that exhibit extensive smearing (2NF3, 2NF8, 1NF12). So, I think all samples will be fine for proceeding with bisulfite conversion and subsequent library construction.

However, I should re-run this gel using equalized DNA quantities for all samples…

 

DNA Isolation – Oly gDNA for BS-seq

Need DNA to prep our own libraries for bisulfite-treated high-throughput sequencing (BS-seq).

Isolated gDNA from the following tissue samples stored in RNAlater (tissue was not weighed) using DNAzol:

2NF1
2NF2
2NF3
2NF4
2NF5
2NF6
2NF7
2NF8
1NF11
1NF12
1NF13
1NF14
1NF15
1NF16
1NF17
1NF18

The sample coding breaks down as follows (see the project wiki for a full explanation):

2NF#

2 = Oysters outplanted in Fidalgo Bay

NF = Broodstock originated in Fidalgo Bay

# = Sample number

1NF#

1 = Oysters outplanted in Oyster Bay

NF = Broodstock originated in Fidalgo Bay

# = Sample number

 

DNA was isolated in the following manner:

  • Homogenized tissues in 500μL of DNAzol (Molecular Research Center; MRC).
  • Added additional 500μL of DNAzol.
  • Added 10μL of RNase A (10mg/mL, ThermoFisher); incubated 10mins @ RT.
  • Added 300μL of chloroform and mixed moderately fast by hand.
  • Incubated 5mins @ RT.
  • Centrifuged 12,000g, 10mins, RT.
  • Transferred aqueous phase to clean tube.
  • Added 500μL of 100% EtOH and mixed by inversion.
  • Pelleted DNA 5,000g, 5mins @ RT.
  • Performed 3 washes w/70% EtOH.
  • Dried pellet 3mins.
  • Resuspended in 100μL of Buffer EB (Qiagen).
  • Centrifuged 12,000g, 10mins, RT to pellet insoluble material.
  • Transferred supe to clean tube.

The samples were quantified using the Qubit dsDNA BR reagents (Invitrogen) according to the manufacturer’s protocol and used 1μL of sample for measurement.

Results:

Qubit data (Google Sheet): 20151216_Oly_gDNA_qubit_quants

SAMPLE CONCENTRATION (ng/μL)
2NF1 76.4
2NF2 175
2NF3 690
2NF4 11.7
2NF5 142
2NF6 244
2NF7 25
2NF8 456
1NF11 182
1NF12 432
1NF13 155
1NF14 21
1NF15 244
1NF16 112
1NF17 25.2
1NF18 278

 

Will run samples on gel tomorrow to evaluate gDNA integrity.

Sample Submission – 2bRAD Libraries for Genewiz

Pooled the libraries into a single sample for sequencing on Illumina HiSeq2500 by Genewiz.

Here’s the list of samples that were pooled:

SAMPLE NAME LIBRARY NAMES INDEX 1 (i7) LENGTH (bp) INDEX 2 (i5) LENGTH
SJW_Oly_2bRAD Oly RAD 02 CGTGAT 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 03 ACATCG 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 04 GCCTAA 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 06 TGGTCA 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 07 CACTGT 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 08 ATTGGC 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 14 GATCTG 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 17 TCAAGT 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 23 CTGATC 6 ATGCAT 6
SJW_Oly_2bRAD Oly RAD 30 AAGCTA 6 ATGCAT 6

Combined 40ng of all samples, except Oly RAD 30. Used only 20ng of Oly RAD 30 because it was the only sample that produced a single peak in qPCR melt curve analysis (i.e. no primer dimer). As such, it’s a rough assumption that the qPCR quantitation of all the other samples is twice as high as they should be due to the contribution of primer dimer amplification.

Calculations for pooling can be seen here (Google Sheet): 20151117_RAD_qPCR_data

The full list of samples (and the individual samples/libraries/indexes) submitted to Genewiz for this project by Katherine Silliman & me can be seen here (Google Sheet): White_BS1511196_R2_barcodes

 

Samples Received – C.gigas Tissue & DNA from Katie Lotterhos

Received 6 samples from Katie today. The box was labeled and stored @ -20C.

 

Here description of the samples, via email:

Lotterhos samples (gigas) arriving tomorrow

1) Mantle tissue samples of C. gigas were collected on 20140705 (source: Pipestem Inlet) by KEL
2) Extraction on 20141028 by VG using Qiagen DNAeasy Blood and Tissue Kit
3) Beadwash on 20150720 by VG using homemade sera-mag speed beads
4) Qubit 3.0 quantification on 20151206 by KEL and the following amounts were sent:

M1: 13 uL of 386 ug/mL
M2: 13.8 uL of 326 ug/mL
M3: 13.15 uL of 380 ug/mL (solution looked cloudy)

 

Sample Submission – Olympia oyster MBD-enriched DNA to ZymoResearch

We opted to go with ZymoResearch for this project because they could do the bisulfite treatment and finish the sequencing by the end of January.

Submitted the following 18 Ostrea lurida MBD-enriched gDNA samples to ZymoResearch for bisulfite treatment and subsequent Illumina sequencing (50bp, single read):

hc1_2B
hc1_4B
hc2_15B
hc2_17
hc3_1
hc3_10
hc3_11
hc3_5
hc3_7
hc3_9
ss2_14B
ss2_18B
ss2_9B
ss3_14B
ss3_15B
ss3_16B
ss3_20
ss3_3B
ss3_4B
ss5_18

The samples will be bisulfite treated, Illumina libraries constructed, multiplexed, and this multiplexed library will be sequenced three times.