Limited genetic diversity and high differentiation among the remnant adder (Vipera berus) populations in the Swiss and French Jura Mountains

Although the adder (Vipera berus) has a large distribution area, this species is particularly threatened in Western Europe due to high habitat fragmentation and human persecution. We developed 13 new microsatellite markers in order to evaluate population structure and genetic diversity in the Swiss and French Jura Mountains, where the species is limited to only a few scattered populations. We found that V. berus exhibits a considerable genetic differentiation among populations (global F-ST = 0.269), even if these are not geographically isolated. Moreover, the genetic diversity within populations in the Jura Mountains and in the less perturbed Swiss Alps is significantly lower than in other French populations, possibly due to post-glacial recolonisation processes. Finally, in order to minimize losses of genetic diversities within isolated populations, suggestions for the conservation of this species in fragmented habitats are proposed.

Spatial variability of the chemical, physical and biological properties in lowland cultivated with irrigated rice

In the areas where irrigated rice is grown in the south of Brazil, few studies have been carried out to investigate the spatial variability structure of soil properties and to establish new forms of soil management as well as determine soil corrective and fertilizer applications. In this sense, this study had the objective of evaluating the spatial variability of chemical, physical and biological soil properties in a lowland area under irrigated rice cultivation in the conventional till system. For this purpose, a 10 x 10 m grid of 100 points was established, in an experimental field of the Embrapa Clima Temperado, in the County of Capão do Leão, State of Rio Grande do Sul. The spatial variability structure was evaluated by geostatistical tools and the number of subsamples required to represent each soil property in future studies was calculated using classical statistics. Results showed that the spatial variability structure of sand, silt, SMP index, cation exchange capacity (pH 7.0), Al3+ and total N properties could be detected by geostatistical analysis. A pure nugget effect was observed for the nutrients K, S and B, as well as macroporosity, mean weighted diameter of aggregates, and soil water storage. The cross validation procedure, based on linear regression and the determination coefficient, was more efficient to evaluate the quality of the adjusted mathematical model than the degree of spatial dependence. It was also concluded that the combination of classical with geostatistics can in many cases simplify the soil sampling process without losing information quality.

Genetic differentiation for size at first reproduction through male versus female functions in the widespread Mediterranean tree Pinus pinaster.

BACKGROUND AND AIMS: The study of local adaptation in plant reproductive traits has received substantial attention in short-lived species, but studies conducted on forest trees are scarce. This lack of research on long-lived species represents an important gap in our knowledge, because inferences about selection on the reproduction and life history of short-lived species cannot necessarily be extrapolated to trees. This study considers whether the size for first reproduction is locally adapted across a broad geographical range of the Mediterranean conifer species Pinus pinaster. In particular, the study investigates whether this monoecious species varies genetically among populations in terms of whether individuals start to reproduce through their male function, their female function or both sexual functions simultaneously. Whether differences among populations could be attributed to local adaptation across a climatic gradient is then considered. METHODS: Male and female reproduction and growth were measured during early stages of sexual maturity of a P. pinaster common garden comprising 23 populations sampled across the species range. Generalized linear mixed models were used to assess genetic variability of early reproductive life-history traits. Environmental correlations with reproductive life-history traits were tested after controlling for neutral genetic structure provided by 12 nuclear simple sequence repeat markers. KEY RESULTS: Trees tended to reproduce first through their male function, at a size (height) that varied little among source populations. The transition to female reproduction was slower, showed higher levels of variability and was negatively correlated with vegetative growth traits. Several female reproductive traits were correlated with a gradient of growth conditions, even after accounting for neutral genetic structure, with populations from more unfavourable sites tending to commence female reproduction at a lower individual size. CONCLUSIONS: The study represents the first report of genetic variability among populations for differences in the threshold size for first reproduction between male and female sexual functions in a tree species. The relatively uniform size at which individuals begin reproducing through their male function probably represents the fact that pollen dispersal is also relatively invariant among sites. However, the genetic variability in the timing of female reproduction probably reflects environment-dependent costs of cone production. The results also suggest that early sex allocation in this species might evolve under constraints that do not apply to other conifers.

Local influence for spatial analysis of soil physical properties and soybean yield using student's t-distribution

