An approximate Markov model for the wright-fisher diffusion and its application to time series data

The joint and accurate inference of selection and demography from genetic data is considered a particularly challenging question in population genetics, since both process may lead to very similar patterns of genetic diversity. However, additional information for disentangling these effects may be obtained by observing changes in allele frequencies over multiple time points. Such data is common in experimental evolution studies, as well as in the comparison of ancient and contemporary samples. Leveraging this information, however, has been computationally challenging, particularly when considering multi-locus data sets. To overcome these issues, we introduce a novel, discrete approximation for diffusion processes, termed mean transition time approximation, which preserves the long-term behavior of the underlying continuous diffusion process. We then derive this approximation for the particular case of inferring selection and demography from time series data under the classic Wright- Fisher model and demonstrate that our approximation is well suited to describe allele trajectories through time, even when only a few states are used. We then develop a Bayesian inference approach to jointly infer the population size and locus-specific selection coefficients with high accuracy, and further extend this model to also infer the rates of sequencing errors and mutations. We finally apply our approach to recent experimental data on the evolution of drug resistance in Influenza virus, identifying likely targets of selection and finding evidence for much larger viral population sizes than previously reported

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

9 pages, 2 figures, 4 additional files.[Background] It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.[Results] Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates.[Conclusion] The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.This research was supported by "Fundación Genoma España" (proyectos piloto CEGEN 2004–2005), Dirección General de Investigación, Ministerio de Educación y Ciencia of Spain (grants BFU2005-00243, BFU2006-01235, BFU2006-15413-CO2-01, SEJ2006-13537) and Direcció General de Recerca, Generalitat de Catalunya (2005SGR00608). SNP genotyping services were provided by the Spanish "Centro Nacional de Genotipado"Peer reviewe

Identifying the favored mutation in a positive selective sweep.

Most approaches that capture signatures of selective sweeps in population genomics data do not identify the specific mutation favored by selection. We present iSAFE (for "integrated selection of allele favored by evolution"), a method that enables researchers to accurately pinpoint the favored mutation in a large region (∼5 Mbp) by using a statistic derived solely from population genetics signals. iSAFE does not require knowledge of demography, the phenotype under selection, or functional annotations of mutations

Secretor Genotype (FUT2 gene) Is Strongly Associated with the Composition of Bifidobacteria in the Human Intestine

Intestinal microbiota plays an important role in human health, and its composition is determined by several factors, such as diet and host genotype. However, thus far it has remained unknown which host genes are determinants for the microbiota composition. We studied the diversity and abundance of dominant bacteria and bifidobacteria from the faecal samples of 71 healthy individuals. In this cohort, 14 were non-secretor individuals and the remainders were secretors. The secretor status is defined by the expression of the ABH and Lewis histo-blood group antigens in the intestinal mucus and other secretions. It is determined by fucosyltransferase 2 enzyme, encoded by the FUT2 gene. Non-functional enzyme resulting from a nonsense mutation in the FUT2 gene leads to the non-secretor phenotype. PCR-DGGE and qPCR methods were applied for the intestinal microbiota analysis. Principal component analysis of bifidobacterial DGGE profiles showed that the samples of non-secretor individuals formed a separate cluster within the secretor samples. Moreover, bifidobacterial diversity (p<0.0001), richness (p<0.0003), and abundance (p<0.05) were significantly reduced in the samples from the non-secretor individuals as compared with those from the secretor individuals. The non-secretor individuals lacked, or were rarely colonized by, several genotypes related to B. bifidum, B. adolescentis and B. catenulatum/pseudocatenulatum. In contrast to bifidobacteria, several bacterial genotypes were more common and the richness (p<0.04) of dominant bacteria as detected by PCR-DGGE was higher in the non-secretor individuals than in the secretor individuals. We showed that the diversity and composition of the human bifidobacterial population is strongly associated with the histo-blood group ABH secretor/non-secretor status, which consequently appears to be one of the host genetic determinants for the composition of the intestinal microbiota. This association can be explained by the difference between the secretor and non-secretor individuals in their expression of ABH and Lewis glycan epitopes in the mucosa

Small-scale spatial variation in infection risk shapes the evolution of a Borrelia resistance gene in wild rodents

This is the author accepted manuscript. The final version is available from Wiley via the DOI in this record.Spatial variation in pathogen-mediated selection is predicted to influence the evolutionary trajectory of host populations and lead to spatial variation in their immunogenetic composition. However, to date few studies have been able to directly link small-scale spatial variation in infection risk to host immune gene evolution in natural, non-human populations. Here we use a natural rodent-Borrelia system to test for associations between landscape-level spatial variation in Borrelia infection risk along replicated elevational gradients in the Swiss Alps and Toll-like receptor 2 (TLR2) evolution, a candidate gene for Borrelia resistance, across bank vole (Myodes glareolus) populations. We found that Borrelia infection risk (i.e. the product of Borrelia prevalence in questing ticks and the average tick load of voles at a sampling site) was spatially variable and significantly negatively associated with elevation. Across sampling sites, Borrelia prevalence in bank voles was significantly positively associated with Borrelia infection risk along the elevational clines. We observed a significant association between naturally occurring TLR2 polymorphisms in hosts and their Borrelia infection status. The TLR2 variant associated with a reduced likelihood of Borrelia infection was most common in rodent populations at lower elevations that face a high Borrelia infection risk, and its frequency changed in accordance with the change in Borrelia infection risk along the elevational clines. These results suggest that small-scale spatial variation in parasite-mediated selection affects the immunogenetic composition of natural host populations, providing a striking example that the microbial environment shapes the evolution of the host’s immune system in the wild.The study was financially supported by the Swiss National Science Foundation (PP00P3_128386 and PP00P3_157455), the University of Zurich Research Priority Program ‘Evolution in Action: from Genomes to Ecosystems’, the Faculty of Science of the University of Zurich, the Baugarten Stiftung, the Stiftung für wissenschaftliche Forschung an der Universität Zürich (STWF 17_027) and the Georges und Antoine Claraz-Schenkung

