Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Bioclimatic constraints to Andean cat distribution: a modelling application for rare species

AimTo identify the bioclimatic niche of the endangered Andean cat (Leopardus jacobita), one of the rarest and least known felids in the world, by developing a species distribution model.LocationSouth America, High Andes and Patagonian steppe. Peru, Bolivia, Chile, Argentina.MethodsWe used 108 Andean cat records to build the models, and 27 to test them, applying the Maxent algorithm to sets of uncorrelated bioclimatic variables from global databases, including elevation. We based our biogeographical interpretations on the examination of the predicted geographic range, the modelled response curves and latitudinal variations in climatic variables associated with the locality data.ResultsSimple bioclimatic models for Andean cats were highly predictive with only 3-4 explanatory variables. The climatic niche of the species was defined by extreme diurnal variations in temperature, cold minimum and moderate maximum temperatures, and aridity, characteristic not only of the Andean highlands but also of the Patagonian steppe. Argentina had the highest representation of suitable climates, and Chile the lowest. The most favourable conditions were centrally located and spanned across international boundaries. Discontinuities in suitable climatic conditions coincided with three biogeographical barriers associated with climatic or topographic transitions.Main conclusionsSimple bioclimatic models can produce useful predictions of suitable climatic conditions for rare species, including major biogeographical constraints. In our study case, these constraints are also known to affect the distribution of other Andean species and the genetic structure of Andean cat populations. We recommend surveys of areas with suitable climates and no Andean cat records, including the corridor connecting two core populations. The inclusion of landscape variables at finer scales, crucially the distribution of Andean cat prey, would contribute to refine our predictions for conservation applications.

Power transformations in correspondence analysis

Power transformations of positive data tables, prior to applying the correspondence analysis algorithm, are shown to open up a family of methods with direct connections to the analysis of log-ratios. Two variations of this idea are illustrated. The first approach is simply to power the original data and perform a correspondence analysis this method is shown to converge to unweighted log-ratio analysis as the power parameter tends to zero. The second approach is to apply the power transformation to thecontingency ratios, that is the values in the table relative to expected values based on the marginals this method converges to weighted log-ratio analysis, or the spectral map. Two applications are described: first, a matrix of population genetic data which is inherently two-dimensional, and second, a larger cross-tabulation with higher dimensionality, from a linguistic analysis of several books.

Evolution in subdivided plant populations: concepts, recent advances and future directions.

SUMMARY: Research into the evolution of subdivided plant populations has long involved the study of phenotypic variation across plant geographic ranges and the genetic details underlying that variation. Genetic polymorphism at different marker loci has also allowed us to infer the long- and short-term histories of gene flow within and among populations, including range expansions and colonization-extinction dynamics. However, the advent of affordable genome-wide sequences for large numbers of individuals is opening up new possibilities for the study of subdivided populations. In this review, we consider what the new tools and technologies may allow us to do. In particular, we encourage researchers to look beyond the description of variation and to use genomic tools to address new hypotheses, or old ones afresh. Because subdivided plant populations are complex structures, we caution researchers away from adopting simplistic interpretations of their data, and to consider the patterns they observe in terms of the population genetic processes that have given rise to them; here, the genealogical framework of the coalescent will continue to be conceptually and analytically useful.

Local acquisition of methicillin-resistance in predominant Staphylococcus aureus clones of western Switzerland.

