Chemical Pollution Levels in a River Explain Site-Specific Sensitivities to Micropollutants within a Genetically Homogeneous Population of Freshwater Amphipods.

Anthropogenic micropollutants alter chemical and ecological conditions of freshwater ecosystems and impact aquatic species that live along the pollution gradient of a river. Species sensitivity to micropollutants depends on the site-specific exposure; however, it remains unclear to what degree this sensitivity relates to the species' genetic structure. Here, we explored the relationship between the toxic sensitivity and genetic structure of the amphipod species Gammarus pulex (Linnaeus, 1758) along an organic micropollutant gradient in the Holtemme River in central Germany. We determined the river's site-specific micropollutant patterns and analyzed the genetic structure of G. pulex using nuclear and mitochondrial genetic markers. Furthermore, we examined the exposure sensitivities and bioaccumulation of the commonly detected insecticide imidacloprid in G. pulex from different sites. Our results show that throughout the Holtemme River, G. pulex forms a well-connected and homogeneous population with no observable pollution-related differences in the genetic structure. However, G. pulex from polluted sites responded more sensitively to imidacloprid; survival times for half of the amphipods were up to 54% shorter, the percentage of immobile individuals increased up to 65%, and the modeled imidacloprid depuration rate was lower in comparison to amphipods from non-polluted sites. Altogether, these results suggest that the level of sensitivity of G. pulex amphipods to micropollutants in the river depends on the degree of pollution: amphipods may thrive in food-rich but polluted habitats; yet, their sensitivity is increased when chronically exposed to organic micropollutants.

RNA-Seq analysis of soft rush (Juncus effusus): transcriptome sequencing, de novo assembly, annotation, and polymorphism identification.

BACKGROUND: Juncus effusus L. (family: Juncaceae; order: Poales) is a helophytic rush growing in temperate damp or wet terrestrial habitats and is of almost cosmopolitan distribution. The species has been studied intensively with respect to its interaction with co-occurring plants as well as microbes being involved in major biogeochemical cycles. J. effusus has biotechnological value as component of Constructed Wetlands where the plant has been employed in phytoremediation of contaminated water. Its genome has not been sequenced. RESULTS: In this study we carried out functional annotation and polymorphism analysis of de novo assembled RNA-Seq data from 18 genotypes using 249 million paired-end Illumina HiSeq reads and 2.8 million 454 Titanium reads. The assembly comprised 158,591 contigs with a mean contig length of 780 bp. The assembly was annotated using the dammit! annotation pipeline, which queries the databases OrthoDB, Pfam-A, Rfam, and runs BUSCO (Benchmarking Single-Copy Ortholog genes). In total, 111,567 contigs (70.3%) were annotated with functional descriptions, assigned gene ontology terms, and conserved protein domains, which resulted in 30,932 non-redundant gene sequences. Results of BUSCO and KEGG pathway analyses were similar for J. effusus as for the well-studied members of the Poales, Oryza sativa and Sorghum bicolor. A total of 566,433 polymorphisms were identified in transcribed regions with an average frequency of 1 polymorphism in every 171 bases. CONCLUSIONS: The transcriptome assembly was of high quality and genome coverage was sufficient for global analyses. This annotated knowledge resource can be utilized for future gene expression analysis, genomic feature comparisons, genotyping, primer design, and functional genomics in J. effusus.

Isolation and characterization of eleven novel microsatellite markers for fine-scale population genetic analyses of Gammarus pulex (Crustacea: Amphipoda).