The genomic basis of eco-evolutionary dynamics.

Recent recognition that ecological and evolutionary processes can operate on similar timescales has led to a rapid increase in theoretical and empirical studies on eco-evolutionary dynamics. Progress in the fields of evolutionary biology, genomics and ecology is greatly enhancing our understanding of rapid adaptive processes, the predictability of adaptation and the genetics of ecologically important traits. However, progress in these fields has proceeded largely independently of one another. In an attempt to better integrate these fields, the centre for 'Adaptation to a Changing Environment' organized a conference entitled 'The genomic basis of eco-evolutionary change' and brought together experts in ecological genomics and eco-evolutionary dynamics. In this review, we use the work of the invited speakers to summarize eco-evolutionary dynamics and discuss how they are relevant for understanding and predicting responses to contemporary environmental change. Then, we show how recent advances in genomics are contributing to our understanding of eco-evolutionary dynamics. Finally, we highlight the gaps in our understanding of eco-evolutionary dynamics and recommend future avenues of research in eco-evolutionary dynamics

Different level of population differentiation among human genes

<p>Abstract</p> <p>Background</p> <p>During the colonization of the world, after dispersal out of African, modern humans encountered changeable environments and substantial phenotypic variations that involve diverse behaviors, lifestyles and cultures, were generated among the different modern human populations.</p> <p>Results</p> <p>Here, we study the level of population differentiation among different populations of human genes. Intriguingly, genes involved in osteoblast development were identified as being enriched with higher <it>F</it>_STSNPs, a result consistent with the proposed role of the skeletal system in accounting for variation among human populations. Genes involved in the development of hair follicles, where hair is produced, were also found to have higher levels of population differentiation, consistent with hair morphology being a distinctive trait among human populations. Other genes that showed higher levels of population differentiation include those involved in pigmentation, spermatid, nervous system and organ development, and some metabolic pathways, but few involved with the immune system. Disease-related genes demonstrate excessive SNPs with lower levels of population differentiation, probably due to purifying selection. Surprisingly, we find that Mendelian-disease genes appear to have a significant excessive of SNPs with high levels of population differentiation, possibly because the incidence and susceptibility of these diseases show differences among populations. As expected, microRNA regulated genes show lower levels of population differentiation due to purifying selection.</p> <p>Conclusion</p> <p>Our analysis demonstrates different level of population differentiation among human populations for different gene groups.</p

Geographic differences in allele frequencies of susceptibility SNPs for cardiovascular disease

<p>Abstract</p> <p>Background</p> <p>We hypothesized that the frequencies of risk alleles of SNPs mediating susceptibility to cardiovascular diseases differ among populations of varying geographic origin and that population-specific selection has operated on some of these variants.</p> <p>Methods</p> <p>From the database of genome-wide association studies (GWAS), we selected 36 cardiovascular phenotypes including coronary heart disease, hypertension, and stroke, as well as related quantitative traits (eg, body mass index and plasma lipid levels). We identified 292 SNPs in 270 genes associated with a disease or trait at <it>P </it>< 5 × 10^-8. As part of the Human Genome-Diversity Project (HGDP), 158 (54.1%) of these SNPs have been genotyped in 938 individuals belonging to 52 populations from seven geographic areas. A measure of population differentiation, <it>F</it>_ST, was calculated to quantify differences in risk allele frequencies (RAFs) among populations and geographic areas.</p> <p>Results</p> <p>Large differences in RAFs were noted in populations of Africa, East Asia, America and Oceania, when compared with other geographic regions. The mean global <it>F</it>_ST(0.1042) for 158 SNPs among the populations was not significantly higher than the mean global <it>F</it>_STof 158 autosomal SNPs randomly sampled from the HGDP database. Significantly higher global <it>F</it>_ST(<it>P </it>< 0.05) was noted in eight SNPs, based on an empirical distribution of global <it>F</it>_STof 2036 putatively neutral SNPs. For four of these SNPs, additional evidence of selection was noted based on the integrated Haplotype Score.</p> <p>Conclusion</p> <p>Large differences in RAFs for a set of common SNPs that influence risk of cardiovascular disease were noted between the major world populations. Pairwise comparisons revealed RAF differences for at least eight SNPs that might be due to population-specific selection or demographic factors. These findings are relevant to a better understanding of geographic variation in the prevalence of cardiovascular disease.</p