Population Genetics Part II: Tips and Tricks, Multiplex PCR and Workflow of Microsatellites- the cheap way

No this blog is not a belated April fools joke… there really is a method to save thousands of dollars on microsatellite marker multiplex and genotyping! If you are just reading my blog for the first time, this is Part II following up on my last blog: How-to use microsatellites for population genetics, Part I: Study Design, DNA extraction, Microsatellite Marker Design/Outsourcing.

When I first set out to work with microsatellites, I was on a budget and I had never had experience with multiplex PCR, genotyping or fluorescent markers before- so there was definitely an uphill learning curve.

To start (assuming you already have your microsatellite markers- see previous blog for this)- the next step is to order your fluorescent markers to see how things work in multiplex PCR. FYI-Don’t order your Liz Size Standard until you are done troubleshooting and checking things on gels because it expires kinda quick! Plus it ships really fast (at least if you are in the USA). 

What are fluorescent markers you ask? If you look at the microsatellite genotyping peaks above, the different color peaks correspond to the different fluorescent tails that are ‘attached’ to the microsatellite primers. FAM= blue, PET = red, VIC =green, and NED = yellow. This way, you can see the different colors and know what markers they correspond to. The orange peaks are from the 600 Liz size standard which enables you to actually calibrate the size of the peaks. This is very easy if you use a program like Geneious as they have a microsatellite plugin.  I will detail more of the data handling in Part III of this blog series.

Why multiplex PCR? Simple.. it is faster (& cheaper if you troubleshoot things quick).PCR Reaction for Multiplex PCR of Microsatellite Markers

Many folks use multiplex genotyping, where you do singleplex PCR with all the separate microsatellite markers, and then you add for example 4 non-overlapping microsatellite marker amplified products (from your singleplex PCR) into a well together for downstream genotyping (more on that process later). This process saves money, since instead of genotyping a single marker for 1$.. you can genotype 4 markers for 1$! However- singleplex PCR takes forever if you have a lot of samples and markers!!!! If you have 400 individual DNA extracts, and each sample requires 10 markers of genotyping.. this means 4000 PCR REACTIONS- YIKES- that would result in carpal tunnel in a heartbeat!

Multiplex PCR in contrast allows you to add 2-4, (or more) microsatellite markers with fluorescent tails into the PCR mix, making sure that markers with the same fluorescent tails don’t overlap in size (ie- a FAM marker amplified product of 100-250 bp compared to another FAM marker amplified product of ~300-400 bp size range- should be fine to put together in a mixture). Whereas you don’t need to worry as much if they are similar size but have different fluorescent markers (such as FAM (blue) versus NED (yellow)). There are pull-up issues, and inhibition issues.. but that is why you will need to test everything out first anyhow before running your final assays.

As for myself, two papers were key in learning how to streamline microsatellite multiplex and genotyping: Blacket et al. 2012 and Culley et al. 2013. Both studies utilized four universal fluorescent tails (different ones in the different studies)- so that all you need to do is to add the non-fluorescent tail (just the ATCGs) of the corresponding fluorescent tail (the ATCGs + the FAM, VIC, NED or PET fluorescent marker) to your forward primer, and then use pig-tails on your reverse primers (such as GT, GTT, GTTT- depending on your reverse primer). Cullen et al. 2013 has an appendix which actually walks you through every step, including the reaction concentration of each Forward primer +tail, Reverse primer +pig-tail, and Fluorescent marker +tail in the final multiplex mix. I ended up using the four universal tails in Blacket et al. 2012, and then used Culley’s reaction mixes.

Screenshot 2018-04-01 12.55.20
Universal Tails used with PCR fluorophores, from Blacket, M.J., Robin, C., Good, R.T., Lee, S.F., Miller, A.D. 2012. Universal primers for fluorescent labelling of PCR fragments–an efficient and cost-effective approach to genotyping by fluorescence. Molecular ecology resources 12, 456-463.

Below is a great pictorial image of how this works (from Blacket et al. 2012)

Screenshot 2018-04-01 12.54.34
Multiplex PCR with universal tails: Process from Blacket, M.J., Robin, C., Good, R.T., Lee, S.F., Miller, A.D. 2012. Universal primers for fluorescent labelling of PCR fragments–an efficient and cost-effective approach to genotyping by fluorescence. Molecular ecology resources 12, 456-463.

