Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

Background Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest. Methods Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea. Results We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics. Conclusions The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

Climatic change and contemporaneous land-use phases north and south of the Alps 2300 BC to 800 AD

Fluctuations in the Δ14C curve and subsequent gaps of archaeological findings at 800–650 and 400–100 BC in western and central Europe may indicate major climate-driven land-abandonment phases. To address this hypothesis radiocarbon-dated sediments from four lakes in Switzerland were studied palynologically. Pollen analysis indicates contemporaneous phases of forest clearances and of intensified land-use at 1450–1250 BC, 650–450 BC, 50 BC–100 AD and around 700 AD. These land-use expansions coincided with periods of warm climate as recorded by the Alpine dendroclimatic and Greenland oxygen isotope records. Our results suggest that harvest yields would have increased synchronously over wide areas of central and southern Europe during periods of warm and dry climate. Combined interpretation of palaeoecological and archaeological findings suggests that higher food production led to increased human populations. Positive long-term trends in pollen values of Cerealia and Plantago lanceolata indicate that technical innovations during the Bronze and Iron Age (e.g. metal ploughs, scythes, hay production, fertilising methods) gradually increased agricultural productivity. The successful adoption of yield-increasing advances cannot be explained by climatic determinism alone. Combined with archaeological evidence, our results suggest that despite considerable cycles of spatial and demographic reorganisation (repeated land abandonments and expansions, as well as large-scale migrations and population decreases), human societies were able to shift to lower subsistence levels without dramatic ruptures in material culture. However, our data imply that human societies were not able to compensate rapidly for harvest failures when climate deteriorated. Agriculture in marginal areas was abandoned, and spontaneous reforestations took place on abandoned land south and north of the Alps. Only when the climate changed again to drier and warmer conditions did a new wide-spread phase of forest clearances and field extensions occur, allowing the reoccupation of previously abandoned areas. Spatial distribution of cereal cultivation and growth requirements of Cerealia species suggest that increases in precipitation were far more decisive in driving crop failures over central and southern Europe than temperature decreases.

Impact of oceanic circulation changes on atmospheric δ¹³ CO₂

δ¹³ CO₂ measured in Antarctic ice cores provides constraints on oceanic and terrestrial carbon cycle processes linked with millennial-scale changes in atmospheric CO₂. However, the interpretation of δ¹³ CO₂ is not straight-forward. Using carbon isotope-enabled versions of the LOVECLIM and Bern3D models, we perform a set of sensitivity experiments in which the formation rates of North Atlantic Deep Water (NADW), North Pacific Deep Water (NPDW), Antarctic Bottom Water (AABW), and Antarctic Intermediate Water (AAIW) are varied. We study the impact of these circulation changes on atmospheric δ¹³ CO₂ as well as on the oceanic δ¹³ CO₂ distribution. In general, we find that the formation rates of AABW, NADW, NPDW, and AAIW are negatively correlated with changes in δ¹³ CO₂: namely, strong oceanic ventilation decreases atmospheric δ¹³ CO₂. However, since large-scale oceanic circulation reorganizations also impact nutrient utilization and the Earth’s climate, the relationship between atmospheric δ¹³ CO₂ levels and ocean ventilation rate is not unequivocal. In both models atmospheric δ¹³ CO₂ is very sensitive to changes in AABW formation rates: increased AABW formation enhances the transport of low δ¹³ CO₂ waters to the surface and decreases atmospheric δ¹³ CO₂. By contrast, the impact of NADW changes on atmospheric δ¹³ CO₂ is less robust and might be model dependent. This results from complex interplay between global climate, carbon cycle, and the formation rate of NADW, a water body characterized by relatively high δ¹³ CO₂.

Millennial to orbital-scale variations of drought intensity in the Eastern Mediterranean

Millennial to orbital-scale rainfall changes in the Mediterranean region and corresponding variations in vegetation patterns were the result of large-scale atmospheric reorganizations. In spite of recent efforts to reconstruct this variability using a range of proxy archives, the underlying physical mechanisms have remained elusive. Through the analysis of a new high-resolution sedimentary section from Lake Van (Turkey) along with climate modeling experiments, we identify massive droughts in the Eastern Med- iterranean for the past four glacial cycles, which have a pervasive link with known intervals of enhanced North Atlantic glacial iceberg calving, weaker Atlantic Meridional Overturning Circulation and Dansgaard-Oeschger cold conditions. On orbital timescales, the topographic effect of large Northern Hemisphere ice sheets and periods with minimum insolation seasonality further exacerbated drought intensities by suppressing both summer and winter precipitation.

