High-throughput Method for Determination of Seed Paternity by Microsatellite Markers

Materials and Reagents

1. Agencourt Chloropure Kit (this kit is discontinued by Beckman Coulter but can be ordered) or Quiagen MagAttract 96 DNA plant core kit (QIAGEN, catalog number: 67163 )
   
2. Qiagen Type-it Multiplex PCR kit (QIAGEN, catalog number: 206243 )
   
3. QIAquick PCR purification kit (QIAGEN, catalog number: 28106 )
   
4. Liquid nitrogen
   
5. 100% Isopropanol, ultrapure
   
6. Freshly prepared 70% ethanol made with nuclease free water
   Note: 70% ethanol is hygroscopic, always make fresh for optimal result.
   
7. Agarose
   
8. Ethydium bromide solution
   
9. Nuclease free water
   
10. Pure molecular biology grade ethanol (96– 100%)
    Note: Calculate the volume of ethanol needed before you start. A new kit of 24 plates needs about 5 L of ethanol.
    
11. GENSCAN 500 ROX (Applied Biosystems, catalog number: 401734 )
12. 1x TAE (Tris-Acetate-EDTA) buffer (see Recipes)
    

Equipment

1. 2.2 ml Ritter Deep-well plates (ABGene, catalog number: AB-6661 )
   
2. Adhesive plate film (ABGene, catalog number: AB-0558 , AB-626 , AB-662 )
   
3. 4 mm steel balls (SPEX SamplePrep, catalog number: 2150 )
   
4. Pipette tips
   
5. Processsing Plate: Deep-well titer plate with a 96-well format (SPEX SamplePrep, catalog number: 2210 )
   
6. Destination Plate: 300 μl round bottom microtiter plate (96 well, 300 μl well capacity, round bottom) (Corning, Costar^®, catalog number: 07-200-105 )
   
7. Desktop centrifuge for 96 well plates
   
8. Reagent reservoir
   
9. Genogrinder (SPEX SamplePrep, catalog number: 2010 )
   
10. Cooling block (SPEX SamplePrep, catalog number: 2665 )
    
11. Steel ball dispenser (SPEX SamplePrep, catalog number: 2100 )
    
12. Multi-channel pipettes (0.1-10, 10-100, 30-300 μl)
    
13. Agencourt Supermagnet magnetic Plate (Agencourt SPRIPlate 96R-Ring Magnet Plate, catalog number: A29164 )
    
14. BioRad Gel casting tray, running tray, power pack etc (Bio-Rad Laboratories)
    
15. Nanodrop spectrophotometer (Nanodrop)
    
16. Pipette (0.2-2 μl), pipette tips, soft tissue paper
    
17. ABI 3100 Genetic Analyzer (Applied Biosysyems)
    

Procedure

1. Sample grinding (Day 1)
   1. Take up to 40 mg fresh tissue sample (leaf punches/seedling) in 2.2 ml plate. Take care to group similar samples/treatments/experimental designs together for easier data handling during sequencing and genotype analysis.
      
   2. Dispense 2 steel balls in each well with dispenser.
      
   3. Balance two plates exactly for grinding in Genogrinder.
      
   4. Seal the plate(s) with adhesive plate film.
      
   5. Cool the plate(s) in liquid nitrogen.
      
   6. Cool the cooling racks along with the plates.
      
   7. Grind in Genogrinder at Speed 250, 1x, for 1 min. If not crushed completely, again for 1 min.
      
      
