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

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

Genomic structure and insertion sites of Helicobacter pylori prophages from various geographical origins

We present the full genomic sequences, insertion sites and phylogenetic analysis of 28 prophages found in H. pylori isolates from patients of distinct disease types, ranging from gastritis to gastric cancer, and geographic origins, covering most continents. The gentic diversity of H pylori is known to be influenced by these genomic elements including prophages who’s geneomes range from 22.6 to 33.0 Kbp. There was a high conservation of integration site shared in over 50% of cases with greater than 40% or prophage genomes harbouring insertion sequences (IS). Furthermore prophage genomes present a robust phylogeographic pattern, revealing four distinct clusters: one African, one Asian and two European prophage populations. There was evidence of recombination within the genome of some prophages, which resulted in genome mosaics composed by different populations, which may yield additional H. pylori phenotypes

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

Recent acquisition of Helicobacter pylori by Baka Pygmies

Both anatomically modern humans and the gastric pathogen Helicobacter pylori originated in Africa, and both species have been associated for at least 100,000 years. Seven geographically distinct H. pylori populations exist, three of which are indigenous to Africa: hpAfrica1, hpAfrica2, and hpNEAfrica. The oldest and most divergent population, hpAfrica2, evolved within San hunter-gatherers, who represent one of the deepest branches of the human population tree. Anticipating the presence of ancient H. pylori lineages within all hunter-gatherer populations, we investigated the prevalence and population structure of H. pylori within Baka Pygmies in Cameroon. Gastric biopsies were obtained by esophagogastroduodenoscopy from 77 Baka from two geographically separated populations, and from 101 non-Baka individuals from neighboring agriculturalist populations, and subsequently cultured for H. pylori. Unexpectedly, Baka Pygmies showed a significantly lower H. pylori infection rate (20.8%) than non-Baka (80.2%). We generated multilocus haplotypes for each H. pylori isolate by DNA sequencing, but were not able to identify Baka-specific lineages, and most isolates in our sample were assigned to hpNEAfrica or hpAfrica1. The population hpNEAfrica, a marker for the expansion of the Nilo-Saharan language family, was divided into East African and Central West African subpopulations. Similarly, a new hpAfrica1 subpopulation, identified mainly among Cameroonians, supports eastern and western expansions of Bantu languages. An age-structured transmission model shows that the low H. pylori prevalence among Baka Pygmies is achievable within the timeframe of a few hundred years and suggests that demographic factors such as small population size and unusually low life expectancy can lead to the eradication of H. pylori from individual human populations. The Baka were thus either H. pylori-free or lost their ancient lineages during past demographic fluctuations. Using coalescent simulations and phylogenetic inference, we show that Baka almost certainly acquired their extant H. pylori through secondary contact with their agriculturalist neighbors

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

Phenotypic and Genetic Divergence among Poison Frog Populations in a Mimetic Radiation

The evolution of Müllerian mimicry is, paradoxically, associated with high levels of diversity in color and pattern. In a mimetic radiation, different populations of a species evolve to resemble different models, which can lead to speciation. Yet there are circumstances under which initial selection for divergence under mimicry may be reversed. Here we provide evidence for the evolution of extensive phenotypic divergence in a mimetic radiation in Ranitomeya imitator, the mimic poison frog, in Peru. Analyses of color hue (spectral reflectance) and pattern reveal substantial divergence between morphs. However, we also report that there is a “transition-zone� with mixed phenotypes. Analyses of genetic structure using microsatellite variation reveals some differentiation between populations, but this does not strictly correspond to color pattern divergence. Analyses of gene flow between populations suggest that, while historical levels of gene flow were low, recent levels are high in some cases, including substantial gene flow between some color pattern morphs. We discuss possible explanations for these observations

Contrasting demographic history and gene flow patterns of two mangrove species on either side of the Central American Isthmus

Comparative phylogeography offers a unique opportunity to understand the interplay between past environmental events and life-history traits on diversification of unrelated but co-distributed species. Here, we examined the effects of the quaternary climate fluctuations and palaeomarine currents and present-day marine currents on the extant patterns of genetic diversity in the two most conspicuous mangrove species of the Neotropics. The black (Avicennia germinans, Avicenniaceae) and the red (Rhizophora mangle, Rhizophoraceae) mangroves have similar geographic ranges but are very distantly related and show striking differences on their life-history traits. We sampled 18 Atlantic and 26 Pacific locations for A.germinans (N=292) and R.mangle (N=422). We performed coalescence simulations using microsatellite diversity to test for evidence of population change associated with quaternary climate fluctuations. In addition, we examined whether patterns of genetic variation were consistent with the directions of major marine (historical and present day) currents in the region. Our demographic analysis was grounded within a phylogeographic framework provided by the sequence analysis of two chloroplasts and one flanking microsatellite region in a subsample of individuals. The two mangrove species shared similar biogeographic histories including: (1) strong genetic breaks between Atlantic and Pacific ocean basins associated with the final closure of the Central American Isthmus (CAI), (2) evidence for simultaneous population declines between the mid-Pleistocene and early Holocene, (3) asymmetric historical migration with higher gene flow from the Atlantic to the Pacific oceans following the direction of the palaeomarine current, and (4) contemporary gene flow between West Africa and South America following the major Atlantic Ocean currents. Despite the remarkable differences in life-history traits of mangrove species, which should have had a strong influence on seed dispersal capability and, thus, population connectivity, we found that vicariant events, climate fluctuations and marine currents have shaped the distribution of genetic diversity in strikingly similar ways