Hierarchical analysis of linkage disequilibrium in Rhizobium populations: evidence for sex?

Population genomics of Vibrionaceae isolated from an endangered oasis reveals local adaptation after an environmental perturbation

Background

In bacteria, pan-genomes are the result of an evolutionary “tug of war” between selection and horizontal gene transfer (HGT). High rates of HGT increase the genetic pool and the effective population size (N_e), resulting in open pan-genomes. In contrast, selective pressures can lead to local adaptation by purging the variation introduced by HGT and mutation, resulting in closed pan-genomes and clonal lineages. In this study, we explored both hypotheses, elucidating the pan-genome of Vibrionaceae isolates after a perturbation event in the endangered oasis of Cuatro Ciénegas Basin (CCB), Mexico, and looking for signals of adaptation to the environments in their genomes.

Results

We obtained 42 genomes of Vibrionaceae distributed in six lineages, two of them did not showed any close reference strain in databases. Five of the lineages showed closed pan-genomes and were associated to either water or sediment environment; their high N_e estimates suggest that these lineages are not from a recent origin. The only clade with an open pan-genome was found in both environments and was formed by ten genetic groups with low N_e, suggesting a recent origin. The recombination and mutation estimators (r/m) ranged from 0.005 to 2.725, which are similar to oceanic Vibrionaceae estimations. However, we identified 367 gene families with signals of positive selection, most of them found in the core genome; suggesting that despite recombination, natural selection moves the Vibrionaceae CCB lineages to local adaptation, purging the genomes and keeping closed pan-genome patterns. Moreover, we identify 598 SNPs associated with an unstructured environment; some of the genes associated with these SNPs were related to sodium transport.

Conclusions

Different lines of evidence suggest that the sampled Vibrionaceae, are part of the rare biosphere usually living under famine conditions. Two of these lineages were reported for the first time. Most Vibrionaceae lineages of CCB are adapted to their micro-habitats rather than to the sampled environments. This pattern of adaptation is concordant with the association of closed pan-genomes and local adaptation.

FINE-SCALE GENETIC AND PHENOTYPIC STRUCTURE IN NATURAL POPULATIONS OF BACILLUS SUBTILIS AND BACILLUS LICHENIFORMIS: IMPLICATIONS FOR BACTERIAL EVOLUTION AND SPECIATION

The genetic and phenotypic structure of sympatric populations of wild bacteria traditionally identified as Bacillus subtilis and B. licheniformis was analyzed. Small soil samples were taken from a single, tiny site in the Sonoran Desert of Arizona, USA, to provide a true population analysis, in contrast to many analyses of genetic structure using bacterial strain collections of widely heterogeneous origin. Genetic analyses of isolates used multilocus enzyme electrophoresis, mismatches in restriction fragment length polymorphism, and variants from Southern hybridization with B. subtilis DNA probes. Phenotypic analyses of isolates used the API test system for detection of growth and acid production on specific carbon sources. The two species were distinct both phenotypically and genetically, despite their known potential for genetic exchange in laboratory experiments. Genic and genotypic diversity were high in both species, and only 16% of observed allozyme variants might possibly be common to both species. Hence, there is probably modest genetic exchange, if any, between the species in nature. Clear hierarchies of population-genetic structure were found for both species. Different types of genetic data yield concordant population structures for B. subtilis. For both species, two-locus and multilocus statistical analyses of linkage demonstrated modest to strong disequilibrium at the species level but truly panmictic subunits within each species. The evidence for extensive genetic recombination within these fine-scale subdivisions is unequivocal, indicating that the sexuality of these bacteria can be well expressed in nature. The relation of these results to processes of bacterial evolution and speciation is discussed.

LONG-TERM EXPERIMENTAL EVOLUTION IN ESCHERICHIA COLI. VII. MECHANISMS MAINTAINING GENETIC VARIABILITY WITHIN POPULATIONS

Six replicate populations of the bacterium Escherichia coli were propagated for more than 10,000 generations in a defined environment. We sought to quantify the variation among clones within these populations with respect to their relative fitness, and to evaluate the roles of three distinct population genetic processes in maintaining this variation. On average, a pair of clones from the same population differed from one another in their relative fitness by approximately 4%. This within-population variation was small compared with the average fitness gain relative to the common ancestor, but it was statistically significant. According to one hypothesis, the variation in fitness is transient and reflects the ongoing substitution of beneficial alleles. We used Fisher's fundamental theorem to compare the observed rate of each population's change in mean fitness with the extent of variation for fitness within that population, but we failed to discern any correspondence between these quantities. A second hypothesis supposes that the variation in fitness is maintained by recurrent deleterious mutations that give rise to a mutation-selection balance. To test this hypothesis, we made use of the fact that two of the six replicate populations had evolved mutator phenotypes, which gave them a genomic mutation rate approximately 100-fold higher than that of the other populations. There was a marginally significant correlation between a population's mutation rate and the extent of its within-population variance for fitness, but this correlation was driven by only one population (whereas two of the populations had elevated mutation rates). Under a third hypothesis, this variation is maintained by frequency-dependent selection, whereby genotypes have an advantage when they are rare relative to when they are common. In all six populations, clones were more fit, on average, when they were rare than when they were common, although the magnitude of the advantage when rare was usually small (~1% in five populations and ~5% in the other). These three hypotheses are not mutually exclusive, but frequency-dependent selection appears to be the primary force maintaining the fitness variation within these experimental populations.

Genetically Determined Response to Artemisinin Treatment in Western Kenyan Plasmodium falciparum Parasites