The freshwater amphipod species Gammarus pulex (Linnaeus, 1758) is widespread across Europe and Asia and is able to live in a broad range of environmental conditions. Yet, it is of great interest to which degree it is able to tolerate and adapt to the current rapid anthropogenic environmental changes affecting its habitat, such as pollution, changes in river morphology, and invasions of alien gammarid species. Microsatellite primers for genetic population studies with G. pulex have been developed but due to the existence of several genetically different lineages within the species, the application of these primers is not always successful. In order to investigate the impacts of anthropogenic environmental changes on the spatio-genetic patterns of G. pulex lineage E in streams in the Saale river catchment in Germany, we designed eleven novel polymorphic microsatellites for this lineage using a high-throughput sequencing approach. These microsatellites enabled highly specific characterization of three closely related populations. The results show genetically distinct populations reflected by both a principal coordinates analysis and an analysis of molecular variance. Several of the newly designed microsatellite primers also enabled successful cross-amplification of the respective microsatellites in specimens of G. pulex lineage C, while only two microsatellites were amplified successfully and showed polymorphisms for all of the analyzed specimens of G. fossarum Koch, 1836. The microsatellites identified here are suitable for future assessments of micro-evolutionary dynamics of G. pulex from central Germany.

Functional and phylogenetic diversity of woody plants drive herbivory in a highly diverse forest.

Biodiversity loss may alter ecosystem processes, such as herbivory, a key driver of ecological functions in species-rich (sub)tropical forests. However, the mechanisms underlying such biodiversity effects remain poorly explored, as mostly effects of species richness - a very basic biodiversity measure - have been studied. Here, we analyze to what extent the functional and phylogenetic diversity of woody plant communities affect herbivory along a diversity gradient in a subtropical forest. We assessed the relative effects of morphological and chemical leaf traits and of plant phylogenetic diversity on individual-level variation in herbivory of dominant woody plant species across 27 forest stands in south-east China. Individual-level variation in herbivory was best explained by multivariate, community-level diversity of leaf chemical traits, in combination with community-weighted means of single traits and species-specific phylodiversity measures. These findings deviate from those based solely on trait variation within individual species. Our results indicate a strong impact of generalist herbivores and highlight the need to assess food-web specialization to determine the direction of biodiversity effects. With increasing plant species loss, but particularly with the concomitant loss of functional and phylogenetic diversity in these forests, the impact of herbivores will probably decrease - with consequences for the herbivore-mediated regulation of ecosystem functions.

Tree diversity promotes functional dissimilarity and maintains functional richness despite species loss in predator assemblages.

The effects of species loss on ecosystems depend on the community's functional diversity (FD). However, how FD responds to environmental changes is poorly understood. This applies particularly to higher trophic levels, which regulate many ecosystem processes and are strongly affected by human-induced environmental changes. We analyzed how functional richness (FRic), evenness (FEve), and divergence (FDiv) of important generalist predators-epigeic spiders-are affected by changes in woody plant species richness, plant phylogenetic diversity, and stand age in highly diverse subtropical forests in China. FEve and FDiv of spiders increased with plant richness and stand age. FRic remained on a constant level despite decreasing spider species richness with increasing plant species richness. Plant phylogenetic diversity had no consistent effect on spider FD. The results contrast with the negative effect of diversity on spider species richness and suggest that functional redundancy among spiders decreased with increasing plant richness through non-random species loss. Moreover, increasing functional dissimilarity within spider assemblages with increasing plant richness indicates that the abundance distribution of predators in functional trait space affects ecological functions independent of predator species richness or the available trait space. While plant diversity is generally hypothesized to positively affect predators, our results only support this hypothesis for FD-and here particularly for trait distributions within the overall functional trait space-and not for patterns in species richness. Understanding the way predator assemblages affect ecosystem functions in such highly diverse, natural ecosystems thus requires explicit consideration of FD and its relationship with species richness.

Mycorrhizas in the Central European flora: relationships with plant life history traits and ecology.

