Construction of individuals ddRADseq libraries for macro-algae (Kelp) V.3

Stéphane Mauger, Komlan Avia

Published: 2024-07-30 DOI: 10.17504/protocols.io.rm7vzjo68lx1/v3

Abstract

This protocol describes a double digested restriction-site associated DNA (ddRADseq) procedure, that is a variation on the original RAD sequencing method (Davey & Blaxter 2011), which is used for de novo SNP discovery and genotyping.

This protocol differs from the original ddRADseq protocol (Peterson et al 2012), in which the samples are pooled just after the ligation to adaptors (i.e. before size selection and PCR). This protocol is an update of the protocol from Claire Daguin Thiebaut et al. (dx.doi.org/10.17504/protocols.io.bv4tn8wn) adapted for macro-algae.

The following protocol is intended for the construction of individual ddRADseq libraries from genomic DNA of various macro-algae samples (Kelp). In the present protocol, we added a genomic DNA purification step to eliminate the inhibitors of PCR and Ligation present in macro-algae (polysaccharides). Moreover, all samples are treated separately until final PCR amplification (Rad Taq enrichment step) performed before pooling.

Despite being slightly more costly and time-consuming in the lab, it allows for fine adjustement of each sample representation in the final library pool ensuring similar number of reads between samples. Finally, we have defined new P1 adapters (barcodes) with variable sequences and variable sizes (6bp to 13bp) to increase the efficiency of the Illumina sequencing.

Briefly, purified genomic DNA from the samples are individually digested with 2 restriction enzymes PstI/HhaI PstI/HhaI or PstI/MesI PstI/MesI (one rare-cutter and one more frequent cutter) then ligated to a barcoded adaptor (among 24 available) at one side, and a single adaptor at the other side, purified with magnetic beads, and PCR-amplified allowing the addition of a Illumina index (among 24 available) for multiplexing a maximum of 576 samples per library . Samples are then pooled in equimolar conditions after visualisation on an agarose gel. Purification and size selection is then performed before final quality control of the library and sequencing.

This protocol has proven its effectiveness in several genetic studies of marco-algae populations.

Overview of ddRADseq libraries preparation
Overview of ddRADseq libraries preparation
Diagram of oligos and adaptaters; final library; BRIEF GLOSSARY
Diagram of oligos and adaptaters; final library; BRIEF GLOSSARY
Final ddRAD library construction
Final ddRAD library construction

Before start

  1. Prepare all buffers and solutions in advance (see Step 1 to Step 5)

  2. If not using Retschcopyright Mixer Mill MM 301 (or equivalente) and Grinding ball for the sample grinding, you can use Lysing Matrix H tube with FastPrep-24TM Classic or manual grinding as a last resort.

Steps

Solutions and buffers preparations

1.

5 M sodium chloride solution (NaCl)

29.2g

Dissolve the slat in MilliQ water and fill up to 100mL.

Autoclave.

Store at Room temperature

2.

Annealing buffer stock (10x)

Annealing buffer composed 100 mM Tris-HCl, pH8; 500 mM NaCl and 10 mM EDTA

5mL

5mL

1mL

39mL

Homogenize and autoclave.

Store at Room temperature

Preparation of double-stranded barcoded P1 adaptors 4µM

3.

Single-stranded oligos NGS grade P1 need to be annealed with their appropriate partner before ligation. We provide sequences for 48 uniquely barcoded adapter P1 oligo pairs (oligos P1_PstI_x.F and P1_PstI_x.R), see the Barcoded_P1_adaptors.xlsx file below .

To create Adapter P1, combine each oligo Forward with its complementary oligo Reverse in a 1:1 ratio in working strength annealing buffer (final buffer concentration 1x) for a total annealed adapter

concentration of 4µM.

In house barcoded P1 adaptors sequences (NGS grade needed) :

Barcoded_P1_adaptors.xlsx

3.1.

In a PCR plate wells, combine each oligo P1_PstI_x.F with its complementary oligo P1_PstI_x.R :

4µL

4µL

10µL

82µL

3.2.
ABCDEFGH
P1_PstI_01P1_PstI_09P1_PstI_17
P1_PstI_02P1_PstI_10P1_PstI_18
P1_PstI_03P1_PstI_11P1_PstI_19
P1_PstI_04P1_PstI_12P1_PstI_20
P1_PstI_05P1_PstI_13P1_PstI_21
P1_PstI_06P1_PstI_14P1_PstI_22
P1_PstI_07P1_PstI_15P1_PstI_23
P1_PstI_08P1_PstI_16P1_PstI_24