In addition to reading these papers (and their supplementary material) thoroughly, I recommend the following: 

  1. Talk to as many people as you can before you start/ while you are getting started -you always learn fabulous tips and tricks as well as what not to do!
  2. Use the multiplex manager program -this will help you simulate what markers are compatible with each other based on the estimated product sizes, the melting temperatures (Tm) and the specific tails and flourophores you want to use. This will help you think about the different multiplex reactions that you can use.
  3.  Order the Qiagen Multiplex Plus Kit– this will streamline everything! I accidentally ordered the regular multiplex kit.. which is an older version and slower- so I had to stick with it once I got started. However- the Plus version enables you to use a faster PCR protocol! This kit is the same thing as their old ‘Type-it kit”, just better.
  4. Order a set of just the universal tails without the fluorophores attached first, in addition to the forward primers with the universal tails, and the reverse primers with the pig tails. This will let you try running all the multiplex reactions out and test them on a gel to make sure you have everything working before you waste your precious flourophores which are expensive. This is also a cheap thing- 4x $6 max.. 24$ to try out a bunch of stuff before spending the big money is well worth it! By the way  you will order all of your flourophores (reporter dyes) attached to the universal tail or the forward primer (the latter is a more expensive technique) from Thermo Fisher Scientific as they have patents on all of them except for FAM, which you can buy cheaper from Sigma or other companies (IDTdna, Elim BioPharm…etc)
  5. When you have everything working (bands are where you expect them to be, and no large gaps between bands which indicate that the msat markers are targeting multiple regions of the DNA sequence rather than one region)- Then order black, sterile micro centrifuge tubes (this link is just an example- but any sterile brand will work)- this will be what you make all your fluorescent primer ‘party’ mixes in, which will protect the fluorescence from degradation- preventing you from having to order more and save you money and time!
  6. Color code everything -from the tops of your individual fluorescent forward-reverse-fluorescent primer mixes (FRT: Forward-tail +Reverse-pigtail +Fluorophore-Tail) to your tubes with primer party mixes (all four FRT), to your  multiplex pcr reaction mixes and to your excel files for the pcr-plates and genotyping plates.
  7. Be organized when you pipette things onto your pcr-plate, and into your genotyping plate, see below image on how I organize pipetting into a pcr-plate, with the samples in the large tubes being moved in the order that I add things to the plate. I close the lid of each tube and move it to the upper tray after I add it to the plate so I don’t lose my place. I also use my tips from the box in order so that I can look at my tip-box as well to see where I should be. IMG_7752
  8. Talk to the genotyping facility about if they permit a ‘troubleshooting’ run (free-of-charge) so that you can test the amount of final pcr product to add into each genotyping well. I used 0.5 ul pcr product (that had amplified products of four markers) with 11 ul genotyping mix (Liz size standard in Hi-di formamide)

Random Tips/Interesting findings:

  • The Liz size-standard has strict “keep in dark and don’t-freeze’ instructions when it is shipped to you. Be sure to put it in the fridge at 4 degrees and NOT in the freezer! Unfortunately it comes in a styrofoam box w/ ice-packs with no-outside labels instructing to not freeze…and so the mailing department or your lab technician might accidentally put it in the freezer (speaking from personal experience…)- The company (Thermo Fisher Scientific) was very nice in shipping me a replacement because of this issue.  However… due to this occurrence I had the opportunity to test whether or not the Liz -size standard would still work when frozen for 6 hours, and when frozen for 24 hours….Results: The size standard still works great when frozen for 6 hours and for 24 hours ! With that said.. obviously don’t purposely test this, but if it is accidentally frozen- chances are you are still ok!
  • I also found that the Liz size-standard works great and is consistent for at least 2 months beyond its written expiry date…
    • For both of these tests I had unexpired and unfrozen Liz-size standard to compare these tests to. No p-value available.. just my experience 😉
  • As for the genotyping- I used less Liz Size Standard than recommended (and so did everyone that I talked to). My specific reaction mixes were the following: 0.5 ul PCR product + 0.5 ul Liz Size standard, and 10.5 ul Hi-Di Formamide per genotyping reaction well. I know a lot of folks that use 0.5 ul PCR product + 0.2 ul Liz-size standard+ 9.3 ul Hi-Di Formamide.. and they have great success as well. I tried the latter mix  and it worked, but because my pcr products had such high fluorescence (and I was over troubleshooting my primer fluorophore mix concentrations- I decided to  instead increase my size-standard so that I could better separate the noise from the signals). My genotyping mix and final primer ‘party’ mixes result in the initial genotype peaks image of this blog- so you can see what I mean by high sample peaks compared to the size-standard.
  • As for the Hi-Di Formamide- I noticed that this does not have good results when you leave it in the fridge overnight and try to use it the next day for a second genotyping plate. However- I had good results with using Hi-Di Formamide that underwent 1-3 freeze-thaw cycles. Thus my advice is that you never leave it at room temperature or in the fridge if you have extra, but also to avoid too many freeze-thaws.
  • Additionally- Im sure you will find this out- but NEVER freeze your pcr-products after the multiplex pcr with the flourescent markers, or after you add the liz-size standard and hi-di formamide. If you can’t get your samples to the genotyping facility right away, then be sure to do a quick denature (95 degrees for 5 min) post combining the pcr product w/ the Liz size standard and hi-di formamide and then just keep in the fridge (and in the dark!) until the next day or two.

