Estimation of allele frequency and association mapping using next-generation sequencing data

<p>Abstract</p> <p>Background</p> <p>Estimation of allele frequency is of fundamental importance in population genetic analyses and in association mapping. In most studies using next-generation sequencing, a cost effective approach is to use medium or low-coverage data (e.g., < 15<it>X</it>). However, SNP calling and allele frequency estimation in such studies is associated with substantial statistical uncertainty because of varying coverage and high error rates.</p> <p>Results</p> <p>We evaluate a new maximum likelihood method for estimating allele frequencies in low and medium coverage next-generation sequencing data. The method is based on integrating over uncertainty in the data for each individual rather than first calling genotypes. This method can be applied to directly test for associations in case/control studies. We use simulations to compare the likelihood method to methods based on genotype calling, and show that the likelihood method outperforms the genotype calling methods in terms of: (1) accuracy of allele frequency estimation, (2) accuracy of the estimation of the distribution of allele frequencies across neutrally evolving sites, and (3) statistical power in association mapping studies. Using real re-sequencing data from 200 individuals obtained from an exon-capture experiment, we show that the patterns observed in the simulations are also found in real data.</p> <p>Conclusions</p> <p>Overall, our results suggest that association mapping and estimation of allele frequencies should not be based on genotype calling in low to medium coverage data. Furthermore, if genotype calling methods are used, it is usually better not to filter genotypes based on the call confidence score.</p

Rural Cooperatives Magazine, January/February 2013

Features - Creating a safety culture; Here today, here tomorrow; Rocky Ford Renaissance; USDA-supported programs help Mali farmers adapt in hard times; Bio-energy impact, base capital financing among topics at Farmer Co-op Conference; Antitrust challenges facing farmers and their cooperatives; Co-op conversions help bring security to manufactured housing owners; The Co-op Nature of (the Affordable Care Act) CO-OPs; Most Co-op Network members say 2012 was a solid economic yea

Genetic diversity in Cypripedium calceolus (Orchidaceae) with a focus on north-western Europe, as revealed by plastid DNA length polymorphisms

Background and Aims: Cypripedium calceolus, although widespread in Eurasia, is rare in many countries in which it occurs. Population genetics studies with nuclear DNA markers on this species have been hampered by its large nuclear genome size. Plastid DNA markers are used here to gain an understanding of variation within and between populations and of biogeographical patterns. Methods: Thirteen length-variable regions (microsatellites and insertions/deletions) were identified in non-coding plastid DNA. These and a previously identified complex microsatellite in the trnL-trnF intergenic spacer were used to identify plastid DNA haplotypes for European samples, with sampling focused on England, Denmark and Sweden. Key Results: The 13 additional length-variable regions identified were two homopolymer (polyA) repeats in the rps16 intron and a homopolymer (polyA) repeat and ten indels in the accD-psa1 intergenic spacer. In accD-psa1, most of these were in an extremely AT-rich region, and it was not possible to design primers in the flanking regions; therefore, the whole intergenic spacer was sequenced. Together, these new regions and the trnL-trnF complex microsatellite allowed 23 haplotypes to be characterized. Many were found in only one or a few samples (probably due to low sampling density), but some commoner haplotypes were widespread. Most of the genetic variation was found within rather than between populations (83 vs. 18%, respectively). Two haplotypes occurred from the Spanish Pyrenees to Sweden. Conclusions: Plastid DNA data can be used to gain an understanding of patterns of genetic variation and seed-mediated gene flow in orchids. Although these data are less information-rich than those for nuclear DNA, they present a useful option for studying species with large genomes. Here they support the hypothesis of long-distance seed dispersal often proposed for orchids