Example of a plate map for barcoded P1 adaptors. Allow enough space betweeen the rows to avoid cross-contaminations between barcodes.

3.3.

The reaction is performed in a thermocycler with the following PCR cycling conditions :

ABCD
Cycle stepTemperatureTimeCycles
Initial Denaturation97.5°C2.5 min1
Annealing96°C (-3°C per cycle)1 min25
Hold4°C

PCR cycling conditions

Store at4°C (or at-20°C for a long-term storage)

Preparation of double-stranded P2 adaptors 40µM

4.

Single-stranded oligos NGS grade P2 need to be annealed with their appropriate partner before PCR. We provide sequences for 4 uniquely adapter P2 oligo pairs (oligos P2_HhaI.F and P2_HhaI.R or P2_MseI.F and P2_MseI.R), see the No-Barcoded_P2_adaptors.xlsx file below .

To create Adapter P2, combine each oligo Forward with its complementary oligo Reverse in a 1:1 ratio in working strength annealing buffer (final buffer concentration 1x) for a total annealed adapter

concentration of 40µM.

No-barcoded P2 adaptors sequences (NGS grade needed) :

No-Barcoded_P2_adaptors.xlsx

4.1.

In 1.5mL microtube, combine oligo P2_HhaI.F with its complementary oligo P2_HhaI.F (or P2_MseI.F and P2_MseI.R)

400µL

400µL

100µL

100µL

Then aliquot this volume into 125µL in each well of a 8- PCR tube strip.

4.2.

The reaction is performed in a thermocycler with the following PCR cycling conditions :

ABCD
Cycle stepTemperatureTimeCycles
Initial Denaturation97.5°C2.5 min1
Annealing96°C (-3°C per cycle)1 min25
Hold4°C

PCR cycling conditions

Pool all reaction in a1.5mL tube.

Store at 4°C (or at -20°C for a long-term storage).

Preparation of Illumina indexed primers mix (5µM)

5.

In 24 1.5mL microtubes, combine each of the 24 Illumina indexed reverse primers ILLPCR2_ind01 to ILLPCR2_ind27 (no primer numbers ind17, ind24 and ind26) with the Illumina no-indexed forward primer ILLPRC1, see the Illumina_indexed_primers.xlsx file below .

5µL ILLPCR1 oligo forward (100µM)

5µL ILLPCR2 oligo reverse (100µM) ind01 to ind27 (one per tube)

90µL nuclease free water and mix by pipetting

Store at 4°C (or at -20°C for a long-term storage)

Ilumina indexed primers sequences (NGS grade needed) :

Illumina_indexed_primers.xlsx

Genomic DNA extraction and purification

6.

Genomic DNA extraction

Upon collection, a piece of tissue was cut out from a spot that was free of algal and animal epiphytes and stored in silica gel. Total genomic DNA was extracted from 15 to 20 mg of grinded dry tissue using the Nucleospin 96 plant kit (Macherey-Nagel, Germany).

Example of micro-algae stored in silica gel
Example of micro-algae stored in silica gel
6.1.

In a Rack of Tube Strips (consumable of NucleoSpin 96 Plant II kit) add15mg to 20mg of each sample with one 3 mm grinding ball stainless steel. Close the Tubes Strips with Cap Strips.

Example of dry algae with 3 mm steel ball before and after grinding process in individual tube
Example of dry algae with 3 mm steel ball before and after grinding process in individual tube
6.2.

Grind dry tissues using Mixer Mill MM400 Retch using 2 cycles of 0h 2m 0s at maximum frequency.

Example of grinded samples in 96 wells plate format
Example of grinded samples in 96 wells plate format
6.3.

The extraction was performed according to the manufacturer's instructions using the PL1 lysis buffer except that we added one wash step with PW1 buffer (2 times PW1 washes in total) and one wash step with PW2 buffer (3 times PW2 washes in total). The extracted DNA was eluted into 120 μL (2 x 60 µL) of the supplied elution buffer.

Instruction-NucleoSpin-96-Plant-II.pdf

Store at 4°C (or at -20°C for a long-term storage)

7.

Genomic DNA purification

The genomic DNA extracts were purified using the NucleoSpin gDNA Clean-up XS, Micro kit for DNA clean up and concentration (Macherey-Nagel, Germany).