References

Blacket, M.J., Robin, C., Good, R.T., Lee, S.F., Miller, A.D. 2012. Universal primers for fluorescent labelling of PCR fragments–an efficient and cost-effective approach to genotyping by fluorescence. Molecular ecology resources 12, 456-463.

Culley, T.M., Stamper, T.I., Stokes, R.L., Brzyski, J.R., Hardiman, N.A., Klooster, M.R., Merritt, B.J. 2013. An efficient technique for primer development and application that integrates fluorescent labeling and multiplex PCR. Appl Plant Sci 1.

 

 

 

 

 

 

 

 

How-to use microsatellites for population genetics, Part I: Study Design, DNA extraction, Microsatellite Marker Design/Outsourcing

So… you want to use microsatellite markers to assess the genetic variation and population structure of your focal study organism? Well if you are anything like me two years ago.. then you have no idea where to start. Otherwise- congratulations if you are already an expert- in which case you probably don’t need to read on 🙂

SeeHearSpeak
“See No Weevil, Hear No Weevil, Speak No Weevil”                                                                          Illustration by Jacki Whisenant, contracted by Julie Hopper. Copyright 2017.

Two years ago, I was just like you (and these weevils above), and felt a bit overwhelmed and lost in undertaking the large task of designing microsatellite markers and genotyping these markers for the two weevils species (Neochetina bruchi and N. eichhorniae) that I have discussed in previous posts. 

Very briefly to recap on my work:  these two weevil species are used all over the world for the biological control of the invasive water hyacinth, including the Sacramento-San Joaquin River Delta, California. They have had variable success, with notable reduction of biomass and cover of water hyacinth in warmer climates compared to more temperate climates such as the Delta. Although temperature plays a large role in their success, I am also investigating the role of genetic variation and particularly whether there is lower genetic diversity and heterozygosity in the Delta compared to the native origin of these weevils (Uruguay and Argentina).

In Part I- (this blog), I will detail the how-to’s of sampling design and strategy, and the development of (or outsourcing) microsatellite markers.

In Part II- (next blog) I will discuss how to make your final microsatellite marker selections, and the workflow of multiplex PCR and genotyping.

In Part III- (come back in a month!) I will detail how to analyze the data with various R-packages and other computer programs, and how to format the data files correctly for these programs.

On this note, please research your study system thoroughly, as every organism is different and may require different sampling strategies and methods than I detail here for two diploid beetle species (Insecta). Additionally.. my overview below on Part I- is very brief and I definitely skip small steps to be succinct. Also my suggestions are not the only way to do things and below this blog, I post links to several other great resources. Lastly- This work is currently in prep for publication and I will post an update again after publication.


Part I: 

 

Figure from: Grunwald et al. 2017, Phytopathology
Figure from: Grunwald, N.J., Everhart, S.E., Knaus, B.J., Kamvar, Z.N. 2017. Best Practices for Population Genetic Analyses. Phytopathology 107, 1000-1010.

Sampling Design and Strategy:

First before you start sampling or ordering primers- make sure that you have a solid study question with a testable hypothesis, and a good study framework.

Next: all of the power in your genetic analyses (aka, accuracy and ability to detect differentiation among populations, etc.) depend on: 1) your sample quality (aka DNA quality), the number of samples (replicates) per treatment or location, 2) the number of high quality microsatellite markers (e.g.quality relating to two important characteristics: markers are polymorphic -having 2 or more alleles per locus-with more being better, and the markers lack true null alleles), 3) the robustness of your PCR  – whether the PCR conditions are truly suitable for your markers, and whether they can result in reproducible data, 4) the assumptions of the data and 5) the choice of statistical tests and whether the tests are truly suitable for the data.

I will cover the latter (regarding statistical tests) in a future blog, but for today I would like to focus on the ideal # of samples and the # of polymorphic markers. There has been debate about how many samples and how many markers are necessary for robust studies, and if you study an endangered species -sometimes you just have to work with what you got!

