Ryan N Gutenkunst

PMID: 17140603;Abstract:

We use observed movement tracks of Atlantic bluefin tuna in the Gulf of Maine and mathematical modeling of this movement to identify possible resource patches. We infer bounds on the overall sizes and distribution of such patches, even though they are difficult to quantify by direct observation in situ. To do so, we segment individual fish tracks into intervals of distinct motion types based on the ratio of net displacement to length of track (Δ D / Δ L) over a time window Δ t. To find the best segmentation, we optimize the fit of a random-walk movement model to each motion type. We compare results from two distinct movement models: biased turning and biased speed, to check the model-dependence of our inferences, and find that uncertainty in choice of movement model dominates the uncertainties of our conclusions. We find that our data are best described using two motion types: "localized" (Δ D / Δ L small) and "long-ranged" (Δ D / Δ L large). The biased turning model leads to significantly better resolution of localized movement intervals than the biased speed model. We hypothesize that localized movement corresponds to exploitation of resource patches. Comparison with visual behavior observations made during tracking suggests that many inferred intervals of localized motion do indeed correspond to feeding activity. From our analysis, we estimate that, on average, bluefin tuna in the Gulf of Maine encounter a resource patch every 2 h, that those patches have an average radius of 0.7-1.2 km, and that, overall, there are at most 5-9 such patches per 100 km² in the region studied. © 2006 Elsevier Ltd. All rights reserved.

The distribution of mutational effects on fitness is central to evolutionary genetics. Typical univariate distributions, however, cannot model the effects of multiple mutations at the same site, so we introduce a model in which mutations at the same site have correlated fitness effects. To infer the strength of that correlation, we developed a diffusion approximation to the triallelic frequency spectrum, which we applied to data from Drosophila melanogaster We found a moderate positive correlation between the fitness effects of nonsynonymous mutations at the same codon, suggesting that both mutation identity and location are important for determining fitness effects in proteins. We validated our approach by comparing it to biochemical mutational scanning experiments, finding strong quantitative agreement, even between different organisms. We also found that the correlation of mutational fitness effects was not affected by protein solvent exposure or structural disorder. Together, our results suggest that the correlation of fitness effects at the same site is a previously overlooked yet fundamental property of protein evolution.

Tumor genome sequencing offers great promise for guiding research and therapy, but spurious variant calls can arise from multiple sources. Mouse contamination can generate many spurious calls when sequencing patient-derived xenografts (PDXs). Paralogous genome sequences can also generate spurious calls when sequencing any tumor. We developed a BLAST-based algorithm, MAPEX, to identify and filter out spurious calls from both these sources.

PMID: 21730125;PMCID: PMC3142009;Abstract:

High-throughput sequencing technology enables population-level surveys of human genomic variation. Here, we examine the joint allele frequency distributions across continental human populations and present an approach for combining complementary aspects of whole-genome, low-coverage data and targeted highcoverage data. We apply this approach to data generated by the pilot phase of the Thousand Genomes Project, including wholegenome 2-4x coverage data for 179 samples from HapMap European, Asian, and African panels as well as high-coverage target sequencing of the exons of 800 genes from 697 individuals in seven populations. We use the site frequency spectra obtained from these data to infer demographic parameters for an Out-of-Africa model for populations of African, European, and Asian descent and to predict, by a jackknife-based approach, the amount of genetic diversity that will be discovered as sample sizes are increased. We predict that the number of discovered nonsynonymous coding variants will reach 100,000 in each population after ∼1,000 sequenced chromosomes per population, whereas ∼2,500 chromosomes will be needed for the same number of synonymous variants. Beyond this point, the number of segregating sites in the European and Asian panel populations is expected to overcome that of the African panel because of faster recent population growth. Overall, we find that the majority of human genomic variable sites are rare and exhibit little sharing among diverged populations. Our results emphasize that replication of disease association for specific rare genetic variants across diverged populations must overcome both reduced statistical power because of rarity and higher population divergence.

Rice is a staple crop that has undergone substantial phenotypic and physiological changes during domestication. Here we resequenced the genomes of 40 cultivated accessions selected from the major groups of rice and 10 accessions of their wild progenitors (Oryza rufipogon and Oryza nivara) to >15 × raw data coverage. We investigated genome-wide variation patterns in rice and obtained 6.5 million high-quality single nucleotide polymorphisms (SNPs) after excluding sites with missing data in any accession. Using these population SNP data, we identified thousands of genes with significantly lower diversity in cultivated but not wild rice, which represent candidate regions selected during domestication. Some of these variants are associated with important biological features, whereas others have yet to be functionally characterized. The molecular markers we have identified should be valuable for breeding and for identifying agronomically important genes in rice.