Stand scale variability of topsoil organic matter composition in a high-elevation Norway spruce forest ecosystem

Our knowledge about the effect of single-tree influence areas on the physicochemical properties of the underlying mineral soil in forest ecosystems is still limited. This restricts our ability to adequately estimate future changes in soil functioning due to forest management practices. We studied the stand scale spatial variation of different soil organic matter species investigated by 13C NMR spectroscopy, lignin phenol and neutral sugar analysis under an unmanaged mountainous high-elevation Norway spruce (Picea abies L.) forest in central Europe. Multivariate geostatistical approaches were applied to relate the spatial patterns of the different soil organic matter species to topographic parameters, bulk density, oxalate- and dithionite-extractable iron, pH, and the impact of tree distribution. Soil samples were taken from the mineral top soil. Generally, the stand scale distribution patterns of different soil organic matter compounds could be divided into two groups: Those compounds, which were significantly spatially correlated with topography/altitude and those with small scale spatial pattern (range ≤ 10 m) that was closely related to tree distribution. The concentration of plant-derived soil organic matter components, such as lignin, at a given sampling point was significantly spatially related to the distance of the nearest tree (p ≤ 0.05). In contrast, the spatial distribution of mainly microbial-derived compounds (e.g. galactose and mannose) could be attributed to the dominating impact of small-scale topography and the contribution of poorly crystalline iron oxides that were significantly larger in the central depression of the study site compared to crest and slope positions. Our results demonstrate that topographic parameters dominate the distribution of overall topsoil organic carbon (OC) stocks at temperate high-elevation forest ecosystems, particularly in sloped terrain. However, trees superimpose topography-controlled OC biogeochemistry beneath their crown by releasing litter and changing soil conditions in comparison to open areas. This may lead to distinct zones with different mechanisms of soil organic matter degradation and also stabilization in forest stands.

Complete genome of a Puumala virus strain from Central Europe

Puumala virus (PUUV) is one of the predominant hantavirus species in Europe causing mild to moderate cases of haemorrhagic fever with renal syndrome. Parts of Lower Saxony in north-western Germany are endemic for PUUV infections. In this study, the complete PUUV genome sequence of a bank vole-derived tissue sample from the 2007 outbreak was determined by a combined primer-walking and RNA ligation strategy. The S, M and L genome segments were 1,828, 3,680 and 6,550 nucleotides in length, respectively. Sliding-window analyses of the nucleotide sequences of all available complete PUUV genomes indicated a non-homogenous distribution of variability with hypervariable regions located at the 3′-ends of the S and M segments. The overall similarity of the coding genome regions to the other PUUV strains ranged between 80.1 and 84.7 % at the level of the nucleotide sequence and between 89.5 and 98.1 % for the deduced amino acid sequences. In comparison to the phylogenetic trees of the complete coding sequences, trees based on partial segments revealed a general drop in phylogenetic support and a lower resolution. The Astrup strain S and M segment sequences showed the highest similarity to sequences of strains from geographically close sites in the Osnabrück Hills region. In conclusion, a primer-walking-mediated strategy resulted in the determination of the first complete nucleotide sequence of a PUUV strain from Central Europe. Different levels of variability along the genome provide the opportunity to choose regions for analyses according to the particular research question, e.g., large-scale phylogenetics or within-host evolution.

Biotic homogenization can decrease landscape-scale forest multifunctionality

Many experiments have shown that local biodiversity loss impairs the ability of ecosystems to maintain multiple ecosystem functions at high levels (multifunctionality). In contrast, the role of biodiversity in driving ecosystem multifunctionality at landscape scales remains unresolved. We used a comprehensive pan-European dataset, including 16 ecosystem functions measured in 209 forest plots across six European countries, and performed simulations to investigate how local plot-scale richness of tree species (α-diversity) and their turnover between plots (β-diversity) are related to landscape-scale multifunctionality. After accounting for variation in environmental conditions, we found that relationships between α-diversity and landscape-scale multifunctionality varied from positive to negative depending on the multifunctionality metric used. In contrast, when significant, relationships between β-diversity and landscape-scale multifunctionality were always positive, because a high spatial turnover in species composition was closely related to a high spatial turnover in functions that were supported at high levels. Our findings have major implications for forest management and indicate that biotic homogenization can have previously unrecognized and negative consequences for large-scale ecosystem multifunctionality.