The modeling and estimation of the parameters that define the spatial dependence structure of a regionalized variable by geostatistical methods are fundamental, since these parameters, underlying the kriging of unsampled points, allow the construction of thematic maps. One or more atypical observations in the sample data can affect the estimation of these parameters. Thus, the assessment of the combined influence of these observations by the analysis of Local Influence is essential. The purpose of this paper was to propose local influence analysis methods for the regionalized variable, given that it has n-variate Student's t-distribution, and compare it with the analysis of local influence when the same regionalized variable has n-variate normal distribution. These local influence analysis methods were applied to soil physical properties and soybean yield data of an experiment carried out in a 56.68 ha commercial field in western Paraná, Brazil. Results showed that influential values are efficiently determined with n-variate Student's t-distribution.

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

The soil CO2 emission has high spatial variability because it depends strongly on soil properties. The purpose of this study was to (i) characterize the spatial variability of soil respiration and related properties, (ii) evaluate the accuracy of results of the ordinary kriging method and sequential Gaussian simulation, and (iii) evaluate the uncertainty in predicting the spatial variability of soil CO2 emission and other properties using sequential Gaussian simulations. The study was conducted in a sugarcane area, using a regular sampling grid with 141 points, where soil CO2 emission, soil temperature, air-filled pore space, soil organic matter and soil bulk density were evaluated. All variables showed spatial dependence structure. The soil CO2 emission was positively correlated with organic matter (r = 0.25, p < 0.05) and air-filled pore space (r = 0.27, p < 0.01) and negatively with soil bulk density (r = -0.41, p < 0.01). However, when the estimated spatial values were considered, the air-filled pore space was the variable mainly responsible for the spatial characteristics of soil respiration, with a correlation of 0.26 (p < 0.01). For all variables, individual simulations represented the cumulative distribution functions and variograms better than ordinary kriging and E-type estimates. The greatest uncertainties in predicting soil CO2 emission were associated with areas with the highest estimated values, which produced estimates from 0.18 to 1.85 t CO2 ha-1, according to the different scenarios considered. The knowledge of the uncertainties generated by the different scenarios can be used in inventories of greenhouse gases, to provide conservative estimates of the potential emission of these gases.

SPATIAL CORRELATION BETWEEN PHYSICAL PROPERTIES OF SOIL AND WEEDS IN TWO MANAGEMENT SYSTEMS

The spatial correlation between soil properties and weeds is relevant in agronomic and environmental terms. The analysis of this correlation is crucial for the interpretation of its meaning, for influencing factors such as dispersal mechanisms, seed production and survival, and the range of influence of soil management techniques. This study aimed to evaluate the spatial correlation between the physical properties of soil and weeds in no-tillage (NT) and conventional tillage (CT) systems. The following physical properties of soil and weeds were analyzed: soil bulk density, macroporosity, microporosity, total porosity, aeration capacity of soil matrix, soil water content at field capacity, weed shoot biomass, weed density, Commelina benghalensis density, and Bidens pilosa density. Generally, the ranges of the spatial correlations were higher in NT than in CT. The cross-variograms showed that many variables have a structure of combined spatial variation and can therefore be mapped from one another by co-kriging. This combined variation also allows inferences about the physical and biological meanings of the study variables. Results also showed that soil management systems influence the spatial dependence structure significantly.

The effect of innovation and sex-specific migration on neutral cultural differentiation

Studies of behaviour are increasingly focusing on acquisition of traits through cultural inheritance. Comparison of patterns of spatial population structure (FST) between neutral genetic loci and behavioural or cultural traits can been used to test hypotheses about demography, life history, and the mechanisms of inheritance/transmission of these traits in humans, chimpanzees and other animals. Here, we develop analytical expectations to show how FST in cultural traits can differ strongly from that measured at neutral genetic markers if migration is largely restricted to one sex but social learning is predominantly modelled on the other (e.g. males migrate, females serve as models for cultural traits), if one individual is the learning model for many, or if rates of innovation (individual learning) are high or rates of social learning are low. We discuss how comparisons of FST between genetic loci and behavioural traits can be applied to evaluate the importance of innovation in shaping patterns of cultural differentiation, as even low rates of innovation can considerably reduce FST, relative to observed structure at neutral genetic loci. Our results also suggest that differentiation in neutral cultural traits should occur over much smaller scales in species with male migration and female enculturation (or the reverse).