2. DNA isolation (Day 1)
   Assembly in step B-1 is performed once for each new Agencourt Chloropure kit. If you have already made the following preparations for a previous experiment, please skip ahead to step B-2. Standard RNase treatment can be included either at the step B-5 with the homogenized lysate (requires large amount of enzyme, but convenient for high-throughput) or selectively after determination of the degree of RNA contamination at "quality check".
   
   1. Add 80 ml of 100% Isopropanol to the wash buffer bottle provided with the kit. After addition of Isopropanol invert the bottle to mix. Once the solution has been thoroughly mixed, store at room temperature.
      
   2. Prepare bind buffer: Combine 6 μl bind buffer with 150 μl of 100% isopropanol for each individual isolation in a nuclease free vessel of suitable size (for example: for 10 isolations, add 60 μl of bind buffer to 1.5 ml of 100% isopropanol in a 15 ml conical tube). Vortex bind buffer bottle thoroughly before combining. Unused combined solution should be discarded.
      
   3. Homogenize each sample in 300 μl of lysis buffer, provided with the kit. (sample input should not exceed three 6 mm lyophilized leaf punches or 40 mg ground seeds or fresh leaf material.)
      
   4. Centrifuge lysate for 10 min at 1,100 RCF at room temp to pellet debris.
      
   5. Transfer 150 μl of homogenized lysate to 1.2 ml processing plate.
      
   6. Pipette 150 μl bind buffer (prepared in step B-2), mix by slowly pipetting 5 times and incubate at room temperature for 5 min.
      
   7. Move the plate onto the Agencourt Supermagnet and separate for 2-4 min. Wait for the solution to clear before proceeding to the next step.
      
   8. Slowly aspirate the cleared solution from the plate and discard.
      Note: Aspirate from the top down to avoid disturbing the pellet. This step must be performed while the plate is situated on the magnet. Do not disturb the separated magnetic beads. If beads are drawn out, leave a few microliters of supernatant behind.
      
   9. Remove the plate from the magnet and add 300 μl of Wash buffer. Pipette mix 10 times and incubate for 1 min at RT.
      
   10. Return plate to the magnet and separate for 2-4 min. Wait for the solution to clear before proceeding to the next step.
       
   11. Slowly aspirate the cleared solution from the plate and discard. This step must be performed while the plate is situated on the magnet.
       
   12. Remove the plate from the magnet and add 300 μl of 70% ethanol. Pipette mix 10 times to re-suspend the beads.
       
   13. Return plate to the magnet and separate for 2 min. Wait for the solution to clear before proceeding to the next step.
       
   14. Slowly aspirate the cleared solution from the plate and discard.
       
   15. Repeat steps l-n for a total of 2 ethanol washes.
       
   16. Let the plate air dry for 5 min at room temperature. The plate should air-dry until the last visible traces of ethanol evaporate. Over-drying the sample may result in a lower recovery.
       
   17. Remove the plate from the magnet and add 50 μl of nuclease-free water. Re-suspend the beads by pipette mixing 10 times.
       Note: Smaller or larger elution volumes can be used for more or less concentrated product; however the minimum elution volume should be 40 μl to ensure complete elution. Optimal elution volumes need to be experimentally determined.
       
       1. DNA quality/quantity check
          1) Quality
          a) Prepare 0.8% agarose gel with 1x TAE buffer. Add 2 μl Ethydium Bromide before casting.
          b) Load 2 μl DNA+1 μl loading dye
          c) Add 1 kb or 100 bp ladder for reference
          d) Run at constant 100 V for 30 min
          e) Check DNA for quality, presence of RNA and quantity in gel doc system and record.
          
       
       
       Figure 1. DNA quality check on 0.8% agarose gel. Lane 1-6, 8-13 represent extracted DNA samples, lane 7 represents DNA marker.
       
       2) Quantity

       a)	Follow the procedure to measure concentration of DNA in Nanodrop.
       b)	Export to excel file for reference.

       2. Dilution of DNA sample
          a) Dilute the samples to 50 ng/μl as working stock.
          b) Dispense them in 3 PCR plates (for 3 multiplex PCR reactions) similarly and maintaining the sample organization of working DNA plate. Make a record of the samples.
          
          