In a perfect world– you will want to make up for what you lack in samples with microsatellite markers (loci) and vice versa. So if you have a lower end of replicates, then you will want a higher number of microsatellite markers (# of loci, and more important is to have polymorphic loci with 2 or more alleles/locus) to test for each individual (replicate), and again vice-versa. There are a couple great papers that discuss sampling strategies and study design that you should definitely check out, particularly the one noted in the figure above (Grunwald et al. 2017), as well as Hale et al. 2012 which states that 25-30 individuals per population should be sufficient to accurately estimate allele frequencies given population (with some caveats). Caveats being that obviously, 25-30 individuals per population would likely NOT be enough if you only have four microsatellite markers, particularly if these markers are not polymorphic or very variable (variability referencing to the # of alleles per locus- the more the better!).. so keep this in mind. In general, with that many samples- 10-15 polymorphic markers should be fine (although the more the better), but again this depends on your study question and study system. Also, more samples might be necessary if you are interested in population differentiation (population genetic structure). In fact, in a landscape genetics study, Landguth et al. 2012 demonstrated that increasing the number of loci (and particularly having more variable loci) is more likely to increase the power of population genetic inferences compared to increasing the number of individuals.

You can also test your samples with genotype accumulation curves to see if you have captured the majority of genetic variation (I used the poppr package in R for this and will discuss more on poppr and its primer in Part III of this blog series).

With that said.. If I would have known 1 year ago what I know now…. I would have asked for folks around the world to collect more weevils for me, and I would have extracted more DNA!  Just remember.. not all of your DNA extractions are going to end up working out..due to various human error and/or preservation issues. Thus its always good to add at least 10-20 more samples than you think you need!

map_with_labels_pop_gen
Sampling locations of Neochetina bruchi and N. eichhorniae individuals that I used for the focal population genetics study (Hopper et al. In Prep). Thanks to all those who sent me weevils!

Designing or Outsourcing Microsatellite Marker Design: 

  • Marker Outsource Options: I want to first be upfront in that I actually ended up outsourcing this component of my study as I was going through a tough time and taking care of my dad who had metastatic cancer via at-home hospice care in Columbus, Ohio for two months. Needless to say- I was working remotely then, which made the decision to outsource this part of the lab work an easy decision. I researched a lot of outsource options and in the end I went with the cheaper and most recommended option by several colleagues- the Savannah River Ecology Lab at the University of Georgia. In the end I have mixed opinions on their work and please email me if you would like more info and I will detail the ups and downs.
  • Brief Workflow for designing microsatellite markers: 
    1. First! Check the literature to make sure microsatellite markers have not already been developed for your species or a sister species (the latter of which will sometimes work). Using previously developed markers is obviously the easiest and cheapest route!
    2. If the markers have not already been developed: Obtain high quality and high molecular weight DNA Extractions. I love doing 5% Chelex DNA extractions, but the resulting DNA can be full of PCR inhibitors- so I always use the second half of the DNAeasy kit to purify and clean up my DNA samples. You can also buy replacement spin columns for these kits way cheaper from Epoch Life Science. Then quantify them on a nano-drop or a similar DNA quantification instrument and additionally run them on a gel to make sure that you have ≥100 uL of ≥50 ng/uL of >10kb DNA per sample.
    3. Send to a sequencing facility (Illumina with paired ends >150bp preferred)
    4. Clean up sequences/fix Errors and Run a program called “Pal_finder”, or use a similar program. Pal_finder can analyze 454 or paired-end Illumina sequences ( ~150bp from each end).  This program sends possible primers to Primer3 for primer design and searches for how often each primer and primer pair occur.

    5. Filter the resulting data set by only including: a) sequences for which primers can be designed (e.g. enough flanking sequence) and b) primer pairs that occurred 1-3 times. Then, sort by motif length (di, tri, tetra, etc.) to quickly find tri or tetra nucleotide repeats and look to see if the motif was found in both directions of the sequence (which can be bad as they typically end up being smaller PCR products, but this depends on your goals). Finally, order a bunch of the primers that look promising-say 48 primer pairs to start, and test them out on a subset of 24 individuals, with an equal distribution of these individuals across all your study locations, or select individuals that you think will have a lot of variation. See Initial PCR testing in the next Blog. 

To be continued…

References

Grunwald, N.J., Everhart, S.E., Knaus, B.J., Kamvar, Z.N. 2017. Best Practices for Population Genetic Analyses. Phytopathology 107, 1000-1010.

Hale, M.L., Burg, T.M., Steeves, T.E. 2012. Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PloS one 7, e45170.

Landguth, E.L., Fedy, B.C., Oyler-McCance, S.J., Garey, A.L., Emel, S.L., Mumma, M., Wagner, H.H., Fortin, M.-J., Cushman, S.A. 2012. Effects of sample size, number of markers, and allelic richness on the detection of spatial genetic pattern. Molecular ecology resources 12, 276-284.

Helpful Resources on Getting Started for Part I

Lecture on Intro to Microsatellites