The purifications were performed according to the manufacturer's instructions with elution into 30 μL (2 x 15 µL) of the supplied elution buffer.

7.1.

Transferring the 120µL into 1.5 mL microtubes. Add nuclease free water to fill up to 400µL.

7.2.

The purifications were performed according to the manufacturer's instructions with elution into 30 μL (2 x 15 µL) of the supplied elution buffer.

Instruction-NucleoSpin-gDNA-Clean-up-XS.pdf

7.3.

Removal of residual ethanol and concentration were performed by incubation 0h 15m 0s at 70°C

7.4.

The purified gDNA of each sample was transferred in a 96 wells PCR plate.

Store at 4°C (or at -20°C for a long-term storage)

Genomic DNA Quality and Quantification

8.

Quality Control of Genomic DNA

Optional: load 1µL to 3µL of the Genomic DNA extract on an agarose gel to evaluate its quantity and quality using electrophoresis.

Example of agarose gel picture of genomic DNA
Example of agarose gel picture of genomic DNA
9.

Quantification of Genomic DNA ( Preparation for one 96 wells PCR plate) using PicoGreen‱ TM

Quantify Genomic DNA extract using PicoGreenTM.

Qant-iT™ PicoGreen® dsDNA Quantification

9.1.

Take out all reagents from the fridge and bring them to room temperature.

Take out the DNA samples from the freezer. DNA samples should be slowly thawed on ice

Note
Quant-iTTMPicoGreen® dsDNA reagent is dissolved in dimethylsulfoxide (DMSO), which freezes below 19 °C. The reagent must be completely thawed before using it by bringing it to room temperature. After the reagent thawed, it is advisable to briefly vortex the tube to make sure it is adequately mixed and to spin it down in a centrifuge or a mini centrifuge.

Note
Quant-iTTMPicoGreen® dsDNA reagent is light sensitive and should be protected from light at all times.

9.10.

Place the plate in a plate reader and measure the fluorescence according to the following parameters:

Excitation ~480 nm

Emission ~520 nm

Integration time 40 s

Lag time 0 s

Gain Optimal

Number of flashes 10

Calculated well highest standard

Shaking 5 s

Equipment

ValueLabel
Synergy 2NAME
absorbance microplate readerTYPE
BioTekBRAND
Synergy2SKU

Equipment

ValueLabel
SPARKNAME
Microwell plate readerTYPE
TECANBRAND
SPARKSKU
9.11.

Plot the measured fluorescent values of the standard samples against their known concentrations and fit a linear curve using linear regression. Make sure that the coefficient of determination (R2) is close to 1 (typically > 0.99). Calculate the DNA concentrations in the unknown samples using the slope and intercept parameters of the linear equation. Output values you obtained are in pg/µl, assuming 1 µl of each sample was used.

9.2.

Preparation of of 1X TE buffer 11mL of 1X TE buffer

In 15mL

550µL

10.450mL

Mix by inverting the tube several times.

9.3.

Preparation of DNA solution at 5000 pg/µl (for 3 ranges)

In 0.5mL

4µL

76µL

Mix by inverting the tube several times.

9.4.

Preparation of the standard range 0 pg/µl to 1000 pg/µl

Prepare the following standard mixture in 8 0.5mL

ABCDE
TubesStandard DNA solution concentration (pg/µL)Standard DNA solution volume (µL)1X TE (µL)Final DNA concentration (pg/µL)
15000421681000
2500021189500
350009171250
4100021189100
55002118950
61001816210
750181625
8001800

Standard DNA solutions preparation Pipette 50µL of each standard mixture in the first two columns of the black, sterile, 96-well plate :

ABCDEFGH
1000 pg/µL1000 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
500 pg/µL500 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
250 pg/µL250 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
100 pg/µL100 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
50 pg/µL50 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
10 pg/µL10 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
5 pg/µL5 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA
0 pg/µL0 pg/µLunknown DNAunknown DNAunknown DNAunknown DNAunknown DNAunknown DNA

Exemple of map plate for PicoGreen‱ quantification

9.5.

Pipette 49 µl of 1X TE buffer in the remaining wells.

49µL

9.6.

Pipette 1 µl of the unknown DNA samples in the remaining wells.

1µL

9.7.

Prepare PicoGreen‱ work solution ® work solution

In 10mL

25µL

4.975mL