Sex-linked genetic markers in the shrews of the Sorex araneus group : from free gene flow to speciation

SUMMARY: The shrews of the Sorex araneus group are morphologically .very similar, but have undergone a spectacular chromosomal evolution. Altogether, the shrews of this group present a complete array of every possible level of chromosomal and genetic differentiation. In South-Western Europe, four species are recognised: S. antiriorii, S. araneus, S. coronatus and S. granarius, which differ essentially by the amount and the composition of Robertsonian metacentric chromosomés. Additionally, several chromosome races of S. araneus are also present in the same region (i.e. Bretolet, Carlit, Cordon, Jura and Vaud). The objective of this thesis was to examine the genetic relationships between populations, races and /or species of the Sorex araneus group with a special emphasis onsex-specific markers (mtDNA and Y chromosome). We first investigate the evolutionary history of the shrews of the Sorex araneus group distributed in the South-Western Europe. The results of. these analyses confirmed the difficulty to draw a single dichotomic tree within this group. Incongruent mtDNA and Y chromosome phylogenies suggest further that genetic and chromosomal evolution are in this group partially independent processes and that the evolutionary history of the south-western European populations of the S. araneus group can only be understood if we consider secondary contacts between taxa, after their divergence (with genetic exchanges by means of hybridization and / or introgression). Using one male-inherited, one female inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii). Overall there results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps after the last glacial period and came into contact in the Rhône Valley. After the Valais shrew (Sorex antinorii) reached the Swiss Alps, it came into contact with the common shrew (Sorex araneus). When two species come into contact and hybridize, endogenous counter-selection of hybrids is usually first expressed as a reduced fertility or viability in hybrids of the heterogametic sex, a mechanism know as Haldane's rule (Haldane 1922). We first evaluated the extent of introgression for Y chromosome, mtDNA and autosomal markers in a hybrid zone between S. antinoriii and S. araneus. The overall level of genetic and karyotypic differentiation between the two species must be strong .enough to allow the detection asymmetric introgression. Secondly, we compared the levels of gene flow between chromosome common to both species and chromosome differently rearranged in each of them. We detected a significantly stronger genetic structure in rearranged chromosomes. Over a 10-year period, we even observed a decrease of genetic structure for common chromosomes. These results strongly support the role of chromosomal rearrangements in the reproductive barrier between S. araneus and S. anfinorii. Overall, this thesis underlines the need to use different inherited (paternally, maternally and / or biparentally) and chromosomally located (on common vs. on rearranged chromosomes) markers to obtain more accurate pictures of genetic relationships between populations or species. RÉSUMÉ: Les musaraignes du groupe Sorex araneus sont morphologiquement très proches, mais ont connu une spectaculaire évolution chromosomique. Prises dans leur ensemble, les musaraignes de ce groupe présentent tous les nivaux possibles de différenciation génétique et chromosomique. Dans le sud-ouest de l'Europe, quatre espèces appartenant à ce groupe sont présentes : S. antinorii, S. araneus, S. coronatus et S. granarius. Celles-ci diffèrent essentiellement par leur caryotype dont la variabilité est principalement due à des fusions Robertsoniennes. De plus, plusieurs races chromosomiques appartenant à S. araneus sont aussi présentes dans la même région (i.e. les races Bretolet, Carlit, Cordon, Jura et Vaud). L'objectif de cette thèse était d'examiner les relations génétiques entre populations, races et/ou espèces du groupe S. araneus, en utilisant particulièrement des marqueurs liés aux sexes (ADN mitochondrial et Chromosome Y). Nous avons dans un premier temps retracé l'histoire évolutive des musaraignes de ce groupe dans le sud-ouest de l'Europe. Les résultats dé ces analyses confirment qu'il est difficile de tracer un simple arbre dichotomique au sein de ce groupe. Les arbres phylogénétiques obtenus sur l'ADN mitochondrial et le chromosome Y sont incongruents et suggèrent de plus que l'évolution génétique et chromosomique sont des processus indépendants. L'histoire évolutive -des populations de ce groupe ne peut. être comprise qu'en considérant des contacts secondaires entre taxa postérieure à leur divergence et induisant des échanges génétiques par hybridation et/ou introgression. Par la suite, nous avons examiné la structure génétique des populations de la musaraigne du Valais, S. antinorii, en utilisant un marqueur transmis par les mâles, un marqueur transmis par les femelles et huit marqueurs transmis par les 2 sexes. Nos résultats suggèrent que deux lignées génétiquement bien différenciées aient colonisé les Alpes Suisses, après les dernières glaciations et entrent en contact dans là Vallée du Rhône. Après avoir franchi les Alpes Suisses, la musaraigne du Valais est entrée en contact avec là musaraigne commune (S. araneus). Lorsque deux espèces entrent en contact et s'hybrident, la sélection contre les hybrides implique habituellement une baisse de fertilité ou de viabilité des hybrides du sexe hétérogamétique (i.e. les mâles XY chez les mammifères). Ce mécanisme est connu sous le nom de règle de Haldane (Haldane 1922) et implique une plus forte structuration génétique de marqueurs males - spécifiques que des marqueurs femelles spécifiques. Nous avons donc évalué le degré d'introgression des marqueurs situés sur le chromosome Y, sur l'ADN mitochondrial et sur des autosomes dans une zone hybride entre S. araneus et S. antinorii. Le niveau de différenciation chromosomique et génétique entre les 2 espèces doit être suffisamment fort pour ne pas permettre la détection d'une introgression asymétrique entre les sexes. Dans un second temps, nous avons comparé les niveaux de flux de gênes mesurés à l'échelle du chromosome, pour des chromosomes communs aux deux espèces et pour des chromosomes différemment arrangées dans chacune des deux espèces. Nous avons détecté une structure génétique significativement plus forte sur les chromosomes réarrangés et comme la zone hybride a été étudiée à dix années d'intervalle, nous observons même une diminution de la structure génétique pour les chromosomes communs au cours du temps.. Ces résultats soutiennent fortement l'hypothèse d'un rôle des réarrangements chromosomiques dans l'établissement d'une barrière reproductive entre S. araneus et S. antinorii. Ainsi cette thèse souligne l'utilité d'utiliser des marqueurs génétiques avec différents modes de transmission. (par les mâles, par les femelles et/ou par les 2 sexes) ou localisés au niveau du chromosome (chromosomes communs vs chromosomes réarrangés) afin d'obtenir une image plus juste ou du moins plus complète des relations génétiques entre populations ou espèces.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

