• evo视讯(中国)

    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
    AI 解读
    浏览 1510

    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 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
    扫码咨询