Mix and protect from light.

9.8.

Pipette 50µL in each well, including the standard and unknown sample wells.

9.9.

Protect the 96-well plate from light and incubate for 0h 5m 0s at room temperature.

Genomic DNA preparation

10.

In a PCR plate, put around 100ngof genomic DNA in a volume of 40µL (in nuclease free water or Tris-HCl 5mM pH 8.5) for each sample. If possible, randomize the location of samples in the microplate. Keep a few empty wells for negative controls.

Double digestion

11.

Double digest around 100ng of high quality genomic DNA with selected restriction enzymes 50µLreaction volume. Use a digestion buffer appropriate for both enzymes.  Here, we will do the protocol for the PstI and HhaI couple of enzymes but it's same with PstI and MseI couple. Both couple of enzyme woks well for micro-algae but it's possible to test double digestion on few sample to select the best couple of enzymes. The best couple given large smear with size range 100 bp to 1000 pb.

11.1.

Vortex all reagents, except enzymes (stored at -20°C), for approximately 0h 0m 5s

Spin down all reagents for approximately 0h 0m 5s and placeOn ice.

In a microtube, prepare the digestion mix, according to the following table for a total volume of 50µL:

ABCDE
Initial concentrationFinal concentrationn=1n=100 (1 plate)
Genomic DNA~100 ng40 µL
Cutsmart buffer10X1X5 µL500 µL
Enzyme 1 (PstI HF)20 u/µL10U0.5 µL50 µL
Enzyme 2 (HhaI or MseI)20 u/µL10U0.5 µL50 µL
nuclease-free water4 µL400 µL
TOTAL50 µL1000 µL

Digestion master mix composition

11.2.

Vortex the master mix and spin down.

Aliquot 125µL of the digestion master mix in each well of a 8-PCR tube strip.

In the DNA plate (containing 40µL per well), add10µL of digestion master mix with a x8 multichannel pipette and mix by pipetting, seal PCR plate and spin down.

Incubate at 37°C 0h 0m 5s

Then store at 4°C

Check digestion on an agarose gel

12.

Check the efficiency of the digestion by electrophoresis of 5µL of digested DNA in a 1.5 % agarose gel (standard quality). High molecular weigh DNA should no longer be visible.

Example of agarose gel picture after digestion of genomic DNA. Size marker is a 100bp ladder.
Example of agarose gel picture after digestion of genomic DNA. Size marker is a 100bp ladder.

Bead purification (96‑well plate format)

13.

This protocol can be used to remove contaminants, unligated adapters, enzymes, buffer additives, salts... and short DNA fragments. The method utilizes a single-size selection step : After adding the appropriate volume of Bead Suspension to the DNA sample, beads will bind larger fragments. The supernatant contains smaller fragments and contaminants that are discarded. For most NGS sequencing applications it is optimal to remove all fragments below 100 bp. This can be achieved by using a volume ratio (bead suspension to sample) of 1:1, which is described in the following protocol.

NucleoMag kit for clean up and size selection Workflow (Macherey-Nagel)
NucleoMag kit for clean up and size selection Workflow (Macherey-Nagel)
14.

Before starting

Prepare 50mL of fresh 80% Molecular Biology Grade Ethanol

40mL

10mL

Remove the NucleoMag® NGS Bead Suspension from the fridge. Let for approximately 30 min to bring the bead suspension to Room temperature.

Then, vortex this Bead Suspension stock solution carefully until homogenized and put in a reagent reservoir.

15.

Binding

This step binds DNA fragments 100 bp and larger to the magnetic beads.

Pipette 45µL of NGS Beads suspension with x8 multichannel pipette and transfer in digestion plate (plate with 45µL of digested template DNA for each sample), carefully mix by pipetting up and down 10 times.

Incubate 0h 5m 0s atRoom temperature

16.

Separation

Place the purification plate onto the 96-well magnetic separator.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

The supernatant contains unwanted low molecular weight contaminants and unwanted smaller DNA fragments.

Remove and discard the supernatant (~90 µl) by pipetting.

Note
Do not disturb the attracted beads while aspirating the supernatant. Remove the supernatant with the multichannel from the opposite side of the well.

Example of separation process where beads have been attracted by the magnets
Example of separation process where beads have been attracted by the magnets
17.

1st wash with 80 % ethanol

Place 80% ethanol in a reagent reservoir.

