Honors & Awards

  • Gerald J. Lieberman Fellowship, VPGE (2016-2017)
  • Graduate Research Fellowship, National Science Foundation (2013-2016)
  • Thouron Award, Thouron Award (2012-2013)

Education & Certifications

  • BA, University of Pennsylvania, Mathematics (2012)
  • MSc (Res), University of Oxford, Statistics (2013)

Stanford Advisors

All Publications

  • More effective drugs lead to harder selective sweeps in the evolution of drug resistance in HIV-1. eLife Feder, A. F., Rhee, S., Holmes, S. P., Shafer, R. W., Petrov, D. A., Pennings, P. S. 2016; 5


    In the early days of HIV treatment, drug resistance occurred rapidly and predictably in all patients, but under modern treatments, resistance arises slowly, if at all. The probability of resistance should be controlled by the rate of generation of resistance mutations. If many adaptive mutations arise simultaneously, then adaptation proceeds by soft selective sweeps in which multiple adaptive mutations spread concomitantly, but if adaptive mutations occur rarely in the population, then a single adaptive mutation should spread alone in a hard selective sweep. Here, we use 6717 HIV-1 consensus sequences from patients treated with first-line therapies between 1989 and 2013 to confirm that the transition from fast to slow evolution of drug resistance was indeed accompanied with the expected transition from soft to hard selective sweeps. This suggests more generally that evolution proceeds via hard sweeps if resistance is unlikely and via soft sweeps if it is likely.

    View details for DOI 10.7554/eLife.10670

    View details for PubMedID 26882502

  • The population genetics of drug resistance evolution in natural populations of viral, bacterial and eukaryotic pathogens. Molecular ecology Wilson, B. A., Garud, N. R., Feder, A. F., Assaf, Z. J., Pennings, P. S. 2016; 25 (1): 42-66


    Drug resistance is a costly consequence of pathogen evolution and a major concern in public health. In this review we show how population genetics can be used to study the evolution of drug resistance and also how drug resistance evolution is informative as an evolutionary model system. We highlight five examples from diverse organisms with particular focus on: 1) identifying drug resistance loci in the malaria parasite Plasmodium falciparum using the genomic signatures of selective sweeps, 2) determining the role of epistasis in drug resistance evolution in influenza, 3) quantifying the role of standing genetic variation in the evolution of drug resistance in HIV, 4) using drug resistance mutations to study clonal interference dynamics in tuberculosis, and 5) analyzing the population structure of the core and accessory genome of Staphylococcus aureus to understand the spread of methicillin resistance. Throughout this review we discuss the uses of sequence data and population genetic theory in studying the evolution of drug resistance. This article is protected by copyright. All rights reserved.

    View details for DOI 10.1111/mec.13474

    View details for PubMedID 26578204

  • Identifying Signatures of Selection in Genetic Time Series. Genetics Feder, A. F., Kryazhimskiy, S. A., Plotkin, J. B. 2013


    Both genetic drift and natural selection cause the frequencies of alleles in a population to vary over time. Discriminating between these two evolutionary forces, based on a time series of samples from a population, remains an outstanding problem with increasing relevance to modern data sets. Even in the idealized situation when the sampled locus is independent of all other loci this problem is difficult to solve, especially when the size of the population from which the samples are drawn is unknown. A standard χ(2)-based likelihood ratio test was previously proposed to address this problem. Here we show that the χ(2) test of selection substantially underestimates the probability of Type I error, leading to more false positives than indicated by its P-value, especially at stringent P-values. We introduce two methods to correct this bias. The empirical likelihood ratio test (ELRT) rejects neutrality when the likelihood ratio statistic falls in the tail of the empirical distribution obtained under the most likely neutral population size. The frequency increment test (FIT) rejects neutrality if the distribution of normalized allele frequency increments exhibits a mean that deviates significantly from zero. We characterize the statistical power of these two tests for selection, and we apply them to three experimental data sets. We demonstrate that both ELRT and FIT have power to detect selection in practical parameter regimes, such as those encountered in microbial evolution experiments. Our analysis applies to a single diallelic locus, assumed independent of all other loci, which is most relevant to full-genome selection scans in sexual organisms, and also to evolution experiments in asexual organisms as long as clonal interference is weak. Different techniques will be required to detect selection in time series of co-segregating linked loci.

    View details for DOI 10.1534/genetics.113.158220

    View details for PubMedID 24318534

  • LDx: Estimation of Linkage Disequilibrium from High-Throughput Pooled Resequencing Data PLOS ONE Feder, A. F., Petrov, D. A., Bergland, A. O. 2012; 7 (11)


    High-throughput pooled resequencing offers significant potential for whole genome population sequencing. However, its main drawback is the loss of haplotype information. In order to regain some of this information, we present LDx, a computational tool for estimating linkage disequilibrium (LD) from pooled resequencing data. LDx uses an approximate maximum likelihood approach to estimate LD (r(2)) between pairs of SNPs that can be observed within and among single reads. LDx also reports r(2) estimates derived solely from observed genotype counts. We demonstrate that the LDx estimates are highly correlated with r(2) estimated from individually resequenced strains. We discuss the performance of LDx using more stringent quality conditions and infer via simulation the degree to which performance can improve based on read depth. Finally we demonstrate two possible uses of LDx with real and simulated pooled resequencing data. First, we use LDx to infer genomewide patterns of decay of LD with physical distance in D. melanogaster population resequencing data. Second, we demonstrate that r(2) estimates from LDx are capable of distinguishing alternative demographic models representing plausible demographic histories of D. melanogaster.

    View details for DOI 10.1371/journal.pone.0048588

    View details for Web of Science ID 000312272600012

    View details for PubMedID 23152785

  • Natural Selection Affects Multiple Aspects of Genetic Variation at Putatively Neutral Sites across the Human Genome PLOS GENETICS Lohmueller, K. E., Albrechtsen, A., Li, Y., Kim, S. Y., Korneliussen, T., Vinckenbosch, N., Tian, G., Huerta-Sanchez, E., Feder, A. F., Grarup, N., Jorgensen, T., Jiang, T., Witte, D. R., Sandbaek, A., Hellmann, I., Lauritzen, T., Hansen, T., Pedersen, O., Wang, J., Nielsen, R. 2011; 7 (10)


    A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.

    View details for DOI 10.1371/journal.pgen.1002326

    View details for Web of Science ID 000296665400027

    View details for PubMedID 22022285