Plant traits have been widely used to characterize different aspects of the ecology of plant species. Despite its wide distribution and its proven significance at the level of individuals, communities, and populations, the ability to form mycorrhizal associations has been largely neglected in these studies so far. Analyzing plant traits associated with the occurrence of mycorrhizas in plants can therefore enhance our understanding of plant strategies and distributions. Using a comparative approach, we tested for associations between mycorrhizal status and habitat characteristics, life history traits, and plant distribution patterns in 1752 species of the German flora (a major part of the Central European flora). Data were analyzed using log-linear models or generalized linear models, both accounting for phylogenetic relationships. Obligatorily mycorrhizal (OM) species tended to be positively associated with higher temperature, drier habitats, and higher pH; and negatively associated with moist, acidic, and fertile soils. Competitive species were more frequently OM, and stress tolerators were non-mycorrhizal (NM), while ruderal species did not show any preference. Facultatively mycorrhizal (FM) species showed the widest geographic and ecological amplitude. Indigenous species were more frequently FM and neophytes (recent aliens) more frequently OM than expected. FM species differed markedly from OM and NM species in almost all analyzed traits. Specifically, they showed a wider geographic distribution and ecological niche. Our study of the relationships between mycorrhizal status and other plant traits provides a comprehensive test of existing hypotheses and reveals novel patterns. The clear distinction between FM and OM + NM species in terms of their ecology opens up a new field of research in plant-mycorrhizal ecology.

Two common, often coexisting grassland plant species differ in their evolutionary potential in response to experimental drought.

For terrestrial plant communities, the increase in frequency and intensity of drought events is considered as one of the most severe consequences of climate change. While single-species studies demonstrate that drought can lead to relatively rapid adaptive genetic changes, the evolutionary potential and constraints to selection need to be assessed in comparative approaches to draw more general conclusions. In a greenhouse experiment, we compare the phenotypic response and evolutionary potential of two co-occurring grassland plant species, Bromus erectus and Trifolium pratense, in two environments differing in water availability. We quantified variation in functional traits and reproductive fitness in response to drought and compared multivariate genetic variance-covariance matrices and predicted evolutionary responses between species. Species showed different drought adaptation strategies, reflected in both their species-specific phenotypic plasticity and predicted responses to selection indicating contrasting evolutionary potential under drought. In T. pratense we found evidence for stronger genetic constraints under drought compared to more favourable conditions, and for some traits plastic and predicted evolutionary responses to drought had opposing directions, likely limiting the potential for adaptive change. Our study contributes to a more detailed understanding of the evolutionary potential of species with different adaptive strategies in response to climate change and may help to inform future scenarios for semi-natural grassland ecosystems.

Plant traits affecting herbivory on tree recruits in highly diverse subtropical forests.

Differences in herbivory among woody species can greatly affect the functioning of forest ecosystems, particularly in species-rich (sub)tropical regions. However, the relative importance of the different plant traits which determine herbivore damage remains unclear. Defence traits can have strong effects on herbivory, but rarely studied geographical range characteristics could complement these effects through evolutionary associations with herbivores. Herein, we use a large number of morphological, chemical, phylogenetic and biogeographical characteristics to analyse interspecific differences in herbivory on tree saplings in subtropical China. Unexpectedly, we found no significant effects of chemical defence traits. Rather, herbivory was related to the plants' leaf morphology, local abundance and climatic niche characteristics, which together explained 70% of the interspecific variation in herbivory in phylogenetic regression. Our study indicates that besides defence traits and apparency to herbivores, previously neglected measures of large-scale geographical host distribution are important factors influencing local herbivory patterns among plant species.