With a x8 multichannel pipette, dispense200µL of 80% ethanol into the purification plate without disturbing the bead pellet.

Incubate the purification plate at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

18.

2nd wash with 80 % ethanol

With a x8 multichannel pipette, dispense200µL of 80% ethanol into the purification plate without disturbing the bead pellet.

Incubate the purification plate at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

19.

Dry the beads

Let the purification plate on the magnetic separator and incubate at Room temperature for maximum 0h 5m 0s in order to allow the remaining traces of ethanol to evaporate.

Note
Take care not to over dry the bead pellet (bead pellet appears cracked in this case) as this will significantly decrease elution efficiency.

20.

Elute DNA fragments

Take the purification plate from the magnetic stand, and add 40µL of nuclease-free water with a x8 multichannel pipette to resuspend the bead pellet by pipetting up and down 10 times.

Incubate the purification plate at Room temperature for 0h 5m 0s .

Separate the magnetic beads against the side of the wells by placing the 96-well plate on the magnetic separator.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

Transfer 35µL of the supernatant containing the digested purified template DNA to a new 96‑well plate. Be careful to avoid pipeting beads during this step .

Seal the plate and store at 4°C (or store -20°C for a long-term storage) until adaptor ligation.

Adaptor ligation

21.

For each sample of one line of the digested purified plate (with 35µL of digested purified template DNA) add 5µL of double-stranded barcoded P1 adaptors at 4 µM. Use one double-stranded barcoded P1 adaptors per line .

Map of plate with the P1 adaptors lines
Map of plate with the P1 adaptors lines

Note
Depending on the number of samples you want to put into the library (maximum 576 samples), you can mix the number of barcoded P1 adaptors (one per lines) and the number of Illumina indexes (one per columns).

Number and map of plates with combinaison of barcoded P1 adaptors and Illumina Indexes depending of the number of samples required into the library.
Number and map of plates with combinaison of barcoded P1 adaptors and Illumina Indexes depending of the number of samples required into the library.
22.

Vortex all reagents, except enzymes (stored at -20°C), for approximately 0h 0m 5s

Spin down all reagents for approximately 0h 0m 5s and placeOn ice.

In a microtube, prepare the ligation mix, according to the following table for a total volume of 60µL:

ABCDE
Initial concentrationFinal concentrationn=1n=100 (1 plate)
Digested purified template DNA + P1 adaptor40 µL
P2 adaptor (HhaI or MseI)40 µM330 nM0.5 µL50 µL
T4 ligase buffer10X1X6 µL600 µL
T4 ligase400 u/µL160U0.4 µL40 µL
nuclease-free water13.1 µL1310 µL
TOTAL60 µL2000 µL

Ligation master mix composition

23.

Vortex the master mix and spin down.

Aliquot 125µL of the ligation master mix in each well of two 8-PCR tube strip.

In the digested purified plate (containing 35µL of digested purified template DNA and 5µLof barcoded P1 adaptors ), add20µL of ligation master mix with a x8 multichannel pipette and mix by pipetting, seal PCR plate and spin down.

Incubate at 16°C 0h 0m 5s

Then store at 4°C or at -20°C if not performing the bead purification the day after.

24.

Before starting

Prepare 50mL of fresh 80% Molecular Biology Grade Ethanol

40mL

10mL

Remove the NucleoMag® NGS Bead Suspension from the fridge. Let for approximately 30 min to bring the bead suspension to Room temperature.

Then, vortex this Bead Suspension stock solution carefully until homogenized and put in a reagent reservoir.

25.

Binding

This step binds DNA fragments 100 bp and larger to the magnetic beads.

Pipette 60µL of NGS Beads suspension with x8 multichannel pipette and transfer in adaptor-ligated plate (plate with 60µL of digested and adaptor-ligated template DNA for each sample), carefully mix by pipetting up and down 10 times.

Incubate 0h 5m 0s atRoom temperature

26.

Separation

Place the purification plate onto the 96-well magnetic separator.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

The supernatant contains unwanted low molecular weight contaminants and unwanted smaller DNA fragments.

Remove and discard the supernatant (~120 µl) by pipetting.

Note
Do not disturb the attracted beads while aspirating the supernatant. Remove the supernatant with the multichannel from the opposite side of the well.

Example of separation process where beads have been attracted by the magnets
Example of separation process where beads have been attracted by the magnets
27.

1st wash with 80 % ethanol

