Sampling distribution for microsatellites amplified by PCR: mean field approximation and its applications to genotyping

Microsatellite evolution: Markov transition functions for a suite of models

This paper takes from the collection of models considered by Whittaker et al. [2003. Likelihood-based estimation of microsatellite mutation rates. Genetics 164, 781-787] derived from direct observation of microsatellite mutation in parent-child pairs and provides analytical expressions for the probability distributions for the change in number of repeats over any given number of generations. The mathematical framework for this analysis is the theory of Markov processes. We find these expressions using two approaches, approximating by circulant matrices and solving a partial differential equation satisfied by the generating function. The impact of the differing choice of models is examined using likelihood estimates for time to most recent common ancestor. The analysis presented here may play a role in elucidating the connections between these two approaches and shows promise in reconciling differences between estimates for mutation rates based on Whittaker's approach and methods based on phylogenetic analyses.

Identification, cloning, and characterization of microsatellite DNA in Euglena gracilis

Microsatellite DNA has been developed into one of the most popular genetic markers. We have identified and cloned microsatellite loci in the genome of a free-living protozoan Euglena gracilis FACHB-848, using the random amplified microsatellites method (RAMS). The digoxigenin-labelled oligonucleotides (CT)10 and (GT)10 served as probes to detect complementary sequences in the randomly amplified polymorphic DNA (RAPD) fingerprints produced by means of Southern blotting. Subsequently, positive RAPD fragments were cloned. From a total of 31 RAPD primer profiles, eight microsatellite loci of E. gracilis were detected and characterized. Further, six sites (i.e. EGMS1, EGMS3, EGMS4, EGMS5, EGMS6, and EGMS7) showed polymorphisms. We found a GT or CT microsatellite every 10.5 kb in the genome of E. gracilis, and similar to animal genomes, the (GT)(n) motif was much more abundant than the (CT)(n) motif. These polymorphic microsatellite DNA will serve as advantageous molecular markers for studying the genetic diversity and molecular ecology of Euglena.

A graphical simulation model of the entire DNA process associated with the analysis of short tandem repeat loci

The use of expert systems to interpret short tandem repeat DNA profiles in forensic, medical and ancient DNA applications is becoming increasingly prevalent as high-throughput analytical systems generate large amounts of data that are time-consuming to process. With special reference to low copy number (LCN) applications, we use a graphical model to simulate stochastic variation associated with the entire DNA process starting with extraction of sample, followed by the processing associated with the preparation of a PCR reaction mixture and PCR itself. Each part of the process is modelled with input efficiency parameters. Then, the key output parameters that define the characteristics of a DNA profile are derived, namely heterozygote balance (Hb) and the probability of allelic drop-out p(D). The model can be used to estimate the unknown efficiency parameters, such as π_extraction. ‘What-if’ scenarios can be used to improve and optimize the entire process, e.g. by increasing the aliquot forwarded to PCR, the improvement expected to a given DNA profile can be reliably predicted. We demonstrate that Hb and drop-out are mainly a function of stochastic effect of pre-PCR molecular selection. Whole genome amplification is unlikely to give any benefit over conventional PCR for LCN.