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.

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis).

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.

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.

Genes mirror geography within Europe.

Understanding the genetic structure of human populations is of fundamental interest to medical, forensic and anthropological sciences. Advances in high-throughput genotyping technology have markedly improved our understanding of global patterns of human genetic variation and suggest the potential to use large samples to uncover variation among closely spaced populations. Here we characterize genetic variation in a sample of 3,000 European individuals genotyped at over half a million variable DNA sites in the human genome. Despite low average levels of genetic differentiation among Europeans, we find a close correspondence between genetic and geographic distances; indeed, a geographical map of Europe arises naturally as an efficient two-dimensional summary of genetic variation in Europeans. The results emphasize that when mapping the genetic basis of a disease phenotype, spurious associations can arise if genetic structure is not properly accounted for. In addition, the results are relevant to the prospects of genetic ancestry testing; an individual's DNA can be used to infer their geographic origin with surprising accuracy-often to within a few hundred kilometres.

Characterization of 13 polymorphic microsatellite loci for two land snail species, Trochulus villosus and T. sericeus (Gastropoda: Pulmonata: Hygromiidae).

Thirteen new microsatellite loci were isolated and tested on two land snail species, Trochulus villosus and T. sericeus (Pulmonata: Hygromiidae), resulting in a set of eight polymorphic markers for each species. The expected heterozygosity was high for all loci and species (between 0.616 and 0.944). Such levels of variability will allow detailed insights into the population genetic structure of some Trochulus species.

Staphylococcus aureus clfB and spa alleles of the repeat regions are segregated into major phylogenetic lineages.

To reliably differentiate among Staphylococcus aureus isolates we recently developed the Double Locus Sequence Typing (DLST) based on the analysis of partial sequences of clfB and spa genes. This method is highly discriminatory and gives unambiguous definition of types. The highly clonal population structure of S. aureus suggests that isolates with identical clfB or spa alleles belong to the same clonal complex (CC) defined by Multi-Locus Sequence Typing (MLST). To test this hypothesis as well as to investigate putative intra-CC genetic structure, we analyzed a total of 289 isolates (186 MSSA and 103 MRSA) with DLST-, spa- and MLST-typing. Among the 289 strains, 242 were clustered into 7 major MLST CCs, 40 into minor CCs and 7 were not grouped into CCs. A total of 205 DLST- and 129 spa-types were observed. With one exception, all DLST-clfB, DLST-spa and spa-type alleles were segregated into CCs. DLST-types sharing an identical allele (clfB or spa) were clustered using eBURST. Except for one strain, all isolates from each DLST cluster belonged to the same CC. However, using both DLST- and spa-typing we were not able to disclose a clear intra-CC structure. Nevertheless, the high diversity of these loci confirmed that they are good markers for local epidemiological investigations.

Ecology and evolution of social-organization: insights from fire ants and other highly eusocial insects

Social organisms exhibit conspicuous intraspecific variation in all facets of their social organization. A prominent example of such variation in the highly eusocial Hymenoptera is differences in the number of reproductive queens per colony, Differences in queen number in ants are associated with differences in a host of reproductive and social traits, including queen phenotype and breeding strategy, mode of colony reproduction, and pattern of sex allocation. We examine the causes and consequences of changes in colony queen number and associated traits using the fire ant Solenopsis invicta as a principal model. Ecological constraints on mode of colony founding may act as important selective forces causing the evolution of queen number in this and many other ants, with social organization generally perpetuated across generations by means of the social environment molding appropriate queen phenotypes and reproductive strategies. Shifts in colony queen number have profound effects on genetic structure within nests and may also influence genetic structure at higher levels (aggregations of nests or local demes) because of the association of queen number with particular mating and dispersal habits. Divergence of breeding habits between populations with different social organizations has the potential to promote genetic differentiation between these social variants. Thus, evolution of social organization can be important in generating intrinsic selective regimes that channel subsequent social evolution and in initiating the development of significant population genetic structure, including barriers to gene flow important in cladogenesis.

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.

