Refine
Year of publication
Language
- English (37)
Is part of the Bibliography
- yes (37)
Keywords
- Larix (4)
- climate change (3)
- metabarcoding (3)
- pollen (3)
- sedimentary ancient DNA (3)
- treeline (3)
- Environmental DNA (2)
- Holocene (2)
- Holocene environmental history (2)
- Larix gmelinii (2)
Institute
- Institut für Biochemie und Biologie (37) (remove)
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
Pollen-based quantitative reconstructions of past climate variables is a standard palaeoclimatic approach. Despite knowing that the spatial extent of the calibration-set affects the reconstruction result, guidance is lacking as to how to determine a suitable spatial extent of the pollen-climate calibration-set. In this study, past mean annual precipitation (P-ann) during the Holocene (since 11.5 cal ka BP) is reconstructed repeatedly for pollen records from Qinghai Lake (36.7 degrees N, 100.5 degrees E; north-east Tibetan Plateau), Gonghai Lake (38.9 degrees N, 112.2 degrees E; north China) and Sihailongwan Lake (42.3 degrees N, 126.6 degrees E; north-east China) using calibration-sets of varying spatial extents extracted from the modern pollen dataset of China and Mongolia (2559 sampling sites and 168 pollen taxa in total). Results indicate that the spatial extent of the calibration-set has a strong impact on model performance, analogue quality and reconstruction diagnostics (absolute value, range, trend, optimum). Generally, these effects are stronger with the modern analogue technique (MAT) than with weighted averaging partial least squares (WA-PLS). With respect to fossil spectra from northern China, the spatial extent of calibration-sets should be restricted to radii between ca. 1000 and 1500 km because small-scale calibration-sets (<800 km radius) will likely fail to include enough spatial variation in the modern pollen assemblages to reflect the temporal range shifts during the Holocene, while too broad a scale calibration-set (>1500 km radius) will include taxa with very different pollen-climate relationships. (C) 2017 Elsevier Ltd. All rights reserved.
How fast the Northern Hemisphere (NH) forest biome tracks strongly warming climates is largely unknown. Regional studies reveal lags between decades and millennia. Here we report a conundrum: Deglacial forest expansion in the NH extra-tropics occurs approximately 4000 years earlier in a transient MPI-ESM1.2 simulation than shown by pollen-based biome reconstructions. Shortcomings in the model and the reconstructions could both contribute to this mismatch, leaving the underlying causes unresolved. The simulated vegetation responds within decades to simulated climate changes, which agree with pollen-independent reconstructions. Thus, we can exclude climate biases as main driver for differences. Instead, the mismatch points at a multi-millennial disequilibrium of the NH forest biome to the climate signal. Therefore, the evaluation of time-slice simulations in strongly changing climates with pollen records should be critically reassessed. Our results imply that NH forests may be responding much slower to ongoing climate changes than Earth System Models predict. <br /> Deglacial forest expansion in the Northern Hemisphere poses a conundrum: Model results agree with the climate signal but are several millennia ahead of reconstructed forest dynamics. The underlying causes remain unsolved.
Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.
Changes in species’ distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species’ range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.
The late-Quaternary climate history of monsoonal Central Asia was inferred from 75 palaeoclimatic records which provide information about moisture conditions in the last 50 ka (or part of this period). Wet conditions occurred during middle and late Marine Isotope Stage 3, while the Last Glacial Maximum (LGM) was characterized by dry climate conditions in the region. A stepwise climate amelioration is suggested by the climate records following the LGM. Several climate signals of this period, which were reported from high-latitude ice core records, are preserved in archives from monsoonal Central Asia as well. During the early Holocene, high effective moisture was inferred from most records from the area dominated by the Indian Monsoon (e.g. the Tibetan Plateau) suggesting that Holocene optimal climate conditions occurred there during this period. In contrast, areas which are dominated by the South-East Asian monsoon (SE Monsoon) and the Westerlies (in north-western and north-central China, Mongolia) do not uniformly show an early Holocene climate optimum. For this area optimal conditions prevailed during the mid-Holocene. These apparent contradictions can possibly be explained by the regional uplift and descent of air masses in the Holocene. During the early Holocene, strengthened insolation possibly caused an enhanced low-level convergence over the Tibetan Plateau which led to the intensification of the summer monsoon. The strong air uplift caused intensified precipitation and air divergence in the upper troposphere over the Tibetan Plateau. The areas adjacent to the north therefore experienced an intensified descent of air masses and consequently increased aridity. The majority of the palaeoclimatic records suggest reduced effective moisture since the late Holocene in the region.
Future precipitation levels remain uncertain because climate models have struggled to reproduce observed variations in temperature-precipitation correlations. Our analyses of Holocene proxy-based temperature-precipitation correlations and hydrological sensitivities from 2,237 Northern Hemisphere extratropical pollen records reveal a significant latitudinal dependence and temporal variations among the early, middle, and late Holocene. These proxy-based variations are largely consistent with patterns obtained from transient climate simulations (TraCE21k). While high latitudes and subtropical monsoon areas show mainly stable positive correlations throughout the Holocene, the mid-latitude pattern is temporally and spatially more variable. In particular, we identified a reversal from positive to negative temperature-precipitation correlations in the eastern North American and European mid-latitudes from the early to mid-Holocene that mainly related to slowed down westerlies and a switch to moisture-limited convection under a warm climate. Our palaeoevidence of past temperature-precipitation correlation shifts identifies those regions where simulating past and future precipitation levels might be particularly challenging.
Recent and subfossil pollen spectra from the Alashan Plateau are presented in order to provide information on desert plant representation and on recent changes in vegetation and climate in this remote area in northern China. The desert vegetation composition is faithfully represented by the surface pollen spectra. The comparison of the desert plant species to the related pollen taxa yielded the following sequence from over-representation to under- representation: Chenopodiaceae, Artemisia, Ephedra fragilis-type s.l., Reaumuria, Nitraria and Calligonum. A 72 cm long sediment record from a small hydrologically-closed inter-dune lake (SE Badan Jilin Sand Sea, southern Alashan Plateau) covering the past similar to 160 years (dated by(137)Cs) was analysed palynologically. Intervals of denser Artemisia coverage on the sand dunes around the lake, indicating wetter climate, occurred from the mid-1850s to the mid-1870s, during the first two decades of the 20th century and from the late 1930s to the beginning of the 1960s
A 13.94-m-long sediment core, collected from a medium-sized lake in the Qilian Mountains (NE Tibetan Plateau, China), was analysed palynologically at 81 horizons. The interpretation of indicator taxa yielded various vertical shifts of the vegetation belts. These palaeovegetation results have been checked with lake surface pollen spectra from 8 takes representing different altitudinal vegetation belts. Our main findings are the following: A short period of the late Marine Isotope Stage 3 (around similar to 46,000 yr ago) was characterized by interglacial temperature conditions with a tree line above its present-day altitude. During the LGM, the vicinity of the lake was not covered by ice but by sparse alpine vegetation and alpine deserts, indicating that the climate was colder by similar to 4-7 degrees C than today Markedly higher temperatures were inferred from higher arboreal pollen frequencies between similar to 13,000 and similar to 7000 yr ago with a Holocene temperature optimum and a maximal Picea-Betula mixed-forest expansion between similar to 9000 and similar to 7000 yr ago, when temperatures exceeded the present-day conditions by at least 1-2 degrees C. Alpine steppes and meadows and sub-alpine shrub vegetation dominated around the lake since the middle Holocene, suggesting that vegetation and climate conditions were exceptionally stable in comparison to previous periods.
Lakes cover large parts of the climatically sensitive Arctic landscape and respond rapidly to environmental change. Arctic lakes have different origins and include the predominant thermokarst lakes, which are small, young and highly dynamic, as well as large, old and stable glacial lakes. Freshwater diatoms dominate the primary producer community in these lakes and can be used to detect biotic responses to climate and environmental change. We used specific diatom metabarcoding on sedimentary DNA, combined with next-generation sequencing and diatom morphology, to assess diatom diversity in five glacial and 15 thermokarst lakes within the easternmost expanse of the Siberian treeline ecotone in Chukotka, Russia. We obtained 163 verified diatom sequence types and identified 176 diatom species morphologically. Although there were large differences in taxonomic assignment using the two approaches, they showed similar high abundances and diversity of Fragilariceae and Aulacoseiraceae. In particular, the genetic approach detected hidden within-lake variations of fragilarioids in glacial lakes and dominance of centric Aulacoseira species, whereas Lindavia ocellata was predominant using morphology. In thermokarst lakes, sequence types and valve counts also detected high diversity of Fragilariaceae, which followed the vegetation gradient along the treeline. Ordination analyses of the genetic data from glacial and thermokarst lakes suggest that concentrations of sulfate (SO42-), an indicator of the activity of sulfate-reducing microbes under anoxic conditions, and bicarbonate (HCO3-), which relates to surrounding vegetation, have a significant influence on diatom community composition. For thermokarst lakes, we also identified lake depth as an important variable, but SO42- best explains diatom diversity derived from genetic data, whereas HCO3- best explains the data from valve counts. Higher diatom diversity was detected in glacial lakes, most likely related to greater lake age and different edaphic settings, which gave rise to diversification and endemism. In contrast, small, dynamic thermokarst lakes are inhabited by stress-tolerant fragilarioids and are related to different vegetation types along the treeline ecotone. Our study demonstrated that genetic investigations of lake sediments can be used to interpret climate and environmental responses of diatoms. It also showed how lake type affects diatom diversity, and that such genetic analyses can be used to track diatom community changes under ongoing warming in the Arctic.
Variations in regional temperature have widespread implications for society, but our understanding of the amplitude and origin of long-term natural variability is insufficient for accurate regional projections. This is especially the case for terrestrial temperature variability, which is currently thought to be weak over long timescales. By performing spectral analysis on climate reconstructions, produced using sedimentary pollen records from the Northern Hemisphere over the last 8,000 years, coupled with instrumental data, we provide a comprehensive estimate of regional temperature variability from annual to millennial timescales. We show that short-term random variations are overprinted by strong ocean-driven climate variability on multi-decadal and longer timescales. This may cause substantial and potentially unpredictable regional climatic shifts in the coming century, in contrast to the relatively muted and homogeneous warming projected by climate models. Due to the marine influence, regions characterized by stable oceanic climate at sub-decadal timescales experience stronger long-term variability, and continental regions with higher sub-decadal variability show weaker long-term variability. This fundamental relationship between the timescales provides a unique insight into the emergence of a marine-driven low-frequency regime governing terrestrial climate variability and sets the basis to project the amplitude of temperature fluctuations on multi-decadal timescales and longer.
Temperature variability over land is enhanced by ocean temperature fluctuations on millennial timescales, with implications for regional-scale climate change, according to an analysis of Northern Hemisphere proxy records and observations.
Alpine ecosystems on the Tibetan Plateau are being threatened by ongoing climate warming and intensified human activities. Ecological time-series obtained from sedimentary ancient DNA (sedaDNA) are essential for understanding past ecosystem and biodiversity dynamics on the Tibetan Plateau and their responses to climate change at a high taxonomic resolution. Hitherto only few but promising studies have been published on this topic. The potential and limitations of using sedaDNA on the Tibetan Plateau are not fully understood. Here, we (i) provide updated knowledge of and a brief introduction to the suitable archives, region-specific taphonomy, state-of-the-art methodologies, and research questions of sedaDNA on the Tibetan Plateau; (ii) review published and ongoing sedaDNA studies from the Tibetan Plateau; and (iii) give some recommendations for future sedaDNA study designs. Based on the current knowledge of taphonomy, we infer that deep glacial lakes with freshwater and high clay sediment input, such as those from the southern and southeastern Tibetan Plateau, may have a high potential for sedaDNA studies. Metabarcoding (for microorganisms and plants), metagenomics (for ecosystems), and hybridization capture (for prehistoric humans) are three primary sedaDNA approaches which have been successfully applied on the Tibetan Plateau, but their power is still limited by several technical issues, such as PCR bias and incompleteness of taxonomic reference databases. Setting up high-quality and open-access regional taxonomic reference databases for the Tibetan Plateau should be given priority in the future. To conclude, the archival, taphonomic, and methodological conditions of the Tibetan Plateau are favorable for performing sedaDNA studies. More research should be encouraged to address questions about long-term ecological dynamics at ecosystem scale and to bring the paleoecology of the Tibetan Plateau into a new era.
Arctic treelines are facing a strong temperature increase as a result of recent global warming, causing possible changes in forest extent, which will alter vegetation-climate feedbacks. However, the mode and strength of the response is rather unclear, as potential changes are happening in areas that are very remote and difficult to access, and empirical data are still largely lacking. Here, we assessed the current population structure and genetic differentiation of Larix Mill. tree stands within the northernmost latitudinal treeline reaching ~ 72° N in the southern lowlands of the Taymyr Peninsula (~ 100° E). We sampled 743 individuals belonging to different height classes (seedlings, saplings, trees) at 11 locations along a gradient from ‘single tree’ tundra over ‘forest line’ to ‘dense forest’ stands and conducted investigations applying eight highly polymorphic nuclear microsatellites. Results suggest a high diversity within sub-populations (HE = 0.826–0.893), coupled, however, with heterozygote deficits in all sub-populations, but pronounced in ‘forest line’ stands. Overall, genetic differentiation of sub-populations is low (FST = 0.005), indicating a region-wide high gene flow, although ‘forest line’ stands harbour few rare and private alleles, likely indicating greater local reproduction. ‘Single tree’ stands, located beyond the northern forest line, are currently not involved in treeline expansion, but show signs of a long-term refuge, namely asexual reproduction and change of growth-form from erect to creeping growth, possibly having persisted for thousands of years. The lack of differentiation between the sub-populations points to a sufficiently high dispersal potential, and thus a rapid northward migration of the Siberian arctic treeline under recent global warming seems potentially unconstrained, but observations show it to be unexpectedly slow.
A strong temperature increase in the Arctic is expected to lead to latitudinal treeline shift. This tundra-taiga turnover would cause a positive vegetation-climate feedback due to albedo decrease. However, reliable estimates of tree migration rates are currently lacking due to the complex processes involved in forest establishment, which depend strongly on seed dispersal. We aim to fill this gap using LAVESI, an individual-based and spatially explicit Larix vegetation simulator. LAVESI was designed to simulate plots within homogeneous forests. Here, we improve the implementation of the seed dispersal function via field-based investigations. We inferred the effective seed dispersal distances of a typical open-forest stand on the southern Taymyr Peninsula (northern central Siberia) from genetic parentage analysis using eight nuclear microsatellite markers. The parentage analysis gives effective seed dispersal distances (median similar to 10 m) close to the seed parents. A comparison between simulated and observed effective seed dispersal distances reveals an overestimation of recruits close to the releasing tree and a shorter dispersal distance generally. We thus adapted our model and used the newly parameterised version to simulate south-to-north transects; a slow-moving treeline front was revealed. The colonisation of the tundra areas was assisted by occasional long-distance seed dispersal events beyond the treeline area. The treeline (similar to 1 tree ha(-1)) advanced by similar to 1.6 m yr(-1), whereas the forest line (similar to 100 trees ha(-1)) advanced by only similar to 0.6 m yr(-1). We conclude that the treeline in northern central Siberia currently lags behind the current strong warming and will continue to lag in the near future.
The biodiversity of tundra areas in northern high latitudes is threatened by invasion of forests under global warming. However, poorly understood nonlinear responses of the treeline ecotone mean the timing and extent of tundra losses are unclear, but policymakers need such information to optimize conservation efforts. Our individual-based model LAVESI, developed for the Siberian tundra-taiga ecotone, can help improve our understanding. Consequently, we simulated treeline migration trajectories until the end of the millennium, causing a loss of tundra area when advancing north. Our simulations reveal that the treeline follows climate warming with a severe, century-long time lag, which is overcompensated by infilling of stands in the long run even when temperatures cool again. Our simulations reveal that only under ambitious mitigation strategies (relative concentration pathway 2.6) will ~30% of original tundra areas remain in the north but separated into two disjunct refugia.
The occurrence of refugia beyond the arctic treeline and genetic adaptation therein play a crucial role of largely unknown effect size. While refugia have potential for rapidly colonizing the tundra under global warming, the taxa may be maladapted to the new environmental conditions. Understanding the genetic composition and age of refugia is thus crucial for predicting any migration response.
Here, we genotype 194 larch individuals from an similar to 1.8 km(2)area in northcentral Siberia on the southern Taimyr Peninsula by applying an assay of 16 nuclear microsatellite markers. For estimating the age of clonal individuals, we counted tree rings at sections along branches to establish a lateral growth rate that was then combined with geographic distance.
Findings reveal that the predominant reproduction type is clonal (58.76%) by short distance spreading of ramets. One outlier of clones 1 km apart could have been dispersed by reindeer. In clonal groups and within individuals, we find that somatic mutations accumulate with geographic distance. Clonal groups of two or more individuals are observed. Clonal age estimates regularly suggest individuals as old as 2,200 years, which coincides with a major environmental change that forced a treeline retreat in the region.
We conclude that individuals with clonal growth mode were naturally selected as it lowers the likely risk of extinction under a harsh environment. We discuss this legacy from the past that might now be a maladaptation and hinder expansion under currently strongly increasing temperatures.
We present a chronology framework named LegacyAge 1.0 containing harmonized chronologies for 2831 pollen records (downloaded from the Neotoma Paleoecology Database and the supplementary Asian datasets) together with their age control points and metadata in machine-readable data formats.
All chronologies use the Bayesian framework implemented in Bacon version 2.5.3. Optimal parameter settings of priors (accumulation.shape, memory.strength, memory.mean, accumulation.rate, and thickness) were identified based on information in the original publication or iteratively after preliminary model inspection.
The most common control points for the chronologies are radiocarbon dates (86.1 %), calibrated by the latest calibration curves (IntCal20 and SHCal20 for the terrestrial radiocarbon dates in the Northern Hemisphere and Southern Hemisphere and Marine20 for marine materials).
The original publications were consulted when dealing with outliers and inconsistencies. Several major challenges when setting up the chronologies included the waterline issue (18.8% of records), reservoir effect (4.9 %), and sediment deposition discontinuity (4.4 %).
Finally, we numerically compare the LegacyAge 1.0 chronologies to those published in the original publications and show that the reliability of the chronologies of 95.4% of records could be improved according to our assessment.
Our chronology framework and revised chronologies provide the opportunity to make use of the ages and age uncertainties in synthesis studies of, for example, pollen-based vegetation and climate change.
The LegacyAge 1.0 dataset, including metadata, datings, harmonized chronologies, and R code used, is openaccess and available at PANGAEA (https://doi.org/10.1594/PANGAEA.933132; Li et al., 2021) and Zenodo (https://doi.org/10.5281/zenodo.5815192; Li et al., 2022), respectively.
Pollen-based biome reconstruction on the Qinghai-Tibetan Plateau during the past 15,000 years
(2022)
Reconstruction of past vegetation change is critical for better understanding the potential impact of future global change on the fragile alpine ecosystems of the Qinghai-Tibetan Plateau (QTP). In this paper, pollen assemblages comprising 58 records from the QTP, spanning the past 15 kyrs, were collected to reconstruct biome compositions using a standard approach. Six forest biomes were identified mainly on the southeastern plateau, exhibiting a pattern of gradual expansion along the eastern margin during early to mid-Holocene times. The alpine meadow biome was separately identified based on an updated scheme, and showed notable westward expansions towards lower latitudes and higher altitudes during early Holocene times. Consistent patterns of migration could also be identified for the alpine steppe biome, which moved eastward during the late Holocene after 4 ka. As the dominant biome type, temperate steppe was distributed widely over the QTP with minor migration patterns, except for a progressive expansion to lower altitudes in the late Holocene times. The desert biome was inferred mainly as covering the northwestern plateau and the Qaidam Basin, in relatively restricted areas. The spatial distribution of the reconstructed biomes represent the large-scale vegetation gradient on the QTP. Monsoonal precipitation expressed predominant controls on the development of alpine ecosystems, while the variations in desert vegetation responded to regional moisture brought by the mid-latitude Westerlies. Temperature changes played relatively minor roles in the variations of alpine vegetation, but exerted more significant impacts on the forest biomes.
Although sedimentary ancient DNA (sedaDNA) has been increasingly used to study paleoecological dynamics (Schulte et al., 2020), the approach has rarely been compared with the traditional method of pollen analysis for investigating past changes in the vegetation composition and diversity of Arctic treeline areas. Here, we provide a history of latitudinal floristic composition and species diversity based on a comparison ofsedaDNA and pollen data archived in three Siberian lake sediment cores spanning the mid-Holocene to the present (7.6-0 cal ka BP), from northern typical tundra to southern open larch forest in the Omoloy region. Our results show that thesedaDNA approach identifies more plant taxa found in the local vegetation communities, while the corresponding pollen analysis mainly captures the regional vegetation development and has its limitations for plant diversity reconstruction. Measures of alpha diversity were calculated based onsedaDNA data recovered from along a tundra to forest tundra to open larch forest gradient. Across all sites,sedaDNA archives provide a complementary record of the vegetation transition within each lake's catchment, tracking a distinct latitudinal vegetation type range from larch tree/alder shrub (open larch forest site) to dwarf shrub-steppe (forest tundra) to wet sedge tundra (typical tundra site). By contrast, the pollen data reveal an open landscape, which cannot distinguish the temporal changes in compositional vegetation for the open larch forest site and forest-tundra site. IncreasingLarixpollen percentages were recorded in the forest-tundra site in the last millenium although noLarixDNA was detected, suggesting that thesedaDNA approach performs better for tracking the local establishment ofLarix. Highest species richness and diversity are found in the mid-Holocene (before 4.4 ka) at the typical tundra site with a diverse range of vegetational habitats, while lowest species richness is recorded for the forest tundra where dwarf-willow habitats dominated the lake's catchment. During the late Holocene, strong declines in species richness and diversity are found at the typical tundra site with the vegetation changing to relatively simple communities. Nevertheless, plant species richness is mostly higher than at the forest-tundra site, which shows a slightly decreasing trend. Plant species richness at the open larch forest site fluctuates through time and is higher than the other sites since around 2.5 ka. Taken together, there is no evidence to suggest that the latitudinal gradients in species diversity changes are present at a millennial scale. Additionally, a weak correlation between the principal component analysis (PCA) site scores ofsedaDNA and species richness suggests that climate may not be a direct driver of species turnover within a lake's catchment. Our data suggest thatsedaDNA and pollen have different but complementary abilities for reconstructing past vegetation and species diversity along a latitude.
Patterns of past vegetation changes over time and space can help facilitate better understanding of the interactions among climate, ecosystem, and human impact. Biome changes in China over the last 22,000 yr (calibrated radiocarbon date, a BP) were numerically reconstructed by using a standard approach of pollen-plant functional type-biome assignment (biomization). The biomization procedure involves pollen data from 2434 surface sites and 228 fossil sites with a high quality of pollen count and C-14 dating, 51 natural and three anthropogenic plant functional types (PFTs), as well as 19 natural and one anthropogenic biome. Surface pollen-based reconstruction of modern natural biome patterns is in good agreement (74.4%) with actual vegetation distribution in China. However, modem large-scale anthropogenic biome reconstruction has not been successful based on the current setup of three anthropogenic PFTs (plantation, secondary, and disturbed PFT) because of the limitation of non-species level pollen identification and the difficulty in the clear assignment of disturbed PFTs. The non-anthropogenic biome distributions of 44 time slices at 500-year intervals show large-scale discrepant and changed vegetation patterns from the last glacial maximum (LGM) to the Holocene throughout China. From 22 ka BP to 19 ka BP, temperate grassland, xerophytic shrubland, and desert dominated northern China, whereas cold or cool forests flourished in central China. Warm-temperate evergreen forests were restricted to far southern China, and tropical forests were absent During 18.5 ka BP to 12 ka BP, cold, cool, and dry biomes extended to some parts of northern, westem, and eastern China. Warm-temperate evergreen and mixed forests gradually expanded to occupy the whole of southern China. A slight northward shift of forest biomes occurred from 15 ka BP to 12 lea BP. During 11.5 ka BP to 9 ka BP, temperate grassland and shrubland gradually stretched to northern and western China. Cold and cool forests widely expanded into northern and central China, as well as in the northern margin of South China along with temperate deciduous forest. Since the early mid-Holocene (approximately 8.5 ka BP to 5.5 ka BP), all forest biomes shifted northward at the expense of herbaceous and shrubby biomes. Simultaneously, cold and cool forest biomes occupied the marginal areas of the Tibetan Plateau and the high mountains in western China. During the middle to late Holocene, from 5 ka to the present, temperate grassland and xerophytic shrubland expanded to the south and east, whereas temperate deciduous forests slightly shifted southward. After 3 lea BP, forest biomes were absent in western China and on the Tibetan plateau surface. Dramatic biome shifts from the LGM to the Holocene were observed in the forest-grassland ecotone and transitional zones between temperate and subtropical climates, between subtropical and tropical regions, and in the mountainous margins of the eastern Tibetan Plateau. Evidence showed more human disturbances during the late Holocene. More pollen records and historical documents are therefore further needed to understand fully the human disturbance-induced large-scale forest changes. In addition, more classifications of anthropogenic biome or land cover, more distinct assignment of pollen taxa to anthropogenic PFTs, and more effective numerical and/or mechanistic techniques in building large-scale human disturbances are required. (C) 2014 Elsevier B.V. All rights reserved.
Simulating biome distribution on the Tibetan Plateau using a modified global vegetation model
(2011)
We used a regionally modified global vegetation model (BIOME4-Tibet) to simulate biome distribution on the Tibetan Plateau under current climate conditions derived from regional meteorological observations. The bioclimatic limits (mean temperatures of the coldest and warmest months, minimum temperature, growing degree-days on 5 degrees C and 0 degrees C bases) for some key alpine plant functional types (temperate deciduous and conifer trees, boreal deciduous and conifer trees, desert woody plants, tundra shrubs, cold herbaceous plants, and lichens/forbs) were redefined based on regional vegetation-climate relationships. Modern vegetation maps confirmed that the BIOME4-Tibet model does a better job of simulating biome patterns on the plateau (gridcell agreement 52%) than the original BIOME4 model (35%). This improved model enhanced our ability to simulate temperate conifer forest, cool conifer and mixed forest, evergreen taiga, temperate xerophytic shrubland, temperate grassland and desert, and steppe and shrub tundra biomes, but made a negligible or reduced difference to the prediction of temperate deciduous forest, warm-temperate mixed forest, and three tundra biomes (erect dwarf-shrub tundra, prostrate dwarf-shrub tundra, and cushion forb, lichen, and moss tundra). Future modification of the vegetation model, by increasing the number of shrub and herb plant functional types, re-parameterization of more precise bioclimatic constraints, and improved representation of soil, permafrost, and snow processes, will be needed to better characterize the distribution of alpine vegetation on the Tibetan Plateau.
Reliable information on past and present vegetation is important to project future changes, especially for rapidly transitioning areas such as the boreal treeline. To study past vegetation, pollen analysis is common, while current vegetation is usually assessed by field surveys. Application of detailed sedimentary DNA (sedDNA) records has the potential to enhance our understanding of vegetation changes, but studies systematically investigating the power of this proxy are rare to date. This study compares sedDNA metabarcoding and pollen records from surface sediments of 31 lakes along a north-south gradient of increasing forest cover in northern Siberia (Taymyr peninsula) with data from field surveys in the surroundings of the lakes. sedDNA metabarcoding recorded 114 plant taxa, about half of them to species level, while pollen analyses identified 43 taxa, both exceeding the 31 taxa found by vegetation field surveys. Increasing Larix percentages from north to south were consistently recorded by all three methods and principal component analyses based on percentage data of vegetation surveys and DNA sequences separated tundra from forested sites. Comparisons of the ordinations using procrustes and protest analyses show a significant fit among all compared pairs of records. Despite similarities of sedDNA and pollen records, certain idiosyncrasies, such as high percentages of Alnus and Betula in all pollen and high percentages of Salix in all sedDNA spectra, are observable. Our results from the tundra to single-tree tundra transition zone show that sedDNA analyses perform better than pollen in recording site-specific richness (i.e., presence/absence of taxa in the vicinity of the lake) and perform as well as pollen in tracing vegetation composition.
Increasing arctic coastal erosion rates imply a greater release of sediments and organic matter into the coastal zone. With 213 sediment samples taken around Herschel Island-Qikiqtaruk, Canadian Beaufort Sea, we aimed to gain new insights on sediment dynamics and geochemical properties of a shallow arctic nearshore zone. Spatial characteristics of nearshore sediment texture (moderately to poorly sorted silt) are dictated by hydrodynamic processes, but ice-related processes also play a role. We determined organic matter (OM) distribution and inferred the origin and quality of organic carbon by C/N ratios and stable carbon isotopes delta C-13. The carbon content was higher offshore and in sheltered areas (mean: 1.0 wt.%., S.D.: 0.9) and the C/N ratios also showed a similar spatial pattern (mean: 11.1, S.D.: 3.1), while the delta C-13 (mean: -26.4 parts per thousand VPDB, S.D.: 0.4) distribution was more complex. We compared the geochemical parameters of our study with terrestrial and marine samples from other studies using a bootstrap approach. Sediments of the current study contained 6.5 times and 1.8 times less total organic carbon than undisturbed and disturbed terrestrial sediments, respectively. Therefore, degradation of OM and separation of carbon pools take place on land and continue in the nearshore zone, where OM is leached, mineralized, or transported beyond the study area.
We investigated 4.84-m-long sediment record spanning over the Late Glacial and Holocene from Lake Donggi Cona to be able to reconstruct circulation pattern on the Tibetan Plateau (TP). Presently, Lake Donggi Cona is located at the boundaries of Westerlies and Asian monsoon circulations in the northeastern TP. However, the exact timing and stimulating mechanisms for climatic changes and monsoon shifts in this region are still debated. We used a 19-ka-long stable isotope record of sedimentary n-alkanes to address this discrepancy by providing insights into paleohydrological conditions. The SD of nC(23) is influenced by lake water evaporation; the BD. values of sedimentary nC(29) are mainly controlled by moisture source and temperature changes. Long-chain n-alkanes dominate over the core whereas three mean clusters (i.e. microbial, aquatic and terrestrial) can be inferred. Multi-proxies suggest five major episodes in the history of Lake Donggi Cona. The Lake Donggi Cona record indicates that the Late Glacial(18.4-14.8 cal ka BP) was dominated by low productivity of mainly microbial and aquatic organisms. Relatively low delta D values suggest low temperatures and moist conditions eventually caused by stronger Westerlies, winter monsoon and melt-water influence. Likely, the shift (similar to 17.9 cal ka BP) from microbial to enhanced aquatic input suggests either a change from deep to shallow water lake or a break in local stratification. Between 14.8 and 13.0 cal ka BP, variable climatic conditions prevailed. Although the Westerlies weekend, the increase in temperature enhanced the permafrost and snow melting (displayed by a high sedimentary accumulation rate). Higher delta D values indicate increasingly arid conditions with higher temperatures which eventually lead to high evaporative conditions and lowest lake levels. Low vegetation cover and high erosion rates led to high sediment accumulation resulting in stratification followed by anoxia in the terminal lake. From 13.0 to 9.2 cal ka BP, lowered values of 813 along with high contents of terrestrial organic matter marked the early-Holocene warming indicating a further strengthening of summer precipitation and higher lake levels. A cooling trend was observed in the mid-Holocene between 9.2 and 3.0 cal ka BP accompanied by higher moisture availability (displayed by lowered SD values) caused by reduced evaporative conditions due to a drop in temperature and recovering Westerlies. After 3.0 cal ka BP, a decrease in lake productivity and cold and semi-arid conditions prevailed suggesting lower lake levels and reduced moisture from recycled air masses and Westerlies. We propose that the summer monsoon was the predominant moisture source during the Belling-Allered warm complex and early -Holocene followed by Westerlies in mid-to-late Holocene period. Stable carbon isotope values-32%o indicate the absence of C-4 -type vegetation in the region contradicting with their presence in the Lake Qinghai record. The 81) record from lake Donggi Cona highlights the importance of the interplay between Westerlies and summer monsoon circulation at this location, which is highly dynamic in northeastern plateau compared to the North Atlantic circulation and insolation changes. Consequently lake Donggi Cona might be an important anchor point for environmental reconstructions on the Tibetan Plateau. (C) 2017 The Authors. Published by Elsevier Ltd.
Siberian larch (Larix Mill.) forests dominate vast areas of northern Russia and contribute important ecosystem services to the world. It is important to understand the past dynamics of larches in order to predict their likely response to a changing climate in the future. Sedimentary ancient DNA extracted from lake sediment cores can serve as archives to study past vegetation. However, the traditional method of studying sedimentary ancient DNA-metabarcoding-focuses on small fragments, which cannot resolve Larix to species level nor allow a detailed study of population dynamics. Here, we use shotgun sequencing and hybridization capture with long-range PCR-generated baits covering the complete Larix chloroplast genome to study Larix populations from a sediment core reaching back to 6700 years from the Taymyr region in northern Siberia. In comparison with shotgun sequencing, hybridization capture results in an increase in taxonomically classified reads by several orders of magnitude and the recovery of complete chloroplast genomes of Larix. Variation in the chloroplast reads corroborates an invasion of Larix gmelinii into the range of Larix sibirica before 6700 years ago. Since then, both species have been present at the site, although larch populations have decreased with only a few trees remaining in what was once a forested area. This study demonstrates for the first time that hybridization capture applied directly to ancient DNA of plants extracted from lake sediments can provide genome-scale information and is a viable tool for studying past genomic changes in populations of single species, irrespective of a preservation as macrofossil.
Aim: The continental-scale distribution of plant functional types, such as evergreen and summergreen needle-leaf forest, is assumed to be determined by contemporary climate. However, the distribution of summergreen needle-leaf forest of larch (Larix Mill.) differs markedly between the continents, despite relatively similar climatic conditions. The reasons for these differences are little understood. Our aim is to identify potential triggers and drivers of the current distribution patterns by comparing species' bioclimatic niches, glacial refugia and postglacial recolonization patterns.
Location: Northern hemisphere.
Taxon: Species of the genus Larix (Mill.).
Methods: We compare species distribution and dominance using species ranges and sites of dominance, as well as their occurrence on modelled permafrost extent, and active layer thickness (ALT). We compare the bioclimatic niches and calculate the niche overlap between species, using the same data in addition to modern climate data. We synthesize pollen, macrofossil and ancient DNA palaeo-evidence of past Larix occurrences of the last 60,000 years and track differences in distribution patterns through time.
Results: Bioclimatic niches show large overlaps between Asian larch species and American Larix laricina. The distribution across various degrees of permafrost extent is distinctly different for Asian L. gmelinii and L. cajanderi compared to the other species, whereas the distribution on different depths of ALT is more similar among Asian and American species. Northern glacial refugia for Larix are only present in eastern Asia and Alaska.
Main Conclusion: The dominance of summergreen larches in Asia, where evergreen conifers dominate most of the rest of the boreal forests, is dependent on the interaction of several factors which allows Asian L. gmelinii and L. cajanderi to dominate where these factors coincide. These factors include the early postglacial spread out of northern glacial refugia in the absence of competitors as well as a positive feedback mechanism between frozen ground and forest.
Diatom diversity in lakes of northwest Yakutia (Siberia) was investigated by microscopic and genetic analysis of surface and cored lake sediments, to evaluate the use of sedimentary DNA for paleolimnological diatom studies and to identify obscure genetic diversity that cannot be detected by microscopic methods. Two short (76 and 73 bp) and one longer (577 bp) fragments of the ribulose 1,5-bisphosphate carboxylase/oxygenase (rbcL) gene, encoding the large subunit of the rbcL, were used as genetic markers. Diverse morphological assemblages of diatoms, dominated by small benthic fragilarioid taxa, were retrieved from the sediments of each lake. These minute fragilarioid taxa were examined by scanning electron microscopy, revealing diverse morphotypes in Staurosira and Staurosirella from the different lakes. Genetic analyses indicated a dominance of haplotypes that were assigned to fragilarioid taxa and less genetic diversity in other diatom taxa. The long rbcL_577 amplicon identified considerable diversification among haplotypes clustering within the Staurosira/Staurosirella genera, revealing 19 different haplotypes whose spatial distribution appears to be primarily related to the latitude of the lakes, which corresponds to a vegetation and climate gradient. Our rbcL markers are valuable tools for tracking differences between diatom lineages that are not visible in their morphologies. These markers revealed putatively high genetic diversity within the Staurosira/Staurosirella species complex, at a finer scale than is possible to resolve by microscopic determination. The rbcL markers may provide additional reliable information on the diversity of barely distinguishable minute benthic fragilarioids. Environmental sequencing may thus allow the tracking of spatial and temporal diversification in Siberian lakes, especially in the context of diatom responses to recent environmental changes, which remains a matter of controversy.
Genetic investigations on eukaryotic plankton confirmed the existence of modern biogeographic patterns, but analyses of palaeoecological data exploring the temporal variability of these patterns have rarely been presented. Ancient sedimentary DNA proved suitable for investigations of past assemblage turnover in the course of environmental change, but genetic relatedness of the identified lineages has not yet been undertaken. Here, we investigate the relatedness of diatom lineages in Siberian lakes along environmental gradients (i.e. across treeline transects), over geographic distance and through time (i.e. the last 7000 years) using modern and ancient sedimentary DNA. Our results indicate that closely-related Staurosira lineages occur in similar environments and less-related lineages in dissimilar environments, in our case different vegetation and co-varying climatic and limnic variables across treeline transects. Thus our study reveals that environmental conditions rather than geographic distance is reflected by diatom-relatedness patterns in space and time. We tentatively speculate that the detected relatedness pattern in Staurosira across the treeline could be a result of adaptation to diverse environmental conditions across the arctic boreal treeline, however, a geographically-driven divergence and subsequent repopulation of ecologically different habitats might also be a potential explanation for the observed pattern.
Shrub encroachment has far-reaching ecological and economic consequences in many ecosystems worldwide. Yet, compositional changes associated with shrub encroachment are often overlooked despite having important effects on ecosystem functioning. We document the compositional change and potential drivers for a northern Namibian Combretum woodland transitioning into a Terminalia shrubland. We use a multiproxy record (pollen, sedimentary ancient DNA, biomarkers, compound-specific carbon (delta C-13) and deuterium (delta D) isotopes, bulk carbon isotopes (delta(13)Corg), grain size, geochemical properties) from Lake Otjikoto at high taxonomical and temporal resolution. We provide evidence that state changes in semiarid environments may occur on a scale of one century and that transitions between stable states can span around 80 years and are characterized by a unique vegetation composition. We demonstrate that the current grass/woody ratio is exceptional for the last 170 years, as supported by n-alkane distributions and the delta C-13 and delta(13)Corg records. Comparing vegetation records to environmental proxy data and census data, we infer a complex network of global and local drivers of vegetation change. While our delta D record suggests physiological adaptations of woody species to higher atmospheric pCO(2) concentration and drought, our vegetation records reflect the impact of broad-scale logging for the mining industry, and the macrocharcoal record suggests a decrease in fire activity associated with the intensification of farming. Impact of selective grazing is reflected by changes in abundance and taxonomical composition of grasses and by an increase of nonpalatable and trampling-resistant taxa. In addition, grain-size and spore records suggest changes in the erodibility of soils because of reduced grass cover. Synthesis. We conclude that transitions to an encroached savanna state are supported by gradual environmental changes induced by management strategies, which affected the resilience of savanna ecosystems. In addition, feedback mechanisms that reflect the interplay between management legacies and climate change maintain the encroached state.
Understanding vegetation-modern pollen relationships is essential to provide confidence in fossil pollen reconstructions of long-term vegetation changes in savanna ecosystems. In this paper we compare the taxonomical composition and the diversity (Hill NO, N1, N2) of vegetation and modern pollen along precipitation and local grazing-intensity gradients in Namibian savannas. Modern pollen was extracted from surface soil samples collected from 5 x 5 m plots distributed along four 500 m gradients. Vegetation was surveyed in each plot. The results show a high correspondence between vegetation and pollen data in terms of composition. Precipitation and grazing explain a significant although low proportion of compositional change in the vegetation and pollen spectra. We identified pollen taxa as indicators of grazing pressure such as Limeum, Alternanthera, and particularly Tribulus. Correspondence between vegetation and pollen data in terms of taxa richness (NO) is limited, probably because of the influence of landscape heterogeneity and openness, as well as low pollen concentrations. In contrast, the effective numbers of common and dominant taxa (N1, N2) are consistent among the different datasets. We conclude that in spite of limitations, modern pollen assemblages can reflect changes in vegetation composition, richness and diversity patterns along precipitation and grazing gradients in savanna environments. (C) 2018 Elsevier B.V. All rights reserved.
Arctic and alpine treelines worldwide differ in their reactions to climate change. A northward advance of or densification within the treeline ecotone will likely influence climate-vegetation feedback mechanisms. In our study, which was conducted in the Taimyr Depression in the North Siberian Lowlands, w present a combined field-and model-based approach helping us to better understand the population processes involved in the responses of the whole treeline ecotone, spanning from closed forest to single-tree tundra, to climate warming. Using information on stand structure, tree age, and seed quality and quantity from seven sites, we investigate effects of intra-specific competition and seed availability on the specific impact of recent climate warming on larch stands. Field data show that tree density is highest in the forest-tundra, and average tree size decreases from closed forest to single-tree tundra. Age-structure analyses indicate that the trees in the closed forest and forest-tundra have been present for at least similar to 240 yr. At all sites except the most southerly ones, past establishment is positively correlated with regional temperature increase. In the single-tree tundra, however, a change in growth form from krummholz to erect trees, beginning similar to 130 yr ago, rather than establishment date has been recorded. Seed mass decreases from south to north, while seed quantity increases. Simulations with LAVESI (Larix Vegetation Simulator) further suggest that relative density changes strongly in response to a warming signal in the forest-tundra while intra-specific competition limits densification in the closed forest and seed limitation hinders densification in the single-tree tundra. We find striking differences in strength and timing of responses to recent climate warming. While forest-tundra stands recently densified, recruitment is almost non-existent at the southern and northern end of the ecotone due to autecological processes. Palaeo-treelines may therefore be inappropriate to infer past temperature changes at a fine scale. Moreover, a lagged treeline response to past warming will, via feedback mechanisms, influence climate change in the future.
Mosses are a major component of the arctic vegetation, particularly in wetlands. We present C / N atomic ratio, delta C-13 and delta N-15 data of 400 brown-moss samples belonging to 10 species that were collected along hydrological gradients within polygonal mires located on the southern Taymyr Peninsula and the Lena River delta in northern Siberia. Additionally, n-alkane patterns of six of these species (16 samples) were investigated. The aim of the study is to see whether the inter-and intraspecific differences in C / N, isotopic compositions and n-alkanes are indicative of habitat, particularly with respect to water level. Overall, we find high variability in all investigated parameters for two different moisture-related groups of moss species. The C / N ratios range between 11 and 53 (median: 32) and show large variations at the intraspecific level. However, species preferring a dry habitat (xero-mesophilic mosses) show higher C / N ratios than those preferring a wet habitat (meso-hygrophilic mosses). The delta C-13 values range between 37.0 and 22.5% (median D 27.8 %). The delta N-15 values range between 6.6 and C 1.7%(median D 2.2 %). We find differences in delta C-13 and delta N-15 compositions between both habitat types. For some species of the meso-hygrophilic group, we suggest that a relationship between the individ-ual habitat water level and isotopic composition can be inferred as a function of microbial symbiosis. The n-alkane distribution also shows differences primarily between xeromesophilic and meso-hygrophilic mosses, i. e. having a dominance of n-alkanes with long (n-C29, n-C31 /and intermediate (n-C25 /chain lengths, respectively. Overall, our results reveal that C / N ratios, isotopic signals and n-alkanes of studied brown-moss taxa from polygonal wetlands are characteristic of their habitat.
Larix populations at the tundra-taiga ecotone in northern Siberia are highly under-represented in population genetic studies, possibly due to the remoteness of these regions that can only be accessed at extraordinary expense. The genetic signatures of populations in these boundary regions are therefore largely unknown. We aim to generate organelle reference genomes for the detection of single nucleotide polymorphisms (SNPs) that can be used for paleogenetic studies. We present 19 complete chloroplast genomes and mitochondrial genomic sequences of larches from the southern lowlands of the Taymyr Peninsula (northernmost range of Larix gmelinii (Rupr.) Kuzen.), the lower Omoloy River, and the lower Kolyma River (both in the range of Larix cajanderi Mayr). The genomic data reveal 84 chloroplast SNPs and 213 putatively mitochondrial SNPs. Parsimony-based chloroplast haplotype networks show no spatial structure of individuals from different geographic origins, while the mitochondrial haplotype network shows at least a slight spatial structure with haplotypes from the Omoloy and Kolyma populations being more closely related to each other than to most of the haplotypes from the Taymyr populations. Whole genome alignments with publicly available complete chloroplast genomes of different Larix species show that among official plant barcodes only the rcbL gene contains sufficient polymorphisms, but has to be sequenced completely to distinguish the different provenances. We provide 8 novel mitochondrial SNPs that are putatively diagnostic for the separation of L. gmelinii and L. cajanderi, while 4 chloroplast SNPs have the potential to distinguish the L. gmelinii/ L. cajanderi group from other Larix species. Our organelle references can be used for a targeted primer and probe design allowing the generation of short amplicons. This is particularly important with regard to future investigations of, for example, the biogeographic history of Larix by screening ancient sedimentary DNA of Larix.
Ecosystem boundaries, such as the Arctic-Boreal treeline, are strongly coupled with climate and were spatially highly dynamic during past glacial-interglacial cycles. Only a few studies cover vegetation changes since the last interglacial, as most of the former landscapes are inundated and difficult to access. Using pollen analysis and sedimentary ancient DNA (sedaDNA) metabarcoding, we reveal vegetation changes on Bol’shoy Lyakhovsky Island since the last interglacial from permafrost sediments. Last interglacial samples depict high levels of floral diversity with the presence of trees (Larix, Picea, Populus) and shrubs (Alnus, Betula, Ribes, Cornus, Saliceae) on the currently treeless island. After the Last Glacial Maximum, Larix re-colonised the island but disappeared along with most shrub taxa. This was probably caused by Holocene sea-level rise, which led to increased oceanic conditions on the island. Additionally, we applied two newly developed larch-specific chloroplast markers to evaluate their potential for tracking past population dynamics from environmental samples. The novel markers were successfully re-sequenced and exhibited two variants of each marker in last interglacial samples. SedaDNA can track vegetation changes as well as genetic changes across geographic space through time and can improve our understanding of past processes that shape modern patterns.
The Fram Strait is an area with a relatively low and irregular distribution of diatom microfossils in surface sediments, and thus microfossil records are scarce, rarely exceed the Holocene, and contain sparse information about past richness and taxonomic composition. These attributes make the Fram Strait an ideal study site to test the utility of sedimentary ancient DNA (sedaDNA) metabarcoding. Amplifying a short, partial rbcL marker from samples of sediment core MSM05/5-712-2 resulted in 95.7% of our sequences being assigned to diatoms across 18 different families, with 38.6% of them being resolved to species and 25.8% to genus level. Independent replicates show a high similarity of PCR products, especially in the oldest samples. Diatom sedaDNA richness is highest in the Late Weichselian and lowest in Mid- and Late Holocene samples. Taxonomic composition is dominated by cold-water and sea-ice-associated diatoms and suggests several reorganisations - after the Last Glacial Maximum, after the Younger Dryas, and after the Early and after the Mid-Holocene. Different sequences assigned to, amongst others, Chaetoceros socialis indicate the detectability of intra-specific diversity using sedaDNA. We detect no clear pattern between our diatom sedaDNA record and the previously published IP25 record of this core, although proportions of pennate diatoms increase with higher IP25 concentrations and proportions of Nitzschia cf. frigida exceeding 2% of the assemblage point towards past sea-ice presence.