Genetic evidence for female-biased dispersal and gene flow in a polygynous primate.

Many models of sex-biased dispersal predict that the direction of sex-bias depends upon a species' mating system. In agreement with this, almost all polygynous mammals show male-biased dispersal whereas largely monogamous birds show female-biased dispersal (FBD). The hamadryas baboon (Papio hamadryas hamadryas) is polygynous and so dispersal is predicted to be male biased, as is found in all other baboon subspecies, but there are conflicting field data showing both female and male dispersal. Using 19 autosomal genetic markers genotyped in baboons from four Saudi Arabian populations, we found strong evidence for FBD in post-dispersal adults but not, as expected, in pre-dispersal infants and young juveniles, when we compared male and female: population structure (F(st)), inbreeding (F(is)), relatedness (r), and the mean assignment index (mAIc). Furthermore, we found evidence for female-biased gene flow as population genetic structure (F(st)), was about four times higher for the paternally inherited Y, than for either autosomal markers or for maternally inherited mtDNA. These results contradict the direction of sex-bias predicted by the mating system and show that FBD has evolved recently from an ancestral state of male-biased dispersal. We suggest that the cost-benefit balance of dispersal to males and females is tightly linked to the unique hierarchical social structure of hamadryas baboons and that dispersal and social organization have coevolved.

Genetic differentiation of common shrew Sorex araneus populations among different alpine valleys revealed by microsatellites