Objectives: Recent population genetic studies suggest that the Staphylococcal Chromosome Cassettes mec (SCCmec) was acquired at a global scale much more frequently than previously thought. We hypothesized that such acquisitions can also be observed at a local level. In the present study, we aimed at investigating the diversity of SCCmec in a local MRSA population, where the dissemination of four MRSA clones has been observed (JCM 2007, 45: 3729). Methods: All the MRSA isolates (one per patient) recovered in the Vaud canton of Switzerland from January 2005 to December 2008 were analyzed in this study. We used the Double Locus Sequence Typing (DLST) method, based on clfB and spa loci, and the e-BURST algorithm to group the types with one allele in common (i.e. clone). To increase the discriminatory power of the DLST method, a third polymorphic marker (clfA) was further analyzed on a sub-sample of isolates. The SCCmec type of each isolate was determined with the first two PCRs of the Kondo scheme. Results: DLST analysis indicated that 1884/2036 isolates (92.5%) belong to the four predominant clones. A majority of isolates in each clone harboured an identical SCCmec type: 61/64 (95%) isolates to DLST clone 1−1 SCCmec IV, 1282/1323 (97%) to clone 2−2 SCCmec II, 237/288 (82%) to clone 3−3 SCCmec IV, and 192/209 (92%) to clone 4−4 SCCmec I. Unexpectedly, different SCCmec types were present in a single predominant DLST clone: SCCmec V plus one unusual type in 3 isolates of clone 1−1; SCCmec I, IV, V, VI plus two unusual types in 41 isolates of clone 2−2; SCCmec I, II, VI plus three unusual types in 51 isolates of clone 3−3; and SCCmec II, IV, V plus one unusual type in 17 isolates of clone 4−4. Interestingly, adding a third locus generally did not change the classification of incongruent SCCmec types, suggesting that these SCCmec elements have been acquired locally during the dissemination of the clones. Conclusion: Although the SCCmec diversity within clones was relatively low at a local level, a significant proportion of isolates with different SCCmec have been identified in the four major clones. This suggests that the local acquisition of SCCmec elements is not a rare event and illustrates the great capacity of S. aureus to quickly adapt to its environment by acquiring new genetic elements.

The interplay between relatedness and horizontal gene transfer drives the evolution of plasmid-carried public goods.

Plasmids carry a wide range of genes that are often involved in bacterial social behaviour. The question of why such genes are frequently mobile has received increasing attention. Here, we use an explicit population genetic approach to model the evolution of plasmid-borne bacterial public goods production. Our findings highlight the importance of both transmission and relatedness as factors driving the evolution of plasmid-borne public goods production. We partition the effects of plasmid transfer of social traits into those of infectivity and the effect of increased relatedness. Our results demonstrate that, owing to its effect on relatedness, plasmid mobility increases the invasion and stability of public goods, in a way not seen in individually beneficial traits. In addition, we show that plasmid transfer increases relatedness when public goods production is rare but this effect declines when production is common, with both scenarios leading to an increase in the frequency of plasmid-borne public goods. Plasmids remain important vectors for the spread of social genes involved in bacterial virulence thus an understanding of their dynamics is highly relevant from a public health perspective.

Climate oscillations and species interactions: large-scale congruence but regional differences in the phylogeographic structures of an alpine plant and its monophagous insect

Aim. To predict the fate of alpine interactions involving specialized species, using a monophagous beetle and its host-plant as a case study. Location. The Alps. Methods. We investigated genetic structuring of the herbivorous beetle Oreina gloriosa and its specific host-plant Peucedanum ostruthium. We used genome fingerprinting (in the insect and the plant) and sequence data (in the insect) to compare the distribution of the main gene pools in the two associated species and to estimate divergence time in the insect, a proxy for the temporal origin of the interaction. We quantified the similarity in spatial genetic structures by performing a Procrustes analysis, a tool from the shape theory. Finally, we simulated recolonization of an empty space analogous to the deglaciated Alps just after ice retreat by two lineages from two species showing unbalanced dependence, to examine how timing of the recolonization process, as well as dispersal capacities of associated species, could explain the observed pattern. Results. Contrasting with expectations based on their asymmetrical dependence, patterns in the beetle and plant were congruent at a large scale. Exceptions occurred at a regional scale in areas of admixture, matching known suture zones in Alpine plants. Simulations using a lattice-based model suggested these empirical patterns arose during or soon after recolonization, long after the estimated origin of the interaction c. 0.5 million years ago. Main conclusions. Species-specific interactions are scarce in alpine habitats because glacial cycles have limited opportunities for coevolution. Their fate, however, remains uncertain under climate change. Here we show that whereas most dispersal routes are paralleled at large scale, regional incongruence implies that the destinies of the species might differ under changing climate. This may be a consequence of the host-dependence of the beetle that locally limits the establishment of dispersing insects.