Place 80% ethanol in a reagent reservoir.

With a x8 multichannel pipette, dispense200µL of 80% ethanol into the purification plate without disturbing the bead pellet.

Incubate the purification plate at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

28.

2nd wash with 80 % ethanol

With a x8 multichannel pipette, dispense200µL of 80% ethanol into the purification plate without disturbing the bead pellet.

Incubate the purification plate at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

29.

Dry the beads

Let the purification plate on the magnetic separator and incubate at Room temperature for maximum 0h 5m 0s in order to allow the remaining traces of ethanol to evaporate.

Note
Take care not to over dry the bead pellet (bead pellet appears cracked in this case) as this will significantly decrease elution efficiency.

30.

Elute DNA fragments

Take the purification plate from the magnetic stand, and add 40µL of nuclease-free water with a x8 multichannel pipette to resuspend the bead pellet by pipetting up and down 10 times.

Incubate the purification plate at Room temperature for 0h 5m 0s .

Separate the magnetic beads against the side of the wells by placing the 96-well plate on the magnetic separator.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

Transfer 35µL of the supernatant containing the adaptor-ligated purified template DNA to a new 96‑well plate. Be careful to avoid pipeting beads during this step .

Seal the plate and store at 4°C (or store -20°C for a long-term storage) until PCR amplification.

Rad Tag enrichment (PCR)

31.

PCR amplification to generate Illumina sequencing indexed libraries :

In this PCR, Illumina indexed primers are incorporated in order to produce fragments compatible with Illumina sequencing, and to insert an index allowing multiplexing of barcoded samples. This index will be read during the sequencing run.

This PCR is expected to have a homogenizing effect. Primers are thus included in limiting quantity, in order to produce equalized amounts of PCR fragments among samples. The number of cycles is limited to a maximum of 15 (optimal with 12). After those cycles, a final PCR cycle is then performed after addition of primers in large excess.

The Reaction mixture for a total volume of 40µL is :

ABCD
Initial concentrationFinal concentrationn=1
Adaptor-ligated purified template DNA10 µL
Primer mix (ILLPCR1 and ILLPCR2ind)5 µM each0.17 µM1.36 µL
Q5 buffer5X1X8 µL
High GC enhancer5X1X8 µL
dNTP mix25mM each0.20 µM0.32 µL
Q5 hotstart hifi polymerase2 u/µL0.8 U0.40 µL
nuclease-free water11.92 µL
Total mix30 µL
TOTAL reaction40 µL

PCR mixture composition

We need to prepare one PCR mixture per index, i.e. 12 PCR mixtures for one plate.

Map of plate with the Illumina indexes columns
Map of plate with the Illumina indexes columns
32.

First PCR mix preparation (with primers in limiting quantity)

Defreeze and vortex all reagents, except enzymes (stored at -20°C), for approximately0h 0m 5s

Spin down all reagents for approximately0h 0m 5s and place On ice.

In 12 0.5mL microtubes, prepare the 1st mix according to the following table (one mix per column) :

ABCDE
Initial concentrationFinal concentrationn=1n=10 (one column of 1 plate)
Adaptor-ligated purified template DNA10 µL
Primer mix (ILLPCR1 and ILLPCR2ind)5 µM each0.17 µM1.36 µL13.6 µL
Q5 buffer5X1X8 µL80 µL
High GC enhancer5X1X8 µL80 µL
dNTP mix25mM each0.20 µM0.32 µL3.2 µL
Q5 hotstart hifi polymerase2 u/µL0.8 U0.40 µL4 µL
nuclease-free water11.92 µL119.2 µL
TOTAL30 µL300 µL

PCR mix composition

Vortex mix all reagents in the mix and spin down.

In a new PCR plate, dispense 30µL of 1st mix in each column.

33.

DNA and mix combination

Spin down the adaptor-ligated purified template DNA plate.

With a multichannel pipette, transfer 10µL of adaptor-ligated purified template DNA into the PCR plate and mix by pipetting.

Finally, aliquot the 40µL of the total mix by dispensing 20µL into 1 additional new empty PCR plates.

Seal the 2 PCR plates and spin down.

The 2 PCR will be performed in parallel in 2 different thermal cyclers, in order to reduce the PCR bias.

ABCD
Cycle stepTemperatureTimeCycles
Hot start initial denaturation98°C30 sec1
Denaturation98°C20 sec15
Annealing60°C30 sec15
Extension72°C40 sec15
Final extension72°C10 min1
Hold4°C