Geographical barriers may affect the genetic structure of populations by reducing gene exchanges among them. In Switzerland, the common shrew Sorer araneus Linnaeus, 1758 is mostly confined to mountainous areas because of a competing sister species, Millet's shrew S. coronatus Millet, 1828, which occupies most of the Swiss lowlands. The structure of common shrew populations found in different alpine valleys may therefore be affected by the topography. Using microsatellites, genetic structuring of seven shrew populations is investigated among four different valleys of, the Swiss Alps. Using the exact G-test, significant genetic structuring is detected between several valleys. Isolation by distance does not fully explain our results. It appears that high mountain ridges (> 2400 m) can significantly reduce gene flow. F- and R-statistics are estimated and compared to the exact G-tests results. Mantel tests show that F-ST, unlike R-ST, is significantly correlated with differentiation. F-ST remains however low even at high differentiation levels, while R-ST has a high variance. We discuss how these results may have wider implications with regards the interpretation of microsatellite data. Finally, a new microsatellite locus, L99, appears to discriminate S. araneus of the Vaud and Cordon races from both S. araneus Valais and S. coronatus.

Genetic variability of the Brazilian hair sheep breeds

The objectives of this work were to investigate the genetic structure of the Brazilian hair sheep breeds and to determine the origin of the Santa Inês breed. Molecular similarity was determined using Randomly Amplified Polymorphic DNA - Polymerase Chain Reaction markers in 238 individuals from five naturalized sheep breeds: Santa Inês (48 animals), Rabo Largo (48), Somali (48), Morada Nova (48) and Bergamasca (46), collected in Goiás, Sergipe, Bahia, and Ceará States as well as in the Federal District. Fifty-four loci were selected from 19 primers, after a pilot test using 140 primers. Qualitative analyses indicate diagnostic markers for all breeds. All breeds were significantly different from each other. Interbreed differences were explained by 14.92% of the total variation. Santa Inês clustered with Bergamasca (97% bootstrap) and with Rabo Largo, composing the third member of the group (81% bootstrap) while Morada Nova and Somali breeds clustered separately. Each breed should be considered as a separate management and conservation unit, and special care should be taken with Rabo Largo, Morada Nova and Somali breeds, represented by small herds in Brazil.

Mites as biological tags of their hosts.

Movements and spatial distribution of host populations are expected to shape the genetic structure of their parasite populations. Comparing the genetic patterns of both interacting species may improve our understanding of their evolutionary history. Moreover, genetic analyses of parasites with horizontal transmission may serve as indicators of historical events or current demographic processes that are not apparent in the genetic signature of their hosts. Here, we compared mitochondrial variation in populations of the ectoparasitic mite Spinturnix myoti with the genetic pattern of its host, the Maghrebian bat Myotis punicus in North Africa and in the islands of Corsica and Sardinia. Mite mitochondrial differentiation among populations was correlated with both host mitochondrial and nuclear differentiation, suggesting spatial co-differentiation of the lineages of the two interacting species. Therefore our results suggest that parasite dispersal is exclusively mediated by host movements, with open water between landmasses as a main barrier for host and parasite dispersal. Surprisingly the unique presence of a continental European mite lineage in Corsica was inconsistent with host phylogeographical history and strongly suggests the former presence of European mouse-eared bats on this island. Parasites may thus act as biological tags to reveal the presence of their now locally extinct host.

Genetic characterization of ninety elite soybean cultivars using coefficient of parentage

The objective of this work was to estimate the coefficient of parentage and to understand the genetic structure of 90 elite soybean cultivars, which are adapted to different Brazilian environments. A total of 4,005 coefficients of parentage (f) were obtained and used to group the cultivars by UPGMA method. The constructed dendrogram showed several related cultivar groups which shared similar ancestors and clearly showed the genetic structure of the main Brazilian cultivars. Effective population sizes (Ne) were also estimated for cultivars in different generations. The average f = 0.2124 value, obtained from cultivars classified into four decades according to the release year, suggested effective soybean population sizes of 11 and 13 calculated using arithmetic and weighted means, respectively. The relatively small Ne and the high parentage coefficient support the conclusion that there is a high similarity degree among the main soybean cultivars in Brazil.

Ecology and life history affect different aspects of the population structure of 27 high-alpine plants.

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics.