Population genetics models of local ancestry

Migrations have played an important role in shaping the genetic diversity of human populations. Understanding genomic data thus requires careful modeling of historical gene flow. Here we consider the effect of relatively recent population structure and gene flow, and interpret genomes of individuals that have ancestry from multiple source populations as mosaics of segments originating from each population. We propose general and tractable models for describing the evolution of these patterns of local ancestry and their impact on genetic diversity. We focus on the length distribution of continuous ancestry tracts, and the variance in total ancestry proportions among individuals. The proposed models offer improved agreement with Wright-Fisher simulation data when compared to state-of-the art models, and can be used to infer various demographic parameters in gene flow models. Considering HapMap African-American (ASW) data, we find that a model with two distinct phases of `European' gene flow significantly improves the modeling of both tract lengths and ancestry variances.Comment: 25 pages with 7 figures; Genetics: Published online before print April 4, 201

The Landscape of Realized Homologous Recombination in Pathogenic Bacteria

Recombination enhances the adaptive potential of organisms by allowing genetic variants to be tested on multiple genomic backgrounds. Its distribution in the genome can provide insight into the evolutionary forces that underlie traits, such as the emergence of pathogenicity. Here, we examined landscapes of realized homologous recombination of 500 genomes from ten bacterial species and found all species have “hot” regions with elevated rates relative to the genome average. We examined the size, gene content, and chromosomal features associated with these regions and the correlations between closely related species. The recombination landscape is variable and evolves rapidly. For example in Salmonella, only short regions of around 1 kb in length are hot whereas in the closely related species Escherichia coli, some hot regions exceed 100 kb, spanning many genes. Only Streptococcus pyogenes shows evidence for the positive correlation between GC content and recombination that has been reported for several eukaryotes. Genes with function related to the cell surface/membrane are often found in recombination hot regions but E. coli is the only species where genes annotated as “virulence associated” are consistently hotter. There is also evidence that some genes with “housekeeping” functions tend to be overrepresented in cold regions. For example, ribosomal proteins showed low recombination in all of the species. Among specific genes, transferrin-binding proteins are recombination hot in all three of the species in which they were found, and are subject to interspecies recombination

Environmental factors influence both abundance and genetic diversity in a widespread bird species.

Genetic diversity is one of the key evolutionary variables that correlate with population size, being of critical importance for population viability and the persistence of species. Genetic diversity can also have important ecological consequences within populations, and in turn, ecological factors may drive patterns of genetic diversity. However, the relationship between the genetic diversity of a population and how this interacts with ecological processes has so far only been investigated in a few studies. Here, we investigate the link between ecological factors, local population size, and allelic diversity, using a field study of a common bird species, the house sparrow (Passer domesticus). We studied sparrows outside the breeding season in a confined small valley dominated by dispersed farms and small-scale agriculture in southern France. Population surveys at 36 locations revealed that sparrows were more abundant in locations with high food availability. We then captured and genotyped 891 house sparrows at 10 microsatellite loci from a subset of these locations (N = 12). Population genetic analyses revealed weak genetic structure, where each locality represented a distinct substructure within the study area. We found that food availability was the main factor among others tested to influence the genetic structure between locations. These results suggest that ecological factors can have strong impacts on both population size per se and intrapopulation genetic variation even at a small scale. On a more general level, our data indicate that a patchy environment and low dispersal rate can result in fine-scale patterns of genetic diversity. Given the importance of genetic diversity for population viability, combining ecological and genetic data can help to identify factors limiting population size and determine the conservation potential of populations

PLoS Genet.

Our understanding of basic evolutionary processes in bacteria is still very limited. For example, multiple recent dating estimates are based on a universal inter-species molecular clock rate, but that rate was calibrated using estimates of geological dates that are no longer accepted. We therefore estimated the short-term rates of mutation and recombination in Helicobacter pylori by sequencing an average of 39,300 bp in 78 gene fragments from 97 isolates. These isolates included 34 pairs of sequential samples, which were sampled at intervals of 0.25 to 10.2 years. They also included single isolates from 29 individuals (average age: 45 years) from 10 families. The accumulation of sequence diversity increased with time of separation in a clock-like manner in the sequential isolates. We used Approximate Bayesian Computation to estimate the rates of mutation, recombination, mean length of recombination tracts, and average diversity in those tracts. The estimates indicate that the short-term mutation rate is 1.4x10(-6) (serial isolates) to 4.5x10(-6) (family isolates) per nucleotide per year and that three times as many substitutions are introduced by recombination as by mutation. The long-term mutation rate over millennia is 5-17-fold lower, partly due to the removal of non-synonymous mutations due to purifying selection. Comparisons with the recent literature show that short-term mutation rates vary dramatically in different bacterial species and can span a range of several orders of magnitude

Stratified dispersal and increasing genetic variation during the invasion of Central Europe by the western corn rootworm, Diabrotica virgifera virgifera

Invasive species provide opportunities for investigating evolutionary aspects of colonization processes, including initial foundations of populations and geographic expansion. Using microsatellite markers and historical information, we characterized the genetic patterns of the invasion of the western corn rootworm (WCR), a pest of corn crops, in its largest area of expansion in Europe: Central and South-Eastern (CSE) Europe. We found that the invaded area probably corresponds to a single expanding population resulting from a single introduction of WCR and that gene flow is geographically limited within the population. In contrast to what is expected in classical colonization processes, an increase in genetic variation was observed from the center to the edge of the outbreak. Control measures against WCR at the center of the outbreak may have decreased effective population size in this area which could explain this observed pattern of genetic variation. We also found that small remote outbreaks in southern Germany and north-eastern Italy most likely originated from long-distance dispersal events from CSE Europe. We conclude that the large European outbreak is expanding by stratified dispersal, involving both continuous diffusion and discontinuous long-distance dispersal. This latter mode of dispersal may accelerate the expansion of WCR in Europe in the future

Bayesian modeling of recombination events in bacterial populations

Background: We consider the discovery of recombinant segments jointly with their origins within multilocus DNA sequences from bacteria representing heterogeneous populations of fairly closely related species. The currently available methods for recombination detection capable of probabilistic characterization of uncertainty have a limited applicability in practice as the number of strains in a data set increases. Results: We introduce a Bayesian spatial structural model representing the continuum of origins over sites within the observed sequences, including a probabilistic characterization of uncertainty related to the origin of any particular site. To enable a statistically accurate and practically feasible approach to the analysis of large-scale data sets representing a single genus, we have developed a novel software tool (BRAT, Bayesian Recombination Tracker) implementing the model and the corresponding learning algorithm, which is capable of identifying the posterior optimal structure and to estimate the marginal posterior probabilities of putative origins over the sites. Conclusion: A multitude of challenging simulation scenarios and an analysis of real data from seven housekeeping genes of 120 strains of genus Burkholderia are used to illustrate the possibilities offered by our approach. The software is freely available for download at URL http://web.abo.fi/fak/ mnf//mate/jc/software/brat.html

Introgression of exotic <i>Cervus</i> (<i>nippon</i> and <i>canadensis</i>) into red deer (<i>Cervus elaphus</i>) populations in Scotland and the English Lake District

Since the mid-19th century, multiple introductions of Japanese sika deer (Cervus nippon nippon) and North American wapiti (C. canadensis) have taken place in the British Isles. While wapiti have generally been unsuccessful, sika have been very successful, especially in Scotland where they now overlap at least 40% of the range of native red deer (C. elaphus). Hybridization between these two species and red deer has been demonstrated in captivity and in the wild. Using a panel of 22 microsatellite loci that are highly diagnostic between red deer and sika, and moderately diagnostic between red deer and wapiti, we investigated the extent of introgression between these species in 2,943 deer sampled from around Scotland and from the English Lake District using the Bayesian clustering software STRUCTURE. We also used a diagnostic mitochondrial marker for red deer and sika. Our survey extends previous studies indicating little introgression of wapiti nuclear alleles into red deer, in particular in Northern Scotland, Kintyre, and the Lake District. We found a new area of extensive sika introgression in South Kintyre. In the North Highlands, we show for the first time geographically scattered evidence of past hybridization followed by extensive backcrossing, including one red-like individual with sika introgression, two sika-like individuals with red deer introgression, and six individuals that were apparently pure sika at the nuclear markers assessed but which carried red deer mitochondria. However, there has not been a collapse of assortative mating in this region. Similarly, in the English Lake District red deer, we found only traces of past sika introgression. No sika alleles were detected in the Central Highlands or the Hebridean red deer refugia. We make suggestions for management to prevent further spread of sika alleles into red deer and vice versa.Peer Reviewe

Secondary contact and admixture between independently invading populations of the Western corn rootworm, diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a Bayesian analysis of the population structure and in an approximate Bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed

Novel associations for hypothyroidism include known autoimmune risk loci

Hypothyroidism is the most common thyroid disorder, affecting about 5% of the general population. Here we present the first large genome-wide association study of hypothyroidism, in 2,564 cases and 24,448 controls from the customer base of 23andMe, Inc., a personal genetics company. We identify four genome-wide significant associations, two of which are well known to be involved with a large spectrum of autoimmune diseases: rs6679677 near _PTPN22_ and rs3184504 in _SH2B3_ (p-values 3.5e-13 and 3.0e-11, respectively). We also report associations with rs4915077 near _VAV3_ (p-value 8.3e-11), another gene involved in immune function, and rs965513 near _FOXE1_ (p-value 3.1e-14). Of these, the association with _PTPN22_ confirms a recent small candidate gene study, and _FOXE1_ was previously known to be associated with thyroid-stimulating hormone (TSH) levels. Although _SH2B3_ has been previously linked with a number of autoimmune diseases, this is the first report of its association with thyroid disease. The _VAV3_ association is novel. These results suggest heterogeneity in the genetic etiology of hypothyroidism, implicating genes involved in both autoimmune disorders and thyroid function. Using a genetic risk profile score based on the top association from each of the four genome-wide significant regions in our study, the relative risk between the highest and lowest deciles of genetic risk is 2.1