Genetically determined artemisinin resistance in Plasmodium falciparum has been described in Southeast Asia. The relevance of recently described Kelch 13-propeller mutations for artemisinin resistance in Sub-Saharan Africa parasites is still unknown. Southeast Asia parasites have low genetic diversity compared to Sub-Saharan Africa, where parasites are highly genetically diverse. This study attempted to elucidate whether genetics provides a basis for discovering molecular markers in response to artemisinin drug treatment in P. falciparum in Kenya. The genetic diversity of parasites collected pre- and post- introduction of artemisinin combination therapy (ACT) in western Kenya was determined. A panel of 12 microsatellites and 91 single nucleotide polymorphisms (SNPs) distributed across the P. falciparum genome were genotyped. Parasite clearance rates were obtained for the post-ACT parasites. The 12 microsatellites were highly polymorphic with post-ACT parasites being significantly more diverse compared to pre-ACT (p < 0.0001). The median clearance half-life was 2.55 hours for the post-ACT parasites. Based on SNP analysis, 15 of 90 post-ACT parasites were single-clone infections. Analysis revealed 3 SNPs that might have some causal association with parasite clearance rates. Further, genetic analysis using Bayesian tree revealed parasites with similar clearance phenotypes were more closely genetically related. With further studies, SNPs described here and genetically determined response to artemisinin treatment might be useful in tracking artemisinin resistance in Kenya.

Population genomics of the symbiotic plasmids of sympatric nitrogen-fixing Rhizobium species associated with Phaseolus vulgaris

Cultivated common beans are the primary protein source for millions of people around the world who subsist on low-input agriculture, enabled by the symbiotic N2 -fixation these legumes perform in association with rhizobia. Within a single agricultural plot, multiple Rhizobium species can nodulate bean roots, but it is unclear how genetically isolated these species remain in sympatry. To better understand this issue, we sequenced and compared the genomes of 33 strains isolated from the rhizosphere and root nodules of a particular bean variety grown in the same agricultural plot. We found that the Rhizobium species we observed coexist with low genetic recombination across their core genomes. Accessory plasmids thought to be necessary for the saprophytic lifestyle in soil show similar levels of genetic isolation, but with higher rates of recombination than the chromosomes. However, the symbiotic plasmids are extremely similar, with high rates of recombination and do not appear to have co-evolved with the chromosome or accessory plasmids. Therefore, while Rhizobium species are genetically isolated units within the microbial community, a common symbiotic plasmid allows all Rhizobium species to engage in symbiosis with the same host in a single agricultural plot.

Genetic variability and population structure of Plasmodium falciparum parasite populations from different malaria ecological regions of Kenya

Transmission intensity, movement of human and vector hosts, biogeographical features, and malaria control measures are some of the important factors that determine Plasmodium falciparum parasite genetic variability and population structure. Kenya has different malaria ecologies which might require different disease intervention methods. Refined parasite population genetic studies are critical for informing malaria control and elimination strategies. This study describes the genetic diversity and population structure of P. falciparum parasites from the different malaria ecological zones in Kenya. Twelve multi-locus microsatellite (MS) loci previously described were genotyped in 225 P. falciparum isolates collected between 2012 and 2013 from five sites; three in lowland endemic regions (Kisumu, Kombewa, and Malindi) and two in highland, epidemic regions (Kisii and Kericho). Parasites from the lowland endemic and highland epidemic regions of western Kenya had high genetic diversity compared to coastal lowland endemic region of Kenya [Malindi]. The Kenyan parasites had a mean genetic differentiation index (FST) of 0.072 (p=0.011). The multi-locus genetic analysis of the 12 MS revealed all the parasites had unique haplotypes. Significant linkage disequilibrium (LD) was observed in all the five parasite populations. Kisumu had the most significant index of association values (0.16; p<0.0001) whereas Kisii had the least significant index of association values (0.03; p<0.0001). Our data suggest high genetic diversity in Kenyan parasite population with the exception of parasite from Malindi where malaria has been on the decline. The presence of significant LD suggests that there is occurrence of inbreeding in the parasite population. Parasite populations from Kisii showed the strongest evidence for epidemic population structure whereas the rest of the regions showed panmixia. Defining the genetic diversity of the parasites in different ecological regions of Kenya after introduction of the artemether-lumefantrine is important in refining the spread of drug resistant strains and malaria transmission for more effective control and eventual elimination of malaria in Kenya.

Genomic characterization of Sinorhizobium meliloti AK21, a wild isolate from the Aral Sea Region

The symbiotic, nitrogen-fixing bacterium Sinorhizobium meliloti has been widely studied due to its ability to improve crop yields through direct interactions with leguminous plants. S. meliloti AK21 is a wild type strain that forms nodules on Medicago plants in saline and drought conditions in the Aral Sea Region. The aim of this work was to establish the genetic similarities and differences between S. meliloti AK21 and the reference strain S. meliloti 1021. Comparative genome hybridization with the model reference strain S. meliloti 1021 yielded 365 variable genes, grouped into 11 regions in the three main replicons in S. meliloti AK21. The most extensive regions of variability were found in the symbiotic plasmid pSymA, which also contained the largest number of orthologous and polymorphic sequences identified by suppression subtractive hybridization. This procedure identified a large number of divergent sequences and others without homology in the databases, the further investigation of which could provide new insight into the alternative metabolic pathways present in S. meliloti AK21. We identified a plasmid replication module from the repABC replicon family, together with plasmid mobilization-related genes (traG and a VirB9-like protein), which suggest that this indigenous isolate harbors an accessory plasmid. Furthermore, the transcriptomic profiles reflected differences in gene content and regulation between S. meliloti AK21 and S. meliloti 1021 (ExpR and PhoB regulons), but provided evidence for an as yet unknown, alternative mechanism involving activation of the cbb3 terminal oxidase. Finally, phenotypic microarrays characterization revealed a greater versatility of substrate use and chemical degradation than for S. meliloti 1021.

Microsatellite markers reveal low levels of population sub-structuring of Plasmodium falciparum in southwestern Nigeria

BACKGROUND

Genetic diversity studies provide evidence of Plasmodium falciparum differentiation that could affect fitness and adaptation to drugs and target antigens for vaccine development. This study describes the genetic structure of P. falciparum populations in urban and rural sites from southwestern Nigeria.

METHODOLOGY