PCR program for the Illumina indexing PCR

After PCR, pool back the 2 PCR plates into a single plate with a multichannel

34.

Final cycle (with primers in large excess)

In a 12-tube PCR strip, prepare the 2nd mix according to the following table (one mix per column) :

ABCDE
Initial concentrationFinal concentrationn=1n=10 (one colomn of 1 plate)
Primer mix (ILLPCR1 and ILLPCR2ind)5 µM each3.35 µM2.68 µL26.8 µL
Q5 buffer5X1X0.80 µL8 µL
dNTP mix25mM each0.20 µM0.32 µL3.2 µL
nuclease-free water0.20 µL2 µL
TOTAL4 µL40 µL

Final cycle PCR mix composition

Mix all reagents by pipetting and spin down.

Dispense 4µL of final cycle mix in each line of the PCR plate with a 12 multichannel pipette.

Seal the PCR plate and spin down.

In a thermocycler, run the final cycle as follows :

ABCD
Cycle stepTemperatureTimeCycles
Denaturing98°C3 min1
Annealing60°C2 min1
Extension72°C12 min1
Hold12°C

PCR program for the final cycle of the illumina PCR

After PCR, place the plate at4°C (or -20°C for a long-term storage).

Check PCR on an agarose gel

35.

Check the efficiency of the PCR by electrophoresis of 5µL of digested DNA in a 1.5 % agarose gel (standard quality).

Example of smears obtained after the Illumina PCR. Size marker is a 100bp ladder.
Example of smears obtained after the Illumina PCR. Size marker is a 100bp ladder.

Sample pooling (equimolar)

36.

Each barcoded and indexed individual can now be pooled in a single tube, in equimolar conditions.

After the normalizing PCR, all smears should have similar intensity on the agarose gel. In this case, pool 5µL of all individuals in a single low binding 1.5mL microtube.

If not, normalization can be made at this step. For this, roughly estimate the concentration of fragments from the gel picture, and pools accordingly. It can be efficient to make intermediate pools for example, one pool for the low, one for the medium, and another one for the high intensity samples in 3 low binding 1.5mLmicrotubes .

Vortex mix and spin down.

Store at4°C (or -20°C for a long-term storage) until bead purification.

Bead purification (microtube format)

37.

Before starting

Prepare 10mLof fresh 80% Molecular Biology Grade Ethanol

8mL

2mL

Remove the NucleoMag® NGS Bead Suspension from the fridge. Let for approximately 30 min to bring the bead suspension to Room temperature.

Then, vortex this Bead Suspension stock solution carefully until homogenized and put in a reagent reservoir.

38.

Binding

This step binds DNA fragments 100 bp and larger to the magnetic beads.

Pipette a volume of NGS Beads suspension to have a ratio 1:1 with the sample pooling volume, and transfer in the pooling sample tube(s).

Carefully mix by pipetting up and down 10 times.

Incubate 0h 5m 0s atRoom temperature

39.

Separation

Place the purification tube(s) onto the magnetic microtube stand.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

The supernatant contains unwanted low molecular weight contaminants and unwanted smaller DNA fragments.

Remove and discard the supernatant by pipetting.

Note
Do not disturb the attracted beads while aspirating the supernatant. Remove the supernatant with the pipette from the opposite side of the well.

40.

1st wash with 80 % ethanol

Dispense1mL of 80% ethanol into the purification tube(s) without disturbing the bead pellet.

Incubate the purification tube(s) at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

41.

2nd wash with 80 % ethanol

Dispense1mL of 80% ethanol into the purification tube(s) without disturbing the bead pellet.

Incubate the purification tube(s) at room temperature for at least 0h 0m 30s

Carefully and completely remove and discard ethanol by pipetting.

42.

Dry the beads

Let the purification tube(s) on the magnetic separator and incubate at Room temperature for maximum 0h 5m 0s in order to allow the remaining traces of ethanol to evaporate.

Note
Take care not to over dry the bead pellet (bead pellet appears cracked in this case) as this will significantly decrease elution efficiency.

43.

Elute DNA fragments

Take the purification tube(s) from the magnetic stand, and add 100µL of nuclease-free water with a pipette to resuspend the bead pellet by pipetting up and down 10 times.

Incubate the purification tube(s) at Room temperature for 0h 5m 0s .