Identification and characterization of microsatellite loci in the rush Juncus effusus (Juncaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for the common wetland plant Juncus effusus to assess levels of within- and among-population genetic diversity. METHODS AND RESULTS: Using a next-generation sequencing approach we identified new microsatellite loci, from which 23 were further characterized in a large population of J. effusus. Expected heterozygosity ranged from 0 to 0.64 with a mean of 0.407. Inbreeding coefficients (mean = 0.984) were very high, confirming earlier suggestions of an inbred mating system. Cross-amplification was tested in J. conglomeratus, J. inflexus, J. drummondii, and J. filiformis, with best results in the former. CONCLUSIONS: The described microsatellite markers will be helpful in assisting studies on, e.g., reproduction, taxonomy, and phylogeography in J. effusus and related species.

Reduced genetic diversity of freshwater amphipods in rivers with increased levels of anthropogenic organic micropollutants.

Anthropogenic chemicals in freshwater environments contribute majorly to ecosystem degradation and biodiversity decline. In particular anthropogenic organic micropollutants (AOM), a diverse group of compounds, including pesticides, pharmaceuticals, and industrial chemicals, can significantly impact freshwater organisms. AOM were found to impact genetic diversity of freshwater species; however, to which degree AOM cause changes in population genetic structure and allelic richness of freshwater macroinvertebrates remains poorly understood. Here, the impact of AOM on genetic diversity of the common amphipod Gammarus pulex (Linnaeus, 1758) (clade E) was investigated on a regional scale. The site-specific AOM levels and their toxic potentials were determined in water and G. pulex tissue sample extracts for 34 sites along six rivers in central Germany impacted by wastewater effluents and agricultural run-off. Population genetic parameters were determined for G. pulex from the sampling sites by genotyping 16 microsatellite loci. Genetic differentiation among G. pulex from the studied rivers was found to be associated with geographic distance between sites and to differences in site-specific concentrations of AOM. The genetic diversity parameters of G. pulex were found to be related to the site-specific AOM levels. Allelic richness was significantly negatively correlated with levels of AOM in G. pulex tissue (p < 0.003) and was reduced by up to 22% at sites with increased levels of AOM, despite a positive relationship of allelic richness and the presence of waste-water effluent. In addition, the inbreeding coefficient of G. pulex from sites with toxic AOM levels was up to 2.5 times higher than that of G. pulex from more pristine sites. These results indicate that AOM levels commonly found in European rivers significantly contribute to changes in the genetic diversity of an ecologically relevant indicator species.

Pollination mode and life form strongly affect the relation between mating system and pollen to ovule ratios.

Pollen to ovule (P : O) ratios have been hypothesized to correlate with the degree of outcrossing and thus with the mating system of a plant. Also, P : O ratios are likely to vary with respect to pollination mode (i.e. wind pollination or animal pollination). Furthermore, constraints on the evolution of mating systems depending on life form may affect P : O ratios. We compiled P : O ratios and outcrossing rates for 107 angiosperm species and analyzed the relation between these traits considering pollination mode, life form and phylogenetic relatedness among species. In general, P : O ratios correlated significantly with outcrossing rates. However, when taking additional factors into account, the relation became ambiguous. The correlation was significantly positive in wind-pollinated species, but only marginally so in animal-pollinated species. Wind-pollinated species had higher P : O ratios than animal-pollinated taxa. In woody perennials, outcrossing was the predominant mating system and outcrossing rates did not correlate with P : O ratios. The results were not altered by accounting for phylogenetic relatedness among species. The results indicate that P : O ratios vary more strongly with pollination mode and life form than with the mating system.

Trait expression and signatures of adaptation in response to nitrogen addition in the common wetland plant Juncus effusus.

Wetland ecosystems are known to mitigate high nutrient loadings and thus can improve water quality and prevent potential biodiversity loss caused by eutrophication. Plant traits affect wetland processes directly through effects on accumulation or metabolization of substances, and indirectly by affecting microbial transformation processes in the soil. Understanding the causes and consequences of intraspecific variation in plant functional traits and associated ecosystem processes can aid applied ecological approaches such as wetland restoration and construction. Here we investigated molecular variation and phenotypic variation in response to three levels of nitrogen availability for a regional set of populations of the common wetland plant Juncus effusus. We asked whether trait expression reveals signatures of adaptive differentiation by comparing genetic differentiation in quantitative traits and neutral molecular markers (QST-FST comparisons) and relating trait variation to soil conditions of the plant's origin. Molecular analyses showed that samples clustered into three very distinct genetic lineages with strong population differentiation within and among lineages. Differentiation for quantitative traits was substantial but did not exceed neutral expectations when compared across treatments or for each treatment and lineage separately. However, variation in trait expression could be explained by local soil environmental conditions of sample origin, e.g. for aboveground carbon-to-nitrogen (C:N) ratios, suggesting adaptive differentiation to contribute to trait expression even at regional level.

Synchronous pulsed flowering: analysis of the flowering phenology in Juncus (Juncaceae).

BACKGROUND AND AIMS: The timing of flowering within and among individuals is of fundamental biological importance because of its influence on total seed production and, ultimately, fitness. Traditional descriptive parameters of flowering phenology focus on onset and duration of flowering and on synchrony among individuals. These parameters do not adequately account for variability in flowering across the flowering duration at individual and population level. This study aims to analyse the flowering phenology of wind-pollinated Juncus species that has been described as temporally highly variable ('pulsed flowering'). Additionally, an attempt is made to identify proximate environmental factors that may cue the flowering, and ultimate causes for the flowering patterns are discussed. METHODS: Flowering phenology was examined in populations of nine Juncus species by estimating flowering synchrony and by using the coefficient of variation (CV) to describe the temporal variation in flowering on individual and population levels. Phenologies were compared with null models to test which patterns deviate from random flowering. All parameters assessed were compared with each other and the performance of the parameters in response to randomization and varying synchrony was evaluated using a model population. Flowering patterns were correlated with temperature and humidity. KEY RESULTS: Most flowering patterns of Juncus were best described as synchronous pulsed flowering, characterized as population-wide concerted flowering events separated by days with no or few open flowers. Flowering synchrony and variability differed from a random pattern in most cases. CV values in combination with a measure of synchrony differentiated among flowering patterns found. Synchrony varied among species and was independent from variability in flowering. Neither temperature nor humidity could be determined as potential cues for the flowering pulses. CONCLUSIONS: The results indicate that selection may act independently on synchrony and variability. We propose that synchronous pulsed flowering in Juncus is an evolved strategy that provides selective benefits by increasing outcrossing and by spreading the risk of reproductive failure.

Phylogenetic turnover during subtropical forest succession across environmental and phylogenetic scales.

Although spatial and temporal patterns of phylogenetic community structure during succession are inherently interlinked and assembly processes vary with environmental and phylogenetic scales, successional studies of community assembly have yet to integrate spatial and temporal components of community structure, while accounting for scaling issues. To gain insight into the processes that generate biodiversity after disturbance, we combine analyses of spatial and temporal phylogenetic turnover across phylogenetic scales, accounting for covariation with environmental differences. We compared phylogenetic turnover, at the species- and individual-level, within and between five successional stages, representing woody plant communities in a subtropical forest chronosequence. We decomposed turnover at different phylogenetic depths and assessed its covariation with between-plot abiotic differences. Phylogenetic turnover between stages was low relative to species turnover and was not explained by abiotic differences. However, within the late-successional stages, there was high presence-/absence-based turnover (clustering) that occurred deep in the phylogeny and covaried with environmental differentiation. Our results support a deterministic model of community assembly where (i) phylogenetic composition is constrained through successional time, but (ii) toward late succession, species sorting into preferred habitats according to niche traits that are conserved deep in phylogeny, becomes increasingly important.

Trait variation in response to varying winter temperatures, diversity patterns and signatures of selection along the latitudinal distribution of the widespread grassland plant Arrhenatherum elatius.

Across Europe, genetic diversity can be expected to decline toward the North because of stochastic and selective effects which may imply diminished phenotypic variation and less potential for future genetic adaptations to environmental change. Understanding such latitudinal patterns can aid provenance selection for breeding or assisted migration approaches. In an experiment simulating different winter temperatures, we assessed quantitative trait variation, genetic diversity, and differentiation for natural populations of the grass Arrhenatherum elatius originating from a large latitudinal gradient. In general, populations from the North grew smaller and had a lower flowering probability. Toward the North, the absolute plastic response to the different winter conditions as well as heritability for biomass production significantly declined. Genetic differentiation in plant height and probability of flowering were very strong and significantly higher than under neutral expectations derived from SNP data, suggesting adaptive differentiation. Differentiation in biomass production did not exceed but mirrored patterns for neutral genetic differentiation, suggesting that migration-related processes caused the observed clinal trait variation. Our results demonstrate that genetic diversity and trait differentiation patterns for A. elatius along a latitudinal gradient are likely shaped by both local selection and genetic drift.

Separation in flowering time contributes to the maintenance of sympatric cryptic plant lineages.

Sympatric cryptic lineages are a challenge for the understanding of species coexistence and lineage diversification as well as for management, conservation, and utilization of plant genetic resources. In higher plants studies providing insights into the mechanisms creating and maintaining sympatric cryptic lineages are rare. Here, using microsatellites and chloroplast sequence data, morphometric analyses, and phenological observations, we ask whether sympatrically coexisting lineages in the common wetland plant Juncus effusus are ecologically differentiated and reproductively isolated. Our results show two genetically highly differentiated, homoploid lineages within J. effusus that are morphologically cryptic and have similar preference for soil moisture content. However, flowering time differed significantly between the lineages contributing to reproductive isolation and the maintenance of these lineages. Furthermore, the later flowering lineage suffered less from predispersal seed predation by a Coleophora moth species. Still, we detected viable and reproducing hybrids between both lineages and the earlier flowering lineage and J. conglomeratus, a coexisting close relative. Flowering time differentiation between the lineages can be explained by neutral divergence alone and together with a lack of postzygotic isolation mechanisms; the sympatric coexistence of these lineages is most likely the result of an allopatric origin with secondary contact.

Interactive effects of landscape history and current management on dispersal trait diversity in grassland plant communities.

Plant communities and their ecosystem functions are expected to be more resilient to future habitat fragmentation and deterioration if the species comprising the communities have a wide range of dispersal and persistence strategies. However, the extent to which the diversity of dispersal and persistence traits in plant communities is determined by the current and historical characteristics of sites and their surrounding landscape has yet to be explored.Using quantitative information on long-distance seed dispersal potential by wind and animals (dispersal in space) and on species' persistence/longevity (dispersal in time), we (i) compared levels of dispersal and persistence trait diversity (functional richness, FRic, and functional divergence, FDiv) in seminatural grassland plant communities with those expected by chance, and (ii) quantified the extent to which trait diversity was explained by current and historical landscape structure and local management history - taking into account spatial and phylogenetic autocorrel.Null model analysis revealed that more grassland communities than expected had a level of trait diversity that was lower or higher than predicted, given the level of species richness. Both the range (FRic) and divergence (FDiv) of dispersal and persistence trait values increased with grassland age. FDiv was mainly explained by the interaction between current grazing intensity and the amount of grassland habitat in the surrounding landscape in 1938.Synthesis. The study suggests that the variability of dispersal and persistence traits in grassland plant communities is driven by deterministic assembly processes, with both history and current management (and their interactions), playing a major role as determinants of trait diversity. While a long continuity of grazing management is likely to have promoted the diversity of dispersal and persistence traits in present-day grasslands, communities in sites that are well grazed at the present day, and were also surrounded by large amounts of grassland in the past, showed the highest diversity of dispersal and persistence strategies. Our results indicate that the historical context of a site within a landscape will influence the extent to which current grazing management is able to maintain a diversity of dispersal and persistence strategies and buffer communities (and their associated functions) against continuing habitat fragmentation.

Wolbachia infections mimic cryptic speciation in two parasitic butterfly species, Phengaris teleius and P. nausithous (Lepidoptera: Lycaenidae).

Deep mitochondrial divergence within species may result from cryptic speciation, from phylogeographic isolation or from endosymbiotic bacteria like Wolbachia that manipulate host reproduction. Phengaris butterflies are social parasites that spend most of their life in close relationship with ants. Previously, cryptic speciation has been hypothesised for two Phengaris species based on divergent mtDNA sequences. Since Phengaris species are highly endangered, the existence of cryptic species would have drastic consequences for conservation and management. We tested for cryptic speciation and alternative scenarios in P. teleius and P. nausithous based on a comprehensive sample across their Palaearctic ranges using COI gene sequences, nuclear microsatellites and tests for Wolbachia. In both species a deep mitochondrial split occurring 0.65-1.97 myrs ago was observed that did not correspond with microsatellite data but was concordant with Wolbachia infection. Haplotypes previously attributed to cryptic species were part of the Wolbachia-infected clades. In both species remaining phylogeographic structure was largely consistent between mitochondrial and nuclear genomes. In P. teleius several mitochondrial and nuclear groups were observed in East Asia while a single haplogroup and nuclear cluster prevailed across continental Eurasia. Neutrality tests suggested rapid demographic expansion into that area. In contrast, P. nausithous had several mitochondrial and nuclear groups in Europe, suggesting a complex phylogeographic history in the western part of the species range. We conclude that deep intraspecific divergences found in DNA barcode studies do not necessarily need to represent cryptic speciation but instead can be due to both infection by Wolbachia and phylogeographic structure.

Isolation by elevation: genetic structure at neutral and putatively non-neutral loci in a dominant tree of subtropical forests, Castanopsis eyrei.

BACKGROUND: The distribution of genetic diversity among plant populations growing along elevational gradients can be affected by neutral as well as selective processes. Molecular markers used to study these patterns usually target neutral processes only, but may also be affected by selection. In this study, the effects of elevation and successional stage on genetic diversity of a dominant tree species were investigated controlling for neutrality of the microsatellite loci used. METHODOLOGY/PRINCIPAL FINDINGS: Diversity and differentiation among 24 populations of Castanopsis eyrei from different elevations (251-920 m) and successional stages were analysed by eight microsatellite loci. We found that one of the loci (Ccu97H18) strongly deviated from a neutral model of differentiation among populations due to either divergent selection or hitchhiking with an unknown selected locus. The analysis showed that C. eyrei populations had a high level of genetic diversity within populations (A(R) = 7.6, H(E) = 0.82). Genetic variation increased with elevation for both the putatively selected locus Ccu97H18 and the neutral loci. At locus Ccu97H18 one allele was dominant at low elevations, which was replaced at higher elevations by an increasing number of other alleles. The level of genetic differentiation at neutral loci was similar to that of other Fagaceae species (F(ST) = 0.032,  = 0.15). Population differentiation followed a model of isolation by distance but additionally, strongly significant isolation by elevation was found, both for neutral loci and the putatively selected locus. CONCLUSIONS/SIGNIFICANCE: The results indicate higher gene flow among similar elevational levels than across different elevational levels and suggest a selective influence of elevation on the distribution of genetic diversity in C. eyrei. The study underlines the importance to check the selective neutrality of marker loci in analyses of population structure.

High selfing and high inbreeding depression in peripheral populations of Juncus atratus.

The mating system of a plant is the prime determinant of its population genetic structure. However, mating system effects may be modified by postzygotic mechanisms like inbreeding depression. Furthermore, historical as well as contemporary ecological factors and population characteristics, like the location within the species range can contribute to genetic variability. Using microsatellite markers we assessed the population genetic structure of the wind-pollinated Juncus atratus in 16 populations from peripheral and nearly central areas of the distribution range and studied the mating system of the species. In three peripheral populations, outcrossing rates at seeds stage were low (mean t(m) = 5.6%), suggesting a highly autogamous mating system. Despite this fact, on adult stage both individual heterozygosity (mean H(O) = 0.48) and gene diversity (mean H(E) = 0.58) were high even in small populations. Inbreeding coefficients were consistently low among all populations (mean F(IS) = 0.15). Within the three peripheral populations indirect estimates of lifetime inbreeding depression were surprisingly high (delta(eq) = 0.96) and inbreeding depression could be shown to act mostly on early seedling establishment. Similar conditions of autogamy combined with high inbreeding depression are typical for plants with a large lifetime genomic mutation rate that cannot avoid selfing by geitonogamy. However, the results presented here are unexpected for small-statured, herbaceous plants. Substantial genetic differentiation among all populations was found (mean F(ST) = 0.24). An isolation-by-distance pattern was apparent on large scale but not on local scale suggesting that the overall pattern was largely influenced by historical factors, e.g. colonization, whereas locally genetic drift was of greater importance than gene flow. Peripheral populations exhibited lower genetic diversity and higher inbreeding coefficients when compared with subcentral populations.