Ten neutral microsatellite loci were genotyped in 196 P. falciparum infections from three localities: Aramoko-Ekiti, a rural community; Lekki, an urban location and Badagry, a peri-urban border settlement. Analysis was performed on the genetic diversity, linkage disequilibrium, population structure and inter-population differentiation.

RESULTS

Allelic diversity values were similar across all populations, with mean expected heterozygosity (HE) values between 0.65 and 0.79. No matching multilocus haplotypes were found and analysis of multilocus LD showed no significant index of association. Genetic differentiation between populations was low (ΦPT = 0.017).

CONCLUSION

The absence of detectable population structure of P. falciparum in southwestern Nigeria is evident in the lack of significant differentiation between populations separated by about 200 km. This implies that a fairly uniform malaria control strategy may be effective over a wide geographic range in this highly endemic region. However, more wide-scale survey across the country will be required to inform malaria control in this large and densely populated endemic region.

Travel, sex, and food: what's speciation got to do with it?

We discuss the potential interactions among travel (dispersal and gene flow), bacterial "sex" (mainly as horizontal gene transfer), and food (metabolic plasticity and responses to nutrient availability) in shaping microbial communities. With regard to our work at a unique desert oasis, the Cuatro Ciénegas Basin in Coahuila, Mexico, we propose that diversification and low phosphorus availability, in combination with mechanisms for nutrient recycling and community cohesion, result in enhanced speciation through reproductive as well as geographic isolation. We also discuss these mechanisms in the broader sense of ecology and evolution. Of special relevance to astrobiology and central to evolutionary biology, we ask why there are so many species on Earth and provide a working hypothesis and a conceptual framework within which to consider the question. Key Words: Microbial ecology-Microbial mats-Evolution-Horizontal gene transfer-Metabolism.

The impact of interspecific competition on lineage evolution and a rapid peak shift by interdemic genetic mixing in experimental bacterial populations

Epistatic interactions between genes in the genome constrain the accessible evolutionary paths of lineages. Two factors involving epistasis that can affect the evolutionary path and fate of lineages were investigated. The first factor concerns the impact of competition with another species lineage that has different epistatic constraints. Five enteric bacterial populations were evolved by point mutation in medium containing a single limiting resource. Single-species and two-species cultures were used to determine whether different asexual lineages have different capacities for producing variants due to epistatic constraints, and whether their survival is determined by local inter-lineage competition with different species. Local inter-lineage competition quickly resulted in one successful lineage, with another lineage becoming extinct before finding a higher peak. The second factor concerns a peak-shifting process, and whether the sexual recombination between different demes can cause peak shifts was investigated. An Escherichia coli population consisting of a male (Hfr) and female strain (F(-)) was evolved in a single limiting resource and compared to evolving populations containing the male or female strain alone. The E. coli sexual lineage was successful due to its ability to escape lower peaks and reach a higher peak, not because of a rapid approach to the nearest local peak the male or female asexual lineage could reach. The data in this study demonstrate that lineage survivability can be determined by the ability to produce beneficial mutations and checked by local competition between lineages of different species. Interspecific competition may prevent a population from evolving through crossing fitness valleys or adaptive ridges if it requires many generations to achieve peak shifts. The data also show that genomic recombination between different conspecific lineages can rapidly carry the combined lineage to a higher peak.

Genomic lineages of Rhizobium etli revealed by the extent of nucleotide polymorphisms and low recombination

Background

Most of the DNA variations found in bacterial species are in the form of single nucleotide polymorphisms (SNPs), but there is some debate regarding how much of this variation comes from mutation versus recombination. The nitrogen-fixing symbiotic bacteria Rhizobium etli is highly variable in both genomic structure and gene content. However, no previous report has provided a detailed genomic analysis of this variation at nucleotide level or the role of recombination in generating diversity in this bacterium. Here, we compared draft genomic sequences versus complete genomic sequences to obtain reliable measures of genetic diversity and then estimated the role of recombination in the generation of genomic diversity among Rhizobium etli.

Results

We identified high levels of DNA polymorphism in R. etli, and found that there was an average divergence of 4% to 6% among the tested strain pairs. DNA recombination events were estimated to affect 3% to 10% of the genomic sample analyzed. In most instances, the nucleotide diversity (π) was greater in DNA segments with recombinant events than in non-recombinant segments. However, this degree of recombination was not sufficiently large to disrupt the congruence of the phylogenetic trees, and further evaluation of recombination in strains quartets indicated that the recombination levels in this species are proportionally low.

Conclusion

Our data suggest that R. etli is a species composed of separated lineages with low homologous recombination among the strains. Horizontal gene transfer, particularly via the symbiotic plasmid characteristic of this species, seems to play an important role in diversity but the lineages maintain their evolutionary cohesiveness.

Local expansion of a panmictic lineage of water bloom-forming cyanobacterium Microcystis aeruginosa

In previous studies, we have demonstrated that the population structure of the bloom-forming cyanobacterium Microcystis aeruginosa is clonal. Expanded multilocus sequence typing analysis of M. aeruginosa using 412 isolates identified five intraspecific lineages suggested to be panmictic while maintaining overall clonal structure probably due to a reduced recombination rate between lineages. Interestingly, since 2005 most strains belonging to one of these panmictic clusters (group G) have been found in a particular locality (Lake Kasumigaura Basin) in Japan. In this locality, multiple, similar but distinct genotypes of this lineage predominated in the bloom, a pattern that is unprecedented for M. aeruginosa. The population structure underlying blooms associated with this lineage is comparable to epidemics of pathogens. Our results may reveal an expansion of the possible adaptive lineage in a localized aquatic environment, providing us with a unique opportunity to investigate its ecological and biogeographical consequences.

The population genomics of begomoviruses: global scale population structure and gene flow

BACKGROUND

The rapidly growing availability of diverse full genome sequences from across the world is increasing the feasibility of studying the large-scale population processes that underly observable pattern of virus diversity. In particular, characterizing the genetic structure of virus populations could potentially reveal much about how factors such as geographical distributions, host ranges and gene flow between populations combine to produce the discontinuous patterns of genetic diversity that we perceive as distinct virus species. Among the richest and most diverse full genome datasets that are available is that for the dicotyledonous plant infecting genus, Begomovirus, in the Family Geminiviridae. The begomoviruses all share the same whitefly vector, are highly recombinogenic and are distributed throughout tropical and subtropical regions where they seriously threaten the food security of the world's poorest people.

RESULTS

We focus here on using a model-based population genetic approach to identify the genetically distinct sub-populations within the global begomovirus meta-population. We demonstrate the existence of at least seven major sub-populations that can further be sub-divided into as many as thirty four significantly differentiated and genetically cohesive minor sub-populations. Using the population structure framework revealed in the present study, we further explored the extent of gene flow and recombination between genetic populations.

CONCLUSIONS

Although geographical barriers are apparently the most significant underlying cause of the seven major population sub-divisions, within the framework of these sub-divisions, we explore patterns of gene flow to reveal that both host range differences and genetic barriers to recombination have probably been major contributors to the minor population sub-divisions that we have identified. We believe that the global Begomovirus population structure revealed here could facilitate population genetics studies into how central parameters of population genetics namely selection, recombination, mutation, gene flow, and genetic drift shape the global begomovirus diversity.

Genetic relationships between clinical and non-clinical strains of Yersinia enterocolitica biovar 1A as revealed by multilocus enzyme electrophoresis and multilocus restriction typing

Background

Genetic relationships among 81 strains of Y. enterocolitica biovar 1A isolated from clinical and non-clinical sources were discerned by multilocus enzyme electrophoresis (MLEE) and multilocus restriction typing (MLRT) using six loci each. Such studies may reveal associations between the genotypes of the strains and their sources of isolation.

Results

All loci were polymorphic and generated 62 electrophoretic types (ETs) and 12 restriction types (RTs). The mean genetic diversity (H) of the strains by MLEE and MLRT was 0.566 and 0.441 respectively. MLEE (DI = 0.98) was more discriminatory and clustered Y. enterocolitica biovar 1A strains into four groups, while MLRT (DI = 0.77) identified two distinct groups. BURST (Based Upon Related Sequence Types) analysis of the MLRT data suggested aquatic serotype O:6,30-6,31 isolates to be the ancestral strains from which, clinical O:6,30-6,31 strains might have originated by host adaptation and genetic change.

Conclusion

MLEE revealed greater genetic diversity among strains of Y. enterocolitica biovar 1A and clustered strains in four groups, while MLRT grouped the strains into two groups. BURST analysis of MLRT data nevertheless provided newer insights into the probable evolution of clinical strains from aquatic strains.

Genetic and metabolic divergence within a Rhizobium leguminosarum bv. trifolii population recovered from clover nodules

Rhizobia are able to establish symbiosis with leguminous plants and usually occupy highly complex soil habitats. The large size and complexity of their genomes are considered advantageous, possibly enhancing their metabolic and adaptive potential and, in consequence, their competitiveness. A population of Rhizobium leguminosarum bv. trifolii organisms recovered from nodules of several clover plants growing in each other's vicinity in the soil was examined regarding possible relationships between their metabolic-physiological properties and their prevalence in such a local population. Genetic and metabolic variability within the R. leguminosarum bv. trifolii strains occupying nodules of several plants was of special interest, and both types were found to be considerable. Moreover, a prevalence of metabolically versatile strains, i.e., those not specializing in utilization of any group of substrates, was observed by combining statistical analyses of Biolog test results with the frequency of occurrence of genetically distinct strains. Metabolic versatility with regard to nutritional requirements was not directly advantageous for effectiveness in the symbiotic interaction with clover: rhizobia with specialized metabolism were more effective in symbiosis but rarely occurred in the population. The significance of genetic and, especially, metabolic complexity of bacteria constituting a nodule population is discussed in the context of strategies employed by bacteria in competition.

Phenotypic and genetic diversity in Sinorhizobium meliloti and S. medicae from drought and salt affected regions of Morocco

BACKGROUND

Sinorhizobium meliloti and S. medicae are symbiotic nitrogen fixing bacteria in root nodules of forage legume alfalfa (Medicago sativa L.). In Morocco, alfalfa is usually grown in marginal soils of arid and semi-arid regions frequently affected by drought, extremes of temperature and soil pH, soil salinity and heavy metals, which affect biological nitrogen fixing ability of rhizobia and productivity of the host. This study examines phenotypic diversity for tolerance to the above stresses and genotypic diversity at Repetitive Extragenic Pallindromic DNA regions of Sinorhizobium nodulating alfalfa, sampled from marginal soils of arid and semi-arid regions of Morocco.

RESULTS

RsaI digestion of PCR amplified 16S rDNA of the 157 sampled isolates, assigned 136 isolates as S. meliloti and the rest as S. medicae. Further phenotyping of these alfalfa rhizobia for tolerance to the environmental stresses revealed a large degree of variation: 55.41%, 82.16%, 57.96% and 3.18% of the total isolates were tolerant to NaCl (>513 mM), water stress (-1.5 MPa), high temperature (40 degrees C) and low pH (3.5), respectively. Sixty-seven isolates of S. meliloti and thirteen isolates of S. medicae that were tolerant to salinity were also tolerant to water stress. Most of the isolates of the two species showed tolerance to heavy metals (Cd, Mn and Zn) and antibiotics (chloramphenicol, spectinomycin, streptomycin and tetracycline). The phenotypic clusters observed by the cluster analysis clearly showed adaptations of the S. meliloti and S. medicae strains to the multiple stresses. Genotyping with rep-PCR revealed higher genetic diversity within these phenotypic clusters and classified all the 157 isolates into 148 genotypes. No relationship between genotypic profiles and the phenotypes was observed. The Analysis of Molecular Variance revealed that largest proportion of significant (P < 0.01) genetic variation was distributed within regions (89%) than among regions (11%).

CONCLUSION

High degree of phenotypic and genotypic diversity is present in S. meliloti and S. medicae populations from marginal soils affected by salt and drought, in arid and semi-arid regions of Morocco. Some of the tolerant strains have a potential for exploitation in salt and drought affected areas for biological nitrogen fixation in alfalfa.

A comparison of homologous recombination rates in bacteria and archaea

It is a standard practice to test for the signature of homologous recombination in studies examining the genetic diversity of bacterial populations. Although it has emerged that homologous recombination rates can vary widely between species, comparing the results from different studies is made difficult by the diversity of estimation methods used. Here, Multi Locus Sequence Typing (MLST) datasets from a wide variety of bacteria and archaea are analyzed using the ClonalFrame method. This enables a direct comparison between species and allows for a first exploration of the question whether phylogeny or ecology is the primary determinant of homologous recombination rate.

Microbial endemism: does phosphorus limitation enhance speciation?

There is increasing evidence for the existence of unique ecosystems that are dominated by locally adapted microbiota which harbour distinct lineages and biological capabilities, much like the macrobiota of Darwin's Galapagos Islands. As a primary example of such a system, we highlight key discoveries from the Cuatro Ciénegas basin in Mexico. We argue that high microbial endemism requires a combination of geographical isolation, long-term continuity and mechanisms for reducing the intensity of horizontal gene transfer (HGT). We also propose that strong phosphorus limitation has an important role in microbial diversification by reducing the intensity of HGT.

Contrasting genetic structure in Plasmodium vivax populations from Asia and South America

Populations of Plasmodium falciparum show striking differences in linkage disequilibrium, population differentiation and diversity, but only fragmentary data exists on the genetic structure of Plasmodium vivax. We genotyped nine tandem repeat loci bearing 2-8 bp motifs from 345 P. vivax infections collected from three Asian countries and from five locations in Colombia. We observed 9-37 alleles per locus and high diversity (He=0.72-0.79, mean=0.75) in all countries. Numbers of multiple clone infections varied considerably: these were rare in Colombia and India, but > 60% of isolates carried multiple alleles in at least one locus in Thailand and Laos. However, only one or two of the nine loci show >1 allele in many samples, suggesting that mutation within infections may result in overestimation of true multiple carriage rates. Identical nine-locus genotypes were frequently found in Colombian populations, contributing to strong linkage disequilibrium. These identical genotypes were strongly clustered in time, consistent with epidemic transmission of clones and subsequent breakdown of allelic associations, suggesting high rates of inbreeding and low effective recombination rates in this country. In contrast, identical genotypes were rare and loci were randomly associated in all three Asian populations, consistent with higher rates of outcrossing and recombination. We observed low but significant differentiation between different Asian countries (standardized FST = 0.13-0.45). In comparison, we see greater differentiation between collection locations within Colombia (standardized FST = 0.4-0.7), and strong differentiation between continents (standardized FST = 0.48-0.79). The observed heterogeneity in multiple clone carriage rates, linkage disequilibrium and population differentiation are similar in some, but not all, respects to those observed in P. falciparum, and have important implications for the design of association mapping studies, and interpretation of P. vivax epidemiology.

Biogeography of actinomycete communities and type II polyketide synthase genes in soils collected in New Jersey and Central Asia

Soil microbial communities are believed to be comprised of thousands of different bacterial species. One prevailing idea is that "everything is everywhere, and the environment selects," implying that all types of bacteria are present in all environments where their growth requirements are met. We tested this hypothesis using actinomycete communities and type II polyketide synthase (PKS) genes found in soils collected from New Jersey and Uzbekistan (n = 91). Terminal restriction fragment length polymorphism analysis using actinomycete 16S rRNA and type II PKS genes was employed to determine community profiles. The terminal fragment frequencies in soil samples had a lognormal distribution, indicating that the majority of actinomycete phylotypes and PKS pathways are present infrequently in the environment. Less than 1% of peaks were detected in more than 50% of samples, and as many as 18% of the fragments were unique and detected in only one sample. Actinomycete 16S rRNA fingerprints clustered by country of origin, indicating that unique populations are present in North America and Central Asia. Sequence analysis of type II PKS gene fragments cloned from Uzbek soil revealed 35 novel sequence clades whose levels of identity to genes in the GenBank database ranged from 68 to 92%. The data indicate that actinomycetes are patchily distributed but that distinct populations are present in North American and Central Asia. These results have implications for microbial bioprospecting and indicate that the cosmopolitan actinomycete species and PKS pathways may account for only a small proportion of the total diversity in soil.

Biogeography: an emerging cornerstone for understanding prokaryotic diversity, ecology, and evolution

New questions about microbial ecology and diversity combined with significant improvement in the resolving power of molecular tools have helped the reemergence of the field of prokaryotic biogeography. Here, we show that biogeography may constitute a cornerstone approach to study diversity patterns at different taxonomic levels in the prokaryotic world. Fundamental processes leading to the formation of biogeographic patterns are examined in an evolutionary and ecological context. Based on different evolutionary scenarios, biogeographic patterns are thus posited to consist of dramatic range expansion or regression events that would be the results of evolutionary and ecological forces at play at the genotype level. The deterministic or random nature of those underlying processes is, however, questioned in light of recent surveys. Such scenarios led us to predict the existence of particular genes whose presence or polymorphism would be associated with cosmopolitan taxa. Furthermore, several conceptual and methodological pitfalls that could hamper future developments of the field are identified, and future approaches and new lines of investigation are suggested.

Recombination and selection shape the molecular diversity pattern of nitrogen-fixing Sinorhizobium sp. associated to Medicago

We investigate the genetic structure and molecular selection pattern of a sympatric population of Sinorhizobium meliloti and Sinorhizobium medicae. These bacteria fix nitrogen in association with plants of the genus Medicago. A set of 116 isolates were obtained from a soil sample, from root nodules of three groups of plants representing among-species, within-species and intraline diversity in the Medicago genus. Bacteria were characterized by sequencing at seven loci evenly distributed along the genome of both Sinorhizobium species, covering the chromosome and the two megaplasmids. We first test whether the diversity of host plants influence the bacterial diversity recovered. Using the same data set, we then analyse the selective pattern at each locus. There was no relationship between the diversity of Medicago plants that were used for sampling and the diversity of their symbionts. However, we found evidence of selection within each of the two main symbiotic regions, located on the two different megaplasmids. Purifying selection or a selective sweep was found to occur in the nod genomic region, which includes genes involved in nodulation specificity, whereas balancing selection was detected in the exo region, close to genes involved in exopolysaccharide production. Such pattern likely reflects the interaction between host plants and bacterial symbionts, with a possible conflict of interest between plants and cheater bacterial genotypes. Recombination appears to occur preferentially within and among loci located on megaplasmids, rather than within the chromosome. Thus, recombination may play an important role in resolving this conflict by allowing different selection patterns at different loci.

Large-scale genetic variation of the symbiosis-required megaplasmid pSymA revealed by comparative genomic analysis of Sinorhizobium meliloti natural strains

Background

Sinorhizobium meliloti is a soil bacterium that forms nitrogen-fixing nodules on the roots of leguminous plants such as alfalfa (Medicago sativa). This species occupies different ecological niches, being present as a free-living soil bacterium and as a symbiont of plant root nodules. The genome of the type strain Rm 1021 contains one chromosome and two megaplasmids for a total genome size of 6 Mb. We applied comparative genomic hybridisation (CGH) on an oligonucleotide microarrays to estimate genetic variation at the genomic level in four natural strains, two isolated from Italian agricultural soil and two from desert soil in the Aral Sea region.

Results

From 4.6 to 5.7 percent of the genes showed a pattern of hybridisation concordant with deletion, nucleotide divergence or ORF duplication when compared to the type strain Rm 1021. A large number of these polymorphisms were confirmed by sequencing and Southern blot. A statistically significant fraction of these variable genes was found on the pSymA megaplasmid and grouped in clusters. These variable genes were found to be mainly transposases or genes with unknown function.

Conclusion

The obtained results allow to conclude that the symbiosis-required megaplasmid pSymA can be considered the major hot-spot for intra-specific differentiation in S. meliloti.

Recombination shapes the natural population structure of the hyperthermophilic archaeon Sulfolobus islandicus

Although microorganisms make up the preponderance of the biodiversity on Earth, the ecological and evolutionary factors that structure microbial populations are not well understood. We investigated the genetic structure of a thermoacidophilic crenarchaeal species, Sulfolobus islandicus, using multilocus sequence analysis of six variable protein-coding loci on a set of 60 isolates from the Mutnovsky region of Kamchatka, Russia. We demonstrate significant incongruence among gene genealogies and a lack of association between alleles consistent with recombination rates greater than the rate of mutation. The observation of high relative rates of recombination suggests that the structure of this natural population does not fit the periodic selection model often used to describe populations of asexual microorganisms. We propose instead that frequent recombination among closely related individuals prevents periodic selection from purging diversity and provides a fundamental cohesive mechanism within this and perhaps other archaeal species.

Marine actinobacteria: perspectives, challenges, future directions

In this paper we evaluate the current state of research on the biology and biotechnology of marine actinobacteria. The topics covered include the abundance, diversity, novelty and biogeographic distribution of marine actinobacteria, ecosystem function, bioprospecting, and a new approach to the exploration of actinobacterial taxonomic space. An agenda for future marine actinobacterial research is suggested based upon consideration of the above issues.

The ecology and genetics of microbial diversity

Natural communities of microbes are often diverse, a fact that is difficult to reconcile with the action of natural selection in simple, uniform environments. We suggest that this apparent paradox may be resolved by considering the origin and fate of diversity in an explicitly ecological context. Here, we review insights into the ecological and genetic causes of diversity that stem from experiments with microbial populations evolving in the defined conditions of the laboratory environment. These studies highlight the importance of environmental structure in governing the fate of diversity and shed light on the genetic mechanisms generating diversity. We conclude by emphasizing the importance of placing detailed molecular-level studies within the context of a sound ecological and evolutionary framework.

Role of HLA class II alleles in susceptibility to and protection from localized cutaneous leishmaniasis

Localized cutaneous leishmaniasis (LCL) is the prevalent form of leishmaniasis in Mexico. It is limited to the skin; reversible upon treatment and the host cellular immune response is intact. Several genes that influence the expression of LCL have been described in the mouse. In humans, we, as well as others, have demonstrated that HLA-DQ3 antigens seem to play some role in host susceptibility. We therefore analyzed at the DNA level, the class II loci of the same patients that were previously studied by serology. The purpose of this study was to assess the contribution of HLA DR, DQ, and DP genes in the protection and/or the susceptibility to LCL. Sixty-five patients with LCL from Comalcalco, state of Tabasco, were recruited and 100 healthy controls were included for comparison. All were Mexican Mestizos. DRB1, DQA1, DQB1, DPA1, and DPB1 alleles were typed using two different methods: PCR-SSO and PCR-SSP. Results indicate that class II genes are relevant for the expression of LCL and several loci contribute independently and sinergically. DRB1*0407 participates in susceptibility with an etiological fraction (EF) of 20% and an odds ratio (OR) of 2.92. Two additional susceptibility genes were found. These are located to the DP locus: DPA1*0401 (OR = 10.07; EF=7%) and DPB1*0101 (OR = 5.99 EF = 13%). Resistance was found associated to DPB1*0401, thus *0401 "motif" could be an ideal candidate for the development of a vaccine. DR2 (DRB1*1500+DRB1*1600) has also a significant p for protection, suggesting that the sequence common to this group of antigens may anchor parasite peptides which trigger a protective response.

Genetic diversity, recombination and cryptic clades in Pseudomonas viridiflava infecting natural populations of Arabidopsis thaliana

Species-level genetic diversity and recombination in bacterial pathogens of wild plant populations have been nearly unexplored. Pseudomonas viridiflava is a common natural bacterial pathogen of Arabidopsis thaliana, for which pathogen defense genes and mechanisms are becoming increasing well known. The genetic variation contained within a worldwide sample of P. viridiflava collected from wild populations of A. thaliana was investigated using five genomic sequence fragments totaling 2.3 kb. Two distinct and deeply diverged clades were found within the P. viridiflava sample and in close proximity in multiple populations, each genetically diverse with synonymous variation as high as 9.3% in one of these clades. Within clades, there is evidence of frequent recombination within and between each sequenced locus and little geographic differentiation. Isolates from both clades were also found in a small sample of other herbaceous species in Midwest populations, indicating a possibly broad host range for P. viridiflava. The high levels of genetic variation and recombination together with a lack of geographic differentiation in this pathogen distinguish it from other bacterial plant pathogens for which intraspecific variation has been examined.

Biogeography, evolution, and diversity of epibionts in phototrophic consortia

Motile phototrophic consortia are highly regular associations in which numerous cells of green sulfur bacteria surround a flagellated colorless beta-proteobacterium in the center. To date, seven different morphological types of such consortia have been described. In addition, two immotile associations involving green sulfur bacteria are known. By employing a culture-independent approach, different types of phototrophic consortia were mechanically isolated by micromanipulation from 14 freshwater environments, and partial 16S rRNA gene sequences of the green sulfur bacterial epibionts were determined. In the majority of the lakes investigated, different types of phototrophic consortia were found to co-occur. In all cases, phototrophic consortia with the same morphology from the same habitat contained only a single epibiont phylotype. However, morphologically indistinguishable phototrophic consortia collected from different lakes contained different epibionts. Overall, 19 different types of epibionts were detected in the present study. Whereas the epibionts within one geographic region were very similar (Dice coefficient, 0.582), only two types of epibionts were found to occur on both the European and North American continents (Dice coefficient, 0.190). None of the epibiont 16S rRNA gene sequences have been detected so far in free-living green sulfur bacteria, suggesting that the interaction between epibionts and chemotrophic bacteria in the phototrophic consortia is an obligate interaction. Based on our phylogenetic analysis, the epibiont sequences are not monophyletic. Thus, the ability to form symbiotic associations either arose independently from different ancestors or was present in a common ancestor prior to the radiation of green sulfur bacteria and the transition to the free-living state in independent lineages. The present study thus demonstrates that there is great diversity and nonrandom geographical distribution of phototrophic consortia in the natural environment.

An ecological perspective on bacterial biodiversity

Bacteria may be one of the most abundant and species-rich groups of organisms, and they mediate many critical ecosystem processes. Despite the ecological importance of bacteria, past practical and theoretical constraints have limited our ability to document patterns of bacterial diversity and to understand the processes that determine these patterns. However, recent advances in molecular techniques that allow more thorough detection of bacteria in nature have made it possible to examine such patterns and processes. Here, we review recent studies of the distribution of free-living bacterial diversity and compare our current understanding with what is known about patterns in plant and animal diversity. From these recent studies a preliminary picture is emerging: bacterial diversity may exhibit regular patterns, and in some cases these patterns may be qualitatively similar to those observed for plants and animals.

Multilocus sequence typing for comparison of veterinary and human isolates of Campylobacter jejuni

Multilocus sequence typing (MLST) has been applied to 266 Campylobacter jejuni isolates, mainly from veterinary sources, including cattle, sheep, poultry, pigs, pets, and the environment, as well as isolates from human cases of campylobacteriosis. The populations of veterinary and human isolates overlap, suggesting that most veterinary sources should be considered reservoirs of pathogenic campylobacters. There were some associations between source and sequence type complex, indicating that host or source adaptation may exist. The pig isolates formed a distinct group by MLST and may well represent a potential pig-adapted clone of C. jejuni. A subset (n = 82) of isolates was reanalyzed with a second MLST scheme which provided a unique set of isolates that had been analyzed at a total of 12 loci. The distribution of isolates among the complexes in each of the two schemes was similar but not identical. In addition to isolates from human outbreaks, one group of isolates that were not epidemiologically linked was also identical at all 12 loci. This group of isolates is believed to represent another stable strain of C. jejuni.

Genetic and ecological structure of Hafnia alvei in Australia

Multilocus enzyme electrophoresis of 161 Hafnia alvei isolates from 158 hosts and 3 water column samples collected in Australia revealed that this species consists of two genetically distinct groups. The two groups of H. alvei differed significantly in their genetic structure and host distribution. The taxonomic class of the host but not geographic locality explained a significant proportion of the observed genetic and biochemical variation among strains within each genetic group.

Detection of Genetic Variation in Taenia solium

Genetic variability among Taenia solium isolates was studied in 160 cysticerci from 6 pigs, 4 from Mexico, 1 from Honduras, and 1 from Tanzania. Random amplified polymorphic DNA (RAPD) analysis performed with 4 commercial primers showed 88% polymorphic loci and an average heterozygosity of 0.077; however, several alleles were fixed within each isolate. Linkage disequilibrium analysis indicated that 3 of the 6 isolates had a random association of alleles, whereas the other 3 had a clonal structure. These results suggest the existence of local lineages in T. solium, with events of genetic recombination within them.

Clonal population structure of Pseudomonas avellanae strains of different origin based on multilocus enzyme electrophoresis

To assess the genetic diversity and genetic relationships ofPseudomonas avellanae, the causative agent of hazelnut decline, a total of 102 strains, obtained from central Italy (provinces of Viterbo and Rome) and northern Greece, were studied using multilocus enzyme electrophoresis (MLEE). Their allelic variation in 10 loci was determined. All loci were polymorphic and 53 electrophoretic types (ETs) were identified from the total sample. The mean genetic diversity (H) was 0·65 and this value ranged from 0·37 for the least polymorphic to 0·82 for the most polymorphic locus. The dendrogram originated from MLEE data indicated two main groups of ETs, A and B. The groups do not appear to be correlated to the geographic origin of the strains, although all the ETs from northern Greece clustered into subgroup B1.Pseudomonas syringaepv.actinidiaeandP. syringaepv.theae, included in the analysis as outgroups, clustered apart. The index of association (IA) forP. avellanaewas 0·90. TheIAvalues were always significantly different from zero for the population subsets studied and no epidemic structure was found. These results would indicate that the population structure ofP. avellanaeis clonal either in northern Greece or in central Italy. The recent outbreaks of the bacterium in new areas of hazelnut cultivation would explain the current clonal structure that is persisting over decades.

Rhizobium etli and Rhizobium gallicum nodulate common bean (Phaseolus vulgaris) in a traditionally managed milpa plot in Mexico: population genetics and biogeographic implications

The stability of the genetic structure of rhizobial populations nodulating Phaseolus vulgaris cultivated in a traditionally managed milpa plot in Mexico was studied over three consecutive years. The set of molecular markers analyzed (including partial rrs, glnII, nifH, and nodB sequences), along with host range experiments, placed the isolates examined in Rhizobium etli bv. phaseoli and Rhizobium gallicum bv. gallicum. Cluster analysis of multilocus enzyme electrophoresis and plasmid profile data separated the two species and identified numerically dominant clones within each of them. Population genetic analyses showed that there was high genetic differentiation between the two species and that there was low intrapopulation differentiation of the species over the 3 years. The results of linkage disequilibrium analyses are consistent with an epidemic genetic structure for both species, with frequent genetic exchange taking place within conspecific populations but not between the R. etli and R. gallicum populations. A subsample of isolates was selected and used for 16S ribosomal DNA PCR-restriction fragment length polymorphism analysis, nifH copy number determination, and host range experiments. Plasmid profiles and nifH hybridization patterns also revealed the occurrence of lateral plasmid transfer among distinct multilocus genotypes within species but not between species. Both species were recovered from nodules of the same plants, indicating that mechanisms other than host, spatial, or temporal isolation may account for the genetic barrier between the species. The biogeographic implications of finding an R. gallicum bv. gallicum population nodulating common bean in America are discussed.

Linkage disequilibrium between minisatellite loci supports clonal evolution of Mycobacterium tuberculosis in a high tuberculosis incidence area

Deciphering the structure of pathogen populations is instrumental for the understanding of the epidemiology and history of infectious diseases and for their control. Although Mycobacterium tuberculosis is the most widespread infectious agent in humans, its actual population structure has remained hypothetical until now because: (i) its structural genes are poorly polymorphic; (ii) adequate samples and appropriate statistics for population genetic analysis have not been considered. To investigate this structure, we analysed the statistical associations (linkage disequilibrium) between 12 independent M. tuberculosis minisatellite-like loci by high-throughput genotyping within a model population of 209 isolates representative of the genetic diversity in an area with a very high incidence of tuberculosis. These loci contain variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs). Highly significant linkage disequilibrium was detected among the MIRU-VNTR loci in this model. This linkage disequilibrium was also evident when the MIRU-VNTR types were compared with the IS6110 restriction fragment length polymorphism types. These results support a predominant clonal evolution of M. tuberculosis.

Identification of complex composition, strong strain diversity and directional selection in local Pseudomonas stutzeri populations from marine sediment and soils

Members of Pseudomonas stutzeri have been isolated world-wide from various habitats including aquatic and terrestrial ecosystems. The global population has a clonal structure, is of exceptionally high genetic diversity and has been grouped into eight genomovars. We have analysed four local populations (n = 89-125) from three geographically separated habitats (two from a marine sediment and two from different soils) by random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR), restriction fragment length polymorphism (RFLP) of the rpoB gene and 16S rDNA sequences in order to quantify the influence of evolutionary forces on closely related groups of proliferating cells in situ. All populations consisted of a complex structure of genomic subgroups with variable numbers of members. The analyses revealed that the two populations from marine sediment were rather similar. At least three of the populations were influenced by migrational input as concluded from the presence of members from different genomovars. All populations showed very high strain diversity suggesting strong influence of mutability. Neutrality tests indicated that two or possibly three of the populations were shaped by directional selection. Thus, the local populations of P. stutzeri reflect already the high genetic diversity of the global population and are influenced, to different extents, by migration, mutation and directional selection.

Statistics of selectively neutral genetic variation

Random models of evolution are instrumental in extracting rates of microscopic evolutionary mechanisms from empirical observations on genetic variation in genome sequences. In this context it is necessary to know the statistical properties of empirical observables (such as the local homozygosity, for instance). Previous work relies on numerical results or assumes Gaussian approximations for the corresponding distributions. In this paper we give an analytical derivation of the statistical properties of the local homozygosity and other empirical observables assuming selective neutrality. We find that such distributions can be very non-Gaussian.

Genetic differentiation between sympatric populations of Bacillus cereus and Bacillus thuringiensis

Little is known about genetic exchanges in natural populations of bacteria of the spore-forming Bacillus cereus group, because no population genetics studies have been performed with local sympatric populations. We isolated strains of Bacillus thuringiensis and B. cereus from small samples of soil collected at the same time from two separate geographical sites, one within the forest and the other at the edge of the forest. A total of 100 B. cereus and 98 B. thuringiensis strains were isolated and characterized by electrophoresis to determine allelic composition at nine enzymatic loci. We observed genetic differentiation between populations of B. cereus and B. thuringiensis. Populations of a given Bacillus species--B. thuringiensis or B. cereus--were genetically more similar to each other than to populations of the other Bacillus species. Hemolytic activity provided further evidence of this genetic divergence, which remained evident even if putative clones were removed from the data set. Our results suggest that the rate of gene flow was higher between strains of the same species, but that exchanges between B. cereus and B. thuringiensis were nonetheless possible. Linkage disequilibrium analysis revealed sufficient recombination for B. cereus populations to be considered panmictic units. In B. thuringiensis, the balance between clonal proliferation and recombination seemed to depend on location. Overall, our data indicate that it is not important for risk assessment purposes to determine whether B. cereus and B. thuringiensis belong to a single or two species. Assessment of the biosafety of pest control based on B. thuringiensis requires evaluation of the extent of genetic exchange between strains in realistic natural conditions.