Bacteria diversity and microbial biomass in forest, pasture and fallow soils in the southwestern Amazon basin

It is well-known that Amazon tropical forest soils contain high microbial biodiversity. However, anthropogenic actions of slash and burn, mainly for pasture establishment, induce profound changes in the well-balanced biogeochemical cycles. After a few years the grass yield usually declines, the pasture is abandoned and is transformed into a secondary vegetation called "capoeira" or fallow. The aim of this study was to examine how the clearing of Amazon rainforest for pasture affects: (1) the diversity of the Bacteria domain evaluated by Polymerase Chain Reaction and Denaturing Gradient Gel Electrophoresis (PCR-DGGE), (2) microbial biomass and some soil chemical properties (pH, moisture, P, K, Ca, Mg, Al, H + Al, and BS), and (3) the influence of environmental variables on the genetic structure of bacterial community. In the pasture soil, total carbon (C) was between 30 to 42 % higher than in the fallow, and almost 47 % higher than in the forest soil over a year. The same pattern was observed for N. Microbial biomass in the pasture was about 38 and 26 % higher than at fallow and forest sites, respectively, in the rainy season. DGGE profiling revealed a lower number of bands per area in the dry season, but differences in the structure of bacterial communities among sites were better defined than in the wet season. The bacterial DNA fingerprints in the forest were stronger related to Al content and the Cmic:Ctot and Nmic:Ntot ratios. For pasture and fallow sites, the structure of the Bacteria domain was more associated with pH, sum of bases, moisture, total C and N and the microbial biomass. In general microbial biomass in the soils was influenced by total C and N, which were associated with the Bacteria domain, since the bacterial community is a component and active fraction of the microbial biomass. Results show that the genetic composition of bacterial communities in Amazonian soils changed along the sequence forest-pasture-fallow.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species.

BACKGROUND: Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications. FINDINGS: We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement. CONCLUSIONS: The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

Origin and expansion of the allotetraploid Aegilops geniculata, a wild relative of wheat.

*This study reconstructs the phylogeography of Aegilops geniculata, an allotetraploid relative of wheat, to discuss the impact of past climate changes and recent human activities (e.g. the early expansion of agriculture) on the genetic diversity of ruderal plant species. *We combined chloroplast DNA (cpDNA) sequencing, analysed using statistical parsimony network, with nonhierarchical K-means clustering of amplified fragment length polymorphism (AFLP) genotyping, to unravel patterns of genetic structure across the native range of Ae. geniculata. The AFLP dataset was further explored by measurement of the regional genetic diversity and the detection of isolation by distance patterns. *Both cpDNA and AFLP suggest an eastern Mediterranean origin of Ae. geniculata. Two lineages have spread independently over northern and southern Mediterranean areas. Northern populations show low genetic diversity but strong phylogeographical structure among the main peninsulas, indicating a major influence of glacial cycles. By contrast, low genetic structuring and a high genetic diversity are detected in southern Mediterranean populations. Finally, we highlight human-mediated dispersal resulting in substantial introgression between resident and migrant populations. *We have shown that the evolutionary trajectories of ruderal plants can be similar to those of wild species, but are interfered by human activities, promoting range expansions through increased long-distance dispersal and the creation of suitable habitats.

Mating frequency of ant queens with alternative dispersal strategies, as revealed by microsatellite analysis of sperm