Diagnostic techniques applied in geostatistics for agricultural data analysis

The structural modeling of spatial dependence, using a geostatistical approach, is an indispensable tool to determine parameters that define this structure, applied on interpolation of values at unsampled points by kriging techniques. However, the estimation of parameters can be greatly affected by the presence of atypical observations in sampled data. The purpose of this study was to use diagnostic techniques in Gaussian spatial linear models in geostatistics to evaluate the sensitivity of maximum likelihood and restrict maximum likelihood estimators to small perturbations in these data. For this purpose, studies with simulated and experimental data were conducted. Results with simulated data showed that the diagnostic techniques were efficient to identify the perturbation in data. The results with real data indicated that atypical values among the sampled data may have a strong influence on thematic maps, thus changing the spatial dependence structure. The application of diagnostic techniques should be part of any geostatistical analysis, to ensure a better quality of the information from thematic maps.

O-chromosome lethal frequencies in Serbian and Montenegrin Drosophila subobscura populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.

Chromosome localization of microsatellite markers in the shrews of the Sorex araneus group.

The extremely high rate of karyotypic evolution that characterizes the shrews of the Sorex araneus group makes this group an exceptionally interesting model for population genetics and evolutionary studies. Here, we attempted to map 46 microsatellite markers at the chromosome arm level using flow-sorted chromosomes from three karyotypically different taxa of the Sorex araneus group (S. granarius and the chromosome races Cordon and Novosibirsk of S. araneus). The most likely localizations were provided for 35 markers, among which 25 were each unambiguously mapped to a single locus on the corresponding chromosomes in the three taxa, covering the three sexual chromosomes (XY1Y2) and nine of the 18 autosomal arms of the S. araneus group. The results provide further evidence for a high degree of conservation in genome organization in the S. araneus group despite the presence of numerous Robertsonian rearrangements. These markers can therefore be used to compare the genetic structure among taxa of the S. araneus group at the chromosome level and to study the role of chromosomal rearrangements in the genetic diversification and speciation process of this group.

Isolation and characterization of highly polymorphic microsatellite loci in the bladder campion, Silene vulgaris (Caryophyllaceae)

This study reports the isolation and characterization of seven highly polymorphic microsatellite loci in Silene vulgaris (Caryophyllaceae). The loci were isolated from two libraries constructed from genomic DNA enriched for CA and GA repeats. These markers yielded nine to 40 alleles per locus (mean 22.1) in a survey of 45 individuals from a single population located in the western Swiss Alps. Average observed heterozygosity ranged from 16.2 to 77.4%. These microsatellite loci should be valuable tools for studying fine-scale genetic structure.

O-chromosome lethal frequencies in Serbian and Montenegrin Drosophila subobscura populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.

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.

Sex-specific estimates of dispersal show female philopatry and male dispersal in a promiscuous amphibian, the alpine salamander (Salamandra atra).

Amphibians display wide variations in life-history traits and life cycles that should prove useful to explore the evolution of sex-biased dispersal, but quantitative data on sex-specific dispersal patterns are scarce. Here, we focused on Salamandra atra, an endemic alpine species showing peculiar life-history traits. Strictly terrestrial and viviparous, the species has a promiscuous mating system, and females reproduce only every 3 to 4 years. In the present study, we provide quantitative estimates of asymmetries in male vs. female dispersal using both field-based (mark-recapture) and genetic approaches (detection of sex-biased dispersal and estimates of migration rates based on the contrast in genetic structure across sexes and age classes). Our results revealed a high level of gene flow among populations, which stems exclusively from male dispersal. We hypothesize that philopatric females benefit from being familiar with their natal area for the acquisition and defence of an appropriate shelter, while male dispersal has been secondarily favoured by inbreeding avoidance. Together with other studies on amphibians, our results indicate that a species' mating system alone is a poor predictor of sex-linked differences in dispersal, in particular for promiscuous species. Further studies should focus more directly on the proximate forces that favour or limit dispersal to refine our understanding of the evolution of sex-biased dispersal in animals.