3. PCR amplification (Day 2)
   1. Multiplex groups
      Primer groups are selected based on primer sequences using software MultiPLX version 2.0 (available at http://bioinfo.ut.ee/?page_id=167, Ref: doi: 10.1093/bioinformatics/bti219) to select multiplex grouping. Medium stringency was chosen to select groups.
      Example of 3 multiplex groups of 8 primer pairs:
      
      Table 1. Example of 3 multiplex groups of 8 primer pairs:
      
      
      
   2. 10x Primer mix
      1. First make 100 μM primer stock (if using 50 μM primer stock, adjust the volume of the table below accordingly).
         
      2. Mix all primers according to the table below to get 3 primer-mix groups (sufficient for 1x PCR plate i.e. 96 samples).
         
         Table 2. Schematic representation of preparation of 3 primer mix groups:
         
         
         
   3. PCR reaction mix
      1. Thaw the 2x Type-it Multiplex PCR Master Mix (MM) (if stored at -20 °C), template DNA (already dispensed in PCR plates), RNase-free water, and the primer mix. Mix the solutions completely before use.
         Note: It is important to mix the solutions completely before use to avoid localized concentrations of salts.
         
      2. Prepare a reaction mix according to the table below. Note: The reaction mix typically contains all the components required for multiplex PCR except the template DNA. Prepare a volume of reaction mix 10% greater than that required for the total number of reactions to be performed. For reaction volumes less than 25 μl, the 1:1 ratio of Type-it Multiplex PCR Master Mix to primer mix and template should be maintained.
         Note: Starting with an initial Mg²⁺ concentration of 3 mM is recommended as provided by the 2x Type-it Multiplex PCR MM.
         
         Table 3. Reaction Mix preparation for one 96 well plate samples
         (modified after manufacturer's instruction)
         
         
         
      3. Mix the reaction mix thoroughly and dispense appropriate volumes into PCR tubes or plates. Note: Mix gently, for example, by pipetting the reaction mix up and down a few times. Due to the hot start, it is not necessary to keep samples on ice during reaction setup.
         
   4. PCR condition
      
      Table 4. PCR condition (according to manufacturer's instruction)
      
      
      
   5. PCR programme:
      Ctrl. Tube
      Lid- 105 °C
      (1)95 °C-05:00 min
      (2)95 °C-00:30 min
      (3)60°C- 01:30 min
      (4)72 °C-00:30 min
      (5)Go to step (2), repeat 29
      (6)60°C-30:00 min
      (7)20 °C- ∞
      (8)End
      (9)Store PCR products in -20 °C, till further processing.
      
   6. PCR check
      Check PCR products randomly on 1.8% agarose gel (this is just to check if everything is working before proceeding to next steps. So choose random samples). Load 4 μl product and 100 bp ladder.
      Note: Use only Xylene Cyanol loading dye as Bromophenol Blue dye interfere with the amplified band range.
      
      
4. PCR purification (Day 3)
   1. Add 75 μl nuclease free water to each well of PCR plate (containing 25 μl PCR product) to make it 100 μl.
      
   2. In a 96 well plate (at-least 1 ml capacity) dispense 300 μl Buffer PM. Add 100 μl PCR sample and mix by pipetteing.
      
   3. Apply the samples to the wells of the QIAquick plate while it is fitted with the collecting deepwell plate.
      
   4. Centrifuge at 1,100 RCF (g) for 2 min.
      
   5. Discard flow-through.
      Note: If the sample volume is more than 600 μl apply the remaining samples to the wells of the QIAquick plate and repeat step 4.
      
   6. Wash wells of QIAquick plate by adding 900 μl of Buffer PE to each well and centrifuge at 1,100 RCF (g) for 2 min.
      
   7. Repeat step 4-f.
      
   8. Discard flow-through and gently soak any buffer remaining on the tip of column with soft tissue.
      
   9. Centrifuge at 1,500 RCF (g) for 4 min.
      Important: This step removes residual Buffer PE from the membrane.
      
   10. Let the plate dry for 5 min or until all the ethanol (of Buffer PE) evaporates. Place it under fume hood for 10 min is sufficient.
       
   11. To elute, place the plate on collection tubes (provided), add 50 μl of RNase-free water (provided) to the center of each well of the QIAquick 96 plate, incubate for 2 min, and centrifuge at 1,500 RCF (g) for 2 min.
       Important: Ensure that the elution buffer is dispensed directly onto the center of QIAquick membrane for complete elution of bound DNA. Please note that the average eluate volume is 60 μl from 80 μl elution buffer volume, and 40 μl from 60 μl elution buffer volume.
       
       1. Quality check
          Check Purified PCR product by nanodrop or 1.8% agarose gel.
          Note: For nanodrop analysis, centrifuge the samples at 1,500 RCF for 2 min and estimate from the upper layer (to precipitate any remaining silica particles from the column, which sometimes interfere with sensitive nanodrop estimation). Imporant: It is better to check randomly the purified product on 1.8% agarose gel to be sure of quality. As this is a multiplex PCR, do not expect a clear band. The gel will appear as smear of many bands.
          
          
          Figure 2. PCR amplification check on 1.8% agarose gel. Lanes 1-15 represent PCR products (with a couple of failed amplification on lane 7 and 14), lane 16 represents 100bp DNA marker.
          
          
       2. Dilution
          It is not needed to further dilute the samples if you start with 50 ng DNA sample and followed the procedure exactly. At this step the concentration of samples should be around 250 ng/μl.
          For the size sequencing set up, it is recommended to dilute samples in this step to facilitate sample handling.
          Add 200 μl nuclease free water to each purified sample tube (1/5th dilution). So the concentration becomes 50 ng/μl.
          
          
5. Size sequencing (Day 4)
   1. Sequencing mix
      The size sequence mix for each sample contains
      
      1. 5 μl purified PCR product (250 ng total)
         
      2. 4.5 μl nuclease free water
         
      3. 0.5 μl GENSCAN 500 ROX marker
         For a single sequence plate (prepare for 105 samples), set up MM as below:
         
      1. Mix 472.5 μl nuclease free water with 52.5 μl of ROX marker.
         
      2. Dispense 5 μl of above MM to each of sequencing plate.
         
      3. Mix in 5 μl purified PCR samples.
         
      4. Do not leave any sample well blank. Fill it with Nuclease free water.
         
      5. Cover the plate with adhesive seal mat.
         
      6. Centrifuge the plates briefly at 1,500 RCF.
         
   2. Prepare sequencing table of samples.
      Prepare the sample list conveniently for recognition later. e.g. Sample name_primer mix group_Sample number_any other info.
      
      Table 5. Sample list example
      
      
      

6. Analysis (Day 5)
   Before starting analysis, import the data files generated by the sequencer to the computer with GeneMapper software. It is recommended to read the manual if you want to modify and fine-tune the marker settings and allele size binning for genotyping.
   1. Genotype table
      1. Genotyping in GeneMapper
         1)  Go to 'file' menu and create new project. Name it conveniently.
         2)  Import your samples that belong to one primer mix group.
         3)  Create marker panel and bin set for each multiplex primer groups.
         Refer to the manual for modifying or creating a new marker panel for new set of markers/primer mix in microsatellite. This involves:a) Marker set creation
         b)  Autobinning of alleles
         c)  Manual binning of alleles
         d) Create new marker panel and bin set.
         4)  Use predefined marker panel for respective samples
         
         1. Go to 'file' menu and create new project. Name it conveniently.
            
         2. Import your samples that belong to one primer mix group.
            
         3. Create marker panel and bin set for each multiplex primer groups
            
         4. Refer to the manual for modifying or creating a new marker panel for new set of markers/primer mix in microsatellite. This involves:
            
            1. Marker set creation
               
            2. Autobinning of alleles
               
            3. Manual binning of alleles
               
            4. Create new marker panel and bin set.
               
         5. Use predefined marker panel for respective samples
            
      2. Analysis parameters
         1)  Select the table setting as ‘Microsatellite default’
         2)  Select the analysis method as ‘Microsatellite analysis method’
         3)  Set the predefined panel as EV1 for primer mix 1 and so on.
         4)  Set size standard as GEN500.
         5)  Analyze
         6)  If some samples fail to match size standards, go to edit size standards and override size. If it does not improve the sizing, then there was some problem with the sample.
         a)  Open the genotype Tab. Do not bother about the GQ (genotype quality) column as it is set to a very high stringency.
         b)  Go to tools-panel manager-microsatellite kit-attenuata panel-EV1 binset. Select ‘show project alleles’. If you see any ‘*’ out of the grey bars (defines the range of allele bins), then either modify the range of existing bins (click on the grey bar and you’ll get handles for modification), or add a new bin (click on the blank space near your odd allele (that is not binned) and name it. Save the panel. Again analyze the samples with the modified binset.
         c)  If you want to see data for only two alleles and want to modify the way you want your data, modify the table output style. However, manual inspection of allele size, height etc is sometimes necessary to ensure optimal result than automation.
         d)  Go to file and export the genotype table in .csv or other excel compatible format.
         
   2. Paternity analysis in COLONY
      For paternity/maternity analysis, i.e. to determine the parent of an offspring the software COLONY (http://www.zsl.org/science/research/software/colony,1154,AR.html).
      COLONY is a Fortran program written by Jinliang Wang . It implements a maximum likelihood method to assign sibship and parentage jointly, using individual multilocus genotypes at a number of codominant or dominant marker loci. (Jones and Wang, 2010).
      1. Formatting of genotype table.
         1)  Keep only the sample name, two allele sizes for each marker for all the marker.
         2)  Compile data for all the marker sets into a single datasheet for all offspring and all marker.
         
      2. Marker table (save as tab delimited .txt file and rename 'marker').
         
         Table 6. Marker details
         
         
         
      3. Parent table (save as tab delimited .txt file and rename 'father/mother')
         
         Table 7. Paternal allele list
         
         
         Table 8. Maternal allele list
         
         
         
      4. Offspring table (save as tab delimited .txt file and rename 'offspring')
         
         Table 9. Offspring allele list
         
         
         
      5. Analysis parameters
         Feed all the information step by step as asked by the software and brose for the files to upload.
         A typical parameter is as follows:
         Number of loci:16
         Number of offspring in the sample:
         Number of male candidates:0
         Number of female candidates:
         Number of known paternal sibships:0
         Number of known maternal sibships:0
         Number of offspring with excluded fathers:0
         Number of offspring with excluded mothers:0
         Male mating system:Polygamy/monogamy
         Female mating system:Polygamy/monogamy
         Number of threads : 1
         Number of Excluded Paternal Sibships:0
         Number of Excluded Maternal Sibships:0
         Seed for random number generator:1,234
         Allele frequency:No updating by accounting for the inferred relationship
         Species:Diploid/haploid
         Sibship size prior:No
         Known population allele frequency:No
         Number of run:1
         Length of run : short/medium/long
         Monitor intermiediate results by:Every 1 second
         Prob. a mum is included in the female candidates:0.5
         Project data input produced:
         NOTE to the Project:
         
      6. Analysis outputs

         1)	Select Show result-'best configuration'.
         2)	Copy the table of assigned father ID/mother ID to each sample.
         3)	Sort parents and count offsprings sired by them in excel. Population analysis in GenAlEx

   3. Population analysis in GenAlEx
      When comparing populations for genotypes, or simply want to look at the genotypic diversity within a population, apart from parental analysis, it is better to use the excel addon 'GenAlEx' (http://biology.anu.edu.au/GenAlEx/Welcome.html). (Ref: Peakall R. and Smouse P.E. 2006. Genalex 6: Genetic analysis in Excel. Population genetic software for teaching and research, Molecular Ecology Notes, Volume 6, Issue 1, pages 288–295).
      This can analyze:
      1. Heterozygosity, F-statistics and Polymorphism
         
      2. Allelic Patterns
         
      3. Allele List
         
      4. Allele Frequencies and Principal Component Analysis (PCA).
         

Recipes

1. 1x TAE (Tris-Acetate-EDTA) buffer
   1. Prepare a stock solution of EDTA (ethylenediamine tetraacetic acid)
      For a 500 ml stock solution of 0.5 M EDTA, dissolve 93.05 g EDTA disodium salt (FW = 372.2) in 400 ml deionized water and adjust the pH with NaOH to 8.0. Make up the final volume to 500 ml.
      
   2. Prepare a 50x Stock Solution of TAE
      Dissolve 242 g Tris base (FW = 121.14) in approximately 750 ml deionized water. Add 57.1 ml glacial acid and 100 ml of 0.5 M EDTA (pH 8.0) and adjust the solution to a final volume of 1 L. This stock solution can be stored at room temperature for a long time. The pH of this buffer needs no further re-adjustment and remains about 8.5.
      
   3. Prepare a Working Solution of 1x TAE
      Dilute the 50X stock solution by 1:50 with deionized water. Final solute concentrations are 40 mM Tris acetate and 1mM EDTA.
      

Acknowledgments

This work was supported by the Max Planck Gesellschaft. The protocol was adapted from the publication: Kessler et al. (2012).

References

1. Bindler, G., van der Hoeven, R., Gunduz, I., Plieske, J., Ganal, M., Rossi, L., Gadani, F. and Donini, P. (2007). A microsatellite marker based linkage map of tobacco. Theor Appl Genet 114(2): 341-349.
   
2. Jones, O. R. and Wang, J. (2010). COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour 10(3): 551-555.
   
3. Kessler, D., Bhattacharya, S., Diezel, C., Rothe, E., Gase, K., Schöttner, M. and Baldwin, I. T. (2012) Unpredictability of nectar nicotine promotes outcrossing by hummingbirds in Nicotiana attenuata. Plant J 71(4): 529-538.