Separate the magnetic beads against the side of the tube by placing the tube(s) on the magnetic separator.

Wait at least 0h 5m 0s until all the beads have been attracted by the magnets or until the liquid appears clear.

Transfer 90µL of the supernatant in a new(s) low binding 1.5mL microtube(s). Be careful to avoid pipeting beads during this step .

Store at 4°C (or store -20°C for a long-term storage) until size selection.

Intermediate pools quantification and final pooling (Optional)

44.

If you made intermediate pools (low, medium and high intensity on agarose gel), estimate the double strand DNA concentration in the pools by fluorimetry with a Qubit equipment.

Qubit_dsDNA_HS_Assay_UG.pdf

Pool in equimolar concentration the 3 intermediate pools in a single low binding 1.5mLmicrotube with a minimum concentration of 20nanomolar (nM) and minimum volume of 30µL.

Vortex mix and spin down.

Store at4°C (or -20°C for a long-term storage) until size selection.

Size selection with sage science Pippin-Prep

45.

Perform the size selection of fragments between 300 and 800 bp using a 1,5% DF marker K agarose gel cassette, according to the Pippin prep manufacturer's instructions : 

Pippin-prep-Quick-Guide-CDF1510-marker-K3.pdf

Note
In the case you do not have access to a Pippin prep, you can alternatively perform a double size-selection with beads (see the bead manufacturer's instructions for details), or by smear excision and purification from an agarose gel. From our experience, size selection with the Pippin prep is the most accurate and repeatable method.

Quality control of the libraries

46.

Control the quality of the library with a Bioanalyzer (Agilent) (or equivalent equipment) in a High Sensitivity DNA chip. Dilute your pool 1:2 or more and load 1µl of the pool before and after size selection, according to the manufacturer's instructions :

Agilent_high_SensitivityDNA_KG_EN.pdf

Example of a Bioanalyzer profile obtained for a ddRADseq library before size selection (DNA HS kit)
Example of a Bioanalyzer profile obtained for a ddRADseq library before size selection (DNA HS kit)
Example of a Bioanalyzer profile obtained for a ddRADseq library after size selection (DNA HS kit)
Example of a Bioanalyzer profile obtained for a ddRADseq library after size selection (DNA HS kit)
47.

Fluorimetric estimation of the dsDNA concentration in the library.

Perform a quick estimation of the DNA concentration with a fluorimetric assay, in a QubitTM apparatus or equivalent, with the Qubit ds 1X DNA HS assay kit, according to the manufacturer's instructions:

Qubit_dsDNA_HS_Assay_UG.pdf

48.

qPCR quantification

In the case you need an accurate estimation of the DNA concentration in your library, perform a qPCR quantification with the NEBNext Library Quant Kit for Illumina, or equivalent, which uses P5 and P7 illumina primers to target the double stranded DNA fragments in the library. Follows the kit's user guide and perform your quantitative qPCR in a qPCR thermocycler (e.g. LightCycler 480, Roche).

NEBNext_manualE7630.pdf

Contrarily to the fluorimetric method (Qubit), the qPCR estimation will only consider dsDNA fragments starting with P5 and ending with P7 illumina sequences, that will be effectively amplified onto the flowcell of the Illumina sequencer.

48.1.

Suggestions to prepare library dilutions for qPCR

In 0.5mL low binding microtube, prepare 1:1 000 dilution of library with buffer supplied in the qPCR kit. Then, prepare the 3 library dilutions (1:10 000 to 1:30 000) to be used on triplicate for qPCR analysis.

1:100 : 1µL of library + 99µL of 1X buffer

1:1 000 : 10µL of library + 90µL of 1X buffer

1:10 000 : 20µL of 1:1 000 dilution + 180µL of 1X buffer

1:20 000 : 50µL of 1:10 000 dilution + 50µL of 1X buffer

1:30 000 : 50µL of 1:10 000 dilution + 100µL of 1X buffer

You should get more than 10 nM, that is the library concentration usually required by the sequencing platform facilities.

Use the average size of the library size range as estimated from the Bioanalyzer profile to convert DNA concentration from nM to ng/µL using the attach file below :

nM_ngµL_Conversion_Calculator.xlsx

49.

Library ready for sequencing

The library is now ready for sequencing in single read or paired-end 150 bases in an Illumina sequencer, with one index read.

推荐阅读

Nature Protocols
Protocols IO
Current Protocols
扫码咨询