In social Hymenoptera (ants, bees, and wasps), the number of males that mate with the same queen affects social and genetic organization of the colony. However, the selective forces leading to single mating in certain conditions and multiple mating in others remain enigmatic. In this study, I investigated whether queens of the wood ant Formica paralugubris adopting different dispersal strategies varied in their mating frequency (the number of males with whom they mated). The frequency of multiple mating was determined by using microsatellite markers to genotype the sperm stored in the spermatheca of queens, and the validity of this method was confirmed by analysing mother-offspring combinations obtained from experimental single-queen colonies. Dispersing queens, which may found new colonies, did not mate with more males than queens that stayed within polygynous colonies, where the presence of numerous reproductive individuals ensured a high level of genetic diversity. Hence, this study provides no support to the hypotheses that multiple mating is beneficial because it increases genetic variability within colonies. Most of the F. paralugubris queens mated with a single male, whatever their dispersal strategy and life history. Moreover, multiple mating had little effect on colony genetic structure: the effective mating frequency was 1.11 when calculated from within-brood relatedness, and 1.13 when calculated from the number of mates detected in the sperm. Hence, occasional multiple mating by F. paralugubris queens may have no adaptive significance.

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones.

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.

Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae).

BackgroundGenetic, phenotypic and ecological divergence within a lineage is the result of past and ongoing evolutionary processes, which lead ultimately to diversification and speciation. Integrative analyses allow linking diversification to geological, climatic, and ecological events, and thus disentangling the relative importance of different evolutionary drivers in generating and maintaining current species richness.ResultsHere, we use phylogenetic, phenotypic, geographic, and environmental data to investigate diversification in the Spanish sand racer (Psammodromus hispanicus). Phylogenetic, molecular clock dating, and phenotypic analyses show that P. hispanicus consists of three lineages. One lineage from Western Spain diverged 8.3 (2.9-14.7) Mya from the ancestor of Psammodromus hispanicus edwardsianus and P. hispanicus hispanicus Central lineage. The latter diverged 4.8 (1.5-8.7) Mya. Molecular clock dating, together with population genetic analyses, indicate that the three lineages experienced northward range expansions from southern Iberian refugia during Pleistocene glacial periods. Ecological niche modelling shows that suitable habitat of the Western lineage and P. h. edwardsianus overlap over vast areas, but that a barrier may hinder dispersal and genetic mixing of populations of both lineages. P. h. hispanicus Central lineage inhabits an ecological niche that overlaps marginally with the other two lineages.ConclusionsOur results provide evidence for divergence in allopatry and niche conservatism between the Western lineage and the ancestor of P. h. edwardsianus and P. h. hispanicus Central lineage, whereas they suggest that niche divergence is involved in the origin of the latter two lineages. Both processes were temporally separated and may be responsible for the here documented genetic and phenotypic diversity of P. hispanicus. The temporal pattern is in line with those proposed for other animal lineages. It suggests that geographic isolation and vicariance played an important role in the early diversification of the group, and that lineage diversification was further amplified through ecological divergence.

The evolution and consequences of sex-specific reproductive variance.

Natural selection favors alleles that increase the number of offspring produced by their carriers. But in a world that is inherently uncertain within generations, selection also favors alleles that reduce the variance in the number of offspring produced. If previous studies have established this principle, they have largely ignored fundamental aspects of sexual reproduction and therefore how selection on sex-specific reproductive variance operates. To study the evolution and consequences of sex-specific reproductive variance, we present a population-genetic model of phenotypic evolution in a dioecious population that incorporates previously neglected components of reproductive variance. First, we derive the probability of fixation for mutations that affect male and/or female reproductive phenotypes under sex-specific selection. We find that even in the simplest scenarios, the direction of selection is altered when reproductive variance is taken into account. In particular, previously unaccounted for covariances between the reproductive outputs of different individuals are expected to play a significant role in determining the direction of selection. Then, the probability of fixation is used to develop a stochastic model of joint male and female phenotypic evolution. We find that sex-specific reproductive variance can be responsible for changes in the course of long-term evolution. Finally, the model is applied to an example of parental-care evolution. Overall, our model allows for the evolutionary analysis of social traits in finite and dioecious populations, where interactions can occur within and between sexes under a realistic scenario of reproduction.